---
title: Virga Cloud Modeling Framework
url: https://www.emergentmind.com/topics/virga-cloud-modeling-framework
type: topic
---

# Virga Cloud Modeling Framework

The Virga Cloud Modeling Framework is a Python-based, 1D equilibrium condensation suite that models cloud vertical structure, particle size distributions, aggregate morphology, and spectral properties in exoplanet and substellar atmospheres. It formalizes the steady-state balance of upward turbulent mixing and downward sedimentation for condensates, grounded in the Ackerman & Marley (2001) EddySed theory and extended to capture variable sedimentation efficiency, fractal aggregate aerosols, and multi-species cloud compositions. Virga is widely employed to interpret atmospheric spectra, benchmark microphysical models, and retrieve physically motivated cloud parameters from observations ranging from HST to JWST [2508.15102, 2603.13167, 2509.06708, 2512.04186, 2110.05903, 2202.01355].

## 1. Governing Physics and Core Equations

Virga operationalizes the vertical mass continuity for condensates as a balance between eddy diffusion and gravitational settling:
\[
-K_{zz} \frac{\partial q_t}{\partial z} - f_{\rm sed} w_* q_c = 0
\]
where $K_{zz}$ is the eddy diffusivity, $q_t$ the total (vapor + condensate) mass mixing ratio, $q_c$ the condensate mixing ratio, $w_*$ the convective velocity scale ($w_*=K_{zz}/H$), and $f_{\rm sed}$ is the dimensionless sedimentation efficiency [2603.13167, 2508.15102]. The first term quantifies upward turbulent transport and the second represents downward particle flux.

Steady-state closure integrates sources and sinks across vertical layers under local phase equilibrium, relating $q_v$ (vapor), $q_{\rm sat}$ (saturation, via Clausius-Clapeyron), and $q_c$ by
\[
q_v(z) = \min[q_t(z), q_{\rm sat}(z)], \quad q_c(z) = q_t(z) - q_v(z)
\]

Virga adopts a log-normal particle size distribution:
\[
\chi(r) = \frac{1}{\sqrt{2\pi} r \ln \sigma_g} \exp\left[ -\frac{1}{2} \left( \frac{\ln(r/r_g)}{\ln \sigma_g} \right)^2 \right]
\]
characterized by mean radius $r_g$ and width $\sigma_g$ [2603.13167].

The particle fall speed $v_{\rm sed}(r,z)$ is computed using physical drag laws across all Reynolds and Knudsen regimes [2508.15102], with transition between Stokes, slip-corrected, and kinetic drag formulations for spheres; aggregates follow specific free-molecular and continuum scaling [2509.06708].

## 2. Sedimentation Efficiency and Altitude Dependence

The sedimentation efficiency parameter $f_{\rm sed}$ is central to Virga's predictive power. It encodes the mass-weighted mean particle settling velocity relative to the convective velocity scale:
\[
f_{\rm sed}(z) = \frac{\int v_{\rm sed}(r,z)\, \chi(r) \, dr }{w_*(z)}
\]
Physically, larger $f_{\rm sed}$ yields rapid settling, large particles, and geometrically thin clouds; small $f_{\rm sed}$ yields slow settling, small particles, and vertically extended clouds [2603.13167, 2110.05903].

Virga supports both constant and variable altitude-dependent $f_{\rm sed}(z)$ parameterizations. The latter are necessary to reproduce cloud profiles from comprehensive microphysical models (e.g., CARMA). Example forms include exponential density dependence
\[
f_{\rm sed}(z) = \alpha \exp\left[ \frac{z-z_T}{6\beta H_0} \right] + \epsilon
\]
where $\alpha$, $\beta$, and $\epsilon$ are tunable parameters enforcing normalization and physical limits [2110.05903, 2202.01355]. The analytic structure of the governing equation is retained for computational efficiency.

Comparison with microphysical models shows that no constant-$f_{\rm sed}$ profile can reproduce the curvature and vertical extent of clouds under heterogeneous nucleation, necessitating flexible forms for retrievals and for jointly interpreting different spectral regions [2110.05903, 2202.01355, 2603.13167].

## 3. Treatment of Particle Morphology: Fractal Aggregates

Virga v2.0 and later generalize beyond compact spheres to arbitrary fractal aggregates [2509.06708, 2512.04186]. Particle morphology is controlled via the fractal dimension $D_f$ (1 for chains, 3 for spheres) and monomer radius $r_{\rm mon}$ or number $N_{\rm mon}$:
\[
N_{\rm mon} = k_0 \left( \frac{R_g}{r_{\rm mon}} \right)^{D_f}
\]
with $k_0$ a prefactor (Tazaki 2021). Aggregate growth is parametrized either with fixed $r_{\rm mon}$ or fixed $N_{\rm mon}$ [2509.06708].

Aggregate settling velocities are derived for both free-molecular and continuum drag regimes, explicitly accounting for the reduced density and cross-sectional scaling of fluffy structures. The closure ensures mass conservation and permits robust solution of the equilibrium fall-speed-radius condition.

Virga's optical property module for aggregates leverages the Modified Mean Field (MMF) theory as implemented in Optool to produce extinction, scattering, and asymmetry parameters across parameter grids of $D_f$, $r_{\rm mon}$, and $N_{\rm mon}$ [2509.06708, 2512.04186].

## 4. Optical Properties and Radiative Transfer Coupling

Virga computes wavelength-dependent cloud optical properties using Mie theory for spheres and MMF theory for aggregates [2508.15102, 2509.06708, 2512.04186]. Layer-by-layer extinction optical depth, single-scattering albedo, and asymmetry parameters are integrated over the local particle size distribution and passed to external radiative transfer codes (e.g., PICASO) for spectral synthesis [2508.15102].

For a fixed cloud mass, the spectral impact of morphology is regime dependent:
- **Rayleigh regime ($2\pi r/\lambda \ll 1$):** fluffier aggregates are more opaque.
- **Geometric regime ($2\pi r/\lambda \gg 1$):** compact particles dominate extinction.

Morphology alters transmission slope, NIR/IR band muting, and the strengths of diagnostic features such as the $10\,\mu$m silicate band [2512.04186, 2509.06708]. The MMF implementation yields efficiency within $\sim$1–2% of DDA models but with $\sim$150× computational speedup [2509.06708].

## 5. Input Parameters, Model Workflow, and Species Support

Virga requires as input vertical profiles of pressure, temperature, and eddy diffusivity. Cloud species selection is flexible, with support for at least 18 condensates (volatile ices, sulfides, silicates, metals, oxides) along with their corresponding vapor pressure parameterizations, refractive indices, and densities [2508.15102]. Users can add new species by supplying optical and vapor-pressure data.

A typical workflow involves:
1. Preparing atmospheric $P$–$T$ and $K_{zz}(z)$ profiles.
2. Selecting condensate species and cloud/morphology parameters ($f_{\rm sed}$, $D_f$, $r_{\rm mon}$).
3. Solving the mass-balance ODE to obtain $q_c(z)$ and particle size distributions.
4. Computing opacity profiles and synthesizing spectra using RT solvers [2508.15102].

Model output includes vertical distributions of $q_c(z)$, $r_g(z)$, number densities, optical depths $\tau(\lambda,z)$, and diagnostic properties for comparison with observed spectra.

## 6. Validation, Applications, and Limitations

Virga has been benchmarked against both observed exoplanet atmospheres (e.g., WASP-17b SiO$_2$ clouds) and microphysics models (CARMA, Sonora grid) [2508.15102, 2202.01355, 2509.06708]. For retrieval studies, $f_{\rm sed}$ and aggregate parameters are tuned to match spectral features, such as the strength and shape of $10\,\mu$m silicate absorption seen in JWST data. Very small $f_{\rm sed}$ and low sticking coefficients $s$ are required to reproduce high-altitude, small-grain clouds as implied by observations of WASP-107b, VHS-1256b, and YSES-1c [2603.13167].

Morphology critically impacts spectra, especially for broadband and short-wavelength datasets. Virga's ability to flexibly explore $f_{\rm sed}(z)$ and morphology parameter spaces makes it amenable to joint retrieval and forward-modeling efforts, but it should be cross-validated against full microphysics (e.g., growth and coagulation) where possible [2110.05903, 2603.13167].

The current model is parametric and does not resolve explicit kinetic microphysics (nucleation, coagulation cascades, break-up, or full radiative feedback), but it is extensible and computationally efficient for retrieval and sensitivity studies [2512.04186, 2509.06708]. Guidance is provided for physical parameter ranges and consistency checks to prevent unphysical solutions (e.g., porous aggregate compressibility limits).

## 7. Future Development and Recommendations

Ongoing and planned Virga developments include support for aggregate/hollow-core particles, coupling to 2-D GCM outputs for spatially resolved clouds, and retrieval-mode interfaces [2508.15102]. Community contributions are encouraged via the open-source code repository. Constraints from laboratory sticking coefficients for silicates, together with panchromatic spectra, are refining microphysics-to-observable connections.

Modelers are encouraged to utilize shape and $f_{\rm sed}(z)$ as free parameters, benchmark against laboratory and microphysical models, and leverage the rapid analytic solutions of Virga for comprehensive parameter space exploration and joint retrievals on multi-instrument datasets [2512.04186, 2603.13167, 2110.05903].

Source: https://www.emergentmind.com/topics/virga-cloud-modeling-framework