Fully-relativistic evolution of vacuum tensor inhomogeneities during inflation

문제 정의
표준적인 단일 장(single-field) 완만한 굴러가기(slow-roll) 인플레이션 모델은 CMB 위상 일관성과 가우시안성을 설명하는 데 성공적이지만, BiCEP 및 Planck와 같은 실험으로부터 엄격한 제약을 받고 있습니다. 이는 원시 확률적 중력파 배경(PGWB) 및 고차 통계량인 원시 이스펙트럼(bispectrum)을 예측하는 비표준 인플레이션 모델에 대한 관심을 다시 불러일으켰습니다. 역사적으로 이러한 스펙트럼에 대한 예측은 선형 운동 방정식이나 비선형 물질 역학을 해결하지만 선형 중력 섹터를 가정하는 격자 시뮬레이션에 의존해 왔습니다. 그러나 다중 장 인플레이션이나 강한 비선형 중력 결합(주요 차수를 넘어서는)이 포함된 시나리오에서는, 메트릭 섭동이 물질 섹터에 미치는 역작용(backreaction)과 그 반대의 과정이 관측 가능한 스펙트럼의 형태를 유의미하게 변화시킬 수 있습니다. 인플레이션 시기 자체를 시뮬레이션하면서 이러한 비선형 중력 효과를 통합하여 차세대 CMB 조사를 위한 정밀한 예측을 생성할 수 있는 견고하고 완전한 상대론적 도구가 부족한 실정입니다.

방법론
저자들은 중력 역작용을 포함하는 완전한 상대론적 수치 시뮬레이션 내에서 1차 인플레이션 텐서 섭동을 초기화하고 추출하기 위한 완전한 프레임워크를 제시합니다. 이 연구는 수치 상대론(Numerical Relativity, NR) 코드인 GRChombo(및 GRTeclyn으로의 포팅 진행 중)를 활용합니다.

주요 방법론적 구성 요소는 다음과 같습니다:

1. CPT-BSSN 대응 관계: 저자들은 1차 우주론적 섭동 이론(CPT) 변수들을 NR에서 사용되는 완전 상대론적 BSSN(Baumert-Shapiro-Shibata-Nakamura) 변수로 번역하는 게이지 불변(gauge-agnostic) 사전(dictionary)을 구축합니다. 이를 통해 선형 섭동 해를 비선형 시뮬레이션을 위한 초기 데이터로 변환할 수 있습니다.
2. 확률적 초기 조건: 격자 우주론의 표준적인 스칼라 초기화에서 영감을 얻은 확률적 텐서 초기 조건 생성 방법을 개발했습니다. 이 방법은 Mukhanov-Saasaki 방정식을 사용하여 Bunch-Davies 진공에서의 모드 함수를 도출합니다. 결정적으로, 저자들은 단순한 "비결맞음(decoherent)" 방식(크기와 위상만을 무작위화하는 방식)과 구별되는, 필드와 그 속도의 실수부와 허수부 사이의 $k$ 의존적 위상 관계를 보존하는 "모드 함수 방법"을 구현합니다.
3. 윈도우 함수: 작은 스케일에서의 큰 공간 미분으로 인한 수치적 불안정성을 완화하기 위해, 초기 데이터에 tanh 형태의 윈도우 함수를 적용하여 고차 모드로의 전력 누설(power leakage)을 억팡합니다.
4. 추출 알고리즘: 진화된 BSSN 변수로부터 가로-트레이스 없는(transverse-traceless, TT) 텐서 섭동을 추출하는 절차를 정의합니다. 섭동적 영역에서는 단순한 차분법($h_{ij} = \tilde{\gamma}_{ij} - \delta_{ij}$)이 사용되지만, 저자들은 향후 비섭동적 응용을 위해 트레이스 및 가로 성분을 분리하는 더 견고한 분해 방법을 개략적으로 설명합니다.

주요 기여

• 포멀리즘: 격자 계산과 섭동 방법 간의 비교를 용이하게 하는, CPT 변수와 BSSN 변수 사이의 완전한 게이지 불변 대응 관계를 유도하였습니다.
• 구현: 깊은 인플레이션 영역을 위해 설계된 특정 GRChombo 구현을 개발하였으며, 여기에는 Mukhanov-Saasaki 해를 통한 양자-고전 전이를 존중하는 확률적 텐서 초기 조건 생성이 포함됩니다.
• 검증: 선형 이론에 대한 엄격한 검증을 수행하였습니다. 저자들은 섭동 한계 내에서 배경의 완전 비선형 진화와 텐서 파워 스펙트럼이 선형ized Mukhanov-Saasaki 해와 일치함을 입증하였습니다.
• 제약 조건 충족: 에너지 및 운동량 제약 조건이 2차 또는 그 이상으로 만족됨을 확인하여, 시뮬레이션 내에서 1차 결과의 타당성을 확인하였습니다.

결과

• 배경 진화: 초지평선(super-horizon) 및 지평선 통과(horizon-crossing) 영역 모두에서, GRChombo에서 추출된 척도 인자($a$)와 허블 매개변수($H$)의 진화는 프라이드만 방정식과 높은 정밀도로 일치합니다.
• 파워 스펙트럼 회복: 시뮬레이션은 하부 지평선($1/k$)에서 상부 지평선($1/k^3$) 스케일링으로의 기대되는 스펙트럼 전환(turnover)을 성공적으로 회복합니다. "모드 함수" 초기화 방법은 필드 성분 간의 올바른 상대적 위상을 유지함으로써, "비결맞음" 방법보다 하부 지평선 영역에서의 매끄러운 전력 감쇄를 정확하게 재현한다는 점에서 우월함이 입증되었습니다.
• 비가우시안성: 왜도(skewness)와 첨도(kurtosis)를 통해 텐서 비가우시안성을 탐색하는 예비 결과를 제시하며, 모드가 지평선을 통과함에 따라 편광장(polarization fields)에서 구조가 동적으로 출현함을 보여줍니다.
• 수렴성: 수렴 테스트를 통해 제약 조건 위반이 억제되고 시간이 경과함에 따라 성장하지 않음을 확인하였으며, 이는 안정적인 수치 해임을 나타냅니다.

의의
본 논문은 이 계산 파이프라인이 메트릭과의 비선형 중력 결합 효과를 포함하는 원시 텐서 이스펙트럼 및 교차 스펙트럼의 출현을 연구하기 위한 필수적인 도구를 제공한다고 주장합니다. 선형 이론과 완전 상대론적 시뮬레이션 사이의 대응 관계를 검증함으로써, 저자들은 텐서-스칼라 비(tensor-to-scalar ratio)를 전례 없는 정밀도로 측정하고자 하는 차세대 CMB 조사(Simons Observatory 및 LiteBIRD 등)의 데이터를 분석하기 위한 도구로서의 토대를 마련하였습니다. 이 연구는 인플레이션의 섭동적 양자 처리와 완전 비선형 수치 상대론 사이의 간극을 메우며, 선형 이론이 불충분할 수 있는 시나리오를 연구할 수 있게 합니다.