[논문 리뷰] The halo mass function from the excursion set method. II. The diffusing barrier

Michele Maggiore, Antonio Riotto|arXiv (Cornell University)|2009. 03. 06.

Galaxies: Formation, Evolution, Phenomena인용 수 4

한 줄 요약

이 논문은 충돌 임계값을 고정된 값이 아니라 확률적이고 확산되는 장벽으로 모델링함으로써 추적 집합 이론을 확장하여, 헬로 질량 함수의 정확도를 향상시킨다. 결과적으로 유도된 질량 함수는 exp{−aδ_c²/(2σ²)} 형태의 수정된 지수 인자로 구성되며, 여기서 a = 1/(1 + D_B)이며, D_B는 장벽의 확산 계수이다. 이는 현실적인 헬로 형성 조건에서 N-체 시뮬레이션 결과와 더 잘 일치한다.

ABSTRACT

In excursion set theory the computation of the halo mass function is mapped into a first-passage time process in the presence of a barrier, which in the spherical collapse model is a constant and in the ellipsoidal collapse model is a fixed function of the variance of the smoothed density field. However, N-body simulations show that dark matter halos grow through a mixture of smooth accretion, violent encounters and fragmentations, and modeling halo collapse as spherical, or even as ellipsoidal, is a significant oversimplification. We propose that some of the physical complications inherent to a realistic description of halo formation can be included in the excursion set theory framework, at least at an effective level, by taking into account that the critical value for collapse is not a fixed constant $\delta_c$, as in the spherical collapse model, nor a fixed function of the variance $\sigma$ of the smoothed density field, as in the ellipsoidal collapse model, but rather is itself a stochastic variable, whose scatter reflects a number of complicated aspects of the underlying dynamics. Solving the first-passage time problem in the presence of a diffusing barrier we find that the exponential factor in the Press-Schechter mass function changes from $\exp\{-\delta_c^2/2\sigma^2\}$ to $\exp\{-a\delta_c^2/2\sigma^2\}$, where $a=1/(1+D_B)$ and $D_B$ is the diffusion coefficient of the barrier. The numerical value of $D_B$, and therefore the corresponding value of $a$, depends among other things on the algorithm used for identifying halos. We discuss the physical origin of the stochasticity of the barrier and we compare with the mass function found in N-body simulations, for the same halo definition.[Abridged]

연구 동기 및 목표

• 표준 추적 집합 이론에서 고정된 충돌 임계값의 과도한 단순화를 해결하기 위해.
• 헬로 형성의 물리적 복잡성—예를 들어 부드러운 첨가 및 폭력적인 이완—을 이론적 프레임워크에 통합하기 위해.
• 헬로 붕괴의 동적이고 비구형 성격을 더 잘 반영하는 확률적 장벽 모델을 개발하기 위해.
• N-체 시뮬레이션 결과와 더 밀접게 일치하는 수정된 헬로 질량 함수를 유도하기 위해.

제안 방법

• 헬로 형성의 역동적 불확실성을 반영하기 위해 확산 계수 D_B를 갖는 확률적 변수로 임계 충돌 임계값 δ_c를 모델링한다.
• 추적 집합 형식에서 시간에 따라 변화하는 확산 장벽을 가진 첫 번째 통과 시간 문제로 헬로 질량 함수를 설정한다.
• 확산 장벽을 가정할 때 첫 번째 통과 시간 문제를 해석적으로 해결하여 수정된 지수 감쇠 인자를 도출한다.
• Press-Schechter 유사 질량 함수의 지수 항에서 효과적 지수 a = 1/(1 + D_B)를 유도한다.
• 확산 계수 D_B를 헬로 식별 알고리즘과 첨가 동역학과 같은 물리적 과정과 연결한다.
• 동일한 헬로 정의를 사용하여 유도된 질량 함수를 N-체 시뮬레이션과 비교함으로써 모델을 검증한다.

실험 결과

연구 질문

• RQ1 충돌 임계값을 확률적이고 확산되는 장벽으로 모델링할 경우, 고정된 임계값 모델 대비 헬로 질량 함수 예측에 어떤 개선이 이루어지는가?
• RQ2 충돌 임계값의 확률성의 물리적 기원은 무엇이며, 현실의 헬로 형성 과정을 어떻게 반영하는가?
• RQ3 장벽의 확산 계수 D_B가 헬로 질량 함수의 형태에 어떤 영향을 미치는가?
• RQ4 a = 1/(1 + D_B)를 갖는 수정된 질량 함수는 N-체 시뮬레이션 결과와 어느 정도 일치하는가?
• RQ5 헬로 찾기 알고리즘의 선택은 효과적 확산 계수 D_B와 결과로 유도된 질량 함수에 어떤 영향을 미치는가?

주요 결과

• 헬로 질량 함수의 지수 인자는 exp{−δ_c²/(2σ²)}에서 exp{−aδ_c²/(2σ²)}로 수정되었으며, 여기서 a = 1/(1 + D_B)이며, D_B는 장벽의 확산을 정량화한다.
• D_B의 값은 헬로 형성의 물리적 특성과 헬로 찾기 알고리즘에 따라 달라지며, 이는 다양한 시뮬레이션 설정에 적응 가능한 모델임을 의미한다.
• 확률적 장벽 모델은 비구형이고 일정하지 않은 붕괴 역학을 반영함으로써 N-체 시뮬레이션 결과와의 일치도를 향상시킨다.
• 효과적 지수 a < 1은 표준 Press-Schechter 형태에 비해 질량 함수의 尾부가 덜 감쇠됨을 의미한다.
• 이 모델은 질량 함수가 구형 붕괴 예측에서 관측된 편차를 물리적으로 설명할 수 있는 근거를 제공한다.
• 이 프레임워크는 분석적 취급이 가능함을 유지하면서도 추적 집합 형식 내에서 복잡한 헬로 성장 과정을 체계적으로 통합할 수 있다.