Orbital stability close to asteroid 624 Hektor using the polyhedral model

We investigate the orbital stability close to the unique L4-point Jupiter binary Trojan asteroid 624 Hektor. The gravitational potential of 624 Hektor is calculated using the polyhedron model with observational data of 2038 faces and 1021 vertexes. Previous studies have presented three different density values for 624 Hektor. The equilibrium points in the gravitational potential of 624 Hektor with different density values have been studied in detail. There are five equilibrium points in the gravitational potential of 624 Hektor no matter the density value. The positions, Jacobian, eigenvalues, topological cases, stability, as well as the Hessian matrix of the equilibrium points are investigated. For the three different density values the number, topological cases, and the stability of the equilibrium points with different density values are the same. However, the positions of the equilibrium points vary with the density value of the asteroid 624 Hektor. The outer equilibrium points move away from the asteroid’s mass center when the density increases, and the inner equilibrium point moves close to the asteroid’s mass center when the density increases. There exist unstable periodic orbits near the surface of 624 Hektor. We calculated an orbit near the primary’s equatorial plane of this binary Trojan asteroid; the results indicate that the orbit remains stable after 28.8375?d.     

Low-altitude,near-polar and near-circular orbits around Europa

The dynamics of orbits around planetary satellites, taking into account the gravitational attraction of a third-body and the non-uniform distribution of mass of the planetary satellite, is studied. The Hamiltonian considered is explicitly time-dependent. Conditions for frozen orbits are presented. Low-altitude, near-polar orbits, very desirable for scientific missions to study planetary satellites such as the Jupiter’s moon Europa, are analyzed. Lifetimes for these orbits are computed through the single and double averaged method. Comparison between the results obtained by the single and double averaged method is presented. The single-averaged model is more realistic, since it does not eliminate the term due to the equatorial ellipticity of the planetary satellite as done by the double-averaged problem. Considering the single-averaged method, we found unstable frozen orbits where the satellite does not impact with the surface of Europa for at least 200 days. We present an approach using the unaveraged disturbing potential to analyze the effects of these terms in the amplitude of the eccentricity.     

The problems of simulation of planetary systems accumulation processes

Evolution of a system consisting of a great number of bodies that are gravitationally interacting and aggregating in contacts is considered. Body motions take place in the gravitational field of a central massive body (Sun or planet) in the same plane and at the initial time of system evolution orbits of all bodies are circular. It is shown that during evolution of the protoplanetary cloud, ring zones of matter rarefaction and condensation develop. Development of the condensation zones leads to the formation of planets, the most part of which acquire a direct rotation about their axes. In the case under consideration, approximate agreement between the law of planetary distances and that of geometric progression takes place as it is observed in planetary and satellite systems. The formation of the terrestrial planets and Jovian planets has been simulated. The principal numerical results have been obtained through digital simulation of planetary accumulation.     

Symplectic integration of space debris motion considering several Earth’s shadowing models

In this work, we present a symplectic integration scheme to numerically compute space debris motion. Such an integrator is particularly suitable to obtain reliable trajectories of objects lying on high orbits, especially geostationary ones. Indeed, it has already been demonstrated that such objects could stay there for hundreds of years. Our model takes into account the Earth’s gravitational potential, luni-solar and planetary gravitational perturbations and direct solar radiation pressure. Based on the analysis of the energy conservation and on a comparison with a high order non-symplectic integrator, we show that our algorithm allows us to use large time steps and keep accurate results. We also propose an innovative method to model Earth’s shadow crossings by means of a smooth shadow function. In the particular framework of symplectic integration, such a function needs to be included analytically in the equations of motion in order to prevent numerical drifts of the energy. For the sake of completeness, both cylindrical shadows and penumbra transitions models are considered. We show that both models are not equivalent and that big discrepancies actually appear between associated orbits, especially for high area-to-mass ratios.     

Entry of dust particles into planetary magnetospheres

While interplanetary dust constitutes a primary source of cosmic particulate matter in planetary magnetospheres, the debris produced by its impact with small satellites and ring material provides an important secondary source. Internal processes, such as volcanic activity, particularly in the smaller satellites, could result in a third source. In the case of the terrestrial magnetosphere there are also artificial (internal) sources: 1–10μ sized A?₂O₃ particles injected by solid rocket mortar burns between near earth and geosynchronous orbit constitute one such source, while the fragments of larger bodies (artificial satellites) due to explosions (e.g., “killer satellites”) and collisions constitute another. Finally, if we include the purely induced cometary magnetosphere among planetary magnetospheres, the injection of cometary dust into it due to entrainment by the outflowing gases constitutes another source.As a result of being immersed in a radiative and plasma environment these dust grains get electrically charged up to some potential (positive or negative). Particularly in those regions where the magnetospheric plasma is hot and dense and their own spatial density is low, the dust grains could get charged to numerically large negative potentials.While this charging may have physical consequences for the larger grains, such as electrostatic erosion (“chipping”) and disruption, it also can effect the dynamics of the smaller grains. Indeed, the small but finite capacitance of these grains, which leads to a phase lag in the gyrophase oscillation of the grain potential, could even lead to the permanent magneto-gravitational capture of interplanetary grains within planetary magnetospheres in certain situations. Here we will review the sources of dust in planetary magnetospheres and discuss their physics and their dynamics under the combined action of both planetary gravitational and magnetospheric electromagnetic forces.     

Planetary protection considerations for MarsNet and Mars sample return missions.

The ESA MarsNet mission proposal consists most probably of a trio of Mars landers. These landers each contain a variety of scientific equipment. The network of stations demands for a definition of its planetary protection requirements. With respect to the MarsNet mission only forward contamination problems will be considered. Future involvement of European efforts in planetary exploration including sample returns will also raise the problem of back contamination. A tradeoff study for the overall scientific benefit with respect to the approximative cost is necessary. Planetary protection guide-lines will be proposed by an interdisciplinary and international board of experts working in the fields of both biology and planetary science. These guide-lines will have to be flexible in order to be modified with respect to new research results, e.g. on adaptation of microorganisms to extreme (space) conditions. Experiments on the survival of microorganisms at conditions of simulated Mars surface and subsurface will have to be conducted in order to obtain a baseline data collection as a reference standard for future guide-lines.     

Mass loading of planetary atmospheres by rocky satellites—II: Transport and enhanced lifetimes of satellite ejecta in planetary magnetospheres

Extensive studies have been conducted concerning individual mass, temporal and positional distribution of submicron rocky ejecta existing in the satellite-planetary gravitational sphere of influence. The transit time of the major portion of the ejecta that is transported from the satellite's gravitational sphere of influence to the planetary magnetopause is about one week and represents a mass loading pulse occurring each satellite orbit. The mass-flux distributions of lunar ejecta at the surface of the magnetopause for a complete lunar orbit are presented. Spatial mass densities of lunar ejecta in specific zones of the magnetosphere provide a means to compare sporadic interplanetary dust spatial mass densities in the same zones.     

Effects of electrostatic forces on the vertical structure of planetary rings

We have calculated the vertical structure of planetary dust rings as it results from a balance between an electrostatic force on the dust grains and the vertical component of the gravitational force from the central planet. The electrostatic force results from the charging of the dust grains by the ambient plasma and a large scale electric field due to a shielding electric field and the resulting vertical dust distribution are strongly dependent on dust size, dust and plasma density, plasma temperature and plasma ion type. The dust density distribution has a different dependence on these parameters in tenuous and in dense dust rings. We solve the relevant equations numerically and also by linearization in the limiting cases of tenuous or dense rings. Our results indicate that the effects treated in this paper may be important in both Jupiter's and Saturn's rings.     

Gravitational wave detection on the Moon and the moons of Mars

The Moon and the moons of Mars should be extremely quiet seismically and could therefore become sensitive gravitational wave detectors, if instrumented properly. Highly sensitive displacement sensors could be deployed on these planetary bodies to monitor the motion induced by gravitational waves. A superconducting displacement sensor with a 10-kg test mass cooled to 2 K will have an intrinsic instrument noise of 10⁻¹⁶ m Hz^−1/2. These sensors could be tuned to the lowest two quadrupole modes of the body or operated as a wideband detector below its fundamental mode. An interesting frequency range is 0.1–1 Hz, which will be missed by both the ground detectors on the Earth and LISA and would be the best window for searching for stochastic background gravitational waves. Phobos and Deimos have their lowest quadrupole modes at 0.2–0.3 Hz and could offer a sensitivity h_min ? 10⁻²² Hz^−1/2 within their resonance peaks, which is within two orders of magnitude from the goal of the Big Bang Observer (BBO). The lunar and Martian moon detectors would detect many interesting foreground sources in a new frequency window and could serve as a valuable precursor for BBO.     

强不规则天体引力场中的动力学研究进展

小行星探测与彗星探测是深空探测的重要方面。一般来说,小行星和彗星因质量都不足以使得万有引力克服应力达到流体静力学平衡,而具有强不规则的外形。研究强不规则天体引力场中的动力学行为及其内在机制,是探测器被不规则天体捕获并对其形成近距离探测轨道的基础。从引力场模型和动力学行为两个方面综述了强不规则天体引力场中动力学的研究进展,在引力场模型的研究方面介绍了强不规则天体引力场建模的球谐函数摄动展开模型、简单特殊体模型及多面体模型的研究现状,在动力学机制的研究方面介绍了强不规则天体引力场中的周期轨道和拟周期轨道、平衡点、流形、分岔与共振以及混沌运动的研究现状,指出了这些方面研究的重点与难点。分析了强不规则体引力场中动力学的研究趋势。     

Collinear equilibrium points in the elliptic R3BP with oblateness and radiation

We have studied the positions and stability of the collinear equilibrium points, L_1,2,3 of an infinitesimal body in the elliptic restricted three-body problem (ER3BP) when both primaries of the system are luminous and oblate spheroids moving in elliptic orbits around their common center of mass. We observe that their positions are affected by the radiation pressure forces, oblateness and the eccentricity of the orbits, but the stability character remains unchanged and are unstable. The effects of the parameters involved on the collinear points, in particular for the binary systems Achird, Luyten 726-8, Kruger 60, Alpha Centauri AB and Xi Bootis, and their stability in general have been investigated numerically using the analytical results obtained.     

Radiation chemistry of ices of planetological interest at low temperature.

We present the results of recent experiments on some physico-chemical effects induced by fast ion colliding with solids of relevance for the physics of planetary objects. The production of molecular solids, polymer-like materials and amorphous carbon by irradiation of frozen hydrocarbons and pentacene is discussed. We also report on a set of experimental results obtained irradiating methanol and water-methanol mixtures. Because of bombardment different species form. The techniques used for the analysis are "in situ" infrared (IR) and Raman spectroscopy. The experimental results are finally discussed in the light of their relevance for planetary physics.     

Asteroid interception and disruption for terminal planetary defense

Modern techniques for planetary defense from comets and asteroids involve the deflection of the bolide via kinetic, gravitational, ablative, or radiative means. While potentially effective, none of these methods are capable of operating in a terminal interdiction mode wherethe threat is discovered with little time prior to impact. We present a practical and effective method for planetary defense which enables extremely short interdiction time scales, but can also operate within longer time scales and can be effective for extremely large threats. Called PI (“Pulverize It”), the method makes use of an array of hypervelocity penetrators which uses the kinetic energy of the asteroid or comet to disrupt it. In the terminal interdiction mode, the fragments of maximum $～10$ m diameter disperse laterally as they continue towards the Earth, and then enter the Earth’s atmosphere where they burn up as a series of airburst events which spatially and temporally de-correlate the energy of the original parent bolide for any arbitrary observer on the ground in the form of acoustical shockwaves and optical pulses. We show that terminal interdiction modes ranging from 2 minutes prior to impact for 20-meter class bolides (such as the Chelyabinsk asteroid), 1 day prior to impact for 100 m-class asteroids, 10 days prior to impact for Apophis-class asteroids ($～370$ m), and even 60 days prior to impact for 1 km-class threats are all possible, though longer warning times are always preferred. Using only technologies readily available today, the PI method allows for a cost-effective and practical roadmap towards robust planetary defense capability.     

Detecting gravitational waves from accreting neutron stars

The gravitational waves emitted by neutron stars carry unique information about their structure and composition. Direct detection of these gravitational waves, however, is a formidable technical challenge. In a recent study we quantified the hurdles facing searches for gravitational waves from the known accreting neutron stars, given the level of uncertainty that exists regarding spin and orbital parameters. In this paper we reflect on our conclusions, and issue an open challenge to the theoretical community to consider how searches should be designed to yield the most astrophysically interesting upper limits. With this in mind we examine some more optimistic emission scenarios involving spin-down, and show that there are technically feasible searches, particularly for the accreting millisecond pulsars, that might place meaningful constraints on torque mechanisms. We finish with a brief discussion of prospects for indirect detection.     

Stability and mass flow in Saturn's rings

Besides gravitational effects, interesting electrodynamical processes could also take place in the vicinity of the rings of Saturn. In part, this is because of the electrostatic charging of the ring particles by the magnetospheric and ionospheric plasma, and in part, the generation of impact plasma by meteoroid bombardment at the ring plane could lead to strong coupling between the rings and the ionosphere via a variety of current systems. The mass transport and angular momentum transfer in association with the ring-ionosphere coupling may cause quite large changes in the ring configuration over the age of the solar system. The presence of the sharp boundary between the B and the C rings perhaps is a good example. To highlight these new developments, we describe several of the electrodynamical mechanisms (with emphasis on their corresponding electric fields and current systems) which have been postulated to be of importance in determining the mass transport of the ring system. Further points are made that, besides mass exchange between the rings and the planetary atmosphere, the mass injection from the rings could also have significant effect on the mass and energy budget of the magnetosphere, maintenance of the E ring, the Titan hydrogen torus as well as aeronomic process in the upper atmosphere of Titan.     

Dynamical characterization of the last prolonged solar minima

The planetary hypothesis of the solar cycle is an old idea in which the gravitational influence of the planets has a non-negligible effect on the causes of the solar magnetic cycle. The advance of this hypothesis is based on phenomenological correlations between dynamical parameters of the Sun’s movement around the barycentre of the Solar System and sunspots time series; and more especially, identifying relationships linking solar barycentric dynamics with prolonged minima (especially Grand Minima events). However, at present there is no clear physical mechanism relating these phenomena. The possible celestial influence on solar cycle modulation is of great importance not only in solar physics but also in Earth sciences, because prolonged solar minima have associated important climatic and telluric variations, in particular, during the Maunder and Dalton Minimum. In this work we looked for a possible causal link in relation with solar barycentric dynamics and prolonged minima events. We searched for particular changes in the Sun’s acceleration and concentrated on long-term variations of the solar cycle. We show how the orbital angular momentum of the Sun evolves and how the inclination of the solar barycentric orbit varies during the epochs of orbital retrogressions. In particular, at these moments, the radial component of the Sun’s acceleration (i.e., in the barycentre-Sun direction) had an exceptional magnitude. These radial impulses occurred at the very beginning of the Maunder Minimum, during the Dalton Minimum and also at the maximum of cycle 22 before the present extended minimum. We also found a strong correlation between the planetary torque and the observed sunspots international number around that maximum. We apply our results in a novel theory of Sun–planets interaction that it is sensitive to Sun barycentric dynamics and found a very important effect on the Sun’s capability of storing hypothetical reservoirs of potential energy that could be released by internal flows and might be related to the solar cycle. This process begins about 40 years before the solar angular momentum inversions, i.e., before Maunder Minimum, Dalton Minimum, and before the present extended minimum. Our conclusions suggest a dynamical characterization of peculiar prolonged solar minima. We discuss the possible implications of these results for the solar cycle including the present extended minimum.     

Alternative transfer to the Earth–Moon Lagrangian points L4 and L5 using lunar gravity assist

Lagrangian points L4 and L5 lie at 60° ahead of and behind the Moon in its orbit with respect to the Earth. Each one of them is a third point of an equilateral triangle with the base of the line defined by those two bodies. These Lagrangian points are stable for the Earth–Moon mass ratio. As so, these Lagrangian points represent remarkable positions to host astronomical observatories or space stations. However, this same distance characteristic may be a challenge for periodic servicing mission. This paper studies elliptic trajectories from an Earth circular parking orbit to reach the Moon’s sphere of influence and apply a swing-by maneuver in order to re-direct the path of a spacecraft to a vicinity of the Lagrangian points L4 and L5. Once the geocentric transfer orbit and the initial impulsive thrust have been determined, the goal is to establish the angle at which the geocentric trajectory crosses the lunar sphere of influence in such a way that when the spacecraft leaves the Moon’s gravitational field, its trajectory and velocity with respect to the Earth change in order to the spacecraft arrives at L4 and L5. In this work, the planar Circular Restricted Three Body Problem approximation is used and in order to avoid solving a two boundary problem, the patched-conic approximation is considered.     

Robotic Instrument for Grinding Rocks Into Thin Sections (GRITS)

We have developed a rock grinding and polishing mechanism for in situ planetary exploration based on abrasive disks, called Grinding Rocks Into Thin Sections (GRITS). Performance characteristics and design considerations of GRITS are presented. GRITS was developed as part of a broader effort to develop an in situ automated rock thin section (ISARTS) instrument. The objective of IS-ARTS was to develop an instrument capable of producing petrographic rock thin sections on a planetary science spacecraft. GRITS may also be useful to other planetary science missions with in situ instruments in which rock surface preparation are necessary.     

Integral inversion of GRAIL inter-satellite gravitational accelerations for regional recovery of the lunar gravitational field

We present an integral-based approach for high-resolution regional recovery of the gravitational field in this article. We derive rigorous remove-compute-restore integral estimators relating the line-of-sight gravitational acceleration to an arbitrary order radial derivative of the gravitational potential. The integral estimators are composed of three terms, i.e., the truncated integration, the low-frequency line-of-sight gravitational acceleration, and the high-frequency truncation error (effect of the distant zones). We test the accuracy of the integral transformations and of the integral estimators in a closed-loop simulation over the Montes Jura region on the nearside of the Moon. In this way, we determine optimal sizes of integration radii and grid discretisation. In addition, we investigate the performance of the regional integral inversion with synthetic and realistic GRAIL observations. We demonstrate that the regional inversion results of the disturbing gravitational potential and its first order radial derivative in the Montes Jura mountain range are less contaminated by high-frequency noise than the global spherical harmonic models.     

Space telescope studies of Jupiter and Saturn simulated by Voyager observations

Space Telescope (ST) observations of Jupiter and Saturn will offer a unique opportunity for monitoring their changing meteorological characteristics. They will provide higher spatial and temporal resolution for composition and vertical structure studies than have been available to date. We have simulated the planetary camera observations of Jupiter and Saturn by Voyager images of the appropriate spatial scale. With this data set we have investigated the meteorological properties of these atmospheres which can be studied at these scales. In addition we have considered the advances obtainable with the high resolution spectrometer on ST compared with observations from ground-based and other Earth-orbiting satellites. These studies will provide insight into the scientific gain and possible problems in the use of ST for planetary studies.