Averaging on the motion of a fast revolving body. Application to the stability study of a planetary system

Exploring the global dynamics of a planetary system involves computing integrations for an entire subset of its parameter space. This becomes time-consuming in presence of a planet close to the central star, and in practice this planet will be very often omitted. We derive for this problem an averaged Hamiltonian and the associated equations of motion that allow us to include the average interaction of the fast planet. We demonstrate the application of these equations in the case of the μ Arae system where the ratio of the two fastest periods exceeds 30. In this case, the effect of the inner planet is limited because the planet’s mass is one order of magnitude below the other planetary masses. When the inner planet is massive, considering its averaged interaction with the rest of the system becomes even more crucial.     

Analytical theory of a lunar artificial satellite with third body perturbations

We present here the first numerical results of our analytical theory of an artificial satellite of the Moon. The perturbation method used is the Lie Transform for averaging the Hamiltonian of the problem, in canonical variables: short-period terms (linked to l, the mean anomaly) are eliminated first. We achieved a quite complete averaged model with the main four perturbations, which are: the synchronous rotation of the Moon (rate ), the oblateness J ₂ of the Moon, the triaxiality C ₂₂ of the Moon ( ) and the major third body effect of the Earth (ELP2000). The solution is developed in powers of small factors linked to these perturbations up to second-order; the initial perturbations being sorted ( is first-order while the others are second-order). The results are obtained in a closed form, without any series developments in eccentricity nor inclination, so the solution apply for a wide range of values. Numerical integrations are performed in order to validate our analytical theory. The effect of each perturbation is presented progressively and separately as far as possible, in order to achieve a better understanding of the underlying mechanisms. We also highlight the important fact that it is necessary to adapt the initial conditions from averaged to osculating values in order to validate our averaged model dedicated to mission analysis purposes.     

Chirikov diffusion in the asteroidal three-body resonance (5, −2, −2)

The theory of diffusion in many-dimensional Hamiltonian system is applied to asteroidal dynamics. The general formulation developed by Chirikov is applied to the Nesvorny-Morbidelli analytic model of three-body (three-orbit) mean-motion resonances (Jupiter-Saturn-asteroid). In particular, we investigate the diffusion along and across the separatrices of the (5, −2, −2) resonance of the (490) Veritas asteroidal family and their relationship to diffusion in semi-major axis and eccentricity. The estimations of diffusion were obtained using the Melnikov integral, a Hadjidemetriou-type sympletic map and numerical integrations for times up to 10⁸ years.     

Homographic solutions in the planar <Emphasis Type="Italic">n</Emphasis> + 1 body problem with quasi-homogeneous potentials

We prove that for generalized forces which are function of the mutual distance, the ring n + 1 configuration is a central configuration. Besides, we show that it is a homographic solution. We apply the above results to quasi-homogeneous potentials.     

The Lambert <Emphasis Type="Italic">W</Emphasis> function and solutions to Kepler’s equation

This paper presents a method for the truncation of infinite Fourier–Bessel representations for functions requiring a solution to Kepler’s equation. Use is made of the Lambert W function to solve for the desired index that bounds the remainder terms of the series, within the prescribed tolerance. The enforcement of a maximum on the number of Bessel functions is also useful in truncating the Bessel functions themselves, resulting in an analytical representation of the solution to a desired tolerance, without the use of infinite series.     

Fourier techniques for an analysis of eclipsing binary light curves. VIa

The aim of the present paper is to deduce relations between the integral transformsA _2m, B_2m,andF _1,2 of the light curves of eclipsing binary systems. The integral transformsA _2m, B_2m,andF _1,2 have been related to one another by means of finite or rapidly converging infinite summations obtained by integrations of the series expansions of trigonometric functions.     

Franck-Condon factors and <Emphasis Type="Italic">r</Emphasis>-centroids for certain band systems of astrophysical molecules SrF and ScF

The Franck-Condon factors and r-centroids, which are very closely related to vibrational transition probabilities, have been evaluated by the more reliable numerical integration procedure for the bands of B ²∑⁺ − X ²∑⁺, F ²∑⁺ − X ²∑⁺ systems of SrF and C ¹∑⁺ − X ¹∑⁺, G ¹Π − X ¹∑⁺ systems of ScF molecules of astrophysical interest, using a suitable potential.     

Numerical Integration of Regular Equations of the Planar Motion of Terrestrial Bodies

The parameters of L matrices are applied to the numerical integration of regular equations describing the motion of minor bodies in the Solar System. The problem of the optimal choice of the regularizing change of variables is formulated in the context of the numerical integration of the equations of motion using the Runge–Kutta–Fehlberg method. Arbitrary perturbations are taken into account. This problem is completely solved in the case of planar motion. The solution of the optimization problem reduces the amount of computations needed to determine the vector of perturbing accelerations. Results of numerical integrations are given.     

Seasonal monsoon forecast for the years 1987 and 1988

The results of several seasonal integrations with an atmospheric global circulation model with prescribed “perfect” sea surface temperatures are presented. These experiments illustrate the results of seasonal simulations for the years 1987 and 1988. These were a dry and a wet monsson year, respectively, when compared to the average. The integrations cover the period from June 1 through August 31 for both years and were carried out at two horizontal resolutions, T42 and T106, of a global model containing two different parameterizations of surface hydrology. The seasonal differences of the motion fields, divergent circulations and rainfall distributions for these respective experiments are compared with the corresponding observed fields.The sensitivity of seasonal simulations to the initial state is explored with integrations starting on two successive dates. In these experiments we diagnose differences of the simulated time mean states from residue free budgets of the complete vorticity equation.     

Eccentricity Generation in Hierarchical Triple Systems with Non-coplanar and Initially Circular Orbits

In a previous paper, we developed a technique for estimating the inner eccentricity in coplanar hierarchical triple systems on initially circular orbits, with comparable masses and with well-separated components, based on an expansion of the rate of change of the Runge-Lenz vector. Now, the same technique is extended to non-coplanar orbits. However, it can only be applied to systems with I ₀ < 39.23° or I ₀ > 140.77°, where I is the inclination of the two orbits, because of complications arising from the so-called ‘Kozai effect’. The theoretical model is tested against results from numerical integrations of the full equations of motion. This revised version was published online in July 2006 with corrections to the Cover Date.     

Non-Existence Of <Emphasis Type="Italic">N</Emphasis>-Dimensional Static Plane Symmetric Solutions In Bimetric Relativity Theory

A problem of static plane symmetric metric in the perfect fluid, the mesonic massive scalar field and in their coupling is studied in Rosen’s (1973) bimetric theory of relativity. It was found that the matter field like either perfect fluid or mesonic massive scalar field or their coupling does not survive in bimetric theory of gravitation when the space–time is governed by n-dimensional static plane symmetric metric.     

The N-gon is not a local minimum of U 2 I for N > 7

A theorem of Palmore's concerning coplanar central configurations of equal mass bodies was shown to be false for all even N 6 by Slaminka and Woerner. Using a variation of that argument I prove that Palmore's Theorem is false for all N 6.Northwestern University     

Comparison between Different Models of Galactic Tidal Effects on Cometary Orbits

Different models of the action of the galactic tide are compared. Each model is a substitute for direct numerical integrations allowing a drastic decrease of the computation time. The models are built using two different techniques, (i) averaging of the fast variable (the mean anomaly) over one cometary period and (ii) fixing the comet in its aphelion direction. Moreover, we consider two different formalisms (Lagrangian and Hamiltonian) and also two different sets of variables. As expected, we find that the model results are independent of the formalism and the set of variables considered, and are highly accurate, whereas mathematical technique leads to poor results. In order to further reduce the computation time, mappings are built from the development of the solution of the models. We show that for these mappings, the set of variables giving the most accurate results is strongly dependent on the cometary eccentricity, e, and semimajor axis, a.     

Fitting a line to a sine

In the numerical integrations of the outer planets, when the inner planets are not considered, the initial conditions are determined by least square adjustment with a complete ephemeris. In the present work, we compute the resulting uncertainty on the initial conditions by the analysis of the very simple problem which consists in fitting by least squares a straight line to a sine curve. Explicit formulae for the computation of the error due to the fitting process are given. In particular, fitting over an integer number of periods is very close to the worst case. The determination of the initial conditions in the LONGSTOP project (Milani,et al., 1987) is analysed as an example. It is shown that the uncertainty on the main frequencies of the secular system due to the determination of the initial conditions is of the order of 0.01 arcsec/year.     

The ideal resonance problem a comparison of two formal solutions I

Garfinkel's solution of the Ideal Resonance problem derived from a Bohlin-von Zeipel procedure, and Jupp's solution, using Poincaré's action and angle variables and an application of Lie series expansions, are compared. Two specific Hamiltonians are chosen for the comparison and both solutions are compared with the numerical solutions obtained from direct integrations of the equations of motion. It is found that in deep resonance the second-mentioned solution is generally more accurate, while in the classical limit the first solution gives excellent agreement with the numerical integrations.This article represents a summary of a much more extensive programme of research, the complete results of which will be published in a future article.     

Disintegration process of hierarchical triple systems. I. Small-mass planet orbiting equal-mass binary

The stability of hierarchical triple system is studied in the case of an extrasolar planet or a brown dwarf orbiting a pair of main sequence stars. The evolution of triple system is well modelled by random walk (RW) diffusion, particularly in the cases where the third body is small and tracing an orbit with a large eccentricity. A RW model neglects the fact that there are many periodic orbits accompanied by stability islands, and hence inherently overestimates the instability of the system. The present work is motivated by the hope to clarify how far the RW model is applicable. Escape time and the surface section technique are used to analyse the outcome of numerical integrations. The analysis shows that the RW-like model explains escape of the third body if the initial configuration is directly outside of the KAM tori. A small gap exists in (q ₂/a ₁, e ₂)-plane between locations of the stability limit curves based on our numerical study and on RW-model (the former is shifted by –1.4 in q ₂/a ₁ direction from the latter).     

Orbital and temporal distributions of meteorites originating in the asteroid belt

Abstract— The recent discovery of the importance of Sun-grazing phenomena dramatically changed our understanding of the dynamics of objects emerging from the asteroid belt via resonant phenomena. The typical lifetimes of such objects are now expected to be <10 Ma, thus demanding a reassessment of our general picture of the meteorite delivery process. By analysing direct numerical integrations of ～2000 test particles beginning in the v₆, 3:1, and 5:2 resonances in the main belt, we have reexamined the orbital and temporal distribution of meteoroids that journey to Earth. Comparing the results with fireball data, we find that the orbital distribution of Earth-impacting chondrites is consistent with a steady-state injection of meteoroids into the 3:1 and v₆, resonances. Because this is the most complete and unbiased data set concerning Earth-impacting meteoroids, the agreement leads us to believe that our model is accurate. The simulations predict a P.M. fall ratio for chondrites ～14% lower than the observed value of ～68%, which argues for a moderate bias being present in this statistic. Most interestingly, the typical meteorite transfer times predicted by our models are several factors lower than the typical chondrite exposure ages, which implies that these meteorites acquired most of their exposure in the main belt before entering the resonances. We discuss some processes that would allow such preexposure. The case of achondrites and iron meteorites is also briefly discussed.     

Critical accretion flow of gas onto compact stars

The critical accretion flow of gas onto compact stars with mass of 0.6M is investigated by numerical integrations of the time-dependent hydrodynamic equations in the sphericallysymmetric and optically thick case. For the compact stars surrounded by such a dense cloud of gas, the radiation pressure force decelerates the infall gas significantly and free fall regime of the gas is not at all attained. This results in incident low velocities at the standing shock front close to the stellar surface, low temperatures of the gas around the compact stars, and no X-ray in white dwarfs but soft X-rays in neutron stars, respectively. Some applications of the results to the X-ray sources are discussed.     

Cartography of the <Emphasis Type="Italic">b</Emphasis>-plane of a close encounter I: semimajor axes of post-encounter orbits

Close planetary encounters play an important role in the evolution of the orbits of small Solar system bodies and are usually studied with the help of numerical integrations. Here we study close encounters in the framework of an analytic theory, focusing on the so-called b-plane, which is the plane centred on the planet and perpendicular to the planetocentric velocity at infinity of the small body. As shown in previous papers, it is possible to identify the initial conditions on the b-plane that lead to post-encounter orbits of given semimajor axis. In this paper we exploit analytical relationships between b-plane coordinates and pre-encounter orbital elements and compute the probability of transition to these post-encounter states, and numerically check the validity of the analytic approach.     

Orbital evolution of giant comet-like objects

We investigated by numerical integrations the long-term orbital evolution of four giant comets or comet-like objects. They are Chiron, P/Schwassmann-Wachmann 1 (SW1), Hidalgo, and 1992AD (5145), and their orbits were traced for 100–200 thousand years (kyr) toward both the past and the future. For each object, 13 orbits were calculated, one for the nominal orbital elements and other 12 with slightly modified elements based on the rms residual of the orbit determination and on the number of observations. As past studies indicate, their orbital evolution is found to be very chaotic, and thus can be described only in terms of probability. Plots of the semi-major axis (a) and perihelion distance (q) of the objects treated here seem to cross each other frequently, suggesting a possibility of their common evolutionary paths. About a half of all the calculated orbits showedq- ora-decreasing evolution. This indicates that, at least on the time scale in question, the giant comet-like objects are possibly on a dynamical track that can lead to capture from the outer solar system. We could hardly find the orbits with perihelia far outside the orbit of Saturn (q>15 AU). This is perhaps because the evolution of the orbits beyond Saturn is so slow that substantial orbital changes do not take place within 100–200 kyr.