BeeTheory – 기초 – 기술 노트 XII

공식화:
은하계 규모의 벌 이론 계산

이 노트는 원반 은하에 적용될 때 BeeTheory 프레임워크를 공식화합니다. 여기에는 관측 입력, 바이리오닉 분포의 기하학적 분해, 각 구성 요소의 파장을 정의하는 적분 방정식, 예측된 회전 곡선을 산출하는 일련의 연산이 명시되어 있습니다. 이 절차는 엄격하게 단방향으로 진행되며, 관찰된 바이리오닉 구조가 파장을 결정하고 회전 곡선을 결정하는 것이지 그 반대의 경우는 절대 없습니다.

1. 하나의 다이어그램에서 계산

단방향 체인

관측된 광도 $\;\longrightarrow\;$ 바리오닉 분해 $(\rho_\text{bar})$
$\big\downarrow$
파장-장 컨볼루션 $\;\longrightarrow\;$ 파동 밀도 $(\rho_\text{wave})$
$\big\downarrow$
질량 적분 $\;\longrightarrow\;$ 밀폐된 파동 질량 $(M_\text{wave})$
$\big\downarrow$
뉴턴 관계 $\;\longrightarrow\;$ 예측된 회전 곡선 $(V_c)$

단계가 반전되지 않습니다. 회전 곡선 $V_c(R)$ 은 입력으로 사용되지 않습니다.

2. 관찰 입력

각 은하에 대해 계산하려면 공개된 5개의 관측값이 필요합니다. 이것들은 유일한 은하별 수량이며, 다른 모든 수치는 이로부터 계산됩니다. 이 단계에서는 회전 곡선에 대한 피팅은 수행되지 않습니다.

기호	수량	출처
$T$	허블 형태 유형	카탈로그(드 보쿨레르 외. 1991, SPARC)
$R_d$	스텔라 디스크 스케일 길이(kpc)	스피처 3.6µm 광도 측정(SPARC)
시그마_d$	중앙 디스크 표면 밝기($L_\odot/\text{pc}^2$)	스피처 3.6µm 광도 측정(SPARC)
$M_\text{HI}$	총 원자 수소 질량($M_odot$)	21cm 전파 관측(SPARC)
$\업실론_\별$	3.6µm의 뛰어난 질량 대 빛 비율	유니버설 고정: $0.5\,M_\odot/L_\odot$ (McGaugh 2014)

은하수의 경우, $R_d$, $Sigma_d$, $M_text{HI}$는 내부 항성 조사(Bovy & Rix 2013) 및 21cm 지도로부터 결정된 유사 값으로 대체됩니다. 동일한 5량 입력 벡터가 사용됩니다.

3. 바리오닉 분해 – 5가지 기하학적 구성 요소

5개의 관측 입력에서 바이리오닉 질량은 5개의 서로 다른 기하학적 구성 요소로 분할됩니다. 각 구성 요소에는 고유한 밀도 프로파일과 특징적인 스케일이 있습니다.

3.1 총 항성 및 가스 질량

$$M_\star \;=\; 2\pi\,R_d^2\,\Sigma_d\,\Upsilon_\star$$

$$M_\text{gas} \;=\; 1.33\,M_\text{HI} \q쿼드 \text{(He 보정; Arnett 1996)}$$

3.2 구성 요소 질량 및 스케일

구성 요소	질량	규모	활성화
벌지	$M_b = 0.20\,M_\star$	r_b = \max(0.5\,R_d,\,0.3\text{ kpc})$.	만약 $T \leq 4$라면
씬 디스크	$M_\text{thin} = 0.75\,(M_\star – M_b)$	$R_d$	항상
두꺼운 디스크	$M_\text{두께} = 0.25\,(M_\star – M_b)$	1.5\,R_d$	항상
가스 링	$M_\text{gas} = 1.33\,M_\text{HI}$	$R_g = 1.7\,R_d$ (Broeils & Rhee 1997)	항상
나선형 암	$M_\text{arm} = 0.10\,M_\text{thin}$ (유효)	R_d$ (씬 디스크 뒤따름)	항상

3.3 밀도 프로파일

벌지(3D 헤른퀴스트)

$$\rho_b(r) \;=\; \frac{M_b\,r_b}{2\pi\,r\,(r + r_b)^3}$$

얇고 두꺼운 항성 원반(2D 지수)

$$\Sigma_\text{thin}(R) \;=\; \frac{M_\text{thin}}{2\pi\,R_d^2}\,e^{-R/R_d}$$

$$\Sigma_\text{thick}(R) \;=\; \frac{M_\text{thick}}{2\pi\,(1.5\,R_d)^2}\,e^{-R/(1.5R_d)}$$

가스 링(중앙 구멍이 있는 2D 지수)

$$\Sigma_\text{gas}(R) \;=\; \frac{M_\text{gas}}{2\pi\,R_g^2}\,\exp\!\left(-\frac{R_\text{hole}}{R} – \frac{R}{R_g}\right), \quad R_\text{hole} = 0.5\,R_g$$.

나선형 암 초과(2D, 얇은 디스크 뒤따름)

$$\Sigma_\text{arm}(R) \;=\; 0.10\;\Sigma_\text{thin}(R)$$

4. 웨이브 커널

각 바이리오닉 질량 원소는 비이론 파장을 생성합니다. 소스 원소에 의해 $vec{r},’$에서 $D = |vec{r} – vec{r},’|$로 구분된 한 점 $vec{r}$에서 생성되는 필드는 주 I의 정규화된 파동 함수에서 파생된 유카와형 커널의 지배를 받습니다:

비이론 파동 커널

$$\mathcal{K}_i(D) \;=\; K_0\,\frac{(1 + \alpha_i\,D)\,e^{-\alpha_i\,D}}{D^2}, \qquad \alpha_i \;=\; \frac{1}{\ell_i}$$.

여기서 $\K_0$은 보편적인 파동 질량 진폭(단일 차원이 없는 수)이고 $\ell_i$는 성분 $i$의 코히어런스 길이입니다. 커널은 짧은 간격에서 준 뉴턴식 $1/D^2$ 동작을 인코딩하며, $\ell_i$ 이상의 스케일에서는 지수 컷오프에 의해 변조됩니다. (1 + \알파 D)\,e^{-\알파 D}$ 형식은 무한대에서의 연속성과 유한한 총 밀폐 질량을 보장합니다.

4.1 컴포넌트 일관성 길이

각 컴포넌트의 일관성 길이는 자연스러운 기하학적 눈금에 해당 차원에 특정한 무차원 상수를 곱하여 설정됩니다:

구성 요소	일관성 길이	기하 상수
벌지(3D 구체)	$\ell_b = c_\text{sph}\,r_b$	$c_\TEXT{SPH}$
씬 디스크(2D)	$\ell_\text{thin} = c_\text{디스크}\,R_d$	$c_\text{디스크}$
두꺼운 디스크(2D)	$\ell_\text{두께} = c_\text{디스크}\,(1.5\,R_d)$.	$c_\text{디스크}$
가스 링(2D)	$\ell_\text{gas} = c_\text{디스크}\,R_g$	$c_\text{디스크}$
나선형 암(2D, 방위각으로 집중)	$\ell_\text{arm} = c_\text{arm}\,R_d$	$c_\text{arm}$

세 가지 기하학적 상수 $(c_\text{sph},\,c_\text{디스크},\,c_\text{arm})$는 은하마다 다르지 않은 보편적인 상수입니다. 전역 파동 질량 진폭 $K_0$ 및 파동-장 결합 $\lambda$와 함께 이론 수준의 전체 매개변수 집합을 구성합니다.

5. 파장 컨볼루션 – 구성 요소별 적분 방정식

위치 $\vec{r}$에서의 파장 밀도는 파동 커널을 가진 바이리오닉 소스 분포의 컨볼루션입니다. 은하 대칭 시스템(축 대칭, 단극 근사)의 경우 각 바이리오닉 구성 요소는 추가적으로 기여합니다:

반경 $r$에서의 총 파장 밀도

$$\rho_\text{wave}(r) \;=\; \lambda \;\sum_{i \in \{\text{얇은, 두꺼운, 기체, 팔, 돌출}\}} \rho_\text{wave}^{(i)}(r)$$

아래에 5개의 적분은 구성 요소당 하나씩 적혀 있습니다. 각 적분은 동일한 공간 지점에서 바이리온 질량 분포를 파장 질량 분포로 변환합니다.

5.1 벌지 – 3D 셸 통합

$$\rho_\text{wave}^{(b)}(r) \;=\; \int_0^{r_\text{max}} \rho_b(r’)\;\mathcal{K}_b\!\left(\sqrt{r^2 + r’^2}\right)\;4\pi r’^2\,dr’$$

적분은 반지름 $r’$의 동심 구형 껍질에 대해 이루어집니다. 중심에서 반지름 $r$ 의 필드 포인트는 단극 근사법에서 각 쉘을 유효 간격 $D = \sqrt{r^2 + r’^2}$ 에서 보게 됩니다. 적분은 $r_\text{최대} = 6\,r_b$까지 확장되며, 그 이상에서는 벌지 밀도가 수치적으로 무시할 수 있는 수준입니다.

5.2 씬 디스크 – 2D 링 통합

$$\rho_\text{wave}^{(\text{thin})}(r) \;=\; \int_0^{R_\text{max}} \시그마_\text{thin}(R’)\;\mathcal{K}_\text{thin}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$

디스크는 반지름 $R’$와 무한대 폭 $dR’$의 동심원 고리로 분해되며, 각 고리는 표면 질량 $\Sigma_\text{thin}(R’)\,2\pi R’\,dR’$을 갖습니다. 동일한 단극 근사가 적용됩니다: 중심으로부터 반경 $r$ 의 파장은 각 고리로부터 유효 이격 $D = \sqrt{r^2 + R’^2}$ 의 기여를 받습니다. 적분 범위는 $R_\text{최대} = 8\,R_d$입니다.

5.3 두꺼운 디스크 – 2D 링 통합

$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{R_\text{max}} \시그마_\text{두께}(R’)\;\mathcal{K}_\text{두께}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$

얇은 디스크 통합과 동일하며, 소스 밀도는 $\Sigma_\text{thick}(R’)$이고 커널 파라미터는 $\alpha_\text{thick} = 1/(c_\text{disk}\,\cdot 1.5\,R_d)$로 설정합니다. 두꺼운 디스크의 방사형 범위가 넓을수록 파동 일관성 범위가 약간 더 넓어집니다.

5.4 가스 링 – 중앙 디플로이먼트와 2D 링 통합

$$\rho_\text{wave}^{(\text{gas})}(r) \;=\; \int_{R_\text{hole}}^{R_\text{max}} \Sigma_\text{gas}(R’)\;\mathcal{K}_\text{gas}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$

기체 분포에는 중앙 구멍이 있으며, 적분 하한에서 $R_\text{hole} = 0.5\,R_g$에서 방사형 컷오프에 의해 포착됩니다. 이 컷오프 밖에서 기체는 항성 원반보다 더 멀리 확장되며, 이는 더 큰 특성 스케일 $R_g = 1.7\,R_d$에 반영되어 일관성 길이 $\ell_\text{gas} = c_\text{disk}\,R_g$에 입력됩니다.

5.5 스파이럴 암 초과 – 진폭이 감소된 2D 링 통합

$$\rho_\text{wave}^{(\text{arm})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{arm}(R’)\;\mathcal{K}_\text{arm}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$

나선형 암은 10\%$ 수준에서 얇은 디스크 표면 밀도의 축 방향 평균 향상으로 처리되며, 자체 일관성 길이 $\ell_\text{arm} = c_\text{arm}\,R_d$를 갖습니다. 따라서 커널은 나선형 구조의 방위각 집중도를 반영하여 얇은 디스크 커널보다 더 좁습니다.

6. 동봉된 파동 질량 및 예측 회전 곡선

총 파장 밀도 $\rho_\text{wave}(r)$가 알려지면 반경 $R$의 구 내에서 둘러싸인 파장 질량은 방사형 적분을 통해 구할 수 있습니다:

밀폐된 파장 질량

$$M_\text{wave}(R) \;=\; \int_0^{R} 4\pi\,r^2\,\rho_\text{wave}(r)\,dr$$

그런 다음 반경 $R$에서 예측된 원주 속도는 뉴턴 관계에서 구적법에서 바이론 및 파장 기여도를 결합하여 구합니다:

예측된 원주 속도

$$V_c^2(R) \;=\; V_\text{bar}^2(R) \;+\; \frac{G\,M_\text{wave}(R)}{R}$$

바이리오닉 속도 $V_\text{bar}(R)$ 는 그 자체로 4개의 디스크형 구성 요소(각 지수 프로파일에 대한 Freeman 1970 공식)와 벌지(Hernquist 동봉 질량 공식)의 기여도의 이차 합입니다:

$$V_\text{bar}^2(R) \;=\; V_\text{bulge}^2 + V_\text{thin}^2 + V_\text{thick}^2 + V_\text{gas}^2$$.

여기서 각 $V_i(R)$ 은 해당 질량 분포의 표준 뉴턴 원형 속도입니다.

7. 이론 수준 매개 변수

은하에 적용되는 전체 BeeTheory 프레임워크는 이론 수준에서 다섯 가지 매개 변수를 포함합니다. 이 매개변수는 은하마다 다르지 않은 보편적인 것입니다.

기호	의미	역할
$K_0$	파동 질량 진폭	웨이브 커널의 무차원 스케일을 설정합니다.
$c_\TEXT{SPH}$	3D 기하 상수	구형 소스의 비율 $\ell/r_\text{scale}$ (벌지)
$c_\text{디스크}$	2D 기하 상수	디스크 및 링 소스의 비율 $\ell/R_\text{scale}$
$c_\text{arm}$	나선형 기하 상수	방위각으로 집중된 팔 초과에 대한 $\ell/R_d$ 비율
$\lambda$	글로벌 파장 필드 커플링	총 파장 밀도 스케일 조정

매개변수의 범용성

다섯 가지 파라미터는 모두 전역적입니다. 은하, 왜성, 거대 나선에 동일한 수치값이 적용됩니다. 은하별 정보는 5개의 관측 입력 $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\-star)$를 통해서만 입력됩니다. 이 모델에는 은하별 조정 가능한 매개변수가 포함되어 있지 않습니다.

8. 계산의 단방향 특성

오픈 체인 – 피드백 없음

전체 계산은 입력에서 출력으로 한 방향으로 흐릅니다. 광도계와 21cm 관측을 통해 바이리온 분해를 결정합니다. 바이리오닉 분해는 파장 밀도를 결정합니다. 파장 밀도는 둘러싸인 파동 질량을 결정합니다. 동봉된 파동 질량은 예측된 회전 곡선을 결정합니다. 어떤 시점에서도 회전 곡선은 계산의 이전 단계에 영향을 미치지 않습니다.

이러한 단방향성에는 세 가지 중요한 결과가 있습니다.

(a) 다섯 가지 이론 수준의 매개변수가 고정되면 회전 곡선은 적합이 아닌 엄격한 예측입니다. 관찰된 회전 곡선과의 비교는 보정이 아니라 테스트입니다.

(b) 이 모델에는 은하별 조정 메커니즘이 없습니다. 회전 곡선 예측의 모든 수정은 입력 벡터 $(T,\,R_d,\,\Sigma_d,\,M_\text{HI},\,\Upsilon_\star)$의 수정 또는 우주 이론 수준의 매개 변수 $(K_0,\,c_\text{sph},\,c_\text{disk},\,c_\text{arm},\,\lambda)$의 변경에서 이루어져야 합니다.

(c) 기준 은하에서 $람다$를 보정하는 것은 해당 은하의 자전 곡선에 맞추는 것과는 다릅니다. 보정을 통해 하나의 전역 수치가 결정되고, 기준 은하의 다른 모든 반경에서의 회전 곡선과 다른 모든 은하의 회전 곡선은 보정된 프레임워크에 대한 엄격한 예측이 됩니다.

9. 중앙 표면 밀도의 역할(참고 XI 개정)

참고 XI의 진단을 통해 잔여 예측 오차는 디스크 스케일 길이 $\R_d$와 무관하게 중심 바이론 표면 밀도 $\Sigma_d$와 강한 상관관계가 있음을 확인했습니다. 위에 제시된 공식화는 이 발견이 통합되기 전의 모델 버전으로, 일관성 길이 표현식 $\ell_i = c_i\,R_d$에서 $\R_d$만 사용합니다.

세분화가 들어가는 곳

개선된 모델에서 일관성 길이 $\ell_i$ 는 $\R_d$ 와 $\Sigma_d$ 모두에 의존하게 되며, 엄격한 선형 관계 $\ell_i = c_i\,R_d$ 를 노트 XI에서 확인된 잔차를 흡수하는 함수 $\ell_i = c_i\,R_d\,\phi(\Sigma_d/\Sigma_\text{ref})$ 로 대체하게 됩니다. 후속 노트에서 $\phi$의 함수 형태와 그 파라미터를 먼저 22은하 보정 세트에서 결정한 다음 나머지 SPARC 샘플에 대한 블라인드 예측으로 검증할 것입니다.

계산의 단방향 구조는 이 세분화로 인해 유지됩니다. $\Sigma_d$는 관측 입력이고, 수정된 일관성 길이는 동일한 컨볼루션 적분으로 피드되며, 회전 곡선은 이전과 같이 나타납니다. 단 하나의 연산 링크, 즉 두 번째 관측 가능성에 대한 $\ell_i$의 종속성만 추가됩니다.

10. 방법론 요약

1. 입력값. 은하당 다섯 가지 관측값: 허블 유형 $T$, 원반 규모 $R_d$, 표면 밝기 $\Sigma_d$, HI 질량 $M_\text{HI}$, 우주 항성 질량 대 빛의 비율 $\Upsilon_\-star$.

2. 바리오닉 분해. 다섯 가지 구성 요소: 벌지($T \leq 4$인 경우), 얇은 디스크, 두꺼운 디스크, 가스 링, 나선형 암 초과. 각각은 분석 밀도 프로파일을 가지고 있습니다.

3. 웨이브 커널. 각 구성 요소의 기하학적 범위에 따라 결정되는 일관성 길이 $\ell_i = c_i\,R_\text{scale}$를 갖는 범용 유카와형 형식 $\mathcal{K}_i(D) = K_0\,(1 + \alpha_i D)\,e^{-\alpha_i D}/D^2$입니다.

4. 컨볼루션. 각 컴포넌트는 링(2D 컴포넌트) 또는 쉘(3D 벌지)에 대한 1차원 적분을 통해 파장 밀도를 생성합니다. 총 파장 밀도는 다섯 가지 구성 요소의 합으로, 글로벌 커플링 $\람다$에 의해 스케일링됩니다.

5. 출력. 동봉된 파동 질량 $M_\text{wave}(R)$을 적분하고 바이론 속도 $V_\text{bar}(R)$와 결합하여 예측된 회전 곡선 $V_c(R)$을 구합니다.

6. 이론 수준 매개변수. (K_0,\,c_\text{sph},\,c_\text{디스크},\,c_\text{암},\,\람다)$ – 범용, 은하별 튜닝 없음. 연구 중인 개선안은 $\Sigma_d$에 대한 종속성을 추가할 것입니다.

7. 방향. 입력 → 바이리온 → 파장 → 회전 곡선. 피드백 없음. 회전 곡선은 맞춤이 아닌 예측입니다.

참고 문헌. Lelli, F., McGaugh, S. S., Schombert, J. M. – SPARC: 스피처 광도계와 정확한 회전 곡선을 이용한 175개 원반 은하의 질량 모델, AJ 152, 157 (2016). – Freeman, K.C. – 나선 은하와 S0 은하의 원반에서, ApJ 160, 811 (1970). – Hernquist, L. – 구형 은하와 벌지에 대한 분석 모델, ApJ 356, 359 (1990). – Broeils, A. H., Rhee, M.-H. – 나선 은하와 불규칙 은하에 대한 21cm의 짧은 WSRT 관측, A&A 324, 877 (1997). – 맥거, S. S. – 은하 자전의 세 번째 법칙, Galaxies 2, 601 (2014). – 보비, J., 릭스, H.-W. – 은하수 원반 표면 밀도 프로파일의 직접 동역학 측정, ApJ 779, 115 (2013). – Arnett, D. – 초신성과 핵 합성, 프린스턴 (1996). – 두테르트르, X. – 꿀벌 이론™: 중력의 파동 기반 모델링, v2, BeeTheory.com (2023).

공식화:은하계 규모의 벌 이론 계산