[논문 리뷰] Dynamical Rotational Instability at Low $T/W$

Joan Centrella, Kimberly C. B. New외 2명|2000. 10. 27.

Pulsars and Gravitational Waves Research인용 수 29

한 줄 요약

이 논문은 $T/|W| \gtrsim 0.14$인 비균일하게 회전하는 별핵에서 $m=1$ 모드가 지배하는 동역학적 회전 불안정성이 발생할 수 있으며, 이로 인해 중성자별 밀도로 빠르게 붕괴됨을 보여준다. 3D 뉴턴 유체역학 시뮬레이션을 통해 $N=3.33$ 다항체에서 $j$-일정 회전 법칙을 적용한 결과, 이 불안정성은 LIGO 또는 LISA에 의해 감지 가능한 중력파를 유도할 수 있다. 소스 질량 척도에 따라 감지 가능성이 달라진다.

ABSTRACT

Dynamical instability is shown to occur in differentially rotating polytropes with N = 3.33 and $T/|W| \gtrsim 0.14$. This instability has a strong m=1 mode, although the m=2, 3, and 4 modes also appear. Such instability may allow a centrifugally-hung core to begin collapsing to neutron star densities on a dynamical timescale. The gravitational radiation emitted by such unstable cores may be detectable with advanced ground-based detectors, such as LIGO II. If the instability occurs in a supermassive star, it may produce gravitational radiation detectable by the space-based detector LISA.

연구 동기 및 목표

• 비균일하게 회전하는 별핵에서 이전에 알려진 것보다 낮은 $T/|W|$ 값에서도 동역학적 회전 불안정성이 발생할 수 있는지 조사하기 위해.
• 이러한 불안정성이 원심력에 의해 끌어진 핵에서 중성자별 밀도로의 빠른 붕괴를 유도할 수 있는지 확인하기 위해.
• 이러한 불안정핵의 중력파 방출 특성을 분석하여 LIGO-II와 같은 고감도 지상 기반 및 LISA와 같은 궤도 기반 탐지기에서 감지 가능한지 평가하기 위해.
• $m=1$ 모드 지배의 불안정성 기구에서의 역할과 핵 진화에 대한 영향을 검토하기 위해.

제안 방법

• 뉴턴 체계에서 중력과 역반작용 없이, $N=3.33$ 및 $\Gamma=1.3$를 가진 비균일하게 회전하는 다항체 별핵의 수치 시뮬레이션을 수행하였다.
• 초기 모델은 $j$-일정 회전 법칙 $\Omega^2 = j_o^2 / (d^2 + \varpi^2)^2$ 를 사용한 자기일관성 필드 방법으로 구성되었으며, $d=0.2$를 사용하였다.
• 두 개의 독립된 3D 유체역학 코드를 사용하여 시뮬레이션을 수행하였다: 원통좌표계 $(\varpi,z,\phi)$ 그리드 코드와 직교좌표계 $(x,y,z)$ 그리드 코드로, 각각 $64\times64\times128$ 또는 $128\times128\times128$ 격자 영역을 사용하였다.
• 불안정성 성장의 씨앗이 되기 위해 1%의 무작위 밀도 편향을 도입하였으며, 격자 경계에서 질량 손실이 두드러지게 발생할 경우 시뮬레이션을 중단하였다.
• 가짜 $m=1$ 신호를 배제하기 위해 질량중심 운동을 모니터링하였으며, 수치 정확도 평가를 위해 각운동량 손실도 추적하였다.
• 초기 불안정 단계에서 중력파 진폭과 주파수를 추정하였으며, $\beta=0.14$일 경우 $h \sim 10^{-24} r_{20}^{-1}$, $\beta=0.18$일 경우 $h \sim 10^{-23} r_{20}^{-1}$로 추정하였다.

실험 결과

연구 질문

• RQ1 비균일하게 회전하는 별핵에서 $T/|W| \lesssim 0.27$ 에서도 동역학적 회전 불안정성이 발생할 수 있는가?
• RQ2 낮은 $T/|W|$ 조건에서 불안정성의 주요 애자이멀 모드($m$)는 무엇인가?
• RQ3 이러한 불안정성이 중성자별 밀도로의 빠른 붕괴를 유도할 수 있는가?
• RQ4 이러한 불안정성의 중력파 특징은 무엇이며, LIGO-II 또는 LISA에 의해 감지 가능한가?
• RQ5 질량 중심 이동이나 각운동량 손실과 같은 수치적 오류는 시뮬레이션 결과에 어떤 영향을 미치는가?

주요 결과

• $N=3.33$ 다항체에서 $T/|W| \gtrsim 0.14$일 때 동역학적 회전 불안정성이 발생하며, 이는 이론적으로 알려진 뉴턴 체계의 임계값 $\sim 0.27$보다 훨씬 낮은 값이다.
• 불안정성은 $m=1$ 모드에 의해 지배되며, $m=2$, $m=3$, $m=4$ 모드도 존재함을 확인하여 비축대칭적이고 비대칭적인 붕괴임을 시사한다.
• 불안정한 핵은 단일 밀도 높은 덩어리로 진화하여, 동역학적 시간 척도 내에서 중성자별 밀도로의 빠른 붕괴를 나타낸다.
• 200 km 반경의 $1.4 M_\odot$ 핵에서 중력파 진폭은 $\beta=0.14$일 경우 $h \sim 10^{-24} r_{20}^{-1}$, $\beta=0.18$일 경우 $h \sim 10^{-23} r_{20}^{-1}$로 추정되며, 이는 LIGO-II에 의해 감지 가능할 수 있다.
• 질량 중심 이동이 격자 경계에서 질량 손실로 인해 발생했지만, 두 코드 모두에서 1개 격자 이내의 수준으로 관측되었으며, 이는 $m=1$ 모드 식별에 영향을 주지 않았다.
• 초기 불안정 단계에서 중력파 진폭과 주파수를 추정한 결과, $\beta=0.14$일 경우 $h \sim 10^{-24} r_{20}^{-1}$, $\beta=0.18$일 경우 $h \sim 10^{-23} r_{20}^{-1}$로 추정하였다.