Decoding the Early-Time Light Curves of Type Ia Supernovae. II. Population Parameters of One Thousand ZTF Supernovae

Resumo Técnico: Decodificando as Curvas de Luz de Tempo Inicial de Supernovas do Tipo Ia. II. Parâmetros Populacionais de Mil Supernovas ZTF

Definição do Problema
As supernovas do tipo Ia (SNe Ia) são indicadores de distância cosmológica padronizáveis, embora seus mecanismos de explosão e sistemas progenitores subjacentes permaneçam debatidos. Enquanto as curvas de luz de tempo tardio próximas ao pico de brilho são bem caracterizadas, a emissão de tempo inicial (dentro da primeira semana) codifica informações críticas sobre o ejecta externo, a distribuição de $^{56}\text{Ni}$ radioativo e potenciais interações com material circunestelar (CSM) ou companheiras. Historicamente, estudos de SNe Ia de tempo inicial foram limitados por amostras pequenas, vieses de seleção (ex: viés de Malmquist) e pela agregação de parâmetros ajustados individualmente, o que pode introduzir correlações espúrias e sistemáticas. Além disso, a relação entre a morfologia de ascensão de tempo inicial e o parâmetro de estiramento (stretch) global da curva de luz ($x_1$) não foi rigorosamente quantificada através de uma população estatisticamente completa. Este artigo aborda a necessidade de uma caracterização robusta, em nível populacional, de SNe Ia de tempo inicial usando uma amostra volumétrica completa para quebrar degenerações entre canais de explosão e propriedades dos progenitores.

Metodologia
Os autores aplicam um arcabouço bayesiano hierárquico, previamente desenvolvido no "Paper I" desta série, a uma amostra volumétrica completa de 972 SNe Ia espectroscópicamente normais do Zwicky Transient Facility Data Release 2 (ZTF DR2, $z < 0.06$). Esta amostra é uma ordem de magnitude maior do que os conjuntos de dados anteriores usados para análises semelhantes.

Componentes metodológicos principais incluem:

• Modelo de Ascensão de Lei de Potência: As curvas de luz de tempo inicial são modeladas usando uma ascensão de lei de potência empírica ($F \propto (t - t_{\text{fl}})^\alpha$) até um limiar de truncamento (30%, 40% ou 50% do fluxo de pico). O modelo inclui parâmetros para o tempo de primeira luz ($t_{\text{fl}}$), amplitude e um índice de ascensão específico do filtro ($\alpha$).
• Inferência Hierárquica: Em vez de ajustar SNe individuais e agregar resultados, o arcabouço infere simultaneamente parâmetros individuais e hiperparâmetros de nível populacional (médias, dispersões e covariâncias). Esta abordagem naturalmente reduz o peso de eventos com amostragem esparsa ou ruidosa, mitigando vieses de seleção e o viés de Malmquist sem exigir cortes de dados rigorosos.
• Tratamento de Incertezas: A função de verossimilhança incorpora um fator de escala ($\beta_d$) para corrigir incertezas fotométricas subestimadas no pipeline do ZTF e propaga incertezas do tempo de pico ($t_0$) inferido a partir de ajustes de curva de luz completos (SALT2).
• Análise de Correlação: O estudo investiga correlações entre parâmetros de tempo inicial ($t_{\text{rise}}$, $\alpha$, evolução de cor $\alpha_g - \alpha_r$) e o parâmetro de estiramento SALT2 ($x_1$) usando um modelo de regressão de erros-nas-variáveises bayesiano com um ajuste linear quebrado para identificar quebras estruturais.
• Comparação Teórica: As correlações observadas são comparadas contra curvas de luz sintéticas de dois arcabouços teóricos: modelos 1D de massa próxima a Chandrasekhar (perto de-$M_{\text{Ch}}$) com perfis de densidade e mistura de $^{56}\text{Ni}$ variáveis (código TURTLS) e modelos 2D de dupla detonação de sub-Chandrasekhar ($M_{\text{Ch}}$) com variações de ângulo de visão.

Resultados Principais

1. Parâmetros Populacionais: O estudo fornece restrições robustas sobre a população de tempo inicial:

   • Tempo de ascensão médio ($t_{\text{rise}}$): $18,55 \pm 0,08$ dias.
   • Índice de ascensão médio na banda $r$ ($\alpha_r$): $2,10 \pm 0,04$ (mais íngreme que o modelo de fireball $\alpha=2$).
   • Evolução de cor média ($\alpha_g - \alpha_r$): $0,20 \pm 0,02$, indicando que o fluxo de tempo inicial evolui mais estritamente na banda $g$ do que na banda $r$ (evolução mais azul).
   • Existe dispersão intrínseca significativa em todos os parâmetros ($\sigma_{t_{\text{rise}}} \approx 1,42$ dias, $\sigma_{\alpha_r} \approx 0,48$).

2. Bifurcação com Estiramento ($x_1$): A relação entre a morfologia de tempo inicial e o estiramento da curva de luz exibe uma quebra estrutural distinta perto de $x_1 \approx 0,2\text{--}0,5$:

   • População de alto estiramento ($x_1 > x_{\text{brk}}$): Mostra tendências claras onde um $x_1$ mais alto correlaciona-se com ascensões mais rasas (menor $\alpha$), cores azuis mais persistentes (menor $\alpha_g - \alpha_r$) e um achatamento da correlação $t_{\text{rise}}$ vs. $x_1$.
   • População de baixo estiramento ($x_1 < x_{\text{brk}}$): Carece de correlações significativas entre parâmetros de tempo inicial e $x_1$.
   • Um "tempo de ascensão ajustado" (medido a partir de um epoch de fração de fluxo em vez de extrapolado para a primeira luz) revela uma correlação contínua e forte com $x_1$, sem quebra, sugerindo que a quebra no $t_{\text{rise}}$ padrão é impulsionada por variações na evolução do fluxo inicial ($t_{X\%}$).

3. Anomalias e Subclasses Peculiares:

   • Eventos tipo 03fg: Exibem tempos de ascensão significativamente mais longos e índices de ascensão mais íngremes, consistentes com excessos de fluxo de curta duração de interações ejecta-CSM em fusões violentas de duplas anãs brancas.
   • Eventos tipo 02cx e 91bg: Mostram tempos de ascensão mais curtos e índices mais rasos, consistentes com deflagrações parciais e baixas massas de ejecta.
   • Anomalias de SNe Ia Normais: Os autores identificam SNe Ia normais com tempos de ascensão incomumente longos (potencialmente excessos de curta duração) e um subconjunto de eventos tipo 91T de alto estiramento com ascensões quase lineares ($\alpha \approx 1$). Estes "ascensores lineares" parecem ocupar a cauda da distribuição de eventos normais luminosos e de evolução lenta, em vez de formar uma subclasse distinta, sugerindo uma mistura outward de $^{56}\text{Ni}$ generalizada na população de alto estiramento.

4. Comparação de Modelos:

   • Modelos Próximos a $M_{\text{Ch}}$: Embora capazes de produzir diversas curvas de luz, modelos 1D com variações de densidade e mistura de $^{56}\text{Ni}$ falham em reproduzir as correlações quebradas entre parâmetros de tempo inicial e estiramento.
   • Duplas Detonações de Sub-$M_{\text{Ch}}$: Modelos 2D com ejecta assimétricos e efeitos de ângulo de visão reproduzem com sucesso uma correlação quebrada na evolução de cor ($\alpha_g - \alpha_r$), mas falham em corresponder aos tempos de ascensão observados (prevendo ascensões mais curtas) e não reproduzem a quebra na relação do índice de ascensão.

Significância e Alegações
O artigo afirma fornecer a primeira caracterização abrangente, em nível populacional, de SNe Ia de tempo inicial usando uma amostra volumétrica completa, indo além das limitações de amostras pequenas e enviesadas. Ao empregar um arcabouço bayesiano hierárquico, os autores demonstram que cortes de qualidade estritos podem introduzir vieses de seleção (ex: enviesando amostras para eventos de alto estiramento) e que inferências populacionais robustas podem ser feitas mesmo com dados de amostragem esparsa.

A principal significância reside na descoberta de uma dicotomia no comportamento de tempo inicial das SNe Ia baseada no estiramento da curva de luz. A quebra observada sugere que os mecanismos físicos que impulsionam a emissão inicial diferem fundamentalmente entre as populações de baixo e alto estiramento. Especificamente, a prevalência de ascensões rasas e cores azuis persistentes na população de alto estiramento aponta para uma mistura outward generalizada de $^{56}\text{Ni}$ nas camadas externas do ejecta.

Os autores concluem que os modelos teóricos atuais (tanto 1D próximo a $M_{\text{Ch}}$ quanto 2D sub-$M_{\text{Ch}}$) são insuficientes para reproduzir totalmente o espaço de parâmetros observado, particularmente as correlações específicas com o estiramento. Eles argumentam que modelos de explosão multidimensionais com configurações de progenitores mais realistas e tratamentos da distribuição de $^{56}\text{Ni}$ e assimetria do ejecta são necessários para explicar a diversidade observada. O artigo posiciona estas restrições empíricas como um marco crítico para desenvolvimentos teóricos futuros e destaca o potencial de próximos grandes levantamentos (ex: Observatório Rubin) para refinar ainda mais estas demografias populacionais.