Self-Consistent Evolution Models Show Weak Double-Diffusive Mixing in Jupiter and Saturn

Resumo Técnico: Modelos de Evolução Autoconsistentes Mostram Mistura Doble-Difusiva Fraca em Júpiter e Saturno

Problema
A estrutura interna dos planetas gigantes, particularmente a existência de núcleos "difusos" ou estendidos em Júpiter e Saturno, permanece um desafio significativo para a ciência planetária. Embora a convecção de dupla difusão tenha sido proposta como um mecanismo para redistribuir elementos pesados de um núcleo denso para o envelope sobrejacente, sua eficiência em contextos evolutivos realistas é incerta. Estudos anteriores basearam-se em estruturas composicionais idealizadas e não trataram o transporte de dupla difusão de forma autoconsistente dentro de cálculos de evolução planetária. Além disso, há uma discrepância entre simulações hidrodinâmicas locais, que frequentemente preveem a fusão rápida de camadas e mistura eficiente, e os requisitos de modelos evolutivos globais que devem reproduzir as restrições observacionais atuais.

Metodologia
Os autores implementam um tratamento autoconsistente do transporte de dupla difusão dentro do código de evolução planetária APPLE.

• Condições Iniciais: O estudo utiliza modelos de interior pós-formação para Júpiter (baseados em Stevenson et al. 2022) e Saturno (baseados em Bodenheimer et al. 2025). Esses modelos apresentam gradientes composicionais estendidos ("núcleos difusos") resultantes da disrupção parcial e evaporação de planetesimais acretados durante a formação, em vez de fronteiras de núcleo nítidas.
• Física de Transporte: Os autores incorporam a prescrição de convecção em camadas derivada por Spruit (2013). Esta abordagem modela o transporte através de escadas convectivas como um processo de difusão aprimorado. As difusividades efetivas composicional ($D_{\text{eff}}$) e térmica ($\kappa_{T,\text{eff}}$) são derivadas com base em regimes assintóticos onde o transporte é independente da espessura da camada, escalando com a razão de estabilidade ($R_\rho$) e difusividades microscópicas.
• Simulação Evolutiva: Os modelos são evoluídos por 4,56 Gyr (a idade do Sistema Solar). O código resolve simultaneamente a evolução térmica e o transporte composicional usando uma grade de massa Lagrangiana. Regiões instáveis à convecção de reversão (critério de Ledoux) são tratadas com a teoria de comprimento de mistura padrão, enquanto regiões estáveis de Ledoux com $1 < R_\rho < \tau^{-1/2}$ utilizam a prescrição de dupla difusão.
• Análise de Sensibilidade: Para testar a robustez dos resultados, os autores realizam simulações onde as difusividades efetivas são aumentadas em fatores de 100 e 1000 em relação aos valores fiduciais de Spruit (2013).

Resultos Principais

• Mistura Ineficiente: As simulações revelam que a convecção de dupla difusão produz uma mistura muito limitada entre o interior profundo e o envelope. Ao longo da história de resfriamento de 4,56 Gyr, menos de $\sim 1 M_\oplus$ de material pesado é redistribuído em Júpiter e Saturno.
  • No modelo fiducial, apenas $\approx 0,11 M_\oplus$ (Júpiter) e $\approx 0,07 M_\oplus$ (Saturno) são redistribuídos.
  • Mesmo com difusividades aumentadas por um fator de 1000, a massa redistribuída permanece modesta ($\approx 2 M_\oplus$ para Júpiter e $\approx 1,3 M_\oplus$ para Saturno).
• Mecanismo de Autolimitação: A ineficiência surge de um ciclo de feedback negativo. O transporte térmico aprimorado modifica a estratificação térmica, reduzindo a razão de estabilidade ($R_\rho$) e, consequentemente, suprimindo as difusividades efetivas. Isso evita uma mistura descontrolada, garantindo que mesmo coeficientes de transporte irrealisticamente grandes falhem em homogeneizar os gradientes composicionais primordiais.
• Restrição Energética: Os autores identificam uma restrição energética fundamental. Na aproximação de Boussinesq usada em muitas simulações hidrodinâmicas, a mistura é energeticamente barata porque a espessura da camada é pequena em comparação com a escala térmica. No entanto, em interiores de planetas gigantes compressíveis, a região estável abrange uma fração substancial do raio planetário ($H \sim 0,25\text{--}0,5 R_p$), tornando a energia necessária para misturar o gradiente composicional comparável ao orçamento de energia térmica disponível. Consequentemente, o trabalho de flutuabilidade disponível para impulsionar o transporte composicional é severamente limitado.

Significância e Alegações
O artigo afirma que a convecção de dupla difusão, quando tratada de forma autoconsistente em um contexto evolutivo global, não pode erodir significativamente os gradientes composicionais gerados durante a formação planetária. As distribuições de elementos pesados observadas em Júpiter e Saturno provavelmente requerem mecanismos de transporte adicionais ou vias de formação que produzam uma mistura inicial mais ampla do que os modelos de acreção padrão preveem.

Os autores sugerem que os núcleos "difusos" podem ser um produto direto do processo de formação (por exemplo, mistura inicial mais extensa durante a ablação de planetesimais) ou o resultado de grandes eventos colisionais, em vez de um produto da mistura de dupla difusão pós-formação. O estudo destaca uma limitação crítica na aplicação de resultados hidrodinâmicos locais (frequentemente baseados na aproximação de Boussinesq) à evolução global de planetas gigantes, enfatizando que o orçamento de energia do interior profundo é o principal regulador da eficiência de mistura.