Bayesian Inference of Neutron Star Properties in $f(Q)$ Gravity Using NICER Observations

Resumo Técnico: Inferência Bayesiana de Propriedades de Estrelas de Neutrons em Gravidade $f(Q)$ usando Observações do NICER

Definição do Problema
Embora a Relatividade Geral (RG) descreva com sucesso fenômenos gravitacionais em várias escalas, ela enfrenta desafios conceituais relativos a singularidades, matéria escura e energia escura, o que motiva a exploração de teorias de Gravidade Modificada (GM). Entre elas, a gravidade teleparalela simétrica (STEGR) e sua extensão, a gravidade $f(Q)$, oferecem uma formulação geométrica distinta onde a gravidade é codificada na não-metricidade do espaço-tempo, em vez da curvatura ou torção. Embora os modelos $f(Q)$ tenham sido estudados na cosmologia, sua aplicação ao regime de campo forte de estrelas de neutrones (ENs) permanece amplamente inexplorada. Especificamente, há uma carência de estudos sistemáticos de inferência bayesiana para restringir os parâmetros dos modelos $f(Q)$ usando dados observacionais, visto que análises anteriores frequentemente dependiam de escolhas fenomenológicas de parâmetros sem um confronto estatístico rigoroso com os dados.

Metodologia
Os autores investigam a estrutura interna de estrelas de neutrones dentro do arcabouço da gravidade $f(Q)$, considerando três formas funcionais representativas: modelos lineares, logarítmicos e exponenciais. Para isolar os efeitos da gravidade modificada das incertezas da física da matéria densa, o estudo adota a equação de estado (EoS) realista e dependente da densidade DDME2 de campo médio relativístico.

A metodologia procede da seguinte forma:

1. Arcabouço Teórico: Os autores derivam as equações de Tolman–Oppenheimer–Volkoff (TOV) modificadas para espaços-tempos estáticos e esfericamente simétricos em gravidade $f(Q)$. Estas equações governam o equilíbrio estelar e são resolvidas numericamente utilizando um esquema de Runge-Kutta de quarta ordem (RK4).
2. Inferência Bayesiana: Um arcabouço bayesiano é empregado para restringir os parâmetros livres ($\alpha, \beta, Q_0$) dos três modelos $f(Q)$. A análise utiliza distribuições posteriores de massa-raio de observações do NICER de quatro pulsares: PSR J0030+0451, PSR J0740+6620, PSR J0437+4715 e PSR J0614+3329.
3. Análise Estatística: O estudo computa as distribuições posteriores para os parâmetros do modelo e propriedades derivadas das ENs (massa, raio, deformabilidade tidal e compactação). A comparação de modelos é realizada utilizando o fator de Bayes (BF) para determinar qual forma funcional é estatisticamente preferida pelos dados.

Principais Contribuições

• Primeiras Restrições Bayesianas em $f(Q)$: Este trabalho representa a primeira análise sistemática de inferência bayesiana de parâmetros de modelos de gravidade $f(Q)$ utilizando observáveis de estrelas de neutrones, indo além dos testes cosmológicos para o regime de campo forte.
• Análise Comparativa de Modelos: O artigo fornece uma comparação estatística rigorosa entre os modelos $f(Q)$ lineares, logarítmicos e exponenciais, avaliando sua viabilidade frente aos dados de alta precisão do NICER.
• Degenerescência de Parâmetros e Correlação: O estudo mapeia as correlações entre os parâmetros do modelo e as propriedades das ENs, revelando como diferentes formas funcionais influenciam o acoplamento entre as modificações gravitacionais e a estrutura estelar.

Resultos

• Preferência de Modelo: A análise bayesiana, apoiada por comparações de fator de Bayes, favorece estatisticamente o modelo $f(Q)$ exponencial sobre os casos linear e logarítmico. O modelo exponencial exibe parâmetros bem restritos com distribuições posteriores quase gaussianas e unimodais. Em contraste, o modelo linear mostra bimodalidade pronunciada e distribuições de parâmetros amplas, indicando fortes degenerescências, enquanto o modelo logarítmico apresenta incertezas moderadas.
• Propriedades de Estrelas de Neutrons: Para o modelo exponencial preferido, o estudo deriva um raio e uma deformabilidade tidal em $1.4 M_\odot$ de $R_{1.4} = 11.27^{+0.53}_{-0.36}$ km e $\Lambda_{1.4} = 156.95^{+84.02}_{-41.73}$, respectivamente. Estes valores são consistentes com as restrições observacionais atuais.
• Massa Máxima e o Gap de Massa: Um achado significativo é que todos os três modelos restritos preveem massas máximas de estrelas de neutrones atingindo $M_{\text{max}} \simeq 2.98 M_\odot$. As regiões de confiança de 95% para a relação massa-raio estendem-se para o gap de massa inferior ($\sim 2.5\text{--}5 M_\odot$). Esta previsão emerge naturalmente do espaço de parâmetros restringido por via bayesiana, sem ser imposta por construção.
• Correlações: O modelo exponencial demonstra uma forte interdependência estrutural entre a massa máxima e outras propriedades das ENs, ao passo que o modelo linear mostra independência estatística entre a massa máxima e outros observáveis.

Significância
O artigo afirma que estes resultados destacam o potencial das estrelas de neutrones como sondas poderosas da gravidade teleparalela simétrica no regime de campo forte. A emergência de configurações estáveis no gap de massa é apresentada como uma previsão clara e testável do arcabouço $f(Q)$, distinguindo-o de cenários padrão de RG que tipicamente produzem massas máximas mais baixas para a mesma EoS. O estudo conclui que o modelo $f(Q)$ exponencial é a descrição estatisticamente mais suportada da estrutura das estrelas de neutrones entre os três testados, sugerindo que os efeitos da gravidade modificada podem significar um aumento significativo no suporte gravitacional de alta densidade. Os autores observam que medições de massa precisas no futuro na região do gap de massa poderiam servir como um poderoso discriminador entre a RG e cenários de gravidade estendida. Eles também sugerem que trabalhos futuros poderiam envolver a inferência bayesiana conjunta com parametrizações flexíveis de EoS para abordar ainda mais as degenerescências entre a gravidade de campo forte e a microfísica.