UV 광도 함수 변동성: SPICE 시뮬레이션의 초신성 피드백 효과

초록

본 연구는 SPICE 방사선‑수소역학 시뮬레이션에 구현된 세 가지 초신성(SN) 피드백 모델(bursty‑sn, smooth‑sn, hyper‑sn)이 z > 5 시대의 UV 광도 함수(UVLF) 변동성에 미치는 영향을 정량화한다. bursty‑sn 모델은 강하고 간헐적인 SN 폭발로 SFR 변동성을 가장 크게 만들며, UVLF 변동도 최대 2.5 까지 상승한다. smooth‑sn은 완만한 피드백으로 변동성을 최소화하고, hyper‑sn은 중간 정도의 변동성을 보이며 낮은 z 에서는 smooth‑sn과 유사한 경향을 나타낸다. 변동성은 주로 저질량 하얼로에서 크게 나타나며, 모든 모델에서 z 가 낮아질수록 감소한다. 결과는 SN 피드백이 고‑z 시대 밝은 은하 과잉 문제를 완화할 수 있음을 시사한다.

상세 분석

본 논문은 최신 방사선‑수소역학 코드 RAMSES‑RT 기반의 SPICE 시뮬레이션을 활용해 초신성 피드백이 고‑z 은하의 UVLF 변동성에 미치는 물리적 메커니즘을 정밀히 탐구한다. 시뮬레이션은 10 h⁻¹ cMpc 부피에 512³개의 DM 입자를 배치해 질량 해상도 m_dm ≈ 6.4 × 10⁵ M⊙을 달성했으며, AMR 레벨 ℓ_max = 16을 통해 z = 5 시점에 28 pc까지의 공간 해상도를 제공한다. 별 형성은 Kretschmer & Teyssier(2020)의 난류‑기반 서브그리드 모델을 적용해 지역 가스의 Mach 수와 virial 파라미터에 따라 가변적인 효율을 부여한다.

세 가지 피드백 모델은 SN 폭발 시점과 에너지 주입 방식을 차별화한다. bursty‑sn은 모든 별 입자가 10 Myr 후에 단일 SN 이벤트를 겪으며, 에너지는 2 × 10⁵¹ erg로 고정한다. smooth‑sn은 질량에 따라 3–40 Myr 사이에 다중 SN를 발생시키고 동일한 에너지를 주입한다. hyper‑sn은 smooth‑sn과 동일한 타이밍을 유지하되, 에너지를 정규분포(평균 1.2 × 10⁵¹ erg, 범위 10⁵⁰–2 × 10⁵¹ erg)에서 추출하고, 금속도에 비례해 HN(10⁵² erg) 비율을 추가한다. 이러한 차이는 각각 별 형성률(SFR)의 시간적 변동성에 직접적인 영향을 미친다.

시뮬레이션 결과, bursty‑sn은 강력하고 간헐적인 에너지 투입으로 가스 흐름을 급격히 교란시켜 SFR이 높은 진폭의 주기적 상승·하강을 보인다. 이는 M_UV–M_halo 관계에서 큰 산란을 일으키며, UVLF 변동성을 2.5까지 끌어올린다. 반면 smooth‑sn은 연속적인 에너지 주입으로 가스가 비교적 안정적인 흐름을 유지해 SFR 변동폭이 작으며, UVLF 변동성은 1.3 수준에 머문다. hyper‑sn은 에너지 분포가 넓어 중간 정도의 변동성을 보이지만, 낮은 질량 하얼로(10⁹–10¹⁰ M⊙)에서는 여전히 bursty‑sn에 근접한 변동을 나타낸다.

특히, 변동성은 하얼로 질량에 강하게 의존한다. 저질량 하얼로는 얕은 중력우물 때문에 SN 피드백이 가스를 효과적으로 탈출시켜 SFR이 급격히 변동하고, 이에 따라 UVLF의 밝기 구간에서도 큰 스캐터가 발생한다. 반대로 고질량 하얼로(>10¹¹ M⊙)에서는 중력에 의해 피드백 효과가 억제돼 변동성이 감소한다. 시간에 따라, 모든 모델에서 z 가 낮아질수록(즉, 우주가 나이가 들수록) 변동성은 점진적으로 감소한다. 이는 가스 공급이 점차 안정화되고, 별 형성 효율이 평균화되기 때문이다.

UVLF 자체는 BPASSv2.2.1 SED 모델을 이용해 1500 Å에서 직접 계산했으며, RASCAS Monte‑Carlo 방사선 전달 코드를 통해 먼지 감쇠를 적용했다. 결과적으로, M_UV ≈ −15 이하에서는 먼지 감쇠가 눈에 띄게 나타나며, 이는 피드백 모델에 따라 감쇠 시작점이 약간씩 달라진다. 시뮬레이션의 UVLF는 관측 데이터(HST, JWST)와 전반적으로 일치하지만, hyper‑sn 모델은 가장 낮은 밝기에서 관측보다 약간 낮은 수치를 보인다.

결론적으로, 초신성 피드백의 강도와 타이밍이 고‑z 시대 UVLF 변동성의 주된 원천임을 확인했으며, 특히 bursty‑sn과 같은 강력한 피드백이 밝은 은하 과잉 문제를 완화할 수 있는 메커니즘을 제공한다. 이는 JWST가 보고한 z > 10 밝은 은하들의 과잉을 설명하는 데 중요한 물리적 인자를 제시한다.