多电子原子的结构
The Structure of Many-Electron Atoms
通讯作者:
基金资助: |
|
Fund supported: |
国家基础科学人才培养基金. |
论述了不同精度Hamilton算符下多电子原子的结构,并对Slater行列式、电子组态、原子光谱项和光谱支项等相关内容进行了讨论。
关键词:
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.
Keywords:
本文引用格式
陈兰, 孙宏伟, 赖城明.
CHEN Lan, SUN Hong-Wei, LAI Cheng-Ming.
结构化学中无论研究原子还是分子总是从求解定态Schrödinger方程Ĥψ = Eψ入手,并以方程的解E和y表示体系的能量与状态。为什么用y可以表示体系的状态?大多数学生在学习了结构化学之后仍回答不了这个问题;而对多电子原子结构中原子光谱项和光谱支项的掌握则只限于会推求而并不了解它们的意义。本文从为什么求解Schrödinger方程谈起,在采用不同精度的Hamilton算符下,对Slater行列式、电子组态、原子光谱项和光谱支项等相关内容进行了讨论,希望能对学生系统地认识多电子原子的结构有所帮助。
1 定态Schrödinger方程的解ψ( ${\vec{r}}$ )与原子的状态
量子力学中假设微观体系的状态用波函数Ψ(
式(1)中Ĥ为Hamilton算符,它包括动能算符和势能函数两部分。动能算符不包含时间,如果势能函数也与时间无关,那么可以对Ψ(
式(2)左右两边同时除以ψ(
式(3)左边仅与时间有关,而右边仅与坐标有关,只有当左右两边都等于同一常数时,等式才能成立。以E表示这个常数,由式(3)左边等于E有:
对式(4)两边积分,可解得f (t),见式(5),其中C为积分常数:
而由式(1)右边等于E有:
式(6)是Hamilton算符的本征方程,称之为定态Schrödinger方程(或与时间无关的Schrödinger方程)。由该方程的解ψ(
式(5)中的常数C已并入ψ(
所以,状态Ψ(
综上,对势能函数不显含时间的体系,通过求解定态Schrödinger方程,得到E与ψ(
2 采用不同精度的Hamilton算符下多电子原子的结构
求解一个体系的定态Schrödinger方程之前,首先要明确它的Hamilton算符。对于n电子原子来讲,考虑旋轨耦合后的Hamilton算符可以表示为:
其中Ĥ0为中心力场近似下的能量算符;Ĥʹ为非球形作用能算符,如果中心力场选得好,Ĥʹ可以很小;Ĥso为旋轨耦合能算符,对于不太重的原子也比较小。因此,Ĥʹ和Ĥso均可以作为微扰项来处理,哪个更重要就先考虑哪个。一般当原子的核电荷数Z < 40时,Ĥʹ > Ĥso,可以先考虑Ĥʹ后考虑Ĥso,并称之为L-S耦合。以下即采用L-S耦合讨论不同精度的Hamilton算符下多电子原子的结构。
2.1 组态:Ĥ = Ĥ0精度下的多电子原子
中心力场近似的基础是电子独立运动近似,即原子中的每个电子都可以看成是在与原子核及其他电子的作用相等效的势场中作独立的运动。在电子独立运动近似下,原子总的Hamilton算符就等于n个电子的Hamilton算符之和,即Ĥ0 =
式(11)中的Vi(ri)为电子i所处的等效中心势场。因为Vi(ri)只和ri有关,与θi和ϕi都没有关系,所以用分离变量法求解方程(11)时,Vi(ri)只进入R方程,对
这个行列式展开后共有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个,共有C26 = 15种可能。如果是2p+1α、2p+1β,那么C原子的状态则用以下Slater行列式表示:
类似的Slater行列式共有15个,均表示C原子可能的状态。而与组态1s22s22p2相应的能量E0 = 2E1s +2E2s + 2E2p,其能级简并度为15。
一般地,给定电子组态(n 1l1)x1(n 2l2)x2(n 3l3)x3⋯,能量E0 = x1En1l1+ x2En2l2+ x3En3l3+ ⋯,而能级简并度为C2(2l1+ 1)x1C2(2l2+ 2)x2C2(2l3+ 3)x3⋯。
2.2 光谱项:Ĥ = Ĥ0 + Ĥʹ精度下的多电子原子
考虑非球形作用Ĥʹ后,原子中电子所受的力不再是有心力,每个电子的能量、轨道角动量和自旋角动量都不再是常量;但原子的总能量、总轨道角动量及其z分量、总自旋角动量及其z分量都还是守恒的。能够同时确定的物理量有原子的总能量、总轨道角动量的平方及其z分量、总自旋角动量的平方及其z分量。因此,可以用Ĥ、L̂2、L̂z、Ŝ2、Ŝz这5个算符共同的本征函数来描述原子的状态,状态函数由L、ML、S和MS 4个量子数决定;能量则由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原子在Ĥ0 + Ĥʹ精度下的3个可取的能级。根据Hund规则,其中S取最大值的3P谱项的能量最低,因为自旋平行的电子必须占据不同的空间轨道,趋于互相回避的电子之间的排斥作用较小;需要注意不要用Hund规则给其余光谱项的能级顺序排队。进一步的研究表明1D谱项的能量比1S的要低,这可以从L越大电子越偏离球对称分布,电子间的排斥作用越小来解释。
2.3 光谱支项:Ĥ = Ĥ0 + Ĥʹ + Ĥso精度下的多电子原子
非球形作用Ĥʹ和旋轨耦合Ĥso同时考虑后,原子的总轨道角动量和总自旋角动量都不再是常量;只有原子的总能量、总角动量及其z分量是守恒的。能够同时确定的物理量有原子的总能量、总角动量的平方及其z分量。因此,可以用Ĥ、Ĵ2、Ĵz这3个算符共同的本征函数来描述原子的状态,状态函数由L、S、J 和MJ 4个量子数决定;能量则由L、S 和J 3个量子数决定,能级用光谱支项2S+1LJ标记。因为在某个J下MJ共有2J + 1个可取的数值,所以光谱支项2S+1LJ的能级简并度为2J + 1。
对于C原子来讲,考虑旋轨耦合后,其3P光谱项要分裂为3P2(简并度为5)、3P1(简并度为3)和3P0(非简并)3个光谱支项,根据Hund规则,亚层半充满前J越小的能量越低,因此C的光谱基项(基支项)为3P0;光谱项1D (S = 0时没有旋轨耦合)直接对应于光谱支项1D2(简并度为5);光谱项1S直接对应于光谱支项1S0(非简并),这5个光谱支项就是C原子在Ĥ0 + Ĥʹ + Ĥso精度下的5个可取的能级。而由MJ的不同所产生的光谱支项的简并度可以通过外加磁场来消除(Zeeman效应),此时,C原子的能量有15个可以取的数值,每个能级只对应于一个状态。
图1
参考文献
/
〈 |
|
〉 |
