Showing 1–31 of 31 results for author: Tsiganis, K

SELENA: Semi-analytical Integrator for Lunar Artificial Satellites

Authors: Christos Efthymiopoulos, Kleomenis Tsiganis, Ioannis Gkolias, Michalis Gaitanas, Carlos Yanez

Abstract: The present report summarizes the main theory and implementation steps associated with SELENA (SEmi-anaLytical intEgrator for a luNar Artificial satellite), i.e. the semi-analytical propagator for lunar satellite orbits developed in the framework of the the R&T R-S20/BS-0005-062 CNES research activity in collaboration between the University of Padova (UniPd), and the Aristotle University of Thessa… ▽ More The present report summarizes the main theory and implementation steps associated with SELENA (SEmi-anaLytical intEgrator for a luNar Artificial satellite), i.e. the semi-analytical propagator for lunar satellite orbits developed in the framework of the the R&T R-S20/BS-0005-062 CNES research activity in collaboration between the University of Padova (UniPd), and the Aristotle University of Thessaloniki (AUTH), both acting as contractors with CNES. A detailed account of the method, algorithms and symbolic manipulations employed in the derivation of the final theory are described in detail in this report: they invoke the use of canonical perturbation theory in the form of Lie series computed in `closed form', i.e., without expansions in the satellite's orbital eccentricity. These algorithms are provided in the form of a symbolic package accompanying the present report. The package contains symbolic algebra programs, as well as explicit data files containing the final Hamiltonian, equations of motion and transformations (i.e. the coefficients and exponents of each variable in each term) leading to the averaging of the short-periodic terms in the satellite's equations of motion. △ Less

Submitted 21 September, 2023; originally announced September 2023.

Comments: Final Report for the CNES R&T R-S20/BS-0005-062 Research Activity "Semi-analytical theory for the motion of lunar artificial satellites"

The Perturbed Full Two-Body Problem: Application to Post-DART Didymos

Authors: Alex J. Meyer, Harrison F. Agrusa, Derek C. Richardson, R. Terik Daly, Oscar Fuentes-Muñoz, Masatoshi Hirabayashi, Patrick Michel, Colby C. Merrill, Ryota Nakano, Andrew F. Cheng, Brent Barbee, Olivier S. Barnouin, Steven R. Chesley, Carolyn M. Ernst, Ioannis Gkolias, Nicholas A. Moskovitz, Shantanu P. Naidu, Petr Pravec, Petr Scheirich, Cristina A. Thomas, Kleomenis Tsiganis, Daniel J. Scheeres

Abstract: With the successful impact of the NASA DART spacecraft in the Didymos-Dimorphos binary asteroid system, we provide an initial analysis of the post-impact perturbed binary asteroid dynamics. To compare our simulation results with observations, we introduce a set of "observable elements" calculated using only the physical separation of the binary asteroid, rather than traditional Keplerian elements.… ▽ More With the successful impact of the NASA DART spacecraft in the Didymos-Dimorphos binary asteroid system, we provide an initial analysis of the post-impact perturbed binary asteroid dynamics. To compare our simulation results with observations, we introduce a set of "observable elements" calculated using only the physical separation of the binary asteroid, rather than traditional Keplerian elements. Using numerical methods that treat the fully spin-orbit-coupled dynamics, we estimate the system's mass and the impact-induced changes in orbital velocity, semimajor axis, and eccentricity. We find that the changes to the mutual orbit depend strongly on the separation distance between Didymos and Dimorphos at the time of impact. If Dimorphos enters a tumbling state after the impact, this may be observable through changes in the system's eccentricity and orbit period. We also find that any DART-induced reshaping of Dimorphos would generally reduce the required change in orbital velocity to achieve the measured post-impact orbit period and will be assessed by the ESA Hera mission in 2027. △ Less

Submitted 31 July, 2023; originally announced July 2023.

Comments: Accepted for publication in PSJ

Predictions for the Dynamical States of the Didymos System before and after the Planned DART Impact

Authors: Derek C. Richardson, Harrison F. Agrusa, Brent Barbee, William F. Bottke, Andrew F. Cheng, Siegfried Eggl, Fabio Ferrari, Masatoshi Hirabayashi, Özgür Karatekin, Jay McMahon, Stephen R. Schwartz, Ronald-Louis Ballouz, Adriano Campo Bagatin, Elisabetta Dotto, Eugene G. Fahnestock, Oscar Fuentes-Muñoz, Ioannis Gkolias, Douglas P. Hamilton, Seth A. Jacobson, Martin Jutzi, Josh Lyzhoft, Rahil Makadia, Alex J. Meyer, Patrick Michel, Ryota Nakano , et al. (11 additional authors not shown)

Abstract: NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $β$ at a reali… ▽ More NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $β$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $β$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 23 pages, 13 figures, published in PSJ

Journal ref: Planet. Sci. J. 3 157 (2022)

Libration-induced Orbit Period Variations Following the DART Impact

Authors: Alex J. Meyer, Ioannis Gkolias, Michalis Gaitanas, Harrison F. Agrusa, Daniel J. Scheeres, Kleomenis Tsiganis, Petr Pravec, Lance A. M. Benner, Fabio Ferrari, Patrick Michel

Abstract: The Double Asteroid Redirection Test (DART) mission will be the first test of a kinetic impactor as a means of planetary defense. In late 2022, DART will collide with Dimorphos, the secondary in the Didymos binary asteroid system. The impact will cause a momentum transfer from the spacecraft to the binary asteroid, changing the orbit period of Dimorphos and forcing it to librate in its orbit. Owin… ▽ More The Double Asteroid Redirection Test (DART) mission will be the first test of a kinetic impactor as a means of planetary defense. In late 2022, DART will collide with Dimorphos, the secondary in the Didymos binary asteroid system. The impact will cause a momentum transfer from the spacecraft to the binary asteroid, changing the orbit period of Dimorphos and forcing it to librate in its orbit. Owing to the coupled dynamics in binary asteroid systems, the orbit and libration state of Dimorphos are intertwined. Thus, as the secondary librates, it also experiences fluctuations in its orbit period. These variations in the orbit period are dependent on the magnitude of the impact perturbation, as well as the system's state at impact and the moments of inertia of the secondary. In general, any binary asteroid system whose secondary is librating will have a non-constant orbit period on account of the secondary's fluctuating spin rate. The orbit period variations are typically driven by two modes: a long-period and short-period, each with significant amplitudes on the order of tens of seconds to several minutes. The fluctuating orbit period offers both a challenge and an opportunity in the context of the DART mission. Orbit period oscillations will make determining the post-impact orbit period more difficult, but can also provide information about the system's libration state and the DART impact. △ Less

Submitted 1 December, 2021; originally announced December 2021.

The Excited Spin State of Dimorphos Resulting from the DART Impact

Authors: Harrison F. Agrusa, Ioannis Gkolias, Kleomenis Tsiganis, Derek C. Richardson, Alex J. Meyer, Daniel J. Scheeres, Matija Ćuk, Seth A. Jacobson, Patrick Michel, Özgür Karatekin, Andrew F. Cheng, Masatoshi Hirabayashi, Yun Zhang, Eugene G. Fahnestock, Alex B. Davis

Abstract: The NASA Double Asteroid Redirection Test (DART) mission is a planetary defense-driven test of a kinetic impactor on Dimorphos, the satellite of the binary asteroid 65803 Didymos. DART will intercept Dimorphos at a relative speed of ${\sim}6.5 \text{ km s}^{-1}$, perturbing Dimorphos's orbital velocity and changing the binary orbital period. We present three independent methods (one analytic and t… ▽ More The NASA Double Asteroid Redirection Test (DART) mission is a planetary defense-driven test of a kinetic impactor on Dimorphos, the satellite of the binary asteroid 65803 Didymos. DART will intercept Dimorphos at a relative speed of ${\sim}6.5 \text{ km s}^{-1}$, perturbing Dimorphos's orbital velocity and changing the binary orbital period. We present three independent methods (one analytic and two numerical) to investigate the post-impact attitude stability of Dimorphos as a function of its axial ratios, $a/b$ and $b/c$ ($a \ge b \ge c$), and the momentum transfer efficiency $β$. The first method uses a novel analytic approach in which we assume a circular orbit and a point-mass primary that identifies four fundamental frequencies of motion corresponding to the secondary's mean motion, libration, precession, and nutation frequencies. At resonance locations among these four frequencies, we find that attitude instabilities are possible. Using two independent numerical codes, we recover many of the resonances predicted by the analytic model and indeed show attitude instability. With one code, we use fast Lyapunov indicators to show that the secondary's attitude can evolve chaotically near the resonance locations. Then, using a high-fidelity numerical model, we find that Dimorphos enters a chaotic tumbling state near the resonance locations and is especially prone to unstable rotation about its long axis, which can be confirmed by ESA's Hera mission arriving at Didymos in late 2026. We also show that a fully coupled treatment of the spin and orbital evolution of both bodies is crucial to accurately model the long-term evolution of the secondary's spin state and libration amplitude. Finally, we discuss the implications of a post-impact tumbling or rolling state, including the possibility of terminating BYORP evolution if Dimorphos is no longer in synchronous rotation. △ Less

Submitted 29 July, 2021; v1 submitted 16 July, 2021; originally announced July 2021.

Comments: 38 pages, 13 figures, Accepted for publication in Icarus

Secular resonance sweeping and orbital excitation in decaying disks

Authors: A. Toliou, K. Tsiganis, G. Tsirvoulis

Abstract: We revisit the problem of secular resonance sweeping during the dissipation of a protoplanetary disk and its possible role in exciting the orbits of primordial asteroids, in the light of recent models of solar system evolution. We develop an integrator that incorporates the gravitational effect of a uniformly (or not) depleting, axisymmetric disk with arbitrary surface density profile; its perform… ▽ More We revisit the problem of secular resonance sweeping during the dissipation of a protoplanetary disk and its possible role in exciting the orbits of primordial asteroids, in the light of recent models of solar system evolution. We develop an integrator that incorporates the gravitational effect of a uniformly (or not) depleting, axisymmetric disk with arbitrary surface density profile; its performance is verified by analytical calculations. The secular response of fictitious asteroids, under perturbations from Jupiter, Saturn and a decaying disk, is thoroughly studied. Note that the existence of a symmetry plane induced by the disk, lifts the inherent degeneracy of the two-planet system, such that the '$s_5$' nodal frequency can also play a major role. We examine different resonant configurations for the planets (2:1, 3:2, 5:3), disk models and depletion scenarios. For every case we compute the corresponding time paths of secular resonances, which show when and where resonance crossing occurs. The excitation of asteroids, particularly in inclination, is studied in the various models and compared to analytical estimates. We find that inclination excitation in excess of $\sim 10^{\circ}$ is possible in the asteroid belt, but the nearly uniform spread of $Δi\sim 20^{\circ}$ observed calls for the combined action of secular resonance sweeping with other mechanisms (e.g. scattering by Mars-sized embryos) that would be operating during terrestrial planet formation. Our results are also applicable to extrasolar planetary systems. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Chaotic transport of navigation satellites

Authors: Ioannis Gkolias, Jerome Daquin, Despoina K. Skoulidou, Kleomenis Tsiganis, Christos Efthymiopoulos

Abstract: Navigation satellites are known from numerical studies to reside in a dynamically sensitive environment, which may be of profound importance for their long-term sustainability. We derive the fundamental Hamiltonian of GNSS dynamics and show analytically that near-circular trajectories lie in the neighborhood of a Normally Hyperbolic Invariant Manifold (NHIM), which is the primary source of hyperbo… ▽ More Navigation satellites are known from numerical studies to reside in a dynamically sensitive environment, which may be of profound importance for their long-term sustainability. We derive the fundamental Hamiltonian of GNSS dynamics and show analytically that near-circular trajectories lie in the neighborhood of a Normally Hyperbolic Invariant Manifold (NHIM), which is the primary source of hyperbolicity. Quasi-circular orbits escape through chaotic transport, regulated by the NHIM's stable and unstable manifolds, following a power-law escape time distribution $P(t) \sim t^{-α}$, with $α\sim 0.8 - 1.5$. Our study is highly relevant for the design of satellite disposal trajectories, using manifold dynamics. △ Less

Submitted 25 September, 2019; originally announced September 2019.

Comments: Accepted for publication in Chaos: An Interdisciplinary Journal of Nonlinear Science

Dynamical cartography of Earth satellite orbits

Authors: Aaron J. Rosengren, Despoina K. Skoulidou, Kleomenis Tsiganis, George Voyatzis

Abstract: We have carried out a numerical investigation of the coupled gravitational and non-gravitational perturbations acting on Earth satellite orbits in an extensive grid, covering the whole circumterrestrial space, using an appropriately modified version of the SWIFT symplectic integrator, which is suitable for long-term (120 years) integrations of the non-averaged equations of motion. Hence, we charac… ▽ More We have carried out a numerical investigation of the coupled gravitational and non-gravitational perturbations acting on Earth satellite orbits in an extensive grid, covering the whole circumterrestrial space, using an appropriately modified version of the SWIFT symplectic integrator, which is suitable for long-term (120 years) integrations of the non-averaged equations of motion. Hence, we characterize the long-term dynamics and the phase-space structure of the Earth-orbiter environment, starting from low altitudes (400 km) and going up to the GEO region and beyond. This investigation was done in the framework of the EC-funded "ReDSHIFT" project, with the purpose of enabling the definition of passive debris removal strategies, based on the use of physical mechanisms inherent in the complex dynamics of the problem (i.e., resonances). Accordingly, the complicated interactions among resonances, generated by different perturbing forces (i.e., lunisolar gravity, solar radiation pressure, tesseral harmonics in the geopotential) are accurately depicted in our results, where we can identify the regions of phase space where the motion is regular and long-term stable and regions for which eccentricity growth and even instability due to chaotic behavior can emerge. The results are presented in an "atlas" of dynamical stability maps for different orbital zones, with a particular focus on the (drag-free) range of semimajor axes, where the perturbing effects of the Earth's oblateness and lunisolar gravity are of comparable order. In some regions, the overlapping of the predominant lunisolar secular and semi-secular resonances furnish a number of interesting disposal hatches at moderate to low eccentricity orbits. All computations were repeated for an increased area-to-mass ratio, simulating the case of a satellite equipped with an on-board, area-augmenting device. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: 25 pages, 10 figures

Journal ref: Advances in Space Research, Volume 63, Issue 1, 2019, Pages 443-460

Medium Earth Orbit dynamical survey and its use in passive debris removal

Authors: Despoina K. Skoulidou, Aaron J. Rosengren, Kleomenis Tsiganis, George Voyatzis

Abstract: The Medium Earth Orbit (MEO) region hosts satellites for navigation, communication, and geodetic/space environmental science, among which are the Global Navigation Satellites Systems (GNSS). Safe and efficient removal of debris from MEO is problematic due to the high cost for maneuvers needed to directly reach the Earth (reentry orbits) and the relatively crowded GNSS neighborhood (graveyard orbit… ▽ More The Medium Earth Orbit (MEO) region hosts satellites for navigation, communication, and geodetic/space environmental science, among which are the Global Navigation Satellites Systems (GNSS). Safe and efficient removal of debris from MEO is problematic due to the high cost for maneuvers needed to directly reach the Earth (reentry orbits) and the relatively crowded GNSS neighborhood (graveyard orbits). Recent studies have highlighted the complicated secular dynamics in the MEO region, but also the possibility of exploiting these dynamics, for designing removal strategies. In this paper, we present our numerical exploration of the long-term dynamics in MEO, performed with the purpose of unveiling the set of reentry and graveyard solutions that could be reached with maneuvers of reasonable DV cost. We simulated the dynamics over 120-200 years for an extended grid of millions of fictitious MEO satellites that covered all inclinations from 0 to 90deg, using non-averaged equations of motion and a suitable dynamical model that accounted for the principal geopotential terms, 3rd-body perturbations and solar radiation pressure (SRP). We found a sizeable set of usable solutions with reentry times that exceed ~40years, mainly around three specific inclination values: 46deg, 56deg, and 68deg; a result compatible with our understanding of MEO secular dynamics. For DV <= 300 m/s (i.e., achieved if you start from a typical GNSS orbit and target a disposal orbit with e<0.3), reentry times from GNSS altitudes exceed ~70 years, while low-cost (DV ~= 5-35 m/s) graveyard orbits, stable for at lest 200 years, are found for eccentricities up to e~0.018. This investigation was carried out in the framework of the EC-funded "ReDSHIFT" project. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: 39 pages, 23 figures

Dynamical Lifetime Survey of Geostationary Transfer Orbits

Authors: Despoina K. Skoulidou, Aaron J. Rosengren, Kleomenis Tsiganis, George Voyatzis

Abstract: In this paper, we study the long-term dynamical evolution of highly-elliptical orbits (HEOs) in the medium-Earth orbit (MEO) region around the Earth. The real population consists primarily of Geosynchronous Transfer Orbits (GTOs), launched at specific inclinations, Molniya-type satellites and related debris. We performed a suite of long-term numerical integrations (up to 200 years) within a realis… ▽ More In this paper, we study the long-term dynamical evolution of highly-elliptical orbits (HEOs) in the medium-Earth orbit (MEO) region around the Earth. The real population consists primarily of Geosynchronous Transfer Orbits (GTOs), launched at specific inclinations, Molniya-type satellites and related debris. We performed a suite of long-term numerical integrations (up to 200 years) within a realistic dynamical model, aimed primarily at recording the dynamical lifetime of such orbits (defined as the time needed for atmospheric reentry) and understanding its dependence on initial conditions and other parameters, such as the area-to-mass ratio (A/m). Our results are presented in the form of 2-D lifetime maps, for different values of inclination, A/m, and drag coefficient. We find that the majority of small debris (> 70%, depending on the inclination) can naturally reenter within 25-90 years, but these numbers are significantly less optimistic for large debris (e.g., upper stages), with the notable exception of those launched from high latitude (Baikonur). We estimate the reentry probability and mean dynamical lifetime for different classes of GTOs and we find that both quantities depend primarily and strongly on initial perigee altitude. Atmospheric drag and higher A/m values extend the reentry zones, especially at low inclinations. For high inclinations, this dependence is weakened, as the primary mechanisms leading to reentry are overlapping lunisolar resonances. This study forms part of the EC-funded (H2020) "ReDSHIFT" project. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: 16 pages, 10 figures

Journal ref: Celest Mech Dyn Astr (2018) 130: 77

NELIOTA: The wide-field, high-cadence lunar monitoring system at the prime focus of the Kryoneri telescope

Authors: E. M. Xilouris, A. Z. Bonanos, I. Bellas-Velidis, P. Boumis, A. Dapergolas, A. Maroussis, A. Liakos, I. Alikakos, V. Charmandaris, G. Dimou, A. Fytsilis, M. Kelley, D. Koschny, V. Navarro, K. Tsiganis, K. Tsinganos

Abstract: We present the technical specifications and first results of the ESA-funded, lunar monitoring project "NELIOTA" (NEO Lunar Impacts and Optical TrAnsients) at the National Observatory of Athens, which aims to determine the size-frequency distribution of small Near-Earth Objects (NEOs) via detection of impact flashes on the surface of the Moon. For the purposes of this project a twin camera instrume… ▽ More We present the technical specifications and first results of the ESA-funded, lunar monitoring project "NELIOTA" (NEO Lunar Impacts and Optical TrAnsients) at the National Observatory of Athens, which aims to determine the size-frequency distribution of small Near-Earth Objects (NEOs) via detection of impact flashes on the surface of the Moon. For the purposes of this project a twin camera instrument was specially designed and installed at the 1.2 m Kryoneri telescope utilizing the fast-frame capabilities of scientific Complementary Metal-Oxide Semiconductor detectors (sCMOS). The system provides a wide field-of-view (17.0' $\times$ 14.4') and simultaneous observations in two photometric bands (R and I), reaching limiting magnitudes of 18.7 mag in 10 sec in both bands at a 2.5 signal-to-noise level. This makes it a unique instrument that can be used for the detection of NEO impacts on the Moon, as well as for any astronomy projects that demand high-cadence multicolor observations. The wide field-of-view ensures that a large portion of the Moon is observed, while the simultaneous, high-cadence, monitoring in two photometric bands makes possible, for the first time, the determination of the temperatures of the impacts on the Moon's surface and the validation of the impact flashes from a single site. Considering the varying background level on the Moon's surface we demonstrate that the NELIOTA system can detect NEO impact flashes at a 2.5 signal-to-noise level of ~12.4 mag in the I-band and R-band for observations made at low lunar phases ~0.1. We report 31 NEO impact flashes detected during the first year of the NELIOTA campaign. The faintest flash was at 11.24 mag in the R-band (about two magnitudes fainter than ever observed before) at lunar phase 0.32. Our observations suggest a detection rate of $1.96 \times 10^{-7}$ events $km^{-2} h^{-1}$. △ Less

Submitted 3 September, 2018; originally announced September 2018.

Comments: Accepted for publication in A&A

Inclined asymmetric librations in exterior resonances

Authors: G. Voyatzis, K. Tsiganis, K. I. Antoniadou

Abstract: Librational motion in celestial mechanics is generally associated with the existence of stable resonant configurations and signified by the existence of stable periodic solutions and oscillation of critical (resonant) angles. When such an oscillation takes place around a value different than 0 or $π$, the libration is called asymmetric. In the context of the planar circular restricted three-body p… ▽ More Librational motion in celestial mechanics is generally associated with the existence of stable resonant configurations and signified by the existence of stable periodic solutions and oscillation of critical (resonant) angles. When such an oscillation takes place around a value different than 0 or $π$, the libration is called asymmetric. In the context of the planar circular restricted three-body problem (CRTBP), asymmetric librations have been identified for the exterior mean-motion resonances (MMRs) 1:2, 1:3 etc. as well as for co-orbital motion (1:1). In exterior MMRs the massless body is the outer one. In this paper, we study asymmetric librations in the 3-dimensional space. We employ the computational approach of Markellos (1978) and compute families of asymmetric periodic orbits and their stability. Stable, asymmetric periodic orbits are surrounded in phase space by domains of initial conditions which correspond to stable evolution and librating resonant angles. Our computations were focused on the spatial circular restricted three-body model of the Sun-Neptune-TNO system (TNO= trans-Neptunian object). We compare our results with numerical integrations of observed TNOs, which reveal that some of them perform 1:2-resonant, inclined asymmetric librations. For the stable 1:2 TNOs librators, we find that their libration seems to be related with the vertically stable planar asymmetric orbits of our model, rather than the 3-dimensional ones found in the present study. △ Less

Submitted 1 March, 2018; originally announced March 2018.

Comments: Accepted for publication in CeMDA

The rectilinear three body problem as a basis for studying highly-eccentric systems

Authors: George Voyatzis, Kleomenis Tsiganis, Michalis Gaitanas

Abstract: The rectilinear elliptic restricted Three Body Problem (TBP) is the limiting case of the elliptic restricted TBP when the motion of the primaries is described by a Keplerian ellipse with eccentricity $e'=1$, but the collision of the primaries is assumed to be a non-singular point. The rectilinear model has been proposed as a starting model for studying the dynamics of motion around highly eccentri… ▽ More The rectilinear elliptic restricted Three Body Problem (TBP) is the limiting case of the elliptic restricted TBP when the motion of the primaries is described by a Keplerian ellipse with eccentricity $e'=1$, but the collision of the primaries is assumed to be a non-singular point. The rectilinear model has been proposed as a starting model for studying the dynamics of motion around highly eccentric binary systems. Broucke (1969) explored the rectilinear problem and obtained isolated periodic orbits for mass parameter $μ=0.5$ (equal masses of the primaries). We found that all orbits obtained by Broucke are linearly unstable. We extend Broucke's computations by using a finer search for symmetric periodic orbits and computing their linear stability. We found a large number of periodic orbits, but only eight of them were found to be linearly stable and are associated with particular mean motion resonances. These stable orbits are used as generating orbits for continuation with respect to $μ$ and $e'<1$. Also, continuation of periodic solutions with respect to the mass of the small body can be applied by using the general TBP. FLI maps of dynamical stability show that stable periodic orbits are surrounded in phase space with regions of regular orbits indicating that systems of very highly eccentric orbits can be found in stable resonant configurations. As an application we present a stability study for the planetary system HD7449. △ Less

Submitted 16 August, 2017; originally announced August 2017.

Comments: to appear in CMDA

Reconstructing the size distribution of the primordial Main Belt

Authors: Georgios Tsirvoulis, Alessandro Morbidelli, Marco Delbo, Kleomenis Tsiganis

Abstract: In this work we aim to constrain the slope of the size distribution of main-belt asteroids, at their primordial state. To do so we turn out attention to the part of the main asteroid belt between 2.82 and 2.96~AU, the so-called "pristine zone", which has a low number density of asteroids and few, well separated asteroid families. Exploiting these unique characteristics, and using a modified versio… ▽ More In this work we aim to constrain the slope of the size distribution of main-belt asteroids, at their primordial state. To do so we turn out attention to the part of the main asteroid belt between 2.82 and 2.96~AU, the so-called "pristine zone", which has a low number density of asteroids and few, well separated asteroid families. Exploiting these unique characteristics, and using a modified version of the hierarchical clustering method we are able to remove the majority of asteroid family members from the region. The remaining, background asteroids should be of primordial origin, as the strong 5/2 and 7/3 mean-motion resonances with Jupiter inhibit transfer of asteroids to and from the neighboring regions. The size-frequency distribution of asteroids in the size range $17<D(\rm{km})<70$ has a slope $q\simeq-1$. Using Monte-Carlo methods, we are able to simulate, and compensate for the collisional and dynamical evolution of the asteroid population, and get an upper bound for its size distribution slope $q=-1.43$. In addition, applying the same 'family extraction' method to the neighboring regions, i.e. the middle and outer belts, and comparing the size distributions of the respective background populations, we find statistical evidence that no large asteroid families of primordial origin had formed in the middle or pristine zones. △ Less

Submitted 7 June, 2017; originally announced June 2017.

Magnitude and timing of the giant planet instability: A reassessment from the perspective of the asteroid belt

Authors: Athanasia Toliou, Alessandro Morbidelli, Kleomenis Tsiganis

Abstract: It is generally accepted today that our solar system has undergone a phase during which the orbits of the giant planets became very unstable. In recent years, many studies have identified traces of this event and have provided reasonable justification for this occurrence. The magnitude (in terms of orbital variation) and the timing of the instability though (early or late with respect to the dispe… ▽ More It is generally accepted today that our solar system has undergone a phase during which the orbits of the giant planets became very unstable. In recent years, many studies have identified traces of this event and have provided reasonable justification for this occurrence. The magnitude (in terms of orbital variation) and the timing of the instability though (early or late with respect to the dispersal of the gas disk) still remains an open debate. The terrestrial planets seem to set a strict constraint: either the giant planet instability happened early, while the terrestrial planets were still forming, or the orbits of Jupiter and Saturn had to separate from each other impulsively, with a large enough `jump' in semimajor axis (Brasser et al. 2009; Kaib and Chambers 2016) for the terrestrial planets to remain stable. Because a large orbital jump is a low probability event, the early instability hypothesis seems to be favored. However, the asteroid belt would also evolve in a different way, assuming different instability amplitudes. These two constraints need to match each other in order to favor one scenario over the other. Considering an initially dynamically cold disk of asteroids, Morbidelli et al. (2010) concluded that a comparably large jump is needed to reconstruct the current asteroid belt. Here we confirm the same conclusion, but considering an asteroid population already strongly excited in eccentricity, such as that produced in the Grand Tack scenario (Walsh et al. 2011). Because the asteroids existed since the time of removal of the gas disk, unlike the terrestrial planets, this constraint on the width of the giant planet jump is valid regardless of whether the instability happened early or late. Hence, at this stage, assuming an early instability does not appear to provide any advantage in terms of the probabilistic reconstruction of the solar system structure. △ Less

Submitted 14 June, 2016; originally announced June 2016.

Journal ref: A&A 592, A72 (2016)

Diffusive chaos in navigation satellites orbits

Authors: J. Daquin, A. J. Rosengren, K. Tsiganis

Abstract: The navigation satellite constellations in medium-Earth orbit exist in a background of third-body secular resonances stemming from the perturbing gravitational effects of the Moon and the Sun. The resulting chaotic motions, emanating from the overlapping of neighboring resonant harmonics, induce especially strong perturbations on the orbital eccentricity, which can be transported to large values,… ▽ More The navigation satellite constellations in medium-Earth orbit exist in a background of third-body secular resonances stemming from the perturbing gravitational effects of the Moon and the Sun. The resulting chaotic motions, emanating from the overlapping of neighboring resonant harmonics, induce especially strong perturbations on the orbital eccentricity, which can be transported to large values, thereby increasing the collision risk to the constellations and possibly leading to a proliferation of space debris. We show here that this transport is of a diffusive nature and we present representative diffusion maps that are useful in obtaining a global comprehension of the dynamical structure of the navigation satellite orbits. △ Less

Submitted 1 June, 2016; originally announced June 2016.

Comments: 8 pages, 5 figures, conference Chaos, complexity and transport (Marseille, France)

NELIOTA: ESA's new NEO lunar impact monitoring project with the 1.2m telescope at the National Observatory of Athens

Authors: A. Z. Bonanos, M. Xilouris, P. Boumis, I. Bellas-Velidis, A. Maroussis, A. Dapergolas, A. Fytsilis, V. Charmandaris, K. Tsiganis, K. Tsinganos

Abstract: NELIOTA is a new ESA activity launched at the National Observatory of Athens in February 2015 aiming to determine the distribution and frequency of small near-earth objects (NEOs) via lunar monitoring. The project involves upgrading the 1.2m Kryoneri telescope at the National Observatory of Athens, procuring two fast-frame cameras, and developing a software system, which will control the telescope… ▽ More NELIOTA is a new ESA activity launched at the National Observatory of Athens in February 2015 aiming to determine the distribution and frequency of small near-earth objects (NEOs) via lunar monitoring. The project involves upgrading the 1.2m Kryoneri telescope at the National Observatory of Athens, procuring two fast-frame cameras, and developing a software system, which will control the telescope and the cameras, process the images and automatically detect NEO impacts. NELIOTA will provide a web-based user interface, where the impact events will be reported and made available to the scientific community and the general public. The objective of this 3.5 year activity is to design, develop and implement a highly automated lunar monitoring system, which will conduct an observing campaign for 2 years in search of NEO impact flashes on the Moon. The impact events will be verified, characterised and reported. The 1.2m telescope will be capable of detecting flashes much fainter than current, small-aperture, lunar monitoring telescopes. NELIOTA is therefore expected to characterise the frequency and distribution of NEOs weighing as little as a few grams. △ Less

Submitted 13 October, 2015; originally announced October 2015.

Comments: 3 pages, 1 figure, IAU GA Symposium 318 proceedings

The origin of long-lived asteroids in the 2:1 mean-motion resonance with Jupiter

Authors: O. Chrenko, M. Brož, D. Nesvorný, K. Tsiganis, D. K. Skoulidou

Abstract: The 2:1 mean-motion resonance with Jupiter harbours two distinct groups of asteroids. The short-lived population is known to be a transient group sustained in steady state by the Yarkovsky semimajor axis drift. The long-lived asteroids, however, can exhibit dynamical lifetimes comparable to $4\,\mathrm{Gyr}$. They reside near two isolated islands of the phase space denoted $\mathrm{A}$ and… ▽ More The 2:1 mean-motion resonance with Jupiter harbours two distinct groups of asteroids. The short-lived population is known to be a transient group sustained in steady state by the Yarkovsky semimajor axis drift. The long-lived asteroids, however, can exhibit dynamical lifetimes comparable to $4\,\mathrm{Gyr}$. They reside near two isolated islands of the phase space denoted $\mathrm{A}$ and $\mathrm{B}$, with an uneven population ratio $\mathrm{B}/\mathrm{A} \simeq 10$. The orbits of $\mathrm{A}$-island asteroids are predominantly highly inclined, compared to island $\mathrm{B}$. The size-frequency distribution is steep but the orbital distribution lacks any evidence of a collisional cluster. These observational constraints are somewhat puzzling and therefore the origin of the long-lived asteroids has not been explained so far. With the aim to provide a viable explanation, we first update the resonant population and revisit its physical properties. Using an $N$-body model with seven planets and the Yarkovsky effect included, we demonstrate that the dynamical depletion of island $\mathrm{A}$ is faster, in comparison with island $\mathrm{B}$. Then we investigate (i) the survivability of primordial resonant asteroids and (ii) capture of the population during planetary migration, following a recently described scenario with an escaping fifth giant planet and a jumping-Jupiter instability. We also model the collisional evolution of the resonant population over past $4\,\mathrm{Gyr}$. Our conclusion is that the long-lived group was created by resonant capture from a narrow part of hypothetical outer main-belt family during planetary migration. Primordial asteroids surviving the migration were probably not numerous enough to substantially contribute to the observed population. △ Less

Submitted 16 May, 2015; originally announced May 2015.

Comments: Accepted for publication in MNRAS. 20 pages, 17 figures

Vertical instability and inclination excitation during planetary migration

Authors: G. Voyatzis, K. I. Antoniadou, K. Tsiganis

Abstract: We consider a two-planet system, which migrates under the influence of dissipative forces that mimic the effects of gas-driven (Type II) migration. It has been shown that, in the planar case, migration leads to resonant capture after an evolution that forces the system to follow families of periodic orbits. Starting with planets that differ slightly from a coplanar configuration, capture can, also… ▽ More We consider a two-planet system, which migrates under the influence of dissipative forces that mimic the effects of gas-driven (Type II) migration. It has been shown that, in the planar case, migration leads to resonant capture after an evolution that forces the system to follow families of periodic orbits. Starting with planets that differ slightly from a coplanar configuration, capture can, also, occur and, additionally, excitation of planetary inclinations has been observed in some cases. We show that excitation of inclinations occurs, when the planar families of periodic orbits, which are followed during the initial stages of planetary migration, become vertically unstable. At these points, {\em vertical critical orbits} may give rise to generating stable families of $3D$ periodic orbits, which drive the evolution of the migrating planets to non-coplanar motion. We have computed and present here the vertical critical orbits of the $2/1$ and $3/1$ resonances, for various values of the planetary mass ratio. Moreover, we determine the limiting values of eccentricity for which the "inclination resonance" occurs. △ Less

Submitted 14 July, 2014; v1 submitted 28 June, 2014; originally announced June 2014.

Comments: Accepted for publication in Celestial Mechanics and Dynamical Astronomy

Journal ref: Celestial Mechanics and Dynamical Astronomy, Volume 119, Issue 3-4, pp. 221-235 (2014)

NOTE: Explaining why the Uranian satellites have equatorial prograde orbits despite the large planetary obliquity

Authors: Alessandro Morbidelli, Kleomenis Tsiganis, Konstantin Batygin, Aurelien Crida, Rodney Gomes

Abstract: We show that the existence of prograde equatorial satellites is consistent with a collisional tilting scenario for Uranus. In fact, if the planet was surrounded by a proto-satellite disk at the time of the tilting and a massive ring of material was temporarily placed inside the Roche radius of the planet by the collision, the proto-satellite disk would have started to precess incoherently around t… ▽ More We show that the existence of prograde equatorial satellites is consistent with a collisional tilting scenario for Uranus. In fact, if the planet was surrounded by a proto-satellite disk at the time of the tilting and a massive ring of material was temporarily placed inside the Roche radius of the planet by the collision, the proto-satellite disk would have started to precess incoherently around the equator of the planet, up to a distance greater than that of Oberon. Collisional damping would then have collapsed it into a thin equatorial disk, from which the satellites eventually formed. The fact that the orbits of the satellites are prograde requires Uranus to have had a non-negligible initial obliquity (comparable to that of Neptune) before it was finally tilted to 98 degrees. △ Less

Submitted 23 August, 2012; originally announced August 2012.

Journal ref: Icarus, 219, 737 (2012)

Trapping in three-planet resonances during gas-driven migration

Authors: A. -S. Libert, K. Tsiganis

Abstract: We study the establishment of three-planet resonances -similar to the Laplace resonance in the Galilean satellites- and their effects on the mutual inclinations of the orbital planes of the planets, assuming that the latter undergo migration in a gaseous disc. In particular, we examine the resonance relations that occur, by varying the physical and initial orbital parameters of the planets (mass,… ▽ More We study the establishment of three-planet resonances -similar to the Laplace resonance in the Galilean satellites- and their effects on the mutual inclinations of the orbital planes of the planets, assuming that the latter undergo migration in a gaseous disc. In particular, we examine the resonance relations that occur, by varying the physical and initial orbital parameters of the planets (mass, initial semi-major axis and eccentricity) as well as the parameters of the migration forces (migration rate and eccentricity damping rate), which are modeled here through a simplified analytic prescription. We find that, in general, for planetary masses below 1.5 M_J, multiple-planet resonances of the form n3:n2:n1=1:2:4 and 1:3:6 are established, as the inner planets, m1 and m2, get trapped in a 1:2 resonance and the outer planet m3 subsequently is captured in a 1:2 or 1:3 resonance with m2. For mild eccentricity damping, the resonance pumps the eccentricities of all planets on a relatively short time-scale, to the point where they enter an inclination-type resonance (as in Libert & Tsiganis 2011); then mutual inclinations can grow to ~35°, thus forming a "3-D system". On the other hand, we find that trapping of m2 in a 2:3 resonance with m1 occurs very rarely, for the range of masses used here, so only two cases of capture in a respective three-planet resonance were found. Our results suggest that trapping in a three-planet resonance can be common in exoplanetary systems, provided that the planets are not very massive. Inclination pumping could then occur relatively fast, provided that eccentricity damping is not very efficient so that at least one of the inner planets acquires an orbital eccentricity higher than e=0.3. △ Less

Submitted 13 March, 2012; originally announced March 2012.

Comments: 19 pages, 9 figures, published in Celestial Mechanics and Dynamical Astronomy

Journal ref: Celestial Mechanics and Dynamical Astronomy, 111, 201-218 (2011)

Formation of '3D' multiplanet systems by dynamical disruption of multiple-resonance configurations

Authors: A. -S. Libert, K. Tsiganis

Abstract: Assuming that giant planets are formed in thin protoplanetary discs, a '3D' system can form, provided that the mutual inclination is excited by some dynamical mechanism. Resonant interactions and close planetary encounters are thought to be the primary inclination-excitation mechanisms, resulting in a resonant and non-resonant system, respectively. Here we propose an alternative formation scenario… ▽ More Assuming that giant planets are formed in thin protoplanetary discs, a '3D' system can form, provided that the mutual inclination is excited by some dynamical mechanism. Resonant interactions and close planetary encounters are thought to be the primary inclination-excitation mechanisms, resulting in a resonant and non-resonant system, respectively. Here we propose an alternative formation scenario, starting from a system composed of three giant planets in a nearly coplanar configuration. As was recently shown for the case of the Solar system, planetary migration in the gas disc (Type II migration) can force the planets to become trapped in a multiply resonant state. We simulate this process, assuming different values for the planetary masses and mass ratios. We show that such a triple resonance generally becomes unstable as the resonance excites the eccentricities of all planets and planet-planet scattering sets in. One of the three planets is typically ejected from the system, leaving behind a dynamically 'hot' (but stable) two-planet configuration. The resulting two-planet systems typically have large values of semimajor axial ratios (a1/a2 < 0.3), while the mutual inclination can be as high as 70°, with a median of \sim30°. A small fraction of our two-planet systems (\sim5 per cent) ends up in the stability zone of the Kozai resonance. In a few cases, the triple resonance can remain stable for long times and a '3D' system can form by resonant excitation of the orbital inclinations; such a three-planet system could be stable if enough eccentricity damping is exerted on the planets. Finally, in the single-planet resulting systems, which are formed when two planets are ejected from the system, the inclination of the planet's orbital plane with respect to the initial invariant plane -presumably the plane perpendicular to the star's spin axis- can be as large as \sim40°. △ Less

Submitted 13 March, 2012; originally announced March 2012.

Comments: 9 pages, 5 figures, published in MNRAS

Journal ref: MNRAS 412, 2353-2360 (2011)

Interaction of free-floating planets with a star-planet pair

Authors: Harry Varvoglis, Vasiliki Sgardeli, Kleomenis Tsiganis

Abstract: The recent discovery of free-floating planets and their theoretical interpretation as celestial bodies, either condensed independently or ejected from parent stars in tight clusters, introduced an intriguing possibility. Namely, that some exoplanets are not condensed from the protoplanetary disk of their parent star. In this novel scenario a free-floating planet interacts with an already existing… ▽ More The recent discovery of free-floating planets and their theoretical interpretation as celestial bodies, either condensed independently or ejected from parent stars in tight clusters, introduced an intriguing possibility. Namely, that some exoplanets are not condensed from the protoplanetary disk of their parent star. In this novel scenario a free-floating planet interacts with an already existing planetary system, created in a tight cluster, and is captured as a new planet. In the present work we study this interaction process by integrating trajectories of planet-sized bodies, which encounter a binary system consisting of a Jupiter-sized planet revolving around a Sun-like star. To simplify the problem we assume coplanar orbits for the bound and the free-floating planet and an initially parabolic orbit for the free-floating planet. By calculating the uncertainty exponent, a quantity that measures the dependence of the final state of the system on small changes of the initial conditions, we show that the interaction process is a fractal classical scattering. The uncertainty exponent is in the range (0.2-0.3) and is a decreasing function of time. In this way we see that the statistical approach we follow to tackle the problem is justified. The possible final outcomes of this interaction are only four, namely flyby, planet exchange, capture or disruption. We give the probability of each outcome as a function of the incoming planet's mass. We find that the probability of exchange or capture (in prograde as well as retrograde orbits and for very long times) is non-negligible, a fact that might explain the possible future observations of planetary systems with orbits that are either retrograde or tight and highly eccentric. △ Less

Submitted 17 July, 2012; v1 submitted 6 January, 2012; originally announced January 2012.

Comments: 19 pages, 12 figures

Formation and Evolution of Planetary Systems in Presence of Highly Inclined Stellar Perturbers

Authors: Konstantin Batygin, Alessandro Morbidelli, Kleomenis Tsiganis

Abstract: The presence of highly eccentric extrasolar planets in binary stellar systems suggests that the Kozai effect has played an important role in shaping their dynamical architectures. However, the formation of planets in inclined binary systems poses a considerable theoretical challenge, as orbital excitation due to the Kozai resonance implies destructive, high-velocity collisions among planetesimals.… ▽ More The presence of highly eccentric extrasolar planets in binary stellar systems suggests that the Kozai effect has played an important role in shaping their dynamical architectures. However, the formation of planets in inclined binary systems poses a considerable theoretical challenge, as orbital excitation due to the Kozai resonance implies destructive, high-velocity collisions among planetesimals. To resolve the apparent difficulties posed by Kozai resonance, we seek to identify the primary physical processes responsible for inhibiting the action of Kozai cycles in protoplanetary disks. Subsequently, we seek to understand how newly-formed planetary systems transition to their observed, Kozai-dominated dynamical states. We find that theoretical difficulties in planet formation arising from the presence of a distant companion star, posed by the Kozai effect and other secular perturbations, can be overcome by a proper account of gravitational interactions within the protoplanetary disk. In particular, fast apsidal recession induced by disk self-gravity tends to erase the Kozai effect, and ensure that the disk's unwarped, rigid structure is maintained. Subsequently, once a planetary system has formed, the Kozai effect can continue to be wiped out as a result of apsidal precession, arising from planet-planet interactions. However, if such a system undergoes a dynamical instability, its architecture may change in such a way that the Kozai effect becomes operative. The results presented here suggest that planetary formation in highly inclined binary systems is not stalled by perturbations, arising from the stellar companion. Consequently, planet formation in binary stars is probably no different from that around single stars on a qualitative level. Furthermore, it is likely that systems where the Kozai effect operates, underwent a transient phase of dynamical instability in the past. △ Less

Submitted 20 June, 2011; originally announced June 2011.

Comments: 9 pages, 7 figures, accepted for publication in Astronomy and Astrophysics

Evidence from the asteroid belt for a violent past evolution of Jupiter's orbit

Authors: Alessandro Morbidelli, Ramon Brasser, Rodney Gomes, Harold F. Levison, Kleomenis Tsiganis

Abstract: We use the current orbital structure of large (>50km) asteroids in the main asteroid belt to constrain the evolution of the giant planets when they migrated from their primordial orbits to their current ones. Minton & Malhotra (2009) showed that the orbital distribution of large asteroids in the main belt can be reproduced by an exponentially-decaying migration of the giant planets on a time scale… ▽ More We use the current orbital structure of large (>50km) asteroids in the main asteroid belt to constrain the evolution of the giant planets when they migrated from their primordial orbits to their current ones. Minton & Malhotra (2009) showed that the orbital distribution of large asteroids in the main belt can be reproduced by an exponentially-decaying migration of the giant planets on a time scale of tau ~ 0.5My. However, self-consistent numerical simulations show that the planetesimal-driven migration of the giant planets is inconsistent with an exponential change in their semi major axes on such a short time scale (Hahn & Malhotra, 1999). In fact, the typical time scale is tau > 5My. When giant planet migration on this time scale is applied to the asteroid belt, the resulting orbital distribution is incompatible with the observed one. However, the planet migration can be significantly sped up by planet-planet encounters. Consider an evolution where both Jupiter and Saturn have close encounters with a Neptune-mass planet (presumably Uranus or Neptune themselves) and where this third planet, after being scattered inwards by Saturn, is scattered outwards by Jupiter. This scenario leads to a very rapid increase in the orbital separation between Jupiter and Saturn that we show here to have only mild effects on the structure of asteroid belt. This type of evolution is called a jumping-Jupiter case. Our results suggest that the total mass and dynamical excitation of the asteroid belt before migration were comparable to those currently observed. Moreover, they imply that, before migration, the orbits of Jupiter and Saturn were much less eccentric than the current ones. △ Less

Submitted 8 September, 2010; originally announced September 2010.

Comments: AJ, in press

Chaotic Transport and Chronology of Complex Asteroid Families

Authors: Bojan Novakovic, Kleomenis Tsiganis, Zoran Knezevic

Abstract: We present a transport model that describes the orbital diffusion of asteroids in chaotic regions of the 3-D space of proper elements. Our goal is to use a simple random-walk model to study the evolution and derive accurate age estimates for dynamically complex asteroid families. To this purpose, we first compute local diffusion coefficients, which characterize chaotic diffusion in proper eccent… ▽ More We present a transport model that describes the orbital diffusion of asteroids in chaotic regions of the 3-D space of proper elements. Our goal is to use a simple random-walk model to study the evolution and derive accurate age estimates for dynamically complex asteroid families. To this purpose, we first compute local diffusion coefficients, which characterize chaotic diffusion in proper eccentricity (e_p) and inclination (I_p), in a selected phase-space region. Then, a Monte-Carlo-type code is constructed and used to track the evolution of random walkers (i.e. asteroids), by coupling diffusion in (e_p,I_p) with a drift in proper semi-major axis (a_p) induced by the Yarkovsky/YORP thermal effects. We validate our model by applying it to the family of (490) Veritas, for which we recover previous estimates of its age (~8.3 Myr). Moreover, we show that the spreading of chaotic family members in proper elements space is well reproduced in our random-walk simulations. Finally, we apply our model to the family of (3556) Lixiaohua, which is much older than Veritas and thus much more affected by thermal forces. We find the age of the Lixiaohua family to be 155+/-36 Myr. △ Less

Submitted 30 October, 2009; originally announced October 2009.

Comments: 11 pages, 11 figures, Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2010, Vol. 402, Issue 2, pp. 1263-1272

Constructing the secular architecture of the solar system II: The terrestrial planets

Authors: Ramon Brasser, Alessandro Morbidelli, Rodney Gomes, Kleomenis Tsiganis, Harold F. Levison

Abstract: We investigate the dynamical evolution of the terrestrial planets during the planetesimal-driven migration of the giant planets. A basic assumption of this work is that giant planet migration occurred after the completion of terrestrial planet formation, such as in the models that link the former to the origin of the Late Heavy Bombardment. The divergent migration of Jupiter and Saturn causes th… ▽ More We investigate the dynamical evolution of the terrestrial planets during the planetesimal-driven migration of the giant planets. A basic assumption of this work is that giant planet migration occurred after the completion of terrestrial planet formation, such as in the models that link the former to the origin of the Late Heavy Bombardment. The divergent migration of Jupiter and Saturn causes the g5 eigenfrequency to cross resonances of the form g5=gk with k ranging from 1 to 4. Consequently these secular resonances cause large-amplitude responses in the eccentricities of the terrestrial planets. We show that the resonances g5=g_4 and g5=g3 do not pose a problem if Jupiter and Saturn have a fast approach and departure from their mutual 2:1 mean motion resonance. On the other hand, the resonance crossings g5=g2 and g5=g1 are more of a concern as they tend to yield a terrestrial system incompatible with the current one. We offer two solutions to this problem. The first uses the fact that a secular resonance crossing can also damp the amplitude of a Fourier mode if the latter is large originally. A second scenario involves a 'jumping Jupiter' in which encounters between an ice giant and Jupiter, without ejection of the former, cause the latter to migrate away from Saturn much faster than if migration is driven solely by encounters with planetesimals. In this case, the g5=g2 and g5=g1 resonances can be jumped over, or occur very briefly. △ Less

Submitted 10 September, 2009; originally announced September 2009.

Comments: Astronomy & Astrophysics (2009) in press

Constructing the secular architecture of the solar system I: The giant planets

Authors: Morbidelli Alessandro, Ramon Brasser, Kleomenis Tsiganis, Rodney Gomes, Harold F. Levison

Abstract: Using numerical simulations, we show that smooth migration of the giant planets through a planetesimal disk leads to an orbital architecture that is inconsistent with the current one: the resulting eccentricities and inclinations of their orbits are too small. The crossing of mutual mean motion resonances by the planets would excite their orbital eccentricities but not their orbital inclinations… ▽ More Using numerical simulations, we show that smooth migration of the giant planets through a planetesimal disk leads to an orbital architecture that is inconsistent with the current one: the resulting eccentricities and inclinations of their orbits are too small. The crossing of mutual mean motion resonances by the planets would excite their orbital eccentricities but not their orbital inclinations. Moreover, the amplitudes of the eigenmodes characterising the current secular evolution of the eccentricities of Jupiter and Saturn would not be reproduced correctly; only one eigenmode is excited by resonance-crossing. We show that, at the very least, encounters between Saturn and one of the ice giants (Uranus or Neptune) need to have occurred, in order to reproduce the current secular properties of the giant planets, in particular the amplitude of the two strongest eigenmodes in the eccentricities of Jupiter and Saturn. △ Less

Submitted 10 September, 2009; originally announced September 2009.

Comments: Astronomy & Astrophysics (2009) in press

Origin of the Structure of the Kuiper Belt during a Dynamical Instability in the Orbits of Uranus and Neptune

Authors: Harold F. Levison, Alessandro Morbidelli, Christa Van Laerhoven, Rodney Gomes, Kleomenis Tsiganis

Abstract: We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. One characteristic of this model is that the proto-planetary disk must have been tru… ▽ More We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. One characteristic of this model is that the proto-planetary disk must have been truncated at roughly 30 to 35 AU so that Neptune would stop migrating at its currently observed location. As a result, the Kuiper belt would have initially been empty. In this paper we present a new dynamical mechanism which can deliver objects from the region interior to ~35 AU to the Kuiper belt without excessive inclination excitation. Assuming that the last encounter with Uranus delivered Neptune onto a low-inclination orbit with a semi-major axis of ~27 AU and an eccentricity of ~0.3, and that subsequently Neptune's eccentricity damped in ~1 My, our simulations reproduce the main observed properties of the Kuiper belt at an unprecedented level. △ Less

Submitted 4 December, 2007; originally announced December 2007.

Dynamics of the giant planets of the solar system in the gaseous proto-planetary disk and relationship to the current orbital architecture

Authors: Alessandro Morbidelli, Kleomenis Tsiganis, Aurelien Crida, Harold F. Levison, Rodney Gomes

Abstract: We study the orbital evolution of the 4 giant planets of our solar system in a gas disk. Our investigation extends the previous works by Masset and Snellgrove (2001) and Morbidelli and Crida (2007, MC07), which focussed on the dynamics of the Jupiter-Saturn system. The only systems that we found to reach a steady state are those in which the planets are locked in a quadruple mean motion resonanc… ▽ More We study the orbital evolution of the 4 giant planets of our solar system in a gas disk. Our investigation extends the previous works by Masset and Snellgrove (2001) and Morbidelli and Crida (2007, MC07), which focussed on the dynamics of the Jupiter-Saturn system. The only systems that we found to reach a steady state are those in which the planets are locked in a quadruple mean motion resonance (i.e. each planet is in resonance with its neighbor). In total we found 6 such configurations. For the gas disk parameters found in MC07, these configurations are characterized by a negligible migration rate. After the disappearance of the gas, and in absence of planetesimals, only two of these six configurations (the least compact ones) are stable for a time of hundreds of millions of years or more. The others become unstable on a timescale of a few My. Our preliminary simulations show that, when a planetesimal disk is added beyond the orbit of the outermost planet, the planets can evolve from the most stable of these configurations to their current orbits in a fashion qualitatively similar to that described in Tsiganis et al. (2005). △ Less

Submitted 17 July, 2007; v1 submitted 12 June, 2007; originally announced June 2007.

Comments: The Astronomical Journal (17/07/2007) in press

Journal ref: Astron.J.134:1790-1798,2007

Chaos and the Effects of Planetary Migration on the Orbit of S/2000 S5 Kiviuq

Authors: V. Carruba, D. Nesvorny', J. A. Burns, M. 'Cuk, K. Tsiganis

Abstract: Among the many new irregular satellites that have been discovered in the last five years, at least six are in the so-called Kozai resonance. Due to solar perturbations, the argument of pericenter of a satellite usually precesses from 0 to 360 degrees. However, at inclinations higher than 39.3 degrees and lower than 140.7 degrees a new kind of behavior occurs for which the argument of pericenter… ▽ More Among the many new irregular satellites that have been discovered in the last five years, at least six are in the so-called Kozai resonance. Due to solar perturbations, the argument of pericenter of a satellite usually precesses from 0 to 360 degrees. However, at inclinations higher than 39.3 degrees and lower than 140.7 degrees a new kind of behavior occurs for which the argument of pericenter oscillates around +/-90 degrees. In this work we will concentrate on the orbital history of the saturnian satellite S/2000 S5 Kiviuq, one of the satellites currently known to be in such resonance Kiviuq's orbit is very close to the separatrix of the Kozai resonance. Due to perturbations from the other jovian planets, it is expected that orbits near the Kozai separatrix may show significant chaotic behavior. This is important because chaotic diffusion may transfer orbits from libration to circulation, and vice versa. To identify chaotic orbits we used two well-known methods: the Frequency Analysis Method (Laskar 1990) and Maximum Lyapunov Exponents (Benettin et al. 1980). Our results show that the Kozai resonance is crossed by a web of secondary resonances, whose arguments involve combinations of the argument of pericenter, the argument of the Great Inequality, longitude of the node, and other terms related to the secular frequencies g5, g6, and s6. Many test orbits whose precession period is close to the period of the Great Inequality (883 yrs), or some of its harmonics, are trapped by these secondary resonances, and show significant chaotic behavior. Planetary migration, by moving the locations of these secondary resonances, may have depleted an original population of Kozai resonators. △ Less

Submitted 28 April, 2004; originally announced April 2004.

Comments: 47 pages, 18 figures, submitted to the Astronomical Journal