量子多体理论的Hartree.docxVIP

量子多体理论中的Hartree-Fock近似邢玉恒 (扬州大学物理科学与技术学院M14336) 我们知道，Hartree-Fock方程又简称为HF方程是一个应用变分法计算多电子体系波函数的方程，是量子化学中最重要的方程之一，基于分子轨道理论的所有量子化学计算方法都是以HF方程为基础的，鉴于分子轨道理论在现代量子化学中的广泛应用，HF方程可以被称作现代量子化学的基石。HF方程的基本思路为：多电子体系波函数是由体系分子轨道波函数为基础构造的斯莱特行列式，而体系分子轨道波函数是由体系中所有原子轨道波函数经过线性组合构成的，那么不改变方程中的算子和波函数形式，仅仅改变构成分子轨道的原子轨道波函数系数，便能使体系能量达到最低点，这一最低能量便是体系电子总能量的近似，而在这一点上获得的多电子体系波函数便是体系波函数的近似。 在量子力学中，关于相互作用的量子力学系集的讨论的出发点是哈密顿函数 八「八，、/)2V2,、八,、，?…一U=Jd3x中?(x)-%—+UCx)中(X)的基态19＞，在19＞中，N个粒子中每一个占据一个确定的单粒子态，因此它的运动与其他粒子无关。这个情况明显的将被粒子间的相互作用修改；尽管如此，实验事实是，对于许多不同的系统，例如:金属、原子和核子，这种单粒子描述方法是一个非常好的近似。因此一个自然的方法就是保留单粒子图像，而认为各个粒子是在一个单粒子势中运动，这个势来自这个粒子与其它粒子的平均相互作用。于是，这个单粒子能量应该是未微扰能量加上对所有其它粒子占据的态进行了平均的相互作用势能。这样，作为第一级近似，我们可以只保留对正规自能的第一级贡献，但是这样的计算不是完全自洽的。 考虑处于静态的且与自旋无关的外势场U(x)中的系统，例如，金属或原子中的电子，外势场破坏了空间均匀性，在绝热近似下，总的哈密顿量H=Jd3xS?(x)-切2+U(x)S(x)，0aL2m」a百=—Jd3xd3xS?(x)Sf(x)v(x-x)*W(x)V(x)12aPPa 这里为简单起见，认为粒子间相互作用势与自旋无关 心一气、呻，=心一X、%，相互作用的存在给求解带来困难，只能借助近似程序Schrodinger变分原理为在约束?g=1下求〈妇Hg的极值，即8?|HI?=0其解为HI?=EI?E为Lagrange乘子 求变分的严格极值等于解多体问题是无希望的，然而可以在某条件下选定的子空间上求变分极值，虽然结果不是严格的，但在数学上是可行的。 取一组正交归一完备的单粒子态I\，设试探函数?为I?=Ikkk=a+a+a+I0将Hamilton算符2-1也表为k的二次量子化表象H=ZkIuIka+a+1Ekk3(Ix-xI)Ikka+a+aalmlm211mlmlmml 1mlmlm其中在坐标表象中的方2—u=——V2+V(x)2m下面计算?IHI?，先计算矩阵元?Ia+aI?=8?Ia+a+aaI?=88—8811mmllllmmlmml于是〈妇Hg=ZvkIuIk+1￡vkkIvIkk-—^vkkIvIkkll2lmlm2lmml变分约束条件V妇?=1相当于VklIkjmfm上式乘以拉氏乘子，对v?IHI?变分，得Hartree-Fock方程uIb+Z(vbIvIbIb-vbIvIbIb=￡Iblmlmmmlll m换到坐标表象[-方2V2+V(X)]甲(x)+工jdx，甲*(X)v(x-xgkm(x帅ki(x)m=1 +Mjdx处*(x)v(x-x)甲(x)甲(x)=￡甲(x)kmklkmlklm=1 左端第一项是电子的动能加上晶体势，第二项是其它的Hartree-Fock粒子对匕粒子的平均库伦势，称直接库伦作用，第三项是泡利原理引起的交换库伦相互作用。方程是非线性的。非定域的，但是自伴的。在某些条件(如长程势)交换项是不重要的。可以略去，这就是哈特利近似。方程简化为方2[-2-V2+Vf(x)]中kl(x)=￡l^ki(x)其中有效势Vf(x)=V(x)+工ftpk(x，)I2v(x-x，)dx，m=1由于Hartree-Fock方程是自伴的，因而其解是正交的。但是方程并没有给出真正的定态，特别是E=v?IHI?w￡l l=1HI?^EI?参量￡l的意义可以从Koopmans定理看出：七表示从N个粒子体系中取出勺粒子时体系的能量降低值。 对于一般的外势UHartree-Fock方程很难解出，因为单粒子波函数七（x）和能量ej两者都必须自洽的定出，对于均匀系统，U（x）为零，正规自能取形式》（x-x，），这些方程就变的简单的多了，容易验证平面波中（x）=V2ok.x满足k自洽的要求，相应的自洽的单粒子能量成为6k=6o+力》k。 基态能量归结为简单的形式E=（2s+1）V（2兀）Jd3k--hLk）6（k-k）k2_|f 这表明，自能修改了