单电子近似

 一点穿连浩动 两仪内反复阴阳

1.3. 单电子近似#

备注

单电子近似是我觉得大部分误差的罪恶之源,不过也是一切算法的开始。

由于多电子波函数不仅在分析求解时难以处理,在数值求解时也因为维度太高无法难以积分。因此,若想处理大一些的体系,则必须做一假设:多电子波函数可以由单电子波函数的乘积态表示.

\[ \ket{\Psi(1,2,\cdots,n)} = \ket{\phi_1(1)}\ket{\phi_2(2)}\cdots\ket{\phi_n(n)} \]

这个近似约束的范围很大,对比BOA,BOA某种程度上是physical的,但这个近似所引入的误差则必须通过平均场来弥补。因为这个近似假设电子和电子之间是彼此无关的,那么这种缺失的效应则要像BOA一样引入额外的势场来弥补。但是,BOA可以说核的质量很大,电子是一种平均效应,但电子彼此之间可以这样说吗?

目前,几乎所有主流的方法都要引入这个近似,而又通过不同的手段弥补缺失的效应。我希望在某一天能看到,单电子近似可以因为算力的提高而不被使用。

电子密度#

将这个式子代入多电子概率密度的定义式(1.2.1):

\[\begin{split} \begin{split} n(r) &= \sum_i n(r_i) \\ &= \sum_i \int (dr^n/dr_i) \Psi^*(\cdots,r,\cdots)\Psi(\cdots,r,\cdots) \\ &= \sum_i \phi_i^*(r)\phi_i(r) \\ \end{split} \end{split}\]

以上定义是在连续表象下做的,也可以推广到离散表象中\(\braket{\xi_a|\phi_i}=C_{ia}\)1这里认为像坐标表象一样,任何一个粒子的态都可以用这一组基分解。由于电子彼此之间没有什么差别,这种分解是合理的。:

\[ n(a) = \sum_i C_{ia}^*C_{ia} \]

密度矩阵#

同样可以得到密度矩阵(1.2.2)的表达式:

\[\begin{split} \begin{split} \gamma(r, r') &= \sum_i \gamma(r_i, r_i') \\ &= \sum_i \int (dr^n/dr_i) \Psi^*(\cdots,r,\cdots)\Psi(\cdots,r',\cdots) \\ &= \sum_i \phi_i^*(r)\phi_i(r') \end{split} \end{split}\]

在离散表象中:

\[ \gamma(a, b) = \sum_i C_{ia}^*C_{ib} \]