Issue |
A&A
Volume 686, June 2024
|
|
---|---|---|
Article Number | A9 | |
Number of page(s) | 16 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/202348546 | |
Published online | 24 May 2024 |
Rotation curves in protoplanetary disks with thermal stratification
Physical model and observational evidence in MAPS disks
1
Dipartimento di Fisica, Università degli Studi di Milano,
Via Celoria 16,
Milano,
20133,
Italy
2
Leiden Observatory, Leiden University,
PO Box 9513,
2300 RA
Leiden,
The Netherlands
e-mail: martire@strw.leidenuniv.nl
3
Laboratoire Lagrange, Université Côte d’Azur, CNRS,
Observatoire de la Côte d’Azur,
06304
Nice,
France
4
School of Physics and Astronomy, Monash University,
Clayton,
VIC 3800,
Australia
5
Univ. Grenoble Alpes, CNRS, IPAG,
38000
Grenoble,
France
6
Alma Mater Studiorum Università di Bologna, Dipartimento di Fisica e Astronomia (DIFA),
Via Gobetti 93/2,
40129
Bologna,
Italy
7
INAF – Osservatorio Astrofisico di Arcetri,
L.go E. Fermi 5,
50125
Firenze,
Italy
8
Institute of Astronomy, University of Cambridge,
Madingley Road,
Cambridge
CB3 0HA,
UK
e-mail: e-mail: cl2000@cam.ac.uk
Received:
10
November
2023
Accepted:
16
February
2024
Context. In recent years, the gas kinematics probed by molecular lines detected with ALMA has opened a new window into the of study protoplanetary disks. High spatial and spectral resolution observations have revealed the complexity of protoplanetary disk structure. Drawing accurate interpretations of these data allows us to better comprehend planet formation.
Aims. We investigate the impact of thermal stratification on the azimuthal velocity of protoplanetary disks. High-resolution gas observations reveal velocity differences between CO isotopologues, which cannot be adequately explained with vertically isothermal models. The aim of this work is to determine whether a stratified model can explain this discrepancy.
Methods. We analytically solved the hydrostatic equilibrium for a stratified disk and we derived the azimuthal velocity. We tested the model with SPH numerical simulations and then we used it to fit for the star mass, disk mass, and scale radius of the sources in the MAPS sample. In particular, we used 12CO and 13CO datacubes.
Results. When thermal stratification is taken into account, it is possible to reconcile most of the inconsistencies between rotation curves of different isotopologues. A more accurate description of the CO rotation curves offers a deeper understanding of the disk structure. The best-fit values of star mass, disk mass, and scale radius become more realistic and more in line with previous studies. In particular, the quality of the scale radius estimate significantly increases when adopting a stratified model. In light of our results, we computed the gas-to-dust ratio and the Toomre Q parameter. Within our hypothesis, for all the sources, the gas-to-dust ratio appears higher but still close to the standard value of 100 (within a factor of 2). The Toomre Q parameter suggests that the disks are gravitationally stable (Q > 1). However, the systems that show spirals presence are closer to the conditions of gravitational instability (Q ~ 5).
Key words: accretion, accretion disks / hydrodynamics / molecular data / protoplanetary disks
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1 Introduction
Our understanding of the physical properties of protoplanetary disks has improved in recent years thanks to the Atacama Large Millimeter Array (ALMA; ALMA Partnership et al. 2015). High spectral and spatial resolution gas observations enable us to probe density, temperature, and velocity fields of protostellar disks, gaining unique information about their structure (Law et al. 2021; Calahan et al. 2021; Teague et al. 2022; Miotello et al. 2023; Pinte et al. 2023; Lodato et al. 2023). More recently, the large program Molecules with ALMA at Planet-forming Scales (MAPS; Öberg et al. 2021) targeted five protoplanetary disks (MWC 480, IM Lup, GM Aur, HD 163296, and AS 209) in several molecular lines. For optically thick line emission, the gas temperature can be measured along the emission surface directly from the peak surface brightness of the channel maps (Law et al. 2021). Given the varying heights of these emitting layers surfaces, it is possible to infer the thermal structure in disks, proving the existence of a vertical thermal stratification in them (Dartois et al. 2003; Rosenfeld et al. 2013; Pinte et al. 2018), as expected from basic radiative transfer arguments (Chiang & Goldreich 1997; D’Alessio et al. 1998, 1999). Although disk models have usually been considered as vertically isothermal, the vertical gradient of temperature leads to considerable corrections in the calculation of density structure and azimuthal velocity, which results in several percent deviations from the Keplerian velocity (Rosenfeld et al. 2013). Accounting for such differences is important not only to infer stellar masses, but also to accurately constrain the disk pressure structure and disk mass. Such parameters are of great importance with respect to interpreting velocity deviations, which may serve as signposts for planets (Pinte et al. 2018; Rabago & Zhu 2021; Bollati et al. 2021; Izquierdo et al. 2021; Bae et al. 2021), dust trapping (Teague et al. 2018; Rosotti et al. 2020) or disk instabilities (Hall et al. 2020; Terry et al. 2022; Longarini et al. 2021; Barraza-Alfaro et al. 2021).
In this paper, we analytically derive the density and velocity field of protostellar disks with thermal stratification, generalizing the work of Takeuchi & Lin (2002). A similar analysis of MAPS data with vertically isothermal disks was performed by Lodato et al. (2023). We test the model against hydrodynami-cal simulations and we apply it to the whole MAPS sample for 12CO and 13CO data. In Sect. 2, we present the model, solving the vertical hydrostatic equilibrium and obtaining an expression for the azimuthal velocity. In Sect. 3, we present the numerical setup and a comparison between the model and simulations. In Sect. 4, we apply the model and we discuss our findings. Finally, in Sect. 5, we compute the gas-to-dust ratio and the Toomre Q parameter, and we draw our conclusions.
2 Model
2.1 Assumptions
In our analytical calculations, we did not make any assumption on the surface density Σ, considering it as arbitrary. However, to apply the model to observations, we were forced to choose a parameterization for the surface density and we assume that it is described by the self-similar solution from Lynden-Bell & Pringle (1974):
(1)
where Md and Rc are the disk mass and the scale radius respectively; R is the cylindrical radius and γ is a free parameter describing the steepness of the surface density. The disk density at the midplane ρmid is:
(2)
where Hmid = cs,mid/Ωk is the typical scale height of the disk at the midplane, is the sound speed at the disk midplane, kb is the Boltzmann constant, Tmid = Tmid,100(R/100 au)−q is the temperature at midplane, µ is the mean molecular weight (usually assumed to be 2.1), mp is the proton mass, and
is the Keplerian angular velocity (G is the gravitational constant and M★ is the stellar mass).
From the literature (Chiang & Goldreich 1997; Dullemond et al. 2020) and observational data (Rosenfeld et al. 2013; Pinte et al. 2018; Law et al. 2021), we know that protoplanetary disks are thermally stratified. We take this into account by defining a function, f, that describes the dependency of the temperature, T, on height, such that:
(3)
(4)
We underline that the isothermal case can be obtained considering f ≡ 1, thus T = Tmid(R). As for the density, we assume that:
(5)
where g describes how the density changes vertically. We note that in order to smoothly connect the functions above to their value at midplane it is necessary that f(z = 0) = 1 = g(z = 0). Assuming a barotropic fluid, the pressure, P, is given by:
(6)
While the profile of f is arbitrary, this does not hold for g, whose value is set by solving the vertical hydrostatic equilibrium.
2.2 Hydrostatic equilibrium and rotation curve
To compute the vertical density profile we assume a non-self-gravitating disk under the condition of hydrostatic equilibrium in the vertical direction:
(7)
where Φ★ = −GM★/r is the stellar potential ( is the spherical radius). Equation (7) can be written as (for further details, see Appendix A):
(8)
Solving for log(f g), we find:
(9)
and hence the density is given by:
(10)
Assuming the condition of centrifugal balance, the rotation curve is given by the radial component of Navier-Stokes equation:
(11)
The first term in Eq. ((11)) can be written as (for further details, see Appendix A)
(12)
where is the Keplerian velocity. Therefore, the rotation curve is
(14)
In the self-similar case, this becomes:
(15)
where γ' = γ + (3 + q)/2. Each term of Eq. (15) can be easily interpreted: [1 + (z/R)2]−3/2 is the star contribution at the height z, γ′ is the effect of the power law scaling of the pressure, (2 − γ)(R/Rc)2−γ is the effect of the exponential truncation and the logarithmic term is the effect of the vertical stratification. Since the latter is the derivative of a product, we do not know its sign a priori; thus, we also do not know if the rotation is accelerated or slowed down by thermal stratification (see Appendix A). In any case, in all our attempts this term never dominates over the variation of gravity with z. Thus, we found rotation to slow down with z and this effect is more pronounced as compared to the isothermal case when considering the parameters of the MAPS sample.
We underline that for the isothermal case (ƒ ≡ 1), this expression reduces to the one derived and analyzed by Lodato et al. (2023), while Eq. (9) is simplified as:
(16)
Therefore, the density in the isothermal case is given by:
(17)
If the disk is self-gravitating, we should add to the right-hand side of Eq. (11) the self-gravitating term (Bertin & Lodato 1999):
(18)
where K(k) and E(k) are complete elliptic integrals (Abramowitz & Stegun 1970) and k2 = 4Rr/[(R + r)2+ z2].
![]() |
Fig. 1 Comparison between the analytical models and the simulation after eight outer orbits and at z = 20 au. Left panel: density field. Right panel: Pressure gradient term of the rotation curve. The red dots represent the simulation data, while the blue and orange line show the thermally stratified model and the isothermal one, respectively. The model with thermal stratification matches very well the simulation. |
2.3 Temperature prescriptions
The two parameterizations of the vertical temperature more often used are given by Dartois et al. (2003) and Dullemond et al. (2020). In this work we will use the one by Dullemond et al. (2020), which is given by:
(19)
where the atmospheric temperature is parameterized as , Tε is considered as an approximation of the temperature at the midplane Tє ≃ Tmid. Zq(R) is defined as Zq(R) = ζ100(R/100au)β and α is a parameter that describes where the transition from midplane to atmospheric temperature occurs in the vertical direction. We note that in this case, f(R, z = 0) ≠ 1 and, thus, the temperature does not smoothly connect to its value at midplane. We discuss this in Appendix B, but we underline that Eq. (19) is a good approximation for the five disks within the MAPS large program in most of the radial extent of the disk.
Once the function f is defined, Eqs. (10) and (14) can be solved semi-analytically and entirely specify the rotation curve. We have implemented this calculation in DYSC1.
3 Comparison with numerical simulations
In this work, we performed numerical smoothed particle hydrodynamics (SPH) simulations of protostellar disks using the PHANTOM code (Price et al. 2018). This code is widely used in the astrophysical community to study gas and dust dynamics in accretion disks (Dipierro et al. 2015; Ragusa et al. 2017; Curone et al. 2022) and it has recently been employed in kine-matical studies (Pinte et al. 2018; Hall et al. 2020; Terry et al. 2022; Verrios et al. 2022). The aim of this simulation is to test the model before applying it to actual data.
To test the analytical model, we simulated a thermally stratified disk using the parameters of MWC 480 presented in Law et al. (2021). The simulation has been performed with N = 106 gas particles, initially distributed as a tapered power law density profile, smoothed at the inner radius, with γ = 1 and Rc = 150 au, between Rin = 10 au and Rout = 400 au. The mass of the star is 2.1 M⊙. For the temperature structure we used the Dullemond prescription given by Eq. (19), with ζ0 = 7au, a = 2.78,β = −0.05, Tmid,100 = 27 K, q = 0.23, Tatm,100 = 69 K, and qatm = 0.7. The αSS Shakura & Sunyaev (Shakura & Sunyaev 1973) viscosity coefficient has been set to 0.005. No self-gravity or dust have been included in the simulation.
We let the system evolve and reach hydrostatic equilibrium. We observed that after a couple of orbits the system reaches a relaxed state. We decided to analyze the output of the simulation after eight outer orbits (~45 kyr). In Fig. 1, we show a comparison between the density and the velocity of the simulations (red dots) at z = 20 au and both the isothermal and stratified model predictions. The red dots represent the azimuthal average of the respective quantity computed by averaging over all SPH particles within each of the 50 radial bins and the error bar is the corresponding standard deviation. Since we are plotting quantities at z = 20 au, we have excluded the inner points because at those radii (R < 100 au) the disk has a smaller hydrostatic height H, causing numerical issues in the resolution of our simulation. The stratified model perfectly describes the density and the velocity field of the simulation and is a significant improvement over the isothermal one. In particular, in the right panel of Fig. 1, we can see that the difference between the azimuthal velocity and the Keplerian velocity reaches the 10–12% and only the stratified model is able to reproduce it.
4 Applying the model
4.1 The curves
In this section, we describe how we applied the model to the entire sample of disks from the MAPS ALMA Large Program (Öberg et al. 2021). We performed our fits under the assumption of vertically isothermal or stratified disk to then compare the results. For the vertically isothermal model, the thermal structure is defined by the hydrostatic height of the disk at R = 100 au and the power law coefficient of the temperature profile, q. These parameters were taken from Zhang et al. (2021). As for the stratified model, Law et al. (2021) obtained the two-dimensional (2D) temperature structure of the MAPS disks, using the Dullemond et al. (2020) prescription given by Eq. (19). We note that the rotation curve traced by a specific molecule is defined by:
(21)
where z(R) is the height of the emitting layer of the considered molecule. For the emitting layer, we use:
(22)
where the best-fit parameters have been obtained by Izquierdo et al. (2023). All the parameters used are summarized in Table 1.
The rotation curves (Fig. 2) can be obtained through different moment maps, according to the disk emission. We underline that with respect to the rotation curve extraction, we are only interested in measuring velocities from the frontside. Since three of the sources have strong contributions from the disk backside, we used a double-Bell decomposition to distinguish between these two components as introduced (Izquierdo et al. 2022). In this work, we have used an improved algorithm that performs this decomposition based on velocity priors obtained from the DISCMINER models (Izquierdo et al., in prep.).
Velocity extraction method, orientation parameters, thermal parameters and emitting surfaces for 12CO and 13CO data of the MAPS disks.
4.2 Results
The 12CO and 13CO data were simultaneously fitted with both the isothermal and stratified model including the self-gravitating contribution. The results are shown in Figs. 3–7, and the best-fitting parameters reported in Table 2 and Appendix C.
To quantify the importance of thermal stratification, we computed the relative difference between the squares of 12CO and 13CO rotation curves, as shown in Fig. 8. According to the vertical isothermal model, this quantity is:
(23)
which solely depends on the different height of the tracer, since it is assumed that the temperature does not change vertically. As for the stratified model, the expression is more complex, since it involves the evaluation of the term given by Eq. (9) at different heights. In this case, we expect to observe larger differences between the velocity of the two isotopologues, since there is an additional shift caused by the different emission temperature. To determine the importance of vertical stratification, we quantified the maximum value of the velocity shift between 12CO and 13CO, which can be predicted in the isothermal case:
(24)
where we assumed that typically z/R < 0.5. Hence, if the quantity is higher than 5%, the system cannot be described by an isothermal model, while it is likely that vertical stratification plays a significant role. It is important to note that the Eq. (24) depends on the star mass through
. We normalized the squared differences of the velocities by the square of the velocity for 13CO since this quantity is independent of the stellar mass. Figure 8 shows this quantity for the studied systems.
In the following, we will present the results of each disk, along with a discussion of the importance of thermal stratification. To compare the results, we performed our fits for both the vertically isothermal and stratified case. In addition, we computed the dust mass from millimetric emission at 283 GHz, using (Hildebrand 1983):
(25)
where d is the distance, Fv is the flux density in Jy, ĸv = 2.3(v/230 GHz)0.4cm2g−1 is the dust opacity, and Bv is the blackbody spectrum. In our analysis, we assumed T = 20 K and v = 283 GHz, while the flux densities were extracted from MAPS data. We recall that this equation implies that dust emission is optically thin. The results are reported in Table 3.
![]() |
Fig. 2 12CO and 13CO rotation curves of the MAPS disks extracted with DISCMINER. |
4.3 MWC 480
MWC 480 is a ~7 Myr Herbig Ae star located in the Taurus-Aurigae star-forming region at a distance of d = 162 pc (Montesinos et al. 2009). The most recent value of the stellar mass has been derived dynamically by Izquierdo et al. (2023) to be M★ = 1.97 M⊙. Zhang et al. (2021) through 2D thermo-chemical models computed disk mass and scale radius of the MAPS disks. For MWC 480, these values are Md = 0.16 M⊙ and Rc = 200 au.
By inspecting the 12CO and 13CO rotation curves (Fig. 2), there is no evident sign of thermal stratification, since the two curves do not differ significantly. Figure 3 shows that the two models are nearly indistinguishable, but in Fig. 8 we see that the stratified model is better at reproducing the data. When we assume an isothermal model, we obtain M★ = 1.969 ± 0.002 M⊙, Md = 0.201 ± 0.002M⊙ and Rc = 80 ± 1 au, while for the stratified model, it is M★ = 2.027 ± 0.002 M⊙, Md = 0.150 ± 0.002 M⊙ and Rc = 128 ± 1 au. The disk mass obtained with the stratified model is in agreement with the literature value (Zhang et al. 2021). Since the reduced chi-squared is smaller in the stratified case (see Table 2), we adopted it as the best-fit model.
![]() |
Fig. 3 Rotation curves of MWC 480 obtained from CO isotopologues. Left panel: Rotation curve obtained from 12CO data (red points), along with our best-fitting curve for the stratified model (blue line) and the isothermal one (orange line). Right panel: same for the 13CO data. |
4.4 IM Lup
IM Lup is a young pre-main sequence star (~1 Myr) located in the Lupus star-forming region at a distance of 158pc (Gaia Collaboration 2018). The dynamical stellar mass is estimated to be 1.1 M⊙ (Teague et al. 2021), and it hosts an unusually large disk, extending out to ≈300 au in the dust continuum and out to ≈1000 au in the gas (Cleeves et al. 2016). The dust continuum emission shows clear evidence of a spiral morphology, which may be triggered by gravitational instability (Huang et al. 2018). Cleeves et al. (2016) first estimated the disk mass from mm visibilities and found a massive disk of 0.2 M⊙. Verrios et al. (2022) claimed that the spiral structure of IM Lup could be generated by an embedded protoplanet. They performed numerical SPH simulations of planet-disk interaction and then post-processed them to compare their results with CO, dust, and scattered light emission. Interestingly, a high disk mass (~0.1M⊙) is required to match the scattered light image, so that the sub-micron sized grains could remain well coupled in the top layers of the disk. Cleeves et al. (2016) first estimated the disk scale radius Rc = 100 au by comparing SED to a simple tapered power-law density profile. Afterwards, Pinte et al. (2018) analyzed CO data and found that a tapered power law density profile with Rc = 284 au more optimally reproduces the data. They also analyzed the rotation curve of the disk and found that while the inner disk is in good agreement with Keplerian rotation around a 1 ± 0.1 M⊙ star, both the 12CO and the 13CO rotation curves become sub-Keplerian in the outer disk. The authors attributed this effect to the pressure gradient. Lodato et al. (2023) analyzed 12CO and 13CO rotation curves and fitted for star mass, disk mass and scale radius with an isothermal model. In particular, the authors found that for the rotation curves extracted with EDDY the best-fit are M★ = 1.012 ± 0.003 M⊙, Md = 0.096 ± 0.003 M⊙, Rc = 89 ± 1 au and DISCMINER are M★ = 1.02 ± 0.02, Md = 0.10 ± 0.01 M⊙, Rc = 66 ± 1 au. We underline that in this work the rotation curves have been obtained again, and they are different from those of Lodato et al. (2023). This is also true for the case of GM Aur.
Figure 4 shows both the isothermal and stratified fit. While for 12CO, both models describe the rotation curve well, for 13CO the isothermal model fails, since the velocity shift is so high that it cannot be explained just in terms of emitting surface. This difference is clearly visible when considering the , which for the stratified model is considerably smaller. The best-fit parameters for the isothermal model are M★ = 1.055 ± 0.002 M⊙, Md = 0.200 ± 0.003 M⊙, and Rc = 55 ± 1 au; while for the stratified model are M★ = 1.1994 ± 0.002 M⊙, Md = 0.106 ± 0.002 M⊙, and Rc = 115 ± 1 au. The effects of thermal stratification are visible in Fig. 8. At R ~ 250au, the difference in the data between 12CO and13 CO is on the order of ~10% and it significantly increases in the outer part. In that case, the stratified model is not able to explain that difference either. Izquierdo et al. (2023) pointed out that the emission from the outer disk is so diffuse that the retrieval of the emitting surface, as well as the velocity extraction, needs to be taken with care. This is possibly an effect of external photoevaporation. Indeed, despite the very weak external radiation field irradiating IM Lup, Haworth et al. (2017) showed that the disk is sufficiently large that the outer part, which is weakly gravitationally bound, can undergo photoevaporation.
Results of the fitting procedure and reduced chi-squared for the two different models: isothermal and stratified.
4.5 GM Aur
GM Aur is a T-Tauri star in the Taurus-Auriga star-forming region that hosts a transition disk. The stellar mass has been estimated dynamically to be M★ = 1.1 M⊙ by Teague et al. (2021), in agreement with previous measurements (Macías et al. 2018). Its CO morphology is very complex, showing spiral arms, tails, and interactions with the environments (Huang et al. 2021). From thermochemical models of MAPS data, Schwarz et al. (2021) obtained a disk mass of Md = 0.2 M⊙ and a scale radius of Rc = 111 au, making GM Aur a possibly gravitationally unstable disk. Lodato et al. (2023) performed a fitting procedure for the star mass, disk mass, and scale radius, using an isothermal model, finding that for GM Aur, the two CO lines provide inconsistent rotation curves, which cannot be attributed only to a difference in the height of the emitting layer. In addition, the authors provided a simple order-of-magnitude estimate of the expected velocity shift due to thermal stratification, concluding that the difference between the two rotation curves could not be explained by this effect. They drew this conclusion by taking into account the different temperature of the two molecules at their emission height, zi(R), given by Law et al. (2021). However, as shown in Appendix A, in the azimuthal velocity, it is not only essential to know the temperature at (R, z), but also its radial and vertical gradient at that location.
By analyzing the rotation curves of the two CO isotopolgues (Fig. 2), a systematic shift between 12CO and 13CO curves is clearly visible, which may possibly be attributed to thermal stratification. When we perform the fitting with the isothermal model, we obtain (as the best-fit parameters) M★ = 0.872 ± 0.003 M⊙, Md = 0.312 ± 0.003 M⊙, and Rc = 56 ± 1 au. This is in agreement with Lodato et al. (2023) and leads to a high (see Table 2). As a matter of fact, Fig. 5 shows that an isothermal model is not able to reproduce both 12CO and 13CO rotation curves. Conversely, when thermal stratification is taken into account, the two rotation curves are compatible and are in agreement with data, especially for R > 180 au. In this case, the best-fit value for the star mass is M★ = 1.128 ± 0.002 M⊙, which is in line with the literature values (Teague et al. 2021; Macías et al. 2018). As for the disk mass, the best-fit value is Md = 0.118 ± 0.002 M⊙. Finally, the best-fit value for the scale radius is Rc = 96 ± 1au, almost twice the value obtained with the isothermal model and in good agreement with Schwarz et al. (2021). A stratified model reproduces the difference between 12CO and 13CO rotation curves very well, as shown in Fig. 8, which leads to a significant decrease in the
value.
![]() |
Fig. 8 Relative difference between the squares of observed 12CO and 13CO rotation curves predicted by the thermally stratified model (blue line), the isothermal model (orange line), and the data (red dots). Except for AS 209, where this quantity is negative in the inner part, it is clearly visible that data are well reproduced by the stratified model. Indeed, the difference of speed between the two curves cannot be explained just in terms of different height. |
4.6 HD 163296
HD 163296 is one of the most well-studied Herbig Ae star system at millimeter wavelengths thanks to its relative close distance (d = 101pc) and bright disk. The disk presents several features that suggest ongoing planet formation, such as dust rings, deviations from Keplerian velocities due to gas pressure variations, “kinks” in the CO emission, and meridional flows (Isella et al. 2016, 2018; Pinte et al. 2018, 2023; Teague et al. 2018; Izquierdo et al. 2022, 2023; Calcino et al. 2022).
This system has also been extensively studied because there is considerable evidence to support a massive disk. By modeling the dust lines, Powell et al. (2019) found that the disk mass is Md = 0.21 M⊙. As for the scale radius, de Gregorio-Monsalvo et al. (2013) through radiative transfer modeling found that Rc = 125 au is the value that better reproduces dust and CO ALMA observations. Guidi et al. (2016) presented a multiwave-length ALMA and VLA study of the disk and, via a modeling of the visibilities, they found that the best-fit value of the scale radius is Rc = 118 au, in agreement with de Gregorio-Monsalvo et al. (2013).
When we fit data with a vertical isothermal model, we obtain (as the best-fit parameters) M★ = 1.842 ± 0.002 M⊙, Md = 0.124 ± 0.001 M⊙ and Rc = 38 ± 1 au. While the star mass is realistic, the scale radius is unrealistically small compared to the gas emission extent on the order of 400 au (Law et al. 2021). Additionally, the isothermal model is not able to reproduce the difference between the rotation curves of the two CO isotopo-logues (see Fig. 8), resulting in a relatively poor fit with a large . If we include the 2D thermal structure, the quality of the fit increases (see
in Table 2). In this case, the best-fit for stellar mass and disk mass does not change significantly (M★ = 1.948 ± 0.002 M⊙, Md = 0.134 ± 0.001 M⊙), while the scale radius does shift to Rc = 91 ± 1au. Comparing our result for the disk mass to the literature values, we observe that our fit gives a value that is roughly half. Figure 6 shows that both the isothermal and the stratified model aptly describe the rotation curve of 12CO and 13CO. However, the shift between them, presented in Fig. 8, is well recovered only by the stratified model, which offers only a partial explanation for the significant increase of the plotted quantity. The presence of pressure modulated substructures in the rotation curves (Izquierdo et al. 2023) impacts the quality of the fit (clearly visible in Fig. 8). One possible solution would be to include them in the fitting model.
4.7 AS 209
AS 209 is a young T-Tauri star in the Ophiucus star forming region (d ~ 121pc). The most recent stellar mass estimate is M★ = 1.14 M⊙ (Izquierdo et al. 2023). Fedele et al. (2018) gave an estimate for the scale radius Rc = 80au through the modeling of the mm visibilities . Afterwards, through thermochemical modeling, they found a dust mass of Mdust = 3.5 × 10−4 M⊙; with a gas-to-dust ratio of 100, this translates into Md = 0.0035 M⊙, which is in agreement with the recent value Md = 0.0045 M⊙ of Zhang et al. (2021). Interestingly, when inspecting the rotation curves of AS 209 (Fig. 2), the 13CO is slower compared to the 12CO, despite it being closer to the midplane. This trend is observed up to ~125 au. A possible explanation for this is the compactness of the disk, which makes more difficult to extract a precise emitting surface due to beam smearing. Indeed, line centroids from pixels near the center of the disk are an averaged composition of multiple surrounding velocities because of the limited resolution and the steepness of v(r). Since AS209 is the smallest disk in the sample, it is more prone to be affected by this in the largest fraction of its total extent compared to the other sources. When we fit with the isothermal model, we obtained, as the best-fit parameters, M★ = 1.272 ± 0.003 M⊙, Md = 0.042 ± 0.003 M⊙, and Rc = 45 ± 1 au. When we fit with the stratified model, we obtain as the best-fit parameters M★ = 1.311 ± 0.001 M⊙, and Rc = 126 ± 2 au, while for the disk mass we report a 3 − σ upper limit of Md = 0.00025 ± 0.00025 M⊙, since the best-fit parameter is compatible with zero. Both models are shown in Fig. 7. As for 12CO, the two models behave in the same way, showing little difference in the outer edge. Conversely, for 13CO the isothermal model works better in the inner part, where13 CO is slower, while in the outer part the stratified model describes well the rotation curve. According to the , the stratified model describes the data better (see Table 2).
5 Discussion
5.1 Thermal stratification in MAPS disks
Table 2 presents a summary of the findings of this study, comparing the isothermal model with the stratified one. It is evident from the results that the reduced χ2 value consistently decreases when employing the stratified model. This indicates that the inclusion of thermal stratification provides a more effective way of describing the observed data. In this context, MWC 480 is particularly interesting. Despite the small kinematic signatures of thermal stratification, as depicted in Fig. 8, the quality of the stratified fit is higher and it yields more reliable values for star mass, disk mass, and scale radius. On the opposite side, GM Aur is the system that shows the strongest effects of thermal stratification, given that the 12CO and 13CO systematically shifted over all the radial extent of the disk. The introduction of thermal stratification is able to reconcile these differences, reducing the by an order of magnitude. The only case where the stratified model encounters challenges in accurately describing both curves is in AS 209. This system is peculiar because the compactness of the disk influences the extraction of emission surfaces. Consequently, contrary to what expected, we observe that the 13CO rotates slower than the 12CO in the inner part. Despite that, the
is smaller when thermal stratification is taken into account.
5.2 Disk masses
In this paragraph, we aim to contextualize our work within the broader framework of disk mass estimation.
One solid tracer of the disk mass is the carbon dioxide HD, which is a good tracer of the disk gas because it follows the distribution of molecular hydrogen and its emission is sensitive to the total mass. The first detection of HD emission in a protoplanetary disk comes from Bergin et al. (2013) for TW Hya. Afterwards, the detection of HD J = 1–0 line was used to estimate disk mass of GM Aur. The HD based disk mass is 2.5– 20.4 × 10−2 M⊙ (McClure et al. 2016), in line with our estimate of 0.118 M⊙. Finally, the non-detection of HD in HD163296 (Kama et al. 2020) translates into an upper limit for the disk mass of 0.067 M⊙, which is almost half of the value we obtained in this work.
Another reliable method to trace the disk mass uses the N2H+. This molecule is a chemical tracer of CO-poor gas and can be used to measure the CO–H2 ratio and calibrate CO-based gas masses. By combining N2H+ with C18O, Trapman et al. (2022) estimated disk masses of three protoplanetary disks, including GM Aur. The value they obtained, 1.5–9.6 × 10−2 M⊙, is slightly higher compared to our estimate, but it is in an overall good agreement. This method has also been used to probe disk masses of protoplanetary disks in the Lupus star-forming region (Anderson et al. 2022).
In this context, it is worth mentioning observations of the 13C17O, a very rare CO isotopologue. Booth et al. (2019) observed this molecule in HD 163296 and this allows for a precise disk mass measurement to be obtained. These authors found that the disk mass that is better at reproducing observations is Md = 0.31 M⊙, which is discrepant with our inferred value.
As for the dust, its ability to trace the mass is discussed in the next subsection.
5.3 Gas-to-dust ratio
With the knowledge of the disk mass, it is possible to evaluate the gas-to-dust ratio, using Eq. (25) for the dust mass. The results are shown in Table 3. We found values between 100–250, within a factor of 2 from the usually assumed valued of 100. This is surprisingly, due to the assumptions we made to obtain the dust mass. Indeed, as we have already mentioned, the optically thin hypothesis for dust emission could lead to a difference of a more than a factor of 2 in the dust mass calculation (Guidi et al. 2016), underestimating it. In addition, the dust opacity could also vary of a factor of ~10, depending on the grain size and composition. Hence, it is significant overall that the inferred gas-to-dust ratio is so close to the standard value. As for AS 209, we estimated an upper limit for this quantity. Indeed, according to Veronesi et al. (in prep.), the minimum measurable mass with the rotation curve is 5% of the star mass. Taking this value as an upper limit for AS 209 disk mass, it is possible to give an upper limit for the gas-to-dust ratio.
![]() |
Fig. 9 Toomre Q parameter of the MAPS disks with the best-fit parameters of the stratified model. We excluded AS 209 because its best-fit disk mass is compatible with zero. |
5.4 Toomre Q
To investigate the presence of gravitational instability, we used our best-fit parameters for the stratified model to compute the Toomre parameter (Toomre 1964) which, based on the hypothesis of nearly Keplerian disk (κ ≃ Ω), is
(26)
where we used Eq. (1) for the surface density. According to the WKB quadratic dispersion relation (Lin & Shu 1964; Toomre 1964), the onset of the instability happens when Q ~ 1. Figure 9 shows the profile of the Q parameter for the MAPS sample, except for AS209, since its disk mass estimate is compatible with zero. Every disk is gravitationally stable, according to the Toomre criterion, since Q > 1. Interestingly, the two disks that showed spiral structures (IM Lup and GM Aur) have a Toomre profile that is lower than the others, with a minimum value of ~4 for GM Aur and ~6 for IM Lup. Lau & Bertin (1978) showed that a WKB description of gravitational instability can still be obtained under less restrictive conditions compared to the quadratic relation. They showed that disks that are locally stable according to the Q criterion might still generate large scale spiral waves. In general, other mechanisms could increase the critical value of the Toomre parameter, such as external irradiation (Lin & Kratter 2016; Löhnert et al. 2020) or dust-driven gravitational instability (Longarini et al. 2023b,a). Hence, we do not exclude the possibility that gravitational instability is at play in GM Aur and IM Lup.
6 Conclusions
The kinematic data of protoplanetary disks display velocity differences between 12CO and 13CO that cannot be explained through a vertically isothermal model, given the systematic shift between rotation curves of CO isotopologues. In this work, we predict how thermal stratification affects the density and the velocity field of a protoplanetary disk.
We used SPH simulations to test our model, finding excellent agreement, and then we applied it to the MAPS sample. We extracted the rotation curves of CO isotopologues (12CO and 13CO) and we carried out fitting for the star mass, disk mass, and scale radius – both with a vertically isothermal and a stratified model. The quality of the fit significantly improves when thermal stratification is taken into account and the best-fit parameters are more realistic and aligned with literature. All our results are summarized in Table 2.
Typically, when thermal stratification is considered, the best-fit value for the star mass tends to rise. This can be intuitively understood, as an isothermal model would favor a star mass that lies between that of 13CO and 12CO; its underestimation is due to the slower rotation of 12CO. Conversely, the stratified model encapsulates the difference between the two curves, mitigating the underestimation issue and resulting in a more accurate mass estimate. While an isothermal model provides a satisfactory fit at small radii, the fit worsens at large radii where the difference between 12CO and 13CO is larger. The fit tries to compensate for this by increasing the disk mass, most of which resides at large radii, thereby changing the predicted curve only in the outer parts of the disk. Ultimately, a more accurate description of the thermal structure through a stratified model leads to a realistic estimate of the scale radius.
We note that the inclusion of the vertical gradient of temperature into our model results in improved χ2 values across all systems under examination. This work demonstrates the impact of thermal stratification in the disk dynamics and highlights the importance of having a precise knowledge of the disk temperature structure to infer physical quantities such as the stellar mass, the disk mass and the disk scale radius in a meaningful manner. Thanks to the generality of our calculations, we can study the density and velocity profile of thermally stratified disk according different prescription of temperature simply.
Acknowledgements
This paper makes use of the following ALMA data: ADS/JAO.ALMA#2018.1.01055.L. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. This work has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement # 823823 (RISE DUSTBUSTERS project). G.R. is funded by the European Union under the European Union’s Horizon Europe Research & Innovation Programme No. 101039651 (discEvol) and by the Fondazione Cariplo, grant no. 2022-1217. S.F. is funded by the European Union (ERC, UNVEIL, 101076613). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council. Neither the European Union nor the granting authority can be held responsible for them. S.F. acknowledges financial contribution from PRIN-MUR 2022YP5ACE. CH is funded by a Research Training Program Scholarship from the Australian Government, and acknowledges funding from the Australian Research Council via DP220103767. M.B. and J.S. have received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (PROTO-PLANETS, grant agreement No. 101002188). The authors thank Pietro Curone and Claudia Toci for useful discussions.
Appendix A Computing the pressure gradient
Assuming a barotropic fluid, the pressure contribution to the rotation curve is
(A.1)
which can also be expressed as:
(A.2)
Assuming hydrostatic equilibrium in the vertical direction, we can write an explicit expression for the last term in Eq. (A.2):
(A.3)
and, assuming f (z = 0) = 1 (therefore T(R, z = 0) = Tmid(R)),
(A.4)
Deriving with respect to R we obtain:
(A.5)
where · = d/dR. Expanding the logarithmic term is useful since it is now clear that the azimuthal velocity depends both on the temperature (f) and on the temperature gradient along the radial direction. Moreover, it is easy to see the contribution of the vertical thermal stratification: in our model f and
are functions of z and thus they have to be integrated to compute azimuthal velocity; whereas in vertically isothermal models, this does not happen, since f = 1 and
. We note that it is not analytically possible to determine generally the sign of Eq. (A.5) and, thus, the effect that we see on the pressure gradient. However, for all the considered cases it never overcomes the gravitational contribution and it does lead to a faster deceleration in the rotation.
Values of parameters used in the estimate (B.1).
Appendix B Dullemond prescription
According to Law et al. (2021), the temperature prescription given by Dullemond et al. (2020) (Eq.(19)) fits the data well, but we note that in this case, f (R, z = 0) ≠ 1 and the temperature does not smoothly connect to its value at midplane since T(R, z = 0) ≠ Tmid(R). To evaluate this discrepancy, considering Tє (R) ≈ Tmid(R), we compute:
(B.1)
which is strongly dependent on Tmid, Tatm and α. As shown in Table B.1, the deviation of the actual midplane temperature from T(R, z = 0) as computed from the prescription given by Dullemond et al. (2020) is ≲ 10% in our regions of interest (R > 100 au). This discrepancy could have relevance for systems such as GM Aur, HD 163296 and MWC 480. To examine this further, we conducted fits for these systems using both Tmid (R) and T(R, z = 0) as the midplane temperature. It was observed that this choice does not significantly alter the results. Thus, Tmid ≈ T(T, z = 0) and f(R, z = 0) ≈ 1. We can consider Eqs. (19) and (20) as a good parameterization for the temperature.
Appendix C Corner plots
In Figs C.1 and C.2, we present the corner plots of the MCMC fitting procedure for the studied parameters under, respectively, the vertically isothermal and the stratified model.
![]() |
Fig. C.1 Corner plots of the MCMC fitting procedure according the vertically isothermal model. They show the distribution of the three relevant fitting parameters for the five disks of the MAPS large program. From top left to bottom: MWC 480, IM Lup, GM Aur, HD 163296, and AS 209. |
![]() |
Fig. C.2 Corner plots of the MCMC fitting procedure according the stratified model. They show the distribution of the three relevant fitting parameters for the five disks of the MAPS large program. From top left to bottom: MWC 480, IM Lup, GM Aur, HD 163296, and AS 209. |
References
- Abramowitz, M., & Stegun, I. A. 1970, Handbook of Mathematical Functions: with Formulas, Graphs, and Mathematical Tables [Google Scholar]
- ALMA Partnership (Brogan, C. L., et al.) 2015, ApJ, 808, L3 [Google Scholar]
- Anderson, D. E., Cleeves, L. I., Blake, G. A., et al. 2022, ApJ, 927, 229 [NASA ADS] [CrossRef] [Google Scholar]
- Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., & Dullemond, C. P. 2009, ApJ, 700, 1502 [CrossRef] [Google Scholar]
- Andrews, S. M., Huang, J., Pérez, L. M., et al. 2018, ApJ, 869, L41 [NASA ADS] [CrossRef] [Google Scholar]
- Bae, J., Teague, R., Andrews, S. M., et al. 2022, ApJ, 934, L20 [NASA ADS] [CrossRef] [Google Scholar]
- Bae, J., Teague, R., & Zhu, Z. 2021, ApJ, 912, 56 [NASA ADS] [CrossRef] [Google Scholar]
- Barraza-Alfaro, M., Flock, M., Marino, S., & Pérez, S. 2021, A&A, 653, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bergin, E. A., Cleeves, L. I., Gorti, U., et al. 2013, Nature, 493, 644 [Google Scholar]
- Bertin, G., & Lodato, G. 1999, A&A, 350, 694 [NASA ADS] [Google Scholar]
- Bollati, F., Lodato, G., Price, D. J., & Pinte, C. 2021, MNRAS, 504, 5444 [Google Scholar]
- Booth, A. S., Walsh, C., Ilee, J. D., et al. 2019, ApJ, 882, L31 [Google Scholar]
- Calahan, J. K., Bergin, E. A., Zhang, K., et al. 2021, ApJS, 257, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Calcino, J., Hilder, T., Price, D. J., et al. 2022, ApJ, 929, L25 [NASA ADS] [CrossRef] [Google Scholar]
- Chiang, E. I., & Goldreich, P. 1997, ApJ, 490, 368 [Google Scholar]
- Cleeves, L. I., Öberg, K. I., Wilner, D. J., et al. 2016, ApJ, 832, 110 [Google Scholar]
- Curone, P., Izquierdo, A. F., Testi, L., et al. 2022, A&A, 665, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- D’Alessio, P., Cantö, J., Calvet, N., & Lizano, S. 1998, ApJ, 500, 411 [Google Scholar]
- D’Alessio, P., Cantó, J., Hartmann, L., Calvet, N., & Lizano, S. 1999, ApJ, 511, 896 [CrossRef] [Google Scholar]
- Dartois, E., Dutrey, A., & Guilloteau, S. 2003, A&A, 399, 773 [CrossRef] [EDP Sciences] [Google Scholar]
- de Gregorio-Monsalvo, I., Ménard, F., Dent, W., et al. 2013, A&A, 557, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dipierro, G., Price, D., Laibe, G., et al. 2015, MNRAS, 453, L73 [NASA ADS] [CrossRef] [Google Scholar]
- Dullemond, C. P., Isella, A., Andrews, S. M., Skobleva, I., & Dzyurkevich, N. 2020, A&A, 633, A137 [EDP Sciences] [Google Scholar]
- Dutrey, A., Guilloteau, S., Piétu, V., et al. 2008, A&A, 490, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Facchini, S., Birnstiel, T., Bruderer, S., & van Dishoeck, E. F. 2017, A&A, 605, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fedele, D., Tazzari, M., Booth, R., et al. 2018, A&A, 610, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Galloway-Sprietsma, M., Bae, J., Teague, R., et al. 2023, ApJ, 950, 147 [NASA ADS] [CrossRef] [Google Scholar]
- Guidi, G., Tazzari, M., Testi, L., et al. 2016, A&A, 588, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hall, C., Dong, R., Teague, R., et al. 2020, ApJ, 904, 148 [NASA ADS] [CrossRef] [Google Scholar]
- Haworth, T. J., Facchini, S., Clarke, C. J., & Cleeves, L. I. 2017, MNRAS, 468, L108 [CrossRef] [Google Scholar]
- Hildebrand, R. H. 1983, QJRAS, 24, 267 [NASA ADS] [Google Scholar]
- Huang, J., Andrews, S. M., Pérez, L. M., et al. 2018, ApJ, 869, L43 [Google Scholar]
- Huang, J., Bergin, E. A., Öberg, K. I., et al. 2021, ApJS, 257, 19 [NASA ADS] [CrossRef] [Google Scholar]
- Isella, A., Guidi, G., Testi, L., et al. 2016, Phys. Rev. Lett., 117, 251101 [NASA ADS] [CrossRef] [Google Scholar]
- Isella, A., Huang, J., Andrews, S. M., et al. 2018, ApJ, 869, L49 [NASA ADS] [CrossRef] [Google Scholar]
- Izquierdo, A. F., Testi, L., Facchini, S., Rosotti, G. P., & van Dishoeck, E. F. 2021, A&A, 650, A179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Izquierdo, A. F., Facchini, S., Rosotti, G. P., van Dishoeck, E. F., & Testi, L. 2022, ApJ, 928, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Izquierdo, A. F., Testi, L., Facchini, S., et al. 2023, A&A, 674, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kama, M., Trapman, L., Fedele, D., et al. 2020, A&A, 634, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lau, Y. Y., & Bertin, G. 1978, ApJ, 226, 508 [NASA ADS] [CrossRef] [Google Scholar]
- Law, C. J., Teague, R., Loomis, R. A., et al. 2021, ApJS, 257, 4 [NASA ADS] [CrossRef] [Google Scholar]
- Lin, M.-K., & Kratter, K. M. 2016, ApJ, 824, 91 [Google Scholar]
- Lin, C. C., & Shu, F. H. 1964, ApJ, 140, 646 [Google Scholar]
- Lodato, G., Rampinelli, L., Viscardi, E., et al. 2023, MNRAS, 518, 4481 [Google Scholar]
- Löhnert, L., Krätschmer, S., & Peeters, A. G. 2020, A&A, 640, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Long, F., Pinilla, P., Herczeg, G. J., et al. 2018, ApJ, 869, 17 [Google Scholar]
- Longarini, C., Lodato, G., Toci, C., et al. 2021, ApJ, 920, L41 [NASA ADS] [CrossRef] [Google Scholar]
- Longarini, C., Armitage, P. J., Lodato, G., Price, D. J., & Ceppi, S. 2023a, MNRAS, 522, 6217 [NASA ADS] [CrossRef] [Google Scholar]
- Longarini, C., Lodato, G., Bertin, G., & Armitage, P. J. 2023b, MNRAS, 519, 2017 [Google Scholar]
- Lynden-Bell, D., & Pringle, J. E. 1974, MNRAS, 168, 603 [Google Scholar]
- Macías, E., Espaillat, C. C., Ribas, Á., et al. 2018, ApJ, 865, 37 [Google Scholar]
- McClure, M. K., Bergin, E. A., Cleeves, L. I., et al. 2016, ApJ, 831, 167 [Google Scholar]
- Miotello, A., Kamp, I., Birnstiel, T., Cleeves, L. C., & Kataoka, A. 2023, in Astronomical Society of the Pacific Conference Series, Protostars and Planets VII, eds. S. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 534, 501 [NASA ADS] [Google Scholar]
- Molyarova, T., Akimkin, V., Semenov, D., et al. 2017, ApJ, 849, 130 [Google Scholar]
- Montesinos, B., Eiroa, C., Mora, A., & Merín, B. 2009, A&A, 495, 901 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Öberg, K. I., Guzmán, V. V., Walsh, C., et al. 2021, ApJS, 257, 1 [CrossRef] [Google Scholar]
- Piétu, V., Dutrey, A., Guilloteau, S., Chapillon, E., & Pety, J. 2006, A&A, 460, L43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pinte, C., Price, D. J., Ménard, F., et al. 2018, ApJ, 860, L13 [Google Scholar]
- Pinte, C., Teague, R., Flaherty, K., et al. 2023, in Astronomical Society of the Pacific Conference Series, Protostars and Planets VII, eds. S. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 534, 645 [NASA ADS] [Google Scholar]
- Powell, D., Murray-Clay, R., Pérez, L. M., Schlichting, H. E., & Rosenthal, M. 2019, ApJ, 878, 116 [NASA ADS] [CrossRef] [Google Scholar]
- Price, D. J., Wurster, J., Tricco, T. S., et al. 2018, PASA, 35, e031 [Google Scholar]
- Rabago, I., & Zhu, Z. 2021, MNRAS, 502, 5325 [NASA ADS] [CrossRef] [Google Scholar]
- Ragusa, E., Dipierro, G., Lodato, G., Laibe, G., & Price, D. J. 2017, MNRAS, 464, 1449 [Google Scholar]
- Rosenfeld, K. A., Andrews, S. M., Hughes, A. M., Wilner, D. J., & Qi, C. 2013, ApJ, 774, 16 [Google Scholar]
- Rosotti, G. P., Teague, R., Dullemond, C., Booth, R. A., & Clarke, C. J. 2020, MNRAS, 495, 173 [Google Scholar]
- Schwarz, K. R., Calahan, J. K., Zhang, K., et al. 2021, ApJS, 257, 20 [NASA ADS] [CrossRef] [Google Scholar]
- Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
- Sierra, A., Pérez, L. M., Zhang, K., et al. 2021, ApJS, 257, 14 [NASA ADS] [CrossRef] [Google Scholar]
- Simon, M., Guilloteau, S., Beck, T. L., et al. 2019, ApJ, 884, 42 [NASA ADS] [CrossRef] [Google Scholar]
- Takeuchi, T., & Lin, D. N. C. 2002, ApJ, 581, 1344 [NASA ADS] [CrossRef] [Google Scholar]
- Teague, R., Bae, J., Bergin, E. A., Birnstiel, T., & Foreman-Mackey, D. 2018, ApJ, 860, L12 [NASA ADS] [CrossRef] [Google Scholar]
- Teague, R., Bae, J., Aikawa, Y., et al. 2021, ApJS, 257, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Teague, R., Bae, J., Andrews, S. M., et al. 2022, ApJ, 936, 163 [NASA ADS] [CrossRef] [Google Scholar]
- Terry, J. P., Hall, C., Longarini, C., et al. 2022, MNRAS, 510, 1671 [Google Scholar]
- Toci, C., Rosotti, G., Lodato, G., Testi, L., & Trapman, L. 2021, MNRAS, 507, 818 [NASA ADS] [CrossRef] [Google Scholar]
- Toomre, A. 1964, ApJ, 139, 1217 [Google Scholar]
- Trapman, L., Rosotti, G., Bosman, A. D., Hogerheijde, M. R., & van Dishoeck, E. F. 2020, A&A, 640, A5 [EDP Sciences] [Google Scholar]
- Trapman, L., Zhang, K., van’t Hoff, M. L. R., Hogerheijde, M. R., & Bergin, E. A. 2022, ApJ, 926, L2 [NASA ADS] [CrossRef] [Google Scholar]
- Trapman, L., Rosotti, G., Zhang, K., & Tabone, B. 2023, ApJ, 954, 41 [NASA ADS] [CrossRef] [Google Scholar]
- Veronesi, B., Ragusa, E., Lodato, G., et al. 2020, MNRAS, 495, 1913 [Google Scholar]
- Verrios, H. J., Price, D. J., Pinte, C., Hilder, T., & Calcino, J. 2022, ApJ, 934, L11 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, K., Booth, A. S., Law, C. J., et al. 2021, ApJS, 257, 5 [NASA ADS] [CrossRef] [Google Scholar]
The code is publicly available at https://github.com/crislong/DySc
All Tables
Velocity extraction method, orientation parameters, thermal parameters and emitting surfaces for 12CO and 13CO data of the MAPS disks.
Results of the fitting procedure and reduced chi-squared for the two different models: isothermal and stratified.
All Figures
![]() |
Fig. 1 Comparison between the analytical models and the simulation after eight outer orbits and at z = 20 au. Left panel: density field. Right panel: Pressure gradient term of the rotation curve. The red dots represent the simulation data, while the blue and orange line show the thermally stratified model and the isothermal one, respectively. The model with thermal stratification matches very well the simulation. |
In the text |
![]() |
Fig. 2 12CO and 13CO rotation curves of the MAPS disks extracted with DISCMINER. |
In the text |
![]() |
Fig. 3 Rotation curves of MWC 480 obtained from CO isotopologues. Left panel: Rotation curve obtained from 12CO data (red points), along with our best-fitting curve for the stratified model (blue line) and the isothermal one (orange line). Right panel: same for the 13CO data. |
In the text |
![]() |
Fig. 4 Same as Fig. 3 but for IM Lup. |
In the text |
![]() |
Fig. 5 Same as Fig. 3 but for GM Aur. |
In the text |
![]() |
Fig. 6 Same as Fig. 3, but for HD 163296. |
In the text |
![]() |
Fig. 7 Same as Fig. 3, but for AS 209. |
In the text |
![]() |
Fig. 8 Relative difference between the squares of observed 12CO and 13CO rotation curves predicted by the thermally stratified model (blue line), the isothermal model (orange line), and the data (red dots). Except for AS 209, where this quantity is negative in the inner part, it is clearly visible that data are well reproduced by the stratified model. Indeed, the difference of speed between the two curves cannot be explained just in terms of different height. |
In the text |
![]() |
Fig. 9 Toomre Q parameter of the MAPS disks with the best-fit parameters of the stratified model. We excluded AS 209 because its best-fit disk mass is compatible with zero. |
In the text |
![]() |
Fig. C.1 Corner plots of the MCMC fitting procedure according the vertically isothermal model. They show the distribution of the three relevant fitting parameters for the five disks of the MAPS large program. From top left to bottom: MWC 480, IM Lup, GM Aur, HD 163296, and AS 209. |
In the text |
![]() |
Fig. C.2 Corner plots of the MCMC fitting procedure according the stratified model. They show the distribution of the three relevant fitting parameters for the five disks of the MAPS large program. From top left to bottom: MWC 480, IM Lup, GM Aur, HD 163296, and AS 209. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.