Lightcurve Modelling of 2,205 ZTF DR2 Type~Ia Supernovae: Implications for SN Ia Physics and Cosmology

Resumo Técnico: Modelagem de Curvas de Luz de 2.205 SNe Ia do ZTF DR2

Problema e Motivação
As supernovas do tipo Ia (SNe Ia) servem como velas padronizáveis para a cosmologia, baseando-se em relações empíricas como a relação "mais brilhante-mais lenta" (Phillips) e a relação "mais brilhante-mais azul". No entanto, as origens físicas dessas relações e a diversidade dos mecanismos de explosão (por exemplo, progenitores de massa de Chandrasekhar vs. massa sub-Chandrasekhar) permanecem incompletamente compreendidas. Estudos sistemáticos anteriores de parâmetros físicos (ex: massa de ejecta $M_{\rm ej}$, massa de níquel sintetizado $M_{\rm Ni}$) foram limitados por tamanhos de amostra pequenos (ex: Scalzo et al. 2014 analisou 337 SNe) ou basearam-se em relações empíricas em vez de ajuste direto da curva de luz. Os modelos padrão de decaimento radioativo de uma zona frequentemente falham em reproduzir os máximos secundários observados em bandas vermelhas/infravermelho próximo, que são atribuídos à física de recombinação Fe III $\to$ Fe II. Este artigo visa realizar a maior extração sistemática de parâmetros físicos de curvas de luz de SNe Ia até o momento, aproveitando o Zwicky Transient Facility (ZTF) Data Release 2 (DR2) para restringir a física da explosão, os canais de progenitores e a base física para a padronização cosmológica.

Metodologia
Os autores analisaram 2.205 SNe Ia espectroscopicamente confirmadas e de qualidade cosmológica da amostra ZTF DR2. A análise empregou um framework de inferência Bayesiana usando um modelo de decaimento radioativo de uma zona modificado. As principais características metodológicas incluem:

• Opacidade Dependente do Comprimento de Onda: O modelo padrão de opacidade cinza foi estendido fenomenologicamente para incluir a física de recombinação Fe III $\to$ Fe II. A opacidade óptica $\kappa$ foi tornada uma função do comprimento de onda ($\lambda$) e do tempo ($t$), permitindo que o modelo reproduza os máximos secundários em bandas vermelhas ($i, z$) ao reduzir a opacidade em comprimentos de onda mais longos, aumentando assim a temperatura e o fluxo do fotosfera.
• Parâmetros Físicos: O modelo restringiu simultaneamente 10 parâmetros livres por evento, incluindo massa de ejecta ($M_{\rm ej}$), fração de níquel ($f_{\rm Ni} = M_{\rm Ni}/M_{\rm ej}$), velocidade de ejecta ($v_{\rm ej}$), tempo de explosão ($t_0$) e parâmetros de opacidade (opacidade base, época de recombinação, largura e força). A extinção do hospedeiro ($A_V$) também foi modelada.
• Framework de Inferência: Uma verossimilhança de mistura Gaussiana (mixture Gaussian likelihood) foi usada para lidar com outliers e dados dispersos. Os autores utilizaram o amostrador nested sampler PyMultiNest via interface Bilby para amostrar as distribuições posteriores para todos os 2.205 eventos.
• Modelagem Hierárquica: Para abordar os vieses de seleção inerentes à amostra de qualidade cosmológica (que exclui eventos subluminosos como os do tipo 91bg), uma análise Bayesiana hierárquica foi realizada em uma subamostra limitada em redshift ($z \le 0,06$, $N=902$). Esta abordagem modelou as distribuições de nível populacional de $M_{\rm ej}$ e $M_{\rm Ni}$ como Gaussianas truncadas enquanto marginalizava sobre as incertezas individuais dos eventos e uma função de seleção dependente do brilho.

Principais Resultados

1. Distribuições de Parâmetros Físicos:

   • Para a amostra completa, a mediana da massa de ejecta é $M_{\rm ej} \approx 1,40 M_\odot$ e a mediana da massa de níquel é $M_{\rm Ni} \approx 1,19 M_\odot$. O alto $M_{\rm Ni}$ é atribuído ao viés de seleção que favorece eventos luminosos.
   • A distribuição da massa de ejecta é amplamente unimodal, com pico próximo à massa de Chandrasekhar, sem uma bimodalidade clara que sugira canais de progenitores distintos.
   • Subtipos: Eventos do tipo SN 1991T são fisicamente distintos, exibindo massas de ejecta sistematicamente mais altas (mediana de $1,64 M_\odot$ vs. $1,38 M_\odot$ para normais) e produzindo $\sim30\%$ mais níquel. 84% dos eventos do tipo 1991T possuem $M_{\rm ej} > 1,4 M_\odot$.

2. Correlações com Parâmetros Empíricos:

   • Uma forte correlação positiva ($\rho = 0,57$) foi encontrada entre a massa de níquel sintetizado ($M_{\rm Ni}$) e o parâmetro de stretch $x_1$ do SALT2. Isso fornece uma base física para a relação "mais brilhante-mais lenta": SNe mais brilhantes sintetizam mais níquel e, dentro deste modelo, possuem massas de ejecta maiores, levando a tempos de difusão de fótons mais longos.
   • O stretch do SALT2 ($x_1$) correlaciona-se com a massa de ejecta ($M_{\rm ej}$), enquanto a cor do SALT2 ($c$) não mostra correlação com $M_{\rm ej}$, sugerindo que a cor é dominada pela extinção de poeira e efeitos de recombinação secundária, em vez de propriedades de explosão globais.

3. Dependências da Galáxia Hospedeira:

   • SNe Ia em hospedeiros de baixa massa ($\log_{10}(M_*/M_\odot) < 10$) sintetizam $\sim12\%$ mais níquel ($\Delta M_{\rm Ni} = 0,13 M_\odot$) do que aquelas em hospedeiros de alta massa.
   • Esta diferença intrínseca na produção de níquel oferece uma explicação física para o observado "degrau de massa" (mass step) nos resíduos de Hubble. Os autores sugerem que a padronização baseada em parâmetros físicos (como $M_{\rm Ni}$) poderia potencialmente remover este degrau, enquanto a padronização empírica (SALT2) deixa um offset residual.

4. Inferência Populacional Hierárquica:

   • Corrigindo os vieses de seleção na subamostra de baixo redshift ($z \le 0,06$), os autores inferiram as verdadeiras distribuições populacionais.
   • A população é bem descrita por distribuições Gaussianas com $\mu_{\rm ej} = 1,26 \pm 0,01 M_\odot$ e $\mu_{\rm Ni} = 0,64 \pm 0,06 M_\odot$.
   • As frações inferidas são: $43 \pm 2\%$ sub-Chandrasekhar ($M_{\rm ej} < 1,2 M_\odot$), $34 \pm 1\%$ near-Chandrasekhar ($1,2 \le M_{\rm ej} \le 1,5 M_\odot$) e $24 \pm 2\%$ super-Chandrasekhar ($M_{\rm ej} > 1,5 M_\odot$).
   • A distribuição contínua e suave e a divisão semelhante entre eventos sub- e near-Chandrasekhar sugerem um único mecanismo de explosão dominante (provavelmente Detonação Retardada) operando através de uma gama de massas de progenitores, em vez de canais distintos e bimodais.

Significância e Alegações
O artigo afirma fornecer a maior caracterização física homogênea da população local de SNe Ia até o momento. Suas principais contribuições são:

• Base Física para Padronização: Estabelece um link direto entre a relação empírica de Phillips e a física fundamental da explosão (massa de níquel e massa de ejecta), confirmando que a relação "mais brilhante-mais lenta" é uma manifestação da física de difusão.
• Origem do Degrau de Massa: Conecta o degrau de massa da galáxia hospedeira às diferenças intrínsecas na eficiência de produção de níquel, provavelmente impulsionadas pela idade ou metalicidade do progenitor, em vez de puramente de sistemáticas observacionais.
• Restrições de Progenitor: Ao analisar uma amostra duas ordens de magnitude maior do que os esforços anteriores de modelagem física, o estudo sugere que a população de SNe Ia cosmologicamente útil é dominada por explosões de massa próxima a Chandrasekhar com uma distribuição contínua de massas de progenitores, desafiando a noção de canais de explosão distintos e separados para diferentes subtipos.
• Viabilidade Metodológica: O trabalho demonstra que modelos fisicamente motivados e dependentes do comprimento de onda podem ser aplicados a milhares de curvas de luz com eficiência computacional comparável ao ajuste empírico, pavimentando o caminho para a padronização física em futuros levantamentos de grande escala (ex: Rubin, Roman).

Os autores permanecem modestos quanto às incertezas sistemáticas, reconhecendo que seu modelo de uma zona e tratamentos de opacidade simplificados podem introduzir vieses. Eles enfatizam que seus resultados caracterizam a população "cosmologicamente útil" (tendenciosa para eventos luminosos) em vez de toda a diversidade de transientes térmonucleares, e que a modelagem hierárquica é essencial para recuperar as propriedades reais da população.