두 차원 뮐러 알고리즘과 물리학적 적용

초록

**
본 논문은 복소수 두 변수에 대한 비선형 방정식 시스템을 풀기 위해 뮐러 알고리즘을 2차원으로 일반화한 새로운 근 찾기 방법을 제안한다. 제안된 방법은 뉴턴·브로이든 방법과 비교했을 때 수렴 속도와 안정성에서 유사하거나 우수한 성능을 보이며, 특히 Heun 함수와 같이 복잡한 특수함수를 포함하는 시스템에서 뛰어난 효율성을 나타낸다. 이를 검증하기 위해 슈바르츠시드 블랙홀의 Regge‑Wheeler 방정식과 Kerr 블랙홀의 Teukolsky 방정식에서 발생하는 준정상모드(QNM) 계산에 적용하였다. 결과는 기존 문헌의 QNM 주파수와 높은 일치도를 보이며, 새로운 알고리즘이 물리학적 복소수 방정식 해결에 실용적임을 입증한다.

**

상세 분석

**
이 연구는 기존 1차원 뮐러 알고리즘을 2차원 복소수 공간으로 확장함으로써, 두 개의 비선형 방정식을 동시에 해결할 수 있는 프레임워크를 구축한다. 핵심 아이디어는 세 점을 이용해 2차 다항식 근사면을 만든 뒤, 그 근을 복소수 평면에서 반복적으로 업데이트하는 것이다. 이때 각 반복 단계에서 Jacobian 행렬을 직접 계산하지 않으므로, 뉴턴 방법이 요구하는 미분 연산 비용을 크게 절감한다. 또한, 브로이든 방법과 달리 근사식이 2차 형태이기 때문에 수렴 속도가 초기에 급격히 향상될 수 있다. 논문은 알고리즘의 수렴 조건을 상세히 분석하고, 복소수 변수의 실수·허수 부분을 동시에 다루는 복합적인 오류 전파 메커니즘을 제시한다. 특히 Heun 함수와 같이 특수함수의 수치적 평가가 비용이 크고 불안정한 경우, 기존 방법들은 미분값이 발산하거나 수치적 잡음에 민감해 실패하는 경우가 많다. 그러나 2차 뮐러 알고리즘은 함수값만을 이용해 근을 추정하므로, 미분이 필요 없는 상황에서도 안정적인 수렴을 보인다. 실험에서는 다양한 테스트 케이스(다항식, 초월함수, Heun 함수 기반 시스템)를 통해 뉴턴·브로이든과 비교했을 때 평균 반복 횟수와 계산 시간에서 경쟁력을 확인하였다. 특히 Heun 함수가 포함된 시스템에서는 뉴턴이 수렴하지 않거나 브로이든이 발산하는 경우, 2차 뮐러는 성공적으로 근을 찾아내어 그 우수성을 입증한다. 마지막으로, 블랙홀 QNM 계산에 적용한 결과는 기존 문헌값과 오차가 10⁻⁶ 이하로 일치함을 보여, 물리학적 복소수 방정식 해결에 실용적인 도구가 될 가능성을 시사한다.

**