中国科学技术大学：《量子化学 Quantum Chemistry》课程教学资源（参考书籍）《物质结构导论》第二章多电子原子.pdf_大学文库

1 第二章多电子原子 ............................................................................................................................................... 41 2.1 变分法 .......................................................................................................................................................... 41 2.1.1 变分原理 ............................................................................................................................................... 41 2.1.2 变分法 ................................................................................................................................................... 41 2.2 氦原子基态的变分处理 .............................................................................................................................. 42 2.2.1 氦原子的 Schrödinger 方程 .................................................................................................................. 42 2.2.2 原子单位 ............................................................................................................................................... 43 2.2.3 单电子近似 ........................................................................................................................................... 43 2.2.4 反对称波函数和泡利(Pauli)原理 ......................................................................................................... 44 2.2.5 氦原子基态的变分处理 ....................................................................................................................... 44 2.3 自洽场方法 .................................................................................................................................................. 46 2.3.1 氦原子总能量的表达式 ....................................................................................................................... 46 2.3.2 哈特利-福克(Hartree-Fock)方程 .......................................................................................................... 48 2.4 中心力场近似 .............................................................................................................................................. 49 2.4.1 中心力场近似 ....................................................................................................................................... 49 2.4.2 屏蔽常数和轨道指数 ........................................................................................................................... 50 2.5 原子内电子的排布 ...................................................................................................................................... 52 2.5.1 Pauli 原理 .............................................................................................................................................. 52 2.5.2 能量最低原理 ....................................................................................................................................... 53 2.5.3 洪特(Hund)规则 .................................................................................................................................... 54 2.6 原子的状态和原子光谱项 .......................................................................................................................... 54 2.6.1 电子组态与原子状态 ........................................................................................................................... 54 2.6.2 原子光谱项 ........................................................................................................................................... 55 2.6.3 举例说明原子光谱项的写法 ............................................................................................................... 56 2.7 原子光谱 ...................................................................................................................................................... 58 2.7.1 原子发射光谱和原子吸收光谱 ........................................................................................................... 58 2.7.2 原子光谱项所对应的能级 ................................................................................................................... 59 2.7.3 原子光谱的选择定则 ........................................................................................................................... 60 2.8 定态微扰理论 .............................................................................................................................................. 60 2.8.1 非简并情况下的定态微扰理论 ........................................................................................................... 61 2.8.2 简并情况下的定态微扰理论 ............................................................................................................... 63 2.9 定态微扰理论的简单应用 .......................................................................................................................... 63 2.9.1 氦原子基态的微扰处理 ....................................................................................................................... 64 2.9.2 氢原子的一级斯达克(Stark)效应 ........................................................................................................ 65 习题 ....................................................................................................................................................................... 65

d da =0 0a2 Eo'(a1,a2,.,am)>=0 (2-6) -=0 da 由这一方程组可解出最佳参数a1,a2,,am·这组参数所对应的就是最佳的近似基态能量，而这组参数所对应的平o(a1,a2,.,am)就是满足能量最低条件的最佳近似基态波函数.但是用这个波函数计算其它力学量的平均值，其近似程度不一定最好. 如果尝试变分函数平与能量较低的各本征态平，都正交，即 ∫Ψ*Ψndr=0(=0,1,2,,n-1) (2-7) 则平n可用本征函数系{展开平n=∑c4 (2-8) 利用{分的正交归一性，得 c=∫g*yn'dr (2-9) 由(2-7)式得 C1=c2==Cm-1=0 (2-10) 所以展开式(2-8)是从i=n-1开始 4n=∑c4 (2-11) i=n ∫ΨhΨhdr 采用(2-1)式的证明方法，容易得到 KE> 了%4%dn≥E, (2-12) 上式表明，也可以用变分法计算激发态的近似能量和近似波函数。因为精确的华，常常很难得到，所以(2-7)式往往无法直接利用.如果用近似的'来代替平，则(2-12) 式不一定可靠，这使得用变分法计算激发态的近似能量和近似波函数比较困难.在后面我们会看到，在某些情况下，这种计算仍然是可行的： 2.2 氦原子基态的变分处理氦原子有两个电子，是最简单的多电子原子体系.本节所讨论的“原子单位”、“单电子近似”、“反对称波函数”等也适用于其它多电子原子和多电子分子. 2.2.1氦原子的Schrodinger方程 2 象对氢原子的处理一样，对于多电子体系的处理，也普遍采用 Borm-Oppenheimer近似，即把电子和原子核分开处理，我们暂不考虑氦原子核的运动，而仅处理氨原子中的两个电子相对于核的运动.这时描述氨原子状态的波函数F应该是两个电子的波函数，即 F=Fx1,1,1,x2,2,22FFr,2)=F1,2) (2-13) F的这三种形式都是等价的，都表示F是两个电子坐标的函数.F1,2)的物理意义是：电子1在1同时电子2在n的几率密度】把坐标原点选在氦原子核上（见图2.1），则氢原子的Schrodinger方程为图2.1氦原子的坐标【编-益明-琴-跨号 +一]中=E中 r1 r2 r12 (2-14) 式中m为电子的质量，-e为电子电荷，为Planck常数除以2π.上式的前两项分别是电子1和2的动能 42

42 ப ப௔భ =0 ப ப௔మ =0 (2-6) ………… ப ப௔೘ =0 由这一方程组可解出最佳参数 a1，a2，…，am．这组参数所对应的就是最佳的近似基态能量，而这组参数所对应的 Ψ0'(a1，a2，…，am)就是满足能量最低条件的最佳近似基态波函数．但是用这个波函数计算其它力学量的平均值，其近似程度不一定最好．如果尝试变分函数 Ψn'与能量较低的各本征态 Ψi都正交，即 ׬Ψi*Ψn'd߬=0 (i=0, 1, 2, …, nെ1) (2-7) 则 Ψn'可用本征函数系{Ψi}展开 Ψn'= i ciΨi (2-8) 利用{Ψi}的正交归一性，得 ci=׬Ψi*Ψn'd߬ (2-9) 由(2-7)式得 c1=c2=...=cn-1=0 (2-10) 所以展开式(2-8)是从 i=nെ1 开始 Ψn'= in ciΨi (2-11) 采用(2-1)式的证明方法，容易得到 = ೙అ ׬ ᇲ∗ு෡అ೙ ᇲୢఛ ೙అ ׬ ᇲ∗అ೙ ᇲୢఛ ≥En (2-12) 上式表明，也可以用变分法计算激发态的近似能量和近似波函数．因为精确的 Ψi 常常很难得到，所以(2-7)式往往无法直接利用．如果用近似的 Ψi'来代替 Ψi，则 (2-12) 式不一定可靠，这使得用变分法计算激发态的近似能量和近似波函数比较困难．在后面我们会看到，在某些情况下，这种计算仍然是可行的． 2.2 氦原子基态的变分处理氦原子有两个电子，是最简单的多电子原子体系．本节所讨论的“原子单位”、“单电子近似”、“反对称波函数”等也适用于其它多电子原子和多电子分子． 2.2.1 氦原子的 Schrödinger 方程象对氢原子的处理一样，对于多电子体系的处理，也普遍采用 Born-Oppenheimer 近似，即把电子和原子核分开处理，我们暂不考虑氦原子核的运动，而仅处理氦原子中的两个电子相对于核的运动．这时描述氦原子状态的波函数 F 应该是两个电子的波函数，即 F=F(x1, y1, z1, x2, y2, z2)=F(r1, r2)=F(1, 2) (2-13) F 的这三种形式都是等价的，都表示 F 是两个电子坐标的函数．|F(r1, r2)|2 的物理意义是：电子 1 在 r1 同时电子 2 在 r2 的几率密度．把坐标原点选在氦原子核上(见图 2.1)，则氦原子的 Schrödinger 方程为 [െ ԰2 2݉ ׏ଵ ଶ െ ԰2 2݉ ׏ଶ ଶ െ ܼ݁2 ݎ1 െ ܼ݁2 2ݎ + ௘మ ௥భమ ]ߔ=Eߔ) 2-14) 式中 m 为电子的质量，െe 为电子电荷， ԰为 Planck 常数除以 2ߨ.上式的前两项分别是电子 1 和 2 的动能图 2.1 氦原子的坐标 r1 r2 2ߠ r12 1 2 dv2 +

算符，第三项和第四项分别为电子1和2与核的作用势能算符：Z为核电荷数，对于氢原子，Z=2:最后一项是电子间的排斥能算符. 2.2.2原子单位对于多电子体系的计算常采用原子单位a.u.,在原子单位下 h=l.m=l,e=1 (2-15) 即角动量以方为单位，质量以电子质量为单位，电荷以电子电荷的绝对值为单位.而且长度以Bor半径ao h2 为单位 1a0= 0.529A (2-16) a2 能量以哈特利(Hartree)为单位 me4 IHartree-h272eV (2-17) 采用原子单位有两个优点：(1)可以简化数学表达式.(2)随着实验技术的不断提高，物理常数的测定日趋精确，因此采用原子单位可使计算结果不受当时实验水平的限制，在原子单位下，氦原子的Schrodinger方程变为【吸-吸-异-异+品6种 (2-18) 2.2.3单电子近似经Born-Oppenheimer近似后，在原子单位下的氦原子的Schrodinger方程虽然已经简化成(2-l8)式，但这仍然是含有六个坐标变量的偏微分方程，很难精确求解.困难在于在多电子体系的Schrodinger方程中含双电子排斥能项 r12=V(x1-x2)2+y1-y2)2+(Z1-22)7 (2-19) 这一项与两个电子的瞬时坐标有关，不能分离变量而化成象氢原子那样的单电子的Schrodinger方程，所以不得不引入新的近似一单电子近似. 单电子近似，或称为轨道近似，是处理多电子体系的基本近似.这个近似的思想是：每一个电子受其它电子的瞬时作用，看成是其它电子的平均势场的作用.这样，每个电子都在原子核和其它电子的平均势场中运动，因而每个电子都有自身的单电子波函数和单电子能量. 例如氨原子的基态，从无机化学的教科书中可知其电子组态为1s2,这实际上己经采用了单电子近似：每个电子都有其单电子波函数中中1，和单电子能量E1·若第一个电子的波函数记为中1：()，第二个电子的波函数记为中.("2)，则第二个电子的电荷密度p(2)应该等于其几率密度中.(22乘以电子的电荷e。 p(2)=el中1.(2) (2-20) 图2.1中d2内电子2的电荷与电子1的作用势能为 1 1 p(2)d2= 32 二e21φ1.(2)dy2 (2-21) 12 电子2分布在整个空间，所以电子1和电子2间的总排斥势能为 ∫e21s2n (2-22) T12 在原子单位下，e=l.电子1的Schrodinger方程为好-号+∫恤260AeM (2-23) 上式就是在单电子近似下的单电子Schrodinger方程，在单电子近似下，对每个电子都能得到一个单电子方程，这样就把对多电子体系的求解问题简化为对各单电子的求解问题. 由方程(2-23)所解出的单电子波函数中1(1)通常称为轨道，原子的单电子波函数称为原子轨道。分子的

43 算符，第三项和第四项分别为电子 1 和 2 与核的作用势能算符；Z 为核电荷数，对于氦原子，Z=2；最后一项是电子间的排斥能算符． 2.2.2 原子单位对于多电子体系的计算常采用原子单位 a.u．，在原子单位下 ħ=l, m=l, e=1 (2-15) 即角动量以 ħ 为单位，质量以电子质量为单位，电荷以电子电荷的绝对值为单位．而且长度以 Bohr 半径 a0 为单位 1a0= ԰మ ௠௔మ ≈0.529Å (2-16) 能量以哈特利(Hartree)为单位 1Hartree= ௠௘ర ԰మ ≈27.2eV (2-17) 采用原子单位有两个优点：(1)可以简化数学表达式．(2)随着实验技术的不断提高，物理常数的测定日趋精确，因此采用原子单位可使计算结果不受当时实验水平的限制．在原子单位下，氦原子的 Schrödinger 方程变为 [െ 1 2 ׏ଵ ଶ െ 1 2 ׏ଶ ଶ െ ܼ 1ݎ െ ܼ 2ݎ ൅ 1 12ݎ ]ߔ=Eߔ) 2-18) 2.2.3 单电子近似经 Born-Oppenheimer 近似后，在原子单位下的氦原子的 Schrödinger 方程虽然已经简化成(2-18)式，但这仍然是含有六个坐标变量的偏微分方程，很难精确求解．困难在于在多电子体系的 Schrödinger 方程中含双电子排斥能项 r12=ඥሺݔଵ െ ݔଶሻଶ ൅ ሺݕଵ െ ݕଶሻଶ ൅ ሺݖଵ െ ݖଶሻଶ (2-19) 这一项与两个电子的瞬时坐标有关，不能分离变量而化成象氢原子那样的单电子的 Schrödinger 方程，所以不得不引入新的近似——单电子近似．单电子近似，或称为轨道近似，是处理多电子体系的基本近似．这个近似的思想是；每一个电子受其它电子的瞬时作用，看成是其它电子的平均势场的作用．这样，每个电子都在原子核和其它电子的平均势场中运动，因而每个电子都有自身的单电子波函数和单电子能量．例如氦原子的基态，从无机化学的教科书中可知其电子组态为 1s 2 ，这实际上已经采用了单电子近似：每个电子都有其单电子波函数中 ϕ1s和单电子能量∈1s．若第一个电子的波函数记为 ϕ1s(r1)，第二个电子的波函数记为 ϕ1s(r2)，则第二个电子的电荷密度ߩ)2)应该等于其几率密度|ϕ1s(2)| 2 乘以电子的电荷 e。 ߩ)2) =e|ϕ1s(2)|2 (2-20) 图 2.1 中 dv2内电子 2 的电荷与电子 l 的作用势能为 ଵ ௥భమ eߩ)2)dv2= ଵ ௥భమ e 2 |ϕ1s(2)|2 dv2 (2-21) 电子 2 分布在整个空间，所以电子 1 和电子 2 间的总排斥势能为 ׬௘ మ|థభೞሺଶሻ|మ ௥భమ dv2 (2-22) 在原子单位下，e=1．电子 1 的 Schrödinger 方程为 [െ 1 2 ׏ଵ ଶ െ ܼ 1ݎ ൅ ׬ห߶1ݏ൫2൯ห2 12ݎ ]ϕ1s(1)=∈1sϕ1s(1) (2-23) 上式就是在单电子近似下的单电子 Schrödinger 方程，在单电子近似下，对每个电子都能得到一个单电子方程，这样就把对多电子体系的求解问题简化为对各单电子的求解问题．由方程(2-23)所解出的单电子波函数 ϕ1s(1)通常称为轨道，原子的单电子波函数称为原子轨道．分子的

单电子波函数称为分子轨道。所以单电子近似又称为轨道近似.方程(2-23)中的能量本征值1，称为轨道能量. 2.2.4反对称波函数和泡利(Pauli原理电子具有完全相同的质量、电荷等固有性质，在同样的条件下，各电子的行为是完全相同的，电子的这种本质上的不可区分性称为全同性，称电子为全同粒子，同样的，任何一种微观粒子都是全同粒子.按照波函数的统计解释，交换任意两个微观粒子的坐标，将不会影响其几率密度，即 (1,2,,i…,,NΦ(1,2，，，…d,,N (2-24) 由此可以得到 Φ(1,2，，i,…,…,F±(1,2，，j…,i,0 (2-25) 这就是说，由于微观粒子的全同性，对于交换任意两个粒子的坐标，其波函数或者是对称的（正号），或者是反对称的（负号）理论和实验都已证明：具有半整数自旋量子数的微观粒子，如电子（弓和质子（弓等，其波函数是反对称的，这样的粒子称为费米(Fermi)子：具有整数自旋量子数的粒子，如光子(s=l)等，其波函数是对称的，这样的粒子称为玻色(Bose)子. 在轨道近似下，N个电子体系的反对称波函数可用一个N阶行列式表示 pa(1)na(1)中a(2)na(2)…中a(N)na(N) 1 (1,2,,N户中b(1)nb(1)中b(2)n(2)…中b(N)b(N) NI (2-26) pn(1)mm(1)中n(2)nn(2)…中n(N)nn(N)l 这种行列式是由J.C.斯雷特(Slater)首先引入的，故称为Slater行列式，式中1/WW为归一化常数，中为单电子的空间波函数（轨道），η为单电子的自旋波函数，交换任意两个电子的坐标，相当于交换行列式中的任意两列，按行列式的性质，中应改变符号，因而Slater行列式可以作为多电子体系的反对称波函数.Slater 行列式表明，电子是不可区分的，不能指定一个电子占据某个特定的自旋轨道，而只能是所有的电子等全重地占据在这一组自旋轨道上。以后我们将会看到，多电子原子的自旋轨道也须象氢原子那样用四个量子数(，1，m和m,)表示，如中a(1)na(1)实际上是中naom.(1)na(msi).如果Slater行列式中有两行完全相同，则行列式为零，这表明两个电子占据完全相同的自旋轨道（即两个电子的四个量子数都相同）的几率为零，或者说两个电子不能处在同一状态，这就是大家熟悉的Pauli原理. 2.2.5氦原子基态的变分处理在多电子原子中，各电子的一组主量子数和角量子数1值所确定的电子状态称为电子组态，氨原子基态的电子组态为1s2.在轨道近似下，应该解单电子方程(2-23)以确定∈1，和中1。但方程(2-23)中的单电子算符 i--异+∫吧 (2-27) T12 含有待确定的轨道中1s,所以方程 h1φ1.(1Fe1,p1.(1) (2-28) 只能用自洽场方法（我们以后讨论）求解，下面我们先用变分法来处理这一问题考虑到中2是球对称的，所以电子2对电子1的排斥势能算符2业也仅是片的函数（因为电子 T12 2的坐标积分掉了)，这就是说，电子2对电子1的排斥势能也是球对称的.这球对称的排斥作用相当于削弱了原子核对电子1的作用，可看成是电子2对原子核的屏蔽，在这种理解的基础上，可作如下近似 1中2s(2)12 —d2≈- (2-29) T12 44

44 单电子波函数称为分子轨道。所以单电子近似又称为轨道近似．方程(2-23)中的能量本征值∈1s 称为轨道能量． 2.2.4 反对称波函数和泡利(Pauli)原理电子具有完全相同的质量、电荷等固有性质，在同样的条件下，各电子的行为是完全相同的，电子的这种本质上的不可区分性称为全同性，称电子为全同粒子．同样的，任何一种微观粒子都是全同粒子．按照波函数的统计解释，交换任意两个微观粒子的坐标，将不会影响其几率密度，即 |ߔ)1, 2, …, i, …,j, …, N)|2 =| ߔ)1, 2, …, j, …,i, …, N)|2 (2-24) 由此可以得到 ߔ)1, 2, …, i, …,j, …, N)=±ߔ)1, 2, …, j, …,i, …, N) (2-25) 这就是说，由于微观粒子的全同性，对于交换任意两个粒子的坐标，其波函数或者是对称的（正号），或者是反对称的（负号）．理论和实验都已证明：具有半整数自旋量子数的微观粒子，如电子(s= ଵ ଶ )和质子(s= ଵ ଶ )等，其波函数是反对称的，这样的粒子称为费米(Fermi)子；具有整数自旋量子数的粒子，如光子(s=1)等，其波函数是对称的，这样的粒子称为玻色(Bose)子．在轨道近似下，N 个电子体系的反对称波函数可用一个 N 阶行列式表示 ߔ)1, 2, …, N)= ଵ √ே!，，ተ ߶௔ሺ1ሻߟ௔ሺ1ሻ ߶௔ሺ2ሻߟ௔ሺ2ሻ ⋯ ߶௔ሺܰሻߟ௔ሺܰሻ ߶௕ሺ1ሻߟ௕ሺ1ሻ ߶௕ሺ2ሻߟ௕ሺ2ሻ ⋯ ߶௕ሺܰሻߟ௕ሺܰሻ ⋯ ⋯⋯ ⋯ ߶௡ሺ1ሻߟ௡ሺ1ሻ ߶௡ሺ2ሻߟ௡ሺ2ሻ ⋯ ߶௡ሺܰሻߟ௡ሺܰሻተ (2-26) 这种行列式是由 J．C．斯雷特(Slater)首先引入的，故称为 Slater 行列式，式中1⁄√ܰ!为归一化常数，ϕ 为单电子的空间波函数（轨道），ߟ为单电子的自旋波函数，交换任意两个电子的坐标，相当于交换行列式中的任意两列，按行列式的性质，ߔ应改变符号，因而 Slater 行列式可以作为多电子体系的反对称波函数．Slater 行列式表明，电子是不可区分的，不能指定一个电子占据某个特定的自旋轨道，而只能是所有的电子等全重地占据在这一组自旋轨道上。以后我们将会看到，多电子原子的自旋轨道也须象氢原子那样用四个量子数（n，l，m 和 ms）表示，如߶௔ሺ1ሻߟ௔ሺ1ሻ实际上是߶௡ೌ௟ೌ௠ೌሺ1ሻߟ௔ሺ݉௦ଵሻ．如果 Slater 行列式中有两行完全相同，则行列式为零，这表明两个电子占据完全相同的自旋轨道（即两个电子的四个量子数都相同）的几率为零，或者说两个电子不能处在同一状态，这就是大家熟悉的 Pauli 原理． 2.2.5 氦原子基态的变分处理在多电子原子中，各电子的一组主量子数 n 和角量子数 l 值所确定的电子状态称为电子组态，氦原子基态的电子组态为 1s 2 ．在轨道近似下，应该解单电子方程(2-23)以确定∈1s 和߶1s。但方程(2-23)中的单电子算符 ݄෠ ଵ=െ 1 2 ׏ଵ ଶ െ ܼ 1ݎ ൅ ׬ห߶1ݏ൫2൯ห2 12ݎ dv2 (2-27) 含有待确定的轨道߶1s，所以方程 ݄෠ ଵϕ1s(1)=∈1sϕ1s(1) (2-28) 只能用自洽场方法（我们以后讨论）求解，下面我们先用变分法来处理这一问题．考虑到 ϕ1s(2)是球对称的，所以电子 2 对电子 1 的排斥势能算符׬| థሺଶሻ|మ ௥భమ 也仅是 r1 的函数（因为电子 2 的坐标积分掉了），这就是说，电子 2 对电子 1 的排斥势能也是球对称的．这球对称的排斥作用相当于削弱了原子核对电子 1 的作用，可看成是电子 2 对原子核的屏蔽，在这种理解的基础上，可作如下近似 ׬| థమೞሺଶሻ|మ ௥భమ dv2ൎ ఙ ௥భ (2-29)

46 =2∈1s+(ߣെܼ)׬ |థభೞሺଵሻ|మ ௥భ ׬(ܼെߣ)+dv1 |థభೞሺଶሻ|మ ௥మ ׬+dv1 |థభೞሺଵሻ|మ|థభೞሺଶሻ|మ ௥భమ dv1dv2 =2∈1s+2(ߣെܼ)׬ ห߶1ݏሺ1ሻห 2 1ݎ dv1+׬ ห߶1ݏሺ1ሻห 2 ห߶1ݏሺ2ሻห 2 12ݎ dv1dv2 (2-40) 其中积分 ׬ ห߶1ݏሺ1ሻห 2 1ݎ ׬=dv1 1 ߣߨ 3 eെ2ݎߣ1 1ݎ ଵݎ ଶsinߠଵdݎଵdߠଵd߮ଵ=ߣ) 2-41) (2-40)式中的最后一个积分，涉及到对两个电子坐标的积分．对电子 2 坐标的积分可选用 r1 方向为 z 轴，得 ׬ ห߶1ݏሺ1ሻห 2 ห߶1ݏሺ2ሻห 2 12ݎ dv1dv2=׬ dv1 |థభೞሺଵሻ|మ|థభೞሺଶሻ|మ ௥భమ 2ݎ 2 sinߠ2dݎ2dߠ2d߮2 =׬ dv1|߶ଵ௦ሺ1ሻ| ׬ଶ 1 ߣߨ 3 eെ2ݎߣ2 ටݎ1 2൅ݎ2 2െ2ݎ1ݎ2cosߠ2 2ݎ 2sinߠ2dݎ2dߠ2d߮2=׬ dv1|߶ଵ௦ሺ1ሻ| ଶ2ߣଷ 1 1ݎ ׬eିଶఒ௥మr2[r1+r2െ|r1െr2|]dr2 =׬ dv1|߶ଵ௦ሺ1ሻ| ଶ4ߣଷ 1 1ݎ 2ݎ2ݎߣeെ2 ׬] 2 ௥భ ଵ dr2+[׬ eെ2ݎߣ2ݎ2 ஶ 2 ௥భ dr2 1 ׬= ߨ 3ߣ eെ2ݎߣ2 [െ 1 1ݎ eିଶఒ௥మ െ ߣeିଶఒ௥మ+ ଵ ௥భ ଵݎ [ ଶsinߠଵdݎଵdߠଵd߮ଵ= ହ ଼ ߣ) 2-42) 把(2-32)、 (2-41)和(2-42)式代入(2-40)式，得 =െߣଶ+2(ߣ െZ)ߣ+ ହ ଼ +ߣଶെ2Zߣ=ߣ ହ ଼ ߣ) 2-43) 按变分法，应该选择变分参数 l 使取极小值，即满足 பாబ ᇲ డఒ =0 (2-44) 得 ߣ=Zെ ହ ଵ଺ (2-45) 将上式代入(2-43)式，可得氦原子(Z=2)基态近似能量 =െ(Zെ ହ ଵ଺) 2 =െ2.85Hartree=െ74.5eV (2-46) 氦原子基态能量的实验值为െ79.0eV，可见用变分法求得的近似能量确实高于实验值．在前面我们曾提到 l 是有效核电荷，即 ߣ=Zെߪ) 2-47) 其中ߪ为屏蔽常数．由(2-45)式可知，氦原子 1s 电子之间的屏蔽常数为 ହ ଵ଺≈0.3．关于其它电子的屏蔽常数我们将在后面讨论． 2.3 自洽场方法哈特利-福克(Hartree-Fock)自洽场力法是目前计算原子轨道和分子轨道的最精确方法．本节从简单的例子开始，逐渐展开对自洽场方法的讨论． 2.3.1 氦原子总能量的表达式上节中我们曾得到了基态氦原子总能量的表达式(2-37)，为简单起见，令 =E0 െ 1 2 ׏ଵ ଶ െ ߣ 1ݎ =ܪ෡ଵ െ 1 2 ׏ଶ ଶ െ ߣ 2ݎ =ܪ෡ଶ ϕ1s=ϕ1 (2-48)

则(2-37式变为 0o2a,+,+7211b2dnd: =∫中(1)A1中1(1)dv+∫中1(2)A2中1(2)dv2+∫p1(1)冲1(2+一φ1(1)φ1(2)dv1d2 32 (2-49) 上式最后一个积分 aa=可b0r品rde (2-50) 称为Coulomb积分，表示两个电子间的Coulomb作用能，如果氦原子中的一个电子被激发到中2轨道且两个电子的自旋平行（都是α自旋），并设中2=2，则氨原子的该激发态的电子波函数为 1|p1(1)a(1)中1(2)a(2) 1,2o2a)a() φ2(2)a(2) (2-51) 该激发态的电子总能量为 E-J(.2)lH.+Fz+Fo(1.2xdtdrz -JIx(1):(2)-:(1)(2IIR,+Rz+1Jx(1):(2)-:(1):(2)ldv:dvzJla(DFl@(2)dm.dma(2-52) T12 由于自旋波函数是归一化的，所以上式关于自旋部分的积分等于1.把上式的积分展开总共有12项，但原子轨道是相互正交的，即 ∫中1(1)p2(1)dv1=∫中2(2)φ1(2)d2=0 (2-53) 共有四项等于零，而且两个电子是不可区分的，交换两个电子的坐标仅表明两个电子的状态互换，对能量没有影响，所以 ∫(1)A1中(1)dv,=∫(2)A2中1(2)dy2 ∫p2(1)A1p21)dv1=∫p2(2)A2p2(2)dv2 :(:(2Far.dr.-..(Fdrdvs (2-54) r12 这样，(2-52)式的积分展开只剩下四项： E=j中i1a,p1)d+fp:La,1)d+j1中i(1r,2p(2 Fdvdvz--了p11)p1)2pu(2pa2dd (2-55) 12 T12 其中最后一项 Kp:2∫（I1四2(2)pa2ind (2-56) 称为交换积分.在氦原子基态时，电子自旋相反，其总能量E[见(2-49)式]不含变换积分项：在两个电子自旋平行的激发态时，总能量E中含有交换积分.可见交换积分只有在同自旋的电子间才能存在，它表示了同自旋电子间的一种相互作用能一交换能，由(2-55)式可知，交换积分项是负的，即自旋相同电子间的交换作用可使总能量降低，这样，在其它条件允许的情况下，电子将尽量保持自旋平行.由此可解释Hud规则：1相同的简并轨道上的电子将尽量分占不同的轨道而保持自旋平行同自旋电子间的交换作用是一种量子效应，起源于微观粒子的全同性，可用电子的反对称波函数加以解释，对于氢原子的基电子态，两个电子的自旋相反，当两个电子坐标相同时，即=n2,由(2-34)式可知，其波函数为 1中(r1)a(1)中1(r1)a(2)1 1,2疙l0,B,B20 (2-57) 这说明两个自旋相反的电子是可以无限靠近的，对于氢原子的两个电子自旋相同的激发态，如果两个电子的坐标相同，即r=n,由(2-51)式可知 1p1(r1)a(1)中1(r1)a(2)儿1 (12)-(r)@(I)()a(2)(r2(ri)[a(1)a(2)-a(1)a(2)0 (2-58) 即同自旋的电子无限靠近的状态是不存在的，同自旋的电子倾向于互相回避，好象电子周围存在一个小空间(Fermi hole),其它自旋相同的电子无法进入这空间，致使同自旋电子间的静电排斥作用减弱，总能量降 9>

47 则(2-37)式变为 =׬ ϕ1(1)ϕ1(2)[ܪ෡ଵ+ܪ෡ଶ+ ଵ ௥భమ ]ϕ1(1)ϕ1(2)dv1dv2 =׬ ϕ1(1)ܪ෡ଵϕ1(1)dv1+׬ ϕ1(2)ܪ෡ଶϕ1(2)dv2+׬ ϕ1(1)ϕ1(2)+ ଵ ௥భమ ϕ1(1)ϕ1(2)dv1dv2 (2-49) 上式最后一个积分 ܬథభథమ=׬| ϕ1(1)|2 ଵ ௥భమ |ϕ1(2)|2 dv1dv2 (2-50) 称为 Coulomb 积分，表示两个电子间的 Coulomb 作用能．如果氦原子中的一个电子被激发到 ϕ2s 轨道且两个电子的自旋平行（都是ߙ自旋），并设 ϕ2s=ϕ2，则氦原子的该激发态的电子波函数为 ߔ)1, 2)= ଵ √ଶ ฬ ߶ଵሺ1ሻߙሺ1ሻ ߶ଵሺ2ሻߙሺ2ሻ ߶ଶሺ1ሻߙሺ1ሻ ߶ଶሺ2ሻߙሺ2ሻฬ (2-51) 该激发态的电子总能量为 E1=׬ ߔ)1, 2)[ܪ෡ଵ+ܪ෡ଶ+ ଵ ௥భమ ] 2 ߔ)1, 2)d߬ଵd߬ଶ =׬] ϕ1(1)ϕ2(2)െϕ2(1)ϕ1(2)][ܪ෡ଵ+ܪ෡ଶ+ ଵ ௥భమ ][ϕ1(1)ϕ2(2)െϕ2(1)ϕ1(2)]dv1dv2׬| ߙ)1)|2 |ߙ)2)|2 dms1dms2 (2-52) 由于自旋波函数是归一化的，所以上式关于自旋部分的积分等于 1．把上式的积分展开总共有 12 项，但原子轨道是相互正交的，即 ׬ ϕ1(1)ϕ2(1)dv1=׬ϕ2(2)ϕ1(2)dv2=0 (2-53) 共有四项等于零，而且两个电子是不可区分的，交换两个电子的坐标仅表明两个电子的状态互换，对能量没有影响，所以 ׬ ϕ1(1)ܪ෡ଵϕ1(1)dv1=׬ϕ1(2)ܪ෡ଶϕ1(2)dv2 ׬ ϕ2(1)ܪ෡ଵϕ2(1)dv1=׬ϕ2(2)ܪ෡ଶϕ2(2)dv2 ׬| ϕ1(1)|2 ଵ ௥భమ |ϕ2(2)|2 dv1dv2=׬| ϕ1(2)|2 ଵ ௥భమ |ϕ2(1)|2 dv1dv2 (2-54) 这样， (2-52)式的积分展开只剩下四项： E1=׬ ϕ1(1)ܪ෡ଵϕ1(1)dv1+׬ ϕ2(1)ܪ෡ଵϕ2(1)dv1+׬| ϕ1(1)|2 ଵ ௥భమ |ϕ2(2)|2 dv1dv2െ׬ ϕ1(1)ϕ2(1) ଵ ௥భమ ϕ1(2)ϕ2(2)dv1dv2 (2-55) 其中最后一项 ܭథభథమ=׬ ϕ1(1)ϕ2(1) ଵ ௥భమ ϕ1(2)ϕ2(2)dv1dv2 (2-56) 称为交换积分．在氦原子基态时，电子自旋相反，其总能量 E0[见(2-49)式]不含变换积分项；在两个电子自旋平行的激发态时，总能量 E1 中含有交换积分．可见交换积分只有在同自旋的电子间才能存在，它表示了同自旋电子间的一种相互作用能——交换能，由(2-55)式可知，交换积分项是负的，即自旋相同电子间的交换作用可使总能量降低，这样，在其它条件允许的情况下，电子将尽量保持自旋平行．由此可解释 Hund 规则：l 相同的简并轨道上的电子将尽量分占不同的轨道而保持自旋平行．同自旋电子间的交换作用是一种量子效应，起源于微观粒子的全同性，可用电子的反对称波函数加以解释，对于氦原子的基电子态，两个电子的自旋相反，当两个电子坐标相同时，即 r1=r2，由(2-34)式可知，其波函数为 ߔ0(1, 2)= ଵ √ଶ ฬ ߶ଵሺ࢘ଵሻߙሺ1ሻ ߶ଵሺ࢘ଵሻߙሺ2ሻ ߶ଶሺ࢘ଵሻߚሺ1ሻ ߶ଶሺ࢘ଵሻߚሺ2ሻฬ≠0 (2-57) 这说明两个自旋相反的电子是可以无限靠近的．对于氦原子的两个电子自旋相同的激发态，如果两个电子的坐标相同，即 r1=r2，由(2-51)式可知 ߔ1(1,2)= ଵ √ଶ ฬ ߶ଵሺ࢘ଵሻߙሺ1ሻ ߶ଵሺ࢘ଵሻߙሺ2ሻ ߶ଶሺ࢘ଵሻߙሺ1ሻ ߶ଶሺ࢘ଵሻߙሺ2ሻฬ= ଵ √ଶ ߶ଵሺ࢘ଵሻ߶ଶሺ࢘ଵሻ[ߙሺ1ሻߙሺ2ሻ െ ߙሺ1ሻߙሺ2ሻ]=0 (2-58) 即同自旋的电子无限靠近的状态是不存在的，同自旋的电子倾向于互相回避，好象电子周围存在一个小空间（Fermi hole），其它自旋相同的电子无法进入这空间，致使同自旋电子间的静电排斥作用减弱，总能量降

低 2.3.2哈特利-福克(Hartree-.Fock)方程下面把氢原子激发态总能量的表达式(2-55)推广到含有N个电子的原子体系的一般情况，设中和中为自旋轨道，则总能量为四a.①r专花.品4，刨2号20四o,0品Q四b②d (2-59) 其中“*”表示复共轭，最后一项求和只遍及同自旋的电子，不同自旋的电子因自旋正交而积分为零，因子1/2 的出现是由于每对电子的积分重复两次，原子轨道应满足正交归一化条件 ∫中*(1)φ1)dr1=6 (2-60) 设原子的总波函数为(1,2，，W,则总能量E=∫中*A中d (2-61) 如果原子轨道有一微小变化δ中，则中将有一微小变化6中，相应的能量E也会变化δE,即 E+6E-∫（中+6中）*A(中+6p)dr=∫中*A中dr+∫6Φ*H中dr+∫中*A6中dr+∫6中*i6中dr 上式的最后一项比起前三项来要小得多，可以忽略，利用(2-61)一式，得 6E=∫6中*H中dr+∫中*H6中dr=∫6中*A中dr+复共轭 (2-62) 按照变分原理，能量E的极小值应最接近真实的能量，在E取极小值时，d中所对应的dE等于零.但中的变化应受条件(2-60)式的限制，所以这是一个条件极值问题，可用拉格朗日(Lagrange)不定乘子法求E 的极小值，即 6[E+∑∑e⑥-∫中，*(1)4(1)dr1F0 (2-63) 其中E为不定乘子，将(2-59)式代入上式，许可利用(2-62)式的展开规则将(2-63)式展开帅t了64*0A0d+复共银+∑∑J6冲，地 +复共轭-∑∑∫6中，*(1)冲2)一中，*2)中(2)dt1dt2+复共轭-∑∑∈6冲，*1)中(1)dt1+复共轭 (2-64) 12 上式对Coulomb积分和交换积分的变分各有四项，因i和j是等价的而合并成两项，同时去掉了原来的因子 1/2.上式中的6中，(1)是独立变化的，上式等于零则要求6中，(1)的系数等于零，即园+Σ2a2ut,u小-Σ6随2r2et,-∑0 (2-65) T12 这个方程叫做Hartree-Fock方程，简称HF方程，可以证明，其中的∈，所构成的矩阵是Hermite矩阵，总可以找到一个U变换而将其变成对角矩阵，即∈Eδ (2-66) 这种变换是把{中}线性组合起来形成一组新的轨道{中？.在{中}下，∈具有对角形式.我们假定这种变换已经完成，并仍用{中}来表示变换后的轨道{中}，则HF方程为 a+∑422t,10∑4222t,1-e00 (2-66) 上式称为Hartree-Fock方程的正则形式.以后用到的HF方程都是这种正则方程.其中e,就是处在中，态的单电子能量本征值，或称为轨道能量，由方程(2-66)所解出的{中}可使总能量E取极小值，并满足正交归一化的条件(2-60)。 48

48 低． 2.3.2 哈特利-福克(Hartree-Fock)方程下面把氦原子激发态总能量的表达式(2-55)推广到含有 N 个电子的原子体系的一般情况，设 ϕi和 ϕj为自旋轨道，则总能量为 E= N i 1 ௜߶ ׬ ∗ሺ1ሻܪ෡ଵ߶௜ሺ1ሻd߬ଵ+ ଵ ଶ N i 1  N j 1 |ሺ1ሻ |߶௜׬ ଶ 1 ݎ12 ห߶௝ሺ2ሻห ଶ d߬ଵd߬ଶെଵ ଶ N i 1  N j 1 ௜߶ ׬ ∗ሺ1ሻ߶௝ሺ1ሻ 1 12ݎ ߶௝ ∗ሺ2ሻ߶௜ሺ2ሻd߬ଵd߬ଶ (2-59) 其中“*”表示复共轭，最后一项求和只遍及同自旋的电子，不同自旋的电子因自旋正交而积分为零，因子 1/2 的出现是由于每对电子的积分重复两次，原子轨道应满足正交归一化条件 ׬ ϕi*(1)ϕj(1)d߬ଵ=ߜij (2-60) 设原子的总波函数为ߔ)1, 2, …, N)，则总能量 E=׬ ߔ*ܪ෡ߔd߬ (2-61) 如果原子轨道有一微小变化ߜϕ，则ߔ将有一微小变化ߔߜ,相应的能量 E 也会变化ߜE，即 ߬dߔߜ෡ܪ*ߔߜ ׬+߬dߔߜ෡ܪ*ߔ ׬+߬dߔ෡ܪ*ߔߜ ׬+߬dߔ෡ܪ*ߔ ׬=߬d)ߔߜ+ߔ)෡ܪ*(ߔߜ+ߔ) ׬=Eߜ+E 上式的最后一项比起前三项来要小得多，可以忽略，利用(2-61)一式，得 (62-2 (复共轭߬ +dߔ෡ܪ*ߔߜ ׬≡߬dߔߜ෡ܪ*ߔ ׬+߬dߔ෡ܪ*ߔߜ ׬=Eߜ 按照变分原理，能量 E 的极小值应最接近真实的能量．在 E 取极小值时，dϕ 所对应的 dE 等于零．但中的变化应受条件(2-60)式的限制，所以这是一个条件极值问题，可用拉格朗日(Lagrange)不定乘子法求 E 的极小值，即 ߜ]E+ i  j ∈ij(ߜijെ׬ ϕi*(1)ϕj(1)d߬ଵ)]=0 (2-63) 其中∈ij 为不定乘子，将(2-59)式代入上式，许可利用(2-62)式的展开规则将(2-63)式展开 [ߜEെߜ i  j ∈ij׬ϕi*(1)ϕj(1)d߬ଵ]= i ׬ ߜϕi*(1)ܪ෡ଵϕi(1)d߬ଵ+复共轭+ i  j ׬ ߜϕi*(1)ϕi(1) ଵ ௥భమ ϕj*(2)ϕj(2)d߬ଵd߬ଶ +复共轭െ  i  j ׬ ߜϕi*(1)ϕj(2) ଵ ௥భమ ϕj*(2)ϕi(2)d߬ଵd߬ଶ+复共轭െ  i  j ∈ij׬ ߜϕi*(1)ϕi(1)d߬ଵ+复共轭 (2-64) 上式对 Coulomb 积分和交换积分的变分各有四项，因 i 和 j 是等价的而合并成两项，同时去掉了原来的因子 1/2．上式中的ߜϕi*(1)是独立变化的，上式等于零则要求ߜϕi*(1)的系数等于零，即 [ܪ෡ଵ+ i ׬ ϕj*(2) ଵ ௥భమ ϕj(2)d߬ଶ]ϕi(1)െ  j (2*(ϕjߜ ׬ ଵ ௥భమ ϕi(2)d߬ଶϕj(1)= j ∈ijϕj(1) (2-65) 这个方程叫做 Hartree-Fock 方程，简称 HF 方程，可以证明，其中的∈ij所构成的矩阵是 Hermite 矩阵，总可以找到一个 U 变换而将其变成对角矩阵，即 ∈ij=∈iߜij (2-66) 这种变换是把{ϕi}线性组合起来形成一组新的轨道{ϕi'}．在{ϕi'}下，∈ij 具有对角形式．我们假定这种变换已经完成，并仍用{ϕi}来表示变换后的轨道{ϕi'}，则 HF 方程为 [ܪ෡ଵ+ j ׬ ϕj*(2) ଵ ௥భమ ϕj(2)d߬ଶ]ϕi(1)െ  j ׬ ϕj*(2) ଵ ௥భమ ϕi(2)d߬ଶϕj(1)=∈iϕi(1) (2-66) 上式称为 Hartree-Fock 方程的正则形式．以后用到的 HF 方程都是这种正则方程．其中∈i 就是处在 ϕi 态的单电子能量本征值，或称为轨道能量，由方程(2-66)所解出的{ϕi}可使总能量 E 取极小值，并满足正交归一化的条件(2-60)

Coulomb算符j,的定义为 j0bI向e京2t, (2-67) 交换算符R的定义为 R01FU2722tra0 (2-67 Fock算符(I)的定义为 (1)=A+∑01)-（明 (2-68) 利用这三个算符，可将HF方程写成下面两种形式： F(1)b(1）=E中(1) (2-69) A1+∑j(1)-(1)]中(1)=e,φ1) (2-70) Fock算符中含有待确定的轨道{中}，所以HF方程只能用自恰场方法求解：先用一组试探轨道{中，)写出Fock算符)，解H亚方程得出一组新的轨道中，，再用{φ，写出新的四，再解H亚方程得到一组更新的{中，2}·.这种迭代过程反复进行，直到进一步迭代时，轨道的变化已经很小（在给定的阈值范围之内）为止，这时就称这些轨道与它们所产生的势场自洽 2.4 中心力场近似在上一节我们得到了HF方程(2-70).为了便于讨论多电子原子的HF方程求解问题，需要把HF方程写成更简单的形式.HF方程中的 a=-好-号 (2-71) 可令电子1的势能算符为1,0，p1),则 ,0,p1上，0，p- ∑)少-叨 11 (2-72) HF方程变为 [-2+V,0,pb1=e (2-73) 2.4.1中心力场近似对于多电子原子轨道中中(1)和：的计算，直接解这一方程仍然比较困难，需要进一步简化.可以证明，闭壳层组态或半满壳层组态，其电荷分布是球对称的，势场也是球对称的，在这种情况下，HF方程的求解问题要简单得多，对于一般的开壳层组态，势场不是球对称的，为使HF方程求解问题简化，一般都引入中心力场近似，把非球对称势场作球平均 Vi))no,0.dp 。sin81d01dp1 (2-74) 用球对称的)代替HF方程的V(r1,O1,p1),从而得到中心力场近似下的HF方程 [-+中1)=e,() (2-75) 在球对称势场下，中(I)具有径向函数Rnr)和球谐函数Ym(01,p1)乘积的形式 A1)=Rni(r1)Yi(01,1) (2-78) 把上式代入到HF方程(2-75)中，因')与角度变量(，P1)无关，所以可以象解氢原子那样分离变量.其中 Ym(B1,p1)所满足的方程就是氢原子波函数角度部分所满足的方程，Ym(日1，p1)与氢原子波函数角度部分完全相同，径向函数R()所满足的方程与氢原子的径向方程不同，这是因为')与氢原子的不同.在原子 49

49 Coulomb 算符ܬ መ ܬ 的定义为௝ መ ௝(1)ϕi(1)=[׬ ϕj*(2) ଵ ௥భమ ϕj(2)d߬ଶ]ϕi(1) (2-67) 交换算符ܭ௝෡的定义为 ܭ)௝෡1)ϕi(1)= [׬ ϕj*(2) ଵ ௥భమ ϕi(2)d߬ଶ]ϕi(1) (2-67) Fock 算符ܨ෠ ௜(1)的定义为 ܨ෠ ௜(1)=ܪ෡ଵ+ j ܬ] መ ௝(1)െܭ)௝෡1)] (2-68) 利用这三个算符，可将 HF 方程写成下面两种形式： ෠ܨ ௜(1)ϕi(1)=∈iϕi(1) (2-69) ܪ෡ଵ+ j ܬ] መ ௝(1)െܭ)௝෡1)]ϕi(1)=∈iϕi(1) (2-70) Fock 算符中含有待确定的轨道{ϕi}，所以 HF 方程只能用自恰场方法求解：先用一组试探轨道{ϕi (0)}写出 Fock 算符 ܨ෠ ௜ ሺ଴ሻ，解 HF 方程得出一组新的轨道{ϕi (1)}，再用{ϕi (1)}写出新的ܨ෠ ௜ ሺଵሻ，再解 HF 方程得到一组更新的{ϕi (2)}…．这种迭代过程反复进行，直到进一步迭代时，轨道的变化已经很小（在给定的阈值范围之内）为止，这时就称这些轨道与它们所产生的势场自洽． 2.4 中心力场近似在上一节我们得到了 HF 方程(2-70)．为了便于讨论多电子原子的 HF 方程求解问题，需要把 HF 方程写成更简单的形式．HF 方程中的 ܪ෡ଵ=െ 1 2 ׏ଵ ଶ െ ܼ 1ݎ (2-71) 可令电子 1 的势能算符为ܸ෠ ௜(r1, θ1, ߮1)，则 ܸ෠ ௜(r1, θ1, ߮1)=Vi(r1, θ1, ߮1)= െ ܼ 1ݎ j ܬ] መ ௝(1)െܭ)௝෡1)] (2-72) HF 方程变为 [െ 1 2 ׏ଵ ଶ+Vi(r1, θ1, ߮1)]ϕi(1)=∈iϕi(1) (2-73) 2.4.1 中心力场近似对于多电子原子轨道中 ϕi(1)和∈i 的计算，直接解这一方程仍然比较困难，需要进一步简化．可以证明，闭壳层组态或半满壳层组态，其电荷分布是球对称的，势场也是球对称的，在这种情况下，HF 方程的求解问题要简单得多，对于一般的开壳层组态，势场不是球对称的，为使 HF 方程求解问题简化，一般都引入中心力场近似，把非球对称势场作球平均 Vi(r1)=׬ ׬೔௏ ሺ௥భ,ఏభ,ఝభሻ ഏ బ మഏ బ ୱ୧୬ఏభୢఏభୢఝభ ׬ ׬ ୱ୧୬ఏభୢఏభୢఝభ ഏ బ మഏ బ (2-74) 用球对称的 Vi(r1)代替 HF 方程的ܸ௜ሺݎଵ, ߠଵ, ߮ଵሻ，从而得到中心力场近似下的 HF 方程 [െ 1 2 ׏ଵ ଶ+Vi(r1)]ϕi(1)=∈iϕi(1) (2-75) 在球对称势场下，ϕi(1)具有径向函数 Rnl(r1)和球谐函数 Ylm(ߠଵ, ߮ଵ)乘积的形式 ϕi(1)=Rnl(r1)Ylm(ߠଵ, ߮ଵ) (2-78) 把上式代入到 HF 方程(2-75)中，因 Vi(r1)与角度变量(ߠଵ, ߮ଵ)无关，所以可以象解氢原子那样分离变量．其中 Ylm(ߠଵ, ߮ଵ)所满足的方程就是氢原子波函数角度部分所满足的方程，Ylm(ߠଵ, ߮ଵ)与氢原子波函数角度部分完全相同，径向函数 Rnl(r1)所满足的方程与氢原子的径向方程不同，这是因为 Vi(r1)与氢原子的不同．在原子

中国科学技术大学：《量子化学 Quantum Chemistry》课程教学资源（参考书籍）《物质结构导论》第二章 多电子原子

中国科学技术大学：《量子化学 Quantum Chemistry》课程教学资源（参考书籍）《物质结构导论》第二章多电子原子