Young Protoplanetary Disk Research Articles

Angular momentum transport via gravitational instability in the Elias 2–27 disc

Gravitational instability is thought to be one of the main drivers of angular momentum transport in young protoplanetary discs. The disc around Elias 2−27 offers a unique example of gravitational instability at work. It is young and massive, displaying two prominent spiral arms in dust continuum emission and global non-axisymmetric kinematic signatures in molecular line data. In this work, we used archival ALMA observations of 13CO line emission to measure the efficiency of angular momentum transport in the Elias 2−27 system through the kinematic signatures generated by gravitational instability, known as “GI wiggles”. Assuming the angular momentum is transported by the observed spiral structure and leveraging previously-derived dynamical disc mass measurements, the amount of angular momentum transport we found corresponds to an α-viscosity of α = 0.038 ± 0.018. This value implies an accretion rate onto the central star of log10 Ṁ⋆ = −6.99 ± 0.17 M⊙ yr−1, which reproduces the one observed value of log10 Ṁ⋆,obs = −7.2 ± 0.5 M⊙ yr−1 very well. The excellent agreement we have found serves as further proof that gravitational instability is the main driver of angular momentum transport acting in this system.

JWST Reveals Excess Cool Water near the Snow Line in Compact Disks, Consistent with Pebble Drift

Previous analyses of mid-infrared water spectra from young protoplanetary disks observed with the Spitzer-IRS found an anticorrelation between water luminosity and the millimeter dust disk radius observed with ALMA. This trend was suggested to be evidence for a fundamental process of inner disk water enrichment proposed decades ago to explain some properties of the solar system, in which icy pebbles drift inward from the outer disk and sublimate after crossing the snow line. Previous analyses of IRS water spectra, however, were uncertain due to the low spectral resolution that blended lines together. We present new JWST-MIRI spectra of four disks, two compact and two large with multiple radial gaps, selected to test the scenario that water vapor inside the snow line is regulated by pebble drift. The higher spectral resolving power of MIRI-MRS now yields water spectra that separate individual lines, tracing upper level energies from 900 to 10,000 K. These spectra clearly reveal excess emission in the low-energy lines in compact disks compared to large disks, demonstrating an enhanced cool component with T ≈ 170–400 K and equivalent emitting radius R eq ≈ 1–10 au. We interpret the cool water emission as ice sublimation and vapor diffusion near the snow line, suggesting that there is indeed a higher inward mass flux of icy pebbles in compact disks. Observation of this process opens up multiple exciting prospects to study planet formation chemistry in inner disks with JWST.

The rotational excitation of the water isotopologues by molecular hydrogen

ABSTRACT We present cross-sections and rate coefficients for rotational transitions in the water isotopologues D2O, H$_2^{18}$O, and HDO induced by collisions with para-H2(j2 = 0) and ortho-H2(j2 = 1). Quantum scattering calculations are performed at the full close-coupling level with the isotopic variants of an accurate full-dimensional H2O−H2 interaction potential. The D2O, H$_2^{18}$O, and HDO cross-sections are compared to the corresponding cross-sections for H2O. Large isotopic effects are observed in the case of D2O and HDO, in particular for collisions with p-H2(j2 = 0), while the 18O isotopic substitution is found to be negligible. Rate coefficients are provided for rotational transitions among all para-D2O, ortho-D2O, and HDO levels with internal energy below 300 cm−1 and for kinetic temperatures in the range 5–300 K. Non-LTE radiative transfer calculations show that the HDO 225.9 GHz and H$_2^{18}$O 203.4 GHz transitions recently detected with ALMA in the young proto-planetary disc V883 Ori should be inverted (weak masers) in a large fraction of the disc.

Dust Enrichment and Grain Growth in a Smooth Disk around the DG Tau Protostar Revealed by ALMA Triple Bands Frequency Observations

Characterizing the physical properties of dust grains in a protoplanetary disk is critical to comprehending the planet formation process. Our study presents Atacama Large Millimeter/submillimeter Array (ALMA) high-resolution observations of the young protoplanetary disk around DG Tau at a 1.3 mm dust continuum. The observations, with a spatial resolution of ≈0.″04, or ≈5 au, revealed a geometrically thin and smooth disk without substantial substructures, suggesting that the disk retains the initial conditions of the planet formation. To further analyze the distributions of dust surface density, temperature, and grain size, we conducted a multiband analysis with several dust models, incorporating ALMA archival data of the 0.87 and 3.1 mm dust polarization. The results showed that the Toomre Q parameter is ≲2 at a 20 au radius, assuming a dust-to-gas mass ratio of 0.01. This implies that a higher dust-to-gas mass ratio is necessary to stabilize the disk. The grain sizes depend on the dust models, and for the DSHARP compact dust, they were found to be smaller than ∼400 μm in the inner region (r ≲ 20 au) while exceeding larger than 3 mm in the outer part. Radiative transfer calculations show that the dust scale height is lower than at least one-third of the gas scale height. These distributions of dust enrichment, grain sizes, and weak turbulence strength may have significant implications for the formation of planetesimals through mechanisms such as streaming instability. We also discuss the CO snowline effect and collisional fragmentation in dust coagulation for the origin of the dust size distribution.

Characterizing fragmentation and sub-Jovian clump properties in magnetized young protoplanetary discs

ABSTRACT We study the initial development, structure, and evolution of protoplanetary clumps formed in three-dimensional resistive magnetohydrodynamic (MHD) simulations of self-gravitating discs. The magnetic field grows by means of the recently identified gravitational instability dynamo. Clumps are identified and their evolution is tracked finely both backward and forward in time. Their properties and evolutionary path is compared with clumps in companion simulations without magnetic fields. We find that magnetic and rotational energy are important in the clumps’ outer regions, while in the cores, despite appreciable magnetic field amplification, thermal pressure is most important in counteracting gravity. Turbulent kinetic energy is of a smaller scale than magnetic energy in the clumps. Compared with non-magnetized clumps, rotation is less prominent, which results in lower angular momentum in much better agreement with observations. In order to understand the very low sub-Jovian masses of clumps forming in MHD simulations, we revisit the perturbation theory of magnetized sheets finding support for a previously proposed magnetic destabilization in low-shear regions. This can help explaining why fragmentation ensues on a scale more than an order of magnitude smaller than that of the Toomre mass. The smaller fragmentation scale and the high magnetic pressure in clumps’ envelopes explain why clumps in magnetized discs are typically in the super-Earth to Neptune mass regime rather than super-Jupiter as in conventional disc instability. Our findings put forward a viable alternative to core accretion to explain widespread formation of intermediate-mass planets.

Formation of pebbles in (gravito-)viscous protoplanetary disks with various turbulent strengths

Aims. Dust plays a crucial role in the evolution of protoplanetary disks. We study the dynamics and growth of initially submicron dust particles in self-gravitating young protoplanetary disks with various strengths of turbulent viscosity. We aim to understand the physical conditions that determine the formation and spatial distribution of pebbles when both disk self-gravity and turbulent viscosity are at work. Methods. We performed thin-disk hydrodynamics simulations of self-gravitating protoplanetary disks over an initial time period of 0.5 Myr using the FEOSAD code. Turbulent viscosity was parameterized in terms of the spatially and temporally constant α parameter, while the effects of gravitational instability on dust growth were accounted for by calculating the effective parameter αGI. We considered the evolution of the dust component, including the momentum exchange with gas, dust self-gravity, and a simplified model of dust growth. Results. We find that the level of turbulent viscosity strongly affects the spatial distribution and total mass of pebbles in the disk. The α = 10−2 model is viscosity-dominated, pebbles are completely absent, and the dust-to-gas mass ratio deviates from the reference 1:100 value by no more than 30% throughout the extent of the disk. On the contrary, the α = 10−3 model and, especially, the α = 10−4 model are dominated by gravitational instability. The effective parameter α + αGI is now a strongly varying function of radial distance. As a consequence, a bottleneck effect develops in the innermost disk regions, which makes gas and dust accumulate in a ring-like structure. Pebbles are abundant in these models, although their total mass and spatial extent is sensitive to the dust fragmentation velocity and to the strength of gravitoturbulence. The use of the standard dust-to-gas mass conversion is not suitable for estimating the mass of pebbles. Conclusions. Our numerical experiments demonstrate that pebbles can already be abundant in protoplanetary disks at the initial stages of disk evolution. Dust growth models that consider disk self-gravity and ice mantles may be important for studying planet formation via pebble accretion.

The role of the drag force in the gravitational stability of dusty planet forming disc – I. Analytical theory

ABSTRACT Recent observations show that planet formation is already underway in young systems, when the protostar is still embedded into the molecular cloud and the accretion disc is massive. In such environments, the role of self-gravity (SG) and gravitational instability (GI) is crucial in determining the dynamical evolution of the disc. In this work, we study the dynamical role of drag force in self-gravitating discs as a way to form planetesimals in early protoplanetary stages. We obtain the dispersion relation for density-wave perturbations on a fluid composed of two phases (gas and dust) interacting through the common gravitation field and the mutual drag force, and we find that the stability threshold is determined by three parameters: the local dust-to-gas density ratio, the dust relative temperature, and the relevant Stokes number. In a region of parameters space, where young protoplanetary discs are likely to be found, the instability can be dust driven, occurring at small wavelengths. In this regime, the Jeans mass is much smaller than the one predicted by the standard GI model. This mechanism can be a viable way to form planetary cores in protostellar discs, since their predicted mass is about ∼10 M⊕.

Direct Formation of Planetary Embryos in Self-gravitating Disks

Giant planets have been discovered at large separations from the central star. Moreover, a striking number of young circumstellar disks have gas and/or dust gaps at large orbital separations, potentially driven by embedded planetary objects. To form massive planets at large orbital separations through core accretion within the disk lifetime, however, an early solid body to seed pebble and gas accretion is desirable. Young protoplanetary disks are likely self-gravitating, and these gravitoturbulent disks may efficiently concentrate solid material at the midplane driven by spiral waves. We run 3D local hydrodynamical simulations of gravitoturbulent disks with Lagrangian dust particles to determine whether particle and gas self-gravity can lead to the formation of dense solid bodies, seeding later planet formation. When self-gravity between dust particles is included, solids of size St = 0.1–1 concentrate within the gravitoturbulent spiral features and collapse under their own self-gravity into dense clumps up to several M ⊕ in mass at wide orbits. Simulations with dust that drift most efficiently, St = 1, form the most massive clouds of particles, while simulations with smaller dust particles, St = 0.1, have clumps with masses an order of magnitude lower. When the effect of dust backreaction onto the gas is included, dust clumps become smaller by a factor of a few but more numerous. The existence of large solid bodies at an early stage of the disk can accelerate the planet formation process, particularly at wide orbital separations, and potentially explain planets distant from the central stars and young protoplanetary disks with substructures.

Explosive instability of dust settling in a protoplanetary disc

ABSTRACT It is shown that gas-dust perturbations in a disc with dust settling to the disc mid-plane exhibit the non-linear three-wave resonant interactions between streaming dust wave (SDW) and two inertial waves (IW). In the particular case considered in this paper, SDW at the wavenumber k• = 2κ/(g$z$ts), where κ, g$z$, and ts are, respectively, epicyclic frequency, vertical gravitational acceleration, and particle’s stopping time, interacts with two IW at the lower wavenumbers k′ and k″ such that k′ < kDSI < k″ < k•, where kDSI = κ/(g$z$ts) is the wavenumber of the linear resonance between SDW and IW associated with the previously discovered linear dust settling instability. The problem is solved analytically in the limit of the small dust fraction. As soon as the dynamical dust back reaction on gas is taken into account, k•, k′, and k″ become slightly non-collinear and the emerging interaction of waves leads to simultaneous explosive growth of their amplitudes. This growth is explained by the conservative exchange with energy between the waves. The amplitudes of all three waves grow because the negative energy SDW transfers its energy to the positive energy IW. The product of the dimension-less amplitude of initially dominant wave and the time of explosion can be less than Keplerian time in a disc. It is shown that, generally, the three-wave resonance of an explosive type exists in a wide range of wavenumbers 0 < k• ≤ 2κ/(g$z$ts). An explosive instability of gas-dust mixture may facilitate the dust clumping and the subsequent formation of planetesimals in young protoplanetary discs.

Warping Away Gravitational Instabilities in Protoplanetary Discs

Abstract We perform three-dimensional smoothed-particle hydrodynamics simulations of warped, non-coplanar gravitationally unstable discs to show that as the warp propagates through the self-gravitating disk, it heats up the disk rendering it gravitationally stable, thus losing their spiral structure and appearing completely axisymmetric. In their youth, protoplanetary discs are expected to be massive and self-gravitating, which results in nonaxisymmetric spiral structures. However recent observations of young protoplanetary discs with the Atacama Large Millimeter/submillimeter Array have revealed that discs with large-scale spiral structure are rarely observed in the midplane. Instead, axisymmetric discs, with some also having ring and gap structures, are more commonly observed. Our work invloving warps, non-coplanar disk structures that are expected to commonly occur in young discs, potentially resolves this discrepancy between observations and theoretical predictions. We demonstrate that they are able to suppress the large-scale spiral structure of self-gravitating protoplanetary discs.

Carbon monoxide gas produced by a giant impact in the inner region of a young system.

Models of terrestrial planet formation predict that the final stages of planetary assembly-lasting tens of millions of years beyond the dispersal of young protoplanetary disks-are dominated by planetary collisions. It is through these giant impacts that planets like the young Earth grow to their final mass and achieve long-term stable orbital configurations1. A key prediction is that these impacts produce debris. So far, the most compelling observational evidence for post-impact debris comes from the planetary system around the nearby 23-million-year-old A-type star HD 172555. This system shows large amounts of fine dust with an unusually steep size distribution and atypical dust composition, previously attributed to either a hypervelocity impact2,3 or a massive asteroid belt4. Here we report the spectrally resolved detection of a carbon monoxide gas ring co-orbiting with dusty debris around HD 172555 between about six and nine astronomical units-a region analogous to the outer terrestrial planet region of our Solar System. Taken together, the dust and carbon monoxide detections favour a giant impact between large, volatile-rich bodies. This suggests that planetary-scale collisions, analogous to the Moon-forming impact, can release large amounts of gas as well as debris, and that this gas is observable, providing a window into the composition of young planets.

Streaming instability with multiple dust species – II. Turbulence and dust–gas dynamics at non-linear saturation

ABSTRACT The streaming instability is a fundamental process that can drive dust–gas dynamics and ultimately planetesimal formation in protoplanetary discs. As a linear instability, it has been shown that its growth with a distribution of dust sizes can be classified into two distinct regimes, fast- and slow-growth, depending on the dust-size distribution and the total dust-to-gas density ratio ϵ. Using numerical simulations of an unstratified disc, we bring three cases in different regimes into non-linear saturation. We find that the saturation states of the two fast-growth cases are similar to its single-species counterparts. The one with maximum dimensionless stopping time τs,max = 0.1 and ϵ = 2 drives turbulent vertical dust–gas vortices, while the other with τs,max = 2 and ϵ = 0.2 leads to radial traffic jams and filamentary structures of dust particles. The dust density distribution for the former is flat in low densities, while the one for the latter has a low-end cut-off. By contrast, the one slow-growth case results in a virtually quiescent state. Moreover, we find that in the fast-growth regime, significant dust segregation by size occurs, with large particles moving towards dense regions while small particles remain in the diffuse regions, and the mean radial drift of each dust species is appreciably altered from the (initial) drag-force equilibrium. The former effect may skew the spectral index derived from multiwavelength observations and change the initial size distribution of a pebble cloud for planetesimal formation. The latter along with turbulent diffusion may influence the radial transport and mixing of solid materials in young protoplanetary discs.

“Ashfall” Induced by Molecular Outflow in Protostar Evolution

Abstract Dust growth and its associated dynamics play key roles in the first phase of planet formation in young stellar objects. Observations have detected signs of dust growth in very young protoplanetary disks. Furthermore, signs of planet formation, gaps in the disk at a distance of several tens of au from the central protostar, are also reported. From a theoretical point of view, however, planet formation in the outer regions is difficult due to the rapid inward drift of dust, called the radial drift barrier. Here, on the basis of three-dimensional magnetohydrodynamical simulations of disk evolution with dust growth, we propose a mechanism called the “ashfall” phenomenon, induced by a powerful molecular outflow driven by a magnetic field that may circumvent the radial drift barrier. We find that the large dust that grows to a size of about a centimeter in the inner region of a disk is entrained by an outflow from the disk. Then, large dust decoupled from gas is ejected from the outflow due to centrifugal force, enriching the grown dust in the envelope and eventually falls onto the outer edge of the disk. The overall process is similar to the behavior of ashfall from volcanic eruptions. In the ashfall phenomenon, the Stokes number of dust increases by reaccreting to the less dense disk outer edge. This may allow the dust grains to overcome the radial drift barrier. Consequently, the ashfall phenomenon can provide a crucial assist for making the formation of the planetesimals in outer regions of the disk possible, and hence the formation of wide-orbit planets and gaps.

Gravito-turbulence and dynamo in poorly ionized protostellar discs – I. Zero-net-flux case

ABSTRACT In their early stages, protoplanetary discs are sufficiently massive to undergo gravitational instability (GI). This instability is thought to be involved in mass accretion, planet formation via gas fragmentation, the generation of spiral density waves, and outbursts. A key and very recent area of research is the interaction between the GI and magnetic fields in young protoplanetary discs, in particular whether this instability is able to sustain a magnetic field via a dynamo. We conduct 3D, stratified shearing-box simulations using two independent codes, PLUTO and Athena++, to characterize the GI dynamo in poorly ionized protostellar discs subject to ambipolar diffusion. We find that the dynamo operates across a large range of ambipolar Elssaser number Am (which characterizes the strength of ambipolar diffusion) and is particularly strong in the regime Am = 10–100, with typical magnetic to thermal energy ratios of order unity. The dynamo is only weakly dependent on resolution (at least for Am ≲ 100), box size, and cooling law. The magnetic field is produced by the combination of differential rotation and large-scale vertical roll motions associated with spiral density waves. Our results have direct implications for the dynamo process in young protoplanetary discs and possibly some regions of active galactic nucleus discs.

High-resolution Mid-infrared Spectroscopy of GV Tau N: Surface Accretion and Detection of NH3 in a Young Protoplanetary Disk

Abstract Physical processes that redistribute or remove angular momentum from protoplanetary disks can drive mass accretion onto the star and affect the outcome of planet formation. Despite ubiquitous evidence that protoplanetary disks are engaged in accretion, the process(es) responsible remain unclear. Here we present evidence for redshifted molecular absorption in the spectrum of a Class I source that indicates rapid inflow at the disk surface. High-resolution mid-infrared spectroscopy of GV Tau N reveals a rich absorption spectrum of individual lines of C2H2, HCN, NH3, and H2O. From the properties of the molecular absorption, we can infer that it carries a significant accretion rate ∼ 10−8–10−7 M ⊙ yr−1, comparable to the stellar accretion rates of active T Tauri stars. Thus, we may be observing disk accretion in action. The results may provide observational evidence for supersonic “surface accretion flows,” which have been found in MHD simulations of magnetized disks. The observed spectra also represent the first detection of NH3 in the planet formation region of a protoplanetary disk. With NH3 only comparable in abundance to HCN, it cannot be a major missing reservoir of nitrogen. If, as expected, the dominant nitrogen reservoir in inner disks is instead N2, its high volatility would make it difficult to incorporate into forming planets, which may help to explain the low nitrogen content of the bulk Earth.

Protoplanetary Disk Rings as Sites for Planetesimal Formation

Abstract Axisymmetric dust rings are a ubiquitous feature of young protoplanetary disks. These rings are likely caused by pressure bumps in the gas profile; a small bump can induce a traffic-jam-like pattern in the dust density, while a large bump may halt radial dust drift entirely. The resulting increase in dust concentration may trigger planetesimal formation by the streaming instability (SI), as the SI itself requires some initial concentration of dust. Here we present the first 3D simulations of planetesimal formation in the presence of a pressure bump modeled specifically after those seen by Atacama Large Millimeter/submillimeter Array. We place a pressure bump at the center of a large 3D shearing box, along with an initial solid-to-gas ratio of Z = 0.01, and we include both particle back-reaction and particle self-gravity. We consider millimeter-sized and centimeter-sized particles separately. For simulations with centimeter-sized particles, we find that even a small pressure bump leads to the formation of planetesimals via the SI; a pressure bump does not need to fully halt radial particle drift for the SI to become efficient. Furthermore, pure gravitational collapse via concentration in pressure bumps (such as would occur at sufficiently high concentrations and without the SI) is not responsible for planetesimal formation. For millimeter-sized particles, we find tentative evidence that planetesimal formation does not occur. If this result is confirmed at higher resolution, it could put strong constraints on where planetesimals can form. Ultimately, our results show that for centimeter-sized particles planetesimal formation in pressure bumps is extremely robust.

Thermal evolution of protoplanetary disks: from β-cooling to decoupled gas and dust temperatures

Aims. We explore the long-term evolution of young protoplanetary disks with different approaches to computing the thermal structure determined by various cooling and heating processes in the disk and its surroundings. Methods. Numerical hydrodynamics simulations in the thin-disk limit were complemented with three thermal evolution schemes: a simplified β-cooling approach with and without irradiation, where the rate of disk cooling is proportional to the local dynamical time; a fiducial model with equal dust and gas temperatures calculated taking viscous heating, irradiation, and radiative cooling into account; and a more sophisticated approach allowing decoupled dust and gas temperatures. Results. We found that the gas temperature may significantly exceed that of dust in the outer regions of young disks thanks to additional compressional heating caused by the infalling envelope material in the early stages of disk evolution and slow collisional exchange of energy between gas and dust in low-density disk regions. However, the outer envelope shows an inverse trend, with the gas temperatures dropping below that of dust. The global disk evolution is only weakly sensitive to temperature decoupling. Nevertheless, separate dust and gas temperatures may affect the chemical composition, dust evolution, and disk mass estimates. Constant-β models without stellar and background irradiation fail to reproduce the disk evolution with more sophisticated thermal schemes because of the intrinsically variable nature of the β-parameter. Constant-β models with irradiation more closely match the dynamical and thermal evolution, but the agreement is still incomplete. Conclusions. Models allowing separate dust and gas temperatures are needed when emphasis is placed on the chemical or dust evolution in protoplanetary disks, particularly in subsolar metallicity environments.

Polarized emission by aligned grains in the Mie regime: Application to protoplanetary disks observed by ALMA

Context. The azimuthal polarization patterns observed in some protoplanetary disks by the Atacama Large Millimetre Array (ALMA) at millimeter wavelengths have raised doubts about whether they are truly produced by dust grains that are aligned with the magnetic field lines. These conclusions were based on the calculations of dust polarized emission in the Rayleigh regime, that is, for grain sizes that are much smaller than the wavelength. However, the grain size in such disks is typically estimated to be in the range of 0.1−1 mm from independent observations. Aims. We study the dust polarization properties of aligned grains in emission in the Mie regime, that is, when the mean grain size approaches the wavelength. Methods. By using the T-MATRIX and DustEM codes, we computed the spectral dependence of the polarization fraction in emission for grains in perfect spinning alignment for various grain size distributions. We restricted our study to weakly-elongated oblate and prolate grains of astrosilicate composition that have a mean size ranging from 10 μm to 1 mm. Results. In the submillimeter and millimeter wavelength range, the polarization by B-field aligned grains becomes negative for grains larger than ∼250 μm, meaning that the polarization vector becomes parallel to the B-field. The transition from the positive to the negative polarization occurs at a wavelength of λ ∼ 1 mm. The regime of negative polarization does not exist for grains that are smaller than ∼100 μm. Conclusions. When using realistic grain size distributions for disks with grains up to the submillimeter sizes, the polarization direction of thermal emission by aligned grains is shown to be parallel to the direction of the magnetic field over a significant fraction of the wavelengths typically used to observe young protoplanetary disks. This property may explain the peculiar azimuthal orientation of the polarization vectors in some of the disks observed with ALMA and attest to the conserved ability of dust polarized emission to trace the magnetic field in disks.

ALMA observations require slower Core Accretion runaway growth

ABSTRACT Due to recent high-resolution ALMA observations, there is an accumulating evidence for presence of giant planets with masses from ${\sim } 0.01 \, {\rm {M}}_{\rm {J}}$ to a few $\, {\rm {M}}_{\rm {J}}$ with separations up to 100 au in the annular structures observed in young protoplanetary discs. We point out that these observations set unique ‘live’ constraints on the process of gas accretion on to sub-Jovian planets that were not previously available. Accordingly, we use a population synthesis approach in a new way: we build time-resolved models and compare the properties of the synthetic planets with the ALMA data at the same age. Applying the widely used gas accretion formulae leads to a deficit of sub-Jovian planets and an overabundance of a few Jupiter mass planets compared to observations. We find that gas accretion rate on to planets needs to be suppressed by about an order of magnitude to match the observed planet mass function. This slower gas giant growth predicts that the planet mass should correlate positively with the age of the protoplanetary disc, albeit with a large scatter. This effect is not clearly present in the ALMA data but may be confirmed in the near future with more observations.

Gravitoturbulent dynamos in astrophysical discs

The origin of large-scale and coherent magnetic fields in astrophysical discs is an important and long standing problem. Researchers commonly appeal to a turbulent dynamo, sustained by the magneto-rotational instability (MRI), to supply the large-scale field. But research over the last decade in particular has demonstrated that various non-ideal MHD effects can impede or extinguish the MRI, especially in protoplanetary disks. In this paper we propose a new scenario, by which the magnetic field is generated and sustained via the gravitational instability (GI). We use 3D stratified shearing box simulations to characterise the dynamo and find that it works at low magnetic Reynolds number (from unity to ~100) for a wide range of cooling times and boundary conditions. The process is kinematic, with a relatively fast growth rate ($< 0.1 \Omega$), and it shares some properties of mean field dynamos. The magnetic field is generated via the combination of differential rotation and spiral density waves, the latter providing compressible horizontal motions and large-scale vertical rolls. At greater magnetic Reynolds numbers the build up of large-scale field is diminished and instead small-scale structures emerge from the breakdown of twisted flux ropes. We propose that GI may be key to the dynamo engine not only in young protoplanetary discs but also in some AGN and galaxies.