[논문 리뷰] Analytical solution of the Thomas-Fermi equation for atoms

M. Oulne|arXiv (Cornell University)|2005. 11. 02.

Quantum Mechanics and Non-Hermitian Physics인용 수 3

한 줄 요약

이 논문은 톰슨-페르미 방정식에 대한 새로운 3파라미터 해석적 시험 함수를 제안하며, 원자 잠재력과 이온화 에너지를 모델링하는 데 있어 정확도를 크게 향상시킨다. 라이츠 변분법을 사용하여, 0 ≤ x ≤ 10 범위에서 평균 오차가 0.28%에 불과하고, x=0에서 도함수를 2% 이내의 오차로 재현함으로써, 이온화 에너지 예측값이 이전의 어떤 해석적 근사보다 허트리-폭크 결과에 더 가까워진다.

ABSTRACT

An approximate analytical solution of the Thomas-Fermi equation for neutral atoms is obtained, using the Ritz variational method, which reproduces accurately the numerical solution, in the range $0\leq x\leq50$, and its derivative at $x=0$. The proposed solution is used to calculate the total ionization energies of heavy atoms. The obtained results are in good agreement with the Hartree-Fock ones and better than those obtained from previously proposed trial functions by other authors.

연구 동기 및 목표

• 톰슨-페르미 방정식에 대해 수치적 해와 물리적 경계 조건에 더 잘 부합하는 더 정확한 해석적 해를 개발하는 것.
• 기존의 시험 함수에 비해 무거운 원자에 대한 총 이온화 에너지 예측을 향상시키는 것.
• 원점(x=0)에서 잠재력의 도함수를 정확히 재현함으로써 물리적 일관성의 중요한 기준을 충족시키는 것.
• 이전의 1, 2, 3파라미터 형태보다 잠재력의 형태와 이온화 에너지 정확도 측면에서 뛰어난 변분 시험 함수를 제공하는 것.
• 무거운 원자(Z ≥ 92)에 대해 허트리-폭크 데이터를 기준 기준으로 사용하여 방법을 검증하는 것.

제안 방법

• 새로운 3파라미터 시험 함수를 제안: φ(x) = (1 + α√x + βx e^{-γ√x})² e^{-2α√x), x=0에서의 적절한 거동과 무한대에서의 감쇠를 반영하도록 설계됨.
• 라이츠 변분법을 적용하여 라그랑지안 함수 L(φ) = ∫₀^∞ F(φ, φ', x) dx 를 최소화함. F는 (1/2)(dφ/dx)² + (2/5)(φ^{5/2}/√x)로 정의됨.
• 입자 수 조건 ∫ρ dv = Z 는 최소화 과정에서 라그랑주 승수를 통해 강제됨.
• 변분 파라미터 α, β, γ 는 맵러의 리LEASE 9를 사용하여 수치적으로 최적화되어 에너지 함수를 최소화하면서 정규화 조건을 만족함.
• 해결된 해는 톰슨-페르미 방정식의 수치적 해와 x=0에서의 도함수 재현 정확도를 테스트함.
• 이온화 에너지는 유도된 x=0에서의 도함수를 공식 E = (12/7)(2/9π²)^{1/3} (dφ/dx)_{x=0} Z^{7/3} 를 통해 계산됨.

실험 결과

연구 질문

• RQ1 0 ≤ x ≤ 50 범위에서 기존의 수치적 해를 더 정확하게 재현할 수 있는 새로운 해석적 시험 함수를 구성할 수 있는가?
• RQ2 특히 이온화 에너지 계산에 중요한 역할을 하는 x=0에서의 도함수 예측 정확도가 이전의 시험 함수보다 향상되었는가?
• RQ3 무거운 원자(Z ≥ 92)의 총 이온화 에너지 예측 정확도가 허트리-폭크 결과와 비교하여 어떻게 되는가?
• RQ4 기존의 근사치들과 달리 원자 번호가 증가함에 따라 정확도를 유지하거나 향상시킬 수 있는가?
• RQ5 기존의 시험 함수들이 왜 해와 그 도함수를 충분한 정밀도로 재현하지 못하는지, 그리고 새로운 형태가 이러한 한계를 어떻게 극복하는가?

주요 결과

• 제안된 해는 0 ≤ x ≤ 10 범위의 47개 점에서 평균 오차가 0.28%에 불과하며, 이는 2파라미터 함수(Eq. 2)의 1.13% 오차와 비교해 크게 낮고, 이전의 다른 함수들은 2.5% 이상의 오차를 보였다.
• x=0에서의 해의 도함수는 -1.61623647이며, 수치적 값 -1.58807102과 비교해 상대 오차가 2% 미만이다.
• 새로운 함수(Eq. 21)로부터 도출된 이온화 에너지 예측값은 Z=100일 때 오차 4.0%, Z=106일 때 오차 2.0%로, 원자 번호가 증가함에 따라 점차 향상된다.
• 모든 이전의 시험 함수(Eq. 1, 2, 3, 4, 5)보다 우수하며, Eq. 5는 x=0에서 도함수가 무한대이므로 제외됨.
• Z=109일 때 신규 방법은 44,403 하트리의 값을 예측하여, 허트리-폭크 값 44,042보다 0.8% 뿐 높게 나와 초무거운 원자에 대해 매우 높은 정밀도를 보임.
• 잠재력의 형태와 원점에서의 도함수 정확도 측면에서 이 해가 중무거운 원소의 원자 성질 계산에 있어 열등한 해보다 뛰어나다.