Inferring Neutron Star Nuclear Properties from Gravitational-Wave and Gamma-Ray Burst Observations

Resumo Técnico: Inferência de Propriedades Nucleares de Estrelas de Nêutrons a partir de Observações de Ondas Gravitacionais e Explosões de Raios Gama

Enunciado do Problema
Observações recentes de explosões de raios gama de longa duração (lGRBs) acompanhadas de emissão de kilonova (ex: GRB 211211A, GRB 230307A) desafiam o paradigma tradicional de que os lGRBs surgem exclusivamente do colapso de estrelas massivas. Se esses eventos se originam de fusões de binárias de estrelas de nêutrons (BNS) ou de estrela de nêutrons–buraco negro (NSBH), uma questão crítica surge: qual propriedade física governa a duração da explosão? Modelos teóricos sugerem que a massa total da fusão determina o tempo de vida do remanescente da fusão. Fusões mais leves podem formar remanescentes de estrelas de nêutrons (NS) de longa duração que alimentam sGRBs curtos, enquanto fusões mais pesadas resultam em remanescentes de curta duração que colapsam em buracos negros (BHs), potencialmente alimentando lGRBs via discos de acreção massivos. No entanto, os limiares de massa específicos que separam esses resultados e sua relação com propriedades nucleares fundamentais, especificamente a massa de Tolman-Oppenheimer-Volkoff ($M_{\text{TOV}}$), permanecem observacionalmente não restringidos.

Metodologia
Os autores desenvolvem um arcabouço bayesiano para inferir escalas de massa características comparando taxas de fusão de ondas gravitacionais (GW) da colaboração LIGO-Virgo-Kagra (LVK) com as taxas observadas de sGRBs e lGRBs associados a kilonovas.

1. Definições de Massa Característica: O estudo parametriza três limiares de massa crítica em relação à massa TOV ($M_{\text{TOV}}$):

   • $M_{\text{VL}} = k_1 M_{\text{TOV}}$: A fronteira entre remanescentes de vida muito longa e de longa duração.
   • $M_{\text{LS}} = k_2 M_{\text{TOV}}$: A fronteira entre remanescentes de longa duração e de curta duração. Fusões com $M_{\text{VL}} < M_{\text{tot}} \le M_{\text{LS}}$ são hipotetizadas para produzir sGRBs.
   • $M_{\text{PC}} = k_3 M_{\text{TOV}}$: A fronteira entre remanescentes de curta duração e colapso imediato (prompt collapse). Fusões com $M_{\text{LS}} < M_{\text{tot}} \le M_{\text{PC}}$ são hipotetizadas para produzir lGRBs. Fusões com $M_{\text{tot}} > M_{\text{PC}}$ sofrem colapso imediato; lGRBs são produzidos apenas se a razão de massa ($q$) permitir a formação de um disco suficiente ($1,2 \le q \le 3$ para BNS, $q \le 3$ para NSBH).

2. Fontes de Dados:

   • Taxas de GW: O catálogo GWTC-4 da LVK (Abac et al. 2025a) e o modelo de população FullPop-4.0 (Abac et al. 2025b) fornecem taxas de fusão locais como função das massas dos componentes.
   • Taxas de GRB: As taxas intrínsecas para sGRBs ($R_{\text{sGRB}}$) e a razão de taxa $\eta = R_{\text{lGRB}}/R_{\text{sGRB}}$ são tratadas como restrições observacionais com incertezas significativas, modeladas via distribuições normais.

3. Inferência: Usando o teorema de Bayes, os autores inferem as distribuições a posteriori para o conjunto de parâmetros $\vec{\Theta} = (k_1, k_2, k_3, M_{\text{TOV}})$. A função de verossimilhança compara as taxas de GRB induzidas por GW previstas (derivadas dos limiares de massa) contra as taxas observadas de GRB. Os priors para $M_{\text{TOV}}$ e o raio da NS na massa TOV ($R_{\text{TOV}}$) são extraídos de Legred et al. (2021), condicionados a observações de pulsares e GW.

Principais Resultados

• Razão de Massa Característica ($k_2$): A análise produz um valor mediano de $k_2 = 1,36^{+0,08}_{-0,09}$ (intervalo de confiança de 68%), com $k_2 > 1,24$ com 90% de confiança. Isso indica que a transição de remanescentes de longa duração para curta duração ocorre em aproximadamente $1,36 M_{\text{TOV}}$.
• Implicação para Tempos de Vida de Remanescentes: O alto valor de $k_2$ sugere que estrelas de nêutrons massivas podem sobreviver por períodos estendidos (centenas de milissegundos) após a fusão, antes de colapsarem, em vez de colapsarem imediatamente.
• Correlações: O estudo identifica uma forte correlação entre $k_2$ e $k_3$, e uma anti-correlação entre $k_2$ e $M_{\text{TOV}}$. Esta anti-correlação permite que as restrições sobre a razão de massa característica $k_2$ sejam mapeadas diretamente para limites superiores da massa TOV.
• Robustez: A inferência de $k_2$ permanece robusta através de uma ampla gama de taxas de sGRB assumidas ($30$a$2000 \text{ Gpc}^{-3}\text{ano}^{-1}$) e razões de taxa ($\eta$ de $0,01$a$3$). Variações nos limiares de razão de massa ($q_{\text{low}}$ e $q_{\text{high}}$) e a inclusão ou exclusão da fronteira $M_{\text{VL}}$ têm impacto mínimo na inferência de $k_2$.

Significância e Alegações
Este artigo estabelece um método novo e independente para restringir a equação de estado (EoS) da estrela de nêutrons, ligando populações de fusão de GW a durações de GRB. Ao identificar uma escala de massa característica ($1,36 M_{\text{TOV}}$) que separa os progenitores de sGRB e lGRB, os autores fornecem uma restrição observacional direta sobre a massa máxima de estrelas de nêutrons.

Os autores observam que seus achados são consistentes com simulações de relatividade numérica (ex: Perna et al. 2025), que preveem uma transição entre $1,3$e$1,4 M_{\text{TOV}}$. Eles enfatizam que, embora os resultados atuais sejam robustos contra incertezas observacionais, melhorias futuras nas medições de taxas de GW e GRB, particularmente detecções conjuntas de GRBs com eventos de GW, tornarão mais estreita a correlação entre massas características e $M_{\text{TOV}}$, potencialmente reduzindo a incerteza sobre a massa TOV para $\sim 0,1 M_{\odot}$. Esta abordagem complementa as técnicas existentes e oferece um caminho para estabelecer limites mais rigorosos para a EoS de NS.