Fractal Aggregate Aerosols in the Virga Cloud Code I: Model Description and Application to a Benchmark Cloudy Exoplanet

Abstract: We introduce new functionality to treat fractal aggregate aerosol particles within the Virga cloud modeling framework. Previously, the open source cloud modeling code Virga (Batalha et al. 2025), the Python version of EddySed (Ackerman & Marley, 2001), assumed spherical particles to compute particle mass and size distributions throughout the atmosphere. The initial release of Virga also assumed spherical particles to compute Mie scattering properties, which include the single scattering albedo, asymmetry parameter, and optical depth as a function of particle radius and composition. However, extensive evidence from Solar system aerosols, astrophysical disks and dust, and Earth climate studies suggests that non-spherical aggregate particles are common compared to idealized compact spherical particles. Following recent advances in microphysical and opacity modeling, we implement a simple parametrization for dynamical and optical (modified mean field theory) effects of fractal aggregate particles into Virga. We then use this new functionality to perform a case study using basic planetary parameters similar to the well-characterized, aerosol-laden mini-Neptune GJ 1214 b, using KCl clouds made of aggregate particles. We choose KCl to most directly explore comparisons to previous studies. We demonstrate 1) how our method compares to previous fractal aggregate particle treatments and 2) how our new fractal treatment affects theoretical spectra of cloudy atmospheres. Overall, our model is faster and more flexible for a wider range of parameter space than previous studies. We explore the limitations of our modeling set-up and offer guidance for future investigations using our framework.