Slow modes in stellar systems with nearly harmonic potentials: I. Spoke approximation,radial orbit instability

Using a consistent perturbation theory for collisionless disk-like and spherical star clusters, we construct a theory of slow modes for systems having an extended central region with a nearly harmonic potential due to the presence of a fairly homogeneous (on the scales of the stellar system) heavy, dynamically passive halo. In such systems, the stellar orbits are slowly precessing, centrally symmetric ellipses (2: 1 orbits). Depending on the density distribution in the system and the degree of halo inhomogeneity, the orbit precession can be both prograde and retrograde, in contrast to systems with 1: 1 elliptical orbits where the precession is unequivocally retrograde. In the first paper, we show that in the case where at least some of the orbits have a prograde precession and the stellar distribution function is a decreasing function of angular momentum, an instability that turns into the well-known radial orbit instability in the limit of low angular momenta can develop in the system. We also explore the question of whether the so-called spoke approximation, a simplified version of the slow mode approximation, is applicable for investigating the instability of stellar systems with highly elongated orbits. Highly elongated orbits in clusters with nonsingular gravitational potentials are known to be also slowly precessing 2: 1 ellipses. This explains the attempts to use the spoke approximation in finding the spectrum of slow modes with frequencies of the order of the orbit precession rate. We show that, in contrast to the previously accepted view, the dependence of the precession rate on angular momentum can differ significantly from a linear one even in a narrow range of variation of the distribution function in angular momentum. Nevertheless, using a proper precession curve in the spoke approximation allows us to partially “rehabilitate” the spoke approach, i.e., to correctly determine the instability growth rate, at least in the principal (O(α _T^−1/2) order of the perturbation theory in dimensionless small parameter α _T, which characterizes the width of the distribution function in angular momentum near radial orbits.     

Simulations of slow bars in anisotropic disk systems

The instability of anisotropic disk systems with elongated stellar orbits has been investigated. N-body generalized polytropic models of stellar disks have been constructed. They are shown to be unstable with respect to the bar formation at any degree of anisotropy. This result differs from the results of the studies of such models by other authors. The bar pattern speed and amplitude have been found. The initial distribution of precession rates and the adiabatic invariants of stellar orbits have been calculated. A bar is shown to be formed in such systems due to the radial orbit instability.     

A unified theory of the formation of galactic bars

We give arguments for a basically unified formation mechanism of slow (Lynden-Bell) and fast (common) galactic bars. This mechanism is based on an instability that is akin to the well-known instability of radial orbits and is produced by the mutual attraction and alignment of precessing stellar orbits (so far, only the formation of slow bars has been explained in this way). We present a general theory of the low-frequency modes in a disk that consists of orbits precessing at different angular velocities. The problem of determining these modes is reduced to integral equations of moderately complex structure. The characteristic pattern angular velocities Ω_p of the low-frequency modes are of the order of the mean orbital precession angular velocity \(\bar \Omega _{pr}\). Bar modes are also among the low-frequency modes; while \(\Omega _p \approx \bar \Omega _{pr}\) for slow bars, Ω_p for fast bars can appreciably exceed even the maximum orbital precession angular velocity in the disk Ω _pr ^max (however, it remains of the order of these precession angular velocities). The possibility of such an excess of Ω_p over Ω _pr ^max is associated with the effect of “repelling” orbits. The latter tend to move in a direction opposite to the direction in which they are pushed. We analyze the pattern of orbital precession in potentials typical of galactic disks. We note that the maximum radius of an “attracting” circular orbit r_c can serve as a reasonable estimate of the bar length l_b. Such an estimate is in good agreement with the available results of N-body simulations.     

Precession of the Orbit of a Planet around Stars Revolving along: Low-Eccentricity Orbits in a Binary System: Analytical Solution

In our previous paper (hereafter, paper I) we presented analytical results on the non-planar motion of a planet around a binary star for the cases of the circular orbits of the components of the binary. We found that the orbital plane of the planet (the plane containing the “unperturbed” elliptical orbit of the planet), in addition to precessing about the angular momentum of the binary, undergoes simultaneously the precession within the orbital plane. We demonstrated that the analytically calculated frequency of this additional precession is not the same as the frequency of the precession of the orbital plane about the angular momentum of the binary, though the frequencies of both precessions are of the same order of magnitude. In the present paper we extend the analytical results from paper I by relaxing the assumption that the binary is circular – by allowing for a relatively small eccentricity ε of the stars orbits in the binary. We obtain an additional, ε-dependent term in the effective potential for the motion of the planet. By analytical calculations we demonstrate that in the particular case of the planar geometry (where the planetary orbit is in the plane of the stars orbits), it leads to an additional contribution to the frequency of the precession of the planetary orbit. We show that this additional, ε-dependent contribution to the precession frequency of the planetary orbit can reach the same order of magnitude as the primary, ε-independent contribution to the precession frequency. Besides, we also obtain analytical results for another type of the non-planar configuration corresponding to the linear oscillatory motion of the planet along the axis of the symmetry of the circular orbits of the stars. We show that as the absolute value of the energy increases, the period of the oscillations decreases.     

Relativistic Precession of Elliptical Orbits of a Circumbinary Planet can be Cancelled

In publications presenting analytical results on the non-coplanar motion of a circumbinary planet it was shown that the unperturbed elliptical orbit of the planet undergoes simultaneously two kinds of the precession: the precession of the orbital plane and the precession of the orbit in its own plane. It is also well-known that there is also the relativistic precession of the planetary orbit in its own plane. In the present paper we study a combined effect of the all of the above precessions. For the general case, where the planetary orbit is not coplanar with the stars orbits, we analyzed the dependence of the critical inclination angle i_c, at which the precession of the planetary orbit in its own plane vanishes, on the angular momentum L of the planet. We showed that the larger the angular momentum, the smaller the critical inclination angle becomes. We presented the analytical result for i_c(L) and calculated the value of L, for which the critical inclination value becomes zero. For the particular case, where the planetary orbit is not coplanar with the stars orbits, we demonstrated analytically that at a certain value of the angular momentum of the planet, the elliptical orbit of the planet would become stationary: no precession. In other words, at this value of the angular momentum, the relativistic precession of the planetary orbit and its precession, caused by the fact that the planet revolves around a binary (rather than single) star, cancel each other out. This is a counterintuitive result.     

The effect of orbital eccentricity and nodal regression on spin-stabilized satellites

The gyroscopic motion of a spin-stabilized satellite due to gravity gradient torques in a circular orbit has been analyzed to varying degrees in numerous publications. This paper shows that the restriction to a circular orbit is, in fact, not essential and with a slight increase in complexity, noncircular orbits may be treated. More importantly, a uniform regression of the orbital node can also be accounted for. The general results are expressed in closed form using Jacobian elliptic functions. Finally, and this is perhaps most important, certain algebraic integrals of the precession are given which can be used to place limits on the excursions of the spin axis without actually solving for the motion. This allows one to design orientations such that the maximum angle between the orbit normal and spin axis never exceeds a specific amount even though the orbit normal is in precession.     

Free and forced obliquities of the Galilean satellites of Jupiter

The obliquity, or angular separation between orbit normal and spin pole, is an important parameter for the geodynamics of most Solar System bodies. Tidal dissipation has driven the obliquities of the Galilean satellites of Jupiter to small, but non-zero values. We present estimates of the free and forced obliquities of these satellites using a simple secular variation model for the orbits, and spin pole precession rate estimates based on gravity field parameters derived from Galileo spacecraft encounters. The free obliquity values are not well constrained by observations, but are presumed to be very small. The forced obliquity variations depend only on the orbital variations and the spin pole precession rate parameters, which are quite well known. These variations are large enough to influence spatial and temporal patterns of tidal dissipation and tidal stress.     

空间碎片轨道预报中动力学模型基准的选取

在不同的轨道预报场景中, 使用的动力学模型也不同. 例如, 在低轨空间碎片的预报中大气阻力是十分重要的摄动力, 而到了中高轨, 大气阻力就可以忽略不计. 如何为不同轨道类型的空间碎片选择最优(满足精度要求下的最简)动力学模型还没有系统、详尽的研究. 将对不同精度需求、不同轨道类型下的大批量轨道进行预报, 通过比较不同动力学模型下的预报结果, 给出各种预报场景的最优动力学模型建议. 可以为不同轨道类型的空间碎片在轨道预报时选择基准动力学模型提供参考或标准.     

Gravitational loss-cone instability in stellar systems with retrograde orbit precession

We study spherical and disc clusters in a near-Keplerian potential of galactic centres or massive black holes. In such a potential orbit precession is commonly retrograde, that is, the direction of the orbit precession is opposite to the orbital motion. It is assumed that stellar systems consist of nearly-radial orbits. We show that if there is a loss-cone at low angular momentum (e.g. due to consumption of stars by a black hole), an instability similar to loss-cone instability in plasma may occur. The gravitational loss-cone instability is expected to enhance black hole feeding rates. For spherical systems, the instability is possible for the number of spherical harmonics   l ≥ 3  . If there is some amount of counter-rotating stars in flattened systems, they generally exhibit the instability independent of azimuthal number m . The results are compared with those obtained recently by Tremaine for distribution functions monotonically increasing with angular momentum.
The analysis is based on simple characteristic equations describing small perturbations in a disc or a sphere of stellar orbits highly elongated in radius. These characteristic equations are derived from the linearized Vlasov equations (combining the collisionless Boltzmann kinetic equation and the Poisson equation), using the action-angle variables. We use two techniques for analysing the characteristic equations: the first one is based on preliminary finding of neutral modes, and the second one employs a counterpart of the plasma Penrose–Nyquist criterion for disc and spherical gravitational systems.     

Alignment time-scale of the microquasar GRO J1655−40

The microquasar GRO J1655−40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady-state disc warped by Lense–Thirring precession and find the time-scale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically 50 times higher in the latter case but we find that, in both the cases, the lifetime of the mass-transfer state is at most a few times the alignment time-scale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter aligned after its supernova kick but that it is most likely to be close to alignment rather than counter alignment now.     

Periodic orbits in the general three-body problem

Stability regions are identified in the neighborhood of periodic orbits. Features of motion in these regions are investigated. The structure of stability regions in the neighborhood of the Schubart, Moore, and Broucke orbits, the S-orbit, and the Ducati orbit is studied. The following features of motion are identified near these periodic orbits: libration, precession, symmetrization, centralization, bounce (a transition between types of trajectories), ejections, etc.     

Lagrange’s planetary equations for the motion of electrostatically charged spacecraft

A spacecraft that generates an electrostatic charge on its surface in a planetary magnetic field will be subject to a perturbative Lorentz force. Active modulation of the surface charge can take advantage of this electromagnetic perturbation to modify or to do work on the spacecraft’s orbit. Lagrange’s planetary equations are derived using the Lorentz force as the perturbation on a Keplerian orbit, incorporating orbital inclination and true anomaly for the first time for an electrostatically charged vehicle. The planetary equations reveal that orbital inclination is a second-order effect on the perturbation, explaining results found in earlier studies through numerical integration. All of the orbital elements are coupled, but the coupling notably does not depend on the magnitude of the electrostatic charge or on the strength of the magnetic field. Analytical expressions that characterize this coupling are tested with a propellantless escape example at Jupiter. A closed-form solution exists that constrains the set of equatorial orbits for which planetary escape is possible, and a sufficient condition is identified for escape from inclined orbits. The analytical solutions agree with results from the numerically integrated equations of motion to within a fraction of a percent.     

Regular and irregular periodic orbits

We distinguish between regular orbits, that bifurcate from the main families of periodic orbits (those that exist also in the unperturbed case) and irregular periodic orbits, that are independent of the above. The genuine irregular families cannot be made to join the regular families by changing some parameters. We present evidence that all irregular families appear inside lobes formed by the asymptotic curves of the unstable periodic orbits. We study in particular a dynamical system of two degrees of freedom, that is symmetric with respect to the x-axis, and has also a triple resonance in its unperturbed form. The distribution of the periodic orbits (points on a Poincaré surface of section) shows some conspicuous lines composed of points of different multiplicities. The regular periodic orbits along these lines belong to Farey trees. But there are also lines composed mainly of irregular orbits. These are images of the x-axis in the map defined on the Poincaré surface of section. Higher order iterations of this map , close to the unstable triple periodic orbit, produce lines that are close to the asymptotic curves of this unstable orbit. The homoclinic tangle, formed by these asymptotic curves, contains many regular orbits, that were generated by bifurcation from the central orbit, but were trapped inside the tangle as the perturbation increased. We found some stable periodic orbits inside the homoclinic tangle, both regular and irregular. This proves that the homoclinic tangle is not completely chaotic, but contains gaps (islands of stability) filled with KAM curves.     

Secular and tidal evolution of circumbinary systems

We investigate the secular dynamics of three-body circumbinary systems under the effect of tides. We use the octupolar non-restricted approximation for the orbital interactions, general relativity corrections, the quadrupolar approximation for the spins, and the viscous linear model for tides. We derive the averaged equations of motion in a simplified vectorial formalism, which is suitable to model the long-term evolution of a wide variety of circumbinary systems in very eccentric and inclined orbits. In particular, this vectorial approach can be used to derive constraints for tidal migration, capture in Cassini states, and stellar spin–orbit misalignment. We show that circumbinary planets with initial arbitrary orbital inclination can become coplanar through a secular resonance between the precession of the orbit and the precession of the spin of one of the stars. We also show that circumbinary systems for which the pericenter of the inner orbit is initially in libration present chaotic motion for the spins and for the eccentricity of the outer orbit. Because our model is valid for the non-restricted problem, it can also be applied to any three-body hierarchical system such as star–planet–satellite systems and triple stellar systems.     

Regular motions in double bars – I. Double-frequency orbits and loops

Bars in galaxies are mainly supported by particles trapped around stable periodic orbits. These orbits represent oscillatory motion with only one frequency, which is the bar driving frequency, and miss free oscillations. We show that a similar situation takes place in double bars: particles get trapped around parent orbits, which in this case represent oscillatory motion with two frequencies of driving by the two bars, and which also lack free oscillations. Thus the parent orbits, which constitute the backbone of an oscillating potential of two independently rotating bars, are the double-frequency orbits. These orbits do not close in any reference frame, but they map on to closed curves called loops. Trajectories trapped around the parent double-frequency orbit map on to a set of points confined within a ring surrounding the loop.     

High-precision repeat-groundtrack orbit design and maintenance for Earth observation missions

The focus of this paper is the design and station keeping of repeat-groundtrack orbits for Sun-synchronous satellites. A method to compute the semimajor axis of the orbit is presented together with a station-keeping strategy to compensate for the perturbation due to the atmospheric drag. The results show that the nodal period converges gradually with the increase of the order used in the zonal perturbations up to \(J_{15}\). A differential correction algorithm is performed to obtain the nominal semimajor axis of the reference orbit from the inputs of the desired nodal period, eccentricity, inclination and argument of perigee. To keep the satellite in the proximity of the repeat-groundtrack condition, a practical orbit maintenance strategy is proposed in the presence of errors in the orbital measurements and control, as well as in the estimation of the semimajor axis decay rate. The performance of the maintenance strategy is assessed via the Monte Carlo simulation and the validation in a high fidelity model. Numerical simulations substantiate the validity of proposed mean-elements-based orbit maintenance strategy for repeat-groundtrack orbits.     

The 1/1 resonance in extrasolar planetary systems

We study orbits of planetary systems with two planets, for planar motion, at the 1/1 resonance. This means that the semimajor axes of the two planets are almost equal, but the eccentricities and the position of each planet on its orbit, at a certain epoch, take different values. We consider the general case of different planetary masses and, as a special case, we consider equal planetary masses. We start with the exact resonance, which we define as the 1/1 resonant periodic motion, in a rotating frame, and study the topology of the phase space and the long term evolution of the system in the vicinity of the exact resonance, by rotating the orbit of the outer planet, which implies that the resonance and the eccentricities are not affected, but the symmetry is destroyed. There exist, for each mass ratio of the planets, two families of symmetric periodic orbits, which differ in phase only. One is stable and the other is unstable. In the stable family the planetary orbits are in antialignment and in the unstable family the planetary orbits are in alignment. Along the stable resonant family there is a smooth transition from planetary orbits of the two planets, revolving around the Sun in eccentric orbits, to a close binary of the two planets, whose center of mass revolves around the Sun. Along the unstable family we start with a collinear Euler–Moulton central configuration solution and end to a planetary system where one planet has a circular orbit and the other a Keplerian rectilinear orbit, with unit eccentricity. It is conjectured that due to a migration process it could be possible to start with a 1/1 resonant periodic orbit of the planetary type and end up to a satellite-type orbit, or vice versa, moving along the stable family of periodic orbits.