Improving constraints on primordial non-Gaussianity from Quaia with a new cosmological observable: angular redshift fluctuations

기술 요약: 각도 적색편이 변동(ARF)을 이용한 Quaia의 원시 비가우시안성 제약 개선

문제 및 동기
$f_{\text{NL}}$로 매개변수화되는 원시 비가우시안성(PNG)의 검출은 인플레이션 모델에 대한 결정적인 테스트이며, 여기서 0이 아닌 값은 가우시안 초기 조건으로부터의 편차를 나타낸다. 현재 가장 강력한 제약은 CMB 바리스펙트럼(예: Planck 2018)에서 오고 있지만, 대규모 구조(LSS)는 규모 의존적 편향 효과를 통해 상호 보완적인 탐사 도구를 제공한다. 이 효과는 Dalal 등(2008)에 의해 예측되었으며, 국소적 PNG에 대해 은하 편향의 특징적인 $1/k^2$ 시그니처를 도입한다. 그러나 LSS로부터 이를 측정하는 것은 거대 척도에서의 관측 계통 오차를 정밀하게 제어해야 한다는 점에서 도전적이다.

최근 연구들은 $f_{\text{NL}}$을 제약하기 위해 2D 각도 파워 스펙트럼과 CMB 렌징과의 교차 상관 관계를 활용했다. 본 연구는 새로운 우주론적 관측량인 각도 적색편이 변동(ARFs)을 포함함으로써 이러한 제약을 개선할 수 있는 가능성을 다룬다. ARF는 주어진 껍질(shell) 내에서 평균 적색편이의 2D 각도 편차를 측정한다. 이론적으로 ARF는 은하 편향에 민감하며 따라서 PNG에 민감하지만, 실제 데이터에 $f_{\text{NL}}$ 제약을 위해 이를 적용한 사례는 이전에 입증된 바 없다.

방법론
저자들은 Gaia DR3와 unWISE 데이터를 기반으로 한 약 130만 개의 퀘이사(QSO)로 구성된 전천 샘플인 Quaia 카탈로그를 분석하였다. 이 샘-플은 넓은 우주론적 부피를 포괄한다. 분석 파이프라인은 다음 단계로 구성되었다:

1. 데이터 준비: Quaia 샘플을 두 개의 적색편이 빈($\bar{z} \approx 0.97$ 및 $\bar{z} \approx 2.10$)으로 나누었다. 은하 수 밀도(density)와 ARF에 대한 HEALPix 맵($N_{\text{side}}=128$)을 생성하였다. ARF 필드는 서베이의 고유 적색편이 분포에 상응하는 넓은 가우시안 윈도우를 사용하여 픽셀당 가중 평균 적색편이 변동으로 정의되었다.
2. 교차 상관 관계: 연구는 Quaia 밀도 및 ARF 맵을 Planck PR4 CMB 렌징 수렴 맵과 결합하였다. 이를 통해 네 가지 관측량인 밀도 자기 상관 스펙트럼($C_{gg}^\ell$), ARF 자기 상관 스펙트럼($C_{zz}^\ell$), 그리고 CMB 렌징과의 교차 상관 스펙트럼($C_{\kappa g}^\ell$ 및 $C_{\kappa z}^\ell$)을 얻었다.
3. 이론적 모델링: 수정된 버전의 CAMB 볼츠만 코드(ARFCAMB)를 사용하여 이론적 파워 스펙트럼을 계산하였다. 이 코드는 ARF 커널을 포함하며 광도 적색편이 불확실성($\sigma_z \approx 0.06(1+z)$)을 고려한다. 모델은 보편성 관계(universality relation, $p=1$ 가정) 하에서 $f_{\text{NL}}$에 의해 유도되는 규모 의존적 편향과 은하 편향($b(z)$)의 유연한 매개변수화를 포함하였다.
4. 계통 오차 및 공분산: 관측 계통 오차를 해결하기 위해 저자들은 Quaia 템플릿에 대한 선형 회귀를 이용한 계통 탈의(deprojection) 기법을 사용하였다. 데이터 시프트가 통계적 변동과 호환되지 않는 저-$\ell$ 오프셋을 제거하기 위해, 모크 실현(mock realization)과의 비교를 통해 척도 컷(scale cut)을 결정하였다. 구체적으로 밀도에는 $\ell_{\text{min}}=14$, ARF 자기 스펙트럼에는 \ell_{\text{min}}=10}이 적용되었다. 공분산 행렬은 1,000개의 상관된 시뮬레이션으로부터 도출되었으며, 실제 데이터에서 저 멀티폴(low multipoles)에서 관찰된 과잉 파워에 맞추기 위해 모크 자기 스펙트럼에 보정이 적용되었다.
5. 매개변수 추론: 측정된 파워 스펙트럼과 공분산 행렬을 사용하여 가우시안 우도(likelihood)를 구축하였다. MCMC 샘플링(emcee)을 사용하여 $f_{\text{NL}}$과 편향 매개변수($b_0, \alpha$)를 제약하였다.

주요 결과
분석 결과는 국소적 비가우시안성 매개변수 $f_{\text{NL}}$(보편성 관계 $p=1$ 가정 시)에 대해 다음과 같은 제약을 산출하였다:

• ARF 단독: ARF 자기 상관과 ARF-CMB 렌징 교차 상관을 결합했을 때 $f_{\text{NL}} = -43^{+50}_{-44}$를 얻었다.
• 밀도 단독: 밀도 자기 상관과 밀도-CMB 렌징 교차 상관을 결합했을 때 $f_{\text{NL}} = -25^{+19}_{-17}$을 얻었으며, 이는 이전의 Quaia 분석(Fabbian et al. 2026)과 일치한다.
• 결합 분석: 네 가지 관측량 모두($C_{gg}^\ell, C_{\kappa g}^\ell, C_{zz}^\ell, C_{\kappa z}^\ell$)를 결합했을 때 $f_{\text{NL}} = -3 \pm 14$ (68% 신뢰 구간)의 결과를 얻었다.
• 적합도 최적화: 저자들은 모델의 적합도(goodness-of-fit)를 높이기 위해 ARF-CMB 교차 상관($C_{\kappa z}^\ell$)을 포함했을 때 $\chi^2$에서 미세한 긴장(tension)이 있음을 언급하였다. 이 특정 교차 상관을 제외하여 모델의 적합도를 개선했을 때 $f_{\text{NL}} = -5^{+16}_{-15}$가 되었다.

결합 제약은 밀도와 CMB 렌징만을 사용한 이전의 Quaia 측정값보다 약 25% 향상된 것이다. 이 결과는 LSS 2점 통계량을 사용한 지금까지의 $f_{\text{NL}}$ 제약 중 두 번째로 엄격하며(DESI DR1 분광 데이터 이후), 투영된 2D 요약 통계량을 사용하여 달성된 것 중 가장 정밀하다.

의의 및 강건성
본 논문은 ARF를 우주론적 관측량으로 포함하는 것이 실제 데이터에서 $f_{\text{NL}}$ 제약을 개선할 수 있는 측정 가능한 경로를 제공한다고 주장한다. 저자들은 다음 테스트를 통해 결과의 강건성을 검증하였다:

• 편향 모델링: 데드 피드백 모델(Piccirilli et al. 2024 및 Laurent et al. 2017 모델 테스트)을 변경하거나 밀도 및 ARF 필드에 대해 독립적인 편향 매개변수를 허용했을 때도 결과는 안정적이었다.
• 척도 컷: 제약은 최소 멀티폴($\ell_{\text{min}}$) 컷의 변화에 대해 안정적인 것으로 나타났으나, 편향 매개변수에서 일부 긴장이 관찰되었는데, 이는 ARF 커널이 밀도 커널에 비해 비선형 척도($k$-modes)에 더 높은 민감도를 갖기 때문으로 분석되었다.
• 보편성 관계: $p=1$ 가정을 $p=1.6$(일부 문헌에서 QSO에 대해 제안된 값)으로 완화했을 때, 제약은 $f_{\text{NL}} = 2^{+20}_{-21}$로 이동하거나, 모델 독립적 접근 방식에서 $b_\phi f_{\text{NL}} = -63^{+79}_{-89}$가 되었다.

저자들은 관측량 간의 내부적 긴장(모델링이 필요한 비선형 척도에서의 ARF 모델링 문제로 추정됨)이 존재함에도 불구하고, 2D 분석에 ARF를 추가하는 것이 유망한 방안이라고 결론지었다. 그들은 DESI, Euclid, LSST와 같은 대규모 서베이의 향후 분석에 ARF를 포함할 것을 권장하며, ARF가 적색편이 불확실성이 무시할 수 없는 수준인 광도 서베이에서 표준 3D 클러스터링에 대한 강력한 2D 대안이 될 수 있음을 강조하였다.