쌍성의 운명: 트윈 별의 동적 반응과 선호도

초록

강한 1차 상전이가 발생하면 질량‑반지름 관계에 두 개의 안정적인 분기가 생긴다. 저자들은 일반 상대론적 수치 시뮬레이션을 통해 같은 질량을 갖는 하드론(HB)과 트윈(TB) 별이 외부 교란에 어떻게 반응하는지 조사했다. 교란 강도가 임계값 이하이면 원래 분기에서 진동하지만, 임계값을 초과하면 질량을 보존한 채 인접한 분기로 ‘이동’한다. 두 분기의 임계값 차이와 결합에너지 차이를 이용해 어느 분기가 자연계에서 더 선호되는지를 판단할 수 있다.

상세 분석

본 연구는 강한 1차 상전이가 존재하는 EOS를 사용해 트윈 별 현상을 정량적으로 분석한다. 저자들은 다중 구간 폴리트로픽 EOS를 채택했으며, 특히 전이 구간에서 음의 압력 구배를 피하기 위해 압력을 미세하게 증가시켜 음이 아닌 초소형 음속을 유지한다. 이렇게 구성된 EOS는 질량‑반지름(M‑R) 곡선에 두 개의 안정적인 분기, 즉 순수 하드론 분기(HB)와 퀘크코어를 가진 트윈 분기(TB)를 만든다. 동일한 정질량(M_b)을 갖는 두 별은 서로 다른 반지름과 중심 밀도를 가지며, 이는 각각 다른 내부 구조와 결합에너지 차이를 초래한다.

수치 실험은 1차원 GR1D 코드를 기반으로 수행했으며, 필요 시 3차원 WhiskyTHC 코드로 검증하였다. 초기 평형 모델에 일정한 내향 방사속도(v_r = –λ c) 를 부여해 교란을 가하고, λ의 크기를 변화시켜 임계 교란 λ_crit을 탐색한다. 결과는 두 분기 모두 선형 안정성을 보이지만, 비선형 교란에 대해서는 ‘임계 현상’이 존재함을 보여준다. λ < λ_crit이면 별은 원래 분기 내에서 작은 진동을 보이며 점차 감쇠하고, λ > λ_crit이면 별은 정질량을 유지한 채 중심 밀도가 급격히 변하면서 인접한 분기의 새로운 정적 해로 전이한다. 흥미롭게도 같은 정질량에 대해 HB와 TB의 λ_crit 값은 서로 다르며, λ_crit이 더 큰 쪽이 ‘선호된’ 분기로 정의될 수 있다.

또한 저자들은 두 분기의 결합에너지 차이 ΔE_bind을 계산해, ΔE_bind > 0(즉, TB가 더 낮은 결합에너지)인 경우 TB가 자연계에서 더 선호된다는 간단한 판정 기준을 제시한다. 이는 복잡한 시뮬레이션 없이도 EOS와 별 구조만으로 선호도를 예측할 수 있게 한다. 이러한 결과는 트윈 별의 비선형 안정성 문제를 해결하고, 관측적으로는 중성자별 질량‑반지름 측정이나 중력파 신호 해석에 중요한 함의를 제공한다.