Context. In the early evolution of a planetary system, a pair of planets may be captured in a mean motion resonance while still embedded in their nesting circumstellar disk. Aims. The goal is to ...estimate the direction and amount of shift in the semimajor axis of the resonance location due to the disk gravity as a function of the gas density and mass of the planets. The stability of the resonance lock when the disk dissipates is also tested. Methods. The orbital evolution of a large number of systems is numerically integrated within a three-body problem in which the disk potential is computed as a series of expansion. This is a good approximation, at least over a limited amount of time. Results. Two different resonances are studied: the 2:1 and the 3:2. In both cases the shift is inwards, even if by a different amount, when the planets are massive and carve a gap in the disk. For super-Earths, the shift is instead outwards. Different disk densities, Σ, are considered and the resonance shift depends almost linearly on Σ. The gas dissipation leads to destabilization of a significant number of resonant systems, in particular if it is fast. Conclusions. The presence of a massive circumstellar disk may significantly affect the resonant behavior of a pair of planets by shifting the resonant location and by decreasing the size of the stability region. The disk dissipation may explain some systems found close to a resonance but not locked in it.
The dynamics of systems of two and three planets, initially placed on circular and nearly coplanar orbits, is explored in the proximity of their stability limit. The evolution of a large number of ...systems is numerically computed and their dynamical behaviour is investigated with the frequency map analysis as chaos indicator. Following the guidance of this analysis, it is found that for two-planet systems the dependence of the Hill limit on the planet mass, usually made explicit through Hill's radius parametrization, does not appear to be fully adequate. In addition, frequent cases of stable chaos are found in the proximity of the Hill limit. For three-planet systems, the usual approach adopted in numerical explorations of their stability, where the planets are initially separated by multiples of the mutual Hill radius, appears too reducing. A detailed sampling of the parameter space reveals that systems with more packed inner planets are stable well within previous estimates of the stability limit. This suggests that a two-dimensional approach is needed to outline when three-planet systems are dynamically stable.
Planet–planet scattering is a major dynamical mechanism able to significantly alter the architecture of a planetary system. In addition to that, it may also affect the formation and retention of a ...debris disc by the system. A violent chaotic evolution of the planets can easily clear leftover planetesimal belts preventing the ignition of a substantial collisional cascade that can give origin to a debris disc. On the other end, a mild evolution with limited steps in eccentricity and semimajor axis can trigger the formation of a debris disc by stirring an initially quiet planetesimal belt. The variety of possible effects that planet–planet scattering can have on the formation of debris discs is analysed and the statistical probability of the different outcomes is evaluated. This leads to the prediction that systems which underwent an episode of chaotic evolution might have a lower probability of harbouring a debris disc.
ABSTRACT We present two-dimensional hydrodynamic simulations using the Smoothed Particle Hydrodynamic code, VINE, to model a self-gravitating binary system. We model configurations in which a ...circumbinary torus+disk surrounds a pair of stars in orbit around each other and a circumstellar disk surrounds each star, similar to that observed for the GG Tau A system. We assume that the disks cool as blackbodies, using rates determined independently at each location in the disk by the time dependent temperature of the photosphere there. We assume heating due to hydrodynamical processes and to radiation from the two stars, using rates approximated from a measure of the radiation intercepted by the disk at its photosphere. We simulate a suite of systems configured with semimajor axes of either a = 62 AU ("wide") or a = 32 AU ("close"), and with assumed orbital eccentricity of either e = 0 or e = 0.3. Each simulation follows the evolution for ∼6500-7500 yr, corresponding to about three orbits of the torus around the center of mass. Our simulations show that strong, sharply defined spiral structures are generated from the stirring action of the binary and that, in some cases, these structures fragment into 1-2 massive clumps. The torus quickly fragments into several dozen such fragments in configurations in which either the binary is replaced by a single star of equal mass, or radiative heating is neglected. The spiral structures extend inwards to the circumstellar environment as large scale material streams for which most material is found on trajectories that return it to the torus on a timescale of 1-200 yr, with only a small fraction accreting into the circumstellar environment. The spiral structures also propagate outwards through the torus, generating net outwards mass flow, and eventually losing coherence at large distances from the stars. The torus becomes significantly eccentric in shape over most of its evolution. In all configurations, accretion onto the stars occurs at a steady rate of a few ×10−8 M yr−1, with the net result that, without replenishment, the disk lifetimes would be shorter than ∼104 yr. Our simulations show that only wide orbit configurations are able to retain circumstellar disks, by virtue of accretion driven from the robust material streams generated in wide configurations, which are very weak in close configurations. In wide, eccentric orbit configurations, accretion is episodic and occurs preferentially onto the secondary, with rates strongly peaked near the binary periapse. Based on our results, we conclude that the GG Tau A torus is strongly self gravitating and that a major contribution to its thermal energy input is the shock dissipation associated with spiral structures generated both by self gravitating disturbances and by the stirring action of the binary. We interpret the sharply defined features observed in the torus as manifestations of such spiral structures. We interpret the low density disk surrounding it as an excretion disk created by the outward mass flux generated by the spiral arms as they propagate outwards. Typical eccentricities calculated for the shape of the tori modeled in our simulations are large enough to account for the supposed ∼20° mutual inclination between the stellar orbit plane of GG Tau A and its surrounding torus through a degeneracy between the interpretation of inclination of the torus and its eccentricity. We therefore interpret the observations in favor of a coplanar system with an eccentric torus. Because accretion onto the disks occurs at rates sufficient to sustain them only in wide orbit configurations, we conclude that the gas currently resident in the circumstellar disks of the GG Tau A system has been accreted from the torus within the past few thousand years. Although circumstellar disks will persist over time spans long enough to permit planet formation, the overall environment remains unfavorable due to high temperatures and other conditions. Given the presence of circumstellar disks, robust accretion streams, and our interpretation of the GG Tau A stellar orbit plane as coplanar with the torus surrounding it, we conclude that the GG Tau A system is in an eccentric, a ∼ 62 AU orbit, resolving questions in the literature regarding its orbit parameters.
The amount of dust present in circumstellar disks is expected to steadily decrease with age due to the growth from m-sized particles to planetesimals and planets. Mature circumstellar disks, however, ...can be observed to contain significant amounts of dust and possess high dust-to-gas ratios. Using HD 163296 as our case study, we explore how the formation of giant planets in disks can create the conditions for collisionally rejuvenating the dust population, halting or reversing the expected trend. We combine N-body simulations with statistical methods and impact scaling laws to estimate the dynamical and collisional excitation of the planetesimals due to the formation of HD 163296's giant planets. We show that this process creates a violent collisional environment across the disk that can inject collisionally produced second-generation dust into it, significantly contributing to the observed dust-to-gas ratio. The spatial distribution of the dust production can explain the observed local enrichments in HD 163296's inner regions. The results obtained for HD 163296 can be extended to any disk with embedded forming giant planets and may indicate a common evolutionary stage in the life of such circumstellar disks. Furthermore, the dynamical excitation of the planetesimals could result in the release of transient, nonequilibrium gas species like H2O, CO2, NH3, and CO in the disk due to ice sublimation during impacts and, due to the excited planetesimals being supersonic with respect to the gas, could produce bow shocks in the latter that could heat it and cause a broadening of its emission lines.
ABSTRACT
Pairs of planets in a system may end up close to their host star on eccentric orbits as a consequence of planet–planet scattering, Kozai, or secular migration. In this scenario, general ...relativity and secular perturbations have comparable time-scales and may interfere with each other with relevant effects on the eccentricity and pericenter evolution of the two planets. We explore, both analytically and via numerical integration, how the secular evolution is changed by general relativity for a wide range of different initial conditions. We find that when the faster secular frequency approaches the general relativity precession rate, which typically occurs when the outer planet moves away from the inner one, it relaxes to it and a significant damping of the proper eccentricity of the inner planet occurs. The proper eccentricity of the outer planet is reduced as well due to the changes in the secular interaction of the bodies. The lowering of the peak eccentricities of the two planets during their secular evolution has important implications on their stability. A significant number of two-planet systems, otherwise chaotic because of the mutual secular perturbations, are found stable when general relativity is included.
Context. We explore the dynamical evolution of a planet that is embedded in a circumstellar disk, as part of a binary system where the orbital plane of the companion star is significantly tilted with ...respect to the initial disk plane. Aims. Our aim is to test whether the planet remains within the disk and continues to migrate towards the star in a Type I/II mode, in spite of the secular perturbations of the companion star. Our findings, could explain why observed exoplanets have significant inclination in relation to the equatorial plane of their host star. Methods. We used two different smoothed particle hydrodynamic codes, VINE and PHANTOM, to model the evolution of a star+disk+planet system with a companion star over time. Results. After an initial coupled evolution, the inclinations of the disk and the planet begin to differ significantly. The period of oscillation of the disk inclination, in relation to the initial plane, becomes shorter for the planet, which evolves independently after about 104 yr, following a perturbed N-body behaviour. However, the planet continues to migrate towards the star because, during its orbital motion, it crosses the disk plane, and friction with the gas causes angular momentum loss. Conclusions. In a significantly inclined binary system, disks and planets are not dynamically coupled for small binary separations but evolve almost independently. The planet abandons the disk, and because of the onset of a significant mutual inclination, it interacts with the gas only when its orbit intersects the disk plane. The drift of the planet towards the star is not due to Type I/II, where the planet is embedded in the disk, but to the friction with the gas while crossing the disk.
Context. Planet–planet (P–P) scattering is an efficient and robust dynamical mechanism for producing eccentric exoplanets. Coupled to tidal interactions with the central star, this phenomenon can ...also explain close-in giant planets on circularized and potentially misaligned orbits. Aims. We explore scattering events occurring close to the star and test if they can reproduce the main features of the observed orbital distribution of giant exoplanets on tight orbits. Methods. In our modeling we exploited a numerical integration code based on the Hermite algorithm and including the effects of general relativity, dynamical tides, and two-body collisions. Results. We find that P–P scattering events occurring in systems with three giant planets initially moving on circular orbits close to their star produce a population of planets similar to that presently observed, including eccentric and misaligned close-in planets. The contribution of tides and general relativity is relevant in determining the final outcome of the chaotic phase. Conclusions. Even if two-body collisions dominate the chaotic evolution of three planets in crossing orbits close to their star, the final distribution shows a significant number of planets on eccentric orbits. The highly misaligned close-in giant planets are instead produced by systems where the initial semimajor axis of the inner planet was around 0.2 au or beyond.
Context. Meteoroids impacting terrestrial planets at high speed may have different effects. On bodies without atmospheres, such as the Moon and Mercury, they form impact craters and contribute to the ...gardening process through which the surface material is constantly mixed. The interaction of high-speed meteoroids with the atmosphere of Venus, the Earth, and Mars, may lead to the deposition in the ionosphere of species such as neutral Mg or Fe and their ionized atoms, caused by ablation processes during the entry. Aims. In this work we estimate and compare the flux and impact speeds onto the planets of the inner solar system by numerically integrating the orbital evolution of putative dust particles of asteroidal and cometary origin. Methods. The trajectories of dust particles of different sizes are computed with a numerical code that accounts for the gravitational forces due to all planets, the Poynting-Robertson drag and the solar wind drag. The flux of dust grains on each planet is estimated by calibrating the outcome of our model with the flux on the Earth reported previously. Results. We obtain new estimates of the flux and impact velocities for both asteroidal and cometary dust particles on Venus and Mars. For Venus we find that cometary grains enter the planet atmosphere at higher speeds, possibly contributing to the upper layers, while asteroidal grains would be relevant for the lower layers, possibly leading to a compositional gradient. This effect is also present for Mars, but it is less marked. We also find that analytical predictions, not taking radiative forces into account, of both flux and average impact speed are reliable for Mars but fail for Venus because of the complex dynamical evolution of grains in the inner solar system. Conclusions. Our results on the velocity distributions and fluxes of micrometeoroids on the terrestrial planets can be used to put stringent contraints on models that estimate either the superficial material mixing that is due to meteoroid impacts or the formation of ionospheric layers for planets with an atmosphere.
I explore the dynamics of small dust particles in transport-dominated circumstellar debris discs in binary star systems. In these tenuous discs the effects of mutual collisions are negligible and ...their morphology is determined by Poynting-Robertson drag and, possibly, by the strong perturbations induced by the interaction with the interstellar medium (ISM) flux of neutral atoms. The force due to the ISM flux can significantly affect the dynamical behaviour of the dust grains, causing a fast inward drift and a large periodic oscillation of both eccentricity and inclination. If the disc is around a star in a binary system, the gravity of the companion star competes with the ISM force and the dynamics is complex. The balance between the two forces depends strongly on the binary semimajor axis a
B and eccentricity e
B. In a scenario with an ISM flux similar to that observed in the Solar system neighbourhood, the binary secular perturbations, assuming a mass ratio of 0.5, dominate over the ISM force when a
B < 600 au and e
B= 0.2. This occurs when the dust disc is generated by a parent body ring encompassed between 50 and 60 aufrom the primary star. For a larger binary eccentricity e
B= 0.6, the limit moves to a
B < 700 au.
Within these values of a
B, the time-scale of the binary secular perturbations is much shorter than the period of the ISM-induced orbital variations, and the disc shape and density distribution are dominated by the companion gravity. It appears slightly eccentric and, if the binary is coplanar with the disc, only a limited warping due to the ISM perturbations is observed. In this scenario, the strong ISM perturbations, which may significantly affect debris discs around single stars embedded in strong ISM winds, are almost completely silenced.
For larger semimajor axes, the scenario is reversed with the ISM perturbations ruling the dynamics of the dust. The disc develops a large clump oriented at 90° with respect to the direction of the ISM flux and it appears highly warped in space. Since the period of the two perturbations can be analytically estimated, it is possible to deduce by easy computations when a given system is governed either by the binary gravity or by the ISM perturbations.