BeeTheory – 통계 분석 – 2025

핏 읽기

χ²/dof = 0.475

0에서 ∞ 사이의 숫자로, 모델이 데이터에 잘 맞는지, 너무 잘 맞는지, 아니면 전혀 맞지 않는지를 알려줍니다. BeeTheory 은하계 시뮬레이션에서 이 숫자가 의미하는 바는 무엇이며, 진정으로 확실한 결과를 얻기 위해 필요한 것은 무엇일까요?

χ²/dof란 무엇인가요?

모델이 값 _Vmodel(_Ri)를 예측하고 데이터가 측정 불확실성 _σi와 함께 관측 값 _Vobs(_Ri)를 제공하는 경우, 감소된 카이제곱은 다음과 같습니다:

\(\frac{\chi^2}{\mathrm{dof}}=\frac{1}{N-p}\sum_{i=1}^{N}\left(\frac{V_{\mathrm{model}}(R_i)-V_{\mathrm{obs}}(R_i)}{\sigma_i}\right)^2\)

여기서 N은 데이터 포인트의 수이고 p는 자유 매개변수의 수입니다.

여기: N = 16 Gaia 2024 포인트, p = 2 적합 파라미터, K 및 α, 따라서 dof = 14입니다.

합계의 각 항은 측정 불확실성 단위로 표시되는 잔차, 즉 풀입니다.

0.5의 풀은 모델이 0.5 시그마 정도 벗어났음을 의미합니다. 풀이 2.0은 2 시그마 불일치를 의미하므로 조사할 가치가 있습니다.

[라텍스]\frac{\chi^2}{\mathrm{dof}}\approx1\quad\Rightarrow\quad\text{모델이 불확실성 수준에서 데이터에 적합함}[/라텍스] [라텍스]\frac{\chi^2}{\mathrm{dof}}\ll1\quad\Rightarrow\quad\text{모델이 너무 잘 맞음: 불확실성이 과대 추정되었거나 과적합일 수 있음}[/라텍스] [라텍스]\frac{\chi^2}{\mathrm{dof}}\gg1\quad\Rightarrow\quad\text{모델이 잘못되었거나 불확실성이 과소평가됨}[/라텍스]

우리 번호: 0.475

0.475

χ² / dof – N = 16 – p = 2 – dof = 14

통계적으로 우수

모델이 틀린 것이 아닙니다.

불확실성이 약간 클 수 있습니다.

χ²/규모

0.475 ← 우리

1.0 이상적

00.511.522.53+

< 0.3
오버핏

0.3-0.7
우수

0.7-1.3
양호

1.3-2.0
수용 가능

> 2.0
poor

0.475의 값은 우수 영역에 속합니다. 16개 데이터 포인트 전체에서 평균 잔차는 0.69σ에 불과합니다.

물리적으로 이것은 BeeTheory가 거의 모든 관측 반경에서 원주 속도를 측정 불확실도 1 미만 이내로 예측한다는 의미입니다.

하지만 “우수”는 “입증”과는 다릅니다.

χ²/dof = 0.475는 측정 불확실성 _σi가 약간 과대평가되었음을 의미할 수도 있습니다. 실제 오차가 30% 더 작다면 동일한 모델에서 χ²/dof ≈ 1.0이 나올 것입니다.

χ²만으로는 “모델이 매우 우수하다”와 “오차가 약간 크다”를 구분할 수 없습니다. 이것은 표준적인 통계적 모호성입니다.

잔여 사항 – 포인트별

원시 숫자 0.475는 정보를 숨깁니다. 개별 당김을 보면 모델이 어느 지점에 못을 박고 어디에서 어려움을 겪는지 등 핏의 구조를 알 수 있습니다.

풀 플롯: 각 Gaia 2024 데이터 포인트에 대한 [라텍스](V_{\mathrm{model}}-V_{\mathrm{obs}})/\sigma_i[/라텍스][/latex]

R(kpc)	Vobs	σ	V모델	잔여	Pull
잔여물 로드 중…

1σ 이내 15/16점

R = 27.3 kpc를 제외한 모든 데이터 포인트의 |pull|이 1.0 미만입니다. 가우스 오차가 있는 잘 보정된 모델에서는 약 68%의 포인트가 1σ 이내에 있을 것으로 예상되지만, 여기서는 94%입니다.

이는 모델이 우수하거나 내부 포인트의 오차 막대가 약간 너무 크다는 것을 의미합니다.

R = 27.3kpc – pull = +1.53σ의 이상치

이 모델은 _Vc = 195.3km/s로 예측하는 반면, Gaia는 173 ± 17km/s로 측정합니다.

불일치는 22.3km/s, 즉 1.53σ입니다. 이는 통계적으로 놀라운 수준은 아니지만 물리적으로 큰 반경에서 모델이 너무 느리게 감소한다는 점에서 유의미합니다.

최적 맞춤 파라미터 K = 0.01349

결합 상수 K = 0.01349 kpc-¹는 가시 질량 단위당 생성되는 BeeTheory 암흑 질량장의 진폭입니다.

이는 16개의 Gaia 데이터 포인트에서 χ²를 최소화하고 α = 0.0744 kpc-¹를 뉴비 + 가이아 결합 적합도에서 고정하여 결정했습니다.

\(K=0.01349\,\mathrm{kpc}^{-1}\) \(\lambda=K\ell^2=\frac{K}{\alpha^2}=\frac{0.01349}{(0.0744)^2}=2.44\) \(\rho_{\mathrm{dark}}(R_\odot=8\,\mathrm{kpc})=0.37\,\mathrm{GeV/cm^3}\)

관측된 로컬 암흑 물질 밀도는 다음과 같습니다:

\(\rho_{\mathrm{obs}}(R_\odot)=0.39\pm0.03\,\mathrm{GeV/cm^3}\)

무차원 결합 λ = 2.44는 단일 성분, 다중 성분 및 원자 H₂ 비교를 포함한 BeeTheory 적합에서 관찰된 λ ∈ [2, 7] 범위에 속합니다.

펨토미터 단위부터 킬로파섹 단위까지 동일한 규모의 이러한 보편성은 BeeTheory 프레임워크의 가장 강력한 내부 일관성 검사입니다.

K가 단일 디스크보다 더 작은 이유

얇은 원반만을 소스로 사용한 이전 적합도에서는 _Ksingle = 0.038 kpc-¹이 나왔습니다. 6개의 은하 성분이 모두 암흑장에 기여하기 때문에 전체 바이리온 소스는 2.8배 더 큽니다.

따라서 동일한 암흑 질량을 생성하려면 K는 비례적으로 더 작아야 합니다.

\(K_{\mathrm{multi}}\approx\frac{K_{\mathrm{single}}}{2.8}=\frac{0.038}{2.8}=0.0136\,\mathrm{kpc}^{-1}\)

이것이 바로 핏이 제공하는 것입니다. 이는 우연이 아니라 일관성을 확인하는 것입니다.

27.3kpc가 모든 모델에 어려운 이유

가이아 2024의 가장 바깥쪽 데이터 포인트인 _Vc(27.3kpc) = 173 ± 17km/s는 비이론에게만 어려운 것이 아닙니다. 은하수 회전 곡선에 적용되는 모든 암흑 물질 모델에서 가장 어려운 지점입니다.

모델	Vc(27.3kpc) 예측값	잔여	Pull	χ²/dof
BeeTheory	195.3km/s	+22.3	+1.53σ	0.475
NFW 프로필	~198km/s	+25	+1.5σ	0.44
아이나스토 α = 0.91	~191km/s	+18	+1.1σ	0.38
등온 후광	~218km/s	+45	+2.6σ	2.6

표준 두 매개변수 모델은 1σ 이내에서 _Vc(27.3kpc) = 173km/s를 재현하지 못합니다.

세 가지 그럴듯한 설명이 있습니다.

측정 자체: R = 27.3kpc는 가장 먼 지점이며, 운동 추적자가 가장 적고 시스템적 불확실성이 가장 큽니다. Gaia DR4로 인해 값이 바뀔 수 있습니다.
가스 디스크 기여도: HI 디스크는 약 25kpc까지 확장되며 큰 반경에서 바이리오닉 질량에 기여합니다. 이를 별도의 구성 요소로 포함하면 감소가 더 심해질 수 있습니다.
일관성 길이가 작을수록 (α = 0.12kpc-¹ 또는 ℓ = 8.3kpc) 가장 바깥쪽 지점은 더 잘 맞지만 안쪽 고원은 더 나빠집니다.

진정으로 견고한 결과는 어떤 모습일까요?

현재 적합도는 통계적으로 양호합니다. 그러나 2개의 자유 매개 변수와 16개의 데이터 포인트가 있는 모델에서 잘 맞는다는 것은 입증된 것과는 다릅니다.

0.475

현재 χ²/dof, 16포인트, 파라미터 2개

< 0.8

대상: 확장된 데이터 집합에 적합

~1.0

이상적: 은하계 샘플 전반에 걸쳐 보정됨

솔리드 검증을 위한 요구 사항

트레이닝 데이터 세트의 χ²/dof가 양호합니다: 달성됨. Gaia 2024에서 χ²/dof = 0.475.
로컬 암흑 물질 밀도 수정: 달성. 0.37 GeV/cm³ 예측 대 0.39 ± 0.03 관측.
일관된 무차원 커플링: 달성. λ = 2.44, H₂ 캘리브레이션의 λ ≈ 3.5와 비교.
! 가장 바깥쪽 가이아 지점을 1σ 이내로 재현하세요: 완전히 달성하지 못했습니다. R = 27.3kpc에서의 잔차는 +1.53σ입니다.
참고 다른 은하계에 대한 블라인드 예측: 아직 완료되지 않았습니다. SPARC 카탈로그에서 자연스러운 테스트를 제공합니다.
참고 ℓ는 첫 번째 원칙에서 도출한 것입니다: 아직 완료되지 않았습니다. ℓ는 현재 적용되어 있으며 파생되지 않았습니다.
참고 클러스터 규모 유효성 검사: 아직 완료되지 않았습니다. 불릿 클러스터 렌더링은 핵심 테스트입니다.
→ 30kpc 이상의 확장된 회전 곡선: Gaia DR4는 단기적으로 중요한 테스트를 제공해야 합니다.

정직한 과학적 현황

현재 형태의 비이론은 은하수 자전 곡선과 최고의 경험적 암흑 물질 프로파일에 부합하는 물리적 동기가 부여된 현상학적 프레임워크로, 물리적 메커니즘을 가지고 있다는 장점이 있습니다.

K와 ℓ는 도출된 것이 아니라 장착된 것이기 때문에 엄격한 의미에서 아직 이론이 아닙니다.

파동 질량 가정에서 ℓ = f(_Lsource)를 도출하고, SPARC에서 보편적으로 테스트한 다음, Gaia DR4에 대한 블라인드 예측을 수행하면 완전한 이론으로 가는 길은 명확합니다.

Data: Ou, X. et al., MNRAS 528, 2024. Model: BeeTheory v2, Dutertre 2023, 다중 구성 요소 적합. SPARC 참조: Lelli, McGaugh, Schombert, AJ 152, 2016. 불릿 클러스터: Clowe et al., ApJL 648, 2006.