The total number of giant planets in debris discs with central clearings

Infrared spectra from the Spitzer Space Telescope ( SSC ) of many debris discs are well fit with a single blackbody temperature which suggest clearings within the disc. We assume that clearings are caused by orbital instability in multiple planet systems with similar configurations to our own. These planets remove dust-generating planetesimal belts as well as dust generated by the outer disc that is scattered or drifts into the clearing. From numerical integrations, we estimate a minimum planet spacing required for orbital instability (and so planetesimal and dust removal) as a function of system age and planet mass. We estimate that a 10⁸ yr old debris disc with a dust disc edge at a radius of 50 au hosted by an A star must contain approximately five Neptune mass planets between the clearing radius and the iceline in order to remove all primordial objects within it. We infer that known debris disc systems contain at least a fifth of a Jupiter mass in massive planets. The number of planets and spacing required is insensitive to the assumed planet mass. However, an order of magnitude higher total mass in planets could reside in these systems if the planets are more massive.     

Weak lensing signal in unified dark matter models

This paper is an extension of the work done by Pierens & Nelson in which they investigated the behaviour of a two-planet system embedded in a protoplanetary disc. They put a Jupiter mass gas giant on the internal orbit and a lower mass planet on the external one. We consider here a similar problem taking into account a gas giant with mass in the range 0.5 to  1 M _J  and a Super-Earth (i.e. a planet with mass  ≤10 M _⊕  ) as the outermost planet. By changing disc parameters and planet masses, we have succeeded in getting the convergent migration of the planets which allows for the possibility of their resonant locking. However, in the case in which the gas giant has the mass of Jupiter, before any mean-motion first-order commensurability could be achieved, the Super-Earth is caught in a trap when it is very close to the edge of the gap opened by the giant planet. This confirms the result obtained by Pierens & Nelson in their simulations. Additionally, we have found that, in a very thin disc, an apsidal resonance is observed in the system if the Super-Earth is captured in the trap. Moreover, the eccentricity of the small planet remains low, while the eccentricity of the gas giant increases slightly due to the imbalance between Lindblad and corotational resonances. We have also extended the work of Pierens & Nelson by studying analogous systems in which the gas giant is allowed to take sub-Jupiter masses. In this case, after conducting an extensive survey over all possible parameters, we have succeeded in getting the 1:2 mean-motion resonant configuration only in a disc with low aspect ratio and low surface density. However, the resonance is maintained just for a few thousand orbits. Thus, we conclude that for typical protoplanetary discs the mean-motion commensurabilities are rare if the Super-Earth is located on the external orbit relative to the gas giant.     

Exoplanet's Figure and Its Interior

Along with the development of the observing technology, the observation and study on the exoplanets’ oblateness and apsidal precession have achieved significant progress. The oblateness of an exoplanet is determined by its interior density profile and rotation period. Between its Love number k₂ and core size exists obviously a negative correlation. So oblateness and k₂ can well constrain its interior structure. Starting from the Lane-Emden equation, the planet models based on different polytropic indices are built. Then the flattening factors are obtained by solving the Wavre's integro-differential equation. The result shows that the smaller the polytropic index, the faster the rotation, and the larger the oblateness. We have selected 469 exoplanets, which have simultaneously the observed or estimated values of radius, mass, and orbit period from the NASA (National Aeronautics and Space Administration) Exoplanet Archive, and calculated their flattening factors under the two assumptions: tidal locking and fixed rotation period of 10.55 hours. The result shows that the flattening factors are too small to be detected under the tidal locking assumption, and that 28% of exoplanets have the flattening factors larger than 0.1 under the fixed rotation period of 10.55 hours. The Love numbers under the different polytropic models are solved by the Zharkov's approach, and the relation between k₂ and core size is discussed.     

Migration of giant planets in a time-dependent planetesimal accretion disc

In this paper we develop further the model for the migration of planets introduced in Del Popolo et al. We first model the protoplanetary nebula as a time-dependent accretion disc, and find self-similar solutions to the equations of the accretion disc that give us explicit formulae for the spatial structure and the temporal evolution of the nebula. These equations are then used to obtain the migration rate of the planet in the planetesimal disc, and to study how the migration rate depends on the disc mass, on its time evolution and on some values of the dimensionless viscosity parameter α . We find that planets that are embedded in planetesimal discs, having total mass of  10^-4-0.1 M_⊙  , can migrate inward a large distance for low values of α (e.g.,   α ≃10^-3-10^-2)  and/or large disc mass, and can survive only if the inner disc is truncated or because of tidal interaction with the star. Orbits with larger a are obtained for smaller values of the disc mass and/or for larger values of α . This model may explain several orbital features of the recently discovered giant planets orbiting nearby stars.     

Stable satellites around extrasolar giant planets

In this work, we study the stability of hypothetical satellites of extrasolar planets. Through numerical simulations of the restricted elliptic three-body problem we found the borders of the stable regions around the secondary body. From the empirical results, we derived analytical expressions of the critical semimajor axis beyond which the satellites would not remain stable. The expressions are given as a function of the eccentricities of the planet, e _P, and of the satellite, e _sat. In the case of prograde satellites, the critical semimajor axis, in the units of Hill's radius, is given by a _E≈ 0.4895   (1.0000 − 1.0305 e _P− 0.2738 e _sat). In the case of retrograde satellites, it is given by a _E≈ 0.9309  (1.0000 − 1.0764 e _P− 0.9812 e _sat). We also computed the satellite stability region ( a _E) for a set of extrasolar planets. The results indicate that extrasolar planets in the habitable zone could harbour the Earth-like satellites.     

Detailed survey of the phase space around Nix and Hydra

We present a detailed survey of the dynamical structure of the phase space around the new moons of the Pluto–Charon system. The spatial elliptic restricted three-body problem was used as model and stability maps were created by chaos indicators. The orbital elements of the moons are in the stable domain on the semimajor axis, eccentricity and inclination spaces. The structures related to the 4:1 and 6:1 mean motion resonances are clearly visible on the maps. They do not contain the positions of the moons, confirming previous studies. We showed the possibility that Nix might be in the 4:1 resonance if its argument of pericentre or longitude of node falls in a certain range. The results strongly suggest that Hydra is not in the 6:1 resonance for arbitrary values of the argument of pericentre or longitude of node.     

The host stars of extrasolar planets have normal lithium abundances

The lithium abundances of planet-harbouring stars have been compared with the lithium abundances of open clusters and field stars. Young (chromospherically active) and subgiant stars have been eliminated from the comparison because they are at different stages of evolution and Li processing than the planet-harbouring stars, and hence have systematically higher Li abundances. The analysis showed that the Li abundances of the planet-harbouring stars are indistinguishable from those of non-planet-harbouring stars of the same age, temperature and composition. This conclusion is opposite to that arrived at by Gonzalez & Laws; it is believed that the field-star sample used by them contained too wide a range of ages, evolutionary types and temperatures to be accommodated by the model that they adopted to describe the dependence of Li on the parameters. The Li abundance does not appear set to provide key insights into the formation and evolution of planetary systems.     

Numerical simulations of type III planetary migration – III. Outward migration of massive planets

We present a numerical study of rapid, so-called type III migration for Jupiter-sized planets embedded in a protoplanetary disc. We limit ourselves to the case of outward migration, and study in detail its evolution and physics, concentrating on the structure of the corotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. We find that the outward-directed type III migration can be started if the initial conditions support the initial average non-dimensional migration rate bigger than one. Unlike the inward-directed migration, in the outward migration the migration rate increases due to the growing of the volume of the co-orbital region. We find the migration to be strongly dependent on the rate of the mass accumulation in the circumplanetary disc, leading to two possible regimes of migration, fast and slow. The structure of the co-orbital region and the stopping mechanism differs between these two regimes.     

Numerical simulations of type III planetary migration – II. Inward migration of massive planets

We present a numerical study of rapid, so-called type III migration for Jupiter-sized planets embedded in a protoplanetary disc. We limit ourselves to the case of inward migration, and study in detail its evolution and physics, concentrating on the structure of the corotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using the results of global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. The initial conditions are chosen to satisfy the condition for rapid inward migration. We find that type III migration can be divided into two regimes, fast and slow. The structure of the co-orbital region, mass accumulation rate and migration behaviour differ between these two regimes. All our simulations show a transition from the fast to the slow regime, ending type III migration well before reaching the star. The stopping radius is found to be larger for more massive planets and less massive discs. A sharp density drop is also found to be an efficient stopping mechanism. In the fast migration regime the migration rate and induced eccentricity are lower for less massive discs, but almost do not depend on planet mass. Eccentricity is damped on the migration time-scale.     

Jupiter and Super-Earth embedded in a gaseous disc

In this paper we investigate the evolution of a pair of interacting planets – a Jupiter-mass planet and a Super-Earth with a mass of  5.5 M _⊕  – orbiting a Solar-type star and embedded in a gaseous protoplanetary disc. We focus on the effects of type I and II orbital migrations, caused by the planet–disc interaction, leading to the capture of the Super-Earth in first-order mean-motion resonances by the Jupiter. The stability of the resulting resonant system in which the Super-Earth is on the internal orbit relative to the Jupiter is studied numerically by means of full 2D hydrodynamical simulations. Our main aim is to determine the Super-Earth behaviour in the presence of the gas giant in the system. It is found that the Jupiter captures the Super-Earth into the interior 3:2 or 4:3 mean-motion resonance, and that the stability of such configurations depends on the initial positions of the planets and on the evolution of the eccentricity. If the initial separation of the orbits of the planets is larger than or close to that required for the exact resonance, the final outcome is the migration of the pair of planets at a rate similar to that of the gas giant, at least for the time of our simulations. Otherwise, we observe a scattering of the Super-Earth from the disc. The evolution of planets immersed in a gaseous disc is compared with their behaviour in the case of the classical three-body problem when the disc is absent.     

Long-term stability of planetary orbits between Jupiter and Saturn

We extend our two previous studies on the existence of stable orbits in the Solar System by examining the domain between Jupiter and Saturn. We place (1) a massless object, (2) a Moon-mass object, (3) a Mars-mass object, (4) an Earth-mass object, and (5) a Uranus-mass object in the said region. Note that these objects are considered separately in the framework of our simulations. Our goal is to explore the orbital stability of those objects. We employ the Lie-integration method, which is fast and well established, allowing us to solve the respective differential equations for the $N$-body system. Hence, we consider the celestial bodies spanning from Jupiter to Neptune, including the aforementioned test object, the main focus for our model simulations. The integrations indicate that in some models the test objects placed in the region between Jupiter and Saturn reside in that region for more than 600 Myr. Between 5 and 10 au, mean-motion resonances (MMRs) take place acting upon the test objects akin to simulations of Paper I and II. Our models indicate relatively small differences for the long-term stability of the five test objects notwithstanding their vastly different masses. Generally, it is found that between $a_{\text{ini}}=7.04$ and 7.13 au the orbits become unstable mostly within 5 million years and further out, that is, up to $a_{\text{ini}}=7.29$ au, the duration of stability lengthens to up to hundreds of millions of years.     

On the width and shape of the corotation region for low-mass planets

We study the coorbital flow for embedded, low-mass planets. We provide a simple semi-analytic model for the corotation region, which is subsequently compared to high-resolution numerical simulations. The model is used to derive an expression for the half-width of the horseshoe region, x _s, which in the limit of zero softening is given by   x _s/ r _p= 1.68( q / h )^1/2  , where q is the planet to central star mass ratio, h is the disc aspect ratio and   r _p  is the orbital radius. This is in very good agreement with the same quantity measured from simulations. This result is used to show that horseshoe drag is about an order of magnitude larger than the linear corotation torque in the zero-softening limit. Thus, the horseshoe drag, the sign of which depends on the gradient of specific vorticity, is important for estimates of the total torque acting on the planet. We further show that phenomena, such as the Lindblad wakes, with a radial separation from corotation of approximately a pressure scaleheight H can affect x _s, even though for low-mass planets   x _s≪ H   . The effect is to distort streamlines and reduce x _s through the action of a back pressure. This effect is reduced for smaller gravitational softening parameters and planets of higher mass, for which x _s becomes comparable to H .     

Young Jupiters are faint: new models of the early evolution of giant planets

Here we show preliminary calculations of the cooling and contraction of a 2 M_J planet. These calculations, which are being extended to 1–10 M_J, differ from other published “cooling tracks” in that they include a core accretion‐gas capture formation scenario, the leading theory for the formation of gas giant planets.We find that the initial post‐accretionary intrinsic luminosity of the planet is ∼3 times less than previously published models which use arbitrary initial conditions. These differences last a few tens of millions of years. Young giant planets are intrinsically fainter than has been previously appreciated. We also discuss how uncertainties in atmospheric chemistry and the duration of the formation time of giant planets lead to challenges in deriving planetary physical properties from comparison with tabulated model values. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Low-energy Lunar Trajectories with Lunar Flybystwo

The low-energy lunar trajectories with lunar flybys are investigated based on the Sun-Earth-Moon bicircular problem (BCP). The characteristics of the distribution of trajectories in the phase space are summarized. Using the invariant manifolds in the BCP system, the low-energy lunar trajectories with lunar flybys are sought. Then, take time as an augmented dimension in the phase space of a nonautonomous system, we present the state space map and reveal the distribution of these lunar trajectories in the phase space. Consequently, we find that the low-energy lunar trajectories exist as families, and that the every moment in the Sun-Earth-Moon synodic period can be the departure date. Finally, we analyse the velocity increment, transfer duration, and system energy for the different trajectory families, and obtain the velocity-impulse optimal family and the transfer-duration optimal family, respectively.     

Dynamical determination of the mass of the Kuiper Belt from motions of the inner planets of the Solar system

In this paper we determine dynamically the mass of the Kuiper Belt Objects by exploiting the latest least-squares determinations of the extra rates of perihelia of the inner planets of the Solar system. By modelling classical Kuiper Belt Objects as an ecliptic ring of finite thickness, we obtain 0.033 ± 0.115 in units of terrestrial masses. For resonant Kuiper Belt Objects, a two-ring model yields 0.018 ± 0.063. These values are consistent with recent determinations obtained using ground- and space-based optical techniques. Some implications for precise tests of Einsteinian and post-Einsteinian gravity are briefly discussed.     

A new perspective on the irregular satellites of Saturn – II. Dynamical and physical origin

The origin of the irregular satellites of the giant planets has been long debated since their discovery. Their dynamical features argue against an in situ formation suggesting that they are captured bodies, yet there is no global consensus on the physical process at the basis of their capture. In this paper, we explore the collisional capture scenario, where the actual satellites originated from impacts occurred within Saturn's influence sphere. By modelling the inverse capture problem, we estimated the families of orbits of the possible parent bodies and the specific impulse needed for their capture. The orbits of these putative parent bodies are compared to those of the minor bodies of the outer Solar system to outline their possible region of formation. Finally, we tested the collisional capture hypothesis on Phoebe by taking advantage of the data supplied by Cassini on its major crater, Jason. Our results presented a realistic range of solutions matching the observational and dynamical data.