Qualitative features of the evolution of some polar satellite orbits

Two special cases of the problem of the secular perturbations in the orbital elements of a satellite with a negligible mass produced by the joint influence of the oblateness of the central planet and the attraction by its most massive (or main) satellites and the Sun are considered. These cases are among the integrable ones in the general nonintegrable evolution problem. The first case is realized when the plane of the satellite orbit and the rotation axis of the planet lie in its orbital plane. The second case is realized when the plane of the satellite orbit is orthogonal to the line of intersection between the equatorial and orbital planes of the planet. The corresponding particular solutions correspond to those polar satellite orbits for which the main qualitative features of the evolution of the eccentricity and pericenter argument are described here. Families of integral curves have been constructed in the phase plane of these elements for the satellite systems of Jupiter, Saturn, and Uranus.     

Long-term evaluation of three-dimensional heliocentric solar sail trajectories with arbitrary fixed sail setting

The solar radiation effects upon the orbital behaviour of an arbitrarily shaped spacecraft (or a solar sail in particular) in a general fixed orientation with respect to the local coordinate frame are investigated. Through introduction of a quasi-angle in the osculating plane, the motion of the orbital plane becomes uncoupled from the in-plane perturbations. Exact solutions in the form of conic sections and logarithmic spirals can readily be formulated for certain specific initial conditions. An effective out-of-plane spiral transfer trajectory is obtained by reversing the force component normal to the orbital plane at specified positions in the orbit. By choosing the appropriate control angles for the sail orientation, any point in space can be reached eventually. In the case of general initial conditions, the long-term orbital behaviour is assessed asymptotically by means of the two-variable expansion procedure. An implicit expression for the eccentricity is derived and explicit results are established by an iteration scheme. The other orbital elements can be expressed in terms of the eccentricity and their asymptotic series for near-circular initial orbits are also obtained. While equations for the higher-order contributions as well as the periodic parts of their solutions can be formulated readily, their secular terms are determined only for a circular initial orbit.     

Motion of dust in mean motion resonances with planets

Effect of stellar electromagnetic radiation on the motion of spherical dust particle in mean motion orbital resonances with a planet is investigated. Planar circular restricted three-body problem with the Poynting–Robertson (P–R) effect yields monotonic secular evolution of eccentricity when the particle is trapped in the resonance. Planar elliptic restricted three-body problem with the P–R effect enables nonmonotonous secular evolution of eccentricity and the evolution of eccentricity is qualitatively consistent with the published results for the complicated case of interaction of electromagnetic radiation with nonspherical dust grain. Thus, it is sufficient to allow either nonzero eccentricity of the planet or nonsphericity of the grain and the orbital evolutions in the resonances are qualitatively equal for the two cases. This holds both for exterior and interior mean motion orbital resonances. Evolutions of argument of perihelion in the planar circular and elliptical restricted three-body problems are shown. Numerical integrations show that an analytic expression for the secular time derivative of the particle’s argument of perihelion does not exist, if only dependence on semimajor axis, eccentricity and argument of perihelion is admitted. Connection between the shift of perihelion and oscillations in secular eccentricity is presented for the planar elliptic restricted three-body problem with the P–R effect. Period of the oscillations corresponds to the period of one revolution of perihelion. Change of optical properties of the spherical grain with the heliocentric distance is also considered. The change of the optical properties: (i) does not have any significant influence on the secular evolution of eccentricity, (ii) causes that the shift of perihelion is mainly in the same direction/orientation as the particle motion around the Sun. The statements hold both for circular and noncircular planetary orbits.     

Gravitational Radiation Damping and Evolution of the Orbit of Compact Binary Stars (Solution by the Second Perturbation Method)

The influence of the gravitational radiation damping on the evolution of the orbital elements of compact binary stars is examined by using the method of perturbation. The perturbation equations with the true anomaly as an independent variable are given. This effect results in both the secular and periodic variation of the semi-major axis, the eccentricity, the mean longitude at the epoch and the mean longitude. However, the longitude of periastron exhibits no secular variation, but only periodic variation. The effect of secular variation of the orbit would lead to collapse of the system of binary stars. The deduced formulae are applied to the calculation of secular variation of the orbital elements for three compact binary stars: PSR19 13 + 16, PSR J0737-3039 and M33X-7. The results obtained are discussed.     

An analytical solution for Kepler's problem

In this paper, we present a framework which provides an analytical (i.e. infinitely differentiable) transformation between spatial coordinates and orbital elements for the solution of the gravitational two-body problem. The formalism omits all singular variables which otherwise would yield discontinuities. This method is based on two simple real functions for which the derivative rules are only required to be known, all other applications – e.g. calculating the orbital velocities, obtaining the partial derivatives of radial velocity curves with respect to the orbital elements – are thereafter straightforward. As it is shown, the presented formalism can be applied to find optimal instants for radial velocity measurements in transiting explanatory systems to constrain the orbital eccentricity as well as to detect secular variations in the eccentricity or in the longitude of periastron.     

太阳多方模型对行星轨道要素变化的长期摄动影响

本文利用天体力学中的摄动理论和天体物理学中的气体星多方模型理论研究了太阳多方模型对行星轨道要素变化的长期摄动影响。文中给出了太阳日多方指数ｎ＝３的模型由于自转、扁度和内部密度分布等因素对行星轨道要素变化的长期摄动效应的理论结果。研究结果表明：行星轨道要素除长轴、偏心率和轨道倾角不受长期摄动外升交点经度、近日点经度以及平近点角均受长期摄动的影响。最后利用理论结果对行星轨道要素的长期投动效应做了数值计算，数值结果在表１中给出。     

在电离中电感应阻力对带电卫星轨道根数的摄动影响

研究了在高空电离层中运动的带电荷的卫星受电感应阻力后对轨道根数产生的摄动影响。研究结果表明 ,电感应阻力对带电卫星的轨道半长轴、轨道偏心率、近地点赤经、历元平赤经均有周期摄动影响 ,但除对半长轴有长期摄动效应外对其它轨道根数均无长期摄动。轨道倾角和升交点赤经不受摄动影响。文中以飞行在高度 1 50 0km的电离层中的导体卫星作为算例。计算结果显示 :带电导体卫星在高空电离层中带有一定电量时电感应阻力对轨道半长轴的缩短产生显著效应     

Effects of solar radiation on the orbits of small particles

Revised equations of motion are formulated on more general assumptions than hitherto making allowance for some reflection of sunlight by a dust-particle, and from these the secular rates of change of the orbital elements of the particle are obtained. The equation for the eccentricity yields numerical results for the time taken for given changes in this element to occur. Other elements turn out to be expressible in terms of the eccentricity and thence are effectively also known in terms of the time. More general forms of earlier results are found, and some new mathematical results in the theory of the process are derived. The time of infall to the Sun associated with almost circular initial motion of a particle is calculated, and also the time from an orbit of initially high eccentricity. In this latter case, infall takes place much more rapidly than from a circular orbit of radius comparable with the average distance in the eccentric orbit. The effect on a particle of a long-period comet during a single return is negligible compared with the change in its binding-energy to the Sun that will in general result from planetary action. The possible history of a dust-particle from original capture by the Sun to final infall to the solar surface is briefly considered.     

On the stability of the satellites of asteroid 87 Sylvia

The triple asteroidal system (87) Sylvia is composed of a 280-km primary and two small moonlets named Romulus and Remus ( Marchis et al. 2005b ). Sylvia is located in the main asteroid belt, with semi-major axis of about 3.49 au, eccentricity of 0.08 and 11° of orbital inclination. The satellites are in nearly equatorial circular orbits around the primary, with orbital radius of about 1360 km (Romulus) and 710 km (Remus). In this work, we study the stability of the satellites Romulus and Remus. In order to identify the effects and the contribution of each perturber, we performed numerical simulations considering a set of different systems. The results from the three-body problem, Sylvia–Romulus–Remus, show no significant variation of their orbital elements. However, the inclinations of the satellites present a long-period evolution with amplitude of about 20° when the Sun is included in the system. Such amplitude is amplified to more than 50° when Jupiter is included. These evolutions are very similar for both satellites. An analysis of these results shows that Romulus and Remus are librating in a secular resonance and their longitude of the nodes are locked to each other. Further simulations show that the amplitude of oscillation of the satellites' inclination can reach higher values depending on the initial values of their longitude of pericentre. In those cases, the satellites get caught in an evection resonance with Jupiter, their eccentricities grow and they eventually collide with Sylvia. However, the orbital evolutions of the satellites became completely stable when the oblateness of Sylvia is included in the simulations. The value of Sylvia's J ₂ is about 0.17, which is very high. However, even just 0.1 per cent of this value is enough to keep the satellite's orbital elements with no significant variation.     

Influence of fast interstellar gas flow on the dynamics of dust grains

The orbital evolution of a dust particle under the action of a fast interstellar gas flow is investigated. The secular time derivatives of Keplerian orbital elements and the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle’s orbit are derived. The secular time derivatives of the semi-major axis, eccentricity, and of the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle’s orbit constitute a system of equations that determines the evolution of the particle’s orbit in space with respect to the gas flow velocity vector. This system of differential equations can be easily solved analytically. From the solution of the system we found the evolution of the Keplerian orbital elements in the special case when the orbital elements are determined with respect to a plane perpendicular to the gas flow velocity vector. Transformation of the Keplerian orbital elements determined for this special case into orbital elements determined with respect to an arbitrary oriented plane is presented. The orbital elements of the dust particle change periodically with a constant oscillation period or remain constant. Planar, perpendicular and stationary solutions are discussed. The applicability of this solution in the Solar System is also investigated. We consider icy particles with radii from 1 to 10 μm. The presented solution is valid for these particles in orbits with semi-major axes from 200 to 3000 AU and eccentricities smaller than 0.8, approximately. The oscillation periods for these orbits range from 10⁵ to 2 × 10⁶ years, approximately.     

Parameterized post-Newtonian orbital effects in extrasolar planets

Perturbative post-Newtonian variations of the standard osculating orbital elements are obtained by using the two-body equations of motion in the parameterized post-Newtonian theoretical framework. The results obtained are applied to the Einstein and Brans–Dicke theories. As a results, the semi-major axis and eccentricity exhibit periodic variation, but no secular changes. The longitude of periastron and mean longitude at epoch experience both secular and periodic shifts. The post-Newtonian effects are calculated and discussed for six extrasolar planets.     

Nonlinear theory of secular perturbations of satellites of an oblate planet

Anonlinear analytical theory of secular perturbations in the problem of the motion of a systemof small bodies around a major attractive center has been developed. Themutual perturbations of the satellites and the influence of the oblateness of the central body are taken into account in the model. In contrast to the classical Laplace-Lagrange theory based on linear equations for Lagrange elements, the third-degree terms in orbital eccentricities and inclinations are taken into account in the equations. The corresponding improvement of the solution turns out to be essential in studying the evolution of orbits over long time intervals. A program inC has been written to calculate the corrections to the fundamental frequencies of the solution and the third-degree secular perturbations in orbital eccentricities and inclinations. The proposed method has been applied to investigate the motion of the major Uranian satellites. Over time intervals longer than 100 years, allowance for the nonlinear terms in the equations is shown to give corrections to the coordinates of Miranda on the order of the orbital eccentricity, which is several thousand kilometers in linear measure. For other satellites, the effect of allowance for the nonlinear terms turns out to be smaller. Obviously, when a general analytical theory of motion for the major Uranian satellites is constructed, the nonlinear terms in the equations for the secular perturbations should be taken into account.     

Secular influence of change in the heliocentric gravitation constant <Emphasis Type="Italic">GM</Emphasis><Subscript>⊙</Subscript> on evolution of orbits of Meteor Streams

The Secular influence of the change in the heliocentric gravitational constant on the evolution of orbits of Meteor Streams is examined by using the method of celestial mechanics with variable mass and variable gravitational constant. The change in the heliocentric gravitational constant includes the combined changes in the sun’s mass and gravitational constant obtained from the modern observation of planets and spacecraft. The perturbation equations are solved by expanding series with mean anomaly. The solutions of the secular and periodic variation of orbital elements are derived. The theoretical results for the secular variables of the semi-major axes, solar distances at perihelion and orbital periods are given for three Meteor Streams: Dracorids, Quadrantids, and Ursids. The numerical results are shown in Table 2. The discussion and conclusion are drawn.     

Motion near frozen orbits as a means for mitigating satellite relative drift

Generally, any initially-close satellites—chief and deputy—moving on orbits with slightly different orbital elements, will depart each other on locally unbounded relative trajectories. Thus, constraints on the initial conditions must be imposed to mitigate the chief-deputy mutual departure. In this paper, it is analytically proven that choosing the chief’s orbit to be a frozen orbit can mitigate the natural relative drift of the satellites. Using mean orbital element variations, it is proven that if the chief’s orbit is frozen, then the mean differential eccentricity is periodic, leading to a periodic variation of the differential mean argument of latitude. On the other hand, if the chief’s orbit is non-frozen, a secular growth in the differential mean argument of latitude leads to a concomitant along-track separation of the deputy from the chief, thereby considerably increasing the relative distance evolution over time. Long-term orbital simulation results indicate that the effect of choosing a frozen orbit vis-à-vis a non-frozen orbit can reduce the relative distance drift by hundreds of meters per day.     

Improved equations for eccentricity generation in hierarchical triple systems

In a series of papers, we developed a technique for estimating the inner eccentricity in hierarchical triple systems, with the inner orbit being initially circular. However, for certain combinations of the masses and the orbital elements, the secular part of the solution failed. In this paper, we derive a new solution for the secular part of the inner eccentricity, which corrects the previous weakness. The derivation applies to hierarchical triple systems with coplanar and initially circular orbits. The new formula is tested numerically by integrating the full equations of motion for systems with mass ratios from 10⁻³ to 10³. We also present more numerical results for short-term eccentricity evolution, in order to get a better picture of the behaviour of the inner eccentricity.     

Radiation pressure and dust particle dynamics

The dynamics of small dust grains orbiting a planet are investigated when solar radiation pressure forces are added to the planet's gravitational central field. In the first part a set of differential equations is derived in a reference frame linked to the solar motion. The complete solution of these equations is given for particles lying in the planet's orbital plane, and we show that the orbital eccentricity may undergo considerable variation. At the same time the pericenter longitude librates or circulates according to initial conditions. With this result we establish a criterion for any orbiting particle (because of its highly eccentric orbit) to collide with its planet's atmosphere. The case of inclined orbit is studied through a numerical integration and allows us to draw conclusions related to the stability of the orbital plane. All solutions are periodic, with the period being independent of the initial conditions. This last point permits us to investigate the different time scales involved in that problem. Finally, the Poynting-Robertson drag is included, along with the radial radiation pressure forces, and the secular trend is considered. A coupling effect between the two components is ascertained, yielding a systematic behavior in the eccentricity and thus in the pericenter distance. Our solutions generalize the results of S. J. Peale (1966, J. Geophys. Res.71, 911–933) and J. A. Burns, P. Lamy, and S. Soter (1979, Icarus40, 1–48) by allowing eccentricities to be large (of order 1) and inclinations to be nonzero and by considering Poynting-Robertson drag.     

Long‐term orbital period behaviors of the neglected Algol type binaries: CC Herculis and XZ Aquilae

Orbital period variations of two neglected Algol type binaries, CC Her and XZ Aql, are studied based on all available times of minima. In the case of CC Her, it is found that the O –C curve displays a tilted sinusoidal variation with an eccentricity of 0.54 ± 0.03 and a period of 52.4 ± 0.4 yr, which can be explained by the light‐time effect due to the presence of an unseen component. The course of the orbital period change in XZ Aql appears less reliable but its O –C curve can be represented by a periodic variation with a period of 36.7 ± 0.6 yr superimposed on an upward parabola. The parabolic variation indicates a secular period increase with a rate of dP /dt = 7.1 s per century. The corresponding conservative mass transfer from less massive component to the more massive one is about 3.26 × 10^–7 M_⊙ yr^–1. It is interesting to see that the O –C variation of CC Her displays no evidence (as upward parabola) on the mass transfer characteristic for Algols. The periodic change of the orbital period of XZ Aql, like CC Her, may be caused by the presence of the thirdbody. The lower limits of the masses of the hypothetical unseen components for CC Her and XZ Aql are found to be 2.69 M_⊙ and 0.47 M_⊙, respectively. The third body of CC Her should be detectable not only spectroscopically but also photoelectrically, if it exists. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1-5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.     

The high amplitude delta Scuti star SZ Lyncis revisited

Nine new times of light maximum are given for SZ Lyn and the time base of observations is herewith extended to 26 years. 154 light maxima have been used to calculate the orbital elements and the secular change in the period of the pulsating component. The results provided here support the post-Main Sequence, high mass hypothesis for the evolutionary state of this dwarf cepheid.