Poynting-Robertson effect: ‘Circular’ orbit

Time evolution of the interplanetary dust particle under the action of the solar electromagnetic radiation (Poynting-Robertson effect) is investigated. Evolution of the initially circular orbit in terms of the orbital elements present in the standard equations for their secular changes is considered. It is pointed out that the osculating eccentricity is practically constant during the motion in spite of generally accepted opinion that the standard equations for the secular changes of orbital elements represent time evolution of the osculating elements.     

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.     

Perihelion motion of small irregular dust particles due to radiation forces

Perihelion motion, i.e. a secular change of longitude of perihelion, of interplanetary dust particles is investigated under the action of solar gravity and solar electromagnetic radiation. As for spherical particle [Kla?ka, J., 2004. Electromagnetic radiation and motion of a particle. Cel. Mech. Dynam. Astron. 89, 1-61]: (i) perihelion motion is of the order ( is heliocentric velocity of the meteoroid and c is the speed of light in vacuum), if a component of electromagnetic radiation acceleration is considered as a part of central acceleration; (ii) perihelion motion is of the first order in if the total electromagnetic radiation force is considered as a disturbing force. The new facts presented in this paper concern irregular dust particles. Detailed numerical calculations were performed for the grains ejected at aphelion of comet Encke. Perihelion motion for irregular interplanetary dust particles exists already in the first order in for both cases of central accelerations. Moreover, perihelion motion of irregular particles exhibits both positive and negative directions during the particle orbital motion. Irregularity of the grains causes not only perihelion motion, but also dispersion of the dust in various directions, also normal to the orbital plane of the parent body.     

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.     

Perturbation equations of celestial mechanics

Perturbation equations of celestial mechanics in terms of orbital elements are completely derived in application to the motion of interplanetary dust particle in the gravational field of the Sun and under the action of disturbing forces. Consideration of change of mass of interplanetary dust particle is the most important feature of this derivation. The results obtained are completely general in the case of constant masses.     

Temperature-influenced dynamics of small dust particles

The motion of spherical dust particles under the action of gravity, electromagnetic radiation force and Lorentz force (LF) is studied theoretically for materials with temperature-dependent dielectric functions in the visible (VIS) spectral range. Even a weak variation of the optical constants with heliocentric distance may influence predominately a long-term dynamical behaviour of submicron-sized and small micron-sized dust grains. It is shown that the lifetime of carbonaceous or Si particles may change by several tens of per cent because of the temperature dependence of particle refractive indices. The orbital inclination is the most evident difference between the evolution of a dust particle with temperature-dependent optical properties and one without. While carbonaceous 2-μm-sized particles with optical constants independent of temperature may evolve in orbits with inclinations greater than an initial value, grains of the same size with variable refractive indices will be spread along orbits characterized with inclinations lower than the initial one. Here the temperature works as a separation factor for particles having slightly different temperature dependences of the optical constants.     

The Effect of Direct Solar Radiation Pressure on a Spacecraft of Complex Shape

An analytical solution for the joint effects of the Earth oblateness and the direct solar radiation pressure on the motion of an Artificial Earth Satellite of complex shape is constructed. The equations of motion are derived in the previous paper (hereafter refered to as paper I). The solution is effected through two canonical transformations retaining secular and periodic terms up to orders 3 and 2 respectively. The developments stressed on the effects of the radiation pressure and its coupling with the earth's gravity. A procedure for the computation of position and velocity is outlined. The conditions of the resonance are determined and the procedure for the transformations in the case of resonance is outlined. The solution revealed as expected that radiation pressure produced secular effects at the third order resulting from the coupling between periodic terms at lower orders. These affect both the main satellite body and the antenna. This revised version was published online in July 2006 with corrections to the Cover Date.     

Perturbations in close binary systems produced by the tides lagging in latitude

The aim of this investigation is to present the periodic and secular perturbations of the orbital elements of close binary systems due to tidal lag in latitude. The variational equations of the problem of plane motion will be set up in terms of the rectengular componentsR, S, andW of the disturbing accelerations. These equations are highly nonlinear with respect to the orbital elements and we present analytic approximations to the effects produced by the perturbing acceleration due to dynamical tides lagging in latitude. The perturbed elements of the orbit have been expressed by means of Hansen coefficients in the compact form of summations.     

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.     

On the Milankovitch orbital elements for perturbed Keplerian motion

We consider sets of natural vectorial orbital elements of the Milankovitch type for perturbed Keplerian motion. These elements are closely related to the two vectorial first integrals of the unperturbed two-body problem; namely, the angular momentum vector and the Laplace–Runge–Lenz vector. After a detailed historical discussion of the origin and development of such elements, nonsingular equations for the time variations of these sets of elements under perturbations are established, both in Lagrangian and Gaussian form. After averaging, a compact, elegant, and symmetrical form of secular Milankovitch-like equations is obtained, which reminds of the structure of canonical systems of equations in Hamiltonian mechanics. As an application of this vectorial formulation, we analyze the motion of an object orbiting about a planet (idealized as a point mass moving in a heliocentric elliptical orbit) and subject to solar radiation pressure acceleration (obeying an inverse-square law). We show that the corresponding secular problem is integrable and we give an explicit closed-form solution.     

Dynamics of close binary systems with mass exchange

The effects of non-isotropic ejection of mass from either component of a binary system on the orbital elements are studied, for the case of a small initial eccentricity of the relative orbit, when all the ejected mass falls on the other component. The problem is transformed to an equivalent two-body problem with isotropic variation of mass, plus a perturbing force which is a function of the intial conditions of ejection of the particles and their final, positions and velocities when they fall on the surface of the other star. The variation of the orbital elements are derived. It is shown that, to first-order terms in the eccentricity, the secular change of the semimajor axis is equal to the one corresponding to the case of zero initial eccentricity. On the contrary, the secular change of the eccentricity is smaller and it depends on the variations of mass ejection due to the finite eccentricity.     

Radiation Pressure and the Poynting-Robertson Effect for Fluffy Dust Particles

A correct understanding of the dynamical effect of solar radiation exerted on fluffy dust particles can be achieved with assistance of a light scattering theory as well as the equation of motion. We reformulate the equation of motion so that the radiation pressure and the Poynting-Robertson effect on fluffy grains are given in both radial and nonradial directions from the center of the Sun. This allows numerical estimates of these radiation forces on fluffy dust aggregates in the framework of the discrete dipole approximation, in which the first term of the scattering coefficients in Mie theory determines the polarizability of homogeneous spheres forming the aggregates.The nonsphericity in shape turns out to play a key role in the dynamical evolution of dust particles, while its consequence depends on the rotation rate and axis of the grains. Unless a fluffy dust particle rapidly revolves on its randomly oriented axis, the nonradial radiation forces may prevent, apart from the orbital eccentricity and semimajor axis, the orbital inclination of the particle from being preserved in orbit around the Sun. However, a change in the inclination is most probably controlled by the Lorentz force as a consequence of the interaction between electric charges on the grains and the solar magnetic field. Although rapidly and randomly rotating grains spiral into the Sun under the Poynting-Robertson effect in spite of their shapes and structures, fluffy grains drift inward on time scales longer at submicrometer sizes and shorter at much larger sizes than spherical grains of the same sizes. Numerical calculations reveal that the dynamical lifetimes of fluffy particles are determined by the material composition of the grains rather than by their morphological structures and sizes. The Poynting-Robertson effect alone is nevertheless insufficient for giving a satisfactory estimate of lifetimes for fluffy dust grains since their large ratios of cross section to mass would reduce the lifetimes by enhancing the collisional probabilities. We also show that the radiation pressure on a dust particle varies with the orbital velocity of the particle but that this effect is negligibly small for dust grains in the Solar System.     

Interplanetary dust particles and solar radiation

The problem of the action of the solar radiation on the motion of interplanetary dust particle is discussed. Differences between the action of electromagnetic solar radiation and that of the solar wind are explained not only from the point of view of the physical nature of these phenomena but also from the point of view of dust particle's orbital evolution. As for the electromagnetic solar radiation, general equation of motion for the particle is written and the most important consequences are: (i) the process of inspiralling toward the Sun is not the only possible motion - even spiralling from the Sun is also possible, and, (ii) the orbital plane of the particle (its inclination) may change in time. As for the solar wind, the effect corresponding to the fact that solar wind particles spread out from the Sun in nonradial direction causes that the process of inspiralling toward the Sun is in more than 50% less effective than for radial spread out; in the region of the asteroid belt (long period orbits) the process of inspiralling is changed into offspiralling. Also shift in the perihelion of dust particle's orbit exists.     

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.     

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.     

Poynting-Robertson effect and orbital motion

Time evolution of the meteoroid's orbit under the action of the solar electromagnetic radiation is discussed in terms of perihelion and aphelion distances. Perturbation equations for secular changes of orbital elements are written for the most simple case. Initial conditions are formulated for the obtained system of perturbation equations and simple example is presented.     

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.     

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.     

Perturbative effects on a Mercurian orbiter due to the solar radiation pressure,solar wind and the coronal mass ejections

To explore the dynamics of a test particle in the near-Mercury’s environment, the orbital motion of an orbiter around Mercury is considered. Different perturbing forces, namely the Mercurian gravity field, the solar radiation pressure, the solar wind and the coronal mass ejections, are taken into account. The order of magnitude of each perturbing term is assessed. The equations of motion in canonical representation are obtained. The Hamiltonian in terms of Hansen coefficients is expressed. A procedure for solution is presented. The short and long periodic terms are removed from the Hamiltonian and the solution is obtained. Long periodic perturbations on the orbital dynamics of an orbiter around Mercury due to the solar events are found as revealed by Eq. (26) in the text. Resonance cases are discussed and the different resonant inclinations are obtained. A procedure for the computation of the position and velocity is presented.     

Motion of dust particles under the influence of a latitudinally dependent force

The Kelperian motion of dust particles in the solar system is mainly influenced by the electromagnetic and plasma Poynting-Robertson drag. The first force is isotropic while the second one shows latitudinal variations due to the observed differences of the solar wind parameters in the ecliptic plane and over the solar poles. Close to the Sun other effects become important, e.g. sublimation and sputtering, as well as for submicron particles Lorentz scattering has to be taken into account. These forces are very weak for dust grains of moderate size (10–100 µ) not too close (>0.03 AU) to the Sun and are neglected here. Assuming that the general form of the latidudinally dependent force is a series expansion in Legendre polynomials, we have studied the averaged equations of motion for the classical elements and found the first integral of them. The general character of motion is the same as for the classical Poynting-Robertson drag: particles spiral towards the Sun. The new features in the orbital evolution under the latitudinally dependent force as compared with the isotropic Poynting-Robertson drag are:

1. not only the semimajor axisa and the eccentricity ε but also the argument of the perihelion ω varies with time,
2. the rate of change ofa, ε, ω depends on the inclination.

An example of particle trajectories in the phase space of elements is presented.