은하 중심 태양의 움직임: CDM, MOND, MOG 시나리오 비교

초록

최근 측정된 태양의 3차원 운동 정보를 이용해, -4.5 Gyr에서 현재까지 태양의 궤적을 뉴턴 역학 기반 CDM(NFW·로그형), MOND, MOG 네 모델로 수치 적분하였다. LSR 회전속도 Θ₀=268 km s⁻¹와 IAU 권고값 Θ₀=220 km s⁻¹ 두 경우를 비교했으며, Θ₀=268 km s⁻¹에서는 태양이 탄생할 때 12–27 kpc, Θ₀=220 km s⁻¹에서는 8.8–9.3 kpc 정도의 거리에서 시작했음을 보였다. CDM 모델이 가장 넓은 궤적을 보이며, MOND·MOG는 상대적으로 제한된 범위에 머문다.

상세 분석

이 논문은 태양의 과거 궤적을 구하기 위해 은하 중심 좌표계에서 3차원 운동 방정식을 -4.5 Gyr에서 현재까지 직접 적분한다는 점에서 독창적이다. 초기 조건은 Reid et al. (2009)의 최신 관측값을 사용했으며, 특히 LSR 회전속도 Θ₀를 268 km s⁻¹와 220 km s⁻¹ 두 값으로 나누어 분석했다. CDM 모델은 두 가지 형태—Navarro‑Frenk‑White(NFW) 프로파일과 로그형 포텐셜—를 채택했으며, 각각 디스크·버스·다크 매터 성분을 포함한다. NFW 모델은 ρₛ, r_vir, c 등 우주론적 파라미터를 Xue et al. (2008)의 평균값으로 설정했고, 로그형 모델은 Law et al. (2005)의 파라미터를 사용했다. MOND은 Milgrom의 가속도 임계값 A₀≈1.2×10⁻¹⁰ m s⁻²와 보간 함수 μ(X)=X/(1+X)를 적용했으며, 외부장 효과(EFE)를 무시할 정도로 A_e/A₀≈0.01–0.02 수준이라고 가정했다. MOG는 Moffat & Toth(2009)의 Yukawa‑형 수정 중 α와 μ(=λ⁻¹)를 은하 질량 M=6.5×10¹⁰ M☉에 맞춰 계산했다.

수치 적분은 Mathematica의 NDSolve를 기본 옵션으로 수행했으며, 결과는 x‑y, x‑z, y‑z 평면에 투영된 궤적으로 제시된다. Θ₀=268 km s⁻¹일 때, CDM 로그형과 NFW 모델은 태양이 약 25–27 kpc까지 외곽으로 이동한 뒤 다시 내측으로 돌아오는 비원형 궤적을 보인다. MOND은 약 12 kpc, MOG는 약 15 kpc 정도에서 시작해 비교적 원형에 가까운 궤적을 유지한다. 반면 Θ₀=220 km s⁻¹ 경우, 모든 모델이 8.8–9.3 kpc 범위 내에서 시작해 크게 차이나지 않는 궤적을 그린다. 이는 LSR 회전속도가 궤도 반경에 미치는 영향을 명확히 보여준다.

논문은 또한 은하 서식 가능 구역(GHZ) 관점에서 결과를 해석한다. 태양이 12 kpc 이상 외곽에서 탄생했다면 금속 함량이 낮아 행성 형성 및 생명 발달에 불리할 수 있다는 점을 지적한다. 따라서 Θ₀=268 km s⁻¹를 채택한 경우는 GHZ와의 일치도가 낮아, 실제 은하 회전속도에 대한 추가 관측이 필요함을 시사한다.

전반적으로, CDM 모델이 가장 큰 궤도 확장을 보이며, MOND·MOG는 상대적으로 제한된 궤적을 제공한다는 점은 각 이론이 은하 규모 중력에 어떻게 반응하는지를 직관적으로 보여준다. 또한, 초기 속도와 회전속도 선택이 결과에 큰 영향을 미치므로, 관측적 불확실성을 고려한 다중 시나리오 분석이 필수적임을 강조한다.