DD CrB의 새로운 물리·궤도 파라미터와 13년 주기의 행성 질량 제3천체 탐색

초록

본 연구는 B형 아광성 주계열성(sdB)과 M형 주계열성(dM)으로 이루어진 근접 이중성 DD CrB의 다중색 광도곡선, 방사속도곡선, 그리고 TESS 기반 정밀 식이시각을 이용해 시스템 파라미터를 재정밀화하고, 식이시각 변동(ETV)을 제3천체의 광행시효(LiTE) 효과로 해석한다. 로머 지연을 통한 질량비 q = 0.299 ± 0.009를 확보하고, PHOEBE 2.17 모델링에서 2차성의 복사율을 물리적 상한인 A₂ = 1.0으로 고정함으로써 T₂ ≈ 2360 K, R₂ ≈ 0.16 R☉를 도출하였다. O–C 분석 결과 13년 주기의 이심률 e ≈ 0.46인 목성 질량(≈1 M_J) 제3천체가 가장 적합한 설명임을 확인하였다.

상세 분석

DD CrB는 전형적인 HW Vir형 포스트공통외피 이중성(PCEB)으로, sdB 주성분의 고온(≈29 000 K)과 dM 보조성분의 저온(≈2 300 K) 사이에 극심한 온도 차이가 존재한다. 이러한 온도 차이는 강력한 반사효과를 유발해 광도곡선의 2차 최소가 매우 얕고, 반사광의 파장 의존성이 모델링에 큰 영향을 미친다. 기존 연구들(Foret et al. 2010, Lee et al. 2017 등)은 복사율(albedo)을 1을 초과하도록 허용해야만 관측된 광도곡선을 재현할 수 있었으며, 이는 모델의 물리적 한계와 반사광 처리 방식의 부정확성을 시사한다.

본 논문은 로머 지연(Rømer delay) 측정을 통해 1차와 2차 식이시각 사이의 시간 차이를 정밀히 구하고, 이를 질량비 q = Δt · c / a (여기서 a는 반지름) 형태의 식에 적용해 q = 0.299 ± 0.009라는 강력한 제약을 얻었다. 이 값은 이전 연구에서 제시된 q ≈ 0.28–0.30 범위와 일치하지만, 오차가 크게 감소해 후속 모델링의 자유도를 크게 줄인다.

PHOEBE 2.17을 이용한 동시 광·방사속도 모델링에서는 2차성의 복사율을 물리적 상한인 A₂ = 1.0으로 고정하였다. 알베도, 온도, 반지름 사이의 강한 퇴화 상관관계에도 불구하고, 다중색( SDSS‑u,g,r,i,z, Johnson‑VR, Bessel‑R 등) 광도곡선과 고해상도 방사속도 데이터를 동시에 피팅함으로써 T₂ ≈ 2360 K, R₂ ≈ 0.16 R☉라는 값이 도출되었다. 이는 이전 연구에서 제시된 T₂ ≈ 3000 K, R₂ ≈ 0.18 R☉보다 낮으며, sdB‑dM 시스템의 반사광 모델링에 있어 알베도 고정이 실제 물리적 파라미터를 더 정확히 추정하게 함을 보여준다.

ETV 분석에서는 284개의 신뢰성 높은 식이시각(주로 1차 최소)과 TESS 단기·장기 데이터(섹터 24–51, 77, 78)를 활용해 O–C 다이어그램을 구축하였다. LiTE(광행시효) 모델을 적용했을 때, 원형 궤도( a ≈ 0.64 au, M₃ ≈ 4.7 M_J)보다 이심률 e ≈ 0.46, 주기 P₃ ≈ 13 yr, 최소 질량 M₃ ≈ 1 M_J인 모델이 χ²와 BIC 기준에서 현저히 우수하였다. 이는 제3천체가 기존 제안보다 더 작은 질량을 가지고, 약간 타원형 궤도를 따라 움직이고 있음을 의미한다.

또한, 장기적인 식이시각 변동이 10 yr 수준의 주기성을 보이며, 이는 제3천체의 존재 외에도 자기활동에 의한 Applegate 메커니즘이나 공통외피 잔류 물질의 중력적 영향도 배제할 수 없다는 점을 논의한다. 그러나 현재 데이터와 모델링 결과는 목성 질량 제3천체가 가장 설득력 있는 설명임을 강하게 뒷받침한다.

결론적으로, 로머 지연을 통한 질량비 정확도 향상, 알베도 고정에 기반한 반사광 모델링 개선, 그리고 장기적인 TESS·지상 관측을 결합한 ETV 분석은 DD CrB와 유사한 HW Vir형 PCEB 시스템의 물리·역학적 특성을 이해하는 데 새로운 기준을 제시한다.