Resolving core-envelope degeneracies in giant planets and sub-Neptunes: constraining the equivalence in the presence of dilute gradients

Resumo Técnico: Resolvendo Degenerescências de Núcleo-Envelope em Gigantes Planetários e Sub-Netunos

Definição do Problema
A estrutura interna de gigantes planetários e sub-netunos representa um problema altamente degenerado e multidimensional na ciência planetária. O reconhecimento de núcleos "difusos" ou diluídos — onde elementos pesados estão distribuídos em gradientes em vez de fronteiras nítidas — e de mantos voláteis massivos complicou a interpretação dos raios observados ($R_p$) e massas ($M_p$). Superar o "vale do raio" exige distinguir entre planetas com atmosferas espessas de H/He, oceanos de água supercrítica ou composições mistas. Os códigos de evolução atuais (códigos da classe Henyey como MESA, APPLE, CEPAM) são computacionalmente caros, o que muitas vezes impede as varreduras densas de parâmetros necessárias para mapear todo o espaço de degenerescências estruturais estáticas, particularmente em relação à largura dos gradientes composicionais ($\sigma$). Além disso, a distribuição espacial específica de elementos pesados afeta significamente o raio total e os momentos de gravidade ($J_n$), porém a topologia funcional desses gradientes permanece mal restringida apenas pela observação.

Metodologia
Os autores apresentam o fuzzycore, um integrador estrutural hidrostático estático, multi-zona e de código aberto, projetado para a exploração rápida de modelos diretos de degenerescências internas. Diferente dos códigos evolutivos, o fuzzycore desacopla o perfil estrutural do transporte de energia dependente do tempo, da mistura convectiva e das curvas críticas de equilíbrio de fase.

Principais componentes metodológicos incluem:

• Equilíbrio Hidrostático Multi-Zona: O framework resolve equações 1D de simetria esférica para conservação de massa e equilíbrio hidrostático. Ele adapta dinamicamente a resolução da grade para resolver gradientes de densidade íngremes através de interfaces.
• Formulação da Equação de Estado (EOS):
  • Núcleo Misto: Utiliza uma lei de volume aditiva (lei de Amagat) para misturar tabelas ab initio de ferro e rocha multifásica (separando MgSiO3 em MgO e SiO2 sob altas pressões) para garantir uma estratificação de densidade monotônica.
  • Envelope Fluido: Modela o envelope como uma mistura ternária (H, He, elementos pesados) usando mistura linear ideal de volumes específicos. Elementos pesados ($Z$) são representados por uma EOS de H2O de ampla faixa para capturar efeitos de peso molecular médio ($\mu$) e compressibilidade através da transição da atmosfera prístina para as camadas basais enriquecidas.
  • Manto Condensado: Suporta um manto de água com fases distintas (Gelo VII/X ou superiônico) para sub-netunos ricos em voláteis, acoplado termodinamicamente à base do envelope.
• Gradientes Composicionais: Em vez de fronteiras discretas, a fração de massa de elementos pesados $Z(x)$ é modelada via uma função de decaimento Gaussiano contínua. Para lidar com a estabilidade numérica conforme os gradientes se tornam mais íngremes ($\sigma \to 0$), o framework emprega média analítica de sub-grade para conservar o orçamento de elementos pesados e evitar picos de densidade não físicos.
• Integração Numérica: O envelope é discretizado em uma zona convectiva de "escada" totalmente misturada. Um método robusto de passo adiabático (adiabat stepper) usando regressões lineares de mínimos quadrados locais sobre um KD-Tree do espaço de fase 3D ($P, T, S$) integra os estados termodinâmicos.
• Disparo de Fronteira (Boundary Shooting): Um método de disparo de fronteira aninhado (método de Brent) resolve o problema de valor de contorno, iterando sobre pressões de interface para garantir a continuidade de massa, pressão e temperatura entre a atmosfera, o manto e o núcleo. Ele suporta modos de alvo de massa e alvo de gravidade superficial.
• Guardas de Estabilidade: Inclui uma "guarda isotérmica" para evitar artefatos de resfriamento numérico em transições de fase de baixa pressão e amostragem adaptativa para envelopes ultra-finos.

Contribuições Principais e Resultados

1. Robustez da Forma Funcional: O benchmarking contra um perfil de Júpiter derivado de evolução (Vazan et al. 2018) demonstra que, quando as metalicidades de fronteira do envelope e os orçamentos integrados de elementos pesados são correspondidos, o raio macroscópico é robusto à topologia exata do gradiente. Seis famílias paramétricas distintas (Gaussiana, sigmoide, exponencial, bilinear, passo único e escadas hidrodinâmicas digitalizadas) produziram uma dispersão de raio de apenas 2,3%. Isso valida o uso de parametrizações analíticas suaves (como o decaimento Gaussiano) para mapear degenerescências.
2. Atlas de Degenerescência de Sub-Netunos: O framework foi aplicado a um planeta fixo de $8 M_\oplus$ para gerar um "atlas de degenerescência de mundo-água" através de uma grade de $15 \times 15 \times 12$ de massa de água, massa de envelope e largura de diluição ($\sigma$).
   • $\mu$-Compactação: O estudo quantifica como a redistribuição de elementos pesados de um núcleo compacto para um envelope diluído aumenta o peso molecular médio, compactando fisicamente a atmosfera. Esta "$\mu$-compactação" pode suprimir a expansão térmica, reduzendo o raio planetário em até 45% para orçamentos de elementos pesados fixos.
   • Cobertura Térmica (Thermal Blanketing): Envelopes ricos em metais atuam como limitadores térmicos, aumentando a diferença de temperatura entre a superfície e o núcleo em cerca de $\sim 5000$ K, embora este seja um achado estrutural estático, não um histórico de resfriamento.
   • Limites Observacionais: Para um $(M_p, R_p)$ observado fixo, a largura de diluição $\sigma$ é essencialmente não restringida. O atlas revela uma família contínua de arquiteturas estáticas (variando massa de água, massa de envelope e $\sigma$) que reproduzem o mesmo raio, destacando que $(M_p, R_p)$ sozinho não pode resolver a estrutura interna.
3. Aplicação a Alvos Específicos:
   • GJ 1214 b (Rico em voláteis): Modelado como um corpo água-rocha compacto com um envelope fino e altamente enriquecido, convergindo para um raio dentro de 1% das observações.
   • Kepler-11 e (Rico em gás/Super-puff): Modelado como um núcleo massivo com um envelope diluído e fortemente poluído. O modelo demonstra que um gradiente de elementos pesados profundamente suspenso é necessário para "lastrear" o envelope massivo até o raio compacto observado, contrapondo a inflação esperada de um envelope de H/He prístino.
4. Acoplamento Atmosfera-Interior: O framework foi projetado para servir como um backend para recuperações atmosféricas (ex: acoplamento com HADES). Ele pode ligar os priors de metalicidade atmosférica de JWST/Ariel a arquiteturas de interiores profundos ao resolver para gravidade superficial e condições de contorno autossuficientes.

Significância e Escopo
O artigo posiciona o fuzzycore não como um substituto para códigos evolutivos totalmente implícitos, mas como uma ferramenta complementar necessária. Ao remover o custo computacional do misturamento dependente do tempo e do transporte de energia, o fuzzycore permite o mapeamento sistemático e de alta resolução das degenerescências estruturais estáticas que seriam excessivamente caras para explorar com resolvedores da classe Henyey.

Os autores enfatizam que o fuzzycore identifica configurações hidrostáticas mecanicamente permitidas. Ele não determina se essas configurações são termodinamicamente consistentes com a física de desmistura ou hidrodinamicamente estáveis contra erosão turbulenta ao longo de escalas de tempo geológicas; essas questões permanecem no domínio da modelagem evolutiva e de estudos de equilíbrio de fase ab initio. No entanto, o framework fornece o "atlas de degenerescência" essencial necessário para interpretar dados observacionais atuais e futuros, demonstrando que a distribuição espacial de elementos pesados é uma dimensão crítica, porém observacionalmente degenerada, na caracterização planetária.