One-loop self-energy using a numerical Green function

문제 정의
결합 전자 자기 에너지(self-energy, SE)는 원자 에너지 준위에 대한 가장 큰 양자 전기 역학(QED) 보정치를 나타낸다. 수소 유사 계량에서는 고정밀 계산이 존재하지만, 이를 가벼운 비수소 유사 계량(예: 분자 수소 이온)으로 확장하는 것은 상당한 과제를 안겨준다. 수소에 대한 기존의 고정밀 방법들은 해석적 그린 함수(Whittaker 함수)와 부분 파동 분해에 의존하며, 이는 막대한 계산 능력과 복잡한 수렴 가속 기술을 요구한다. 또한, 가벼운 계량에서 널리 사용되는 "포텐셜 전개법(potential-expansion method)"은 낮은 핵 전하($Z$) 값에 대해 수렴이 느리고 파인만 게이지(Feynman gauge)에서 큰 수치적 상쇄가 발생하는 문제가 있다. 분자 계량에서의 주요 장애물은 해석적인 디락-쿨롱 그린 함수의 부재이며, 이로 인해 유한 기저 집합으로부터 구축된 수치적 그린 함수를 사용해야 한다. B-스플라인(B-splines)이나 가우시안 기저를 사용한 이전의 수치적 그린 함수 시도들은 무거운 계량에 국한되었으며 5자리 정밀도를 넘어서지 못했다.

방법론
저자들은 지수 기저 집합 내에서 방사형 디락 방정식을 풀음으로써 얻은 수치적 그린 함수를 사용하여 수소 유사 원자의 1-루프 자기 에너지를 계산한다. 방법론은 다음과 같다:

1. 수치적 그린 함수 구축: 그린 함수는 유한 기저 집합으로부터 유도된 디락 고유 벡터와 고유값의 합으로 표현된다. 저자들은 향후 두 중심 포텐셜로의 확장을 고려하여, 이전에 분자 이온 연구에서 사용했던 것과 유사한 지수 기저 집합($\pi_\mu(x) = x^{\ell+1}e^{-\alpha_\mu x}$)을 채택하였다.
2. 포텐셜 전개 및 차감: SE 보정은 0-, 1-, 그리고 다중 포텐셜 항으로 분리된다. 발산하는 0- 및 1-포텐셜 항은 운동량 공간에서 재규격화되어 평가된다. 가장 수치적으로 까다로운 유한한 다중 포텐셜 항은 수치적 그린 함수를 사용하여 계산된다. 이 항을 격리하기 위해, 저자들은 전체 그린 함수로부터 (해석적으로 또는 동일한 기저를 통해 수치적으로 계산된) 0- 및 1-포텐셜 기여분을 차감한다.
3. 운동 균형(Kinetic Balance) 방식: 저자들은 기저 집합에 대해 세 가지 운동 균형 방식인 NKB(No Kinetic Balance), RKB(Restricted Kinetic Balance), DKB(Dual Kinetic Balance)를 비교한다. 저자들은 DKB가 정밀도를 크게 향방시키고, 특히 $\kappa > 0$인 경우 가상 상태(spurious states)를 제거한다는 것을 발견하였다.
4. 수렴 가속: 부분 파동 전개의 느린 수렴을 해결하기 위해, 저자들은 운동량 공간에서 계산 가능한 근사적인 다중 포텐셜 그린 함수를 차감하고 다시 더하는 차감 기법을 구현한다.
5. 게이지 선택: 계산은 파인만 게이지와 쿨롱 게이지 모두에서 수행된다. 낮은 $Z$에서 문제를 일으키는 수치적 상쇄를 완화하기 위해 쿨롱 게이지가 활용된다.
6. 수치적 구현: 모든 적분(방사형 및 에너지)은 고정밀 산술(MPFUN2020)과 특정 쿼드러처(quadrature) 기법(tanh-sinh 및 exp-sinh)을 사용하여 수치적으로 수행된다. 저자들은 1-포텐셜 그린 함수($G^1_\kappa$)를 계산하는 알고즘을 최적화하여 계산 복잡도를 $O(n_q^3 n_G^3)$에서 $O(\max(n_G, n_q)^4)$로 줄였다.

주요 결과

• 수소 유사 우라늄 ($Z=92$): DKB 지수 기저 집합을 사용하여, 저자들은 파인만 게이지에서 약 $10^{-5}$의 상대 불확도를 달면서 바닥 상태 자기 에너지를 달성하였다. 결과는 0-, 1-, 그리고 결합 상태 기여분 사이의 명확한 오차 상쇄를 보여주며, 이는 수치적 차감 접근법의 타당성을 입증한다. DKB 방식이 NKB 및 RKB보다 우수하여 정밀도를 2자리 수만큼 향상시킨다는 것이 확인되었다.
• 수소 원자 ($Z=1$):
  • 파인만 게이지: 수렴 가속 기법을 적용함으로써, 저자들은 수소 원자로 계산을 확장하여 총 SE 보정에 대해 $10^{-4}$의 상대 불측도를 달성하였다. 주요 제한 사항은 나머지 항의 큰 상대 불확도로 인한 무한 부분 파동으로의 외삽(extrapolation) 불안정성이었다.
  • 쿨롱 게이지: 쿨롱 게이지에서의 계산은 파인만 게이지에서 나타나는 수치적 상쇄를 피함으로써, 개별 부분 파동 기여도에 대해 더 높은 정밀도($10^{-5}$ ~ $10^{-6}$)를 확보하였다. 이는 총 SE 보정의 정밀도를 $10^{-5}$ 상대 불확도로 개선하였다.
• 기저 집합 성능: 지수 기저 집합은 저자들이 비선형 매개변수의 완전한 최적화를 수행하지 않았음에도 불구하고, 이전의 고정밀 연구에서 사용된 가우시안 기저 집합과 대등한 성능을 보였다.

의의 및 주장
본 논문은 순수 수치적 그린 함수 접근법을 사용하여 가벼운 일전자 계량에서 고정밀 자기 에너지 계산을 수행하는 것이 가능함을 입증한다. 저자들은 자신들의 방법이 해석적 그린 함수를 사용할 수 없는 분자 수소 이온과 같은 비수소 유사 계량으로 SE 계산을 확장하는 데 있어 실질적인 진전이라고 주장한다.

본 연구는 다음을 강조한다:

1. 지수 기저 내의 수치적 그린 함수는 무거운 원자에 대한 해석적 방법과 대등한 정밀도를 달성할 수 있으며, 수렴 가속 및 게이지 최적화와 결다면 가벼운 원자에서도 실행 가능하다.
2. 쿨롱 게이지는 낮은 $Z$ 시스템에서 수치적 상쇄를 줄임으로써 정밀도를 높이는 뚜렷한 이점을 제공한다.
3. 현재 접근법의 주요 제한 사항은 무한한 부분 파동으로의 외삽이며, 이것이 여전히 지배적인 불확도 원인으로 남아 있다.

저자들은 쿨롱 게이지의 이점을 완전히 활용하기 위해 기저 집합 매개변수의 추가적인 최적화가 필요하지만, 현재의 결과는 복잡한 분자 계량을 위한 향후 응용에 있어 수치적 프레임워크의 타당성을 검증한다고 결론짓는다.