Cusp条件总结笔记

引言

在量子力学和电子结构计算中，cusp条件用于确保波函数在存在库仑相互作用的系统中具有正确的行为，特别是在粒子非常接近时。这种条件旨在消除波函数在奇点处的发散问题，使波函数在数学上保持有限且物理上合理。

尖点（Cusp）的定义

在数学中，尖点（cusp）通常指函数在某些点处导数的不连续或函数曲线出现尖锐转折的地方。在波函数的讨论中，cusp用于描述在粒子之间相互作用势能存在奇点（如库仑势能）时，波函数在这些点的行为。

物理背景

在量子力学中，电子与原子核或两个电子之间的库仑相互作用在距离很近时表现出 ( \frac{1}{r} ) 的奇点。如果不对波函数进行合理的修改，这种奇点会导致波函数在这些点附近的行为不稳定，不符合实际物理情况。

Cusp条件的形式

电子接近原子核的情况

薛定谔方程在电子接近原子核时简化为：
$\left[- \frac{1}{2L^2} \nabla_e^2 - \frac{Ze^2}{rL}\right] \psi = E \psi \tag{C.1}$

其中，( Z ) 是核电荷，( L ) 是重标距离。将第一个项用球坐标表示，可以得到：
$-\frac{1}{2}\frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial \psi}{\partial r} \right) - \frac{Ze^2}{rL} \psi = E \psi \tag{C.2}$

考虑 ( r ) 很小时的情况，假设波函数的形式为 ( \psi = e^{-cr} )。代入上式：
$-\frac{1}{2} \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial e^{-cr}}{\partial r} \right) - \frac{Ze^2}{rL} e^{-cr} = E e^{-cr}$

计算导数：
$\frac{\partial e^{-cr}}{\partial r} = -ce^{-cr}$
$\frac{\partial}{\partial r} \left( r^2 \frac{\partial e^{-cr}}{\partial r} \right) = \frac{\partial}{\partial r} \left( -cr^2 e^{-cr} \right) = -c e^{-cr} \left( 2r - cr^2 \right)$

因此，有：
$-\frac{1}{2} \frac{1}{r^2} \left( -c e^{-cr} (2r - cr^2) \right) - \frac{Ze^2}{rL} e^{-cr} = E e^{-cr}$

简化上式：
$\frac{1}{2} \frac{c (2r - cr^2)}{r^2} e^{-cr} - \frac{Ze^2}{rL} e^{-cr} = E e^{-cr}$

在 ( r \to 0 ) 时，忽略高阶项：
$\frac{1}{2} \frac{2c}{r} e^{-cr} - \frac{Ze^2}{rL} e^{-cr} \approx 0$

为了消除奇点项，需满足：
$\frac{2c}{r} = \frac{2Ze^2}{rL}$

因此，有：
$c = Z L e^2$

两个电子接近的情况

当两个电子相互接近时，使用相对坐标 ( r_{12} = r_1 - r_2 ) 的薛定谔方程简化为：
$\left[- \frac{\nabla_{12}^2}{L^2} + \frac{e^2}{L r_{12}}\right] \psi = E \psi \tag{C.4}$

假设波函数的形式为 ( \psi = e^{-cr_{12}} )，代入上式：
$-\frac{1}{L^2} \nabla_{12}^2 e^{-cr_{12}} + \frac{e^2}{L r_{12}} e^{-cr_{12}} = E e^{-cr_{12}}$

计算导数：
$\nabla_{12} e^{-cr_{12}} = -c e^{-cr_{12}} \hat{r}_{12}$
$\nabla_{12}^2 e^{-cr_{12}} = c^2 e^{-cr_{12}}$

因此，有：
$-\frac{c^2}{L^2} e^{-cr_{12}} + \frac{e^2}{L r_{12}} e^{-cr_{12}} = E e^{-cr_{12}}$

简化上式：
$-\frac{c^2}{L^2} e^{-cr_{12}} + \frac{e^2}{L r_{12}} e^{-cr_{12}} = E e^{-cr_{12}}$

为了消除奇点项，需满足：
$\frac{e^2}{L r_{12}} = \frac{c^2}{L^2}$

因此，有：
$c = \frac{e^2 L}{2}$

对于反对称的情况（如两个电子的相对 ( p ) 态），cusp条件为：
$c = \frac{e^2 L}{4}$

应用cusp条件的原因

消除奇点：通过应用cusp条件，可以消除波函数在势能奇点处的发散问题，确保波函数在所有空间内是有限且连续的。
提高计算精度：在电子结构计算中，应用cusp条件可以显著提高计算结果的精度，因为它正确描述了电子在非常接近时的行为。
物理合理性：确保波函数在粒子相互作用强烈时的正确行为，使得计算结果更符合实际物理情况。

结论

cusp条件反映了波函数在处理库仑相互作用时的尖点行为。通过应用cusp条件，可以确保波函数在这些点处具有适当的行为，消除奇点引起的发散问题，从而得到物理上合理和数学上稳定的结果。这在量子力学和电子结构计算中具有重要意义。

参考Cusp conditions (attaccalite.com)

Cusp条件