[논문 리뷰] Hyperboloidal discontinuous time-symmetric numerical algorithm with higher order jumps for gravitational self-force computations in the time domain

Lidia J. Gomes Da Silva, Rodrigo Panosso Macedo외 4명|arXiv (Cornell University)|2023. 06. 22.

Radio Astronomy Observations and Technology인용 수 4

한 줄 요약

이 논문은 시간 도메인에서 중력 자기력 계산을 위한 고차수 점프를 가진 초구형 불연속 시간 대칭 수치 알고리즘을 제시한다. 이는 불연속 라그랑주 및 에르미트 보간을 활용하여 입자가 있는 점에서의 특이성에도 불구하고 고차수 정확도를 달성한다. 이 방법은 정확한 복사 플럭스 측정과 입자 근처의 장/도함수 평가를 가능하게 하며, 극단적으로 질량비가 큰 병합 시스템(EMRI)에서의 전체 자기일관된 궤도 진화에 적합함을 보여준다.

ABSTRACT

Within the next decade the Laser Interferometer Space Antenna (LISA) is due to be launched, providing the opportunity to extract physics from stellar objects and systems, such as extit{Extreme Mass Ratio Inspirals}, (EMRIs) otherwise undetectable to ground based interferometers and Pulsar Timing Arrays (PTA). Unlike previous sources detected by the currently available observational methods, these sources can extit{only} be simulated using an accurate computation of the gravitational self-force. Whereas the field has seen outstanding progress in the frequency domain, metric reconstruction and self-force calculations are still an open challenge in the time domain. Such computations would not only further corroborate frequency domain calculations and models, but also allow for full self-consistent evolution of the orbit under the effect of the self-force. Given we have extit{a priori} information about the local structure of the discontinuity at the particle, we will show how to construct discontinuous spatial and temporal discretisations by operating on discontinuous Lagrange and Hermite interpolation formulae and hence recover higher order accuracy. In this work we demonstrate how this technique in conjunction with well-suited gauge choice (hyperboloidal slicing) and numerical (discontinuous collocation with time symmetric) methods can provide a relatively simple method of lines numerical algorithm to the problem. This is the first of a series of papers studying the behaviour of a point-particle prescribing circular geodesic motion in Schwarzschild in the extit{time domain}. In this work we describe the numerical machinery necessary for these computations and show not only our work is capable of highly accurate flux radiation measurements but it also shows suitability for evaluation of the necessary field and it's derivatives at the particle limit.

연구 동기 및 목표

• 극단적으로 질량비가 큰 병합 시스템(EMRI)에서 중력 자기력 계산을 위한 안정적인 시간 도메인 수치 프레임워크를 개발하는 것.
• 주파수 도메인에서의 진전에도 불구하고, 메트릭 복원 및 입자 근처의 장 평가가 여전히 열려 있는 문제인 시간 도메인에서 자기력 계산의 정확성을 해결하는 것.
• 점입자에서의 특이성을 극복하기 위해, 공간 및 시간 이산화 과정에 불연속성의 구조에 대한 사전 지식을 통합하는 것.
• 초구형 절단과 시간 대칭 불연속 콜로케이션을 사용한 선형 방법 기반의 방법으로, 복사 플럭스 및 장 도함수 측정에서 고차수 정확도와 안정성을 입증하는 것.

제안 방법

• 초구형 절단을 사용하여 미래의 영원한 무한대에서 복사파를 자연스럽게 흡수함으로써, 시간 도메인 시뮬레이션에서 장기적 안정성을 향상시킨다.
• 불연속 라그랑주 및 에르미트 보간 공식을 사용하여 입자 근처의 불연속성의 국소적 구조를 포착하는 공간 및 시간 이산화를 구성한다.
• 장의 특이적 행동으로부터 유도된 점프 조건을 보간 프레임워크에 통합함으로써 고차수 정확도를 회복한다.
• 시간 대칭 불연속 콜로케이션 방법을 선형 방법 기반의 수식에 적용하여 안정성 및 수렴 성질을 향상시킨다.
• 기본 모델로 EMRI 웨이브폼에 사용되는 슈바르츠실트 시공간 내 원운동을 하는 점입자에 대해 알고리즘을 구현한다.
• 수치적 검증을 위해 [39]의 주파수 도메인 데이터와 비교하여 시간 도메인 결과의 RWZ 주요 함수 및 그 도함수를 검증하고, 오차 추정치는 식 (76)을 통해 계산한다.

실험 결과

연구 질문

• RQ1 점입자 특이성이 존재하는 상황에서도 시간 도메인 수치 알고리즘이 중력 자기력 계산에서 고차수 정확도를 달성할 수 있는가?
• RQ2 장이 입자 위치에서 점프를 보일 경우, 불연속 공간 및 시간 이산화를 어떻게 구성하여 정확도를 유지할 수 있는가?
• RQ3 초구형 절단이 장기적인 시간 도메인 자기력 시뮬레이션의 안정성과 정확도를 어느 정도 향상시키는가?
• RQ4 이 방법은 플럭스 및 입자 근처의 장 도함수를 정확하게 계산하여 자기일관된 궤도 진화를 가능하게 할 수 있는가?
• RQ5 제안된 알고리즘은 RWZ 주요 함수에 대해 수렴성과 수치 오차 측면에서 주파수 도메인 결과와 어떻게 비교되는가?

주요 결과

• RWZ 주요 함수의 내부 및 외부 해에 대해 고차수 수렴을 달성하였으며, 특정 모드에서는 수치 오차가 $1.7 \times 10^{-9}$까지 낮아졌다.
• 복사 플럭스 측정이 매우 정확했으며, 첫 번째 $l=5$ 모드에서 오차 범위는 $[1.7 \times 10^{-9}, 7.9 \times 10^{-8}]$이었다.
• 입자 위치에서 RWZ 함수의 일阶 시간 및 반경 방향 도함수를 성공적으로 계산하였으며, 대부분의 모드에서 오차는 $10^{-8}$ 이하였다.
• 위상도에서 $(l,m)=(2,2)$ 모드의 경우 $\tau=10,000$에서 해밀토니안의 구조 보존이 관찰되어 장기적 안정성이 확인되었다.
• 수치 수렴 연구를 통해 이론적 기대에 부합하는 오차 스케일링을 보이며 최적의 시간 이산화 단계 선택이 확인되었다.
• 이 방법은 자기력의 입자 위치에서의 평가와의 호환성을 보이며, EMRI 모델에서 자기일관된 궤도 진화에 필수적인 조건을 충족한다.