Attitude stability of a spacecraft on a stationary orbit around an asteroid subjected to gravity gradient torque

Attitude stability of spacecraft subjected to the gravity gradient torque in a central gravity field has been one of the most fundamental problems in space engineering since the beginning of the space age. Over the last two decades, the interest in asteroid missions for scientific exploration and near-Earth object hazard mitigation is increasing. In this paper, the problem of attitude stability is generalized to a rigid spacecraft on a stationary orbit around a uniformly-rotating asteroid. This generalized problem is studied via the linearized equations of motion, in which the harmonic coefficients $C_{20}$ and $C_{22}$ of the gravity field of the asteroid are considered. The necessary conditions of stability of this conservative system are investigated in detail with respect to three important parameters of the asteroid, which include the harmonic coefficients $C_{20}$ and $C_{22}$ , as well as the ratio of the mean radius to the radius of the stationary orbit. We find that, due to the significantly non-spherical shape and the rapid rotation of the asteroid, the attitude stability domain is modified significantly in comparison with the classical stability domain predicted by the Beletskii–DeBra–Delp method on a circular orbit in a central gravity field. Especially, when the spacecraft is located on the intermediate-moment principal axis of the asteroid, the stability domain can be totally different from the classical stability domain. Our results are useful for the design of attitude control system in the future asteroid missions.     

Sun–Earth L_{1} and L_{2} to Moon transfers exploiting natural dynamics

This paper examines the design of transfers from the Sun–Earth libration orbits, at the \(L_{1}\) and \(L_{2}\) points, towards the Moon using natural dynamics in order to assess the feasibility of future disposal or lifetime extension operations. With an eye to the probably small quantity of propellant left when its operational life has ended, the spacecraft leaves the libration point orbit on an unstable invariant manifold to bring itself closer to the Earth and Moon. The total trajectory is modeled in the coupled circular restricted three-body problem, and some preliminary study of the use of solar radiation pressure is also provided. The concept of survivability and event maps is introduced to obtain suitable conditions that can be targeted such that the spacecraft impacts, or is weakly captured by, the Moon. Weak capture at the Moon is studied by method of these maps. Some results for planar Lyapunov orbits at \(L_{1}\) and \(L_{2}\) are given, as well as some results for the operational orbit of SOHO.     

Long term stability of Earth Trojans

We explore the long-term stability of Earth Trojans by using a chaos indicator, the Frequency Map Analysis. We find that there is an extended stability region at low eccentricity and for inclinations lower than about $50^{\circ }$ even if the most stable orbits are found at $i \le 40^{\circ }$ . This region is not limited in libration amplitude, contrary to what found for Trojan orbits around outer planets. We also investigate how the stability properties are affected by the tidal force of the Earth–Moon system and by the Yarkovsky force. The tidal field of the Earth–Moon system reduces the stability of the Earth Trojans at high inclinations while the Yarkovsky force, at least for bodies larger than 10 m in diameter, does not seem to strongly influence the long-term stability. Earth Trojan orbits with the lowest diffusion rate survive on timescales of the order of $10^9$  years but their evolution is chaotic. Their behaviour is similar to that of Mars Trojans even if Earth Trojans appear to have shorter lifetimes.     

High precision predictions for near-Earth asteroids: the strange case of (3908) Nyx

In November 2004 radar delay measurements of near-Earth asteroid (3908) Nyx obtained at the Arecibo radio telescope turned out to be \(7.5\sigma \) away from the orbital prediction. We prove that this discrepancy was caused by a poor astrometric treatment and an incomplete dynamical model, which did not account for nongravitational perturbations. To improve the astrometric treatment, we remove known star catalog biases, apply suitable weights to the observations, and use an aggressive outlier rejection scheme. The main issue related to the dynamical model is having not accounted for the Yarkovsky effect. Including the Yarkovsky perturbation in the model makes the orbital prediction and the radar measurements statistically consistent by both reducing the offset and increasing the prediction uncertainty to a more realistic level. This analysis shows the sensitivity of high precision predictions to the astrometric treatment and the Yarkovsky effect. By using the full observational dataset we obtain a \(5\sigma \) detection of the Yarkovsky effect acting on Nyx corresponding to an orbital drift \(da/dt = (142 \pm 29)\)  m/year. In turn, we derive constraints on thermal inertia and bulk density. In particular, we find that the bulk density of Nyx is around 1 g/cm \(^3\) , possibly less. To make sure that our results are not corrupted by an asteroid impact or a close approach with a perturbing asteroid not included in our dynamical model, we show that the astrometry provides no convincing evidence of an impulsive variation of Nyx’s velocity while crossing the main belt region.     

Evolution of the \mathcal {L}_1 halo family in the radial solar sail circular restricted three-body problem

We present a detailed investigation of the dramatic changes that occur in the \(\mathcal {L}_1\) halo family when radiation pressure is introduced into the Sun–Earth circular restricted three-body problem (CRTBP). This photo-gravitational CRTBP can be used to model the motion of a solar sail orientated perpendicular to the Sun-line. The problem is then parameterized by the sail lightness number, the ratio of solar radiation pressure acceleration to solar gravitational acceleration. Using boundary-value problem numerical continuation methods and the AUTO software package (Doedel et al. in Int J Bifurc Chaos 1:493–520, 1991) the families can be fully mapped out as the parameter \(\beta \) is increased. Interestingly, the emergence of a branch point in the retrograde satellite family around the Earth at \(\beta \approx 0.0387\) acts to split the halo family into two new families. As radiation pressure is further increased one of these new families subsequently merges with another non-planar family at \(\beta \approx 0.289\) , resulting in a third new family. The linear stability of the families changes rapidly at low values of \(\beta \) , with several small regions of neutral stability appearing and disappearing. By using existing methods within AUTO to continue branch points and period-doubling bifurcations, and deriving a new boundary-value problem formulation to continue the folds and Krein collisions, we track bifurcations and changes in the linear stability of the families in the parameter \(\beta \) and provide a comprehensive overview of the halo family in the presence of radiation pressure. The results demonstrate that even at small values of \(\beta \) there is significant difference to the classical CRTBP, providing opportunity for novel solar sail trajectories. Further, we also find that the branch points between families in the solar sail CRTBP provide a simple means of generating certain families in the classical case.     

Cosmic-ray antimatter: A primary origin hypothesis

We examine the possibility that the observed cosmic-ray protons are of primary extragalactic origin. The present \(\bar p\) data are consistent with a primary extragalactic component having \(\bar p\) /p?3.2±0.7 x 10^-4 independent of energy. Following the suggestion that most extragalactic cosmic rays are from active galaxies, we propose that most of the observed \(\bar p\) 's are alos from the same sites. This would imply the possibility of destroying the corresponding \(\bar \alpha \) 'sat the source, thus leading to a flux ratio \(\bar \alpha \) /α< \(\bar p\) /p. We further predict an estimate for \(\bar \alpha \) α～10^-5, within the range of future cosmic-ray detectors. the cosmological implications of this proposal are discussed.     

The broad band spectral index of TeV blazars detected by Fermi LAT

Using γ-ray data detected by Fermi Large Area Telescope (LAT) and multi-wave band data for 35 TeV blazars sample, we have studied the possible correlations between different broad band spectral indices ( $\alpha_{\rm r.ir}$ , $\alpha_{\rm{r.o}}$ , $\alpha_{\rm r.x}$ , $\alpha_{\rm r.\gamma}$ , $\alpha_{\rm{ir.o}}$ , $\alpha_{\rm ir.x}$ , $\alpha_{\rm ir.\gamma}$ , $\alpha_{\rm o.x}$ , $\alpha_{\rm o.\gamma}$ , $\alpha_{\rm r.x}$ , $\alpha_{\rm x.\gamma}$ ) in all states (average/high/low). Our results are as follows: (1) For our TeV blazars sample, the strong positive correlations were found between $\alpha_{\rm r.ir}$ and $\alpha_{\rm{r.o}}$ , between $\alpha_{\rm r.ir}$ and $\alpha_{\rm r.x}$ , between $\alpha_{\rm r.ir}$ and $\alpha_{\rm r.\gamma}$ in all states (average/high/low); (2) For our TeV blazars sample, the strong anti-correlations were found between $\alpha_{\rm r.ir}$ and $\alpha_{\rm x.\gamma}$ , between $\alpha_{\rm{r.o}}$ and $\alpha_{\rm ir.\gamma}$ , between $\alpha_{\rm{r.o}}$ and $\alpha_{\rm o.\gamma}$ , between $\alpha_{\rm{r.o}}$ and $\alpha_{\rm x.\gamma}$ , between $\alpha_{\mathrm{ir.o}}$ and $\alpha_{\rm o.\gamma}$ , between $\alpha_{\rm r.x}$ and $\alpha_{\rm x.\gamma}$ , between $\alpha_{\rm ir.x}$ and $\alpha_{\rm x.\gamma}$ in all states (average/high/low). The results suggest that the synchrotron self-Compton radiation (SSC) is the main mechanism of high energy γ-ray emission and the inverse Compton scattering of circum-nuclear dust is likely to be a important complementary mechanism for TeV blazars. Our results also show that the possible correlations vary from state to state in the same pair of indices, Which suggest that there may exist differences in the emitting process and in the location of the emitting region for different states.     

Secular dynamics of a lunar orbiter: a global exploration using Prony’s frequency analysis

We study the secular dynamics of lunar orbiters, in the framework of high-degree gravity models. To achieve a global view of the dynamics, we apply a frequency analysis (FA) technique which is based on Prony’s method. This allows for an extensive exploration of the eccentricity ( $e$ )—inclination ( $i$ ) space, based on short-term integrations ( $\sim $ 8 months) over relatively high-resolution grids of initial conditions. Different gravity models are considered: 3rd, 7th and 10th degree in the spherical harmonics expansion, with the main perturbations from the Earth being added. Since the dynamics is mostly regular, each orbit is characterised by a few parameters, whose values are given by the spectral decomposition of the orbital elements time series. The resulting frequency and amplitude maps in ( $e_0,i_0$ ) are used to identify the dominant perturbations and deduce the “minimum complexity” model necessary to capture the essential features of the long-term dynamics. We find that the 7th degree zonal harmonic ( $J_7$ term) is of profound importance at low altitudes as, depending on the initial secular phases, it can lead to collision with the Moon’s surface within a few months. The 3rd-degree non-axisymmetric terms are enough to describe the deviations from the 1 degree-of-freedom zonal problem; their main effect is to modify the equilibrium value of the argument of periselenium, $\omega $ , with respect to the “frozen” solution ( $\omega =\pm 90^{\circ }, \forall \Omega $ , where $\Omega $ is the nodal longitude). Finally, we show that using FA on a fine grid of initial conditions, set around a suitably chosen ‘first guess’, one can compute an accurate approximation of the initial conditions of a periodic orbit.     

On the “blue sky catastrophe” termination in the restricted four-body problem

The restricted three-body problem (R3BP) possesses the property that some classes of doubly asymptotic (i.e., homoclinic or heteroclinic) orbits are limit members of families of periodic orbits, this phenomenon has been known as the “blue sky catastrophe” termination principle. A similar case occurs in the restricted four body problem for the collinear equilibrium point $L_{2}$ L 2 . In the restricted four body problem with primaries in a triangle relative equilibrium, we show that the same phenomenon observed in the R3BP occurs. We prove that there exists a critical value of the mass parameter $\mu _{b}$ μ b such that for $\mu =\mu _{b}$ μ = μ b a Hamiltonian Hopf bifurcation takes place. Moreover we show that for $\mu >\mu _{b}$ μ > μ b the stable and unstable manifolds of $L_{2}$ L 2 intersect transversally and the spectrum corresponds to a complex saddle. This proves that Henrard’s theorem applies at least for $\mu $ μ close to $\mu _{b}$ μ b . In particular there exists a family of periodic orbits having the homoclinic orbit as a limit.     

Solving Kepler’s equation via Smale’s \alpha -theory

We obtain an approximate solution $\tilde{E}=\tilde{E}(e,M)$ of Kepler’s equation $E-e\sin (E)=M$ for any $e\in [0,1)$ and $M\in [0,\pi ]$ . Our solution is guaranteed, via Smale’s $\alpha $ -theory, to converge to the actual solution $E$ through Newton’s method at quadratic speed, i.e. the $n$ -th iteration produces a value $E_n$ such that $|E_n-E|\le (\frac{1}{2})^{2^n-1}|\tilde{E}-E|$ . The formula provided for $\tilde{E}$ is a piecewise rational function with conditions defined by polynomial inequalities, except for a small region near $e=1$ and $M=0$ , where a single cubic root is used. We also show that the root operation is unavoidable, by proving that no approximate solution can be computed in the entire region $[0,1)\times [0,\pi ]$ if only rational functions are allowed in each branch.     

Coupling between corotation and Lindblad resonances in the presence of secular precession rates

We investigate the dynamics of two satellites with masses $\mu _s$ and $\mu '_s$ orbiting a massive central planet in a common plane, near a first order mean motion resonance $m+1{:}m$ (m integer). We consider only the resonant terms of first order in eccentricity in the disturbing potential of the satellites, plus the secular terms causing the orbital apsidal precessions. We obtain a two-degrees-of-freedom system, associated with the two critical resonant angles $\phi = (m+1)\lambda ' -m\lambda - \varpi $ and $\phi '= (m+1)\lambda ' -m\lambda - \varpi '$ , where $\lambda $ and $\varpi $ are the mean longitude and longitude of periapsis of $\mu _s$ , respectively, and where the primed quantities apply to $\mu '_s$ . We consider the special case where $\mu _s \rightarrow 0$ (restricted problem). The symmetry between the two angles $\phi $ and $\phi '$ is then broken, leading to two different kinds of resonances, classically referred to as corotation eccentric resonance (CER) and Lindblad eccentric Resonance (LER), respectively. We write the four reduced equations of motion near the CER and LER, that form what we call the CoraLin model. This model depends upon only two dimensionless parameters that control the dynamics of the system: the distance $D$ between the CER and LER, and a forcing parameter $\epsilon _L$ that includes both the mass and the orbital eccentricity of the disturbing satellite. Three regimes are found: for $D=0$ the system is integrable, for $D$ of order unity, it exhibits prominent chaotic regions, while for $D$ large compared to 2, the behavior of the system is regular and can be qualitatively described using simple adiabatic invariant arguments. We apply this model to three recently discovered small Saturnian satellites dynamically linked to Mimas through first order mean motion resonances: Aegaeon, Methone and Anthe. Poincaré surfaces of section reveal the dynamical structure of each orbit, and their proximity to chaotic regions. This work may be useful to explore various scenarii of resonant capture for those satellites.     

An analytical proof on certain determinants connected with the collinear central configurations in the n-body problem

In this paper we give a short analytical proof of the inequalities proved by Albouy–Moeckel through computer algebra, in the cases $n=5$ and $n=6$ . These inequalities guarantee that, in the $n$ -body problem, the family of mass vectors making a given collinear configuration a central configuration is 2-dimensional. The induction techniques here may be used to prove the inequalities for general $n$ with more subtle estimation but currently the inequalities still remains unproved for $n\ge 7$ .     

F and G Taylor series solutions to the Stark and Kepler problems with Sundman transformations

The classic $F$ and $G$ Taylor series of Keplerian motion are extended to solve the Stark problem and to use the generalized Sundman transformation. Exact recursion formulas for the series coefficients are derived, and the method is implemented to high order via a symbolic manipulator. The results lead to fast and accurate propagation models with efficient discretizations. The new $F$ and $G$ Stark series solutions are compared to the Modern Taylor Series (MTS) and 8th order Runge–Kutta–Fehlberg (RKF8) solutions. In terms of runtime, the $F$ and $G$ approach is shown to compare favorably to the MTS method up to order 20, and both Taylor series methods enjoy approximate order of magnitude speedups compared to RKF8 implementations. Actual runtime is shown to vary with eccentricity, perturbation size, prescribed accuracy, and the Sundman power law. The method and results are valid for both the Stark and the Kepler problems. The effects of the generalized Sundman transformation on the accuracy of the propagation are analyzed. The Taylor series solutions are shown to be exceptionally efficient when the unity power law from the classic Sundman transformation is applied. An example low-thrust trajectory propagation demonstrates the utility of the $F$ and $G$ Stark series solutions.     

On some long time dynamical features of the Trojan asteroids of Jupiter

The equation of motion of long periodic libration around the Lagrangian point $L_4$ L 4 in the restricted three-body problem is investigated. The range of validity of an approximate analytical solution in the tadpole region is determined by numerical integration. The predictions of the model of libration are tested on the Trojan asteroids of Jupiter. The long time evolution of the orbital eccentricity and the longitude of the perihelion of the Trojan asteroids, under the effect of the four giant planets, is also investigated and a slight dynamical asymmetry is shown between the two groups of Trojans at $L_4$ L 4 and $L_5$ L 5 .     

Integral wide-field spectroscopy in astronomy: the Imaging FTS solution

Long-slit grating spectrometers in scanning mode and Fabry–Perot interferometers as tunable filters are commonly used to perform integral wide-field spectroscopy on extended astrophysical objects as HII regions and nearby galaxies. The goal of this paper is to demonstrate, by comparison, through a thorough review of the imaging Fourier transform spectrometer (IFTS) properties, that this instrument represents another interesting solution. After a brief recall of the performances, regarding FOV and spectral resolution, of the grating spectrometer, without and with integral field units (IFU), and of the imaging Fabry–Perot, it is demonstrated that for an IFTS the product of the maximum resolution R by the entrance beam étendue U is equal to $2.6\,N\times S_I$ with $N\,\times \,N$ the number of pixels of the detector array and S $_I$ the area of the interferometer beamsplitter. As a consequence, the IFTS offers the most flexible choice of field size and spectral resolution, up to high values for both parameters. It also presents on a wide field an important multichannel advantage in comparison to integral field grating spectrometers, even with multiple IFUs. To complete, the few astronomical IFTSs, built behind ground-based telescopes and in space, for the visible range up to the sub-millimetric domain, are presented. Through two wide-field IFTS projects, one in the visible, the other one in the mid-infrared, the question is addressed of the practical FOV and resolution limits, set by the optical design of the instrument, which can be achieved. Within the 0.3 to $\sim $ 2.5 $\upmu$ m domain, a Michelson interferometer with wide-field diopric collimators provides the easiest solution. This design is illustrated by a $11^{\prime}\times 11^{\prime}$ -field IFTS in the 0.35–0.90 $\upmu$ m range around an off-axis interferometer, called SITELLE, proposed for the 3.6-m CFH Telescope. At longer wavelengths, an all-mirror optics is required, as studied for a spaceborne IFTS, H2EX, for the 8–29 $\upmu$ m range, a $20^{\prime} \times 20^{\prime}$ field, and a high resolution of $\simeq 3\times 10^4$ at 10 $\upmu$ m. To comply with these characteristics, the interferometer is designed with cat’s eye retroreflectors. In the same domain and up to the far infrared, if the instrument aims only at a low spectral resolution (few thousands) and a smaller field (few arcmins $^2$ ), roof-top or corner cube mirrors, as for the IFTS SPIRE on the Herschel space telescope, are usable. At last, perspectives are opened, behind an ELT in the visible and the near infrared with the SITELLE optical combination, in the 2–5 $\upmu$ m on the Antarctic plateau or in space up to longer wavelengths, with the H2EX design, to provide the missing capability of global high spectral resolution studies of extended sources, from comets to distant galaxy clusters.     

Continua of central configurations with a negative mass in the n-body problem

The number of equivalence classes of central configurations of $n \le 4$ bodies of positive mass is known to be finite, but it remains to be shown if this is true for $n \ge 5$ . By allowing one mass to be negative, Gareth Roberts constructed a continuum of inequivalent planar central configurations of $n = 5$ bodies. We reinterpret Roberts’ example and generalize the construction of his continuum to produce a family of continua of central configurations, each with a single negative mass. These new continua exist in even dimensional spaces $\mathbb R ^k$ for $k \ge 4$ .     

High order normal form stability estimates for co-orbital motion

We present estimates of the size of the analytic domain of stability for co-orbital motions obtained by a high order normal form in the framework of the elliptic restricted three body problem. As a demonstration example, we consider the motion of a Trojan body in an extrasolar planetary system with a giant planet of mass parameter $\mu =0.005$ μ = 0.005 and eccentricity $e^{\prime }=0.1$ e ′ = 0.1 . The analysis contains three basic steps: (i) derivation of an accurate expansion of the Hamiltonian, (ii) computation of the normal form up to an optimal order (in the Nekhoroshev sense), and (iii) computation of the optimal size of the remainder at various values of the action integrals (proper elements) of motion. We explain our choice of variables as well as the method used to expand the Hamiltonian so as to ensure a precise model. We then compute the normal form up to the normalisation order $r=50$ r = 50 by use of a computer-algebraic program. We finally estimate the size $||R||$ | | R | | of the remainder as a function of the normalization order, and compute the optimal normalization order at which the remainder becomes minimum. It is found that the optimal value $\log (||R_{opt}||)$ log ( | | R o p t | | ) can serve in order to construct a stability map for the domain of co-orbital motion using only series. This is compared to the stability map found by a purely numerical approach based on chaotic indicators.     

Tidal synchronization of close-in satellites and exoplanets. A rheophysical approach

This paper presents a new theory of the dynamical tides of celestial bodies. It is founded on a Newtonian creep instead of the classical delaying approach of the standard viscoelastic theories and the results of the theory derive mainly from the solution of a non-homogeneous ordinary differential equation. Lags appear in the solution but as quantities determined from the solution of the equation and are not arbitrary external quantities plugged in an elastic model. The resulting lags of the tide components are increasing functions of their frequencies (as in Darwin’s theory), but not small quantities. The amplitudes of the tide components depend on the viscosity of the body and on their frequencies; they are not constants. The resulting stationary rotations (often called pseudo-synchronous) have an excess velocity roughly proportional to $6ne^2/(\chi ^2+\chi ^{-2})$ ( $\chi $ is the mean-motion in units of one critical frequency—the relaxation factor—inversely proportional to the viscosity) instead of the exact $6ne^2$ of standard theories. The dissipation in the pseudo-synchronous solution is inversely proportional to $(\chi +\chi ^{-1})$ ; thus, in the inviscid limit, it is roughly proportional to the frequency (as in standard theories), but that behavior is inverted when the viscosity is high and the tide frequency larger than the critical frequency. For free rotating bodies, the dissipation is given by the same law, but now $\chi $ is the frequency of the semi-diurnal tide in units of the critical frequency. This approach fails, however, to reproduce the actual tidal lags on Earth. In this case, to reconcile theory and observations, we need to assume the existence of an elastic tide superposed to the creeping tide. The theory is applied to several Solar System and extrasolar bodies and currently available data are used to estimate the relaxation factor $\gamma $ (i.e. the critical frequency) of these bodies.     

Retrograde resonance in the planar three-body problem

We continue the investigation of the dynamics of retrograde resonances initiated in Morais and Giuppone (Mon Notices R Astron Soc 424:52–64, doi:10.1111/j.1365-2966.2012.21151.x, 2012). After deriving a procedure to deduce the retrograde resonance terms from the standard expansion of the three-dimensional disturbing function, we concentrate on the planar problem and construct surfaces of section that explore phase-space in the vicinity of the main retrograde resonances (2/ $-$ 1, 1/ $-$ 1 and 1/ $-$ 2). In the case of the 1/ $-$ 1 resonance for which the standard expansion is not adequate to describe the dynamics, we develop a semi-analytic model based on numerical averaging of the unexpanded disturbing function, and show that the predicted libration modes are in agreement with the behavior seen in the surfaces of section.     

Non-minimal derivatively coupled quintessence in the Palatini formalism

The quintessence dark energy model with a kinetic coupling to gravity within the Palatini formalism is studied in this paper. Two different coupling forms: $\hat{R}\partial^{\mu}\phi\partial_{\mu}\phi$ and $\hat {R}_{\mu\nu}\partial^{\mu}\phi\partial^{\nu}\phi$ are analyzed, respectively. We find that both the model with the $\hat{R}\partial^{\mu}\phi\partial_{\mu}\phi$ coupling and the one with the $\hat{R}_{\mu\nu}\partial^{\mu}\phi\partial^{\nu}\phi$ coupling can realize the phantom divide line crossing from phantom to quintessence at late time for its effective equation-of-state. Furthermore, the former can behave like phantom. These features are different from those found in the $\hat {R}\phi^{2}$ coupling case.