多电子原子的结构

陈　兰　孙宏伟　赖城明

（南开大学化学学院，天津300071）



·自学之友·

多电子原子的结构

陈兰孙宏伟*赖城明

（南开大学化学学院，天津300071）

论述了不同精度Hamilton算符下多电子原子的结构，并对Slater行列式、电子组态、原子光谱项和光谱支项等相关内容进行了讨论。

结构化学；多电子原子

结构化学中无论研究原子还是分子总是从求解定态Schrödinger方程Hˆψ=Eψ入手，并以方程的解E和ψ表示体系的能量与状态。为什么用ψ可以表示体系的状态？大多数学生在学习了结构化学之后仍回答不了这个问题；而对多电子原子结构中原子光谱项和光谱支项的掌握则只限于会推求而并不了解它们的意义。本文从为什么求解Schrödinger方程谈起，在采用不同精度的Hamilton算符下，对Slater行列式、电子组态、原子光谱项和光谱支项等相关内容进行了讨论，希望能对学生系统地认识多电子原子的结构有所帮助。

1　定态Schrödinger方程的解ψ?与原子的状态

量子力学中假设微观体系的状态用波函数Ψ（r→，t）描述，而状态随时间的变化则服从含时间的Schrödinger方程［1］：

式（2）左右两边同时除以ψ（）r→f（）

t，有：

式（3）左边仅与时间有关，而右边仅与坐标有关，只有当左右两边都等于同一常数时，等式才能成立。以E表示这个常数，由式（3）左边等于E有：

对式（4）两边积分，可解得f（）t，见式（5），其中C为积分常数：

而由式（1）右边等于E有：

式（6）是Hamilton算符的本征方程，称之为定态Schrödinger方程（或与时间无关的Schrödinger方程）。由该方程的解ψ）乘以f（t）即得到波函数Ψ）：

综上，对势能函数不显含时间的体系，通过求解定态Schrödinger方程，得到E与ψ），E即体系的能量，而ψ）乘以f（t）即体系的状态。因为由Ψ）或ψ）（均已归一化）求得的粒子在空间分布的概率、体系的能量及其他可观测的物理量都是相同的，所以，这类体系的状态也可以简单地用ψ）来表示。这就解释了为什么在研究原子和分子时总是从求解定态Schrödinger方程出发，因为这些体系都可以看成是由原子核和电子组成的，带电粒子之间的相互作用为库仑引力，只与粒子的相对位置有关，体系总的势能函数不显含时间（不显含时间并不表示不随时间变化，势能通常与粒子的坐标有关，坐标随时间变化时势能也会随时间变化，但是当坐标复原时势能也能复原，这就是不显含时间或与时间无关），所以，通过求解定态Schrödinger方程即可确定体系的能量与状态，进而可以研究体系其他的物理量与性质。

2　采用不同精度的Hamilton算符下多电子原子的结构

求解一个体系的定态Schrödinger方程之前，首先要明确它的Hamilton算符。对于n电子原子来讲，考虑旋轨耦合后的Hamilton算符可以表示为：

2.1组态：Hˆ=Hˆ0精度下的多电子原子

这个行列式展开后共有n！项，每一项都是n个自旋轨道的连乘，相应于n个电子在n个自旋轨道中的一种排布方式，而n！项则囊括了n个电子在n个自旋轨道中的所有排布；当2个电子交换时，行列式两列互换，Ψ0相差一个负号，满足波函数对电子交换必须是反对称的要求。而原子的能量E0则为n个电子的能量之和：

所以，当给定每个电子的ni和li（给定电子所处的亚层）时，原子的能量就可以确定下来。而给定每个电子所处的亚层即给定一个电子组态，因此，在中心力场近似下原子的能级由电子组态来决定。

以C原子为例，其基态的电子组态为1s22s22p2，即6个电子中有2个在1s亚层，状态为1sα、1sβ；2个在2s亚层，状态为2sα、2sβ；还有2个在2p亚层，状态可以是2p+1α、2p+1β、2p0α、2p0β、2p-1α、2p-1β6个中的任意2个，共有=15种可能。如果是2p+1α、2p+1β，那么C原子的状态则用以下Slater行列式表示：类似的Slater行列式共有15个，均表示C原子可能的状态。而与组态1s22s22p2相应的能量E0=2E1s+ 2E2s+2E2p，其能级简并度为15。

考虑非球形作用Hˆʹ后，原子中电子所受的力不再是有心力，每个电子的能量、轨道角动量和自旋角动量都不再是常量；但原子的总能量、总轨道角动量及其z分量、总自旋角动量及其z分量都还是守恒的。能够同时确定的物理量有原子的总能量、总轨道角动量的平方及其z分量、总自旋角动量的平方及其z分量。因此，可以用ˆ这5个算符共同的本征函数来描述原子的状态，状态函数由L、ML、S和MS4个量子数决定；能量则由L和S 2个量子数决定，能级用光谱项2S+1L标记。因为在某个L和某个S下（确定能量），ML共有2L+1个可取的数值，MS共有2S+1个可取的数值，那么4个量子数L、ML、S和MS的组合（确定状态）共有（2L+1）（2S+1）个可取的数值，所以，光谱项2S+1L的能级简并度为（2L+1）（2S+1），即光谱项2S+1L包含的线性独立的波函数为（2L+1）（2S+1）个，这些波函数可以由组态波函数的线性组合得到。

例如，在考虑了非球形作用后，C的1s22s22p2组态在中心力场近似下15重简并的能级E0要分裂为3P（简并度为9）、1D（简并度为5）和1S（非简并）3个光谱项，这就是C原子在Hˆ0+Hˆʹ精度下的3个可取的能级。根据Hund规则，其中S取最大值的3P谱项的能量最低，因为自旋平行的电子必须占据不同的空间轨道，趋于互相回避的电子之间的排斥作用较小；需要注意不要用Hund规则给其余光谱项的能级顺序排队。进一步的研究表明1D谱项的能量比1S的要低，这可以从L越大电子越偏离球对称分布，电子间的排斥作用越小来解释。

对于C原子来讲，考虑旋轨耦合后，其3P光谱项要分裂为3P2（简并度为5）、3P1（简并度为3）和3P0（非简并）3个光谱支项，根据Hund规则，亚层半充满前J越小的能量越低，因此C的光谱基项（基支项）为3P0；光谱项1D（S=0时没有旋轨耦合）直接对应于光谱支项1D2（简并度为5）；光谱项1S直接对应于光谱支项1S0（非简并），这5个光谱支项就是C原子在精度下的5个可取的能级。而由MJ的不同所产生的光谱支项的简并度可以通过外加磁场来消除（Zeeman效应），此时，C原子的能量有15个可以取的数值，每个能级只对应于一个状态。

图1　C原子的能级示意图

关于光谱项和光谱支项更为详细的内容，因为超出了本科大纲的要求，我们不作过多的讨论，有兴趣的学生可参考徐光宪先生的《量子化学》中册［2］。最后作为本文的小结，在图1给出了不同精度的Hamilton算符下C原子的能级示意图。

［1］Levine，I.N.Quantum Chemistry.6th ed.北京:世界图书出版公司，2011.

［2］徐光宪，黎乐民，王德民.量子化学——基本原理和从头计算法.第2版.北京:科学出版社，2009.

The Structure of Many-Electron Atoms

CHEN LanSUN Hong-Wei*LAI Cheng-Ming
（College of Chemistry，Nankai University，Tianjin 300071，P.R.China）

In this paper，we represent the structures of many-electron atoms with different approximate Hamiltonian operators.Some related contents such as Slater determinants，electronic configurations，terms and levels are also discussed.

Structural chemistry；Many-electron atoms

O641；G64

10.3866/PKU.DXHX201509011

，Email:sunhw@nankai.edu.cn

国家基础科学人才培养基金（J1103306）