Integrable models of galactic discs with double nuclei

We introduce a new class of 2D mass models, whose potentials are of Stäckel form in elliptic coordinates. Our model galaxies have two separate strong cusps that form double nuclei. The potential and surface density distributions are locally axisymmetric near the nuclei and become highly non-axisymmetric outside the nucleus. The surface density diverges toward the cuspy nuclei with the law     Our model is sustained by four general types of regular orbits: butterfly , nucleophilic banana , horseshoe and aligned loop orbits. Horseshoes and nucleophilic bananas support the existence of cuspy regions. Butterflies and aligned loops control the non-axisymmetric shape of outer regions. Without any need for central black holes, our distributed mass models resemble the nuclei of M31 and NGC 4486B. It is also shown that the self-gravity of the stellar disc can prevent the double nucleus to collapse.     

Regular and chaotic dynamics in 3D reconnecting current sheets

We consider the possibility of particles being injected at the interior of a reconnecting current sheet (RCS), and study their orbits by dynamical systems methods. As an example we consider orbits in a 3D Harris type RCS. We find that, despite the presence of a strong electric field, a 'mirror' trapping effect persists, to a certain extent, for orbits with appropriate initial conditions within the sheet. The mirror effect is stronger for electrons than for protons. In summary, three types of orbits are distinguished: (i) chaotic orbits leading to escape by stochastic acceleration, (ii) regular orbits leading to escape along the field lines of the reconnecting magnetic component, and (iii) mirror-type regular orbits that are trapped in the sheet, making mirror oscillations. Dynamically, the latter orbits lie on a set of invariant KAM tori that occupy a considerable amount of the phase space of the motion of the particles. We also observe the phenomenon of 'stickiness', namely chaotic orbits that remain trapped in the sheet for a considerable time. A trapping domain, related to the boundary of mirror motions in velocity space, is calculated analytically. Analytical formulae are derived for the kinetic energy gain in regular or chaotic escaping orbits. The analytical results are compared with numerical simulations.     

Algorithmic regularization of the few-body problem

Using a modified leapfrog method as a basic mapping, we produce a new numerical integrator for the stellar dynamical few-body problem. We do not use coordinate transformation and the differential equations are not regularized, but the leapfrog algorithm gives regular results even for collision orbits. For this reason, application of extrapolation methods gives high precision. We compare the new integrator with several others and find it promising. Especially interesting is its efficiency for some potentials that differ from the Newtonian one at small distances.     

Orbits supporting bars within bars

High-resolution observations of the inner regions of barred disc galaxies have revealed many asymmetrical, small-scale central features, some of which are best described as secondary bars. Because orbital time-scales in the galaxy centre are short, secondary bars are likely to be dynamically decoupled from the main kiloparsec-scale bars. Here we show that regular orbits exist in such doubly barred potentials, and that they can support the bars in their motion. We find orbits in which particles remain on loops : closed curves which return to their original positions after two bars have come back to the same relative orientation. Stars trapped around stable loops could form the building blocks for a long-lived, doubly barred galaxy. Using the loop representation, we can find which orbits support the bars in their motion, and the constraints on the sizes and shapes of self-consistent double bars. In particular, it appears that a long-lived secondary bar may exist only when an inner Lindblad resonance is present in the primary bar, and that it would not extend beyond this resonance.     

Discrete Virial Theorem

We reexamine the classical virial theorem for bounded orbits of arbitrary autonomous Hamiltonian systems possessing both regular and chaotic orbits. New and useful forms of the virial theorem are obtained for natural Hamiltonian flows of arbitrary dimension. A discrete virial theorem is derived for invariant circles and periodic orbits of natural symplectic maps. A weak and a strong form of the virial theorem are proven for both flows and maps. While the Birkhoff Ergodic Theorem guarantees the existence of the relevant time averages for both regular and chaotic orbits, the convergence is very rapid for the former and extremely slow for the latter. This circumstance leads to a simple and efficient measure of chaoticity. The results are applied to several problems of current physical interest, including the Hénon–Heiles system, weak chaos in the standard map, and a 4D Froeschlé map.     

Stability of axial orbits in analytical galactic potentials

We study the stability of axial orbits in analytical galactic potentials as a function of the energy of the orbit and the ellipticity of the potential. The problem is solved by an analytical method, the validity of which is not limited to small amplitudes. The lines of neutral stability divide the parameter space in regions corresponding to different organizations of the main families of orbits in the symmetry planes.     

Chaos in cuspy triaxial galaxies with a supermassive black hole: a simple toy model

This paper summarises an investigation of chaos in a toy potential which mimics much of the behaviour observed for the more realistic triaxial generalisations of the Dehnen potentials, which have been used to model cuspy triaxial galaxies both with and without a supermassive black hole. The potential is the sum of an anisotropic harmonic oscillator potential, ${\text{V}}_{\text{0}} = \frac{1}{2}\left( {a^2 x^2 + b^2 y^2 + c^2 z^2 } \right)$ , and aspherical Plummer potential, ${\text{V}}_{\text{P}} = M_{BH} /\sqrt {r^2 + \varepsilon ^2 } $ , with $r^2 = x^2 + y^2 + z^2$ . Attention focuses on three issues related tothe properties of ensembles of chaotic orbits which impact on chaotic mixing and the possibility of constructing self-consistent equilibria:(1) What fraction of the orbits are chaotic? (2) How sensitive are the chaotic orbits, that is, how large are their largest (short time) Lyapunov exponents? (3) To what extent is the motion of chaotic orbits impeded by Arnold webs, that is, how 'sticky' are the chaotic orbits? These questions are explored as functions of the axis ratio a: b: c, black hole mass M _BH, softening length ε, and energy E with the aims of understanding how the manifestations of chaos depend onthe shape of the system and why the black hole generates chaos. The simplicity of the model makes it amenable to a perturbative analysis. That it mimics the behaviour of more complicated potentials suggests that much of this behaviour should be generic.     

Fast detection of chaotic behavior in galactic potentials

In the present paper we apply a new method of distinction between ordered and chaotic motion in galactic potentials. The method uses the Fourier Transform of a series of time intervals each one representing the time that elapsed between two successive points on the Poincaré surface of section. Examples of the methods ability to achieve an early and clear detection of an orbit's behavior are provided using two galactic potentials. The new method can also be applied in order to have an early distinction between ordered and sticky orbits. The method is generalized in order to be used in models with more than two dimensions. Finally we have tried to find an one‐number index to give us the nature of the orbit instead of checking by eye the whole spectrum. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Periodic Orbits in Analytical Planar Galactic Potentials

We study the regular families of periodic orbits in an analytical planar galactic potential, using the method of Lindstedt. We obtain analytical expressions describing these orbits, validity of which is not limited to small amplitudes. We can delimit, in the space of the parameters, the domain of existence of each family of orbits. This revised version was published online in July 2006 with corrections to the Cover Date.     

Stable Chaos versus Kirkwood Gaps in the Asteroid Belt: A Comparative Study of Mean Motion Resonances

We have shown, in previous publications, that stable chaos is associated with medium/high-order mean motion resonances with Jupiter, for which there exist no resonant periodic orbits in the framework of the elliptic restricted three-body problem. This topological “defect” results in the absence of the most efficient mechanism of eccentricity transport (i.e., large-amplitude modulation on a short time scale) in three-body models. Thus, chaotic diffusion of the orbital elements can be quite slow, while there can also exist a nonnegligible set of chaotic orbits which are semiconfined (stable chaos) by “quasi-barriers” in the phase space. In the present paper we extend our study to all mean motion resonances of order q≤9 in the inner main belt (1.9-3.3 AU) and q≤7 in the outer belt (3.3-3.9 AU). We find that, out of the 34 resonances studied, only 8 possess resonant periodic orbits that are continued from the circular to the elliptic three-body problem (regular families), namely, the 2/1, 3/1, 4/1, and 5/2 in the inner belt and the 7/4, 5/3, 11/7, and 3/2 in the outer belt. Numerical results indicate that the 7/3 resonance also carries periodic orbits but, unlike the aforementioned resonances, 7/3-periodic orbits belong to an irregular family. Note that the five inner-belt resonances that carry periodic orbits correspond to the location of the main Kirkwood gaps, while the three outer-belt resonances correspond to gaps in the distribution of outer-belt asteroids noted by Holman and Murray (1996, Astron. J.112, 1278-1293), except for the 3/2 case where the Hildas reside. Fast, intermittent eccentricity increase is found in resonances possessing periodic orbits. In the remaining resonances the time-averaged elements of chaotic orbits are, in general, quite stable, at least for times t∼250 Myr. This slow diffusion picture does not change qualitatively, even if more perturbing planets are included in the model.     

A simple algorithm for orbit classification

We describe a simple algorithm for classifying orbits into orbit families. This algorithm works by finding patterns in the sign changes of the principal coordinates. Orbits in the logarithmic potential are studied as an application; we classify orbits into boxlet families and examine the influence of the core radius on the set of stable orbit families.     

The tidal stripping of satellites

We present an improved analytic calculation for the tidal radius of satellites and test our results against N -body simulations.
The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star within the satellite . We demonstrate that this last point is critical and suggest using three tidal radii to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power-law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically.
Over short times (≲1–2 Gyr ∼1 satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.     

On the 3D dynamics and morphology of inner rings

We argue that inner rings in barred spiral galaxies are associated with specific 2D and 3D families of periodic orbits located just beyond the end of the bar. These are families located between the inner radial ultraharmonic 4 : 1 resonance and corotation. They are found in the upper part of a type-2 gap of the x1 characteristic, and can account for the observed ring morphologies without any help from families of the x1-tree. Due to the evolution of the stability of all these families, the ring shapes that are favoured are mainly ovals, as well as polygons with 'corners' on the minor axis, on the sides of the bar. On the other hand, pentagonal rings, or rings of the NGC 7020-type hexagon, should be less probable. The orbits that make the rings belong in their vast majority to 3D families of periodic orbits and orbits trapped around them.     

Finding how many isolating integrals of motion an orbit obeys

The correlation dimension, that is the dimension obtained by computing the correlation function of pairs of points of a trajectory in phase space, is a numerical technique introduced in the field of non-linear dynamics in order to compute the dimension of the manifold in which an orbit moves, without the need of knowing the actual equations of motion that give rise to the trajectory. This technique has been proposed in the past as a method to measure the dimension of stellar orbits in astronomical potentials, that is the number of isolating integrals of motion the orbits obey. Although the algorithm can in principle yield that number, some care has to be taken in order to obtain good results. We studied the relevant parameters of the technique, found their optimal values, and tested the validity of the method on a number of potentials previously studied in the literature, using the Smaller Alignment Index (SALI), Lyapunov exponents and spectral dynamics as gauges.     

Orbital dynamics of three-dimensional bars – I. The backbone of three-dimensional bars. A fiducial case

In this series of papers we investigate the orbital structure of three-dimensional (3D) models representing barred galaxies. In the present introductory paper we use a fiducial case to describe all families of periodic orbits that may play a role in the morphology of three-dimensional bars. We show that, in a 3D bar, the backbone of the orbital structure is not just the x1 family, as in two-dimensional (2D) models, but a tree of 2D and 3D families bifurcating from x1. Besides the main tree we have also found another group of families of lesser importance around the radial 3:1 resonance. The families of this group bifurcate from x1 and influence the dynamics of the system only locally. We also find that 3D orbits elongated along the bar minor axis can be formed by bifurcations of the planar x2 family. They can support 3D bar-like structures along the minor axis of the main bar. Banana-like orbits around the stable Lagrangian points build a forest of 2D and 3D families as well. The importance of the 3D x1-tree families at the outer parts of the bar depends critically on whether they are introduced in the system as bifurcations in z or in   z˙   .     

Dark halos acting as chaos controllers in asymmetric triaxial galaxy models

We study the regular or chaotic character of orbits in a 3D dynamical model,describing a triaxial galaxy surrounded by a spherical dark halo component.Our numerical experiments suggest that the percentage of chaotic orbits decreases exponentially as the mass of the dark halo increases.A linear increase of the percentage of the chaotic orbits was observed as the scale length of the halo component increases. In order to distinguish between regular and chaotic motion,we chose to use the total angular momentum ...