Exponential $f(R)$ cosmology with massive neutrinos as a dynamical dark energy framework

Resumo Técnico: Cosmologia $f(R)$ Exponencial com Neutrinos Massivos

Enunciado do Problema
O modelo cosmológico padrão $\Lambda$CDM, embora bem-sucedido, enfrenta desafios teóricos e observacionais persistentes. Os problemas teóricos incluem os problemas de ajuste fino (fine-tuning) e de coincidência associados à constante cosmológica $\Lambda$, e a falta de uma partícula confirmada de matéria escura. Observacionalmente, existem tensões significativas, notadamente a "tensão de Hubble", uma discrepância de $\sim 5\sigma$ entre medições locais da escada de distâncias ($H_0 \approx 73,04$ km s$^{-1}$ Mpc$^{-1}$) e inferências do Universo primordial a partir da Radiação Cósmica de Fundo (CMB) dentro do arcabouço $\Lambda$CDM ($H_0 \approx 67,36$ km s$^{-1}$ Mpc$^{-1}$). Adicionalmente, dados recentes de Oscilações Acústicas de Bárions (BAO) da colaboração DESI sugerem desvios do $\Lambda$CDM, insinuando uma energia escura dinâmica.

Uma tensão secundária envolve a massa total de neutrinos ($\sum m_\nu$). Experimentos de oscilação de neutrinos estabelecem um limite inferior de $\sum m_\nu \gtrsim 0,06$ eV, enquanto dados cosmológicos dentro do modelo $\Lambda$CDM frequentemente resultam em limites superiores abaixo desse limiar (por exemplo, $\sum m_\nu < 0,064$ eV a partir do DESI). Modelos fenomenológicos como $w_0w_a$CDM podem relaxar esses limites, mas carecem de uma base teórica sólida. Este trabalho investiga se uma modificação da gravidade teoricamente motivada, especificamente a gravidade $f(R)$ exponencial, combinada com neutrinos massivos, pode simultaneamente abordar a tensão de Hubble e as restrições de massa de neutrino, mantendo a consistência com o modelo padrão em altos redshifts.

Metodologia
Os autores empregam um modelo de gravidade $f(R)$ exponencial generalizado definido pelo Lagrangiano:
$$f(R) = R - 2\Lambda \left\{ 1 - \exp\left[ -\beta \left( \frac{R}{2\Lambda} \right)^\alpha \right] \right\}$$
onde $\alpha$ e $\beta$ são parâmetros adimensionais que controlam a transição da Relatividade Geral (GR) em altas curvaturas para a gravidade modificada em baixas curvaturas. O modelo inclui um setor de matéria composto por matéria escura fria, bárions, radiação e neutrinos massivos. A densidade de energia dos neutrinos é computada via uma integral de espaço de fase sobre a distribuição de Fermi–Dirac, contabilizando a transição dos regimes relativístico para o não-relativístico.

A dinâmica cosmológica é derivada das equações de Friedmann modificadas em uma métrica FLRW espacialmente plana. Para facilitar a análise numérica, os autores mapeiam os parâmetros $f(R)$ para parâmetros equivalentes de $\Lambda$CDM ($H_0^*, \Omega_i^*$) usando um parâmetro de Hubble adimensional $E(z) = H(z)/H_0^*$. O sistema de equações diferenciais acopladas para $E(z)$ e o escalar de Ricci adimensional $\mathcal{R}(z)$ é integrado numericamente de uma condição inicial de alto redshift (onde o modelo recupera o $\Lambda$CDM) até $z=0$.

Uma análise de Cadeia de Markov Monte Carlo Bayesiana (MCMC) é realizada utilizando o pacote emcee para restringir o espaço de parâmetros $(H_0^*, \alpha, \beta, \omega_b^*, \omega_{bc}^*, \omega_\nu^*)$. A análise utiliza quatro conjuntos de dados observacionais:

1. Cronômetros Cósmicos (CC): 31 pontos de dados medindo $H(z)$.
2. DESI DR2 BAO: Medições de $D_M(z)/r_d$ e $D_H(z)/r_d$ através de seis bins de redshift.
3. CMB: Verossimilhanças comprimidas da Planck 2018 para a escala acústica $\theta_*$ e densidades físicas $\omega_b, \omega_{bc}$.
4. Pantheon+ (Pan+): 1550 módulos de distância de Supernovas Tipo Ia.

Os resultados são comparados contra o modelo padrão $\Lambda$CDM e o modelo fenomenológico $w_0w_a$CDM usando o Critério de Informação de Akaike (AIC) e o Critério de Informação Bayesiano (BIC) para avaliar a preferência estatística.

Principais Contribuições e Resultados
O estudo fornece a primeira restrição conjunta sobre os parâmetros do modelo $f(R)$ exponencial generalizado ($\alpha, \beta$) e a massa total de neutrino ($\sum m_\nu$) utilizando os dados mais recentes do DESI DR2.

• Viabilidade do Modelo: O modelo $f(R)$ exponencial permanece consistente com as observações atuais. As condições de estabilidade $f_R > 0$ e $f_{RR} > 0$ são satisfeitas em todo o espaço de parâmetros posterior e no intervalo de redshift, confirmando a ausência de instabilidades de tipo ghost e taquionicas.
• Tensão de Hubble: O modelo alivia parcialmente a tensão de Hubble. Ao incluir as supernovas Pantheon+, o $H_0$ inferido para o modelo $f(R)$ é $70,89 \pm 0,19$ km s$^{-1}$ Mpc$^{-1}$. Este valor reduz a tensão com as medições locais comparado ao $\Lambda$CDM ($71,43 \pm 0,13$) e ao $w_0w_a$CDM ($71,13 \pm 1,20$), embora uma discrepância de $2,3\sigma$ permaneça entre os valores do $f(R)$ e do $\Lambda$CDM na combinação total dos dados.
• Restrições de Massa de Neutrino: A inclusão de neutrinos massivos no arcabouço $f(R)$ leva a um relaxamento suave do limite superior de $\sum m_\nu$ em comparação ao $\Lambda$CDM. Com os dados do Pantheon+, o limite superior de $1\sigma$ é $\sum m_\nu < 0,0290$ eV para $f(R)$, comparado a $< 0,0126$ eV para $\Lambda$CDM e $< 0,0613$ eV para $w_0w_a$CDM. Embora este seja um relaxamento, o limite permanece abaixo do limite inferior de oscilação ($\sim 0,06$ eV), indicando que a tensão não foi totalmente resolvida.
• Preferência Estatística:
  • Sem os dados do Pantheon+, o $\Lambda$CDM é estatisticamente preferido.
  • Com a inclusão dos dados do Pantheon+, tanto o $w_0w_a$CDM ($\Delta$AIC = 43,75) quanto o $f(R)$ exponencial ($\Delta$AIC = 20,46) mostram evidência "muito forte" e "forte", respectivamente, sobre o $\Lambda$CDM de acordo com as métricas AIC e BIC.
• Degenerescências de Parâmetros: A análise revela fortes degenerescências entre os parâmetros de modificação ($\alpha, \beta$) e os parâmetros cosmológicos ($H_0, \Omega_m$), indicando que a gravidade modificada pode compensar parcialmente mudanças na densidade de matéria e na história da expansão. A distribuição posterior para $\beta$ é não-gaussiana, refletindo o papel não linear da correção exponencial.

Significância e Alegações
O artigo afirma que a interação entre gravidade modificada e física de neutrinos oferece um arcabouço viável e teoricamente robusto para investigar as tensões cosmológicas. Diferente das parametrizações fenomenológicas de energia escura, o modelo $f(R)$ exponencial fornece uma base teórica sólida enraizada em modificações do Lagrangiano gravitacional.

Os autores enfatizam que, embora o modelo não resolva totalmente a tensão de Hubble ou a discrepância da massa de neutrino, ele demonstra que uma abordagem combinada — modificando o setor gravitacional enquanto inclui neutrinos massivos — pode melhorar o ajuste estatístico aos dados observacionais e relaxar as restrições sobre a massa de neutrino de forma mais eficaz do que o $\Lambda$CDM sozinho. O trabalho destaca a importância de conjuntos de dados de baixo redshift (como SNe Ia) para discriminar entre modelos alternativos e a cosmologia padrão. Os autores concluem que investigações futuras devem incluir efeitos perturbativos (por exemplo, lentes fracas, distorções de espaço de redundância) para avaliar plenamente a viabilidade do modelo, uma vez que a análise atual está restrita à cosmologia de fundo.