The contribution to Galactic Centre {\gamma}-ray excess from cluster-born millisecond pulsars. Constraints from direct N-body simulations

Utilizando simulações N-corpos de alta resolução, o estudo demonstra que uma população de pulsares de milissegundos originada em aglomerados globulares (atuais e disruptados) é capaz de reproduzir as propriedades do excesso de raios gama no Centro Galáctico, favorecendo uma origem astrofísica em detrimento da interpretação por matéria escura.

O essencial

🌌 O Mistério do "Brilho Fantasma" no Centro da Galáxia

Imagine que você está olhando para o centro da nossa galáxia, a Via Láctea. Lá existe um buraco negro supermassivo (Sagittarius A*) e uma densa aglomeração de estrelas. Os astrônomos sabem que essa região deve brilhar em raios gama (uma forma de luz muito energética, invisível aos nossos olhos) devido a colisões de partículas cósmicas.

Mas, quando o telescópio Fermi-LAT olhou para lá, descobriu algo estranho: havia muito mais brilho do que a física padrão previa. Era como se alguém tivesse ligado um holofote extra no centro da cidade, mas ninguém sabia quem era o responsável. Esse fenômeno é chamado de Excesso de Raios Gama do Centro Galáctico (GCE).

Por anos, dois suspeitos principais foram investigados:

1. Matéria Escura: Partículas misteriosas que se aniquilam e liberam energia.
2. Pulsares de Milissegundos (MSPs): "Foguetes" de estrelas mortas que giram super rápido e lançam feixes de luz.

Este novo estudo, feito por uma equipe internacional, decidiu investigar o segundo suspeito: os Pulsares.

🚌 A Analogia dos Ônibus e dos Passageiros

Para entender como os pesquisadores chegaram à conclusão, vamos usar uma analogia de transporte público:

• A Galáxia é uma cidade grande.
• Os Aglomerados Globulares (GCs) são como ônibus antigos e cheios que viajam pela cidade.
• As Estrelas de Nêutrons são os passageiros dentro desses ônibus.
• Os Pulsares são passageiros especiais que, se tiverem sorte, ganham um "superpoder" e começam a brilhar (emitir raios gama).

O Problema

O centro da cidade (o Centro Galáctico) está muito iluminado, mas não sabemos de onde vem essa luz. Será que os "ônibus" (aglomerados globulares) estão soltando passageiros (estrelas de nêutrons) lá no centro, e esses passageiros viram "super-heróis" (pulsares) que iluminam tudo?

A Investigação (Simulações Computacionais)

Os cientistas não puderam esperar 1 bilhão de anos para ver o que acontece. Então, eles criaram um mundo virtual no computador (uma simulação N-corpos).

1. O Cenário: Eles colocaram 12 ônibus reais (aglomerados globulares que ainda existem hoje) e 54 ônibus teóricos (que já foram destruídos há muito tempo) em uma estrada virtual que simula a gravidade da Via Láctea.
2. A Viagem: Eles deixaram esses ônibus viajarem por 8 bilhões de anos no computador.
3. O Efeito: À medida que os ônibus passavam perto do centro da cidade, a gravidade forte "puxava" os passageiros para fora. Alguns passageiros caíam do ônibus e ficavam presos no centro; outros eram jogados para longe.

🔍 O Que Eles Descobriram?

Ao final da simulação, eles contaram quantos "passageiros" (estrelas de nêutrons) ficaram no centro da cidade.

• A Taxa de Conversão: Eles sabiam que nem toda estrela de nêutrons vira um pulsar brilhante. Então, usaram dados reais de aglomerados que já observamos para criar uma "fórmula mágica": Para cada 100 estrelas de nêutrons que caem no centro, quantas viram pulsares brilhantes? Eles calcularam que cerca de 1,1% se tornam pulsares.
• O Resultado:
  • Os ônibus que ainda existem hoje (aglomerados vivos) já soltaram passageiros suficientes para criar um brilho considerável.
  • Mas, quando eles somaram os passageiros dos ônibus destruídos (aqueles que já não existem mais, mas deixaram seus "fantasmas" no centro), o brilho aumentou muito!

🎯 A Conclusão: Quem é o culpado?

O estudo comparou o brilho gerado por esses "ônibus destruídos e vivos" com o brilho real que o telescópio Fermi vê.

• O Veredito: O brilho gerado pelos pulsares vindos desses aglomerados combina perfeitamente com o brilho misterioso que vemos.
• A Forma do Brilho: O estudo também notou que esses pulsares não ficam espalhados aleatoriamente. Eles formam uma forma achatada (como um disco) e também uma forma vertical, o que bate exatamente com o que os telescópios observam no céu.

Em resumo: A equipe concluiu que não precisamos inventar matéria escura para explicar esse brilho. A explicação mais simples e provável é que o Centro Galáctico está cheio de "fósseis" de aglomerados de estrelas destruídos, que deixaram para trás uma multidão de pulsares girando rápido, iluminando a região.

💡 Por que isso é importante?

Imagine que você vê uma luz forte na janela de uma casa. Você pode pensar:

1. "É um fantasma!" (Matéria Escura).
2. "É alguém com uma lanterna potente!" (Pulsares).

Este estudo diz: "Olha, a lanterna faz todo o sentido. A gente sabe que tem muita gente com lanternas lá dentro, e se somarmos todas elas, a luz bate exatamente com o que vemos. Não precisamos inventar fantasmas."

Os autores sugerem que, no futuro, telescópios de rádio (como o FAST ou o SKA) devem procurar por esses "ônibus fantasmas" e seus passageiros no centro da galáxia para confirmar essa teoria. Se encontrarem, teremos resolvido um dos maiores mistérios da astronomia moderna!

Resumo técnico

Resumo Técnico: Contribuição do Excesso de Raios-Gama do Centro Galáctico oriundo de Pulsares de Milissegundos nascidos em Aglomerados

1. O Problema e o Contexto

O Excesso de Raios-Gama do Centro Galáctico (GCE, do inglês Galactic Centre Excess) é um sinal de radiação gama difusa e espacialmente estendida observado pelo telescópio Fermi-LAT na região central da Via Láctea, centrado em Sagitário A*. Este excesso excede as previsões dos modelos astrofísicos padrão de interação de raios cósmicos com o meio interestelar.
As principais hipóteses para a origem do GCE são:

1. Aniquilação de Matéria Escura: Partículas massivas de interação fraca (WIMPs) aniquilando-se no halo de matéria escura do centro galáctico.
2. População de Pulsares de Milissegundos (MSPs): Uma população não resolvida de pulsares de milissegundos, cujas emissões coletivas imitam o espectro e a morfologia observados.

O objetivo deste estudo é reavaliar a interpretação baseada em MSPs, focando especificamente no papel dos Aglomerados Globulares (GCs). A hipótese central é que os GCs, através de sua evolução dinâmica e destruição ao longo do tempo cósmico, depositaram uma grande quantidade de estrelas de nêutrons (NS) no centro galáctico, que posteriormente foram recicladas em MSPs, contribuindo para o sinal de raios gama.

2. Metodologia

Os autores empregaram uma abordagem puramente dinâmica e de alta resolução para rastrear a evolução e a deposição de estrelas de nêutrons:

• Simulações N-corpos Diretas: Utilizaram o código de alta ordem paralelo $\phi$-GPU, que implementa um esquema de integração Hermite de quarta ordem com passos de tempo em blocos hierárquicos. O código inclui evolução estelar atualizada (perda de massa, supernovas, kicks natais).
• Potencial Galáctico: As simulações foram realizadas em um potencial galáctico da Via Láctea variável no tempo, derivado da simulação IllustrisTNG-100, capturando tanto a dinâmica interna dos aglomerados quanto as marés externas.
• Amostras Simuladas:
  1. GCs Atuais: 12 aglomerados globulares reais observados na Via Láctea (ex: Terzan 5, NGC 6642), cujas órbitas foram integradas para trás por 8 Gyr.
  2. GCs Destruídos (Teóricos): Um conjunto de aglomerados teóricos (50 com massa inicial de $6 \times 10^4 M_\odot$ e 4 mais massivos) que foram completamente destruídos ao longo de 5 Gyr, depositando seus restos no centro galáctico.
• Mapeamento NS para MSP: Como as simulações não modelam a evolução binária necessária para "reciclar" estrelas de nêutrons em MSPs, os autores aplicaram um fator de eficiência empírico ($k$). Este fator foi calibrado usando a razão observada entre o número de MSPs e o número de estrelas de nêutrons ligadas em aglomerados reais.
  • Fator mediano encontrado: $k \approx 0,0114$ (ou seja, ~1,14% das NS retidas tornam-se MSPs).
• Cálculo de Fluxo: Converteram a distribuição de NS depositadas no centro (dentro de 1 kpc) em emissão de raios gama, assumindo luminosidades médias por pulsar ($\langle L \rangle$) variando entre $1 \times 10^{33}$e$8 \times 10^{33}$erg s$^{-1}$.

3. Contribuições Principais

• Abordagem Dinâmica Completa: Diferente de estudos anteriores que assumiam frações fixas de retenção de estrelas de nêutrons, este trabalho rastreia individualmente as órbitas de todas as NS formadas, considerando kicks natais, marés e perda de massa, fornecendo uma distribuição de fase completa (posição e velocidade) no centro galáctico.
• Inclusão de Aglomerados Destruídos: O estudo quantifica explicitamente a contribuição de aglomerados que já não existem, mas cujos destroços (incluindo NS) foram acumulados no bojo galáctico, um componente frequentemente negligenciado ou modelado de forma simplista.
• Calibração Estatística Robusta: O uso de uma calibração baseada em dados observacionais (Fermi-LAT e FAST) e previsões teóricas para determinar a taxa de conversão NS-MSP, aplicando-a consistentemente tanto a aglomerados sobreviventes quanto aos destruídos.

4. Resultados Chave

• Distribuição de Estrelas de Nêutrons:
  • A maioria das NS dos aglomerados atuais permanece ligada aos seus aglomerados, mas uma fração significativa é ejetada ou depositada no centro galáctico.
  • As NS provenientes de aglomerados destruídos exibem uma morfologia axissimétrica, com sobredensidades pronunciadas no plano galáctico e uma estrutura menos proeminente perpendicular a ele.
• Reprodução do Excesso de Raios Gama (GCE):
  • Os MSPs derivados apenas dos aglomerados observados atuais já fornecem um sinal significativo.
  • Ao adicionar a contribuição dos aglomerados destruídos e assumir uma luminosidade média de $\langle L \rangle \approx 8 \times 10^{33}$ erg s$^{-1}$, o modelo consegue reproduzir com precisão a amplitude e a concentração central do GCE observado.
  • Para um ajuste ideal com os dados observacionais (subtraindo a contribuição do Buraco Negro Supermassivo), foi necessário aumentar a contribuição dos aglomerados destruídos por um fator de aproximadamente 2.
• Sensibilidade à Luminosidade:
  • Com luminosidades menores ($\langle L \rangle \le 4 \times 10^{33}$ erg s$^{-1}$), o modelo subestima o fluxo observado, sugerindo incertezas na população inicial de aglomerados destruídos ou na eficiência de formação de MSPs.
  • No entanto, mesmo com o número atual de aglomerados destruídos simulados, o uso do coeficiente máximo de conversão ($k$) permite atingir o nível de fluxo observado.

5. Significado e Conclusões

• Favorecimento da Origem Astrofísica: Os resultados favorecem fortemente a origem em MSPs para o GCE em detrimento da aniquilação de Matéria Escura. A contribuição combinada de MSPs entregues por aglomerados sobreviventes e destruídos reproduz naturalmente tanto a intensidade quanto a morfologia do sinal observado sob premissas razoáveis.
• Morfologia do Sinal: A estrutura axissimétrica e as sobredensidades no plano galáctico previstas pela distribuição de NS de aglomerados destruídos são consistentes com características recentes relatadas no GCE, oferecendo uma explicação física para a morfologia do excesso.
• Implicações Futuras: O estudo sugere que buscas profundas por pulsares de rádio no bojo galáctico (com telescópios como FAST, MeerKAT e SKA) ao longo das regiões de sobredensidade previstas poderiam confirmar esta população. Além disso, a modelagem futura deve incorporar a evolução binária completa e atrito dinâmico autoconsistente para refinar as previsões.

Em suma, o trabalho demonstra que a evolução dinâmica de aglomerados globulares é um mecanismo viável e robusto para explicar o Excesso de Raios Gama do Centro Galáctico, sem a necessidade de invocar nova física de partículas.