CARMApy: An Open-Source Python Framework for Simulating Microphysical Clouds in Planetary Atmospheres

Resumen Técnico de CARMApy: Un Marco de Trabajo en Python de Código Abierto para la Simulación de Nubes Microfísicas en Atmósferas Planetarias

Planteamiento del Problema
Las nubes ejercen influencias de primer orden sobre las propiedades atmosféricas planetarias, afectando las inferencias relativas a la composición (por ejemplo, las relaciones C/O y Mg/Si), la estructura térmica y la dinámica. Sin embargo, las incertidumbres en los procesos de formación de nubes se propagan significamente en estas inferencias. Los modelos de nubes existentes suelen depender de suposiciones simplificadoras que limitan su fidelidad física. Los modelos de "nube gris" priorizan la eficiencia sobre la complejidad física, mientras que los modelos de Condensación de Nubes en Equilibrio (ECC) asumen una saturación instantánea y requieren distribuciones de tamaño de partícula prescritas con parámetros ajustados (por ejemplo, $f_{sed}$) de un significado físico incierto. Incluso los modelos de Nubes Cinéticas por Método de Momentos (MMKC), que incluyen más física, todavía dependen de formas funcionales prescritas para las distribuciones de tamaño de partícula. Estas limitaciones dificultan la capacidad de reproducir las características espectrales observadas, particularmente con el advenimiento de los datos de alta resolución del JWST.

Aunque el Community Aerosol and Radiation Model for Atmospheres (CARMA) es un modelo de esquema de celdas (bin-scheme) bien validado, capaz de resolver distribuciones de tamaño de partícula arbitrarias desde primeros principios, su implementación heredada en Fortran presenta una barrera de entrada. El código requiere una experiencia significativa para operar, lo que limita su uso a un pequeño subconjunto de investigadores. Existe la necesidad de una interfaz de usuario amigable y accesible que mantenga el rigor físico de CARMA al tiempo que permite una aplicación más amplia en la comunidad de exoplanetas.

Metodología
Los autores presentan CARMApy, un envoltorio (wrapper) de Python de código abierto para ExoCARMA 2.0, una versión modernizada del código Fortran CARMA 3.0. CARMApy está diseñado para ser amigable al usuario, estar bien documentado y ser capaz de ejecutarse en hardware estándar con soporte de multihilo mediante OpenMP.

El núcleo del marco de trabajo es un modelo microfísico de esquema de celdas con pasos de tiempo que resuelve las distribuciones de tamaño de partícula a través de una rejilla de niveles de presión y masas de partículas definida por el usuario. El modelo simula los siguientes procesos microfísicos desde primeros principios:

• Nucleación: Se modelan tanto la nucleación homogénea (gas a partícula) como la heterogénea (condensación sobre semillas/CCN) utilizando la teoría de nucleación clásica, incluyendo factores de Zeldovich y factores de forma para los casos heterogéneos.
• Crecimiento y Evaporación: Modelados mediante la difusión de moléculas de vapor, corregidos por efectos de no-continuo (número de Knudsen), el efecto Kelvin (curvatura) y factores de ventilación para partículas en movimiento.
• Coagulación: Los núcleos de colisión contabilizan la difusión browniana, la mejora por difusión convectiva y la coalescencia gravitacional (asentamiento diferencial). El modelo asume que las partículas se adhieren o rebotan según eficiencias empíricas, pero actualmente no incluye fragmentación.
• Transporte: El transporte vertical incluye el sedimentación (calculada a través de regímenes laminar, transicional y turbulento) y la difusión de Eddies.

Numéricamente, el código emplea un método de división de operadores (operator-splitting). Dentro de cada paso de tiempo, los procesos se aplican secuencialmente: transporte vertical, coagulación, nucleación/crecimiento/evaporación y, finalmente, agotamiento de gas. El transporte y la difusión se resuelven utilizando el Método Parabólico Segmentado (PPM) para advectar partículas entre celdas. La coagulación y el crecimiento/evaporación se resuelven mediante esquemas semi-implícitos para asegurar la conservación de la masa y concentraciones de partículas no negativas.

CARMApy incluye interfaces integradas para:

• Inicialización: Establecer abundancias de equilibrio químico en la base de la nube utilizando pyFastChem y extender los perfiles atmosféricos de forma adiabática.
• Condensados: Soporte para hasta 10 condensados predeterminados (por ejemplo, TiO$_2$, Mg$_2$SiO$_4$) con propiedades de material especificadas por el usuario (densidad, tensión superficial, presión de vapor de saturación).
• Post-procesamiento: Generar espectros calculando la extinción, dispersión y parámetros de asimetría mediante PyMieScatt e interactuando con PICASO para la transferencia radiativa.
• Simulaciones 2D: Un modo para modelar perfiles de nubes que varían longitudinalmente (por ejemplo, para Júpiteres Calientes) mediante la advección de columnas de nubes alrededor del ecuador.

Resultados Clave
Los autores compararon CARMApy contra el código heredado ExoCARMA 1.0 utilizando un caso de prueba de una enana marrón de 2000 K con condensados de TiO$_2$ y Mg$_2$SiO$_4$.

• Exactitud: Los resultados resultaron consistentes con versiones anteriores. La masa total de la nube y del gas para cada especie concordó dentro de diferencias relativas de $\sim 3 \times 10^{-4}$ y $\sim 3 \times 10^{-5}$, respectivamente. Los espectros de emisión mostraron errores absolutos máximos en la temperatura de brillo de $\sim 0.3$ K (error relativo $\sim 10^{-4}$), atribuidos a la propagación del error de punto flotante a lo largo de decenas de miles de pasos de tiempo.
• Rendimiento: El envoltorio de Python, aprovechando el multihilo de OpenMP, se ejecutó aproximadamente 1.9 veces más rápido en modo de un solo hilo y 3.8 veces más rápido en modo multihilo en comparación con el código heredado ExoCARMA 1.0.
• Salidas: El artículo demuestra la generación de diversos productos de datos, incluyendo perfiles de distribución de tamaño de partícula, mapas de tasas microfísicas (nucleación, crecimiento, evaporación), perfiles de abundancia de gas agotado y espectros de transmisión y emisión con asimetría de limbo para simulaciones 2D.

Significado y Reivindicaciones
El artículo afirma que CARMApy democratiza con éxito el acceso a la modelización microfísica de alta fidelidad. Al envolver un robusto código Fortran de medio siglo de antigüedad en una interfaz Python moderna, los autores permiten que investigadores con mínima experiencia previa en modelización de nubes o en Fortran puedan simular la física compleja de las nubes. El marco se presenta como una herramienta para mejorar la interpretación de observaciones de alta resolución (como las del JWST) al resolver distribuciones de tamaño de partícula arbitrarias y capturar efectos de no-equilibrio que los modelos más simples pasan por alto. Los autores enfatizan que, si bien el código es un paso significativo hacia adelante, se necesita trabajo futuro para restringir las propiedades de los materiales de las partículas de las nubes en regímenes de exoplanetas e incorporar física adicional como agregados fractales y nubes radiativamente activas. El código se libera como software de código abierto con tutoriales y documentación para facilitar la adopción por parte de la comunidad.