Light deflection in the gravitational field of a solar system body with finite distance of source and observer for sub-micro-arcsecond astrometry

기술 요약: 서브 마이크로-아르크초(sub-$\mu$as) 천문측정학을 위한 유한한 거리의 광원 및 관측자를 가진 태양계 천체의 중력장 내 빛의 굴절

문제 정의
고정밀 천문측정학의 주요 목표는 관측된 빛의 신호 방향을 기반으로 천체의 공간적 위치를 결정하는 것이다. Gaia와 같은 현재의 미션들은 마이크로-아크초($\mu$as) 정확도를 달le 성취하였으나, 제안된 미래 미션들(예: GaiaNIR, Theia, TOLIMAN, SHERA)은 약 0.1 $\mu$as를 목표로 하는 서브 마이크로-아크초(sub-$\mu$as) 정밀도를 지향한다. 측정 목표보다 약 10배 더 나은 이론적 모델 정확도(0.01 $\mu$as 또는 10 나노-아크초 임계값)를 달성하기 위해서는 태양계 내 빛의 전파에 대한 상대론적 모델링이 정교해져야 한다.

기존의 빛 궤적에 대한 해석적 해법들은 종종 광원과 관측자의 거리를 무한하다고 가정하거나, 서브-$\mu$as 수준에서는 무시할 수 있는 번거로운 항들을 포함하고 있다는 점에서 중요한 과제가 발생한다. 또한, 예상되는 수십억 개의 관측 데이터를 처리하기 위해 측지선 방정식(geodesic equations)을 수치적으로 적분하는 것은 계산적으로 불가능하다. 본 연구에서 다루는 구체적인 문제는 다음을 고려하여 관측자 위치에서의 단위 접선 벡터(unit tangent vector)에 대한 단순화된 해석적 표현식을 도출하는 것이다:

1. 중력체로부터 광원과 관측자의 유한한 거리.
2. 임의의 천체 형상, 내부 구조 및 회전 운동을 허용하는 질량 다중극(mass-multipoles)과 스핀 다중극(spin-multipoles)의 전체 집합.
3. 빛의 굴절각에 10 나노-아크초 이상의 기여를 하는 항들로만 제한된 축소.

방법론
본 연구는 1.5PN 근사(order $O(c^{-4})$)까지의 포스트-뉴턴(Post-Newtonian, PN) 체계를 활용한다. 고립된 천체의 중력장은 질량 다중극($M_L$)과 스핀 다중극($S_L$)으로 전개된 메트릭 텐서(metric tensor)로 기술된다.

1. 측지선 방정식 및 적분: 빛의 전파에 대한 측지선 방정식을 이전 연구에서 도입된 새로운 변수($\tau$ 및 $\xi$)를 사용하는 고급 적분법을 사용하여 해결한다. 이는 광원($x_0$)과 관측자($x_1$)가 유한한 거리에 있는 경계값 문제(boundary value problem)로 문제를 변환한다.
2. 단위 접선 벡터의 유도: 측지선 방정식을 적분하고 초기 방향($\sigma$)을 관측된 방향($k$)으로 변환함으로써 관측자에서의 단위 접선 벡터($n$)를 유도한다. 이 해법은 배경 방향($k$)과 질량 다중극($\Delta n^{ML}$) 및 스핀 다중극($\Delta n^{SL}$)에 의한 섭동의 합으로 표현된다.
3. 크기 분석 및 단순화: 개별 섭동 항들의 크기를 체계적으로 분석한다. 항들은 충격 매개변수($d_k$), 관측자까지의 거리($r_1$), 그리고 천체의 반지름($P$)에 따른 스케일링을 기준으로 분류된다.
   • 10 나노-아크초 미만의 기여를 하는 항들을 식별하여 제외한다.
   • 두 가지 특정 미분 연산($U^L$ 및 $V^L$)을 분석하며, 특정 성분($U^2_L$ 및 $V^2_L$)이 서브-$\mu$as 천문측정학에서 무시할 수 있음을 보여준다.
4. 축대칭 사례: 천체를 축대칭(예: 태양 및 거대 행성)으로 근사하는 경우, 질량 및 스핀 다중극은 영역 조화 계수($J_l$)를 사용하여 표현된다. 결과적으로 도출된 단위 접선 벡터의 표현식은 제1종 및 제2종 체비쇼프 다항식(Chebyshev polynomials)을 사용하여 나타낼 수 있음을 보여준다.

주요 기여**

• 단순화된 단위 접선 벡터: 본 논문은 정지 상태의 임의의 천체에 대해 "단순화된 단위 접선 벡터"(식 133)를 유도한다. 이 식은 10 나노-아크초 미만의 기여를 하는 항들을 버리면서도 서브-$\mu$as 정확도에 유효한 모든 항을 유지한다. 이는 임의의 유한한 광원 및 관측자 거리에 대해 유효하다.
• 체비쇼프 다항식 표현: 축대칭 천체의 경우, 빛의 굴절 항들을 체비쇼프 다다항식($T_l$ 및 $U_{l-1}$)을 사용하여 간결하게 표현할 수 있음을 입증한다. 이를 통해 임의의 차수의 질량 및 스핀 다중극에 대한 빛의 굴절 효과를 효율적으로 결정할 수 있다.
• 상한선 증명: 광원과 관측자의 거리가 무한할 때 계산된 총 빛 굴절이, 광원과 관측자가 유한한 거리에 있을 때의 빛 굴절 효과에 대한 엄격한 상한선(upper limit)임을 증명한다. 단, 관측자는 지구 근처(예: 라그랑주 점 L2)에 위치한다고 가정한다.
• 다중극 차수 요구사항: 수치 분석을 통해 서브-$\mu$as 정확도를 위해서는 태양계 천체 모델에 대해 질량 다중극은 최대 $l=8$차까지, 스핀 다중극은 최대 $l=3$차까지 포함되어야 함을 확립한다.

결과

• 무시 가능한 항들: 본 연구는 다중극의 고차 미분 및 특정 스핀 다중극 항(예: 스핀의 경우 $l=4, 6, \dots$)이 1 나노-아크초 미만의 기여를 하며 안전하게 무시될 수 있음을 정량화한다.
• 목성 및 토성 접선 광선: 목성과 토성의 접선 광선(grazing rays)에 대해, 무시된 모든 항의 합을 계산하였다. 토성의 경우 그 합은 $\le$ 14.9 nas이며, 목성의 경우 $\le$ 36.2 nas이다. 논문은 목성의 경우 충격 매개변수가 반지름보다 약간 클 경우($d_k \ge 1.7 P$), 무시된 항들이 10 nas 임계값 아래로 떨어진다고 언급한다. 극단적인 목성 접선 사례의 경우, 10 nas 정확도 목표를 유지하기 위해 정확한 질량 사중극(1PN) 항을 별도로 유지해야 할 수도 있다.
• 2PN 및 3PN 항: 본 논문은 고차 포스트-뉴턴 항들을 간략히 평가한다. 2PN 질량 단극(mass-monopole) 및 질량 사중극 항은 수 마이크로-아크초(예: 태양의 경우 3136 $\mu$as)에 달할 수 있으며 반드시 포함되어야 함을 발견했다. 3PN 항은 일반적으로 더 작지만 태양의 경우(11.2 $\mu$as) 무시할 수 없다.
• 수치 표: 본 논문은 태양과 거대 행성에 대한 질량 다중극, 스핀 다중극 및 고차 PN 항에 의한 빛 굴절의 상한치를 나열한 종합적인 표(IV, V, VI)를 제공한다.

의의
본 논문은 유도된 단순화된 단위 접선 벡터가 GaiaNIR 및 Theia와 같은 미래 우주 천문측정 미션의 효과적인 관측 데이터 환원을 위해 필수적이라고 주장한다. 10 나노-아크초 미만의 기여를 무시하면서도 충분히 정확하고 간결한 해석적 해법을 제공함으로써, 본 연구는 수십억 개의 천문측정 데이터를 처리하는 계산적 병목 현상 문제를 해결한다. 체비쇼프 다항식이 축대칭 천체에 대한 굴절을 기술할 수 있다는 데 대한 입증은 태양계 천체의 중력장 내 빛의 전파를 모델링하는 데 있어 계산적으로 효율적인 방법을 제공하며, 마이크로-아크초에서 서브 마이크로-아크초 천문측정으로의 전환을 용이하게 한다.