[논문 리뷰] Simulation of Primordial Black Holes with large negative non-Gaussianity

Albert Escrivà, Yuichiro Tada외 2명|arXiv (Cornell University)|2022. 02. 02.

Cosmology and Gravitation Theories참고 문헌 90인용 수 73

한 줄 요약

이 논문은 국소형 비정규성(fNL가 음수로 파arameterized된)의 레이디에이션 지배 우주에서 초기 블랙홀(PBH) 형성에 대해 고정밀 수치 시뮬레이션을 수행한다. 이전의 분석적 추정치에서 사용된 보편적인 평균 임계 압축 함수(¯Cc = 2/5)와는 달리, 시뮬레이션은 fNL ≲ −0.336일 때조차도 유형 I PBH 형성이 가능함을 확인한다. 이는 특정 비정규성 프로파일에 대해 이 보편적 임계값이 실패함을 시사한다. 대신, q-파aram터 피팅 공식은 fNL ≳ −1인 범위에서 오차가 2% 미만인 더 견고한 분석적 추정치를 제공한다.

ABSTRACT

In this work, we have performed numerical simulations of primordial black hole (PBH) formation in the Friedman-Lema\^itre-Robertson-Walker universe filled by radiation fluid, introducing the local-type non-Gaussianity to the primordial curvature fluctuation. We have compared the numerical results from simulations with previous analytical estimations on the threshold value for PBH formation done in the previous paper arXiv:2109.00791, particularly for negative values of the non-linearity parameter $f_{ m NL}$. Our numerical results show the existence of PBH formation of (the so-called) type I also in the case $f_{ m NL} \lesssim -0.336$, which was not found in the previous analytical expectations using the critical averaged compaction function. In particular, although the universal value for the averaged critical compaction function $\bar{\mathcal{C}}_{c}=2/5$ found previously in the literature is not satisfied for all the profiles considered in this work, an alternative direct analytical estimate has been found to be roughly accurate to estimate the thresholds, which gives the value of the critical averaged density with a few $\%$ deviation from the numerical one for $f_{ m NL}\gtrsim -1$.

연구 동기 및 목표

• 큰 음수 국소형 비정규성(fNL)이 존재할 경우, 보편적인 평균 임계 압축 함수(¯Cc = 2/5)가 PBH 형성에 대해 유효한지 테스트하기 위해.
• 분석적 추정치가 형성이 불가능하다고 예측하는 경우(즉, fNL ≲ −0.336일 경우) 유형 I PBH가 실제로 형성되는지 조사하기 위해.
• 강한 비정규성 하에서 보편적인 ¯Cc = 2/5 가정에 의존하지 않는 대안 분석적 임계값 추정 수단으로서 q-파aram터 피팅 공식의 견고함을 평가하기 위해.
• 음수 fNL에 따른 PBH 질량 함수 및 비율 예측을, 새로 규명된 PBH 형성 가능 영역을 반영하여 업데이트하기 위해.

제안 방법

• 기존 연구들과 비교해 훨씬 향상된 코드를 사용하여, 복사로 지배되는 프리드만–레마트르–로버르손–워커(FLRW) 우주에서 PBH 형성의 수치 시뮬레이션을 수행한다.
• 지역화된 곡률 편미분을 모델링하기 위해 단일 파수 k*를 가진 단색 스펙트럼을 도입한다.
• fNL 파aram터를 통해 국소형 비정규성을 구현하며, 초기 곡률 편미분을 ζ(r) = µ sinc(k*r) + (3/5)fNL µ² sinc²(k*r)로 수정한다.
• 보편적인 ¯Cc = 2/5 가정에 의존하지 않는 수식적 임계값 추정을 위해, q-파aram터 피팅 공식(Eq. 3.8)을 사용한다.
• 다양한 fNL 값에 대해 수치적 PBH 형성 임계값을 평균 압축 함수 접근법과 q-파aram터 피팅 공식 간에 비교한다.
• 피크 이론을 통한 통계적 유도를 통해 PBH 질량 함수를 도출하며, 임계 행동 M ∝ (µ − µc)^γ 및 γ ≈ 0.36을 고려하고 현재 암흑물질 밀도에 정규화한다.

실험 결과

연구 질문

• RQ1 보편적인 평균 임계 압축 함수 ¯Cc = 2/5에 기반한 분석적 예측과 달리, fNL ≲ −0.336일 경우에도 PBH 형성(특히 유형 I)이 발생하는가?
• RQ2 큰 음수 fNL 값에서 PBH 형성 임계값을 추정하는 데 있어 q-파aram터 피팅 공식(δc(q))의 정확도는 어느 정도인가?
• RQ3 비정규성이 PBH 질량 함수 및 총 비율에 미치는 영향은 무엇이며, 특히 fNL ≲ −0.336 영역에서 새로 접근 가능한 영역에서 어떤가?
• RQ4 보편적인 ¯Cc = 2/5 가정의 실패는 특정 프로파일에 의존하는 비정규성 곡률 변동의 행동 탓인가?

주요 결과

• 수치 시뮬레이션은 fNL ≲ −0.336인 경우에도 유형 I PBH 형성이 발생함을 확인하였으며, 이는 보편적인 ¯Cc = 2/5 기준에 기반한 분석적 예측와 정반대이다.
• 평균 임계 압축 함수 ¯Cc는 모든 비정규성 프로파일에 대해 보편적으로 2/5가 아니며, fNL ≲ −0.336일 경우에 크게 이격됨을 확인하여 보편성 가정의 붕괴를 시사한다.
• q-파aram터 피팅 공식(Eq. 3.8)은 수치 결과와 2% 미만의 오차로 일치하는 견고한 분석적 추정치를 제공하며, fNL ≳ −1 범위에서 유용하다.
• fNL ≲ −1인 경우, fNL 전개가 신뢰할 수 없게 되며, q-파aram터 피팅의 정확도도 점차 떨어지므로 그 적용 범위의 한계가 있음을 시사한다.
• PBH 질량 함수는 fNL ≲ −0.336 영역을 포함하도록 업데이트되었으며, fNL = −1일 때 σ²₀를 적절히 조정하면 총 PBH 비율 ftot_PBH가 1에 도달한다.
• 임계점 µc가 경계 (3/5)µfNL = −1/2와 만날 때 fNL ≈ −1.01에서 발생하며, 이는 형성 제도의 전환점으로서의 의미를 가진다.