A Comprehensive Sulfur Chemistry Network Including Excited S(1D) and SO(1{\Delta}) for the XODIAC Photochemical Model: Accounting for Missing Sulfur Processes in Venus and Exo-Venus Analogs

Resumo Técnico: Rede Química de Enxofre Abrangente Incluindo S(¹D) e SO(¹∆) para o Modelo Fotoquímico XODIAC

Título Original: A Comprehensive Sulfur Chemistry Network Including Excited S(¹D) and SO(¹∆) for the XODIAC Photochemical Model: Accounting for Missing Sulfur Processes in Venus and Exo-Venus Analogs
Autores: Priyankush Ghosh et al.
Data: Abril de 2026 (Versão de rascunho)

1. O Problema

A química do enxofre é fundamental para a estrutura atmosférica, formação de nuvens e composição de Vênus e exoplanetas do tipo Vênus (exo-Vênus). No entanto, os modelos fotoquímicos atuais enfrentam limitações críticas:

• Dados Cinéticos Incompletos: Reações envolvendo espécies de enxofre no estado fundamental e excitado sob condições de alta temperatura e rico em CO₂ (como em Vênus) possuem dados cinéticos incertos ou inexistentes.
• Espécies Ignoradas: Espécies excitadas de longa vida, como o enxofre atômico no estado singlete S(¹D) e o monóxido de enxofre no estado excitado SO(¹∆), são frequentemente negligenciadas, apesar de sua potencial reatividade superior.
• Discrepâncias Observacionais: Modelos existentes falham em reproduzir as abundâncias observadas de polissulfetos (especificamente S₃ e S₄) na atmosfera de Vênus, sugerindo a ausência de processos químicos ou fontes de enxofre não contabilizadas.
• Acoplamento com CO₂: A reação entre enxofre e dióxido de carbono (CO₂), o reservatório dominante de carbono em Vênus, carece de dados de taxas de reação precisos derivados de primeiros princípios.

2. Metodologia

Os autores desenvolveram uma abordagem computacional de múltiplas etapas para estender a rede química do modelo fotoquímico unidimensional XODIAC:

• Cálculos de Estrutura Eletrônica:
  • Utilizaram a teoria do funcional de densidade (DFT) com o funcional híbrido WB97X-D/def2-TZVP para otimizar geometrias e frequências.
  • Refinaram energias eletrônicas de ponto único no nível de alta precisão CCSD(T)/aug-cc-pV(T+d)Z.
  • Investigaram superfícies de energia potencial (PES) para as reações S(³P) + CO₂ (estado fundamental) e S(¹D) + CO₂ (estado excitado).
  • Aplicaram correções de acoplamento spin-órbita e energias de excitação baseadas em dados experimentais do NIST para garantir a precisão dos estados excitados.
• Cálculos Cinéticos:
  • Utilizaram o solver de Equação Mestre (MESS) para calcular coeficientes de taxa dependentes da temperatura (150–2000 K) e pressão (1 bar).
  • Empregaram a teoria do estado de transição (TST) e a teoria do espaço de fase (PST) para canais com e sem barreira.
  • Ajustaram os coeficientes de taxa à equação de Arrhenius modificada.
• Termodinâmica:
  • Geraram coeficientes polinomiais da NASA (7 termos) para as espécies excitadas S(¹D) e SO(¹∆) usando o módulo Arkane (RMG), validando-os contra dados do banco de dados Burcat para espécies no estado fundamental.
• Modelagem Fotoquímica:
  • Implementaram a nova rede (XODIAC-2025.v2) no modelo 1D XODIAC.
  • Simularam Vênus e três análogos de exo-Vênus sob diferentes perfis de temperatura-pressão (T-P), fluxos estelares (1x e 1000x solar) e condições de fronteira inferior (com e sem fonte de enxofre atômico de 1 ppm).

3. Contribuições Principais

1. Novos Dados Cinéticos: Determinação de parâmetros cinéticos para as reações de S(³P) + CO₂ e S(¹D) + CO₂, incluindo a formação de SO(³Σ), SO(¹∆) e CO.
2. Expansão da Rede Química: Inclusão sistemática de espécies excitadas (S(¹D), SO(¹∆)) e novos canais fotoquímicos de dissociação que geram esses estados, com taxas de ramificação estimadas estatisticamente.
3. Banco de Dados Termodinâmico: Fornecimento de coeficientes da NASA para S(¹D) e SO(¹∆), permitindo sua integração consistente em modelos de equilíbrio e fotoquímica.
4. Diagnóstico de Fontes Faltantes: Identificação de que a discrepância nos polissulfetos em Vênus não é resolvida apenas pela nova cinética, mas requer uma fonte adicional de enxofre reduzido na baixa atmosfera.

4. Resultados

• Mecanismo de Reação S + CO₂:

  • Estado Fundamental (S(³P)): A reação possui uma barreira de entrada alta (~284 kJ/mol), tornando-a lenta. Ocorre via um complexo intermediário instável [SOCO].
  • Estado Excitado (S(¹D)): A formação do complexo de reação é sem barreira devido à natureza nucleofílica do S(¹D). A reação procede mais rapidamente, embora ainda enfrente barreiras de ativação significativas para a dissociação em produtos.
  • Conclusão Cinética: Sob condições de Vênus, a reação direta S + CO₂ é cineticamente desfavorável devido às altas barreiras de ativação, limitando seu impacto direto na química global em comparação com outros caminhos.

• Simulações em Vênus:

  • A inclusão da nova rede (v2) produziu efeitos menores nas abundâncias acima de 60 km, com pequenas variações fotoquímicas.
  • O modelo ainda falha em reproduzir as abundâncias observadas de S₃ e S₄.
  • Solução Proposta: A introdução de uma fonte de fronteira inferior de 1 ppm de enxofre atômico (representando processos não resolvidos na superfície ou na atmosfera profunda) aumentou as abundâncias de S₃ e S₄ em 1 a 2 ordens de grandeza, trazendo o modelo para dentro das margens de erro observacionais. Isso sugere a existência de um fluxo de enxofre reduzido não contabilizado ou vias de formação mais eficientes.

• Simulações em Exo-Vênus:

  • Condições Isotérmicas/Fluxo Baixo: As mudanças foram modestas.
  • Alta Irradiação + Iso-termia de Alta Altitude: Em atmosferas fortemente irradiadas (1000x fluxo solar) com uma isotermia de alta altitude, a química excitada torna-se mais relevante.
  • Impacto da Fonte de Enxofre: Em análogos de exo-Vênus, a adição de 1 ppm de enxofre atômico na superfície causou alterações pronunciadas não apenas nos polissulfetos, mas também em SO, H₂S, H₂SO₄, S₂O e CO em altitudes médias. O modelo híbrido (M9) mostrou as maiores variações, indicando que a presença de enxofre superficial pode reestruturar significativamente a química atmosférica de exoplanetas rochosos.

5. Significância e Conclusões

• Validação de Modelos: O estudo demonstra que, embora a inclusão de estados excitados e novas cinéticas seja crucial para a completude teórica, ela por si só não resolve as discrepâncias observacionais em Vênus.
• Importância de Fontes Profundas: A necessidade de uma fonte de enxofre atômico na baixa atmosfera para reproduzir S₃ e S₄ em Vênus aponta para lacunas no entendimento dos processos de degaseificação vulcânica, interações superfície-atmosfera ou química em camadas profundas.
• Implicações para Exoplanetas: Para análogos de exo-Vênus, especialmente aqueles sob alta irradiação estelar, a química do enxofre é altamente sensível a fontes de superfície e processos de estado excitado. Ignorar esses fatores pode levar a previsões errôneas sobre a composição atmosférica e a detecção de bioassinaturas ou geoassinaturas.
• Futuro: O trabalho enfatiza a necessidade de modelos fotoquímicos 1D mais completos que integrem fontes de enxofre profundo e vias fotoquímicas para estados excitados, além de mais estudos laboratoriais sob condições extremas de Vênus.

Em suma, este artigo fornece uma base teórica robusta para a química do enxofre em atmosferas ricas em CO₂, destacando que a "química faltante" em Vênus pode residir não apenas em reações gasosas desconhecidas, mas em fluxos de matéria da superfície ou atmosfera profunda que alimentam a formação de polissulfetos.