[논문 리뷰] A thermal-kinetic subgrid model for supernova feedback in simulations of galaxy formation

Evgenii Chaikin, Joop Schaye외 4명|arXiv (Cornell University)|2022. 11. 09.

Astrophysics and Star Formation Studies인용 수 4

한 줄 요약

이 논문은 은하 형성 시뮬레이션에서 초신성 피드백을 위한 열역학적-운동량 하위그리드 모델을 제안한다. 고온의 열 에너지 주입과 저속도 운동량 킥을 조합하여 임의의 복사 손실을 억제하고 피드백 효율을 향상시킨다. 이 모델은 해상도 스케일 전반에서 뛰어난 수렴성을 보이며, 은하수 질량과 소형 은하에서 관측된 항성 형성률, 속도 분산 및 질량 적재 요소를 잘 재현한다.

ABSTRACT

We present a subgrid model for supernova feedback designed for cosmological simulations of galaxy formation that may include a cold interstellar medium (ISM). The model uses thermal and kinetic channels of energy injection, which are built upon the stochastic kinetic and thermal models for stellar feedback used in the OWLS and EAGLE simulations, respectively. In the thermal channel, the energy is distributed statistically isotropically and injected stochastically in large amounts per event, which minimizes spurious radiative energy losses. In the kinetic channel, we inject the energy in small portions by kicking gas particles in pairs in opposite directions. The implementation of kinetic feedback is designed to conserve energy, linear and angular momentum, and is statistically isotropic. To test the model, we run simulations of isolated Milky Way-mass and dwarf galaxies, in which the gas is allowed to cool down to 10 K. Using the thermal and kinetic channels together, we obtain smooth star formation histories and powerful galactic winds with realistic mass loading factors. Furthermore, the model produces spatially resolved star formation rates (SFRs) and velocity dispersions that are in agreement with observations. We vary the numerical resolution by several orders of magnitude and find excellent convergence of the global SFRs and wind mass loading. We show that large thermal-energy injections generate a hot phase of the ISM and modulate the star formation by ejecting gas from the disc, while the low-energy kicks increase the turbulent velocity dispersion in the neutral ISM, which in turn helps suppress star formation.

연구 동기 및 목표

• 제한된 수치 해상도에도 관측된 은하 성질—특히 항성 형성률과 은하 외부 유동—을 정확히 재현할 수 있는 하위그리드 피드백 모델을 개발하는 것.
• 저해상도 시뮬레이션에서 과도한 복사 냉각 손실로 인해 효율이 떨어지는 순수 열 피드백 모델의 문제를 해결하는 것.
• 특히 저속도 킥을 통한 운동량 피드백이 중성 가스 은하간성 매질에서 난류 속도 분산을 증가시켜 항성 형성 억제에 기여하는 방식을 탐색하는 것.
• 다양한 해상도 수준에서 핵심 은하 관측량의 수치 수렴성을 확보하여 천체역학적 시뮬레이션의 강건성을 확보하는 것.

제안 방법

• 복사 에너지 손실을 최소화하기 위해 큰 크기의 등방성 폭발 형태로 에너지를 무작위로 주입하는 열 피드백 채널을 구현하며, 지연 냉각 기법을 사용한다.
• 짝을 이룬 가스 입자에 반대 방향으로 킥을 가하여 운동량을 주입하는 운동량 피드백 채널을 도입하며, 에너지, 선형 운동량, 각운동량을 모두 보존한다.
• 모든 채널에 통계적 등방성 프레임워크를 적용하여 물리적 현실성 확보 및 피드백 효과에 대한 방향성 편향 방지.
• 운동량 주입 비율 $ f_{ m kin} $ 와 특성 킥 속도 $ riangle v_{ m kick} $ 를 캘리브레이션하며, 기준 값으로 $ f_{ m kin} = 0.1 $ 과 $ riangle v_{ m kick} = 50 $ km s$^{-1} $ 을 사용한다.
• 열역학적 냉각이 10 K까지 가능하도록 한 은하수 질량과 소형 은하의 고립 시뮬레이션을 실행하여 피드백 효율성과 수렴성 테스트.
• 관측 데이터와 비교하기 위해 공간적으로 해상도가 높은 항성 형성률, H i 속도 분산 및 외부 유동 질량 적재 요소를 분석한다.

실험 결과

연구 질문

• RQ1 다양한 해상도에서 고립된 은하에서 관측된 케니쿠트-스미스 항성 형성 관계를 열역학적-운동량 피드백 모델이 재현할 수 있는가?
• RQ2 열 에너지와 운동량 에너지 주입 간의 균형이 중성 은하간성 매질에서 난류 속도 분산과 항성 형성 억제에 어떤 영향을 미치는가?
• RQ3 항성 형성률과 외부 유동 질량 적재 요소의 해상도 스케일 전반에서 수렴성을 달성하기 위해 $ f_{ m kin} $ 와 $ riangle v_{ m kick} $ 의 최적 조합은 무엇인가?
• RQ4 저에너지 운동량 킥이 난류 속도 분산을 증가시키고 국소 항성 형성 효율을 감소시키지만 전체 은하 성질을 파괴하지는 않으며, 그 정도는 어느 정도인가?
• RQ5 열 피드백 채널과 운동량 피드백 채널이 서로 어떻게 보완되어 현실적인 은하 외부 유동을 유지하고 밀도가 높은 가스에서 과도한 냉각을 억제하는가?

주요 결과

• 열역학적-운동량 모델은 매끄러운 항성 형성 역사와 관측 결과와 일치하는 질량 적재 요소를 가진 현실적인 은하 외부 유동을 생성한다.
• 이 모델은 수치 해상도가 수십 개의 주요 스케일에 걸쳐 전반적으로 뛰어난 수렴성을 보이며, 전역 항성 형성률과 외부 유동 질량 적재 요소 모두에서 수렴성이 뛰어나다.
• 관측 결과의 기울기를 정확히 재현하기 위해 작은 비율이지만 비제로인 운동량 비율($ f_{ m kin} o 0.1 $)이 필수적이다; 순수 열 또는 순수 운동량 모델은 관측 결과와 일치하지 않는다.
• 저에너지 킥($ riangle v_{ m kick} = 50 $ km s$^{-1} $)은 중성 은하간성 매질의 난류 속도 분산을 몇 배 증가시켜 국소 항성 형성 억제를 유도하지만 총 항성 형성률은 감소시키지 않는다.
• $ riangle v_{ m kick} o 10^2 $ km s$^{-1} $ 일 때, 은하 성질인 항성 형성률과 외부 유동 질량 적재 요소는 킥 속도의 추가 감소에 민감하지 않게 되어 수렴성이 확인된다.
• 모델에서 유도된 공간적으로 해상도가 높은 H i 속도 분산은 Zhou 등 (2017) 및 Law 등 (2022)의 관측 결과와 잘 일치하며, 특히 $ f_{ m kin} = 0.1 $ 과 $ 0.3 $ 일 때 뚜렷한 일치를 보인다.