Convective stability analysis of massive neutron stars formed in binary mergers

Resumo Técnico: Análise de Estabilidade Convectiva de Estrelas de Nêutrons Massivas Formadas em Fusões Binárias

Definição do Problema
A evolução pós-fusão de fusões de estrelas de nêutrons binárias (BNS) é governada pela dinâmica da estrela de nêutrons (NS) massiva resultante. Um aspecto crítico dessa evolução é o potencial desenvolvimento de instabilidades hidrodinâmicas, particularmente a instabilidade convectiva, que pode impulsionar movimentos fluídicos de grande escala, alterar a estrutura térmica e excitar modos de ondas gravitacionais (GW). Estudos anteriores produziram conclusões conflitantes: algumas simulações utilizando o critério de Schwarzschild sugeriram que NSs massivas desenvolvem instabilidade convectiva de grande escala, excitando subsequentemente modos inerciais observáveis [74, 75]. Inversamente, outros estudos incorporando radiação de neutrinos encontraram estratificação estável [55]. Uma limitação fundamental em análises anteriores é a frequente negligência das forças restauradoras de rotação. Dado que os remanescentes de fusão exibem rotação rápida e diferencial, a contribuição rotacional para a estabilidade (frequentemente descrita pelo critério de Solberg-Høiland) é essencial. Este artigo visa resolver essas discrepâncias realizando simulações de hidrodinâmica totalmente relativísticas gerais e aplicando um critério de estabilidade abrangente que considera tanto a flutuabilidade térmica quanto os efeitos rotacionais.

Metodologia
Os autores utilizam o código de relatividade numérica SACRA-MPI para simular a fusão de sistemas BNS de massas iguais durante uma duração de aproximadamente 100 ms pós-fusão.

• Equação de Estado (EOS): Três EOS de politropos segmentados (APR4, SLy4, MPA1) são utilizados, aumentados com um componente térmico de lei-$\Gamma$ ($\Gamma_{\text{th}} = 1.8$) para modelar o aquecimento por choque.
• Configuração Numérica: As simulações utilizam uma grade cartesiana aninhada com refinamento de 2:1 e simetria de espelho através do plano orbital. Três resoluções são testadas para o modelo APR4 para garantir a convergência.
• Análise de Estabilidade: Os autores derivam e aplicam os critérios de estabilidade de Solberg-Høiland relativísticos (Critério I e Critério II) para estrelas quentes e diferencialmente rotativas. Estes critérios incorporam a interação entre a flutuabilidade térmica (discriminante de Schwarzschild) e as forças restauradoras de rotação (gradientes de momento angular específico). A análise é realizada nos perfis médios angulares de entropia específica e momento angular específico.
• Análise de Modos: O estudo analisa a excitação de modos fluídicos e a emissão de GW, focando especificamente no modo de um braço $m=1$ e na potencial excitação de modos inerciais $m=2$ após o amortecimento do modo quadrupolar $f$.

Principais Contribuições

1. Derivação de Critérios de Estabilidade Relativísticos: O artigo fornece uma aplicação rigorosa dos critérios de Solberg-Høiland relativísticos (Eqs. 46 e 47) para NSs pós-fusão massivas, contabilizando explicitamente a força restauradora de rotação que foi negligenciada em estudos que dependiam apenas do critério de Schwarzschild.
2. Estrutura de Simulação Independente: Os autores conduzem uma análise independente usando configurações numéricas (escolhas de EOS, resoluções, esquemas hidrodinâmicos) distintas de estudos conflitantes anteriores, evitando comparações diretas um-para-um para garantir uma investigação autossuficiente da física.
3. Testes de Simetria e Conservação de Momento: O estudo introduz um teste de preservação de simetria para verificar se o código numérico não quebra artificialmente a simetria $\pi$. Além disso, investiga a correlação entre o crescimento da amplitude do modo $m=1$ e a violação da conservação do momento linear, um teste crítico para distinguir instabilidades físicas de artefatos numéricos.

Resultados

• Estabilidade Convectiva: As simulações não mostram evidências de instabilidade convectiva de grande escala. Tanto a entropia específica média angular quanto o momento angular específico aumentam para fora dentro das NSs massivas. Embora regiões locais possam parecer instáveis sob o critério de Schwarzschild isoladamente, a inclusão de efeitos rotacionais (especificamente o aumento externo do momento angular específico) estabiliza essas regiões. Os critérios de estabilidade global (Critério I e II) são satisfeitos para a vasta maioria da estrela.
• Modos Inerciais: Contrariamente aos achados das Refs. [74, 75], nenhuma excitação de modos inerciais observáveis é detectada após o amortecimento dos modos $f$ quadrupolares. O sinal de GW permanece dominado pelo modo $f$ durante toda a simulação, sem a emergência de modos inerciais discretos em frequências de quilohertz.
• Instabilidade $m=1$: Observa-se que um modo de um braço $m=1$ cresce e satura, consistente com a literatura anterior. No entanto, os autores encontram uma forte correlação entre o crescimento da amplitude deste modo e a violação da conservação do momento linear. Isso sugere que a excitação do modo de um braço nessas simulações específicas pode ter uma origem numérica, em vez de ser uma instabilidade puramente física.
• Perfis Térmicos e Rotacionais: A NS massiva estabelece um estado de quase-equilíbrio dentro de 20–30 ms. O núcleo é frio e rotaciona lentamente, enquanto o envelope é aquecido por interações de braços espirais e rotaciona diferencialmente. O momento angular específico aumenta para fora, fornecendo uma forte força centrífuga estabilizadora.

Significância e Alegações
O artigo alega que a inclusão de forças restauradoras de rotação altera fundamentalmente a avaliação de estabilidade de NSs pós-fusão massivas. Os autores afirmam que seus resultados, que mostram uma configuração estável, são incompatíveis com estudos anteriores que previam convecção de grande escala e subsequente excitação de modos inerciais. Eles argumentam que a ausência de instabilidade convectiva em suas simulações explica a falta de modos inerciais excitados.

Em relação ao modo $m=1$, os autores mantêm uma postura modesta, afirmando que "não podem rejeitar a possibilidade de que a excitação do modo de um braço tenha uma origem numérica" devido à sua correlação com a não conservação do momento linear. Eles enfatizam que, embora o código passe em testes de preservação de simetria, a conservação do momento linear em simulações de malha refinada continua sendo um desafio que requer investigação adicional para estabelecer definitivamente a natureza física da instabilidade de um braço.

O estudo conclui que pesquisas futuras incorporando microfísica (gradientes de composição) e campos magnéticos são necessárias para fornecer uma conclusão definitiva, mas dentro do escopo de suas atuais simulações hidrodinâmicas, as NSs massivas formadas em fusões BNS parecem ser convectivamente estáveis.