Accurate modeling for 3$\times$2pt analyses in Roman and Rubin: a study of model approximations

기술 요약: Roman 및 Rubin의 3×2pt 분석을 위한 정확한 모델링

문제 정의
현대 우주론은 은하 클러스터링(galaxy clustering), 은하-은하 렌징(galaxy-galaxy lensing), 우주 전단(cosmic shear)을 결합하여 우주론적 매개변수를 제약하는 "3×2pt" 분석에 크게 의존하고 있습니다. 이 프레임워크는 Stage-III 조사(예: DES, KiDS)에는 성공적으로 적용되었으나, 다가오는 Stage-IV 조사, 특히 Nancy Grace Roman 우주 망원경의 고위도 이미징 조사(HLIS)와 Vera C. Rubin 천문대의 LSST(Legacy Survey of Space and Time)는 전례 없는 통계적 힘을 제공할 것입니다. 저자들은 Stage-III 데이터에는 허용 가능했던 모델링 근사치들이 Stage-IV 제약 조건에서 상당한 편향(bias)을 유발할 수 있다고 주장합니다. 본 논문은 3×2pt 분석에서 흔히 사용되는 세 가지 특정 모델링 근사치를 체계적으로 조사합니다:

1. 3차원 파워 스펙트럼을 각도 척도로 투영하기 위한 **Limber 근사(Limber approximation)**의 적용.
2. 은하 과밀도 장(galaxy overdensity field)에서 **적색편이 공간 왜곡(Redshift-Space Distortions, RSD)**을 무시하는 것.
3. **비선형 물질 파워 스펙트럼(nonlinear matter power spectrum)**에 대한 덜 정확한 모델의 사용(특히 HMCode와 EuclidEmulator2의 비교).

방법론
본 연구는 Roman HLIS 우주론 PIT에서 개발된 파이프라인인 CoCoA(CosmoLike-Cobaya Architecture)로 구현된 시뮬레이션 우도 분석 프레임워크를 활용합니다. 견고성을 보장하기 위해, 저자들은 먼저 Rubin LSST의 DESC(Dark Energy Science Collaboration) 표준인 **CCL(Core Cosmology Library)**을 대상으로 CoCoA를 검증했습니다.

• 설계 조사 구성: 분석은 5년 치 Roman HLIS 데이터셋(medium tier)과 10년 치 Rubin LSST 데이터셋에 대한 제약 조건을 예측합니다. 두 조사 모두 10개의 렌즈-은하 빈(bin)과 다양한 수의 소스-은하 빈(Roman의 경우 8개, Rubin의 경우 5개)을 사용합니다.
• 모델링 프레임워크: 분석은 5개의 자유 매개변수($A_s, n_s, H_0, \Omega_b, \Omega_m$)를 갖는 평탄한 $\Lambda$CDM 우주론을 가정합니다. 섭동 매개변수(nuisance parameters)로는 광학적 적색편이(photo-z) 이동, 전단 교정 편향, 고유 정렬(intrinsic alignment, NLA 모델), 선형 은하 편향이 포함됩니다.
• 스케일 컷(Scale Cuts): 바리온 피드백 및 비선형 은하 편향의 불확실성을 완화하기 위해 저자들은 보수적인 스케일 컷을 적용합니다. 실공간(real space)에서 이는 렌징의 경우 최소 횡방향 반경 6 Mpc/h, 클러스터링의 경우 8 Mpc/h에 해당합니다. 조화 공간(harmonic space)에서는 특정 다중극 범위($20 < \ell < 4000$)로 변환됩니다.
• 실험 설계: 저자들은 근사의 효과를 격리하기 위해 다섯 가지 모델 시나리오를 정의합니다:
  1. Baseline: 비-Limber(Non-Limber) 클러스터링, RSD 포함, EuclidEmulator2를 사용한 비선형 파워 스펙트럼.
  2. Limber: 클러스터링에 Limber 근사를 사용.
  3. noRSD: RSD 보정을 생략.
  4. HMCode: EuclidEmulator2 대신 HMCode를 사용.
  5. Combined: Limber, RSD 생략, 그리고 HMCode를 모두 사용.

주요 결과
본 연구는 이러한 근사치들이 데이터 벡터(via $\Delta\chi^2$)와 그에 따른 우주론적 제약($\Omega_m - S_8$ 평면에서의 편향)에 미치는 영향을 정량화합니다.

1. 코드 검증: CoCoA와 CCL의 비교는 기준 설정에서 $\Delta\chi^2$ 값이 1 미만으로 나타나 매우 우수한 일치를 보였으며, 이는 모델링 파이프라인의 신뢰성을 확인시켜 줍니다.
2. 개별 근사의 영향:
   • Limber 근사: 은하 클러스터링에 대해 비-Limber 계산을 무시하면 Roman의 경우 ~0.7–0.8$\sigma$, Rubin의 경우 ~1.8–2.0$\sigma$의 편향이 발생합니다. 이 효과는 조화 공간에서 영향을 받는 대규모 모드(large-scale modes)가 스케일 컷에 의해 덜 제거되기 때문에 실공간보다 조화 공간에서 더 두드러집니다.
   • 적색편이 공간 왜곡 (RSD): RSD를 무시하면 Limber 근사와 비슷하거나 약간 더 큰 편향을 생성하며, Roman의 경우 ~0.7–0.8$\sigma$, Rubin의 경우 ~2.0–2.3$\sigma$에 달합니다. 저자들은 Limber 근사가 본질적으로 RSD 기여를 억제하기 때문에 RSD가 두 효과 중 더 지배적인 효과라고 언급합니다.
   • 비선형 파워 스펙트럼: EuclidEmulator2에서 HMCode로 전환하면 더 작은 편향(Roman ~0.6$\sigma$, Rubin ~1.0–1.2$\sigma$)이 발생합니다. 이는 두 모델 간의 차이가 주로 소규모 척도에서 발생하며, 이 영역이 기준 스케일 컷에 의해 대부분 제거되었기 때문입니다.
3. 결합 효과: 세 가지 근사를 동시에 적용하면 편향이 누적됩니다. Rubin의 경우 결합된 편향은 실공간에서 ~2.8$\sigma$, 조화 공간에서 ~3.1$\sigma$에 도달합니다. Roman의 경우 결합된 편향은 ~1.0–1.1$\sigma$입니다.
4. 조사 민감도: Rubin의 제약 능력이 더 크기 때문에 편향은 일관되게 Roman보다 Rubin에서 더 크게 나타납니다. 저자들은 Roman에 대해 "medium tier"를 가정했지만, "wide tier" 구성을 사용할 경우 더욱 엄격한 제약이 가능하며 근사치를 사용할 경우 잠재적으로 더 큰 편향이 발생할 수 있다고 언급했습니다.

의의 및 결론
본 논문은 Stage-IV 조사에 있어 Stage-III 분석에서 사용된 표준 근사치들이 불충분하다고 결론짓습니다. 조사된 세 가지 효과 중 어느 하나라도 무시하면 Roman의 경우 1$\sigma$에 근접하거나 초과하는 편향을, Rubin의 경우 2$\sigma$에 달하는 편향을 유발할 수 있습니다.

저자들은 향-Stage-IV 분석을 위한 세 가지 주요 시사점을 강조합니다:

1. 비-Limber 알고리즘: Limber 근사가 고정밀 조사에서의 은하 클러스터링에 더 이상 충분하지 않으므로, 비-Limber 알고리즘이 수치적으로 안정적이고 계산 효율적인지 최적화하는 것이 중요합니다.
2. RSD 모델링: 선형 차수의 RSD 효과는 데이터 벡터와 제약 조건에 큰 영향을 미치며, 이는 더 현실적인 RSD 모델의 개발과 포함을 촉구합니다.
3. 비선형 파워 스펙트럼: 비선형 물질 파워 스펙트럼 모델의 선택은 특정 스케일 컷 전략과 결합될 때 매우 중요합니다.

이 연구는 CCL에 대한 CoCoA 파이프라인의 검증 역할을 수행하며, Roman 및 Rubin 조사의 최종 데이터 릴리스를 위해 더 높은 충실도(fidelity)를 갖춘 모델링이 필요하다는 정량적 근거를 제공합니다. 저자들은 향-연구가 비선형 은하 편향, 더 복잡한 고유 정렬 모델, 바리온 피드백을 포함하도록 확장되어야 하며, Roman의 medium 및 wide tier 조합을 조사해야 한다고 제안합니다.