Exploring the 7:4 mean motion resonance—I: Dynamical evolution of classical transneptunian objects

In the transneptunian classical region (), an unexpected orbital excitation in eccentricity and inclination, dynamically distinct populations and the presence of chaotic regions are observed. For instance, the 7:4 mean motion resonance () appears to have been causing unique dynamical excitation according to observational evidences, namely, an apparent shallow gap in number density and anomalies in the colour distribution, both features enhanced near the 7:4 mean motion resonance location. In order to investigate the resonance dynamics, we present extensive computer simulation results totalizing almost 10,000 test particles under the effect of the four giant planets for the age of the solar system. A chaotic diffusion experiment was also performed to follow tracks in phase space over 4-5 Gyr. The 7:4 mean motion resonance is weakly chaotic causing irregular eccentricity and inclination evolution for billions of years. Most 7:4 resonant particles suffered significant eccentricities and/or inclinations excitation, an outcome shared even by those located in the vicinity of the resonance. Particles in stable resonance locking are rare and usually had 0.25<e<0.3. For other regions, 7:4 resonants had quite large mobility in phase space typically leaving the resonance (and being scattered) after reaching a critical e∼0.2. The escape happened in 10⁸-10⁹ yr time scales. Concerning the inclination dependence for 7:4 resonants, we found strong instability islands for approximately i>10°. Taking into account those particles still locked in the resonance at the end of the simulations, we determined a retainability of 12-15% for real 7:4 resonant transneptunian objects (TNOs). Lastly, our results demonstrate that classical TNOs associated with the 7:4 mean motion resonance have been evolving continuously until present with non-negligible mixing of populations.     

A perturbative treatment of motion near the 3/1 commensurability

A semianalytic perturbation theory for motion near the 3/1 commensurability in the planar elliptic restricted three-body problem is presented. The predictions of the theory are in good agreement with the features found on numerically generated surfaces of section; a global understanding of the phase space is achieved. The unusual features of the motion discovered by J. Wisdom (1982, Astron. J.87, 577–593; 1983a, Icarus56, 51–74) are explained. The principal cause of the large chaotic zone near the 3/1 commensurability is identified, and a new criterion for the existence of large-scale chaotic behavior is presented.     

Numerical simulation of the final stages of terrestrial planet formation

Three representative numerical simulations of the growth of the terrestrial planets by accretion of large protoplanets are presented. The mass and relative-velocity distributions of the bodies in these simulations are free to evolve simultaneously in response to close gravitational encounters and occasional collisions between bodies. The collisions between bodies, therefore, arise in a natural way and the assumption of expressions for the relative velocity distribution and the gravitational collision cross section is unnecessary. These simulations indicate that the growth of bodies with final masses approaching those of Venus and the Earth is possible, at least for the case of a two-dimensional system. Simulations assuming an initial uniform distribution of orbital eccentricities on the interval from 0 to e_max are found to produce final states containing too many bodies with masses which are too small when e_max < 0.10, while simulations with e_max > 0.20 result in too many catastrophic collisions between bodies thus preventing rapid accretion of planetary-size bodies. The e_max = 0.15 simulation ends with a state surprisingly similar to that of the present terrestrial planets and, therefore, provides a rough estimate of the range of radial sampling to be expected for the terrestrial planets.     

Modeling the 3-D secular planetary three-body problem: Discussion on the outer υ Andromedae planetary system

The three-dimensional secular behavior of a system composed of a central star and two massive planets is modeled semi-analytically in the frame of the general three-body problem. The main dynamical features of the system are presented in geometrical pictures allowing us to investigate a large domain of the phase space of this problem without time-expensive numerical integrations of the equations of motion and without any restriction on the magnitude of the planetary eccentricities, inclinations and mutual distance. Several regimes of motion of the system are observed. With respect to the secular angle Δ?, possible motions are circulations, oscillations (around 0° and 180°), and high-eccentricity/inclination librations in secular resonances. With respect to the arguments of pericenter, ω₁ and ω₂, possible motions are direct circulation and high-inclination libration around ±90° in the Lidov-Kozai resonance. The regions of transition between domains of different regimes of motion are characterized by chaotic behavior. We apply the analysis to the case of the two outer planets of the υ Andromedae system, observed edge-on. The topology of the 3-D phase space of this system is investigated in detail by means of surfaces of section, periodic orbits and dynamical spectra, mapping techniques and numerical simulations. We obtain the general structure of the phase space, and the boundaries of the spatial secular stability. We find that this system is secularly stable in a large domain of eccentricities and inclinations.     

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

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

A new dynamical model for the study of galactic structure

In the present article, we present a new gravitational galactic model, describing motion in elliptical as well as in disk galaxies, by suitably choosing the dynamical parameters. Moreover, a new dynamical parameter, the S(g) spectrum, is introduced and used, in order to detect islandic motion of resonant orbits and the evolution of the sticky regions. We investigate the regular or chaotic character of motion, with emphasis in the different dynamical models and make an extensive study of the sticky regions of the system. We use the classical method of the Poincaré r ? p_r phase plane and the new dynamical parameter of the S(g) spectrum. The L.C.E is used, in order to make an estimation of the degree of chaos in our galactic model. In both cases, the numerical calculations, suggest that our new model, displays a wide variety of families of regular orbits, compared to other galactic models. In addition to the regular motion, this new model displays also chaotic regions. Furthermore, the extent of the chaotic regions increases, as the value of the flatness parameter b of the model increases. Moreover, our simulations indicate, that the degree of chaos in elliptical galaxies, is much smaller than that in dense disk galaxies. In both cases numerical calculations show, that the degree of chaos increases linearly, as the flatness parameter b increases. In addition, a linear relationship between the critical value of angular momentum and the b parameter if found, in both cases (elliptical and disk galaxies). Some theoretical arguments to support the numerical outcomes are presented. Comparison with earlier work is also made.     

The range of validity of the two-body approximation in models of terrestrial planet accumulation: II. Gravitational cross sections and runaway accretion

The validity of the two-body approximation in calculating encounters between planetesimals has been evaluated as a function of the ratio of unperturbed planetesimal velocity (with respect to a circular orbit) to mutual escape velocity when their surfaces are in contact (V/V_e). Impact rates as a function of V/V_e are calculated to within ～20% by numerical integration of the equations of motion. It is found that when V/V_e > 0.4, the two-body approximation is a good one. At low velocities (V/V_e < 0.1) two-body “collision-course” trajectories fail to lead to impacts. On the other hand, at these low velocities many impacts result from encounter trajectories with unperturbed separation distances far beyond the two-body gravitational radius. These two effects tend to cancel, and the resulting impact rates remain within a factor of ～3 of the two-body value in spite of these major differences in the nature of the impact trajectories. Therefore, on the average, the two-body approximation is useful well below the value of V/V_e for which it fails to describe individual encounters, and the required corrections are not large. As a consequence of this “anomalous gravitational focusing” planetesimals will continue to interact even when their orbits are noncrossing. This reduces the difficulty with premature isolation of planetesimal embryos during accumulation. Quantitatively, when 0.06 ? V/V_e ? 0.2, the impact rate varies approximately with the fifth power of the radius of the larger body, and is about a factor of 3 above that predicted using the conventional two-body gravitational cross-section formula. At lower values of V/V_e , the impact rate increases less rapidly. Finally, at the lowest values of V/V_e (<.02), the impact rate increases only in proportion to the geometric cross section, as a consequence of the swarm being essentially two dimensional for large unperturbed encounter distances. The gravitational enhancement in effective cross section is thereby limited to a value of about 3000. This leads to an optimal size for growth of planetesimals from a swarm of given eccentricity, and places a limit on the extent of runaway accretion.     

Dynamical portrait of the Lixiaohua asteroid family

In this paper we present a comprehensive analysis of the dynamics in the region of the (3556) Lixiaohua asteroid family. The family lies in a particularly interesting region of the phase space, crossed by several two-body and three-body mean motion resonances. Also, members of this family can have close encounters with large asteroids, such as Ceres. We have identified the mean motion resonances which contribute to the long-term dynamical evolution of the family and our results confirm that the members of this family can be classified into a number of groups, exhibiting different dynamical behavior. We show for the first time that in the Lixiaohua region, apart from the chaotic diffusion in proper eccentricity and inclination (e _p and I _p ), there is at least one extended chaotic zone where several resonances overlap, thus giving rise to chaotic diffusion in proper semi-major axis (a _p ) as well. Using a code of Monte Carlo type, we simulate the evolution of the family, according to the model which combines the chaotic diffusion (in a _p , e _p and I _p ), Yarkovsky/YORP thermal effect and random walk in a _p due to the close encounters with massive asteroids. These simulations show that all these effects should be taken into account in order to accurately explain the observed distribution of family members in the space of proper elements, although a “minimal” model that accounts for chaotic diffusion in (e _p , I _p ), Yarkovsky-induced drift in a _p and random walk in a _p due to the close encounters with the most massive asteroids is enough to grossly characterize the shape of the family.     

Modulational instability of ion-acoustic waves in electron–positron–ion plasmas

The stability of modulation of ion-acoustic waves in a collisionless electron–positron–ion plasma with warm adiabatic ions is studied. Using the Krylov–Bogoliubov–Mitropolosky (KBM) perturbation technique a nonlinear Schrödinger equation governing the slow modulation of the wave amplitude is derived for the system. It is found that for given set of parameters having finite ion temperature ratio (T _i /T _e ) the waves are unstable for the values of k lying in the range k _min<k<k _max. On increasing the ion temperature ratio (T _i /T _e ), it is found that k _min and k _max, both decreases and product PQ increases. The range of unstable region shifts towards the small wave number k, as temperature ratio (T _i /T _e ) increases. The positron concentration and temperature ratio of positron to electron, change the unstable region slightly. As positron concentration increases both k _min and k _max for modulational instability increases and maximum value of the product PQ shifts towards the larger value of k.     

Some properties of a two-body system under the influence of the Galactic tidal field

In this paper the effect of the Galactic tidal field on a Sun–comet pair will be considered when the comet is situated in the Oort cloud and planetary perturbations can be neglected. First, two averaged models were created, one of which can be solved analytically in terms of Jacobi elliptic functions. In the latter system, switching between libration and circulation of the argument of perihelion is prohibited. The non-averaged equations of motion are integrated numerically in order to determine the regions of the ( e ,  i ) phase space in which chaotic orbits can be found, and an effort is made to explain why the chaotic orbits manifest in these regions only. It is evident that for moderate values of semimajor axis, a ∼50 000 au , chaotic orbits are found for ( e <0.15 , 40°≤ i ≤140°) as determined by integrating the evolution of the comet over a period of 10⁴ orbits. These regions of chaos increase in size with increasing semimajor axis. The typical e-folding times for these orbits range from around 600 Myr to 1 Gyr and thus are of little practical interest, as the time-scales for chaos arising from passing stars are much shorter. As a result of Galactic rotation, the chaotic regions in ( e ,  i ) phase space are not symmetric for prograde and retrograde orbits.     

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.     

A Map for a Group of Resonant Cases in a quartic Galactic Hamiltonian

We present a map for the study of resonant motion in a potential made up of two harmonic oscillators with quartic perturbing terms. This potential can be considered to describe motion in the central parts of non-rotating elliptical galaxies. The map is based on the averaged Hamiltonian. Adding on a semi-empirical basis suitable terms in the unperturbed averaged Hamiltonian, corresponding to the 1:1 resonant case, we are able to construct a map describing motion in several resonant cases. The map is used in order to find thex − p _x Poincare phase plane for each resonance. Comparing the results of the map, with those obtained by numerical integration of the equation of motion, we observe, that the map describes satisfactorily the broad features of orbits in all studied cases for regular motion. There are cases where the map describes satisfactorily the properties of the chaotic orbits as well.     

Stable Chaos in High-Order Jovian Resonances

In a previous publication (Tsiganis et al. 2000, Icarus146, 240-252), we argued that the occurrence of stable chaos in the 12/7 mean motion resonance with Jupiter is related to the fact that there do not exist families of periodic orbits in the planar elliptic restricted problem and in the 3-D circular problem corresponding to this resonance. In the present paper we show that nonexistence of resonant periodic orbits, both for the planar and for the 3-D problem, also occurs in other jovian resonances—namely the 11/4, 22/9, 13/6, and 18/7—where cases of real asteroids on stable-chaotic orbits have been identified. This property may provide a “protection mechanism”, leading to semiconfinement of chaotic orbits and extremely slow migration in the space of proper elements, so that diffusion is practically unrelated to the value of the Lyapunov time, T_L, of chaotic orbits. However, we show that, in more complicated dynamical models, the long-term evolution of chaotic orbits initiated in the vicinity of these resonances may also be governed by secular resonances. Finally, we find that stable-chaotic orbits have a characteristic spectrum of autocorrelation times: for the action conjugate to the critical argument the autocorrelation time is of the order of the Lyapunov time, while for the eccentricity- and inclination-related actions the autocorrelation time may be longer than 10³T_L. This behavior is consistent with the trajectory being sticky around a manifold of lower-than-full dimensionality in phase space (e.g., a 4-D submanifold of the 5-D energy manifold in a three-degrees-of-freedom autonomus Hamiltonian system) and reflects the inability of these “flawed” resonances to modify secular motion significantly, at least for times of the order of 200 Myr.     

Exploring the nature of orbits in a galactic model with a massive nucleus

In the present article, we use an axially symmetric galactic gravitational model with a disk–halo and a spherical nucleus, in order to investigate the transition from regular to chaotic motion for stars moving in the meridian (r,z) plane. We study in detail the transition from regular to chaotic motion, in two different cases: the time independent model and the time evolving model. In both cases, we explored all the available range regarding the values of the main involved parameters of the dynamical system. In the time dependent model, we follow the evolution of orbits as the galaxy develops a dense and massive nucleus in its core, as mass is transported exponentially from the disk to the galactic center. We apply the classical method of the Poincaré (r,p_r) phase plane, in order to distinguish between ordered and chaotic motion. The Lyapunov Characteristic Exponent is used, to make an estimation of the degree of chaos in our galactic model and also to help us to study the time dependent model. In addition, we construct some numerical diagrams in which we present the correlations between the main parameters of our galactic model. Our numerical calculations indicate, that stars with values of angular momentum L_z less than or equal to a critical value L_zc, moving near to the galactic plane, are scattered to the halo upon encountering the nuclear region and subsequently display chaotic motion. A linear relationship exists between the critical value of the angular momentum L_zc and the mass of the nucleus M_n. Furthermore, the extent of the chaotic region increases as the value of the mass of the nucleus increases. Moreover, our simulations indicate that the degree of chaos increases linearly, as the mass of the nucleus increases. A comparison is made between the critical value L_zc and the circular angular momentum L_z0 at different distances from the galactic center. In the time dependent model, there are orbits that change their orbital character from regular to chaotic and vise versa and also orbits that maintain their character during the galactic evolution. These results strongly indicate that the ordered or chaotic nature of orbits, depends on the presence of massive objects in the galactic cores of the galaxies. Our results suggest, that for disk galaxies with massive and prominent nuclei, the low angular momentum stars in the associated central regions of the galaxy, must be in predominantly chaotic orbits. Some theoretical arguments to support the numerically derived outcomes are presented. Comparison with similar previous works is also made.     

Orbital evolution of Uranus’s new outer satellites

Data on three recently discovered satellites of Uranus are used to determine basic evolutional parameters of their orbits: the extreme eccentricities and inclinations, as well as the circulation periods of the pericenter arguments and of the longitudes of the ascending nodes. The evolution is mainly investigated by analytically solving Hill’s double-averaged problem for the Uranus-Sun-satellite system, in which Uranus’s orbital eccentricity e _U and inclination i _U to the ecliptic are assumed to be zero. For the real model of Uranus’s evolving orbit with e _U≠0 and i _U≠0, we refine the evolutional parameters of the satellite orbits by numerically integrating the averaged system. Having analyzed the configuration and dynamics of the orbits of Uranus’s five outer satellites, we have revealed the possibility of their mutual crossings and obtained approximate temporal estimates.     

The rhomboidal symmetric four-body problem

We consider the planar symmetric four-body problem with two equal masses m ₁?=?m ₃?>?0 at positions (±x ₁(t),?0) and two equal masses m ₂?=?m ₄?>?0 at positions (0, ±x ₂(t)) at all times t, referred to as the rhomboidal symmetric four-body problem. Owing to the simplicity of the equations of motion this problem is well suited to study regularization of the binary collisions, periodic solutions, chaotic motion, as well as the four-body collision and escape manifolds. Furthermore, resonance phenomena between the two interacting rectilinear binaries play an important role.     

Pulsars and interstellar medium: Multiple regression analysis of related parameters

An empirical relation which relates the 408 MHz galactic continuum background temperature (408GCBT) to dispersion measures, position and radio-luminosity of 325 pulsars is obtained by means of multple stepwise regression analysis. This relation showns that pulsars may be considered as galactic probes for the distribution of 408GCBT and interstellar electron density (IED) in interstellar medium (ISM).Peculiar pulsars (O-C±2.5) point out galactic regions where the observedT ₄₀₈ are higher (or lower) andn _e lower (or higher) than the averaged ones.Normal pulsars (–2.5T ₄₀₈ andn _e are in agreement, on, the average.Standard pulsars (O-C±0.05) show galactic regions where observed and computedT ₄₀₈ andn _e are in good agreement. Recent models of pulsar disk systems, suggested by Michel and Dessler (1981) could justify the conclusions drawn for peculiar pulsars having O-C>2.5.