Hu-Sawicki $f(R)$ Gravity in a Non-Flat Universe: Constraints from DESI-DR2, BBN, and Type Ia Supernovae

기술 요약: 비평탄 우주에서의 Hu-Sawicki $f(R)$ 중력

문제 제기
표준 $\Lambda$CDM 우주론 모델은 성공적이지만, 우주 상수 문제와 지속적인 허블 텐션($H_0$)과 같은 이론적 과제에 직면해 있다. 대안적인 중력 이론, 특히 $f(R)$ 중력은 암흑 에너지 성분을 도입하는 대신 중력 라그랑지안을 수정함으로써 후기 우주의 가속 팽창을 설명하는 메커니즘을 제공한다. 그러나 실행 가능한 대부분의 $f(R)$ 모델 분석은 공간적으로 평탄한 우주($\Omega_k = 0$)를 가정한다. 인플레이션 모델이 열린 우주나 닫힌 우주를 생성할 수 있다는 점을 고려할 때 이러한 가정은 성급할 수 있으며, 최근 Planck 2018 데이터는 $\Lambda$CDM 내에서 미세하게 닫힌 우주를 선호하는 경향을 보였다. 더욱이, 암흑 에너지 파라미터와 공간 곡률 사이에는 강력한 퇴화(degeneracy)가 존재한다. 본 연구는 비평탄 배경에서의 실행 가능한 $f(R)$ 모델에 대한 상세한 분석이 부족하다는 점을 다루며, 특히 공간 곡률을 자유롭게 변수로 허용했을 때 수정 중력 파라미터에 대한 제약이 어떻게 변화하는지, 그리고 현재의 관측 데이터를 사용하여 이 모델들이 $\Lambda$CDM과 어떻게 비교되는지를 조사한다.

방법론
저자들은 다음과 같이 정의된 구체적이고 실행 가능한 일반 상대성 이론의 확장인 Hu-Sawicki $f(R)$ 중력 모델을 분석한다:
$$f(R) = R - c_1 \frac{R_{HS} (R/R_{HS})^p}{c_2 (R/R_{HS})^p + 1}$$
여기서 $p=1$로 고정되며, 모델은 무차원 수정 중력 파라미터 $b$ (원래의 파라미터 $c_1, c_2, R_{HS}$와 관련됨)에 의해 특징지어진다. $b \to 0$ 극한은 $\Lambda$CDM 모델을 회복한다.

본 연구는 수정 중력 파라미터 $b$와 공간 곡률 밀도 파라미터 $\Omega_k$를 모두 자유 파라미터로 취급한다. 분석에는 다음의 관측 데이터셋이 활용된다:

1. DESI-DR2: 적색편프트 $0.295 \le z \le 2.330$ 범위를 커버하는 Dark Energy Spectroscopic Instrument(두 번째 데이터 릴리스)의 바리온 음향 진동(BAO) 측정값.
2. BBN: 원시 중수소 및 헬륨 함량에 기반한 빅뱅 핵합성(BBN) 제약 조건.
3. Type Ia 초신성 (SN Ia): 견고성을 테스트하기 위해 사용된 네 가지 별도의 컴필레이션:
   • PantheonPlus (PP)
   • Cepheid 주계열성 거리로 보정된 PantheonPlus+SH0ES (PPS)
   • Union3
   • DESY5 (Dark Energy Survey Year 5)

이론적 예측은 $f(R)$ 시나리오를 위해 수정된 볼츠만 솔버 CLASS와 MCMC 샘플러 MontePython을 연동하여 생성되었다. 제약 조건은 68% 및 95% 신뢰 수준에서 도출되었다. 모델 비교는 추가적인 자유 파라미터 $b$가 모델의 개선된 적합도를 정당화하는지 평가하기 위해 아카이케 정보 기준(AIC)과 베이지안 정보 기준(BIC)을 사용하여 수행되었다.

주요 기여 및 결과

• 허블 상수 및 물질 밀도: 네 가지 데이터 조합 모두에서 Hu-Sawicki 모델의 허블 상수($H_0$)와 물질 밀도($\Omega_m$)는 $\Lambda$CDM 값과 유사하게 유지된다. 주목할 만한 $H_0$의 변화는 PPS 조합(SH0ES로 보정됨)에서 나타나는데, 이는 $H_0 \approx 71.9$ km/s/Mpc를 기록하며, 이는 $f(R)$ 섹터의 고유한 특징이라기보다 보정 소스의 특성을 반영한다.
• 수정 중력 파라미터($b$)에 대한 제약:
  • 네 가지 조합 중 세 가지(PP, Union3, DESY5yr)에서 데이터는 $b \neq 0$을 $2\sigma$ 이상의 유의도로 선호하며, DESY5yr이 가장 뚜렷한 편차($b = 0.501^{+0.11}_{-0.092}$)를 보인다.
  • PPS 조합은 예외적으로 음수 값($b = -0.26^{+0.18}_{-0.16}$)을 선호하는데, 저자들은 이를 팽창사 요구 사항을 변화시키는 특정 SH0ES 보정 때문이라고 설명한다.
  • $b$의 부호와 크기는 현재 데이터에 의해 유일하게 결정되지 않으며, 사용된 SN Ia 보정에 크게 의존한다.
• 공간 곡률 ($\Omega_k$):
  • $b$를 자유롭게 두었을 때, $\Omega_k$는 네 가지 조합 모두에서 평탄한 우주($\Omega_k = 0$)와 $1\sigma$ 수준에서 일치한다.
  • 이는 동일한 데이터에 대한 $\Lambda$CDM 적합 결과와 극명하게 대조된다. $\Lambda$CDM에서는 PPS가 닫힌 우주($\Omega_k = -0.063 \pm 0.020$)를 약간 선호하고, DESY5yr은 열린 우주를 선호한다.
  • 저자들은 $b, H_0, \Omega_k$ 사이에 강한 양의 상관관계가 있음을 확인했다. 이는 $b$라는 추가적인 자유도가 허용될 때, 수정 중력 섹터로 곡률 신호가 "재흡수"될 수 있음을 시사한다. 즉, $\Lambda$CDM 하에서 탐지된 곡률 징후는 수정 중력 모델에 의해 상쇄될 수 있다.
• 통계적 모델 비교:
  • AIC: 네 가지 사례 중 세 가지에서 Hu-Sawicki 모델 쪽으로 기울어지며 (DESY5yr은 $\Delta$AIC = 7.78로 강한 증거를 보임), 이는 모델이 더 나음을 시사한다.
  • BIC: 추가 파라미터에 대한 더 엄격한 페널티로 인해, BIC는 대부분의 경우 $\Lambda$CDM을 선호한다. 오직 DESY5yr 조합만이 BIC 하에서 Hu-Sawicki 모델에 대한 미미한 선호를 보인다 ($\Delta$BIC = 2.26).
  • PPS 조합은 두 기준 모두에서 일관되게 $\Lambda$CDM을 선호한다.

의의 및 주장
본 논문은 곡률이 제약되지 않은 Hu-Sawicki $f(R)$ 시나리오가 $\Lambda$CDM에 대한 통계적으로 실행 가능한 대안이라고 주장하지만, 현재 데이터에 의해 명확하게 선호되는 것은 아니다. 주요 의의는 곡률 제약이 기초가 되는 중력 모델과 독립적이지 않다는 발견에 있다. $\Lambda$CDM 적합에서 나타나는 비제로 곡률에 대한 선호는, $b$라는 수정 중력 파라미터가 도입될 때 곡률과 중력 수정 사이의 기하학적 퇴화로 인해 완화되거나 제거될 수 있다.

저자들은 현재의 DESI-DR2, BBN, SN Ia 데이터가 모든 데이터셋에 걸쳐 $\Lambda$CDM보다 Hu-Sawicki 모델을 통계적으로 결정적으로 선호한다는 증거를 제공하지 못한다고 결론짓는다. 일부 조합에서 발견된 $b \neq 0$에 대한 미미한 선호는 일반 상대성 이론으로부터의 실제적인 이탈이라기보다는 관측의 한계나 보정 계통 오차를 반영할 수 있음을 강조한다. 관찰된 편차가 물리적인 것인지 판단하기 위해서는 BAO 정밀도의 향상과 초신성 보정 계통 오차에 대한 더 엄격한 통제가 필요하다.