Numerical experiments in the 3/1 andv 6 overlapping resonance region

Numerical experiments of fictitious small bodies with initial eccentricities e=0.1 have been performed in the overlapping region of the 3/1 mean motion resonance and of thev ₆ secular resonance 2.48a2.52AU for different values of the initial inclination 16°i20°. An analysis for thev ₆ secular resonance shows that the topology is different from the one found outside the overlapping region: the critical argument for thev ₆ resonance in the overlapping region rotates in opposite direction as compared to the purev ₆ region. In the 3/1 resonance region the secular resonancev ₅ is dominant, and some secondary secular resonances asv ₆ –v ₁₆ andv ₅ +v ₆ are present.     

Numerical results to the Sitnikov-problem

We present numerical results of the so-called Sitnikov-problem, a special case of the three-dimensional elliptic restricted three-body problem. Here the two primaries have equal masses and the third body moves perpendicular to the plane of the primaries' orbit through their barycenter. The circular problem is integrable through elliptic integrals; the elliptic case offers a surprisingly great variety of motions which are until now not very well known. Very interesting work was done by J. Moser in connection with the original Sitnikov-paper itself, but the results are only valid for special types of orbits. As the perturbation approach needs to have small parameters in the system we took in our experiments as initial conditions for the work moderate eccentricities for the primaries' orbit (0.33e _primaries 0.66) and also a range of initial conditions for the distance of the 3 ^rd body (= the planet) from very close to the primaries orbital plane of motion up to distance 2 times the semi-major axes of their orbit. To visualize the complexity of motions we present some special orbits and show also the development of Poincaré surfaces of section with the eccentricity as a parameter. Finally a table shows the structure of phase space for these moderately chosen eccentricities.     

Air drag effect on the motion of an artificial Earth satellite

Two methods have been used to compute and compare the perturbations in perigee distance for an artificial Earth satellite. The two methods have used different air density models. The first (Helali, 1987) used the TD model, formulated by Sehnel (1986a), which contains terms that describe all the principal changes of the thermospheric density due to solar activity, geomagnetic activity, and the height. The second method (Davis, 1963) used a model of the density which takes into account the rotation of the atmosphere, the bulging atmosphere and the height. For different values of eccentricities from 0.001 to 0.05 we computed the perturbations P _r in the perigee distance at different heights from 200 to 350 km for both methods. The results show a good agreement for the computed values of P _r for different values of e (0 < e 0.02) in both methods at perigee heights from 250 to 350 km. Meanwhile, for perigee heights smaller than about 250 km we found a maximum difference in P _r amounting to 20 metres/revolution for e = 0.005 and 0.01.     

Chaotic behaviour of trajectories for the asteroidal resonances

A systematic study of the main asteroidal resonances of the third and fourth order is performed using mapping techniques. For each resonance one-parameter family of surfaces of section is presented together with a simple energy graph which helps to understand and predict the changes in the surfaces of section within the family. As the truncated Hamiltonian for the planar, elliptic, restricted three-body problem is used for the mapping, the method is expected to fail for high eccentricities. We compared, therefore, the surfaces of section with trajectories calculated by symplectic integrators of the fourth and six order employing the full Hamiltonian. We found a good agreement for small eccentricities but differences for the higher eccentricities (e 0.3).     

Some peculiarities in the Asteroidal Belt

The analysis of the fine structure of the Asteroidal Belt evidenciates a group of asteroids next to the resonance 4/9 with Jupiter. In this group and in other groups associated to the Hirayama families there are indications that their orbital parameters can be represented by quantum numbers as defined here and in two of our previous works. Together with this the distribution of the eccentricities and inclinations of the orbital planes of short period comets and diverse type of asteroids indicates that they can be classified as objects with e > sin i and objects with e > sin i with a limit e = sin i which determinates geometrical properties of the orbits related with discrete states in the solar system. This study lets open the possibility of following studies in order to confirm the quantum characteristics of the Asteroidal Belt being these characteristics common to all the solar system and depending of the same fundamental constant of action per mass unit H ₀ = 1/2 ₀ × T ₀ (potential × time) because only a small part of all the available data in the Asteroid Belt is used here.     

Stability criteria in many-body systems

An expansion of the force function ofn-body dynamical systems, where the equations of motion are expressed in the Jacobian coordinate system, is shown to give rise naturally to a set of (n–1) (n–2) dimensionless parameters ^ki _li {i = 2,...,n;k = 2,...,i – 1 (i 3);l =i + 1,...,n (i n – 1)}, representative of the size of the disturbances on the Keplerian orbits of the various bodies. The expansion is particularized to the casen=3 which involves the consideration of only two parameters ²³ and ₃₂. Further, the work of Szebehely and Zare (1977) is reviewed briefly with reference to a sufficient condition for the stability of corotational coplanar three-body systems, in which two of the bodies form a binary system. This condition is sufficient in the sense that it precludes any possibility of an exchange of bodies, i.e. Hill type stability, however, it is not a necessary condition. These two approaches are then combined to yield regions of stability or instability in terms of the parameters ²³ and ₃₂ for any system of given masses and orbital characteristics (neglecting eccentricities and inclinations) with the following result: that there is a readily applicable rule to assess the likelihood of stability or instability of any given triple system in terms of ²³ and ₃₂.Treating a system ofn bodies as a set of disturbed three-body systems we use existing data from the solar system, known triple systems and numerical experiments in the many-body problem to plot a large number of triple systems in the ²³, ₃₂ plane and show the results agree well with the ²³, ₃₂ analysis above (eccentricities and inclinations as appropriate to most real systems being negligible). We further deal briefly with the extension of the criteria to many-body systems wheren>4, and discuss several interesting cases of dynamical systems.     

Formation of terrestrial planets in a dissipating gas disk with Jupiter and Saturn

We have performed N-body simulations on final accretion stage of terrestrial planets, including the eccentricity and inclination damping effect due to tidal interaction with a gas disk. We investigated the dependence on a depletion time scale of the disk, and the effect of secular perturbations by Jupiter and Saturn. In the final stage, terrestrial planets are formed through coagulation of protoplanets of about the size of Mars. They would collide and grow in a decaying gas disk. Kominami and Ida [Icarus 157 (2002) 43-56] showed that it is plausible that Earth-sized, low-eccentricity planets are formed in a mostly depleted gas disk. In this paper, we investigate the formation of planets in a decaying gas disk with various depletion time scales, assuming disk surface density of gas component decays exponentially with time scale of τ_gas. Fifteen protoplanets with are initially distributed in the terrestrial planet regions. We found that Earth-sized planets with low eccentricities are formed, independent of initial gas surface density, when the condition (τ_cross+τ_growth)/2?τ_gas?τ_cross is satisfied, where τ_cross is the time scale for initial protoplanets to start orbit crossing in a gas-free case and τ_growth is the time scale for Earth-sized planets to accrete during the orbit crossing stage. In the cases satisfying the above condition, the final masses and eccentricities of the largest planets are consistent with those of Earth and Venus. However, four or five protoplanets with the initial mass remain. In the final stage of terrestrial planetary formation, it is likely that Jupiter and Saturn have already been formed. When Jupiter and Saturn are included, their secular perturbations pump up eccentricities of protoplanets and tend to reduce the number of final planets in the terrestrial planet regions. However, we found that the reduction is not significant. The perturbations also shorten τ_cross. If the eccentricities of Jupiter and Saturn are comparable to or larger than present values (∼0.05), τ_cross become too short to satisfy the above condition. As a result, eccentricities of the planets cannot be damped to the observed value of Earth and Venus. Hence, for the formation of terrestrial planets, it is preferable that the secular perturbations from Jupiter and Saturn do not have significant effect upon the evolution. Such situation may be reproduced by Jupiter and Saturn not being fully grown, or their eccentricities being smaller than the present values during the terrestrial planets' formation. However, in such cases, we need some other mechanism to eliminate the problem that numerous Mars-sized planets remain uncollided.     

Rotating Stratified Heterogeneous Oblate Spheroid in Newtonian Physics

By solving the Euler hydrodynamical equations we have obtained closed form solutions for the angular velocities and pressures of a three stratified non-confocal heterogeneous oblate spheroid. Limiting and particular solutions cases, such as a spheroid with N layers, a stratified spheroid with the same eccentricities, as well as confocal layered spheroids are also explicitly written down. As an application, we have numerically estimated planet Earth's outer and inner cores' ellipticities to be _o=1/413.318 and _i=1/424.616, respectively. These Earth's ellipticities values are in good agreement with those found in the literature.     

Solution of canonical equations which result from elimination of short-periodic terms in second-order planetary theory

We solve the first order non-linear differential equation and we calculate the two quadratures to which are reduced the canonical differential equations resulting from the elimination of the short period terms in a second order planetary theory carried out through Hori's method and slow Delaunay canonical variables when powers of eccentricities and the sines of semi-inclinations which are >3 are neglected and the eccentricity of the disturbing planet is identically equal to zero. The procedure can be extended to the case when the eccentricity of the disturbing planet is not identically equal to zero. In this latter general case, we calculatedthe two quadratures expressing angular slow Delaunay canonical variable ₁ of the disturbed planet and angular slow Delaunay canonical variable ₂ of the disturbing planet in terms of timet.     

The expansion of negative powers of mutual distance between two planets

In this part we obtain the expression for ^–s by the application of Smart's method, which involves Taylor's theorem for functions of several variables. We neglected terms of power higher than the fourth with respect to eccentricities and tangents of inclinations.     

On the elimination of the critical terms of a first order theory of Jupiter perturbed by Saturn carried out through Hori's method and Poincaré canonical variables

The elimination of the critical terms inside the Hamiltonian of a first order theory of Jupiter perturbed by Saturn is carried out through the Poincaré canonical variables and the Hori's procedure. Powers of the eccentricities and the sines of inclinations which are>3 are neglected. The Poincaré variablesL ₁,H ₁,P ₁, ₁,K ₁,Q ₁ of Jupiter which result from a previous elimination of the short period terms are expressed in terms of the Poincaré canonical variablesL _u ,H _u ,P _u , _u ,Q _u ;u=1, 2; index 1 Jupiter, index 2 Saturn resulting from the elimination of the short period and critical terms. The differential equations inH _u ,P _u ,K _u ,Q _u are solved through the method of Lagrange and the analytical expressions ofL ₁,H ₁,P ₁, ₁,K ₁,Q ₁ as functions of timet are finally obtained.     

A Symplectic Mapping Model for the Study of 2:3 Resonant Trans-Neptunian Motion

A symplectic mapping is constructed for the study of the dynamical evolution of Edgeworth-Kuiper belt objects near the 2:3 mean motion resonance with Neptune. The mapping is six-dimensional and is a good model for the Poincaré map of the real system, that is, the spatial elliptic restricted three-body problem at the 2:3 resonance, with the Sun and Neptune as primaries. The mapping model is based on the averaged Hamiltonian, corrected by a semianalytic method so that it has the basic topological properties of the phase space of the real system both qualitatively and quantitatively. We start with two dimensional motion and then we extend it to three dimensions. Both chaotic and regular motion is observed, depending on the objects' initial inclination and phase. For zero inclination, objects that are phase-protected from close encounters with Neptune show ordered motion even at eccentricities as large as 0.4 and despite being Neptune-crossers. On the other hand, not-phase-protected objects with eccentricities greater than 0.15 follow chaotic motion that leads to sudden jumps in their eccentricity and are removed from the 2:3 resonance, thus becoming short period comets. As inclination increases, chaotic motion becomes more widespread, but phase-protection still exists and, as a result, stable motion appears for eccentricities up to e = 0.3 and inclinations as high as i = 15°, a region where plutinos exist.     

Secular resonances in the primitive solar nebula

We present here a very simple model that could explain the relatively high eccentricities and inclinations observed in the minor planet belt. This model is based upon the sweeping of the secular resonances ₆ and ₁₆ through the belt due to the gravitational effect of the dissipation of a primitive solar nebula. The sweeping of the ₁₆ secular resonance (responsible for the high inclinations) is very sensitive to the density profile of the nebula. For the model to work we need a density profile proportional to ^–k with between 1.0 and 1.5.     

A discussion of Hill's method of secular perturbations and its application to the determination of the zero-rank effects in non-singular vectorial elements of a planetary motion

Knowledge of the perturbations of zero-rank is essential for the understanding of the behavior of a planetary or cometary orbit over a long interval of time. Recent investigations show that these zero-rank perturbations can cause large oscillations in both the shape and position of the orbit. At present we lack a complete analytical theory of these perturbations that can be applied to cases where either the eccentricity or inclination is large or has large oscillations. For this reason we here develop formulas for the numerical integration of the zero-rank effects, using a modified Hill's theory and suitable vectorial elements. The scalar elements of our theory are the two components of Hamilton's vector in a moving ideal reference frame and the three components of Gibb's rotation vector in an inertial system. The integration step can be taken to be several hundred years in the planetary or cometary case, and a few days in the case of a near-Earth space probe. We re-discuss Hill's method in modern symbolism and by applying the vectorial analysis in a pseudo-euclidean spaceM ₃, we obtain a symmetrical computational scheme in terms of traces of dyadics inM ₃. The method is inapplicable for two orbits too close together. In Hill's method the numerical difficulty caused by such proximity appears in the form of a small divisor, whereas in Halphen's method it appears as a slow convergence of a hypergeometric series. Thus, in Hill's method the difficulty can be watched more directly than in Halphen's method. The methods of numerical averaging have, at the present time, certain advantages over purely analytical methods. They can treat a large range of eccentricities and orbital inclinations. They can also treat the free secular oscillations as well as the forced ones, and together with their mutual cross-effects. At the present time, no analytical theory can do this to the full extent.Basic Notations m the mass of the disturbed body - M the mass of the Sun - f the gravitational constant - f(M+m) - r the heliocentric position vector of the disturbed body - r |r| - r ⁰ the unit vector alongr - n ⁰ the unit vector normal tor and lying in the orbital plane of the disturbed body - a the semi-major axis of the orbit of the disturbed body - e the eccentricity of the orbit of the disturbed body - g the mean anomaly of the disturbed body - the eccentric anomaly of the disturbed body - p a(1–e ²) - P ₁ the unit vector directed from the Sun toward the perihelion of the disturbed body - P ₂ the unit vector normal toP ₁ and lying in the orbital plane of the disturbed body - s - the true orbital longitude of the disturbed body, reckoned from the departure point of the ideal system of coordinates - X the true orbital longitude of the perihelion of the disturbed body in the ideal system of coordinates reckoned from the departure point - the angular distance of the ascending node from the departure point - R ₁,R ₂,R ₃ the unit vectors along the axes of the ideal system of coordinates,R ₁ andR ₂ are in the osculating orbital plane of the disturbed body,R ₃ is normal to this plane. The intersection ofR ₁ with the celestial sphere is the departure point - R ₃ P ₁×P ₂ - S ₁,S ₂,S ₃ the initial values ofR ₁,R ₂,R ₃, respectively - q the Gibb's vector. This vector defines the rotation of the orbital plane of the disturbed body from its initial position to the position at the given timet - m the mass of the disturbing body - r the heliocentric position vector of the disturbing body - a the semi-major axis of the orbit of the disturbing body - e the eccentricity of the orbit of the disturbing body - g the mean anomaly of the disturbing body - the eccentric anomaly of the disturbing body - P₁ the unit vector directed from the Sun toward the perihelion of the disturbing body - P₂ the unit vector normal toP₁ and lying in the orbital plane of the disturbing body - A₁ a P₁ - A₂ - |r–r|     

Cyclotron waves in the solar wind

Cyclotron waves in the solar wind near 1 AU with frequencies well below the electron cyclotron frequency and wavelengths much larger than the electron cyclotron radius but less than the proton cyclotron radius are considered. The cyclotron radii are defined from parallel thermal velocity of electron component and proton component with respect to the interplanetary magnetic field. No LH cyclotron waves are found to propagate for _p < 0, where _p 1 –T _p/T _p is the temperature anisotropy of the proton component with respect to the interplanetary magnetic field. The damping or growth of RH cyclotron waves is found to depend on the frequency range and the temperature anisotropy of the proton component. The RH cyclotron waves are damped in the frequency range _r | _p | _p for _p < 0, where _p is the proton cyclotron frequency. RH cyclotron instabilities occur in the frequency range | _p | _p > _r > | _p | _p /(1– _r ) for _p < 0. The marginal state is at _r =| _p | _p .Abstract presented at theInternational Symposium on Solar-Terrestrial, São Paulo, Brazil, 17–22 June, 1974     

Angular momenta in the solar system

A statistical study of the orbital parameters of comets, asteroids and meteor streams shows that the vectors representing their angular momenta per mass unit (or the average angular momentum for meteor streams) are not arbitrarily distributed in the space: They are clustered around determinated values of angles . This synthesizes the eccentricities and inclinations of the orbital planes in a unique parameter adequated for the statistical purposes of the present work being defined by cos = cos (arc sin e) cos i.The discreteness of the obtained distribution N() and its relation with the components of the angular momenta per mass unit is analysed having this distribution common features for objects of different nature and located in different places in the solar potential well. Some hypotheses concerning to these effects are discussed.     

Period distributions in pulsating and binary stars

We analyze the hypothesis of quantization in bands for the angular momenta of binary systems and for the maount of actionA _c in stable and pulsating stars. This parameter isA _c=Mv _eff R _eff, where the effective velocity corresponds to the kinetic energy in the stellar interior and the effective radius corresponds to the potential energyGM ²/R _eff. Analogous parameters can be defined for a pulsating star withm=M where is the rate of the massm participating in the oscillation to the total massM andv _osc,R _osc the effective velocity and oscillation radius.From an elementary dimensional analysis one has thetA _c (energy x time) (period)^1/3 independently ifA _c corresponds to the angular momentum in a binary system, or to the oscillation in a pulsating star or the inner energy and its time-scaleP _eff in a stable star.From evolving stellar models one has that P _effP _eff(solar)1.22 hr a near-invariant for the Main Sequence and for the range of masses 0.6M <M<1.6M .With this one can give scalesn _k=kn ₁ withk integers andn ₁=(P/P ₁)^1/3 withP ₁=P _eff1.22 hr. In these scales proportional toA _c, one sees that the periods in binary and pulsating stars are clustered in discrete unitsn ₁,n ₂,n ₃, etc.This can be seen in pulsating Scuti, Cephei, RR Lyrae, W Virginis, Cephei, semi-regular variables, and Miras and in binary stars as cataclysmic binaries, W Ursa Majoris, Algols, and Lyrae with the corresponding subgroups in all these materials. Phase functions (n _k) in RR Lyrae and Cephei are also associated with discrete levelsn _k.the suggested scenario is that the potential energies and the amounts of actionE _p(t), A_c(t) are indeed time-dependent, but the stars remain more time in determinated most proble states. The Main Sequence itself is an example of this. These most probable states in binary systems, or pulsating or stable stars, must be associated with velocities sub-multiplesc/ _F , given by the velocity of light and the fine structure constant.Additional tests for such a hypothesis are suggested when the sufficient amount of observational data are available. They can made with oscillation velocities in pulsating stars and velocity differences of pairs of galaxies.     

Green's theorem and Green's functions for the steady-state cosmic-ray equation of transport

Green's Theorem is developed for the spherically-symmetric steady-state cosmic-ray equation of transport in interplanetary space. By means of it the momentum distribution functionF _o(r,p), (r=heliocentric distance,p=momentum) can be determined in a regionr _arr_bwhen a source is specified throughout the region and the momentum spectrum is specified on the boundaries atr _a andr _b . Evaluation requires a knowledge of the Green's function which corresponds to the solution for monoenergetic particles released at heliocentric radiusr _o , Examples of Green's functions are given for the caser _a =0,r _b = and derived for the cases of finiter _a andr _b . The diffusion coefficient is assumed of the form = _o(p)r ^b . The treatment systematizes the development of all analytic solutions for steady-state solar and galactic cosmic-ray propagation and previous solutions form a subset of the present solutions.     

Genericity of the non-periodic solutions of the central force problem

We study the classical problem of two-dimensional motion of a particle in the field of a central force proportional to a real power of the distancer. for negative energy and (0, 2), each energy levelI _h is foliated by the invariant toriI _hc of constant angular momentumc and, by Liouville-Arnold's theorem, the flow on eachI _hc is conjugated to a linear flow of rotation number _h (c).A well-known result asserts that if we require _h (c) to be rational for every value ofh andc, the, must be equal to one (Kepler's problem). In this paper we prove that for almost every (0, 2) _h (c) is a non-constant continuous function ofc, for everyh<0. In particular, we deduce that motion under central potentials is generically non-periodic.Partially supported by CIRIT under grant No. EE88/2.