首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 78 毫秒
1.
On the migration of a system of protoplanets   总被引：1，自引：0，他引：1  
The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter-sized planets is studied numerically. The disc is assumed to be flat and non-self-gravitating; this is modelled by the planar (two-dimensional) Navier–Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass M Jup each, and the radial distances from the star are one and two semimajor axes of Jupiter, respectively.
During the evolution a joint wide annular gap is created by the planets. Both planets increase their mass owing to accretion of gas from the disc: after about 2500 orbital periods of the inner planet it has reached a mass of 2.3  M Jup, while the outer planet has reached a mass of 3.2  M Jup. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc. The semimajor axis of the inner planet remains constant as there is very little gas left in its vicinity to induce any migration. When the distance of close approach eventually becomes smaller than the mutual Hill radius, the eccentricities increase strongly and the system may become unstable.
If disc depletion occurs rapidly enough before the planets come too close to each other, a stable system similar to our own Solar system may remain. Otherwise the orbits may become unstable and produce systems like υ And.  相似文献   

2.
We present the results of hydrodynamic simulations of Jovian mass protoplanets that form in circumbinary discs. The simulations follow the orbital evolution of the binary plus protoplanet system acting under their mutual gravitational forces, and forces exerted by the viscous circumbinary disc. The evolution involves the clearing of the inner circumbinary disc initially, so that the binary plus protoplanet system orbits within a low density cavity. Continued interaction between disc and protoplanet causes inward migration of the planet towards the inner binary. Subsequent evolution can take three distinct paths: (i) the protoplanet enters the 4 : 1 mean motion resonance with the binary, but is gravitationally scattered through a close encounter with the secondary star; (ii) the protoplanet enters the 4 : 1 mean motion resonance, the resonance breaks, and the planet remains in a stable orbit just outside the resonance; (iii) when the binary has initial eccentricity   e bin≥ 0.2  , the disc becomes eccentric, leading to a stalling of the planet migration, and the formation of a stable circumbinary planet.
These results have implications for a number of issues in the study of extrasolar planets. The ejection of protoplanets in close binary systems provides a source of 'free-floating planets', which have been discovered recently. The formation of a large, tidally truncated cavity may provide an observational signature of circumbinary planets during formation. The existence of protoplanets orbiting stably just outside a mean motion resonance (4 : 1) in the simulations indicate that such sites may harbour planets in binary star systems, and these could potentially be observed. Finally, the formation of stable circumbinary planets in eccentric binary systems indicates that circumbinary planets may not be uncommon.  相似文献   

3.
There is evidence for the existence of massive planets at orbital radii of several hundred au from their parent stars where the time-scale for planet formation by core accretion is longer than the disc lifetime. These planets could have formed close to their star and then migrated outwards. We consider how the transfer of angular momentum by viscous disc interactions from a massive inner planet could cause significant outward migration of a smaller outer planet. We find that it is in principle possible for planets to migrate to large radii. We note, however, a number of effects which may render the process somewhat problematic.  相似文献   

4.
The migration and growth of protoplanets in protostellar discs   总被引：1，自引：0，他引：1  
We investigate the gravitational interaction of a Jovian-mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions are investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Non-linear two-dimensional hydrodynamic simulations are employed using three independent numerical codes.
A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near-circular orbit after a characteristic viscous time‐scale of ∼104 initial orbital periods. This is found to be independent of whether the protoplanet is allowed to accrete mass or not. Inward migration is helped by the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion on to the central star. Maximally accreting protoplanets reach about 4 Jovian masses on reaching the neighbourhood of the central star. Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is ∼5 Jupiter masses, if they originate in protoplanetary discs similar to the minimum-mass solar nebula. This is because of the reduced accretion rates obtained for planets of increasing mass.
Assuming that some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets as well as their inward orbital migration can be accounted for by consideration of disc–protoplanet interactions during the late stages of giant planet formation.  相似文献   

5.
We investigate the migration of massive extrasolar planets caused by gravitational interaction with a viscous protoplanetary disc. We show that a model in which planets form at 5 au at a constant rate, before migrating, leads to a predicted distribution of planets that is a steeply rising function of log( a ), where a is the orbital radius. Between 1 and 3 au, the expected number of planets per logarithmic interval in a roughly doubles. We demonstrate that, once selection effects are accounted for, this is consistent with current data, and then extrapolate the observed planet fraction to masses and radii that are inaccessible to current observations. In total, approximately 15 per cent of stars targeted by existing radial velocity searches are predicted to possess planets with masses  0.3< M p sin( i )<10 M J  and radii  0.1< a <5 au  . A third of these planets (around 5 per cent of the target stars) lie at the radii most amenable to detection via microlensing. A further  5–10  per cent of stars could have planets at radii of  5< a <8 au  that have migrated outwards. We discuss the probability of forming a system (akin to the Solar system) in which significant radial migration of the most massive planet does not occur. Approximately  10–15  per cent of systems with a surviving massive planet are estimated to fall into this class. Finally, we note that a smaller fraction of low-mass planets than high-mass planets is expected to survive without being consumed by the star. The initial mass function for planets is thus predicted to rise more steeply towards small masses than the observed mass function.  相似文献   

6.
Planets orbiting a planetesimal circumstellar disc can migrate inward from their initial positions because of dynamical friction between planets and planetesimals. The migration rate depends on the disc mass and on its time evolution. Planets that are embedded in long-lived planetesimal discs, having total mass of 10−4– 0.01 M⊙ , can migrate inward a large distance and can survive only if the inner disc is truncated or as a result of tidal interaction with the star. In this case the semimajor axis, a , of the planetary orbit is less than 0.1 au. Orbits with larger a are obtained for smaller values of the disc mass or for a rapid evolution (depletion) of the disc. This model may explain not only several of the orbital features of the giant planets that have been discovered in recent years orbiting nearby stars, but also the metallicity enhancement found in several stars associated with short-period planets.  相似文献   

7.
Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time-scale that is short compared with the orbital period. We use a combination of hydrodynamic simulations and N -body orbit integrations to study the long-term evolution of a fragmenting disc with an initial mass ratio to the star of   M disc/ M *= 0.1  . For a disc that is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar-mass star) up to  ≈0.01 M⊙  . Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 – in which the companions are close to or beyond the deuterium burning limit – appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ∼1 au should differ from that of stars with lower mass planetary companions.  相似文献   

8.
We investigate the properties of circumplanetary discs formed in three-dimensional, self-gravitating radiation hydrodynamical models of gas accretion by protoplanets. We determine disc sizes, scaleheights, and density and temperature profiles for different protoplanet masses, in solar nebulae of differing grain opacities.
We find that the analytical prediction of circumplanetary disc radii in an evacuated gap  ( R Hill/3)  from Quillen & Trilling yields a good estimate for discs formed by high-mass protoplanets. The radial density profiles of the circumplanetary discs may be described by power laws between   r −2  and   r −3/2  . We find no evidence for the ring-like density enhancements that have been found in some previous models of circumplanetary discs. Temperature profiles follow a  ∼ r −7/10  power law regardless of protoplanet mass or nebula grain opacity. The discs invariably have large scaleheights  ( H / r > 0.2)  , making them thick in comparison with their encompassing circumstellar discs, and they show no flaring.  相似文献   

9.
We investigate the gravitational interaction between a planet and an optically thin protoplanetary disc, performing local three-dimensional hydrodynamical simulations. In the present study, we take account of radiative energy transfer in optically thin discs. Before the stage of planetary accretion, dust opacity is expected to decrease significantly because of grain growth and planetesimal formation. Thus, it would be reasonable to consider optically thin discs in the disc–planet interaction. Furthermore, we focus on small planets that can neither capture disc gas nor open a disc gap. The one-sided torque exerted on a planet by an optically thin disc is examined for various values of the disc optical thickness (<1). In optically thin discs, the temperature behind the density waves is lower than the unperturbed value because of radiative cooling. Heating due to shock dissipation is less effective than radiative cooling. Because of radiative cooling, the density distribution around the planet is not axisymmetric, which exerts an additional torque on the planet. The torque enhancement becomes maximum when the cooling time is comparable with the Keplerian period. The enhancement is significant for low-mass planets. For planets with  3 M⊕  , the additional one-sided torque can be 40 per cent of the torque in the isothermal case. The radiative cooling is expected to change the differential torque and the migration speed of planets, too.  相似文献   

10.
We have performed N -body numerical simulations of the exchange of angular momentum between a massive planet and a 3D Keplerian disc of planetesimals. Our interest is directed at the study of the classical analytical expressions of the lineal theory of density waves, as representative of the dynamical friction in discs 'dominated by the planet' and the orbital migration of the planets with regard to this effect. By means of a numerical integration of the equations of motion, we have carried out a set of numerical experiments with a large number of particles  ( N ≥10 000)  , and planets with the mass of Jupiter, Saturn and one core mass of the giant planets in the Solar system  ( M c=10 M⊕)  . The torque, measured in a phase in which a 'steady forcing' is clearly measurable, yields inward migration in a minimum-mass solar disc  (Σ∼10 g cm-2  ), with a characteristic drift time of ∼ a few 106 yr. The planets predate the disc, but the orbital decay rate is not sufficient to allow accretion in a time-scale relevant to the formation of giant planets. We found reductions of the measured torque on the planet, with respect to the linear theory, by a factor of 0.38 for M c, 0.04 for Saturn and 0.01 for Jupiter, due to the increase in the perturbation on the disc. The behaviour of planets whose mass is larger than M c is similar to the one of type II migrators in gaseous discs. Our results suggest that, in a minimum mass, solar planetesimals disc, type I migrations occur for masses smaller than M c, whereas for this mass value it could be a transition zone between the two types of migration.  相似文献   

11.
We consider the minimum mass planet, as a function of radius, that is capable of opening a gap in an α-accretion disc. We estimate that a half-Jupiter mass planet can open a gap in a disc with accretion rate     for viscosity parameter  α= 0.01  , and solar mass and luminosity. The minimum mass is approximately proportional to     . This estimate can be used to rule out the presence of massive planets in gapless accretion discs. We identify two radii at which an inwardly migrating planet may become able to open a gap and so slow its migration; the radius at which the heating from viscous dissipation is similar to that from stellar radiation in a flared disc, and the radius at which the disc becomes optically thin in a self-shadowed disc. In the inner portions of the disc, we find that the minimum planet mass required to open a gap is only weakly dependent on radius. If a migrating planet is unable to open a gap by the time it reaches either of the transition radii, then it is likely to be lost on to the star. If a gap-opening planet cuts off disc accretion allowing the formation of a central hole or clearing in the disc then we would estimate that the clearing radius would approximately be proportional to the stellar mass.  相似文献   

12.
Recent observations point to the presence of structured dust grains in the discs surrounding young brown dwarfs, thus implying that the first stages of planet formation take place also in the substellar regime. Here, we investigate the potential for planet formation around brown dwarfs and very low-mass stars according to the sequential core accretion model of planet formation. We find that, for a brown dwarf mass 0.05 M⊙, our models predict a maximum planetary mass of  ∼5   M⊕  , orbiting with semimajor axis ∼ 1 au. However, we note that the predictions for the mass–semimajor axis distribution are strongly dependent upon the models chosen for the disc surface density profiles and the assumed distribution of disc masses. In particular, if brown dwarf disc masses are of the order of a few Jupiter masses, Earth-mass planets might be relatively frequent, while if typical disc masses are only a fraction of Jupiter mass, we predict that planet formation would be extremely rare in the substellar regime. As the observational constraints on disc profiles, mass dependencies and their distributions are poor in the brown dwarf regime, we advise caution in validating theoretical models only on stars similar to the Sun and emphasize the need for observational data on planetary systems around a wide range of stellar masses. We also find that, unlike the situation around solar-like stars, Type II migration is totally absent from the planet formation process around brown dwarfs, suggesting that any future observations of planets around brown dwarfs would provide a direct measure of the role of other types of migration.  相似文献   

13.
Using 2D magnetohydrodynamic (MHD) numerical simulations performed with two different finite-difference Eulerian codes, we analyse the effect that a toroidal magnetic field has on low-mass planet migration in non-turbulent protoplanetary discs. The presence of the magnetic field modifies the waves that can propagate in the disc. In agreement with a recent linear analysis, we find that two magnetic resonances develop on both sides of the planet orbit, which contribute to a significant global torque. In order to measure the torque exerted by the disc on the planet, we perform simulations in which the latter is either fixed on a circular orbit or allowed to migrate. For a     planet, when the ratio β between the square of the sound speed and that of the Alfven speed at the location of the planet is equal to 2, we find inward migration when the magnetic field   B φ  is uniform in the disc, reduced migration when   B φ  decreases as   r −1  and outward migration when   B φ  decreases as   r −2  . These results are in agreement with predictions from the linear analysis. Taken as a whole, our results confirm that even a subthermal stable field can stop inward migration of an earth-like planet.  相似文献   

14.
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.  相似文献   

15.
We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼104-yr bursts with ̇ ∼10−5 M⊙ yr−1, separated by quiescent intervals lasting ∼105 yr where ̇ ≈10−8 M⊙ yr−1. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.  相似文献   

16.
We use numerical simulations to model the migration of massive planets at small radii and compare the results with the known properties of 'hot Jupiters' (extrasolar planets with semimajor axes   a < 0.1  au). For planet masses   M pl sin  i > 0.5 M J  , the evidence for any 'pile-up' at small radii is weak (statistically insignificant), and although the mass function of hot Jupiters is deficient in high-mass planets as compared to a reference sample located further out, the small sample size precludes definitive conclusions. We suggest that these properties are consistent with disc migration followed by entry into a magnetospheric cavity close to the star. Entry into the cavity results in a slowing of migration, accompanied by a growth in orbital eccentricity. For planet masses in excess of 1 Jupiter mass we find eccentricity growth time-scales of a few ×105 yr, suggesting that these planets may often be rapidly destroyed. Eccentricity growth appears to be faster for more massive planets which may explain changes in the planetary mass function at small radii and may also predict a pile-up of lower mass planets, the sample of which is still incomplete.  相似文献   

17.
Detectable debris discs are thought to require dynamical excitation ('stirring'), so that planetesimal collisions release large quantities of dust. We investigate the effects of the secular perturbations of a planet, which may lie at a significant distance from the planetesimal disc, to see if these perturbations can stir the disc, and if so over what time-scale. The secular perturbations cause orbits at different semimajor axes to precess at different rates, and after some time   t cross  initially non-intersecting orbits begin to cross. We show that   t cross∝ a 9/2disc/( m pl e pl a 3pl)  , where   m pl, e pl  and   a pl  are the mass, eccentricity and semimajor axis of the planet, and   a disc  is the semimajor axis of the disc. This time-scale can be faster than that for the growth of planetesimals to Pluto's size within the outer disc. We also calculate the magnitude of the relative velocities induced among planetesimals and infer that a planet's perturbations can typically cause destructive collisions out to 100 s of au. Recently formed planets can thus have a significant impact on planet formation in the outer disc which may be curtailed by the formation of giant planets much closer to the star. The presence of an observed debris disc does not require the presence of Pluto-sized objects within it, since it can also have been stirred by a planet not in the disc. For the star ε Eridani, we find that the known radial velocity planet can excite the planetesimal belt at 60 au sufficiently to cause destructive collisions of bodies up to 100 km in size, on a time-scale of 40 Myr.  相似文献   

18.
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 .  相似文献   

19.
We study the viscous evolution of protoplanetary discs driven by the combined action of magnetohydrodynamic turbulence, resulting from the magneto-rotational instability (MRI), self-gravity torques, parametrized in terms of an effective viscosity and an additional viscous agent of unspecified origin. The distribution of torques driving the evolution of the disc is calculated by analysing where in the disc the MRI develops and, to incorporate the effect of self-gravity, calculating the Toomre parameter. We find that, generally, discs rapidly evolve towards a configuration where the intermediate regions, from a fraction of an au to a few au, are stable against the MRI due to their low-ionization degree. As an additional source of viscosity is assumed to operate in those regions, subsequent evolution of the disc is eruptive. Brief episodes of high mass accretion ensue as the criterion for the development of the MRI is met in the low-ionization region. The radial distribution of mass and temperature in the disc differs considerably from disc models with constant α parameter or layered accretion models, with potentially important consequences on the process of planet formation.  相似文献   

20.
We investigate the fast (type III) migration regime of high-mass protoplanets orbiting in protoplanetary discs. This type of migration is dominated by corotational torques. We study the details of flow structure in the planet's vicinity, the dependence of migration rate on the adopted disc model and the numerical convergence of models (independence of certain numerical parameters such as gravitational softening).
We use two-dimensional hydrodynamical simulations with adaptive mesh refinement, based on the flash code with improved time-stepping scheme. We perform global disc simulations with sufficient resolution close to the planet, which is allowed to freely move throughout the grid. We employ a new type of equation of state in which the gas temperature depends on both the distance to the star and planet, and a simplified correction for self-gravity of the circumplanetary gas.
We find that the migration rate in the type III migration regime depends strongly on the gas dynamics inside the Hill sphere (Roche lobe of the planet) which, in turn, is sensitive to the aspect ratio of the circumplanetary disc. Furthermore, corrections due to the gas self-gravity are necessary to reduce numerical artefacts that act against rapid planet migration. Reliable numerical studies of type III migration thus require consideration of both the thermal and the self-gravity corrections, as well as a sufficient spatial resolution and the calculation of disc–planet attraction both inside and outside the Hill sphere. With this proviso, we find type III migration to be a robust mode of migration, astrophysically promising because of a speed much faster than in the previously studied modes of migration.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号