Injection of Oort Cloud comets: the fundamental role of stellar perturbations

Celestial Mechanics and Dynamical Astronomy, 102, pp. 111-132, http://dx.doi.org./10.1007/s10569-008-9140-yInternational audienc

The structure of invariant tori in a 3D galactic potential

We study in detail the structure of phase space in the neighborhood of stable periodic orbits in a rotating 3D potential of galactic type. We have used the color and rotation method to investigate the properties of the invariant tori in the 4D spaces of section. We compare our results with those of previous works and we describe the morphology of the rotational, as well as of the tube tori in the 4D space. We find sticky chaotic orbits in the immediate neighborhood of sets of invariant tori surrounding 3D stable periodic orbits. Particularly useful for galactic dynamics is the behavior of chaotic orbits trapped for long time between 4D invariant tori. We find that they support during this time the same structure as the quasi-periodic orbits around the stable periodic orbits, contributing however to a local increase of the dispersion of velocities. Finally we find that the tube tori do not appear in the 3D projections of the spaces of section in the axisymmetric Hamiltonian we examined.Comment: 26 pages, 34 figures, accepted for publication in the International Journal of Bifurcation and Chao

Application of new dynamical spectra of orbits in Hamiltonian systems

In the present article, we investigate the properties of motion in Hamiltonian systems of two and three degrees of freedom, using the distribution of the values of two new dynamical parameters. The distribution functions of the new parameters, define the S(g) and the S(w) dynamical spectra. The first spectrum definition, that is the S(g) spectrum, will be applied in a Hamiltonian system of two degrees of freedom (2D), while the S(w) dynamical spectrum will be deployed in a Hamiltonian system of three degrees of freedom (3D). Both Hamiltonian systems, describe a very interesting dynamical system which displays a large variety of resonant orbits, different chaotic components and also several sticky regions. We test and prove the efficiency and the reliability of these new dynamical spectra, in detecting tiny ordered domains embedded in the chaotic sea, corresponding to complicated resonant orbits of higher multiplicity. The results of our extensive numerical calculations, suggest that both dynamical spectra are fast and reliable discriminants between different types of orbits in Hamiltonian systems, while requiring very short computation time in order to provide solid and conclusive evidence regarding the nature of an orbit. Furthermore, we establish numerical criteria in order to quantify the results obtained from our new dynamical spectra. A comparison to other previously used dynamical indicators, reveals the leading role of the new spectra.Comment: Published in Nonlinear Dynamics (NODY) journal. arXiv admin note: text overlap with arXiv:1009.1993 by other author

On the connection between the Nekhoroshev theorem and Arnold Diffusion

The analytical techniques of the Nekhoroshev theorem are used to provide estimates on the coefficient of Arnold diffusion along a particular resonance in the Hamiltonian model of Froeschl\'{e} et al. (2000). A resonant normal form is constructed by a computer program and the size of its remainder $||R_{opt}||$ at the optimal order of normalization is calculated as a function of the small parameter $\epsilon$. We find that the diffusion coefficient scales as $D\propto||R_{opt}||^3$, while the size of the optimal remainder scales as $||R_{opt}|| \propto\exp(1/\epsilon^{0.21})$ in the range $10^{-4}\leq\epsilon \leq 10^{-2}$. A comparison is made with the numerical results of Lega et al. (2003) in the same model.Comment: Accepted in Celestial Mechanics and Dynamical Astronom

An Overview of the 13:8 Mean Motion Resonance between Venus and Earth

It is known since the seminal study of Laskar (1989) that the inner planetary system is chaotic with respect to its orbits and even escapes are not impossible, although in time scales of billions of years. The aim of this investigation is to locate the orbits of Venus and Earth in phase space, respectively to see how close their orbits are to chaotic motion which would lead to unstable orbits for the inner planets on much shorter time scales. Therefore we did numerical experiments in different dynamical models with different initial conditions -- on one hand the couple Venus-Earth was set close to different mean motion resonances (MMR), and on the other hand Venus' orbital eccentricity (or inclination) was set to values as large as e = 0.36 (i = 40deg). The couple Venus-Earth is almost exactly in the 13:8 mean motion resonance. The stronger acting 8:5 MMR inside, and the 5:3 MMR outside the 13:8 resonance are within a small shift in the Earth's semimajor axis (only 1.5 percent). Especially Mercury is strongly affected by relatively small changes in eccentricity and/or inclination of Venus in these resonances. Even escapes for the innermost planet are possible which may happen quite rapidly.Comment: 14 pages, 11 figures, submitted to CMD

Order statistics and heavy-tail distributions for planetary perturbations on Oort cloud comets

This paper tackles important aspects of comets dynamics from a statistical point of view. Existing methodology uses numerical integration for computing planetary perturbations for simulating such dynamics. This operation is highly computational. It is reasonable to wonder whenever statistical simulation of the perturbations can be much more easy to handle. The first step for answering such a question is to provide a statistical study of these perturbations in order to catch their main features. The statistical tools used are order statistics and heavy tail distributions. The study carried out indicated a general pattern exhibited by the perturbations around the orbits of the important planet. These characteristics were validated through statistical testing and a theoretical study based on Opik theory.Comment: 9 pages, 12 figures, submitted for publication in Astronomy and Astrophysic

Resonance Zones and Lobe Volumes for Volume-Preserving Maps

We study exact, volume-preserving diffeomorphisms that have heteroclinic connections between a pair of normally hyperbolic invariant manifolds. We develop a general theory of lobes, showing that the lobe volume is given by an integral of a generating form over the primary intersection, a subset of the heteroclinic orbits. Our definition reproduces the classical action formula in the planar, twist map case. For perturbations from a heteroclinic connection, the lobe volume is shown to reduce, to lowest order, to a suitable integral of a Melnikov function.Comment: ams laTeX, 8 figure

The Instability Transition for the Restricted 3-Body Problem. III. The Lyapunov Exponent Criterion

We establish a criterion for the stability of planetary orbits in stellar binary systems by using Lyapunov exponents and power spectra for the special case of the circular restricted 3-body problem (CR3BP). The centerpiece of our method is the concept of Lyapunov exponents, which are incorporated into the analysis of orbital stability by integrating the Jacobian of the CR3BP and orthogonalizing the tangent vectors via a well-established algorithm originally developed by Wolf et al. The criterion for orbital stability based on the Lyapunov exponents is independently verified by using power spectra. The obtained results are compared to results presented in the two previous papers of this series. It is shown that the maximum Lyapunov exponent can be used as an indicator for chaotic behaviour of planetary orbits, which is consistent with previous applications of this method, particularly studies for the Solar System. The chaotic behaviour corresponds to either orbital stability or instability, and it depends solely on the mass ratio of the binary components and the initial distance ratio of the planet relative to the stellar separation distance. Our theoretical results allow us to link the study of planetary orbital stability to chaos theory noting that there is a large array of literature on the properties and significance of Lyapunov exponents. Although our results are given for the special case of the CR3BP, we expect that it may be possible to augment the proposed Lyapunov exponent criterion to studies of planets in generalized stellar binary systems, which is strongly motivated by existing observational results as well as results expected from ongoing and future planet search missions.Comment: 10 pages, 8 figures, 3 tables; accepted by Astronomy and Astrophysic

Interplay Between Chaotic and Regular Motion in a Time-Dependent Barred Galaxy Model

We study the distinction and quantification of chaotic and regular motion in a time-dependent Hamiltonian barred galaxy model. Recently, a strong correlation was found between the strength of the bar and the presence of chaotic motion in this system, as models with relatively strong bars were shown to exhibit stronger chaotic behavior compared to those having a weaker bar component. Here, we attempt to further explore this connection by studying the interplay between chaotic and regular behavior of star orbits when the parameters of the model evolve in time. This happens for example when one introduces linear time dependence in the mass parameters of the model to mimic, in some general sense, the effect of self-consistent interactions of the actual N-body problem. We thus observe, in this simple time-dependent model also, that the increase of the bar's mass leads to an increase of the system's chaoticity. We propose a new way of using the Generalized Alignment Index (GALI) method as a reliable criterion to estimate the relative fraction of chaotic vs. regular orbits in such time-dependent potentials, which proves to be much more efficient than the computation of Lyapunov exponents. In particular, GALI is able to capture subtle changes in the nature of an orbit (or ensemble of orbits) even for relatively small time intervals, which makes it ideal for detecting dynamical transitions in time-dependent systems.Comment: 21 pages, 9 figures (minor typos fixed) to appear in J. Phys. A: Math. Theo

Entangling power of quantized chaotic systems

We study the quantum entanglement caused by unitary operators that have classical limits that can range from the near integrable to the completely chaotic. Entanglement in the eigenstates and time-evolving arbitrary states is studied through the von Neumann entropy of the reduced density matrices. We demonstrate that classical chaos can lead to substantially enhanced entanglement. Conversely, entanglement provides a novel and useful characterization of quantum states in higher dimensional chaotic or complex systems. Information about eigenfunction localization is stored in a graded manner in the Schmidt vectors, and the principal Schmidt vectors can be scarred by the projections of classical periodic orbits onto subspaces. The eigenvalues of the reduced density matrices are sensitive to the degree of wavefunction localization, and are roughly exponentially arranged. We also point out the analogy with decoherence, as reduced density matrices corresponding to subsystems of fully chaotic systems are diagonally dominant.Comment: 21 pages including 9 figs. (revtex