LavAtmos 2.0: Incorporating Volatiles Species in Vaporization Models

Resumen Técnico: LavAtmos 2.0 – Incorporación de especies volátiles en modelos de vaporización

1. Planteamiento del Problema
Los exoplanetas rocosos de periodo ultra-corto (USP) están sometidos a una irradiación estelar extrema que puede fundir sus superficies, creando océanos de lava. La interacción directa entre estos océanos y la atmósfera permite un equilibrio termoquímico que vincula la composición del interior planetario con la atmósfera. Sin embargo, los modelos previos de vaporización de lava (como la versión anterior de LavAtmos) han asumido simplificaciones críticas:

• Han ignorado la presencia de elementos volátiles (Carbono, Hidrógeno, Nitrógeno, Azufre y Fósforo) en la atmósfera, asumiendo que se han perdido debido a la erosión atmosférica.
• Han tratado la adición de volátiles como un proceso pasivo (sumar abundancias calculadas de lava a una atmósfera volátil predefinida), sin considerar cómo los volátiles afectan las propias reacciones de vaporización de la lava.
• Observaciones recientes del telescopio JWST de 55-Cnc e sugieren la existencia de una atmósfera volátil significativa, lo que hace necesario reevaluar cómo estos gases interactúan con la vaporización de la roca fundida.

2. Metodología
Los autores han desarrollado LavAtmos 2.0, una actualización del código de vaporización termoquímica LavAtmos. La metodología clave incluye:

• Integración de FastChem: Se ha integrado el código de equilibrio químico de fase gaseosa FastChem para expandir el número de especies gaseosas consideradas de 31 a 523.
• Cálculo de Equilibrio Acoplado: El código resuelve simultáneamente las leyes de acción de masas y conservación de masa. Calcula la presión parcial de oxígeno ($P_{O_2}$) que satisface tanto las reacciones de vaporización de la lava (líquido-gas) como el equilibrio químico de la atmósfera (gas-gas).
• Algoritmo Iterativo: El modelo utiliza un bucle iterativo para encontrar la abundancia de oxígeno atmosférico que coincide con la $P_{O_2}$ impuesta por la fusión de la lava, asegurando que la ley de acción de masas se cumpla en todo el sistema.
• Escenarios de Prueba:
  • Atmósferas "Puras": Se probaron atmósferas dominadas individualmente por C, H, N, S o P.
  • Atmósferas "Complejas": Mezclas dominadas por un elemento pero conteniendo todos los otros volátiles principales.
  • Caso de Estudio 55-Cnc e: Se aplicó el modelo a dos composiciones atmosféricas propuestas para este exoplaneta (una pobre en carbono y otra rica en carbono) basadas en datos de JWST, comparando los resultados con el método de "vaporización sin volátiles" (el enfoque anterior).

3. Contribuciones Clave

• Desarrollo de LavAtmos 2.0: La primera herramienta que modela de manera autoconsistente la interacción entre la vaporización de un océano de lava y una atmósfera rica en múltiples especies volátiles (C, H, N, S, P).
• Corrección de Sesgos en Modelos Previos: Demuestra que los métodos anteriores (sumar volátiles a la vaporización) subestiman sistemáticamente las presiones parciales de las especies vaporizadas y alteran incorrectamente el estado de oxidación de la atmósfera.
• Nueva Perspectiva sobre 55-Cnc e: Proporciona un marco para interpretar los datos de JWST considerando que la lava puede actuar como un reservorio masivo de oxígeno, influyendo en la relación C/O de la atmósfera.

4. Resultados Principales

• Aumento de la Vaporización: La inclusión de elementos volátiles en los cálculos de equilibrio aumenta significativamente las presiones parciales de las especies vaporizadas de la lava (como SiO, SiO₂, Na y K) en comparación con los modelos sin volátiles.
  • Mecanismo: Los volátiles reaccionan con el $O_2$ para formar especies oxigenadas (ej. $H_2O$, $CO_2$), reduciendo la presión parcial de $O_2$ libre. Según la ley de acción de masas, una menor $P_{O_2}$ desplaza el equilibrio de vaporización de la lava hacia la producción de más especies gaseosas (ej. $SiO_2(l) \leftrightarrow SiO(g) + \frac{1}{2}O_2(g)$).
  • Efecto del Carbono: El Carbono tiene el efecto más fuerte en aumentar la vaporización debido a su alta propensión a formar especies oxigenadas estables ($CO$, $CO_2$), reduciendo drásticamente la $P_{O_2}$. El Nitrógeno tiene el efecto más débil.
• Aumento del Oxígeno Elemental: Contrario a la intuición, la presencia de volátiles conduce a un aumento en la abundancia total de oxígeno elemental en la atmósfera. Esto se debe a que la reducción de la $P_{O_2}$ fomenta la vaporización de óxidos de la lava, liberando más oxígeno al sistema.
• Impacto en Especies Volátiles: Las presiones parciales de especies volátiles clave ($CO_2$, $H_2O$, $SO_2$) difieren en hasta seis órdenes de magnitud entre el nuevo modelo y el método antiguo, especialmente en regímenes donde la atmósfera volátil domina sobre el vapor de roca.
• Relación C/O y 55-Cnc e:
  • El modelo LavAtmos 2.0 predice una relación C/O más baja que los modelos tradicionales cuando hay un océano de lava presente, debido al aporte masivo de oxígeno desde la lava.
  • Para 55-Cnc e, una atmósfera rica en carbono en equilibrio con un océano de lava podría producir una relación C/O que coincide con las observaciones de JWST, lo cual no ocurre con una atmósfera pobre en carbono.
• Detección de Océanos de Lava: Las especies como SiO, SiO₂, Na y K podrían ser más fáciles de detectar en atmósferas volátiles de lo que se pensaba anteriormente, ya que sus presiones parciales aumentan, contrarrestando el efecto de supresión espectral que se creía que tenían los volátiles.

5. Significado e Implicaciones

• Reevaluación de Observaciones: Los resultados sugieren que las firmas espectrales de exoplanetas rocosos calientes (como 55-Cnc e) deben interpretarse considerando la interacción dinámica con océanos de lava. Ignorar este acople lleva a errores en la estimación de la temperatura superficial, la composición atmosférica y la estructura térmica.
• Nuevos Indicadores de Océanos de Lava: Una relación C/O baja en un exoplaneta USP con atmósfera volátil podría servir como un indicador indirecto de la presencia de un océano de lava en la superficie, siempre que el océano sea lo suficientemente profundo para amortiguar la fugacidad de oxígeno del sistema.
• Futuro de la Modelización: El trabajo subraya la necesidad de modelos termoquímicos autoconsistentes que incluyan tanto la química de la roca fundida como la de los volátiles atmosféricos para comprender la historia de formación y evolución química de los planetas rocosos.

En conclusión, LavAtmos 2.0 demuestra que los elementos volátiles no son meros espectadores en la atmósfera de los planetas de lava, sino que activamente modifican la química de vaporización, aumentando la abundancia de especies rocosas detectables y alterando fundamentalmente el estado redox de la atmósfera.