Assessing the stability of ultracompact spinning boson stars with nonlinear evolutions

Título: Avaliação da Estabilidade de Estrelas de Bósons Ultracompactas Giratórias com Evoluções Não Lineares

Autores: Tamara Evstafyeva, Nils Siemonsen e William E. East.
Instituições: Perimeter Institute for Theoretical Physics, Princeton Gravity Initiative, Princeton University.

1. O Problema e o Contexto

O artigo investiga a estabilidade dinâmica de estrelas de bósons (BSs) ultracompactas e giratórias que possuem anéis de luz (light rings - LRs) estáveis.

• Contexto Teórico: Objetos ultracompactos sem horizonte de eventos (como estrelas de bósons) são candidatos a "imitadores de buracos negros". Diferentemente dos buracos negros, que possuem apenas um anel de luz instável, esses objetos podem possuir um anel de luz interno onde geodésicas nulas são estavelmente presas.
• A Hipótese de Instabilidade: Especula-se que o aprisionamento estável de modos de massa zero (como ondas gravitacionais ou escalares) dentro desse anel de luz leve a um acúmulo de energia através de efeitos não lineares fracos, resultando em uma instabilidade não linear do tipo "anel de luz".
• Trabalho Anterior: Estudos anteriores (ex: Ref. [14]) sugeriram que estrelas de bósons ultracompactas poderiam sofrer dessa instabilidade em escalas de tempo de $t\mu \sim 5000$, levando à formação de buracos negros ou migração para modelos menos compactos.
• Objetivo: Reinvestigar essa estabilidade utilizando simulações de relatividade numérica totalmente não lineares, testando se a instabilidade proposta é real ou um artefato numérico.

2. Metodologia

Os autores realizaram simulações de evolução temporal (3+1 e 2+1 dimensões) de soluções de estrelas de bósons estacionárias, perturbando-as para observar se o sistema se desvia do estado estacionário.

• Modelos Estudados: Focaram no modelo S005 (anteriormente reportado como instável) e no modelo S02. Ambos são estrelas de bósons escalares com auto-interações solitônicas e possuem anéis de luz estáveis.
• Códigos Numéricos: Utilizaram dois códigos independentes de relatividade numérica para garantir robustez:
  1. ExoZvezda: Baseado no sistema CCZ4 (Conformal and Covariant Z4) com gauge de "moving puncture" e malha adaptativa (AMR).
  2. GHC (Generalized Harmonic Code): Baseado na formulação Generalized Harmonic com gauge estacionária.
• Tipos de Perturbações:
  1. Erro de Discretização: Deixaram o erro numérico atuar como a única perturbação inicial.
  2. Perturbações Explícitas:
     • Perturbação azimutal não eixo-simétrica: $\phi \to \phi e^{-2\epsilon \cos(n\phi)}$.
     • Perturbação Gaussiana localizada no anel de luz estável: $|\phi| \to |\phi| + \epsilon e^{-(r-r_0)^2/\Delta^2}$.
• Diagnósticos: Monitoraram a estacionaridade através de:
  • Máximo global da derivada temporal da amplitude do campo escalar ($\Phi$).
  • Máximo da derivada temporal do componente $g_{tt}$ da métrica ($G$).
  • Decomposição modal azimutal para detectar modos instáveis.
  • Violações de vínculos (constraints) para distinguir instabilidades físicas de artefatos numéricos.

3. Contribuições e Descobertas Chave

A. Ausência de Instabilidade Física

• Resultado Principal: Não foi encontrada nenhuma evidência de instabilidade (linear ou não linear) para o modelo S005 nas escalas de tempo simuladas ($t\mu \sim 10^4$).
• Comportamento das Perturbações: As perturbações iniciais, mesmo com amplitudes grandes (mas abaixo do limiar de colapso imediato), exibiram um decaimento lento ou permaneceram estáveis, sem levar à formação de buracos negros ou a um crescimento exponencial descontrolado.
• Contraste com Trabalhos Anteriores: Os resultados contradizem as conclusões da Ref. [14], sugerindo que a instabilidade reportada anteriormente pode não ser um fenômeno físico intrínseco a esses objetos, ou que as escalas de tempo necessárias para observá-la (se existirem) são muito maiores do que as simuladas.

B. Identificação de Instabilidades Não Físicas (Artefatos)

O trabalho destaca desafios cruciais na simulação desses objetos extremos:

1. Instabilidades de Gauge: Ao utilizar a formulação Generalized Harmonic com gauge estacionária, os autores observaram um crescimento exponencial de diagnósticos dependentes de gauge (como $\Phi$ e $G$). No entanto, quantidades invariantes de gauge (como a amplitude máxima do campo escalar normalizada) permaneceram constantes. Isso indica que o crescimento era um artefato da escolha de coordenadas, não uma instabilidade física.
2. Importância do Amortecimento de Vínculos (Constraint Damping): Simulações usando a formulação BSSNOK (sem amortecimento de vínculos) mostraram um crescimento gradual de violações de vínculos, levando a instabilidades numéricas. A formulação CCZ4, com amortecimento adequado, foi essencial para manter a estabilidade numérica a longo prazo.

C. Análise de Convergência e Robustez

• Os resultados foram verificados com múltiplas resoluções, mostrando convergência de ordem 3 a 4.
• A estabilidade observada persistiu independentemente do código utilizado (CCZ4 vs. Generalized Harmonic) e da dimensão da simulação (2+1 vs 3+1), desde que as instabilidades de gauge e de vínculos fossem tratadas corretamente.

4. Significância e Conclusões

• Validação de Objetos Exóticos: Os resultados reforçam a possibilidade de que estrelas de bósons ultracompactas giratórias com anéis de luz estáveis podem ser objetos astrofísicos viáveis e dinamicamente estáveis, desafiando a ideia de que o aprisionamento estável de luz inevitavelmente leva a uma instabilidade catastrófica em escalas de tempo acessíveis.
• Cuidados em Relatividade Numérica: O artigo serve como um guia técnico importante, demonstrando como distinguir entre instabilidades físicas reais e artefatos numéricos (especialmente instabilidades de gauge e de vínculos) ao estudar objetos ultracompactos.
• Limitações e Futuro: Os autores enfatizam que não provaram a estabilidade não linear absoluta (pois escalas de tempo muito longas ou tipos específicos de perturbação ainda não foram testados). Se uma instabilidade existir, ela deve ser extremamente fraca, com escalas de tempo muito superiores a $10^4/\mu$.
• Recomendação: Para estudos futuros, sugere-se investigar o impacto do aprisionamento nulo estável começando com configurações de matéria explicitamente presas no anel de luz para estudar a reação gravitacional (back-reaction) de forma mais direta.

Em resumo, o trabalho fornece uma reavaliação rigorosa que, até o momento, não apoia a existência da instabilidade não linear do anel de luz para estrelas de bósons ultracompactas nas condições testadas, atribuindo relatos anteriores a possíveis limitações metodológicas ou escalas de tempo insuficientes.