중력상대론적 회전·축대칭 중성자별의 종횡 순환 모델

초록

이 논문은 일반 상대성 이론 하에서 종횡 순환을 포함한 정적·축대칭 중성자별 모델을 최초로 구축한다. 새로운 코드 GRNS를 이용해 고정점 반복법으로 해를 구하고, 순환 흐름의 와류 수에 따라 구분되는 2차원 모드 집합을 발견하였다. 최대 순환 속도 1000 km/s에서 표면 변형은 약 1 % 수준이며, 회전과 순환이 동시에 존재하는 경우에도 표면 회전이 없고 회전 속도가 순환 속도를 초과하지 않는다.

상세 분석

본 연구는 기존의 회전 중성자별 모델에 종횡(meridional) 순환을 추가함으로써, 중성자별 내부 유체역학을 보다 현실적으로 기술한다. 일반 상대성 이론의 3+1 분해와 축대칭성을 가정하여, 메트릭은 두 개의 잠재함수(중력 퍼텐셜)와 한 개의 회전 퍼텐셜으로 기술된다. 유체는 바리스틱 방정식과 연속 방정식을 만족하며, 순환 흐름은 스트림 함수 ψ로 표현된다. ψ는 라플라시안 형태의 비선형 편미분 방정식을 만족하고, 경계 조건은 중심에서의 정칙성 및 표면에서의 압력 제로 조건으로 설정된다.

코드 GRNS는 이 복합 방정식계를 고정점 반복법으로 풀며, 초기값으로는 기존의 회전 없는 토르소이드 모델을 사용한다. 반복 과정에서 ψ와 메트릭 퍼텐셜을 교대로 업데이트하고, 수렴 기준은 전체 에너지와 질량의 변화율이 10⁻⁸ 이하가 될 때이다. 수치 실험에서는 ψ의 고유 모드가 정수 n에 따라 구분되는 일련의 와류 구조를 형성함을 확인했다. n=1은 단일 와류, n=2 이상은 다중 와류가 표면 근처까지 확장된다.

물리적 결과로는 순환 속도가 10⁸ cm s⁻¹(≈1000 km s⁻¹) 정도일 때, 와류가 표면을 비대칭적으로 팽창·수축시켜 반경 변형이 0.5~1 %에 달한다는 점이다. 변형 정도는 와류가 표면에 가까울수록, 그리고 순환 속도가 클수록 크게 증가한다. 또한 회전과 순환을 동시에 고려한 모델에서는, 회전 속도가 순환 속도를 초과하지 않으며, 표면 자체는 회전하지 않는다. 이는 내부 흐름이 강한 순환에 의해 지배되며, 회전은 주로 내부 코어에서만 나타난다는 해석을 가능하게 한다.

이러한 결과는 중성자별의 비축대칭 변형이 중력파 방출이나 펄서 타이밍 노이즈에 미치는 영향을 재평가할 필요성을 제시한다. 특히, 와류 구조가 장기간 유지될 경우, 관측 가능한 전자기 신호의 주기적 변동이나 중력파 스펙트럼에 미세한 변조를 일으킬 수 있다. 향후 연구에서는 보다 복잡한 방정식 상태(예: 초전도·초유체 핵)와 자기장 효과를 포함시켜, 실제 천체에서의 적용 가능성을 검증해야 할 것이다.