First-Principles Turbulence-Driven Deflagration-to-Detonation Transition Mechanism for Near-Chandrasekhar Mass White Dwarf Progenitors

Resumo Técnico: Mecanismo de Transição Deflagração-Detonção Impulsionado por Turbulência a Partir de Primeiros Princípios para Progenitores de Anãs Brancas com Massa Próxima ao Limite de Chandrasekhar

Declaração do Problema
As supernovas do tipo Ia (SNe Ia) servem como velas padrão cosmologicamente calibráveis críticas e fontes primárias de elementos do grupo do ferro no universo. No entanto, o mecanismo físico que governa sua explosão, especificamente a transição de uma deflagração subsônica para uma detonação supersônica (DDT) em anãs brancas (WDs) com massa próxima ao limite de Chandrasekhar (near-$M_{Ch}$), permaneceu sem resolução. Simulações anteriores basearam-se em prescrições ad hoc, tipicamente acionando a detonação quando a frente de chama atingia um limiar crítico de densidade calibrado empiricamente (por exemplo, $\rho_{DDT} \approx 10^7$ g/cm$^3$). Essa abordagem introduziu incertezas significativas, impedindo uma compreensão a partir de primeiros princípios de como as propriedades dos progenitores (como densidade central e topologia de ignição) mapeiam para resultados observáveis. A falta de um mecanismo físico de DDT dificultou a capacidade de inferir características dos progenitores a partir de curvas de luz e espectros observados, uma limitação cada vez mais problemática para levantamentos de cosmologia de precisão de próxima geração.

Metodologia
Os autores apresentam as primeiras simulações hidrodinâmicas globais e tridimensionais de progenitores de WDs near-$M_{Ch}$ que incorporam um mecanismo ab initio validado em laboratório para DDT impulsionado por turbulência (tDDT), em vez de impor um gatilho de detonação via prescrição.

• Hidrodinâmica e Condições Iniciais: Utilizando o código FLASH 4.3, a equipe simulou seis modelos distintos abrangendo um fator de seis na densidade de ignição central ($\rho_c$ de $1.0 \times 10^9$ a $6.0 \times 10^9$ g/cm$^3$) e variando topologias de ignição (deslocamentos de núcleo único e ignição multiponto com 100 núcleos). Os perfis iniciais de WD foram mapeados a partir de modelos hidrostáticos 1D com gradientes de temperatura adiabáticos representando a fase de cozimento lento (simmering).
• O Mecanismo tDDT: O cerne da metodologia é a implementação do critério de deflagração de Chapman-Jouguet (CJ). Ao contrário de modelos anteriores que exigiam combustão distribuída ou choques com intensidade de detonação, este mecanismo postula que uma detonação inicia quando a velocidade da chama turbulenta ($S_T$) excede a velocidade de deflagração de CJ ($S_{CJ}$). Essa condição é satisfeita quando a turbulência local dobra a superfície da chama suficientemente para criar uma fuga de pressão. O critério é avaliado localmente usando um operador diferencial calibrado ($OP_2$) para extrair a intensidade turbulenta solenoidal ($u'$) do fluxo resolvido. A detonação é acionada quando $u'$ excede uma intensidade crítica ($U_{CJ}$) sobre uma escala de comprimento crítica ($L_{CJ}$), derivada da escala de turbulência de Kolmogorov e propriedades da chama.
• Pós-processamento: As simulações utilizaram partículas traçadoras Lagrangianas passivas para registrar histórias termodinâmicas. Estas foram pós-processadas usando o código Torch (rede de 489 isótopos) para calcular rendimentos de nucleossíntese. Os rendimentos 3D foram mapeados para uma grade de velocidade homóloga 1D para cálculos de transferência radiativa usando o código SuperNu, gerando espectros sintéticos para comparação com observações.

Resultados Chave

• Iniciação Robusta de Detonação: Todos os seis modelos, apesar de abrangerem densidades centrais e geometrias de ignição diversas, alcançaram com sucesso a detonação via mecanismo tDDT. A transição foi imediata, ocorrendo enquanto os penachos de Rayleigh-Taylor (RT) ainda estavam subindo, em vez de após uma pré-expansão extensa.
• Convergência de Resultados Espectrais: Apesar de variações significativas nas condições iniciais, todos os modelos convergiram para espectros sintéticos quase idênticos no pico de luminosidade. A classificação espectral usando o código SNID identificou consistentemente todos os modelos como correspondentes ao subtipo sobre-luminoso SN 1999aa (Ia-99aa), independentemente do deslocamento de ignição ou da densidade central. Isso sugere que o mecanismo tDDT impulsiona progenitores diversos para uma configuração de detonação comum.
• Rendimentos de Nucleossíntese:
  • Produção de $^{56}$Ni: A natureza imediata do tDDT levou a altos rendimentos de $^{56}$Ni ($0.79$–$1.11 M_\odot$), excedendo os de SNe Ia normais ($\sim 0.4$–$0.6 M_\odot$). Os rendimentos diminuíram monotonicamente com o aumento da densidade central do progenitor devido à captura eletrônica aprimorada.
  • Elementos Estáveis do Grupo do Ferro (IGEs): Modelos de alta densidade produziram rendimentos significativamente maiores de IGEs estáveis e ricos em nêutrons (por exemplo, $^{54}$Fe, $^{58}$Ni, $^{55}$Mn) em comparação com modelos de baixa densidade. A estratificação desses elementos mostrou inversões distintas em casos de alta densidade, com espécies ricas em nêutrons dominando o ejetos mais internos.
  • Elementos de Massa Intermediária (IMEs): Os rendimentos de IMEs (Si, S, Ca) foram suprimidos em uma ordem de grandeza em comparação com modelos padrão de detonação atrasada. Isso é atribuído à detonação imediata consumindo combustível em altas densidades, levando a uma incineração quase completa até o Equilíbrio Estatístico Nuclear (NSE) em vez de combustão incompleta.
  • $^{44}$Ti: Os rendimentos permaneceram baixos ($\sim 10^{-5} M_\odot$), distinguindo esses eventos de cenários de dupla detonação com massa sub-Chandrasekhar.
• Morfologia da Chama: As simulações revelaram que a geometria de ignição influencia significativamente a estrutura da chama pré-detonação (por exemplo, espessura do caule, desenvolvimento de penachos RT), mas não altera a configuração final de detonação ou a aparência espectral.

Significado e Alegações
O artigo alega fornecer o primeiro caminho fisicamente motivado e autoconsistente para detonação atrasada em simulações de SNe Ia que não depende de parâmetros ajustáveis. O significado principal reside em demonstrar que o critério CJ impulsionado por turbulência atua como um "mecanismo de convergência". Ele impulsiona progenitores com condições iniciais vastamente diferentes (densidade, topologia de ignição) para um estado de detonação comum, oferecendo assim uma base física para a calibrabilidade empírica das SNe Ia.

Os autores concluem que a uniformidade espectroscópica de SNe Ia near-$M_{Ch}$ não é uma coincidência que requer ajuste fino, mas uma consequência mecânica do gatilho controlado por turbulência. No entanto, o artigo permanece modesto quanto às implicações mais amplas:

• Não alega definitivamente que todas as SNe Ia surgem desse canal específico, observando que mais trabalho é necessário para entender como esse mecanismo pode produzir detonações mais atrasadas ou falhar (potencialmente resultando em SNe Iax).
• Reconhece explicitamente que as simulações são puramente hidrodinâmicas. O papel dos campos magnéticos, que poderiam suprimir a cascata turbulenta ou alterar o limiar de CJ, permanece uma questão em aberto e uma área crítica para futuras simulações de Magnetohidrodinâmica (MHD).
• Os autores alertam contra atribuir imediatamente eventos semelhantes a 91T/99aa exclusivamente a progenitores near-$M_{Ch}$ sem levar em conta potenciais efeitos de campos magnéticos.

Em resumo, este trabalho estabelece uma estrutura de primeiros princípios para DDT em SNe Ia, conectando com sucesso a turbulência hidrodinâmica a características espectrais observáveis sem ajuste empírico, enquanto identifica os campos magnéticos como a próxima grande fronteira para compreender a diversidade completa dos resultados de SNe Ia.