넓은 이중성 3D 궤도 모델링을 통한 중력 상수 베이지안 추정: 일반 알고리즘 및 HARPS 파일럿 연구

초록

이 논문은 Gaia와 HARPS에서 얻은 3차원 상대 위치·속도 데이터를 이용해 넓은 이중성의 궤도를 베이지안 방식으로 전면 분석하는 일반 알고리즘을 제시한다. 32개의 샘플에 적용한 결과, 가속도가 (g_N>10^{-9},\mathrm{m,s^{-2}})인 경우는 뉴턴 중력과 일치하지만, (g_N<10^{-9},\mathrm{m,s^{-2}}) 영역에서는 (\Gamma\equiv\log_{10}\sqrt{G/G_N})가 양의 값을 보여 약 3.5σ 수준의 편차가 관측된다. 다만 이 신호는 한 쌍(HD 189739/HD 189760)에 크게 의존하므로, 더 많은 고정밀 시스템 확보와 삼중성 검증이 필요하다.

상세 분석

본 연구는 넓은 이중성(WSB)의 3차원 상대 변위 (\mathbf r)와 상대 속도 (\mathbf v)를 한 시점에서 정확히 측정하면, 중력 파라미터 (G)가 주어졌을 때 타원 궤도 해를 구할 수 있다는 점에 착안한다. 그러나 궤도 기하학적 파라미터(반장축 (a), 이심률 (e), 근일점 인자 (\phi_0), 위상 (\phi), 경사 (i), 회전 (\theta))와 (G) 사이에 강한 퇴화가 존재하므로, 각 시스템에 대해 베이지안 사후 확률분포를 구하고 이를 통계적으로 합치는 것이 필수적이다.

알고리즘은 다음과 같은 핵심 절차를 포함한다.

1. 관측 입력: Gaia DR3에서 제공하는 정확한 적경·적위 차이 ((x’,y’))와 HARPS에서 얻은 고정밀 방사속도 차이 ((v_{z}’)) 및 적경·적위 방향의 상대 속도 ((v_{x}’,v_{y}’))를 사용한다. 라인오브사이트 거리 차 (\Delta d)는 불확실성이 크지만, 본 파일럿에서는 (\sigma_{\Delta d}<s/\sqrt2)인 경우만 선택하였다.
2. 모델 파라미터: 관측된 ((x’,y’))가 거의 무오차이므로 자유 파라미터는 ({e,,\phi_0,,\phi,,i})와 질량 스케일 팩터 (\log_{10}f_M), 중력 변형 파라미터 (\Gamma)로 제한된다. 반면 (a)와 (\theta)는 식 (2)–(7)에서 ((x’,y’))와 ((e,\phi_0,\phi,i))에 의해 고정된다.
3. 베이지안 프레임워크: 사후 확률 (p(\Theta|{\rm data})\propto L({\rm data}|\Theta),\prod_l f_{\rm pr}(\Theta_l)) 로 정의한다. 여기서 (\Theta={e,i,\phi_0,\Delta\phi,\log_{10}f_M,\Gamma})이며, (\Delta\phi\equiv\phi-\phi_0)이다.
   • 우도는 관측값과 모델값의 차를 정규분포 가정으로 구성한 (\chi^2) 형태(식 12)이다.
   • 사전은 물리적 논리와 경험적 지식을 반영한다. 경사 (i)는 (|\sin i|) 분포, 근일점 인자 (\phi_0)는 균등, 위상 (\Delta\phi)는 케플러 법칙에 따라 (\propto (1-e^2)^{3/2}/(1+e\cos\Delta\phi)^2) 로 설정한다. 질량 스케일 (f_M)는 평균 0, σ = 5 %인 로그정규분포, 이심률 (e)는 열역학적 분포 (f(e)=2e) 혹은 무분포 두 경우를 시험한다. (\Gamma)는 (