Black Hole-Neutron Star Binaries near Neutron Star Disruption Limit in the Mass Regime of Event GW230529

Resumo Técnico: Binários de Buraco Negro-Estrela de Neutrons Próximos ao Limite de Disrupção da Estrela de Neutrons na Região de Massa do Evento GW230529

Definição do Problema
A detecção do evento de ondas gravitacionais GW230529 pela Colaboração LIGO–Virgo–KAGRA (LVK) marcou o primeiro provável de cisalhamento (merger) de buraco negro–estrela de neutrones (BHNS) com um potencial significativo de produzir contrapartidas eletromagnéticas (EM). Este evento envolve um buraco negro primário de 2,5–4,5 $M_\odot$ e um objeto compacto secundário de 1,2–2,0 $M_\odot$, colocando-o em uma região de massa onde a estrela de neutrones (NS) pode sofrer disrupção de maré em vez de mergulhar diretamente no buraco negro (BH). A disrupção de maré é um pré-requisito para gerar sinais EM observáveis, como kilonovas e surtos gama de binários compactos (cbGRBs), que requerem a formação de um disco de acreção e a ejeção de matéria rica em nêutrons.

No entanto, prever se um cisalhamento BHNS específico resultará em disrupção depende fortemente de fórmulas de ajuste analítico, como o modelo FNH18 (Foucart et al. 2018), que estima a massa do remanescente de disco pós-cisalhamento. Esses modelos foram calibrados principalmente em simulações com razões de massa $q \gtrsim 3$ e spins distantes do limiar crítico necessário para a disrupção. Consequentemente, há uma falta de dados de relatividade numérica (NR) de alta fidelidade para sistemas próximos ao limite de disrupção, particularmente para os parâmetros de massa e spin consistentes com o GW230529. Essa lacuna limita a capacidade de prever com precisão a prevalência de sinais EM e de restringir a equação de estado (EOS) da estrela de neutrones a partir de observações EM.

Metodologia
Os autores utilizam o Spectral Einstein Code (SpEC) para realizar simulações de alta resolução de cisalhamentos BHNS. O estudo foca em uma razão de massa de $q=3$ e uma massa de estrela de neutrones de $M_{NS} = 1,35 M_\odot$, parâmetros consistentes com o evento GW230529. Duas equações de estado distintas são empregadas para testar diferentes níveis de compacidade da estrela de neutrones:

1. SFHo: Uma EOS mais "suave" com um raio $R_{NS} \approx 11,9$ km.
2. DD2: Uma EOS mais "rígida" com um raio $R_{NS} \approx 13,2$ km.

As simulações variam a magnitude do spin do buraco negro ($\chi_{BH}$) para sondar a região imediatamente ao redor do spin crítico ($\chi_{crit}$) previsto pelo modelo FNH18 para a disrupção de maré. Para SFHo, $\chi_{crit} \approx 0,037$; para DD2, $\chi_{crit} \approx -0,23$. O estudo inclui simulações com spins tanto ligeiramente abaixo quanto acima desses limiares críticos.

Principais características metodológicas incluem:

• Transporte de Neutrinos por Monte Carlo: Isso marca a primeira aplicação de transporte de radiação de neutrinos por Monte Carlo em simulações de cisalhamento BHNS. Diferente dos modelos de vazamento (leakage models), este método rastreia trajetórias de neutrinos em espaço de posição e momento para modelar com precisão absorção, emissão e espalhamento. Isso é crucial para determinar a composição (fração eletrônica, $Y_e$) do remanescente pós-cisalhamento e dos fluxos (outflows), o que dita o potencial para nucleossíntese de processo r.
• Estudos de Resolução: Testes de convergência foram realizados para o caso SFHo em três resoluções para quantificar incertezas numéricas.
• Condições Iniciais: Os dados iniciais foram gerados usando o solver Spells, com eccentricidades orbitais reduzidas para $e \lesssim 0,009$.

Principais Resultados

1. Massa do Disco vs. Predições Analíticas:

   • EOS SFHo: Os resultados numéricos para a massa total do remanescente geralmente caem dentro dos erros de modelagem esperados da fórmula FNH18. No entanto, as simulações mostram consistentemente massas de remanescente ligeiramente superiores às que o modelo prevê.
   • EOS DD2: Uma discrepância significativa é observada. Para spins próximos ao limite de disrupção, o modelo FNH18 subestima a massa do disco remanescente por um fator que excede as margens de erro declaradas pelo modelo. Por exemplo, em $\chi_{BH} = -0,030$, o modelo prevê uma massa remanescente de $\sim 0,012 M_\odot$, enquanto a simulação produz $\sim 0,049 M_\odot$.
   • Limiar de Spin Crítico: Apesar das discrepâncias de massa, o valor do spin crítico necessário para desencadear a disrupção (onde a massa do remanescente transita de zero para não-zero) é encontrado como consistente entre o modelo analítico e as simulações. A principal diferença reside na inclinação da curva de massa do remanescente ($dM_{rem}/d\chi_{BH}$) próximo ao limiar; as simulações mostram um aumento mais íngreme do que o modelo.

2. Emergência de Discos de Baixa Massa:

   • As simulações revelam a formação de discos de acreção não-nulos mesmo para spins ligeiramente abaixo do limite de disrupção previsto (por exemplo, o caso SFHo com $\chi_{BH} = -0,0130$). Embora esses discos sejam de baixa massa ($\sim 10^{-3} M_\odot$), eles são potencialmente suficientes para alimentar curtos cbGRBs, conforme sugerido por estudos anteriores (ex: Gottlieb et al.).

3. Ejeção de Massa e Composição:

   • Ejecta Dinâmica: A quantidade de massa ejetada dinamicamente durante o cisalhamento é muito baixa, tipicamente $\lesssim 0,004 M_\odot$. Para as configurações de spin mais baixas, a massa não-ligada está próxima dos limites de resolução numérica, com incertezas que podem atingir $\sim 100\%$.
   • Ventos de Disco: O estudo sugere que, para esses sistemas, a massa total ejetada via ventos de disco de tempo posterior (viscosos e magneticamente impulsionados) provavelmente dominará sobre o ejecta dinâmico inicial.
   • Fração Eletrônica ($Y_e$): O ejecta dinâmico é altamente rico em nêutrons ($Y_e < 0,1$). No entanto, os fluxos (outflows) e o material do disco mostram uma gama mais ampla de $Y_e$ (até $\sim 0,3$) devido a interações fracas (aquecimento por neutrinos). A EOS DD2 produz um ejecta mais rico em nêutrons comparado à SFHo.

4. Implicações Observacionais:

   • Qualquer sinal de kilonova proveniente desses cisalhamentos específicos do tipo GW230529 é previsto como tênue, provavelmente dominado pela emissão óptica dos ventos de disco, em vez da emissão infravermelha do ejecta dinâmico.
   • A massa total de ejecta disponível para a nucleossíntese de processo r é $\lesssim 0,02 M_\odot$.

Significância e Alegações
O artigo alega que, embora o modelo analítico FNH18 permaneça amplamente confiável para prever a existência de um limiar de disrupção (o spin crítico), ele subestima sistematicamente a massa do respectivo disco remanescente, particularmente para equações de estado mais rígidas (DD2) perto do limite de disrupção.

Os autores enfatizam que seus resultados destacam os limites da calibração atual quando aplicada a cisalhamentos BHNS próximos ao limiar de disrupção. Especificamente:

• A suposição do modelo de uma massa remanescente "zero" abaixo do spin crítico não é estritamente precisa; discos de baixa massa podem se formar mesmo abaixo do limite previsto.
• Esses discos de baixa massa são fisicamente significativos como potenciais progenitores para surtos gama curtos, mesmo que não produzam kilonovas brilhantes.
• Os achados sugerem que uma recalibração modesta da fórmula de ajuste pode ser necessária à medida que mais simulações na região de baixa razão de massa e spin próximo ao crítico se tornem disponíveis.

O estudo não pretende alterar definitivamente a avaliação de probabilidade da LVK para o GW230529 ser um cisalhamento com disrupção, observando que as diferenças observadas dificilmente mudarão qualitativamente os resultados da análise da LVK. No entanto, ele fornece uma base numérica mais rigorosa para entender o comportamento de binários BHNS neste espaço de parâmetros especificamente crítico para observação.