---
title: 'DiskMINT: A Self-Consistent Disk Modeling Framework'
url: https://www.emergentmind.com/topics/diskmint
type: topic
---

# DiskMINT: A Self-Consistent Disk Modeling Framework

DiskMINT, short for **Disk Model for INdividual Targets**, is an open-source, self-consistent thermochemical disk-modeling framework for inferring the gas and dust properties of protoplanetary disks from continuum and CO isotopologue observations. In its initial form, it was introduced as a Python tool to estimate gas disk masses from optically thin $\mathrm{C^{18}O}$ rotational lines by coupling an SED-derived dust structure to vertical hydrostatic equilibrium, isotope-selective chemistry, and grain-surface conversion of CO into $\mathrm{CO_2}$ ice [2307.02657]. Subsequent developments extended the framework to radially and vertically decoupled gas and dust distributions for individual disks and, later, to DiskMINT-GARDEN, a public grid of 480 self-consistent models plus a fast machine-learning inference tool for large ALMA samples [2509.15487] [2606.26332].

## 1. Origins and scientific scope

DiskMINT was developed to address a persistent interpretation problem in protoplanetary-disk observations: faint CO isotopologue emission had often been taken to imply either very low gas masses or severe CO depletion. The framework instead models the vertical density, temperature, and chemistry together, so that the emitting layer is not imposed ad hoc but derived from the same physical assumptions used to interpret the observables [2307.02657].

The original formulation focused on individual targets and used long-wavelength SED constraints, multiple dust grain sizes, and $\mathrm{C^{18}O}$ line emission to infer disk gas masses. In later work, DiskMINT was extended to infer not only total gas mass but also its radial distribution, while allowing gas and dust surface densities to be spatially decoupled both radially and vertically [2509.15487]. DiskMINT-GARDEN generalized this target-specific framework into a survey-scale inference system by combining the DiskMINT physics with a precomputed model grid and an XGBoost regressor that maps a small set of ALMA observables to $(M_{\rm gas}, \epsilon, R_c)$ [2606.26332].

In the DiskMINT literature, “self-consistent” has a specific meaning. The pressure-supported vertical density profile, the dust and gas temperature fields, and the chemistry that sets the molecular abundances are computed together and mutually constrained; no piece is held fixed or imposed ad hoc. The resulting thermochemical structure places the CO-emitting layers, their temperatures and optical depths, and the continuum optical depths in a single internally consistent model [2606.26332].

## 2. Dynamical and thermal structure

A central component of DiskMINT is vertical hydrostatic equilibrium. In the grid formulation, the standard relation
$$
\frac{dP}{dz} = -\rho \Omega_K^2 z,\qquad \Omega_K=\sqrt{\frac{GM_\star}{r^3}},
$$
is solved together with the energy balance so that the gas and dust temperature fields and the density structure are updated iteratively [2606.26332]. The local scale height follows
$$
H=\frac{c_s}{\Omega_K},\qquad c_s=\sqrt{\frac{k_B T}{\mu m_H}}.
$$
In the extended individual-target framework, DiskMINT solves the same vertical hydrostatic problem with the actual $T_g(r,z)$ rather than assuming vertical isothermality, and evaluates the midplane scale height as $H(r)=c_s(r,\mathrm{mid})/\Omega_K(r)$ [2509.15487].

Disk surface densities are represented with tapered power-law profiles. DiskMINT-GARDEN adopts
$$
\Sigma(r)=\Sigma_1 (r/1\,{\rm au})^{-\gamma}\exp[-(r/R_c)^{2-\gamma}]
$$
with $\gamma=1$, $R_{\rm out}=1000$ au, and the same profile for gas and dust with a constant $\epsilon$ across radii [2606.26332]. In the radially decoupled framework, the corresponding similarity-solution form is
$$
\Sigma(r)=\Sigma_0 (r/R_c)^{-p}\exp[-(r/R_c)^{2-p}],
$$
with $\Sigma_g(r)$ and $\Sigma_d(r)$ treated as independent functions constrained by spatially resolved ALMA radial profiles [2509.15487].

Dust settling is included in the more advanced DiskMINT models through a steady-state turbulent diffusion equation,
$$
\frac{\partial}{\partial z}\left[\ln\left(\frac{\rho_d(a)}{\rho_g}\right)\right]
= - \frac{\Omega^2 \tau_s}{D} z,
$$
with
$$
\tau_s = \frac{\rho_\bullet a}{\rho_g c_s},\qquad D=\frac{\alpha_v c_s h}{Sc},
$$
and $Sc \approx 1$ [2509.15487]. This produces distinct vertical distributions for grains of different sizes, rather than a single well-mixed dust layer.

The thermal treatment differs slightly across published implementations. In the 2023 RU Lup application, the gas temperature was derived from a cross-section-weighted mean of dust temperatures plus viscous heating,
$$
T_g(r,z)=\left[T_{g,d}^4(r,z)+T_{g,v}^4(r,z)\right]^{1/4},
$$
with
$$
T_{g,v}=\left[\frac{3GM_\star \dot{M}_{\rm acc}}{8\pi \sigma_{\rm SB} r^3}
\left(1-\sqrt{\frac{r_\star}{r}}\right)\right]^{1/4}
$$
[2307.02657]. In the IM Lup extension, the gas temperature in the dense molecular layer where $\mathrm{C^{18}O}$ originates is computed to be in thermal equilibrium with the dust grains, which is treated as appropriate for a massive disk [2509.15487].

## 3. Chemistry and radiative transfer

DiskMINT’s reduced CO chemistry is tailored to the dominant pathways controlling CO in disks. The network includes CO, $\mathrm{CO_2}$, $\mathrm{H_2O}$, atomic C/C$^+$, O/O$_2$, $e^-$, $\mathrm{N_2}$, $\mathrm{HCO^+}$, $\mathrm{N_2H^+}$, and hydrocarbons needed to close dominant pathways [2606.26332]. Three processes are treated as critical.

First, CO freeze-out onto grains is included, with the accretion rate per unit volume written as
$$
R_{\rm freeze}=n_{\rm CO} v_{\rm th} S \pi a^2 n_{\rm gr},
$$
where
$$
v_{\rm th}=\sqrt{\frac{8k_B T}{\pi m_{\rm CO}}}.
$$
Thermal desorption, photodesorption, and cosmic-ray–induced desorption are also included [2606.26332].

Second, DiskMINT incorporates grain-surface conversion of CO ice to $\mathrm{CO_2}$ ice on water-rich ice mantles. This process moves the effective CO snow surface upward by tens of K compared to pure CO ice, so the gaseous CO layer becomes thinner and located higher, around $T \sim 35$ K in the self-consistent vertical temperature structure [2606.26332]. In the individual-target studies, this reaction is treated as the principal grain-surface pathway lowering gaseous CO above the midplane snowline [2509.15487].

Third, the framework includes isotope-selective photodissociation. Because self-shielding and line overlap are weaker for rare isotopologues, $\mathrm{C^{18}O}$ becomes selectively underabundant in UV-exposed layers. DiskMINT therefore computes vertical abundance structures in which photodestruction competes with freeze-out and grain-surface chemistry, rather than assuming a fixed isotopologue abundance [2606.26332].

Continuum radiative transfer is solved with RADMC-3D. The transfer equation along a ray is
$$
\frac{dI_\nu}{ds} = -\kappa_\nu \rho I_\nu + \kappa_\nu \rho S_\nu,
$$
and for an isothermal slab the emergent intensity reduces to
$$
I_\nu \approx B_\nu(T_{\rm dust})(1-e^{-\tau_\nu}).
$$
DiskMINT explicitly computes $\tau_\nu$ from the converged vertical structure, rather than relying on the optically thin estimator
$$
M_{\rm dust}\approx \frac{F_\nu d^2}{\kappa_\nu B_\nu(T_{\rm dust})},
$$
which is valid only when $\tau_\nu \ll 1$ [2606.26332].

For line radiative transfer, the individual-target workflows use non-LTE calculations with LIME to generate synthetic cubes and images of CO isotopologues [2307.02657] [2509.15487]. By contrast, DiskMINT-GARDEN uses LTE line luminosities for inference, with LTE verified in the grid because gas densities in the emitting layers exceed $n_{\rm crit}\sim 10^4$–$3\times 10^4\ {\rm cm^{-3}}$ for $\mathrm{C^{18}O}$ $J=2$–1 and $J=3$–2 [2606.26332].

## 4. Observables and inference workflows

DiskMINT supports two inference regimes. The first is a full forward-modeling mode for individual targets. In the 2023 workflow, the long-wavelength SED is fitted first; DiskMINT then constructs a self-consistent dust-plus-gas structure, computes $\mathrm{C^{18}O}$ abundances with the reduced chemistry, performs non-LTE line radiative transfer with LIME, and compares the resulting total luminosity, velocity profile, and radial intensity profile to the data [2307.02657]. In the structured 2025 formulation, $\Sigma_d(r)$ and $\Sigma_g(r)$ are updated iteratively via
$$
\Sigma_{d,i+1}(r)=\Sigma_{d,i}(r)\times R_{{\rm cont},i}(r),\qquad
R_{{\rm cont},i}(r)\equiv \frac{I_{\rm obs}(r)}{I_{{\rm model},i}(r)},
$$
and
$$
\Sigma_{g,i+1}(r)=\Sigma_{g,i}(r)\times R_{{\rm C^{18}O},i}(r),\qquad
R_{{\rm C^{18}O},i}(r)\equiv \frac{I_{\rm obs}(r)}{I_{{\rm model},i}(r)},
$$
while solving vertical hydrostatic equilibrium and dust settling at each step [2509.15487].

The second regime is DiskMINT-GARDEN, which converts the same physics into a survey-scale estimator. The public grid contains 480 self-consistent models spanning stellar mass $0.1$–$2.0\,M_\odot$, gas disk mass $10^{-5}$–$10^{-1}\,M_\star$, dust-to-gas ratio $0.003$–$0.1$, and characteristic radius $10$–$300$ au [2606.26332]. For each model it provides the Band 6 continuum luminosity $L_{\rm mm}$ near 234 GHz, the $\mathrm{C^{18}O}$ $J=2$–1 and $J=3$–2 line luminosities $L_{\rm C^{18}O}$, and the dust size metric $R_{{\rm dust},90}$, the radius enclosing 90% of the continuum flux [2606.26332].

The machine-learning stage uses a gradient-boosted decision tree regressor in log space to map
$$
O=(L_{\rm C^{18}O},L_{\rm mm},R_{{\rm dust},90})
$$
at fixed $M_\star$ to
$$
\Theta=(M_{\rm gas},\epsilon,R_c).
$$
Training uses a 90/10 split with early stopping and regularization; the validation performance is reported as $R^2 \approx 0.99$ with residual scatter $\sim 0.15$ dex across parameters [2606.26332]. Statistical uncertainties are propagated by Monte Carlo sampling in observable space and are typically $\lesssim 5\%$, increasing to $\lesssim 25\%$ at ${\rm S/N}\sim 3$, while the dominant uncertainty is systematic and is conservatively taken as a factor of $\sim 2$ on inferred $M_{\rm gas}$ [2606.26332].

A recurrent inference issue is optical-depth degeneracy. $L_{\rm C^{18}O}$ increases with $M_{\rm gas}$ but saturates when $\tau_{\rm line}$ rises, especially for small $R_c$; $R_{{\rm dust},90}$ is therefore used to distinguish compact optically thick disks from extended optically thin ones, while $L_{\rm mm}$ constrains $M_{\rm dust}$ and, with $\epsilon$, informs shielding and temperature in the emitting layers [2606.26332].

## 5. Applications to RU Lup, IM Lup, and ALMA survey samples

DiskMINT was first demonstrated on RU Lup, a high-accreting star whose disk had previously been inferred to have $M_{\rm gas}\sim 1.5\times 10^{-3}\,M_\odot$ and gas-to-dust ratio $\sim 4$. The best-fit DiskMINT model with vertical hydrostatic equilibrium yielded $M_{\rm gas}\sim 1.2\times 10^{-2}\,M_\odot$, while a parametric Gaussian vertical distribution that better matched the IR SED yielded $M_{\rm gas}\sim 2.1\times 10^{-2}\,M_\odot$ [2307.02657]. Both models reproduced the total $\mathrm{C^{18}O}$ luminosities for $J=2$–1 and $J=3$–2, the $J=2$–1 velocity profile, and the deprojected radial intensity profile. The comparison indicated that the inferred gas mass was about an order of magnitude higher than earlier DALI-grid estimates for the same source [2307.02657].

In IM Lup, DiskMINT was used in its radially and vertically decoupled form. By fitting the multi-wavelength SED, the millimeter continuum, and the $\mathrm{C^{18}O}$ radial emission profiles, the study derived a gas disk mass of $0.02$–$0.08\,M_\odot$, with a well-mixed best fit of $\sim 3.9\times10^{-2}\,M_\odot$ and a structured best fit of $\sim 2.8\times10^{-2}\,M_\odot$ [2509.15487]. The gas disk was found to extend beyond the dust disk, with $R_{\rm gas}\approx 481$ au versus $R_{\rm dust}\approx 264$–$300$ au, and the outer disk was inferred to be drift-dominated, with $\epsilon(r)\approx 0.01$–$0.02$ between $\sim 40$ and $300$ au [2509.15487]. The same work concluded that this dust-to-gas ratio is likely insufficient for strong streaming-instability-driven clumping in the outer disk.

DiskMINT-GARDEN was applied to archival ALMA observations of 34 disks with $>3\sigma$ detections of Band 6 continuum and $\mathrm{C^{18}O}$: 14 from AGE-PRO, 14 large disks from exoALMA plus ancillary measurements, and 6 from DSHARP, MAPS, and other programs [2606.26332]. The inferred gas masses agreed with independent dynamical disk-mass measurements and HD-based masses to within a factor $\lesssim 2$ for reliable cases, consistent with the adopted systematic floor [2606.26332]. When compared with DALI-based chemical modeling, DiskMINT recovered similar disk masses without invoking large uniform CO or elemental C depletion. Instead, grain-surface CO$\rightarrow \mathrm{CO_2}$ conversion and the upward-shifted effective CO snow surface were found to reduce CO columns in the emitting layers naturally [2606.26332]. The survey application therefore concluded that extant data suggest little chemical processing due to disk evolutionary processes.

## 6. Assumptions, limitations, and significance

DiskMINT assumes axisymmetric, flared disks under vertical hydrostatic equilibrium. The survey grid does not include explicit non-axisymmetric substructures or spirals, and warns that models approaching gravitational instability, as flagged by the Toomre parameter $Q$ at $R_{{\rm dust},90}$, require caution [2606.26332]. In compact, high-surface-density disks, both continuum and $\mathrm{C^{18}O}$ can become optically thick, and line luminosities then become less sensitive to mass and more sensitive to temperature, increasing the systematic floor [2606.26332].

The treatment of dust microphysics varies by implementation. The RU Lup study used dsharp\_opac with an ISM-like mixture of 64% astronomical silicates and 36% graphite [2307.02657], whereas the later work adopted the DIANA standard mixture and size distribution, with opacities computed via Optool [2509.15487]. This suggests that DiskMINT’s mass estimates are not tied to a single dust-opacity package, although $L_{\rm mm}$ and $R_{{\rm dust},90}$ remain sensitive to dust assumptions.

The grid-based inference assumes co-spatial gas and dust surface densities with a global $\epsilon$, and does not include explicit turbulent mixing or radial drift, even though targeted DiskMINT runs can treat radially varying $\epsilon$ [2606.26332]. The dense-layer gas temperature is tied closely to dust in the individual-target applications; the IM Lup study explicitly notes that omission of more complete gas heating and cooling introduces at most a factor $\lesssim 2$ mass uncertainty for lower-density disks [2509.15487]. For this reason, lower-mass or more UV-irradiated systems may require fuller thermochemical treatments or additional tracers.

Within those limits, DiskMINT’s principal significance is methodological. It recasts $\mathrm{C^{18}O}$ from a potentially misleading proxy into a robust tracer of total gas mass, and in extended form of radial gas distribution, provided that the line is interpreted in a model where vertical structure, dust temperatures, shielding, freeze-out, and grain-surface chemistry are solved together [2307.02657] [2509.15487]. DiskMINT-GARDEN further shows that the same physics can be compressed into a reproducible survey pipeline based only on widely available ALMA observables, without requiring $\mathrm{N_2H^+}$ [2606.26332].

Planned extensions follow directly from the current limitations. The 2026 framework identifies additional tracers such as HD and $^{13}\mathrm{CO}$, simultaneous inference of vertical CO emitting surfaces, radially varying $\epsilon$ and substructure, and coupling to spatially resolved line moments as natural directions for reducing degeneracies and improving robustness, especially in compact or optically thick disks [2606.26332].

Source: https://www.emergentmind.com/topics/diskmint