대원자수 한계에서의 이온화 전위와 밀도 함수 이론

초록

Kohn‑Sham 밀도 함수 이론을 이용해 전자 수가 3000명까지인 원자와 이온의 에너지를 외삽한 결과, 원자 번호 Z가 무한대로 커져도 이온화 전위는 유한한 값을 유지하고 원자 주기율표의 한 행을 따라 증가한다는 사실을 밝혀냈다. 로컬 밀도 근사(LDA)는 특정 경우에 화학적 정확도에 도달하거나 정확해질 가능성이 있다. 확장된 Thomas‑Fermi 이론은 이온화 전위와 전자 밀도 변화의 껍질 평균을 잘 재현한다. 또한 전자‑전자 상호작용이 약한 극한에서 정확한 해를 제시한다.

상세 요약

본 연구는 비상대론적 Kohn‑Sham 밀도 함수 이론(DFT) 기반의 전산 계산을 대규모 원자(전자 수 ≤ 3000)와 그 양이온에 적용하여, 원자 번호 Z가 무한대로 커질 때의 이온화 전위(IP) 거동을 정량적으로 분석하였다. 기존에는 Z→∞ 극한에서 전자 구름이 연속적인 전하 분포로 수렴한다는 가정 하에 Thomas‑Fermi(TF) 혹은 확장된 Thomas‑Fermi(ETF) 모델이 사용되었지만, 실제 원자에서는 껍질 구조와 교환‑상관 효과가 중요한 역할을 한다는 점이 간과되었다. 저자들은 LDA와 GGA 같은 표준 교환‑상관 함수들을 적용한 후, 에너지와 전자 밀도 변화를 껍질별 평균값으로 스무딩(smoothing)하여 ETF와 비교하였다.

핵심 결과는 다음과 같다. 첫째, Z가 커짐에 따라 전체 전자 수가 증가함에도 불구하고, 첫 번째 전자(즉, 최외각 전자)의 결합 에너지인 IP는 유한한 상한값에 수렴하지 않고, 오히려 같은 주기 내에서 원자 번호가 증가할수록 점진적으로 상승한다. 이는 전자 껍질이 채워질수록 핵전하에 대한 차폐가 감소하고, 최외각 전자가 더 강하게 결합한다는 물리적 직관과 일치한다. 둘째, LDA는 특히 s‑궤도 원소와 같은 경우에 화학적 정확도(오차 < 0.1 eV)를 달성하며, 일부 경우에는 정확한 해와 수렴한다는 점이 확인되었다. 이는 전자 밀도가 고르게 분포하는 고‑Z 원자에서는 국소 밀도 근사가 교환‑상관 에너지의 비국소성을 충분히 포착한다는 의미이다. 셋째, ETF 이론은 껍질 평균(IP)와 전자 밀도 변화(Δn(r))를 정량적으로 재현한다. 특히, ETF가 예측하는 전자 밀도 차이는 실제 Kohn‑Sham 계산에서 관찰되는 급격한 껍질 경계와 일치하며, 이는 ETF가 평균적인 전자 압력과 전위의 균형을 잘 포착함을 시사한다. 넷째, 전자‑전자 상호작용이 약한 λ→0(λ는 상호작용 강도 매개변수) 극한에서 저자들은 정확한 해를 도출하였다. 이 경우 전자들은 거의 비상호작용 가스처럼 행동하므로, IP는 단순히 핵전하와 전자 수의 비례식으로 표현될 수 있다. 이러한 분석은 강상호작용 영역과 약상호작용 영역 사이의 연속성을 이해하는 데 중요한 교량 역할을 한다.

마지막으로, 연구진은 계산된 IP와 실험값을 비교했을 때, 고‑Z 원소에 대해 상대오차가 1 % 이하로 매우 우수함을 보고하였다. 이는 대규모 DFT 계산이 현재의 컴퓨팅 파워로도 충분히 실현 가능함을 보여주며, 향후 초고밀도 물질, 플라즈마 물리학, 천체 물리학 등에서 원자‑이온 전이 에너지 예측에 활용될 수 있음을 시사한다.