On the triple nature of the PSR J0435+3233 system

기술 요약: PSR J0435+3233 시스템의 삼중 성계 특성에 관하여

문제 제기
본 논문은 Commensal Radio Astronomy FAST Survey(Wu et al. 2026, 이하 W26)에 의해 발견된 밀리초 펄사(MSP)인 PSR J0435+3233의 타이밍 에페머리스(ephemeris)에서 나타나는 중대한 이상 현상을 다룬다. W26은 관측된 스핀 다운율($\dot{P}$)이 $4.88 \times 10^{-17}$ s s$^{-1}$라고 보고했는데, 이는 일반적인 은하계 MSP보다 두 자릿수 높은 수치이다. 만약 이것이 고유한 특성이라면, 이는 약 $10^6$년의 특성 연령, $1.26 \times 10^{10}$ G의 표면 자기장, 그리고 약 $5.9 \times 10^{37}$ erg s$^{-1}$의 스핀 다운 출력($\dot{E}$)을 의미하며, 이는 현재의 MSP 형성 및 진화 모델에 도전하는 결과이다. 또한, Fermi LAT(Large Area Telescope)가 감마선 펄세이션을 검출했음에도 불구하고, 그 플럭스(flux)는 W26이 예측한 높은 $\dot{E}$와 일치하지 않았으며, 이는 비정상적으로 낮은 감마선 효율성을 시사한다. 아울러, W26에서 관측된 스핀 주파수 및 궤도 주파수 미분값은 상당한 가변성(예: 5년간 $|\dot{P}|$가 약 60% 증가)과 상관적 진화를 보였는데, 이는 단순한 이중성이나 표준적인 "스파이더(spider)" 펄사 모델로는 설명하기 어렵다.

연구 방법론
저자들은 이러한 비정상적인 타이밍 동작이 고유한 특성이 아니라, 먼 거리에 있는 삼차 동반성(tertiary companion)의 중력장 내에서 내부 이중성계(내부 펄사와 저질량 동반성)가 가속됨으로써 발생하는 현상이라고 제안한다.

1. 주파수 미분 분석: 저자들은 W26에서 보고된 스핀 주파수($f$)와 궤도 주 frequency($f_B$) 미분의 상관적 진화를 분석하였다. 저자들은 $\dot{f}/f \simeq \dot{f}_B/f_B \simeq -a_l/c$ 관계를 따르는 $\dot{f}$와 $\dot{f}_B$의 거의 동일한 변화를 입증함으로써, 이것이 선시선 방향 가속도($a_l$)에 의한 것임을 보여주었다. 이를 통해 원형 궤도(고차 미분에서 부호 반전이 필요함)를 배제하고 편심 궤도를 가진 외곽 궤도를 지지하였다.
2. 케플러 파라미터 추정: 저자들은 W26의 다섯 가지 스핀 주파수 미분값과 세 가지 궤도 주파수 미분값을 사용하여, MCMC(Markov-Chain Monte-Carlo) 샘플러를 통해 외곽 궤도의 다섯 가지 케플러 파라미터(투영 장반경 $x_O$, 주기 $P_{B,O}$, 이심률 $e_O$, 근일점 인수 $\omega_O$, 진근점 이각 $\nu_O$)를 추정하였다.
3. 계층적 삼중성 타이밍 모델: 저자들은 PINT(그래디언트 기반 최적화를 위한 JAX 사용)를 사용하여 계층적 삼중계 시스템 타이밍 모델을 구현하였다. 이 모델은 내부 궤도 위상을 계산하기 전에 외곽 궤도의 로머 지연(Rømer delay)을 고려한다. 저자들은 W26의 바리센터(barycentric) 도착 시간(ToA)에 이 모델을 적합(fit)시켰다.
4. 감마선 통합: 타이밍 기선(baseline)을 정교화하고 확장하기 위해, 저자들은 Fermi-LAT 감마선 광자를 라디오 타이밍 솔루션에 따라 폴딩(folding)하였다. 저자들은 결합 로그 가능도(joint log-likelihood)를 최대화하기 위해 감마선 H-통계량에 기초하여 라디오 타이밍으로부터 얻은 사후 샘플(posterior samples)의 가중치를 재설정함으로써, 타이밍 기선을 2008년부터 2024년까지 확장하였다.
5. 광학적 대조물 식별: 저자들은 펄사 위치를 Gaia DR3 및 2MASS 카탈로그와 교차 대조하여 삼차 동반성의 잠재적 광학적 대조물을 식별하였으며, 해당 천체의 고유 운동, 시차 및 광도 특성을 분석하였다.

주요 결과

• 삼중계 확인: 데이터는 계층적 삼중계 시스템에 의해 가장 잘 설명된다. 내부 이중성은 MSP와 저질량 동반성(아마도 백색 왜성)으로 구성되며 약 8일의 궤도를 가진다. 이 이중성은 공통 질량 중심을 기준으로 긴 주기의 편심 궤도($P_{B,O} \approx 72$ 년, $e_O \approx 0.60$)를 가진 삼차 동반성과 함께 공전한다.
• $\dot{P}$ 이상 현상 해결: 관측된 큰 $\dot{P}$는 주로 외곽 궤도로 인한 변화하는 선시선 방향 가속도의 결과물이다. 고유 스핀 다운율은 관측된 값보다 최소 두 자릿수 낮게 제한되며, 이는 일반적인 은하계 MSP와 일치한다.
• 확장된 감마선 검출: 삼중계 타이밍 모델은 2008년 8월 Fermi-LAT 데이터의 시작 시점까지 펄사 위상을 성공적으로 예측하며, 이는 이전 모델들이 라디오 데이터 범위를 벗어나면 결맞음(coherence)을 잃었던 점을 크게 개선한 것이다. 결합된 라디오/감마선 솔루션은 높은 H-통계량($H=171.9$)을 산출한다.
• 삼차 동반성 식별: 펄사로부터 11 mas 떨어진 곳에 위치한 별 2MASS J04353375+3233080이 삼차 동반성의 광학적 대조물로 식별되었다. Gaia의 시차와 광도 측정을 바탕으로, 이 별은 거리 약 1.5~2.2 kpc에 위치한 주계열성(분광형 G0V–F0V, 질량 약 $1.2 M_\odot$)이다.
• 궤도 기하학: 외곽 궤도의 투영 장반경은 약 2477초이며 궤도 경사각 $i_O \approx 31^\circ$ 또는 $149^\circ$이다. 모델은 외곽 궤도의 근일점 통과가 2036년 초에 일어날 것임을 예측한다.
• 타이밍 잔차: 새로운 모델은 FAST ToA의 잔차 RMS를 W26의 $1.79 \mu$s에서 $1.49 \mu$s로 감소시켰으며, 축소된 $\chi^2$ 값은 1.002이다. 그러나 저자들은 조석 상호작용이나 상대론적 시간 지연과 같은 고차 효과가 여전히 잔차에 존재할 수 있다고 언급하였다.

의의 및 주장
본 논문은 PSR J0435+3233이 계층적 삼중계 시스템임을 결론지으며, 이를 통해 기존의 모순되었던 스핀 다운율 및 감마선 효율 문제를 해결한다. 저자들은 다음과 같이 주장한다:

1. 이 시스템은 이례적인 형성 경로를 가정하거나 극단적인 자기장을 끌어들이지 않고도, 이전에 나타났던 비정상적인 타이밍 미분값들을 자연스럽게 설명해 준다.
2. 삼차 동반성의 식별과 외곽 궤도 파라미터의 결정은 시스템에 대한 일관된 물리적 모델을 가능하게 한다.
3. 이 시스템은 자기 중력 천체에 대한 강한 등가 원리(SEP) 및 자유 낙하의 보편성(UFF)을 테스트하는 데 중요한 잠재력을 가진다. 잘 알려진 삼중계인 PSR J0337+1715와 비교했을 때 더 큰 투영 장반경을 가진 PSR J0435+3233은, 특히 노드트 베크트 효과(Nordtvedt effect)와 관련하여 UFF 위반에 대해 더욱 엄격한 제약을 제공할 수 있다.
4. 삼차 별의 시선 속도 변화를 측정하기 위한 고해解 분광 관측과 2036년 근일점 접근 시까지 지속되는 펄사 타이밍 관측은, 시스템의 질량비, 경사각 및 거리를 정밀하게 결정하여 이러한 중력 테스트를 더욱 정교하게 만들 것이다.