Global solution of the ideal resonance problem

If a dynamical problem ofN degress of freedom is reduced to the Ideal Resonance Problem, the Hamiltonian takes the form 1 $$\begin{array}{*{20}c} {F = B(y) + 2\mu ^2 A(y)\sin ^2 x_1 ,} & {\mu \ll 1.} \\ \end{array} $$ Herey is the momentum-vectory _k withk=1,2?N, x ₁ is thecritical argument, andx _k fork>1 are theignorable co-ordinates, which have been eliminated from the Hamiltonian. The purpose of this Note is to summarize the first-order solution of the problem defined by (1) as described in a sequence of five recent papers by the author. A basic is the resonance parameter α, defined by 1 $$\alpha \equiv - B'/\left| {4AB''} \right|^{1/2} \mu .$$ The solution isglobal in the sense that it is valid for all values of α² in the range 1 $$0 \leqslant \alpha ^2 \leqslant \infty ,$$ which embrances thelibration and thecirculation regimes of the co-ordinatex ₁, associated with α² < 1 and α² > 1, respectively. The solution includes asymptotically the limit α² → ∞, which corresponds to theclassical solution of the problem, expanded in powers of ε ≡ μ², and carrying α as a divisor. The classical singularity at α=0, corresponding to an exact commensurability of two frequencies of the motion, has been removed from the global solution by means of the Bohlin expansion in powers of μ = ε^1/2. The singularities that commonly arise within the libration region α² < 1 and on the separatrix α² = 1 of the phase-plane have been suppressed by means of aregularizing function 1 $$\begin{array}{*{20}c} {\phi \equiv \tfrac{1}{2}(1 + \operatorname{sgn} z)\exp ( - z^{ - 3} ),} & {z \equiv \alpha ^2 } \\ \end{array} - 1,$$ introduced into the new Hamiltonian. The global solution is subject to thenormality condition, which boundsAB″ away from zero indeep resonance, α² < 1/μ, where the classical solution fails, and which boundsB′ away from zero inshallow resonance, α² > 1/μ, where the classical solution is valid. Thedemarcation point 1 $$\alpha _ * ^2 \equiv {1 \mathord{\left/ {\vphantom {1 \mu }} \right. \kern-\nulldelimiterspace} \mu }$$ conventionally separates the deep and the shallow resonance regions. The solution appears in parametric form 1 $$\begin{array}{*{20}c} {x_\kappa = x_\kappa (u)} \\ {y_1 = y_1 (u)} \\ {\begin{array}{*{20}c} {y_\kappa = conts,} & {k > 1,} \\ \end{array} } \\ {u = u(t).} \\ \end{array} $$ It involves the standard elliptic integralsu andE((u) of the first and the second kinds, respectively, the Jacobian elliptic functionssn, cn, dn, am, and the Zeta functionZ (u).     

A note on theH-functions of transfer problems in multiplying media

Some useful results and remodelled representations ofH-functions corresponding to the dispersion function $$T\left( z \right) = 1 - 2z^2 \sum\limits_1^n {\int_0^{\lambda r} {Y_r } \left( x \right){\text{d}}x/\left( {z^2 - x^2 } \right)} $$ are derived, suitable to the case of a multiplying medium characterized by $$\gamma _0 = \sum\limits_1^n {\int_0^{\lambda r} {Y_r } \left( x \right){\text{d}}x > \tfrac{1}{2} \Rightarrow \xi = 1 - 2\gamma _0< 0} $$      

Sur un modele explicitant les differents types morphologiques des galaxies

In the now classical Lindblad-Lin density-wave theory, the linearization of the collisionless Boltzmann equation is made by assuming the potential functionU expressed in the formU=U ₀ + \(\tilde U\) +... WhereU ₀ is the background axisymmetric potential and \(\tilde U<< U_0 \) . Then the corresponding density distribution is \(\rho = \rho _0 + \tilde \rho (\tilde \rho<< \rho _0 )\) and the linearized equation connecting \(\tilde U\) and the component \(\tilde f\) of the distribution function is given by $$\frac{{\partial \tilde f}}{{\partial t}} + \upsilon \frac{{\partial \tilde f}}{{\partial x}} - \frac{{\partial U_0 }}{{\partial x}} \cdot \frac{{\partial \tilde f}}{{\partial \upsilon }} = \frac{{\partial \tilde U}}{{\partial x}}\frac{{\partial f_0 }}{{\partial \upsilon }}.$$ One looks for spiral self-consistent solutions which also satisfy Poisson's equation $$\nabla ^2 \tilde U = 4\pi G\tilde \rho = 4\pi G\int {\tilde f d\upsilon .} $$ Lin and Shu (1964) have shown that such solutions exist in special cases. In the present work, we adopt anopposite proceeding. Poisson's equation contains two unknown quantities \(\tilde U\) and \(\tilde \rho \) . It could be completelysolved if a second independent equation connecting \(\tilde U\) and \(\tilde \rho \) was known. Such an equation is hopelesslyobtained by direct observational means; the only way is to postulate it in a mathematical form. In a previouswork, Louise (1981) has shown that Poisson's equation accounted for distances of planets in the solar system(following to the Titius-Bode's law revised by Balsano and Hughes (1979)) if the following relation wasassumed $$\rho ^2 = k\frac{{\tilde U}}{{r^2 }} (k = cte).$$ We now postulate again this relation in order to solve Poisson's equation. Then, $$\nabla ^2 \tilde U - \frac{{\alpha ^2 }}{{r^2 }}\tilde U = 0, (\alpha ^2 = 4\pi Gk).$$ The solution is found in a classical way to be of the form $$\tilde U = cte J_v (pr)e^{ - pz} e^{jn\theta } $$ wheren = integer,p =cte andJ _v (pr) = Bessel function with indexv (v ² =n ² + α²). By use of the Hankel function instead ofJ _v (pr) for large values ofr, the spiral structure is found to be given by $$\tilde U = cte e^{ - pz} e^{j[\Phi _v (r) + n\theta ]} , \Phi _v (r) = pr - \pi /2(v + \tfrac{1}{2}).$$ For small values ofr, \(\tilde U\) = 0: the center of a galaxy is not affected by the density wave which is onlyresponsible of the spiral structure. For various values ofp,n andv, other forms of galaxies can be taken into account: Ring, barred and spiral-barred shapes etc. In order to generalize previous calculations, we further postulateρ ₀ =kU ₀/r ², leading to Poisson'sequation which accounts for the disc population $$\nabla ^2 U_0 - \frac{{\alpha ^2 }}{{r^2 }}U_0 = 0.$$ AsU ₀ is assumed axisymmetrical, the obvious solution is of the form $$U_0 = \frac{{cte}}{{r^v }}e^{ - pz} , \rho _0 = \frac{{cte}}{{r^{2 + v} }}e^{ - pz} .$$ Finally, Poisson's equation is completely solvable under the assumptionρ =k(U/r ². The general solution,valid for both disc and spiral arm populations, becomes $$U = cte e^{ - pz} \left\{ {r^{ - v} + } \right.\left. {cte e^{j[\Phi _v (r) + n\theta ]} } \right\},$$ The density distribution along the O _z axis is supported by Burstein's (1979) observations.     

Experimental testing of relativistic effects,variability of the gravitational constant and topography of Mercury surface from radar observations 1964–1989

The new analysis of radar observations of inner planets for the time span 1964–1989 is described. The residuals show that Mercury topography is an important source of systematic errors which have not been taken into account up to now. The longitudinal and latitudinal variations of heights of Mercury surface were found and an approximate map of equatorial zone |?|≤120° was constructed. Including three values characterizing global nonsphericity of Mercury surface into the set of parameters under determination allowed to improve essentially all estimates. In particular, the variability of the gravitational constantG was evaluated: $$\dot G/G = (0.47 \pm 0.47) \times 10^{ - 11} yr^{ - 1} $$ . The correction to Mercury perihelion motion: $$\Delta \dot \pi = - 0''.017 \pm 0''.052 cy^{ - 1} $$ and linear combination of the parameters of PPN formalism: $$\upsilon = (2 + 2\gamma - \beta )/3 = 0.9995 \pm 0.0013$$ were determined; they are in a good agreement with General Relativity predictions. The obtained values Δ.π and ν correspond to the negligible solar oblateness, the estimate of solar quadrupole moment being: $$J_2 = ( - 0.13 \pm 0.41) \times 10^{ - 6} $$ .     

Static spherically-symmetric scalar-field theory in general relativity

A spherically-symmetric static scalar field in general relativity is considered. The field equations are defined by $$\begin{gathered} R_{ik} = - \mu \varphi _i \varphi _k ,\varphi _i = \frac{{\partial \varphi }}{{\partial x^i }}, \varphi ^i = g^{ik} \varphi _k , \hfill \\ \hfill \\ \end{gathered} $$ where ?=?(r,t) is a scalar field. In the past, the same problem was considered by Bergmann and Leipnik (1957) and Buchdahl (1959) with the assumption that ?=?(r) be independent oft and recently by Wyman (1981) with the assumption ?=?(r, t). The object of this paper is to give explicit results with a different approach and under a more general condition $$\phi _{;i}^i = ( - g)^{ - 1/2} \frac{\partial }{{\partial x^i }}\left[ {( - g)^{1/2} g^{ik} \frac{\partial }{{\partial x^k }}} \right] = - 4\pi ( -g )^{ - 1/2} \rho $$ where ?=?(r, t) is the mass or the charge density of the sources of the field.     

Changes in binding energy of interacting galaxies

It is shown that the fractional increase in binding energy of a galaxy in a fast collision with another galaxy of the same size can be well represented by the formula $$\xi _2 = 3({G \mathord{\left/ {\vphantom {G {M_2 \bar R}}} \right. \kern-\nulldelimiterspace} {M_2 \bar R}}) ({{M_1 } \mathord{\left/ {\vphantom {{M_1 } {V_p }}} \right. \kern-\nulldelimiterspace} {V_p }})^2 e^{ - p/\bar R} = \xi _1 ({{M_1 } \mathord{\left/ {\vphantom {{M_1 } {M_2 }}} \right. \kern-\nulldelimiterspace} {M_2 }})^3 ,$$ whereM ₁,M ₂ are the masses of the perturber and the perturbed galaxy, respectively,V _p is the relative velocity of the perturber at minimum separationp, and \(\bar R\) is the dynamical radius of either galaxy.     

Comparison of the classical and the global solutions of the Ideal Resonance Problem

The Ideal Resonance Problem is defined by the Hamiltonian $$F = B(y) + 2\varepsilon A(y) \sin ^2 x,\varepsilon \ll 1.$$ The classical solution of the Problem, expanded in powers of ε, carries the derivativeB′ as a divisor and is, therefore, singular at the zero ofB′, associated with resonance. With α denoting theresonance parameter, defined by $$\alpha \equiv - B'/|4AB''|^{1/2} \mu ,\mu = \varepsilon ^{1/2} ,$$ it is shown here that the classical solution is valid only for $$\alpha ^2 \geqslant 0(1/\mu ).$$ In contrast, the global solution (Garfinkelet al., 1971), expanded in powers ofμ=ε^1/2, removes the classical singularity atB′=0, and is valid for all α. It is also shown here that the classical solution is an asymptotic approximation, for largeα ², of the global solution expanded in powers ofα ^?2. This result leads to simplified expressions for resonancewidth and resonantamplification. The two solutions are compared with regard to their general behavior and their accuracy. It is noted that the global solution represents a perturbed simple pendulum, while the classical solution is the limiting case of a pendulum in a state offast circulation.     

Observational constraints on G-corrected holographic dark energy using a Markov chain Monte Carlo method

We constrain holographic dark energy (HDE) with time varying gravitational coupling constant in the framework of the modified Friedmann equations using cosmological data from type Ia supernovae, baryon acoustic oscillations, cosmic microwave background radiation and X-ray gas mass fraction. Applying a Markov Chain Monte Carlo (MCMC) simulation, we obtain the best fit values of the model and cosmological parameters within 1σ confidence level (CL) in a flat universe as: $\varOmega_{b}h^{2}=0.0222^{+0.0018}_{-0.0013}$ , $\varOmega_{c}h^{2}=0.1121^{+0.0110}_{-0.0079}$ , $\alpha_{G}\equiv \dot{G}/(HG) =0.1647^{+0.3547}_{-0.2971}$ and the HDE constant $c=0.9322^{+0.4569}_{-0.5447}$ . Using the best fit values, the equation of state of the dark component at the present time w _d0 at 1σ CL can cross the phantom boundary w=?1.     

Nonlinear newtonian theory and the effect of time dilatation

The fact that the energy density ρ_g of a static spherically symmetric gravitational field acts as a source of gravity, gives us a harmonic function \(f\left( \varphi \right) = e^{\varphi /c^2 } \) , which is determined by the nonlinear differential equation $$\nabla ^2 \varphi = 4\pi k\rho _g = - \frac{1}{{c^2 }}\left( {\nabla \varphi } \right)^2 $$ Furthermore, we formulate the infinitesimal time-interval between a couple of events measured by two different inertial observers, one in a position with potential φ-i.e., dt _φ and the other in a position with potential φ=0-i.e., dt ₀, as $${\text{d}}t_\varphi = f{\text{d}}t_0 .$$ When the principle of equivalence is satisfied, we obtain the well-known effect of time dilatation.     

Ultra-high energy neutrino fluxes from supermassive AGN black holes

We compute the ultra-high energy (UHE) neutrino fluxes from plausible accreting supermassive black holes closely linking to the 377 active galactic nuclei (AGNs). They have well-determined black hole masses collected from the literature. The neutrinos are produced via simple or modified URCA processes, even after the neutrino trapping, in superdense proto-matter medium. The resulting fluxes are ranging from: (1) (quark reactions)— $J^{q}_{\nu\varepsilon}/(\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1})\simeq8.29\times 10^{-16}$ to 3.18×10^?4, with the average $\overline{J}^{q}_{\nu\varepsilon}\simeq5.53\times 10^{-10}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ , where ε _d ～10^?12 is the opening parameter; (2) (pionic reactions)— $J^{\pi}_{\nu\varepsilon} \simeq0.112J^{q}_{\nu\varepsilon}$ , with the average $J^{\pi}_{\nu\varepsilon} \simeq3.66\times 10^{-11}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ ; and (3) (modified URCA processes)— $J^{URCA}_{\nu\varepsilon}\simeq7.39\times10^{-11} J^{q}_{\nu\varepsilon}$ , with the average $\overline{J}^{URCA}_{\nu\varepsilon} \simeq2.41\times10^{-20} \varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ . We conclude that the AGNs are favored as promising pure neutrino sources, because the computed neutrino fluxes are highly beamed along the plane of accretion disk, peaked at high energies and collimated in smaller opening angle θ～ε _d .     

Some identities for spheroidal harmonics

The spheroidal harmonics expressions $$\left[ {P_{2k}^{2s} \left( {i\xi } \right)P_{2k - 2r}^{2s} \left( \eta \right) - P_{2k - 2r}^{2s} \left( {i\xi } \right)P_{2k}^{2s} \left( \eta \right)} \right]e^{i2s\theta } $$ and $$\left[ {\eta ^2 P_{2k}^{2s} \left( {i\xi } \right)P_{2k - 2r}^{2s} \left( \eta \right) + \xi ^2 P_{2k - 2r}^{2s} \left( {i\xi } \right)P_{2k}^{2s} \left( \eta \right)} \right]e^{i2s\theta } $$ , have ξ²+η² as a factor. A method is presented for obtaining for these two expressions the coefficient of ξ²+η² in the form of a linear combination of terms of the formP _2m ^2s (iξ)P _2n ^2s (η)e ^i2sθ. Explicit formulae are exhibited for the casesr=1, 2, 3 and any positive or zero integersk ands. Such identities are useful in gravitational potential theory for ellipsoidal distributions when matching Legendre function expansions are employed.     

Solution of the equation of transfer for interlocked multiplets by the method of Laplace transform and Winer-Hopf technique with planck function as a nonlinear function of optical depth

The equation of transfer for interlocked multiplets has been solved by Laplace transformation and the Wiener-Hopf technique developed by Dasgupta (1978) considering two nonlinear forms of Planck function: i.e., (a) $$B{\text{ }}_{\text{v}} (T) = B(t) = b_0 + b_1 {\text{ }}e^{ - \alpha t} ,$$ (b) $$B{\text{ }}_{\text{v}} (T) = B(t) = b_0 + b_1 t + b_2 E_2 (t).$$ Solutions obtained by Dasgupta (1978) or by Chandrasekhar (1960) may be obtained from our solutions by dropping the nonlinear terms.     

The stability of differential rotation with non-axisymmetric perturbations

A plane-wave analysis on a simplified scheme based on the Boussinesq approximation and shallow convection is used to establish the necessary conditions for stability of a differentiallyrotating, compressible flow between two coaxial cylinders subject to non-axisymmetric perturbations. To test the adequateness of this simplification, the sufficient conditions for stability are again established which agree with those obtained by a normal-mode analysis on an exact scheme in an earlier paper by the author. This model is applicable to stellar models with rotation Ω=Ω(ω), where ω is the radial distance from the axis of rotation (thez-axis). A necessary condition for stability, in the non-dissipative case, is found to be that $$\frac{1}{\varrho }G_\varpi S_\varpi + \frac{{k_z^2 }}{M}\Phi - \frac{1}{4}\frac{{m^2 }}{M}\left( {D\Omega } \right)^2 \geqslant 0$$ everywhere. Here,m andk _z are the wave numbers in the ø- andz-direction, \(M \equiv k_z^2 + m^2 /\varpi ^2 ,D \equiv d/d\varpi ,G_\varpi \equiv - \varrho ^{ - 1} Dp,\varrho \) the density,p the pressure,S _ω and Φ the Schwarzschild and the Rayleigh discriminants defined as \(S_\varpi \equiv \left( {\gamma p/\varrho } \right)^{ - 2} Dp - D\varrho and \Phi \equiv ^{ - 3} d\left( {\varpi ^4 \Omega ^2 } \right)/d\varpi \) respectively, γ the ratio of specific heats. This condition is also a sufficient one. Some conjectures regarding the stabilizing influence of uniform rotation and the destabilizing influence of differential rotation are also verified. The most striking instability mechanism introduced by shear forces and by radiative dissipation is the excitation of the stable motion of small oscillations into that of oscillations with growing amplitude, i. e., overstability. In the case of radiative dissipation and axisymmetric perturbations, the Goldreich-Schubert criterion is only necessary but not sufficient for stability. Instability sets in as soon as the Schwarzschild criterion is violated. When the perturbations are non-axisymmetric, instability always sets in as overstability as long as rotation is differential. This may explain the convective turbulence in the upper atmosphere where the radiation is active.     

A few observations of and remarks on V1500 Cygni

The development of the post-nova light curve of V1500 Cyg inUBV andHβ, for 15 nights in September and October 1975 are presented. We confirm previous reports that superimposed on the steady decline of the light curve are small amplitude cyclic variations. The times of maxima and minima are determined. These together with other published values yield the following ephemerides from JD 2 442 661 to JD 2 442 674: $$\begin{gathered} {\text{From}} 17 {\text{points:}} {\text{JD}}_{ \odot \min } = 2 442 661.4881 + 0_{^. }^{\text{d}} 140 91{\text{n}} \hfill \\ \pm 0.0027 \pm 0.000 05 \hfill \\ {\text{From}} 15 {\text{points:}} {\text{JD}}_{ \odot \max } = 2 442 661.5480 + 0_{^. }^{\text{d}} 140 89{\text{n}} \hfill \\ \pm 0.0046 \pm 0.0001 \hfill \\ \end{gathered} $$ with standard errors of the fits of ±0 _. ^d 0052 for the minima and ±0 _. ^d 0091 for the maxima. Assuming V1500 Cyg is similar to novae in M31, we foundr=750 pc and a pre-nova absolute photographic magnitude greater than 9.68.     

Isovortical orbits in uniformly rotating coordinates

For the conservative, two degree-of-freedom system with autonomous potential functionV(x,y) in rotating coordinates; $$\dot u - 2n\upsilon = V_x , \dot \upsilon + 2nu = V_y $$ , vorticity (v _x -u _y ) is constant along the orbit when the relative velocity field is divergence-free such that: $$u(x,y,t) = \psi _y , \upsilon (x,y,t) = - \psi _x $$ . Unlike isoenergetic reduction using the Jacobi, integral and eliminating the time,non-singular reduction from fourth to second-order occurs when (u,v) are determined explicitly as functions of their arguments by solving for ψ (x, y, t). The orbit function ψ satisfies a second-order, non-linear partial differential equation of the Monge Ampere type: $$2(\psi _{xx} \psi _{yy} - \psi _{xy}^2 ) - 2(\psi _{xx} + \psi _{yy} ) + V_{xx} + V_{yy} = 0$$ . Isovortical orbits in the rotating frame arenot level curves of ψ because it contains time explicitly due to coriolis effects. Rather, (x, y) coordinates along the orbit are obtained, from (u, v) either by numerical integration of the kinematic equations, or by partial differentiation of the Legendre transform ? of ψ. In the latter case, ? is shown to satisfy a non-linear, second-order partial differential equation in three independent variables, derived from the Monge-Ampere Equation. Complete reduction to quadrature is possible when space-time symmetries exist, as in the case of central force motion.     

Recent progress in the theory of Trojan asteroids

In previous publications the author has constructed a long-periodic solution of the problem of the motion of the Trojan asteroids, treated as the case of 1:1 resonance in the restricted problem of three bodies. The recent progress reported here is summarized under three headings:

1. The nature on the long-periodic family of orbits is re-examined in the light of the results of the numerical integrations carried out by Deprit and Henrard (1970). In the vicinity of the critical divisor $$D_k \equiv \omega _1 - k\omega _2 ,$$ not accessible to our solution, the family is interrupted by bifurcations and shortperiodic bridges. Parametrized by the normalized Jacobi constant α², our family may, accordingly, be defined as the intersection of admissible intervals, in the form $$L = \mathop \cap \limits_j \left\{ {\left| {\alpha - \alpha _j } \right| > \varepsilon _j } \right\};j = k,k + 1, \ldots \infty .$$ Here, {α_j(m)} is the sequence of the critical α_j corresponding to the exactj: 1 commensurability between the characteristic frequencies ω₁ and ω₂ for a given value of the mass parameterm. Inasmuch as the ‘critical’ intervals |α?α_j|<ε_j can be shown to be disjoint, it follows that, despite the clustering of the sequence {α_j} at α=1, asj→∞, the family extends into the vicinity of the separatrix α=1, which terminates the ‘tadpole’ branch of the family.
2. Our analysis of the epicyclic terms of the solution, carrying the critical divisorD _k , supports the Deprit and Henrard refutation of the E. W. Brown conjecture (1911) regarding the termination of the tadpole branch at the Lagrangian pointL ₃. However, the conjecture may be revived in a refined form. “The separatrix α=1 of the tadpole branch spirals asymptotically toward a limit cycle centered onL ₃.”
3. The periodT(α,m) of the libration in the mean synodic longitude λ in the range $$\lambda _1 \leqslant \lambda \leqslant \lambda _2$$ is given by a hyperelliptic integral. This integral is formally expanded in a power series inm and α² or \(\beta \equiv \sqrt {1 - \alpha ^2 }\) .

The large amplitude of the libration, peculiar to our solution, is made possible by the mode of the expansion of the disturbing functionR. Rather than expanding about Lagrangian pointL ₄, with the coordinatesr=1, θ=π/3, we have expandedR about the circler=1. This procedure is equivalent to analytic continuation, for it replaces the circle of convergence centered atL ₄ by an annulus |r?1|<ε with 0≤θ<2π.     

Gaussian variational equations for osculating elements of an arbitrary separable reference orbit

If a satellite orbit is described by means of osculating Jacobi α's and β's of a separable problem, the paper shows that a perturbing forceF makes them vary according to $$\dot \alpha _\kappa = {\text{F}} \cdot \partial {\text{r/}}\partial \beta _k {\text{ }}\dot \beta _k = {\text{ - F}} \cdot \partial {\text{r/}}\partial \alpha _k ,{\text{ (}}k = 1,2,3).{\text{ (A1)}}$$ Herer is the position vector of the satellite andF is any perturbing force, conservative or non-conservative. There are two special cases of (A1) that have been previously derived rigorously. If the reference orbit is Keplerian, equations equivalent to (A1), withF arbitrary, were derived by Brouwer and Clemence (1961), by Danby (1962), and by Battin (1964). IfF=?gradV ₁(t), whereV ₁ may or may not depend explicitly on the time, Equations (A1) reduce to the well known forms (e.g. Garfinkel, 1966) $$\dot \alpha _\kappa = {\text{ - }}\partial V_1 {\text{/}}\partial \beta _k {\text{ }}\dot \beta _k = \partial V_1 {\text{/}}\partial \alpha _k ,{\text{ (}}k = 1,2,3).{\text{ (A2)}}$$ holding for all separable reference orbits. Equations (A1) can of course be guessed from Equations (A2), if one assumes that \(\dot \alpha _k (t)\) and \(\dot \beta _k (t)\) depend only onF(t) and thatF(t) can always be modeled instantaneously as a potential gradient. The main point of the present paper is the rigorous derivation of (A1), without resort to any such modeling procedure. Applications to the Keplerian and spheroidal reference orbits are indicated.     

Isovortical orbits of autonomous,conservative, two degree-of-freedom dynamical systems

For an autonomous, conservative, two degree-of-freedom dynamical system, vorticity (the curl of velocity) is constant along the orbit if the velocity field is divergence-free such that: $$u\left( {x, v} \right) - \psi _y , v\left( {x, y} \right) = - \psi _x .$$ Isovortical orbits in configuration space are level curves of a scalar autonomous function Ψ (x, v) satisfying a second-order, non-linear partial differential equation of the Monge-Ampere type: $$2\left( {\psi _{xx} \psi _{yy} - \psi _{xy}^2 } \right) + U_{xx} + U_{yy} = 0,$$ where U(x. y) is the autonomous potential function. The solution Soc the time variable is reduced to a quadrature following determinatio of Ψ. Self-similar solutions of the Monge-Ampere equation under Birkhoff's one-parameter transformation group are derived for homogeneous (power-law) potential functions. It is shown that Keplerian orbits belong to the class of planar isovortical flows.     

Optical spectral decomposition of NGC 4051

Considering the host galaxy contribution, a spectral decomposition method is used to reanalyzed the archive data of optical spectra for a narrow line Seyfert 1 galaxy, NGC 4051. The light curves of the continuum f _λ (5100 Å), and Hβ, He ii, Fe ii emission lines are given. We find strong flux correlations between line emissions of Hβ, He ii, Fe ii and the continuum f _λ (5100 Å). These low-ionization lines (Hβ, Fe ii, He ii) have “inverse” intrinsic Baldwin effects. Using the methods of the cross-correlation function and the Monte Carlo simulation, we find the time delays, with respect to the continuum, are $3.45^{+12.0}_{-0.5}~\mbox{days}$ with the probability of 34 % for the intermediate component of Hβ, $6.45^{+13.0}_{-1.0}~\mbox{days}$ with the probability of 65 % for the intermediate component of He ii. From these intermediate components of Hβ and He ii, the calculated central black hole masses are $0.86^{+4.35}_{-0.33}\times 10^{6}$ and $0.82^{+3.12}_{-0.45}\times 10^{6}~M_{\odot }$ . We also find that the time delays for Fe ii are $9.7^{+3.0}_{-5.0}~\mbox{days}$ with the probability of 36 %, $8.45^{+1.0}_{-2.0}~\mbox{days}$ with the probability of 18 % for the total epochs and “subset 1” data, respectively. It seems that the Fe ii emission region is outside of the Hβ emission region.     

Etude analytique du voisinage de la resonance 4:4:1 dans les systemes a trois degres de liberte

The periodic solutions for an Hamiltonian system with $$H = \frac{1}{2}\mathop \Sigma \limits_1^3 (\dot x_\alpha ^2 + \omega _\alpha ^2 x_\alpha ^2 ) - \varepsilon x_1 x_\alpha ^2 - \eta x_2 x_\alpha ^2 $$ are investigated analytically. The frequencies ω_α, α=1, 2, 3 are assumed near the ratio 4—4—1. We find different families of periodic solutions whose periods are in the vicinity of the period T′=2π/ω₃=2π/ω′. As in the case of the problem with two degrees of liberty, for particular values of ω₁, ω₂, ω₃ and ε, η, we find that the families near the x₃-axis are discontinuous. These families are periodic with periods near the period T′ in a region for ε, η, approximatively [0; 0.4] if we choose \(\omega ' = \sqrt {0.1} \) and h=0.00765.