[论文解读] Node-less atomic wave functions, Pauli repulsion and systematic projector augmentation

Peter E. Blöchl, Clemens J. Först|arXiv (Cornell University)|Oct 22, 2012

Cold Atom Physics and Bose-Einstein Condensates被引用 3

一句话总结

本文提出了一种新颖的方法，用于构建无节点的原子轨道波函数以及能量依赖、节点减少的分波，同时完整保留原子本征态的信息。通过逆解薛定谔方程，推导出描述核心电子泡利排斥作用的有效势能，从而在PAW方法中实现系统化的投影器增强方案，并提供一种透明且具有物理动机的替代传统赝势的方法，显著提升了关联体系的计算精度。

ABSTRACT

A construction of node-less atomic orbitals and energy-dependent, node-reduced partial waves is presented, that contains the full information of the atomic eigenstates and that allows to represent the scattering properties in a transparent manner. By inverting the defining Schrödinger equation, the Pauli repulsion by the core electrons can be represented as effective potential. This construction also provides a description of the Pauli repulsion by an environment. Furthermore, the representation leads to a new systematic scheme for a projector augmentation. The relation to Slater orbitals is discussed.

研究动机与目标

• 开发一种物理上透明的原子波函数表示方法，消除节点的同时完整保留本征态信息。
• 将核心电子引起的泡利排斥作用描述为一种有效、能量依赖、非局域的势能。
• 在投影增强波（PAW）方法中，实现一种系统化且具有物理基础的投影器增强方案。
• 构建一个框架，用于构造紧束缚轨道，并将强电子关联效应整合到基于DFT的模拟中。

提出的方法

• 通过关系式 $(\hat{H} - E_n)|u_n\rangle = -|u_{n-1}\rangle$ 递归构造无节点波函数，从基态开始。
• 利用能量差，推导标准束缚态 $|\psi_n\rangle$ 与无节点函数 $|u_n\rangle$ 之间的变换关系：$|\psi_n\rangle = \sum_{m=1}^n |u_m\rangle \prod_{j=1}^{m-1}(E_j - E_n)$。
• 通过逆解薛定谔方程，将泡利排斥作用表示为有效势能，得到一种非局域的赝势描述。
• 使用能量依赖的节点减少波函数 $|q_n(E)\rangle$，其定义为 $|q_n(E)\rangle = \frac{|q_{n-1}(E)\rangle - |q_{n-1}(E_{n-1})\rangle}{E - E_{n-1}}$，以确保在原点处的正确行为。
• 对非齐次径向薛定谔方程采用 $r$ 的幂级数展开，推导出系数 $a_j$ 的递推关系。
• 推导能量导数关系，表明 $|q_{n+j}(E_n)\rangle = \frac{1}{j!}|q_n^{(j)}(E_n)\rangle$，从而实现高阶投影器的系统化构造。

实验结果

研究问题

• RQ1 如何系统化地构造无节点原子波函数，使其在消除径向节点的同时完整保留原子本征态的全部信息？
• RQ2 能否通过逆解薛定谔方程，将核心电子引起的泡利排斥作用表示为一种有效、能量依赖、非局域的势能？
• RQ3 所得到的波函数如何实现PAW方法中系统化且物理透明的投影器增强方案？
• RQ4 这些无节点波函数与斯莱特型轨道或传统赝势之间，在数学和物理上存在何种联系？

主要发现

• 无节点波函数 $|u_n\rangle$ 通过能量差从束缚态 $|\psi_n\rangle$ 递归构造，可精确重构原始本征态。
• $|q_n(E)\rangle$ 在原点处的主导行为按 $r^{\ell + 2(n-1)}$ 变化，且系数与能量无关，证实了量子数增加时排斥作用增强。
• 在 $E_n$ 处对 $|q_n(E)\rangle$ 求能量导数，可得 $|q_{n+j}(E_n)\rangle = \frac{1}{j!}|q_n^{(j)}(E_n)\rangle$，从而实现系统化、递归式的投影器构造。
• 描述泡利排斥的有效势能直接由薛定谔方程导出，具有非局域性，避免了早期赝势方法中常见的发散问题。
• 该方法为紧束缚轨道提供了自然基底，并可通过具有物理基础的增强方案，将强电子关联效应整合到DFT框架中。
• 该形式化方法可推广至相对论体系，且与PAW方法保持一致，为标准赝势构建提供了一种透明且系统化的替代方案。