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.     

Poynting-Robertson drag and orbital resonance

Both the Poynting-Robertson drag and resonant orbits appear to be very important for the motion of small grains in the early solar system. While orbital resonances are very often stable and tend to force bodies into noncircular orbits, the Poynting-Robertson drag produces secular variations in the semimajor axis and tends to circularize the orbits. We study numerically the competition between the Poynting-Robertson drag and the gravitational interaction of grains with Jupiter near the 2/1 resonance. Computations are based on the plane-restricted problem. Numerical investigations show that the grains always cross the resonance region without any oscillation, except in the special case where the grains were initially inside the resonance. In both cases the variations of the osculating elements exhibit a drastic step, which can be explained by Greenberg's and Schubart's theories.     

磁感应阻力对导体卫星在圆形轨道上的摄动影响

用功能转换原理和摄动理论的两种方法重点研究了导体卫星在地球磁场和有电导率介质空间飞行时磁感应阻力对圆形轨道半径的摄动影响。理论研究表明：导体卫星在圆形轨道上受磁感应阻力后轨道半径有随时间变化的长期摄动效应外还有周期性变化。此外，文中还讨论了磁感应阻力对圆形轨道的其它要素的摄动影响概况。     

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.     

Secular perturbations of asteroids in secular resonance

When the precessional rate of the orbital plane of an asteroid is nearly equal to that of Jupiter, the orbital inclination of the asteroid changes quite largely due to this near equality of their precessional rates, which is called a secular resonance. In the vicinity of the exact resonance the difference of their longitudes of nodes librates with quite a long period of order of 1×10⁶ yr. In this paper we treat this secular resonance by a method of semianalytical secular perturbations with use of numerical averaging for both non-resonant and resonant asteroids and show that the results by the semi-analytical treatment agrees qualitatively with those obtained by direct numerical integrations of asteroid's orbits.     

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.     

太阳质量损失对流星群轨道改变的长期影响

将作者在变质量天体力学所得理论结果应用于太阳质量损失对流星群轨道根数变化的长期效应上。太阳质量损失包括光子辐射和太阳风造成的质量损失。利用G—M型变质量天体轨道根数变化方程的一阶和二阶解对15个流星群轨道半长轴、近日点距离、轨道周期和近日点经度因太阳质量损失造成的每世纪的长期改变效应做了数值计算,并得出计算结果。其计算结果表明,太阳质量损失使流星群轨道半长轴每世纪的改变效应较明显,它们同太阳距离的扩大影响值得关注,但对轨道周期的拉长每世纪的影响甚小,对近日点经度只有量级变化小到可以略而不计。     

Orbital YORP and asteroid orbit evolution, with application to Apophis

Photon thrust from shape alone can produce quasi-secular changes in an asteroid's orbital elements. An asteroid in an elliptical orbit with a north–south shape asymmetry can steadily alter its elements over timescales longer than one orbital trip about the Sun. This thrust, called here orbital YORP (YORP = Yarkovsky–O'Keefe–Radzievskii–Paddack), operates even in the absence of thermal inertia, which the Yarkovsky effects require. However, unlike the Yarkovsky effects, which produce secular orbital changes over millions or billions of years, the change in an asteroid's orbital elements from orbital YORP operates only over the precession timescale of the orbit or of the asteroid's spin axis; this is generally only thousands or tens of thousands of years. Thus while the orbital YORP timescale is too short for an asteroid to secularly journey very far, it is long enough to warrant investigation with respect to 99942 Apophis, which might conceivably impact the Earth in 2036. A near-maximal orbital YORP effect is found by assuming Apophis is without thermal inertia and is shaped like a hemisphere, with its spin axis lying in the orbital plane. With these assumptions orbital YORP can change its along-track position by up to ±245 km, which is comparable to Yarkovsky effects. Though Apophis' shape, thermal properties, and spin axis orientation are currently unknown, the practical upper and lower limits are liable to be much less than the ±245 km extremes. Even so, the uncertainty in position is still likely to be much larger than the 0.5 km “keyhole” Apophis must pass through during its close approach in 2029 in order to strike the Earth in 2036.     

Asteroids with large secular orbital variations

As the classical linear theory of secular perturbations for asteroids is known not to be adequate for computing the perturbations of asteroids with high eccentricities and/or inclinations, a seminumerical method to calculate the secular perturbations by including higher-degree terms in the disturbing function has been developed. It is here applied to asteroids with small values of $\text{(1 ? e}{}^2\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{cos}\text{i}$, since the secular variations as well as their deviations from the results derived by the classical linear theory are generally large for such asteroids. It is found that the arguments of perihelion for five of the numbered asteroids are librating around 90 or 270°. For asteroids with $\text{(1 ? e}{}^2\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext><mtext>cos</mtext><mtext>\; i</mtext>}}$ less than 0.85 the results of the secular variations are tabulated. Also the stability of such orbits is discussed by comparing the orbital properties of short-periodic comets with them. Generally speaking, orbits of the asteroids are more stable than those of the short-periodic comets, and asteroids with librating arguments of perihelion are more stable than those with circular coplanar orbits although their orbital elements are changed more by secular perturbations.     

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.     

Orbital resonance in a dissipative medium

Orbital resonances tend to force bodies into noncircular orbits. If a body is also under the influence of an eccentricity-reducing medium, it will experience a secular change in semimajor axis which may be positive or negative depending on whether its orbit is exterior or interior to that of the perturbing body. Thus a dissipative medium can promote either a loss or a gain in orbital energy. This process may explain the resonant structure of the asteroid belt and of Saturn's rings. For reasonable early solar system parameters, it would clear a gap near the 2:1 resonance with Jupiter on a time scale of a few thousand years; the gap width would be comparable to the Kirkwood gap presently at the location in the asteroid belt. Similarly, a gap comparable in width to Cassini's division would be cleared in Saturn's rings at the 2:1 resonance with Mimas in ～10⁶ yr. Most of the material from the gap would be deposited at the outer edge of ring B. The process would also affect the radial distribution of preplanetary material. Moreover, it provides an explanation for the large amplitude of the Titan-Hyperion libration. Consideration of the effects of dissipation on orbits near the stable L₄ and L₅ points of the restricted three-body problem indicates that energy loss causes particles to move away from these points. This results explains the large amplitude of Trojan asteroids about these points and the possible capture of Trojan into orbit about Jupiter.     

Nodal regression of the quadrantid meteor stream: An analytic approach

The mean orbit of the Quadrantid meteor stream has a high eccentricity and inclination with an aphelion close to the orbit of Jupiter. The nodal regression rate, a quantity which has been well determined from observations, cannot be calculated with sufficient accuracy using standard low-order expansions of the disturbing function. By using a high-order expansion of the disturbing function we show how the behavior of the longitude of ascending node of the Quadrantid stream is a result of both secular and resonant effects. Our analysis illustrates how the proximity of the stream's orbit to the 2: 1 commensurability with Jupiter dominates the short-term variations in orbital elements.     

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.     

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

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

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.     

Periodic motions around pulsating stars

The periodic motion of a test particle (dust, grain, or a larger body) around a pulsating star with a luminosity oscillation of small amplitude (featured by a small parameterB) is being studied. The perturbations of all orbital elements are determined to first order inB, by using Delaunay-type canonical variables and a method whose bases were put forth by von Zeipel. According to the value of the ratio oscillation frequency/dynamic frequency, three possible situations are pointed out: nonresonant (NR), quasi-resonant (QR), and resonant (R). The solution of motion equations shows that only in the (QR) and (R) cases there are orbital parameters (argument of periastron and mean anomaly) affected by secular perturbations. These solutions (which indicate a secularly stable motion in a first approximation) are valid over prediction times of orderB ^–1 in the (NR) case andB ^–1/2 in the (QR) and (R) cases. The theory may be applied to various astronomical situations.     

Io's sodium emission cloud

Strong evidence that Io's sodium emission is due to resonant scattering is given by our observations which show a monotonic increase of emission intensity with residual solar intensity. In addition we detected no emission during three eclipse observations of Io. We propose a resonant scattering model with two spacial components comprising an optically thick atmosphere extending 10³ km above Io's surface surrounded by an optically thin cloud which forms a partial torus around Jupiter. In this model a flux of 10⁷ cm^?2 sec^?1 sodium atoms are sputtered from Io's surface by heavy energetic ions which are accelerated in a plasma sheath around Io. The atoms sputtered from the surface collide with atoms in Io's atmosphere so the equipartition of kinetic energy is established. The total sodium abundance is about 3 × 10¹³ cm^?2. During Io's day, sodium and other atmospheric constituents are ionized, giving rise to the ionosphere observed by Pioneer 10. Atoms escape by means of Jeans escape from the critical level, which is at the top of the atmosphere and the base of the cloud. We have observed sodium emission 6arcsec (6 Io diameters) above and below Io's orbital plane and 23arcsec toward Jupiter in Io's orbital plane. No emission was detected at maximum elongation 180° from Io. We interpret these results to mean that atoms escaping from Io form a partial torus whose thickness is about 12 arcsec and whose length is at least one-fifth of Io's orbital circumference.     

Tidal Evolution of the Kepler Candidate Two-planet Systems

The tidal evolution of ten Kepler candidate two-planet systems is investigated by using the general secular perturbation theory, and then a general picture of tidal evolution for these systems is described. Taking the KOI 1239 system as an example, the tidal effect of the system is studied in detail, the results indicate that the dissipative term of the tidal effect causes the attenuation of planets’ orbital eccentricities, and it plays a dominant role in the process of orbital evolution, however, the conservative term of the tidal effect and the relativistic effect may damp the excitation of the eccentricity of the inner planet under the secular perturbation of the outer planet. In addition, the process of tidal evolution is also affected by both the initial eccentricity of the outer planet and the planet's tidal dissipation coeffcient. At the same time, the numerical simulation on the tidal evolution of the KOI 1239 system is also made, and the numerical results are consistent with those of general secular perturbation theory.     

Disc–Planet Interactions: Migration and Resonances in Extrasolar Planetary Systems

We consider orbital resonances in multiplanet systems. These are expected to arise during or just after formation in a gaseous disc. Disc–planet interaction naturally produces orbital migration and circularization through the action of tidal torques which in turn may lead to an orbital resonance. The mass and angular momentum content of the disc is likely to be comparable to that in the planets so that it is essential to fully incorporate the disc in the analysis.We study the orbital evolution of two planets locked in 2:1 commensurability through migration tidally induced by the disc using both analytic methods and numerical hydrodynamic simulations. The planets are assumed to orbit in an inner cavity containing at most only a small amount of disc material. Results are found to be sensitive to initial surface density profile, planet masses and disc parameters. The evolution may range between attaining and subsequently maintaining a resonance lock with two angles librating to divergent migration with no commensurability formed. In the former case eccentricities increase monotonically with time while the system undergoes inward migration. If the migration is halted by loss of the disc leaving the planets in a final configuration, there is likely to be a low probability of seeing resonant planets at small radii as well as a sensitive dependence on past history.We have also considered a multiplanet system in secular apsidal resonance. We consider the system as being in just one secular normal mode and include the effects of a gaseous disc. It is suggested that a normal mode may be selected by adding in some weak dissipative process in the disc and that it may remain, involving only the planets, when the disc is slowly removed.     

Orbital period changes of contact binary systems: direct evidence for thermal relaxation oscillation theory

Orbital period changes of ten contact binary systems (S Ant, ε CrA, EF Dra, UZ Leo, XZ Leo, TY Men, V566 Oph, TY Pup, RZ Tau and AG Vir) are studied based on the analysis of their     curves. It is discovered that the periods of the six systems, S Ant, ε CrA, EF Dra, XZ Leo, TY Men and TY Pup, show secular increases. For UZ Leo, its secular period increase rate is revised. For the three systems, V566 Oph, RZ Tau and AG Vir, weak evidence is presented that a periodic oscillation (with periods of 20.4, 28.5 and 40.9 yr respectively) is superimposed on a secular period increase. The cyclic period changes can be explained by the presence of an unseen third body in the three systems. All the sample stars studied are contact binaries with     .
Furthermore, orbital period changes of 27 hot contact binaries have been checked. It is found that, apart from AW UMa with the lowest mass ratio     , none shows an orbital period decrease. The relatively weak magnetic activity in the hotter contact binaries means little angular momentum loss (AML) from the systems via magnetic stellar winds. The period increases of these W UMa binaries can be explained by mass transfer from the secondary to the primary components, which is in agreement with the prediction of the thermal relaxation oscillation (TRO) models. This suggests that the evolution of a hotter W UMa star is mainly controlled by TRO. On the other hand, for a cooler W UMa star     , its evolution may be TRO plus AML, which coincides with the recent results of Qian.