Celephais: efficient spectral initial data code for precessing compact binaries

Resumo Técnico: Celephaïs – Dados Iniciais Espectrais Eficientes para Binárias Compactas Precessantes

Definição do Problema
Pesquisas de larga escala em relatividade numérica (NR) exigem a geração de dados iniciais de binárias compactas que sejam tanto espectralmente precisos quanto computacionalmente econômicos. Os desafios atuais incluem a construção de dados para sistemas com massas desiguais e orientações de spin arbitrárias (incluindo spins desalinhados e precessantes) sem impor simetria equatorial. Os métodos espectrais globais tradicionais enfrentam uma tensão entre robustez e custo: a fatoração direta do Jacobiano multidomínio globalmente acoplado exige armazenamento de $O(N_u^2)$, o que se torna proibitivo para sistemas de alta resolução ou de alta razão de massa. O refinamento de grade uniforme introduz graus de liberdade (DOFs) desnecessários em regiões onde as caudas espectrais já estão bem resolvidas, particularmente em binárias de alta razão de massa, onde os requisitos de resolução radial diferem significativamente entre os domínios.

Metodologia
O artigo apresenta o Celephaïs, um código de dados iniciais de objetos compactos construído sobre a biblioteca espectral Kadath. O código resolve as equações da Sanduíche Conformal Estendida (XCTS) para sistemas de estrelas de nêutrons binárias (BNS) e de buraco negro–estrela de nêutron (BHNS). Os principais componentes metodológicos incluem:

1. Construção de Jacobiano Esparso: Ao contrário de abordagens anteriores que dependem de métodos Newton–Krylov sem matriz (matrix-free) com precondicionadores prescritos, o Celephaïs constrói explicitamente o Jacobiano esparso para seu sistema elíptico de forma forte. Ele utiliza um grafo de dependência para identificar zeros estruturais antes da diferenciação, empregando diferenciação automática de modo direto em lote (batched forward-mode) para avaliar apenas as entradas não nulas.
2. Precondicionamento por Separação de Paridade: O código explora a separação aproximada de paridade dos campos sob simetria de reflexão ($y \to -y$). Isso permite que o Jacobiano seja diagonalizado em blocos para os setores par e ímpar. A matriz esparsa é fatorada usando a biblioteca MUMPS.
3. Newton–Krylov sem Matriz com Precondicionamento Atualizado: Os fatores LU do Jacobiano inicial são reutilizados como um precondicionador à direita para as subsequentes iterações Newton–Krylov sem matriz. Esta abordagem captura os acoplamentos complexos do primeiro passo de Newton enquanto mantém a eficiência de memória das soluções sem matriz. O precondicionador é periodicamente atualizado ou remontado se o resíduo crescer bruscamente.
4. Refinamento Adaptativo $hp$: Para lidar com a ineficiência de grades uniformes, o código implementa um esquema de refinamento $hp$ anisotrópico. Ele avalia independentemente as caudas espectrais de cada domínio e direção de coordenada.
   • Movimentos $p$: Aumentam a ordem espectral em direções específicas.
   • Movimentos $h$: Subdividem as cascas radiais.
     O esquema visa as caudas espectrais não resolvidas, concentrando a resolução onde necessário (por exemplo, na direção radial em binárias de alta razão de massa) enquanto mantém a resolução angular fixa quando suficiente.
5. Redução de Excentricidade: O código estende as estimativas de balanço de força pós-newtonianas (PN) para orientações de spin arbitrárias para inicializar a frequência orbital ($\Omega$) e a velocidade de aproximação radial ($\dot{a}$). Estas são refinadas via evoluções curtas para reduzir a excentricidade orbital.

Principais Contribuições e Resultados

• Eficiência e Escalabilidade: Para uma BNS precessante representativa (Tabela I) com uma estrela sem spin e outra com spin, a primeira montagem do Jacobiano leva menos de 10 minutos, e o trabalho completo usa menos de 17 GB de RAM em um MacBook Pro com um processador Apple M4 Max. O método resolve sistemas com sucesso sem requisitos de armazenamento denso.
• Desempenho do Refinamento Adaptativo:
  • Para um benchmark BHNS de razão de massa-20, o esquema $hp$ adaptativo recupera a precisão de restrição de uma grade uniforme usando aproximadamente três vezes menos incógnitas.
  • Para um BHNS de razão de massa-10, o refinamento $p$ anisotrópico atinge a precisão da grade uniforme ($N=15$) com cerca de 1,6 vezes menos incógnitas.
• Validação:
  • Compatibilidade: As soluções geradas pelo código anterior FUKA (para spins alinhados) permanecem soluções válidas no Celephaïs.
  • Consistência PN: As sequências de energia de ligação para sistemas BNS e BHNS concordam com as estimativas pós-newtonianas (incluindo termos de spin-órbita, spin-spin e quadrupolo induzido por spin) nas resoluções de produção.
  • Redução de Excentricidade: Uma única correção baseada em evolução, inicializada pelas estimativas estendidas de PN, reduz a excentricidade medida de uma BNS precessante para abaixo de $10^{-2}$, removendo oscilações visíveis da frequência orbital.
  • Evolução Completa: Uma BNS precessante foi evoluída de forma estável através do inspiral e da fusão usando o código SACRA. A forma de onda extraída permitiu a reconstrução do eixo de precessão.
  • Reconstrução do Eixo de Precessão: O eixo de radiação reconstruído da forma de onda segue de perto um modelo de precessão simples de spin único (ordem 3.5PN) durante a maior parte do inspiral, validando a geometria de spin inicial.

Significância e Alegações
O autor afirma que o Celephaïs estabelece uma rota eficiente para gerar dados iniciais espectralmente resolvidos para sistemas BNS e BHNS com orientações de spin arbitrárias. O código demonstra que:

1. A construção explícita de Jacobiano esparso combinada com a fatoração MUMPS e o precondicionamento atualizado é uma estratégia viável para resolver sistemas elípticos complexos e globalmente acoplados em hardware padrão (ex: laptops).
2. O refinamento adaptativo $hp$ reduz significativamente o custo computacional (em termos de incógnitas) para binárias de alta razão de massa sem sacrificar a precisão das restrições.
3. As estimativas estendidas de balanço de força PN fornecem um ponto de partida robusto para a redução de excentricidade em sistemas precessantes, exigindo apenas uma calibração dinâmica mínima.

O artigo observa limitações, incluindo o fato de que a evolução atual não estabelece a convergência de fase da forma de onda (exigindo estudos de múltiplas resoluções futuros) e que o buraco negro é tratado via excisão em vez de uma formulação de puncture. Trabalhos futuros estão planejados para incluir formulações de puncture e extensões para gravidade escalar-tensor e escalar-Gauss-Bonnet.