The resonant structure of Jupiter's Trojan asteroids – II. What happens for different configurations of the planetary system

In a previous paper, we have found that the resonance structure of the present Jupiter Trojan swarms could be split up into four different families of resonances. Here, in a first step, we generalize these families in order to describe the resonances occurring in Trojan swarms embedded in a generic planetary system. The location of these families changes under a modification of the fundamental frequencies of the planets and we show how the resonant structure would evolve during a planetary migration. We present a general method, based on the knowledge of the fundamental frequencies of the planets and on those that can be reached by the Trojans, which makes it possible to predict and localize the main events arising in the swarms during migration. In particular, we show how the size and stability of the Trojan swarms are affected by the modification of the frequencies of the planets. Finally, we use this method to study the global dynamics of the Jovian Trojan swarms when Saturn migrates outwards. Besides the two resonances found by Morbidelli et al. which could have led to the capture of the current population just after the crossing of the 2:1 orbital resonance, we also point out several sequences of chaotic events that can influence the Trojan population.     

The autoparametric 3:2 resonance in conservative systems

In the resonance model, high‐frequency quasi‐periodic oscillations (QPOs) are supposed to be a consequence of nonlinear resonance between modes of oscillations occurring within the innermost parts of an accretion disk. Several models with a prescribed mode–mode interaction were proposed in order to explain the characteristic properties of the resonance in QPO sources. In this paper, we examine nonlinear oscillations of a system having a quadratic nonlinearity and we show that this case is particularly relevant for QPOs. We present a very convenient way how to study autoparametric resonances of a fully general system using the method of multiple scales. We concentrate to conservative systems and discuss their behavior near the 3:2 parametric resonance. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Study of resonances for the restricted 3‐body problem

Our aim is to identify and classify mean‐motion resonances (MMRs) for the coplanar circular restricted three‐body problem (CR3BP) for mass ratios between 0.10 and 0.50. Our methods include the maximum Lyapunov exponent, which is used as an indicator for the location of the resonances, the Fast Fourier Transform (FFT) used for determining what kind of resonances are present, and the inspection of the orbital elements to classify the periodicity. We show that the 2:1 resonance occurs the most frequently. Among other resonances, the 3:1 resonance is the second most common, and furthermore both 3:2 and 5:3 resonances occur more often than the 4:1 resonance. Moreover, the resonances in the coplanar CR3BP are classified based on the behaviour of the orbits. We show that orbital stability is ensured for high values of resonance (i.e., high ratios) where only a single resonance is present. The resonances attained are consistent with the previously established resonances for the solar system, i.e., specifically, in regards to the asteroid belt. Previous work employed digital filtering and Lyapunov characteristic exponents to determine stochasticity of the eccentricity, which is found to be consistent with our usage of Lyapunov exponents as an alternate approach based on varying the mass ratio instead of the eccentricity. Our results are expected to be of principal interest to future studies, including augmentations to observed or proposed resonances, of extra‐solar planets in binary stellar systems (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analyticity, scaling and renormalization for some complex analytic dynamical systems

We review some results about the analytic structure of Lindstedt series for some complex analytic dynamical systems: in particular, we consider Hamiltonian maps (like the standard map and its generalizations), the semi-standard map and Siegel's problem of the linearization of germs of holomorphic diffeomorphisms of ( ,0). The analytic structure of those series can be studied numerically using Padé approximants, and one can show the existence of natural boundaries for real, diophantine values of the rotation number; by complexifying the rotation number, we show how these natural boundaries arise from the accumulation of singularities due to resonances, providing a new intuitive insight into the mechanism of the break-down of invariant KAM curves. Moreover, we study the Lindstedt series at resonances, i.e. for rational values of the rotation number, by suitably rescaling to 0 the value of the perturbative parameter, and a simple analytic structure emerges. Finally, we present some proofs for the simplest models and relate these results to renormalization ideas.     

Alfvénic resonant cavities in the solar atmosphere: Simple aspects

We investigate the propagation of Alfvén waves in a simple medium consisting of three uniform layers; each layer is characterized by a different value for the Alfvén speed, _A. We show how the central layer can act as a resonant cavity under quite general conditions. If the cavity is driven externally, by an incident wave in one of the outer layers, there result resonant transmission peaks, which allow large energy fluxes to enter the cavity from outside. The transmission peaks result from the destructive interference between a wave which leaks out of the cavity, and a directly reflected wave. We show that there are two types of resonances. The first type occurs when the cavity has the largest (or smallest) of the three Alfvén speeds; this situation occurs on coronal loops. The second type occurs when the cavity Alfvén speed is intermediate between the other two values of _A; this situation may occur on solar spicules. Significant heating of the cavity can occur if the waves are damped. We show that if the energy lost to heat greatly exceeds the energy lost by leakage out of the cavity, then the cavity heating can be independent of the damping rate. This conclusion is shown to apply to coronal resonances and to the spicule resonances. This conclusion agrees with a point made by Ionson (1982) in connection with the coronal resonances. Except for a numerical factor of order unity, we recover Ionson's expression for the coronal heating rate. However, Ionson's qualities are much too large. For solar parameters, the maximum quality is of the order of 100, but the heating is independent of the damping rate only when dissipation reduces the quality to less than about 10.     

Resonant behavior of comet Halley and the Orionid stream

Comet 1P/Halley has the unique distinction of having a very comprehensive set of observational records for almost every perihelion passage from 240 B.C. This has helped to constrain theoretical models pertaining to its orbital evolution. Many previous works have shown the active role of mean motion resonances (MMR) in the evolution of various meteoroid streams. Here, we look at how various resonances, especially the 1:6 and 2:13 MMR with Jupiter, affect comet 1P/Halley and thereby enhance the chances of meteoroid particles getting trapped in resonance, leading to meteor outbursts in some particular years. Comet Halley itself librated in the 2:13 resonance from 240 B.C. to 1700 A.D. and in the 1:6 resonance from 1404 B.C. to 690 B.C., while stream particles can survive for time scales of the order of 10,000 yr and 1,000 yr in the 1:6 and 2:13 resonances, respectively. This determines the long‐term dynamical evolution and stream structure, influencing the occurrence of Orionid outbursts. Specifically, we are able to correlate the occurrence of enhanced meteor phenomena seen between 1436–1440, 1933–1938, and 2006–2010 with the 1:6 resonance and meteor outbursts in 1916 and 1993 with the 2:13 resonance. Ancient as well as modern observational records agree with these theoretical simulations to a very good degree.     

Dynamics of the TrES‐2 system

The TrES‐2 system harbors one planet which was discovered with the transit technique. In this work we investigate the dynamical behavior of possible additional, lower‐mass planets. We identify the regions where such planets can move on stable orbits and show how they depend on the initial eccentricity and inclination. We find, that there are stable regions inside and outside the orbit of TrES‐2b where additional, smaller planets can move. We also show that those planets can have a large orbital inclination which makes a detection with the transit technique very difficult (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hydromagnetic wave coupling in the magnetosphere—Plasmapause effects on impulse-excited resonances

Recent analytical and numerical modelling has demonstrated the possibility that impulsively-stimulated compressional hydromagnetic cavity resonances can drive local field-line resonances in the magnetosphere. This paper extends the modelling to include axisymmetric plasmapause structures with realistic radial variation in the magnetospheric cavity. The results show that: (a) the plasmapause plays an important rôle in determining which cavity resonances are dominant; (b) when the wave fields are significantly non-axisymmetric, additional cavity resonances are evident which are at least partly trapped within the plasmasphere; (c) the position of the plasmapause determines where (and whether) cavity resonances couple significantly to field-line resonances; (d) for the small “azimuthal” wavenumber chosen, there is no evidence of a compressional surface wave on the plasmapause.     

Confined chaotic motion in three-body resonances: trapping of trans-Neptunian material induced by gas-drag

Jiang & Yeh proposed gas-drag-induced resonant capture as a mechanism able to explain the dominant 3:2 resonance observed in the trans-Neptunian belt. Using a model of a disc–star–planet system they concluded that gaseous drag in a protoplanetary disc can trap trans-Neptunian object (TNO) embryos into the 3:2 resonance rather easily although it could not trap objects into the 2:1 resonance. Here we further investigate this scenario using numerical simulations within the context of the planar restricted four-body problem by including both present-day Uranus and Neptune. Our results show that mean motion and corotation resonances are possible and trapping into both the 3:2 and 2:1 resonances as well as other resonances is observed. The associated corotation centres may easily form larger planetesimals from smaller ones. Corotation resonances evolve into pure Lindblad resonances in a time-scale of 0.5 Myr. The non-linear corotation and mean motion resonances produced are very size selective. The 3:2 resonance is dominant for submetric particles but for larger particles the 2:1 resonance is stronger. In summary, our calculations show that confined chaotic motion around the resonances not only increases trapping efficiency but also the orbital eccentricities of the trapped material, modifying the relative abundance of trapped particles in different resonances. If we assume a more compact planetary system, instead of using the present-day values of the orbital elements of Uranus and Neptune, our results remain largely unchanged.     

Capture into first-order resonances and long-term stability of pairs of equal-mass planets

Massive planets form within the lifetime of protoplanetary disks, and therefore, they are subject to orbital migration due to planet–disk interactions. When the first planet reaches the inner edge of the disk, its migration stops and consequently the second planet ends up locked in resonance with the first one. We detail how the resonant trapping works comparing semi-analytical formulae and numerical simulations. We restrict to the case of two equal-mass coplanar planets trapped in first-order resonances, but the method can be easily generalized. We first describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using series expansions up to different orders in eccentricity as well as a non-expanded Hamiltonian. Then we show that during convergent migration the planets evolve along these families of equilibrium points. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict precisely analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that in most cases the resonant configurations achieved by migration become unstable after the removal of the protoplanetary disk. Here we probe the stability of the resonances as a function of planetary mass. For this purpose, we fictitiously increase the masses of resonant planets, adiabatically maintaining the low-amplitude libration regime until instability occurs. We discuss two hypotheses for the instability, that of a low-order secondary resonance of the libration frequency with a fast synodic frequency of the system, and that of minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter nevertheless play the decisive role in the destabilization of resonant pairs. We show evidence that as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.     

Secular spin dynamics of inner main-belt asteroids

Understanding the evolution of asteroid spin states is challenging work, in part because asteroids have a variety of orbits, shapes, spin states, and collisional histories but also because they are strongly influenced by gravitational and non-gravitational (YORP) torques. Using efficient numerical models designed to investigate asteroid orbit and spin dynamics, we study here how several individual asteroids have had their spin states modified over time in response to these torques (i.e., 951 Gaspra, 60 Echo, 32 Pomona, 230 Athamantis, 105 Artemis). These test cases which sample semimajor axis and inclination space in the inner main belt, were chosen as probes into the large parameter space described above. The ultimate goal is to use these data to statistically characterize how all asteroids in the main belt population have reached their present-day spin states. We found that the spin dynamics of prograde-rotating asteroids in the inner main belt is generally less regular than that of the retrograde-rotating ones because of numerous overlapping secular spin-orbit resonances. These resonances strongly affect the spin histories of all bodies, while those of small asteroids (?40 km) are additionally influenced by YORP torques. In most cases, gravitational and non-gravitational torques cause asteroid spin axis orientations to vary widely over short (?1 My) timescales. Our results show that (951) Gaspra has a highly chaotic rotation state induced by an overlap of the s and s₆ spin-orbit resonances. This hinders our ability to investigate its past evolution and infer whether thermal torques have acted on Gaspra's spin axis since its origin.     

Reducing the probability of capture into resonance

A migrating planet can capture planetesimals into mean motion resonances. However, resonant trapping can be prevented when the drift or migration rate is sufficiently high. Using a simple Hamiltonian system for first- and second-order resonances, we explore how the capture probability depends on the order of the resonance, drift rate and initial particle eccentricity. We present scaling factors as a function of the planet mass and resonance strength to estimate the planetary migration rate above which the capture probability drops to less than half. Applying our framework to multiple extrasolar planetary systems that have two planets locked in resonance, we estimate lower limits for the outer planet's migration rate, allowing resonance capture of the inner planet.
Mean motion resonances are comprised of multiple resonant subterms. We find that the corotation subterm can reduce the probability of capture when the planet eccentricity is above a critical value. We present factors that can be used to estimate this critical planet eccentricity. Applying our framework to the migration of Neptune, we find that Neptune's eccentricity is near the critical value that would make its 2 : 1 resonance fail to capture twotinos. The capture probability is affected by the separation between resonant subterms and so is also a function of the precession rates of the longitudes of periapse of both planet and particle near resonance.     

Long-term evolution of the spin of Mercury: I. Effect of the obliquity and core-mantle friction

The present obliquity of Mercury is very low (less than 0.1°), which led previous studies to always adopt a nearly zero obliquity during the planet’s past evolution. However, the initial orientation of Mercury’s rotation axis is unknown and probably much different than today. As a consequence, we believe that the obliquity could have been significant when the rotation rate of the planet first encountered spin-orbit resonances. In order to compute the capture probabilities in resonance for any evolutionary scenario, we present in full detail the dynamical equations governing the long-term evolution of the spin, including the obliquity contribution.The secular spin evolution of Mercury results from tidal interactions with the Sun, but also from viscous friction at the core-mantle boundary. Here, this effect is also regarded with particular attention. Previous studies show that a liquid core enhances drastically the chances of capture in spin-orbit resonances. We confirm these results for null obliquity, but we find that the capture probability generally decreases as the obliquity increases. We finally show that, when core-mantle friction is combined with obliquity evolution, the spin can evolve into some unexpected configurations as the synchronous or the 1/2 spin-orbit resonance.     

A modified peak‐bagging technique for fitting low‐𝓁 solar p‐modes

We introduce a modified version of a standard power spectrum ‘peak‐bagging’ technique which is designed to gain some of the advantages that fitting the entire low‐degree p‐mode power spectrum simultaneously would bring, but without the problems involved in fitting a model incorporating many hundreds of parameters. Employing Monte‐Carlo simulations we show that by using this modified fitting code it is possible to determine the true background level in the vicinity of the p‐mode peaks. In addition to this we show how small biases in other mode parameters, which are related to inaccurate estimates of the true background, are also consequently removed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A planetary system with an escaping Mars

The chaotic behaviour of the motion of the planets in our Solar System is well established. In this work to model a hypothetical extrasolar planetary system our Solar System was modified in such a way that we replaced the Earth by a more massive planet and let the other planets and all the orbital elements unchanged. The major result of former numerical experiments with a modified Solar System was the appearance of a chaotic window at κ _E ∈ (4, 6), where the dynamical state of the system was highly chaotic and even the body with the smallest mass escaped in some cases. On the contrary for very large values of the mass of the Earth, even greater than that of Jupiter regular dynamical behaviour was observed. In this paper the investigations are extended to the complete Solar System and showed, that this chaotic window does still exist. Tests in different ‘Solar Systems’ clarified that including only Jupiter and Saturn with their actual masses together with a more ‘massive’ Earth (4 < κ _E < 6) perturbs the orbit of Mars so that it can even be ejected from the system. Using the results of the Laplace‐Lagrange secular theory we found secular resonances acting between the motions of the nodes of Mars, Jupiter and Saturn. These secular resonances give rise to strong chaos, which is the cause of the appearance of the instability window. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Impulse-excited hydromagnetic cavity and field-line resonances in the magnetosphere

We present detailed numerical results from a model which determines the time development of hydromagnetic waves within a hemi-cylindrical magnetospheric cavity subject to a short-duration compressional stimulus at the magnetopause. The model allows a realistic radial variation of Alfvén speed, arbitrary axial asymmetry, and the inclusion of ionospheric Joule dissipation.The results show the development of a set of compressional cavity resonances, and their coupling to a corresponding set of field-line resonances at positions where the cavity eigenperiods match the uncoupled field-line eigenperiods. Transient solutions having variable period with radial distance are important in establishing the field-line resonances, and can be locally dominant in certain regions of the outer magnetosphere.The variation with axial wavenumber m of the cavity eigenperiods is obtained, and it is shown that the magnetopause impulse drives the field-line resonances most efficiently for m ∼ 3. This suggests that small azimuthal wavenumbers are important in the terrestrial magnetosphere under analogous conditions. A variety of waveform and polarisation diagrams is given to aid in observational identification of such coupled resonances.     

Connection between spectral complexity and X‐ray weakness: testing the reflection interpretation

In recent work it was demonstrated that narrow‐line Seyfert 1 galaxies, which possessed spectral complexity in the 2–10 keV band were at the same time X‐ray weak. In this contribution I show how X‐ray weakness can be understood in the context of reflection and light bending picture. In fact, X‐ray weakness should be expected from objects that are in a reflection dominated state. With simultaneous UV and X‐ray data available with most XMM‐Newton observations, an estimate of the X‐ray weakness is relatively straightforward. As such, it is an easy way to substantiate conclusions of reflection dominated spectra, and we use this method to examine recent claims. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase-space volume of regions of trapped motion: multiple ring components and arcs

The phase-space volume of regions of regular or trapped motion, for bounded or scattering systems with two degrees of freedom respectively, displays universal properties. In particular, sudden reductions in the phase-space volume or gaps are observed at specific values of the parameter which tunes the dynamics; these locations are approximated by the stability resonances. The latter are defined by a resonant condition on the stability exponents of a central linearly stable periodic orbit. We show that, for more than two degrees of freedom, these resonances can be excited opening up gaps, which effectively separate and reduce the regions of trapped motion in phase space. Using the scattering approach to narrow rings and a billiard system as example, we demonstrate that this mechanism yields rings with two or more components. Arcs are also obtained, specifically when an additional (mean-motion) resonance condition is met. We obtain a complete representation of the phase-space volume occupied by the regions of trapped motion.     

Planets in the asteroid belt

Abstract— The main asteroid belt has lost >99.9% of its solid mass since the time at which the planets were forming, according to models for the protoplanetary nebula. Here we show that the primordial asteroid belt could have been cleared efficiently if much of the original mass accreted to form planetsized bodies, which were capable of perturbing one another into unstable orbits. We provide results from 25 N‐body integrations of up to 200 planets in the asteroid belt, with individual masses in the range 0.017–0.33 Earth masses. In the simulations, these bodies undergo repeated close encounters which scatter one another into unstable resonances with the giant planets, leading to collision with the Sun or ejection from the solar system. In response, the giant planets' orbits migrate radially and become more circular. This reduces the size of the main‐belt resonances and the clearing rate, although clearing continues. If ～3 Earth masses of material was removed from the belt this way, Jupiter and Saturn would initially have had orbital eccentricities almost twice their current values. Such orbits would have made Jupiter and Saturn 10–100x more effective at clearing material from the belt than they are on their current orbits. The time required to remove 90% of the initial mass from the belt depends sensitively on the giant planets' orbits, and weakly on the masses of the asteroidal planets. 18 of the 25 simulations end with no planets left in the belt, and the clearing takes up to several hundred million years. Typically, the last one or two asteroidal planets are removed by interactions with planets in the terrestrial region     

Stability of the Triangular Libration Points in the Unrestricted Planar Problem of a Symmetric Rigid Body and a Point Mass

In papers (Godziewski and Maciejewski, 1998a, b, 1999), we investigate unrestricted, planar problem of a dynamically symmetric rigid body and a sphere. Following the original statement of the problem by Kokoriev and Kirpichnikov (1988), we assume that the potential of the rigid body is approximated by the gravitational field of a dumb-bell. The model is described in terms of a 2D Hamiltonian depending on three parameters.In this paper, we investigate the stability of triangular equilibria permissible by the dynamics of the model, under the assumption of low-order resonances. We analyze all resonances of order smaller than four, and we examine the stability with application of theorems by Markeev and Sokolsky. These are the possible following cases: the non-diagonal resonance of the first order with two null characteristic frequencies (unstable); resonances of the first order with one nonzero frequency (diagonal and non-diagonal, stable and unstable); the second-order resonance, which is non-diagonal and stable, and the third-order resonance which is generically unstable, except for three points in the parameters' space, corresponding to stable equilibria.We discuss a perturbed version of Kokoriev and Kirpichnikov model, and we find that if the perturbation is small and depends on the coordinates only, the triangular equilibria persist, except if for the unperturbed equilibria the first-order resonance occurs. We show that the resonances of the order higher than two are also preserved if the perturbation acts.