[논문 리뷰] Numerov numerical method applied to the Schrödinger equation

Francisco Caruso, V. Oguri|arXiv (Cornell University)|2014. 03. 17.

Quantum Mechanics and Non-Hermitian Physics참고 문헌 6인용 수 3

한 줄 요약

이 논문은 양자역학에서 고유값 문제를 해결하기 위해 Numerov 방법을 교육적 응용한 것이다. 특히 1차원 조화진동자와 수소원자에 대해 적용한다. 표준 Numerov 방법을 수정하여 슈뢰딩거 방정식에 존재하는 일阶 도함수 항을 다룰 수 있도록 확장하고, 수치적 결과를 해석적 해와 비교하여 정확도를 검증하며, 2차 유한차분 스킴을 사용하여 에너지 준위와 원형파동함수를 고정밀도로 계산한다.

ABSTRACT

In this paper it is shown how to solve numerically eigenvalue problems associated to second order linear ordinary differential equations, containing also terms which depend on the variable. A didactic presentation of the Numerov Method is given and, in the sequel, it is applied to two quantum non-relativistic problems with well known analytical solutions: the simple harmonic oscillator and the hydrogen atom. The numerical results are compared to those obtained analytically.

연구 동기 및 목표

• 양자역학에서 두 번째 차수 선형 상미분방정식을 해결하기 위한 Numerov 방법에 대한 교육적 소개를 제공하기 위해.
• 표준 Numerov 방법을 파동함수의 일阶 도함수를 포함하는 항을 다룰 수 있도록 확장하기 위해.
• 기본적인 양자역학 시스템인 조화진동자와 수소원자에 대해 정적 슈뢰딩거 방정식을 수치적으로 해결하기 위해.
• 수치적 고유값과 파동함수를 해석적 해와 비교하여 방법의 정확도를 검증하기 위해.
• 현대 양자역학 계산에서 이 방법의 강건성과 일반적 적용 가능성에 대한 시범을 보이기 위해.

제안 방법

• 수치적 미분을 피하기 위해 일阶 도함수를 직접 사용하지 않고, Numerov 방법을 두 번째 차수 유한차분 스킴으로 적용하여 슈뢰딩거 방정식을 이산화한다.
• 파동함수의 일阶 도함수에 의존하는 항을 포함하기 위해 수정된 Numerov 공식을 유도하여, 각운동량이 있는 수소원자에의 적용을 가능하게 한다.
• 세 점 스텐실을 사용하여 초기 조건에서부터 파동함수를 전진 및 후진으로 전파함으로써 안정성과 정확도를 확보한다.
• 공간 영역의 양 끝(예: 원점과 무한대)에 경계 조건을 적용하여 고유값 문제를 해결한다.
• 공개된 링크를 통해 이용 가능한 C++로 구현된 ROOT 프레임워크를 사용한다.
• 수치적 에너지 준위와 원형파동함수를 정확한 해석적 결과와 비교하여 방법을 검증한다.

실험 결과

연구 질문

• RQ1 일阶 도함수 항을 포함하도록 확장된 Numerov 방법이 1차원 조화진동자의 에너지 고유값을 얼마나 정확하게 계산할 수 있는가?
• RQ2 수정된 Numerov 방법이 3차원 수소원자의 원형파동함수와 에너지 준위를 정확히 재현할 수 있는가?
• RQ3 조화진동자와 수소원자에 대해 수치적 고유값과 해석적 해를 비교했을 때 발생하는 수치 오차는 얼마인가?
• RQ4 비영인 각운동량을 포함하는 슈뢰딩거 방정식의 특이 경계값 문제를 해결할 때 이 방법은 어떻게 성능을 발휘하는가?
• RQ5 Numerov 방법은 표준 교과서 사례를 넘어서 현대 양자역학 문제에 얼마나 일반화되고 적용될 수 있는가?

주요 결과

• 1차원 조화진동자와 수소원자(l=1)의 수치적 에너지 준위는 해석적 해와 매우 밀접하게 일치하며, 편차는 2% 미만이다.
• 수소원자(l=1)의 경우 계산된 에너지 준위는 E2 = -3.47 eV(해석적: -3.40 eV), E3 = -1.54 eV(해석적: -1.51 eV), E4 = -0.83 eV(해석적: -0.85 eV)였다.
• 수치적으로 구한 원형파동함수는 형태와 노드 구조에서 해석적 해와 정확히 일치하여 방법의 구현 정확성을 확인한다.
• 효과적 위치에너지와 원형 도함수를 포함시켜, 수정된 Numerov 방법이 수소원자 슈뢰딩거 방정식의 일阶 도함수 항을 성공적으로 처리하였다.
• 이 방법은 큰 공간 영역에서 제로 근처의 특이 행동을 포함하여도 두 번째 차수 수렴성과 안정성을 보였다.
• ROOT 프레임워크를 사용한 C++로의 구현은 교육적 및 연구적 응용에 효율적이고 적합하다.