에너지 모멘텀 제곱 중력에서 보편 관계를 이용한 중성자별 f모드 진동수 제한

초록

본 연구는 에너지‑모멘텀 제곱 중력(EMSG) 이론의 자유 파라미터 α가 중성자별(NS)의 구조와 비방사형 f‑모드 진동수에 미치는 영향을 조사한다. 수정된 토만‑오펜하이머‑볼크(TO‑V) 방정식을 세 가지 대표적인 방정식 상태(EOS) — Stiff, Intermediate, Soft—에 적용해 질량‑반경 관계와 사운드 스피드 프로파일을 계산하였다. Cowling 근사하에 비방사형 진동 방정식을 풀어 f‑모드 주파수를 얻고, 압축도(C)와 f‑모드, 그리고 Love‑수와 f‑모드 사이의 보편 관계(C‑f, f‑Love)를 구축하였다. GW170817·GW190814 등 관측 제약을 이용해 α 범위에 따른 f‑모드 주파수를 제한하고, 양의 α에서는 주파수가 감소, 음의 α에서는 증가함을 확인하였다.

상세 분석

EMSG는 라그랑지안에 α T_{\mu\nu}T^{\mu\nu} 항을 추가함으로써 에너지‑모멘텀 텐서의 제곱을 중력에 직접 결합한다. 이론적 전개에서 저자들은 L_m = P를 선택하고, 유효 에너지 밀도와 압력 E_{\rm eff}, P_{\rm eff} 을 각각 (11)·(12)식으로 정의하였다. 결과적으로 수정된 토만‑오펜하이머‑볼크 방정식(16)·(17)은 α에 따라 비선형적인 항을 포함하게 되며, 이는 고밀도 영역에서 특히 크게 작용한다.

수치해석에서는 세 가지 EOS(Stiff, Intermediate, Soft)를 사용해 중심 밀도 E_c 를 초기조건으로 설정하고, α를 -5×10^{-38}, -2.5×10^{-38}, 0, +2.5×10^{-38}, +5×10^{-38} (cm³/erg) 범위에서 변화를 탐색하였다. 질량‑반경 곡선은 α가 양수일 때 전반적으로 반경이 증가하고 최대 질량이 약간 감소하는 경향을 보였으며, 음수 α에서는 반대로 반경이 축소되고 최대 질량이 약간 상승하였다. 사운드 스피드 c_s² 프로파일은 모든 경우에 0 ≤ c_s² ≤ 1을 만족했으며, 특히 Stiff EOS에서는 높은 중심 밀도에서 c_s²가 1/3 근처의 임계값을 초과하지 않도록 α에 따라 미세 조정되는 ‘phase transition’ 구간이 나타났다. 이는 물리적 타당성을 확보하는 중요한 검증이다.

비방사형 진동은 Cowling 근사(시공간 메트릭 변동 무시) 하에 방정식(20)·(21)으로 기술되었다. 여기서 유효 압력·밀도는 α‑수정 항을 포함한 형태로 대체되었으며, 경계조건(23)·(24)을 적용해 l = 2( quadrupole) 모드의 고유주파수 ω를 구했다. 결과는 질량이 일정할 때 양의 α가 ω를 감소시키고, 음의 α가 ω를 증가시킨다. 특히, α = +5×10^{-38} 일 때 최대 안정 질량 근처에서 f‑모드 주파수가 약 1.8 kHz까지 낮아지는 반면, α = ‑5×10^{-38} 에서는 2.6 kHz 수준으로 상승하였다.

보편 관계 구축에서는 압축도 C = M/R과 f‑모드 주파수 사이의 C‑f 관계, 그리고 Love‑수(Λ)와 f‑모드 사이의 f‑Love 관계를 각각 다항식 형태로 피팅하였다. α에 따른 편차는 전체적으로 5 % 이내에 머물렀으며, 이는 기존 GR 기반 보편 관계와 크게 다르지 않음을 시사한다. 마지막으로 GW170817에서 얻은 Λ ≈ 300 ± 100 제한과 GW190814의 질량‑반경 추정치를 결합해, α가 0에 가까운 경우(즉, GR)에는 표준 f‑모드 주파수 ≈ 2.0 kHz가 가장 가능성이 높으며, α가 양(음)일 경우에는 1.6–2.4 kHz(2.2–2.8 kHz) 범위로 이동한다는 결론을 도출하였다.

전반적으로 이 연구는 EMSG가 고밀도 핵물리와 중성자별 진동학에 미치는 미세한 영향을 정량화하고, 관측 가능한 GW 데이터와 결합해 α의 허용 범위를 좁히는 데 기여한다는 점에서 의미가 크다.