Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe interaction of planets with a disc with MHD turbulence IV: Migration rates of embedded protoplanets
121
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
(Abridged) We present global disc and local shearing box simulations of planets interacting with a MHD turbulent disc. We examine the torque exerted by the disc on the embedded planets as a function of planet mass, and thus make a first study of orbital migration of planets due to interaction with turbulent discs. Global simulations were performed for a disc with H/R=0.07 and planet masses M_p=3,10,30 Earth masses, and 3 Jupiter masses. Shearing box runs were performed for values of (M_p/M_*)/(H/R)^3=0.1,0.3,1.0 and 2.0, M_* being the central mass. These allow embedded and gap forming planets to be examined. In all cases the instantaneous torque exerted on a planet showed strong fluctuations. In the embedded cases it oscillated between negative and positive values, and migration occurs as a random walk, unlike the usual type I migration. Running time averages for embedded planets over 20-25 orbital periods show that large fluctuations occur on longer time scales, preventing convergence of the average torque to well defined values, or even to a well defined sign. Fluctuations become relatively smaller for larger masses, giving better convergence, due to the planets perturbation of the disc becoming larger than the turbulence in its vicinity. Eventually gap formation occurs, with a transition to type II migration. The existence of significant fluctuations occurring in turbulent discs on long time scales is important for lower mass embedded protoplanets. If significant fluctuations occur on the longest disc evolutionary time scales, convergence of torque running averages for practical purposes will not occur, and the migration behaviour of low mass protoplanets considered as an ensemble would be very different from predictions of type I theory for laminar discs.
rate research
Read More
We present a global MHD simulation of a turbulent accretion disc interacting with a protoplanet of 5 Jupiter masses. The disc model had H/r=0.1,and a value of the Shakura & Sunyaev alpha ~ 0.005. The protoplanet opened a gap in the disc, with the interaction leading to inward migration on the expected time scale. Spiral waves were launched by the protoplanet and although they were diffused and dissipated through interaction with the turbulence, they produced an outward angular momentum flow which compensated for a reduced flux associated with the turbulence, so maintaining the gap. When compared with laminar disc models with the same estimated alpha, the gap was found to be deeper and wider indicating that the turbulent disc behaved as if it possessed a smaller alpha. This may arise for two reasons. First, the turbulence does not provide a source of constantly acting friction in the near vicinity of the planet that leads to steady mass flow into the gap region. Instead the turbulence is characterised by large fluctuations in the radial velocity, and time averaging over significant time scales is required to recover the underlying mass flow through the disc. Near the planet the disc material experiences high amplitude perturbations on time scales that are short relative to the time scale required for averaging. The disc response is thus likely to be altered relative to a Navier--Stokes model. Second, the simulation indicates that an ordered magnetic connection between the inner and outer disc can occur enabling angular momentum to flow out across the gap, helping to maintain it independently of the protoplanets tide. This type of effect may assist gap formation for smaller mass protoplanets which otherwise would not be able to maintain them.
We present the results of both global cylindrical disc simulations and local shearing box simulations of protoplanets interacting with a disc undergoing MHD turbulence with zero net flux magnetic fields. We investigate the nature of the disc response and conditions for gap formation. This issue is an important one for determining the type and nature of the migration of the protoplanet, with the presence of a deep gap being believed to enable slower migration. For both types of simulation we find a common pattern of behaviour for which the main parameter determining the nature of the response is $M_p R^3/(M_* H^3)$, with $M_p$, $M_*$, $R$, and $H$ being the protoplanet mass, the central mass, the orbital radius and the disc semi-thickness respectively. We find that as this parameter is increased towards 0.1, the presence of the protoplanet is first indicated by the appearance of the well known trailing wake which, although it may appear erratic on account of the turbulence, appears to be well defined. Once the above parameter exceeds a number around unity a gap starts to develop inside which the magnetic energy density tends to be concentrated in the high density wakes. This gap formation condition can be understood from simple dimensional considerations of the conditions for nonlinearity, and the balance of angular momentum transport due to Maxwell and Reynolds stresses with that due to tidal torques. An important result is that the basic flow morphology in the vicinity of the protoplanet is very similar in both the local and global simulations. This indicates that local shearing box simulations, which are computationally less demanding, capture much of the physics of disc-planet interaction. Thus they may provide a useful tool for studying the local interaction between forming protoplanets and turbulent, protostellar discs.
We predict magnitudes for young planets embedded in transition discs, still affected by extinction due to material in the disc. We focus on Jupiter-size planets at a late stage of their formation, when the planet has carved a deep gap in the gas and dust distributions and the disc starts being transparent to the planet flux in the infrared (IR). Column densities are estimated by means of three-dimensional hydrodynamical models, performed for several planet masses. Expected magnitudes are obtained by using typical extinction properties of the disc material and evolutionary models of giant planets. For the simulated cases located at $5.2$ AU in a disc with local unperturbed surface density of $127$ $mathrm{g} cdot mathrm{cm}^{-2}$, a $1$ $M_J$ planet is highly extincted in J-, H- and K-bands, with predicted absolute magnitudes $ge 50$ mag. In L- and M-bands extinction decreases, with planet magnitudes between $25$ and $35$ mag. In the N-band, due to the silicate feature on the dust opacities, the expected magnitude increases to $40$ mag. For a $2$ $M_J$ planet, the magnitudes in J-, H- and K-bands are above $22$ mag, while for L-, M- and N-bands the planet magnitudes are between $15$ and $20$ mag. For the $5$ $M_J$ planet, extinction does not play a role in any IR band, due to its ability to open deep gaps. Contrast curves are derived for the transition discs in CQ Tau, PDS70, HL Tau, TW Hya and HD163296. Planet mass upper-limits are estimated for the known gaps in the last two systems.
We study how migration affects stars of a galaxy with a thin stellar disc and thicker stellar components. The simulated galaxy has a strong bar and lasting spiral arms. We find that the amplitude of the churning (change in angular momentum) is similar for thin and thick components, and of limited amplitude, and that stars of all components can be trapped at the corotation of the bar. At the exception of those stars trapped at the corotation, we find that stars that are far from their initial guiding radius are more likely so due to blurring rather than churning effects. We compare the simulation to orbits integration with a fixed gravitational potential rotating at a constant speed. In the latter case, stars trapped at corotation are churned periodically outside and inside the corotation radius, with a zero net average. However, as the bar speed of the simulated galaxy decreases and its corotation radius increases, stars trapped at corotation for several Gyrs can be churned outwards on average. We study the location of extreme migrators (stars experimenting the largest churning) and find that extreme migrators come from regions on the leading side of the effective potential local maxima.
As the fundamental physical process with many astrophysical implications, the diffusion of cosmic rays (CRs) is determined by their interaction with magnetohydrodynamic (MHD) turbulence. We consider the magnetic mirroring effect arising from MHD turbulence on the diffusion of CRs. Due to the intrinsic superdiffusion of turbulent magnetic fields, CRs with large pitch angles that undergo mirror reflection, i.e., bouncing CRs, are not trapped between magnetic mirrors, but move diffusively along the magnetic field, leading to a new type of parallel diffusion. This diffusion is in general slower than the diffusion of non-bouncing CRs with small pitch angles that undergo gyroresonant scattering. The critical pitch angle at the balance between magnetic mirroring and pitch-angle scattering is important for determining the diffusion coefficients of both bouncing and non-bouncing CRs and their scalings with the CR energy. We find non-universal energy scalings of diffusion coefficients, depending on the properties of MHD turbulence.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
