cloelike: A Python Library for Cosmological Likelihood Inference in the Euclid Era

cloelike: Euclid 시대의 우주론적 가능도 추론을 위한 Python 라이브러리 기술 요약

문제 제기
ESA Euclid 임무는 고정밀 광학 및 분광 관측을 통해 암흑 에너지, 암흑 물질, 그리고 초기 밀도 요동을 탐구하는 것을 목표로 합니다. 이러한 목표를 달성하기 위해서는 통계적 추론을 위한 견고하고 검증되었으며 재현 가능한 가능도 프레임워크가 필요합니다. 우주론 커뮤니티는 CosmoSIS, Cobaya 와 같은 광범위한 추론 도구를 보유하고 있지만, 기존 공개 가능도 구현체들은 두 가지 한계를 자주 겪습니다. 즉, 특정 관측 대상에 국한되거나, 완전 통합형 종단간 (end-to-end) Euclid 분석에 필요한 모듈성과 네이티브 데이터 형식 호환성이 부족하다는 점입니다. 따라서 외부 베이지안 프레임워크와 원활하게 통합되면서 엄격하게 Euclid 데이터 표준을 준수하는 $\chi^2$ 통계량 계산을 위한 통일되고 경량화된 인터페이스가 필요합니다.

방법론 및 소프트웨어 설계
cloelike는 가우시안 가능도 추론을 위한 모듈식이고 구성 가능한 라이브러리를 제공함으로써 이러한 한계를 해결하도록 설계된 Python 패키지입니다. 그 아키텍처는 엄격한 관심사 분리 원칙에 기반합니다:

1. 핵심 기능: 이 패키지는 표준 다변량 가우시안 로그 가능도 평가에 exclusively 집중합니다:
   $$-2 \ln \mathcal{L} = (\mathbf{d} - \mathbf{t})^T \mathbf{C}^{-1} (\mathbf{d} - \mathbf{t}) + b$$
   여기서 $\mathbf{d}$는 데이터 벡터, $\mathbf{t}$는 이론적 예측, $\mathbf{C}$는 공분산 행렬, $b$는 가산 편향입니다. 이 패키지는 매개변수 샘플링이나 모델 관리를 수행하지 않으며, 이러한 작업은 Cobaya, CosmoSIS, Nautilus 와 같은 외부 엔진에 위임합니다.

2. 모듈식 아키텍처:

   • Base Classes: 광학 관측 대상의 경우, 공유 베이스 클래스가 데이터 벡터 수용, 마스킹, 스케일 컷, 공분산 로딩, 정밀 행렬 처리, 일관성 검사를 처리합니다.
   • Mixins: 이론적 관측량 계산은 독립적이고 관측 대상별인 믹스인 (예: 조화 공간 각 파워 스펙트럼, 구성 공간 상관 함수, COSEBIs) 으로 위임됩니다. 이 설계는 코드 중복을 피하고, 이론적 예측 모듈만 구현함으로써 새로운 관측량을 추가할 수 있게 합니다.
   • 프로토콜 기반 통합: cloelike는 Python 프로토콜을 사용하여 cloelib 이론적 예측 라이브러리와 직접 인터페이스합니다. 이는 가능도 정의와 이론 계산이 분리되도록 보장하여, cloelib 내의 비선형 모델링 (예: 은하 편향, 적색편이 공간 왜곡) 업데이트로부터 cloelike가 즉시 혜택을 받을 수 있게 합니다.

3. 데이터 호환성: 이 라이브러리는 euclidlib 데이터 형식 내에서 구조화된 데이터 벡터, 공분산 행렬, 그리고 이론적 예측을 기반으로 작동합니다. 이는 다음을 지원합니다:

   • 광학 관측 대상: 약한 렌즈 (WL), 광학 은하 군집 (GCph), 은하 - 은하 렌즈 (GGL) 에 대한 조화 공간 각 파워 스펙트럼 및 실공간 2 점 상관 함수를 모든 결합 조합 (3$\times$2pt, 2$\times$2pt) 에서 다룹니다.
   • 분광 관측 대상: 전체 모양 (full-shape) 파워 스펙트럼 다중극 및 중입자 음향 진동 (BAO).

주요 기여

• 통일된 인터페이스: cloelike는 모든 주요 Euclid 관측 대상에 걸쳐 광학 및 분광 데이터를 모두 지원하는 $\chi^2$ 통계량 계산을 위한 최소한의 통일된 Python 인터페이스를 제공합니다.
• 구성 가능성: 이 프레임워크는 사용자들이 기본 인프라를 수정하지 않고도 데이터 컨테이너, 공분산 표현, 관측량 제공자를 결합하여 공동 가능도 (예: WL 과 GCph 결합) 를 구성할 수 있게 합니다.
• 성능 및 확장성: 구현은 벡터화 연산을 활용하고, 중간 양 (예: 역공분산 행렬) 을 캐싱하며, 고차원 매개변수 공간에서의 효율성을 보장하기 위해 병렬 샘플링 코드와 호환됩니다.
• 커뮤니티 채택: 이 소프트웨어는 이전 CLOE 소프트웨어의 자연스러운 진화로, Euclid 컨소시엄 내에서 참조 가능도 구현체로 활발히 사용되고 있습니다.

결과 및 활용
본 논문은 두 가지 대표적인 가능도에 대한 활용 예시를 통해 라이브러리의 유용성을 보여줍니다:

1. 조화 광학 가능도: 우주 전단 각 파워 스펙트럼에 대한 로그 가능도의 초기화 및 평가를 보여주는 예시로, 데이터 사전 구축, 스케일 컷 설정, 그리고 특정 cloelib 배경 및 요동 클래스를 사용한 가능도 인스턴스화를 상세히 설명합니다.
2. 분광 전체 모양 가능도: 분광 은하 군집 전체 모양 파워 스펙트럼에 대한 로그 가능도 계산을 보여주는 예시로, 다중극 ($P_0, P_2, P_4$) 처리, 기준 우주론 추출, 단위 재조정, 그리고 적색편이 구간 전반에 걸친 유해 매개변수 (예: 은하 편향, EFT 매개변수) 정의를 포함합니다.

의의
저자들은 cloelike가 커뮤니티가 공식 Euclid 결과를 재현하는 장벽을 낮추면서도 새로운 체계적 오차와 확장된 우주론 모델에 대한 유연한 탐색을 가능하게 한다고 주장합니다. 기존 종단간 추론 프레임워크를 보완하는 전용 검증 가능도 레이어로 작용함으로써 기능 중복 없이 상호 운용성을 보장합니다. 그 모듈식 설계는 우주 마이크로파 배경과의 교차 상관과 같은 새로운 관측량의 포함을 용이하게 하며, Euclid 관측 데이터의 전체 과학적 잠재력을 추출하는 데 필요한 일관된 다중 관측 대상 분석을 지원합니다. 이 패키지는 커뮤니티 검증과 재현성을 지원하기 위해 MIT 라이선스 하에 공개적으로 출시되었습니다.