Charged particle motion on and off the equatorial plane of a magnetic star in Rosen's bimetric gravity

In this paper we have studied in detail (numerically) the trajectories of charged particles in a magnetic field (dipolar at infinity) associated with a static star in the framework of Rosen's bimetric theory of gravity. It was found that there do exist potential wells that allow possible trapping of particles in stable orbits both on and off the equatorial plane. A particularly interesting feature that has shown up is the fact that the characteristics of the effective potential wellV _eff depend on the ratio of the magnetic field strength parameter λ, and the angular momentumL of the charged particle. For values lower than a critical (λ/L) _c the potential well lies within the regionr≤2m.     

Periodic orbits under combined effects of oblateness and radiation in the restricted problem of three bodies

In this paper, we study the existence of libration points and their linear stability when the three participating bodies are axisymmetric and the primaries are radiating, we found that the collinear points remain unstable, it is further seen that the triangular points are stable for 0<μ<μ _c , and unstable for where , it is also observed that for these points the range of stability will decrease. In addition to this we have studied periodic orbits around these points in the range 0<μ<μ _c , we found that these orbits are elliptical; the frequencies of long and short orbits of the periodic motion are affected by the terms which involve parameters that characterize the oblateness and radiation repulsive forces. The implication is that the period of long periodic orbits adjusts with the change in its frequency while the period of short periodic orbit will decrease.     

Effect of 3rd-degree gravity harmonics and Earth perturbations on lunar artificial satellite orbits

In a previous work we studied the effects of (I) the J ₂ and C ₂₂ terms of the lunar potential and (II) the rotation of the primary on the critical inclination orbits of artificial satellites. Here, we show that, when 3rd-degree gravity harmonics are taken into account, the long-term orbital behavior and stability are strongly affected, especially for a non-rotating central body, where chaotic or collision orbits dominate the phase space. In the rotating case these phenomena are strongly weakened and the motion is mostly regular. When the averaged effect of the Earth’s perturbation is added, chaotic regions appear again for some inclination ranges. These are more important for higher values of semi-major axes. We compute the main families of periodic orbits, which are shown to emanate from the ‘frozen eccentricity’ and ‘critical inclination’ solutions of the axisymmetric problem (‘J ₂ + J ₃’). Although the geometrical properties of the orbits are not preserved, we find that the variations in e, I and g can be quite small, so that they can be of practical importance to mission planning.     

Artificial Martian frozen orbits and Sun-Synchronous orbits using continuous low-thrust control

This paper analyses three types of artificial orbits around Mars pushed by continuous low-thrust control: artificial frozen orbits, artificial Sun-Synchronous orbits and artificial Sun-Synchronous frozen orbits. These artificial orbits have similar characteristics to natural frozen orbits and Sun-Synchronous orbits, and their orbital parameters can be selected arbitrarily by using continuous low-thrust control. One control strategy to achieve the artificial frozen orbit is using both the transverse and radial continuous low-thrust control, and another to achieve the artificial Sun-Synchronous orbit is using the normal continuous low-thrust control. These continuous low-thrust control strategies consider J ₂, J ₃, and J ₄ perturbations of Mars. It is proved that both control strategies can minimize characteristic velocity. Relevant formulas are derived, and numerical results are presented. Given the same initial orbital parameters, the control acceleration and characteristic velocity taking into account J ₂, J ₃, and J ₄ perturbations are similar to those taking into account J ₂ perturbations for both Mars and the Earth. The control thrust of the orbit around Mars is smaller than that around the Earth. The magnitude of the control acceleration of ASFOM-4 (named as Artificial Sun-Synchronous Frozen Orbit Method 4) is the lowest among these strategies and the characteristic velocity within one orbital period is only 0.5219 m/s for the artificial Sun-Synchronous frozen orbit around Mars. It is evident that the relationship among the control thrusts and the primary orbital parameters of Martian artificial orbits is always similar to that of the Earth. Simulation shows that the control scheme extends the orbital parameters’ selection range of three types of orbits around Mars, compared with the natural frozen orbit and Sun-Synchronous orbit.     

Trajectories of satellites under the combined influences of earth oblateness and air drag

Trajectories of satellites under the influences of earth oblateness and air drag are derived by the asymptotic method in nonlinear mechanics. Based on the assumptions: (1) the dominant oblateness factor of the earth is the second harmonic (J ₂), (2) a non-rotating, spherically symmetric atmosphere and an exponential distribution of atmospheric density, (3) original elliptical orbits being of small eccentricity, closed-form solutions for the improved first order approximation are obtained. After finding the osculating orbital elements of the resulting trajectories, we expose the behavior of osculating orbits at various inclinations.     

Trapping and dynamical evolution of interplanetary dust particles in Earth’s quasi-satellite resonance

We used numerical simulations to model the orbital evolution of interplanetary dust particles (IDPs) evolving inward past Earth’s orbit under the influence of radiation pressure, Poynting–Robertson light drag (PR drag), solar wind drag, and gravitational perturbations from the planets. A series of β values (where β is the ratio of the force from radiation pressure to that of central gravity) were used ranging from 0.0025 up to 0.02. Assuming a composition consistent with astronomical silicate and a particle density of 2.5 g cm⁻³ these β values correspond to dust particle diameters ranging from 200 μm down to 25 μm. As the dust particle orbits decay past 1 AU between 4% (for β = 0.02, or 25 μm) and 40% (for β = 0.0025, or 200 μm) of the population became trapped in 1:1 co-orbital resonance with Earth. In addition to traditional horseshoe type co-orbitals, we found about a quarter of the co-orbital IDPs became trapped as so-called quasi-satellites. Quasi-satellite IDPs always remain relatively near to Earth (within 0.1–0.3 AU, or 10–30 Hill radii, R_H) and undergo two close-encounters with Earth each year. While resonant perturbations from Earth halt the decay in semi-major axis of quasi-satellite IDPs their orbital eccentricities continue to decrease under the influence of PR drag and solar wind drag, forcing the IDPs onto more Earth-like orbits. This has dramatic consequences for the relative velocity and distance of closest approach between Earth and the quasi-satellite IDPs. After 10⁴–10⁵ years in the quasi-satellite resonance dust particles are typically less than 10R_H from Earth and consistently coming within about 3R_H. In the late stages of evolution, as the dust particles are escaping the 1:1 resonance, quasi-satellite IDPs can have deep close-encounters with Earth significantly below R_H. Removing the effects of Earth’s gravitational acceleration reveals that encounter velocities (i.e., velocities “at infinity”) between quasi-satellite IDPs and Earth during these close-encounters are just a few hundred meters per second or slower, well below the average values of 2–4 km s⁻¹ for non-resonant Earth-crossing IDPs with similar initial orbits. These low encounter velocities lead to a factor of 10–100 increase in Earth’s gravitationally enhanced impact cross-section (σ_grav) for quasi-satellite IDPs compared to similar non-resonant IDPs. The enhancement in σ_grav between quasi-satellite IDPs and cometary Earth-crossing IDPs is even more pronounced, favoring accretion of quasi-satellite dust particles by a factor of 100–3000 over the cometary IDPs. This suggests that quasi-satellite dust particles may dominate the flux of large (25–200 μm) IDPs entering Earth’s atmosphere. Furthermore, because quasi-satellite trapping is known to be directly correlated with the host planet’s orbital eccentricity the accretion of quasi-satellite dust likely ebbs and flows on 10⁵ year time scales synchronized with Earth’s orbital evolution.     

Analytical Approach to the Motion of a Lunar Artificial Satellite

The canonical equations of motion of an artificial lunar satellite are formulated including the effects of the asphericity of the Moon comprising the harmonics J ₂, J ₂₂, J ₃, J ₃₁, J ₄ andJ ₅, the oblateness of the Earth up to the second zonal harmonic, as well as the disturbing function due to the attractions of the Earth and of the Sun (terms are retained up to order 10^-6 for the higher orbits and 10^-8 for the lower orbits). This revised version was published online in July 2006 with corrections to the Cover Date.     

On the impact of the Yarkovsky effect on Apophis’ orbit

In October 2009, a new set of optical observations of Apophis, a potentially hazardous asteroid, was published. These data have significantly expanded the interval of observations and their total number. In the article we compare the results of refinement of Apophis’ orbit made at the Jet Propulsion Laboratory (JPL, United States), the University of Pisa (Italy), and the Institute of Applied Astronomy (IAA) of the Russian Academy of Sciences with consideration for new observations. New orbits lead to a significant decrease in the probability of Apophis’ collision with the Earth in 2036. As a result of processing a large number of observations of asteroids approaching the Earth and main belt asteroids less than 40 km in size, with a large number of optical and, in many cases, radar observations in different oppositions, one of the authors revealed that additional acceleration affects their motion. This acceleration can be represented by the transversal component A ₂ in the orbital coordinate system. The presence of this acceleration can be interpreted as the Yarkovsky effect. The statistical properties of distribution of A ₂ for asteroids, for which it is determined quite reliably, evidence in favor of this interpretation. The value of additional acceleration for bodies the size of Apophis falls in the range ±10⁻¹³ AU/day². In this paper we have calculated the probability of Apophis colliding with the Earth in 2036 at different values of the transversal component of additional acceleration A ₂. For the resulting points, a plot of the probability of the collision against the A ₂ value has been constructed. At A ₂ = −8.748 × 10⁻¹⁴ AU/day² (and zero values of the radial A ₁ and normal A ₃ components) the nominal solution for Apophis’ orbit on April 13, 2029, is only 90 m from the middle of a “keyhole” 600 m in width, which leads to a collision of Apophis with the Earth in 2036. Since the scattering ellipse in the target plane in 2029 significantly overlaps the keyhole, the probability of collision at the given additional acceleration value is 0.0022. This result has been verified by the Monte Carlo method. Tests of 10000 random sets of orbital elements, which were found taking into account their correlation, have shown that 22 cases have resulted in virtual asteroids colliding with Earth in 2036. A plot of the probability of the collision against the value of A ₂ has been constructed.     

Periodic orbits generated by Lagrangian solutions of the restricted three body problem when one of the primaries is an oblate body

We have studied periodic orbits generated by Lagrangian solutions of the restricted three body problem when one of the primaries is an oblate body. We have determined the periodic orbits for different values of μ, h and A (h is energy constant, μ is mass ratio of the two primaries and A is an oblateness factor). These orbits have been determined by giving displacements along the tangent and normal to the mobile coordinates as defined by Karimov and Sokolsky (Celest. Mech. 46:335, 1989). These orbits have been drawn by using the predictor-corrector method. We have also studied the effect of oblateness by taking some fixed values of μ, A and h. As starters for our method, we use some known periodic orbits in the classical restricted three body problem.     

Third-order secular Solution of the variational equations of motion of a satellite in orbit around a non-spherical planet

We constructed an analytical theory of satellite motion up to the third order relative to the oblateness parameter of the Earth (J ₂). Equations of secular variations was developed for the first three orbital elements (a, e, i) of an artificial satellite. The secular variations are solved in a closed form.     

Constraining the relative inclinations of the planets B and C of the millisecond pulsar PSR B1257+12

We investigate on the relative inclination of the planets B and C orbiting the pulsar PSR B1257+12. First, we show that the third Kepler’s law does represent an adequate model for the orbital periods P of the planets, because other Newtonian and Einsteinian corrections are orders of magnitude smaller than the accuracy in measuring P _B/C. Then, on the basis of available timing data, we determine the ratio sin i _C/ sin i _B = 0.92±0.05 of the orbital inclinations i _B and i _C independently of the pulsar’s mass M. It turns out that coplanarity of the orbits of B and C would imply a violation of the equivalence principle. Adopting a pulsar mass range 1 ≲ M ≲ 3, in solar masses (supported by present-day theoretical and observational bounds for pulsar’s masses), both face-on and edge-on orbital configurations for the orbits of the two planets are ruled out; the acceptable inclinations for B span the range 36 deg ≲ i _B ≲ 66 deg, with a corresponding relative inclination range 6 deg ≲ (i _C − i _B) ≲ 13 deg.     

Three-dimensional time optimal double angular momentum reversal trajectory using solar sails

A new concept of three dimensional non-Keplerian trajectories with double angular momentum reversal is investigated with high performance solar sails. The main discussion of this paper is about such 3D solar inverse orbits with inner constraints. The problem is addressed in a time optimal control framework solved by an indirect method. Two typical solar inverse orbits have been achieved and presented in a 3D non-dimensional dynamic model in the Heliocentric Inertial Frame. Starting from the Earth orbit ecliptic plane, a sailcraft in the inverse orbit exhibits a butterfly shape trajectory. As such, the new orbits are symmetrical with respect to a plane which contains the Sun-perihelion line. The relation of the sail attitude angles between the two symmetrical parts of the orbits are used to reduce the simulation effort. The quasi-heliostationary property at its aphelia is demonstrated with variation of the orbital radius. Evolutions of the orbital velocity and optimal sail orientations are also outlined and discussed to benefit future design work. As is suited for space observation guaranteed by its butterfly shape, the inverse orbits are thoroughly studied in terms of the concerned parameters. The discussion of the parametric influence is ranked in order as perihelion distance r _E , required maximum position z _max, perihelion position z _f and the sail lightness number β. Suitable ranges of each parameter are adopted to illustrate the orbital variation trend. Through numerical simulations the features of such inverse orbits are further emphasized to provide an initial reference for future researchers.     

Stability regions of equilibrium points in restricted four-body problem with oblateness effects

In this paper, we extend the basic model of the restricted four-body problem introducing two bigger dominant primaries m ₁ and m ₂ as oblate spheroids when masses of the two primary bodies (m ₂ and m ₃) are equal. The aim of this study is to investigate the use of zero velocity surfaces and the Poincaré surfaces of section to determine the possible allowed boundary regions and the stability orbit of the equilibrium points. According to different values of Jacobi constant C, we can determine boundary region where the particle can move in possible permitted zones. The stability regions of the equilibrium points expanded due to presence of oblateness coefficient and various values of C, whereas for certain range of t (100≤t≤200), orbits form a shape of cote’s spiral. For different values of oblateness parameters A ₁ (0<A ₁<1) and A ₂ (0<A ₂<1), we obtain two collinear and six non-collinear equilibrium points. The non-collinear equilibrium points are stable when the mass parameter μ lies in the interval (0.0190637,0.647603). However, basins of attraction are constructed with the help of Newton Raphson method to demonstrate the convergence as well as divergence region of the equilibrium points. The nature of basins of attraction of the equilibrium points are less effected in presence and absence of oblateness coefficients A ₁ and A ₂ respectively in the proposed model.     

Sun-synchronous orbits near critical inclination

The long period dynamics of Sun-synchronous orbits near the critical inclination 116.6° are investigated. It is known that, at the critical inclination, the average perigee location is unchanged by Earth oblateness. For certain values of semimajor axis and eccentricity, orbit plane precession caused by Earth oblateness is synchronous with the mean orbital motion of the apparent Sun (a Sun-synchronism). Sun-synchronous orbits have been used extensively in meteorological and remote sensing satellite missions. Gravitational perturbations arising from an aspherical Earth, the Moon, and the Sun cause long period fluctuations in the mean argument of perigee, eccentricity, inclination, and ascending node. Double resonance occurs because slow oscillations in the perigee and Sun-referenced ascending node are coupled through the solar gravity gradient. It is shown that the total number and infinitesimal stability of equilibrium solutions can change abruptly over the Sun-synchronous range of semimajor axis values (1.54 to 1.70 Earth radii). The effect of direct solar radiation pressure upon certain stable equilibria is investigated.     

Acceleration of Relativistic Protons During the 20 January 2005 Flare and CME

The origin of relativistic solar protons during large flare/CME events has not been uniquely identified so far. We perform a detailed comparative analysis of the time profiles of relativistic protons detected by the worldwide network of neutron monitors at Earth with electromagnetic signatures of particle acceleration in the solar corona during the large particle event of 20 January 2005. The intensity – time profile of the relativistic protons derived from the neutron monitor data indicates two successive peaks. We show that microwave, hard X-ray, and γ-ray emissions display several episodes of particle acceleration within the impulsive flare phase. The first relativistic protons detected at Earth are accelerated together with relativistic electrons and with protons that produce pion-decay γ rays during the second episode. The second peak in the relativistic proton profile at Earth is accompanied by new signatures of particle acceleration in the corona within ≈1R _⊙ above the photosphere, revealed by hard X-ray and microwave emissions of low intensity and by the renewed radio emission of electron beams and of a coronal shock wave. We discuss the observations in terms of different scenarios of particle acceleration in the corona.     

Satellite motion in a non-singular gravitational potential

We study the effects of a non-singular gravitational potential on satellite orbits by deriving the corresponding time rates of change of its orbital elements. This is achieved by expanding the non-singular potential into power series up to second order. This series contains three terms, the first been the Newtonian potential and the other two, here R ₁ (first order term) and R ₂ (second order term), express deviations of the singular potential from the Newtonian. These deviations from the Newtonian potential are taken as disturbing potential terms in the Lagrange planetary equations that provide the time rates of change of the orbital elements of a satellite in a non-singular gravitational field. We split these effects into secular, low and high frequency components and we evaluate them numerically using the low Earth orbiting mission Gravity Recovery and Climate Experiment (GRACE). We show that the secular effect of the second-order disturbing term R ₂ on the perigee and the mean anomaly are 4″.307×10⁻⁹/a, and −2″.533×10⁻¹⁵/a, respectively. These effects are far too small and most likely cannot easily be observed with today’s technology. Numerical evaluation of the low and high frequency effects of the disturbing term R ₂ on low Earth orbiters like GRACE are very small and undetectable by current observational means.     

Poynting–Robertson effect in non-singular gravitational potential

We use a non-singular potential that appears in the literature under the influence of which the Poynting-Robertson effect is studied. For that, dust particles originating within the asteroid belt are used, in circular and elliptic orbits, and expressions for the semimajor axis as a function of time are obtained. The derived expressions are written in terms of the two basic dust particle parameters, namely the density and the diameter. In both cases, we obtain expressions for the time that the dust particles take to reach the orbit of Earth under the action of the non-singular potential and solar radiation. For the non-singular potential, dust particles of diameter 10⁻³ m in circular and elliptical orbits require times of the order of 4.058×10⁷ and 2.823×10⁷ y to reach the orbit of the Earth respectively. Finally, the derived expressions and numerical results are compared with those of the Newtonian potential.     

Artificial halo orbits for low-thrust propulsion spacecraft

We consider periodic halo orbits about artificial equilibrium points (AEP) near to the Lagrange points L ₁ and L ₂ in the circular restricted three body problem, where the third body is a low-thrust propulsion spacecraft in the Sun–Earth system. Although such halo orbits about artificial equilibrium points can be generated using a solar sail, there are points inside L ₁ and beyond L ₂ where a solar sail cannot be placed, so low-thrust, such as solar electric propulsion, is the only option to generate artificial halo orbits around points inaccessible to a solar sail. Analytical and numerical halo orbits for such low-thrust propulsion systems are obtained by using the Lindstedt Poincaré and differential corrector method respectively. Both the period and minimum amplitude of halo orbits about artificial equilibrium points inside L ₁ decreases with an increase in low-thrust acceleration. The halo orbits about artificial equilibrium points beyond L ₂ in contrast show an increase in period with an increase in low-thrust acceleration. However, the minimum amplitude first increases and then decreases after the thrust acceleration exceeds 0.415 mm/s². Using a continuation method, we also find stable artificial halo orbits which can be sustained for long integration times and require a reasonably small low-thrust acceleration 0.0593 mm/s².     

Orbital dynamics of the gravitational field in Bardeen space-time

In this paper, the orbital dynamics of the gravitational field in Bardeen space-time are investigated. The equations of the particle’s orbital motion are given by solving the Lagrangian equation, and the stability and types of orbits are studied by means of analysing the effective potential of particles. Also, with the help of phase-plane method, the closed and non-closed orbits of test particle are analysed. We find that the stability and types of orbits in the Bardeen space-time are determined not only by the particles’ energy but also by the angular momentum. And for q=0.5M and b<3.3731M, absorbed by the black hole is the only fate of the test particle. We also find that the position of the innermost stable circular orbit of Bardeen black hole occurs at r _min =5.5722M.     

The Poynting-Robertson effect in the Newtonian potential with a Yukawa correction

We consider a Yukawa-type gravitational potential combined with the Poynting-Robertson effect. Dust particles originating within the asteroid belt and moving on circular and elliptic trajectories are studied and expressions for the time rate of change of their orbital radii and semimajor axes, respectively, are obtained. These expressions are written in terms of basic particle parameters, namely their density and diameter. Then, they are applied to produce expressions for the time required by the dust particles to reach the orbit of Earth. For the Yukawa gravitational potential, dust particles of diameter \(10^{ - 3}\) m in circular orbits require times of the order of \(8.557 \times 10^{6}\) yr and for elliptic orbits of eccentricities \(e =0.1, 0.5\) require times of \(9.396 \times 10^{6}\) and \(2.129 \times 10^{6}\) yr respectively to reach Earth’s orbit. Finally, various cases of the Yukawa potential are studied and the corresponding particle times to reach Earth’s are derived per case along with numerical results for circular and various elliptical orbits.