更新时间:2022-08-25 17:36
赝势(pseudopotential),或有效势(effective potential),是指在对能带结构进行数值计算时所引入的一个虚拟的势。引入有效势有助于实现一个复杂的系统的近似计算。事实上,有效势近似法是正交平面波方法(Orthogonalized Plane Wave method,OPW method)的延伸,其应用范围包括原子物理学和中子散射。“有效势”这个概念是由汉斯·赫尔曼于1934年首先发表的。
在有效势近似中,将原子的核电子(即非价电子)以及原子核共同产生的一个复杂的势置换成一个“有效势”(有效势)之后,薛定谔方程中的库仑势能项会变成一个有利于进行下一步计算的有效势能项。所构造的有效势通过替代原子中所有电子共同产生的势,简化了原子中心部分的态,从而可用包含较少节点的赝波函数来描述价电子。较少的节点意味着可以用较少的傅里叶级数项写出波函数,这也使平面波基组的计算变得实用。通常的计算中只考虑那些有化学反应活性的价电子;核电子则被看作和原子核“冻结”在一起,形成了一个刚性的不可极化的“粒子核”。根据所在的化学环境,自洽地更新有效势是一种修正上述“冻结的核心”的方法;但此做法较少见。
临界半径较大的有效势被称作“软”有效势,具有更快的收敛速度,同时也更难模拟出现实系统的特征。
早期的有效势基于对原子光谱的拟合,并没有取得较大的成功。有效势在如今能获得广泛应用,很大一部分应归功于沃尔特·哈里森(Walter Harrison)在1958年对铝的近自由电子的费米面,以及詹姆斯·C·菲利普斯于同年对硅和锗的共价能隙的成功拟合。后来,菲利普斯及其同事将此工作推广到其他的半导体中,并称其为“半经验有效势”(semiempirical pseudopotential)。
在现代的平面波电子结构数值计算中,范数守恒(Norm-conserving)和超软(Ultrasoft)有效势是两种最常用的有效势。这两种有效势使基组可用较低的截断频率(即傅里叶展开项中的最高频率)来描述电子的波函数,从而在有限的计算资源下达到一定的数值收敛。这些方法的一个变种是线性缀加平面波方法(Linear Augmented Plane Wave,LAPW),即在原子核周围加上一些原子函数作为基组。
范数守恒(Norm-conserving)有效势是由 Hamann,Schlüter 和 Chiang(HSC)于1979年首先提出的。最初的HSC范数守恒有效势的形式如下:
其中将某一单粒子波函数,如科恩-沈吕九轨道,映射至由标记的角动量。是作用在被映射部分的有效势。不同的角动量态会受到不同的有效势作用,也就是说HSC范数守恒有效势是非局域性的;这一点与作用在整个单粒子波函数上的局域性有效势是不同的。
构造的范数守恒有效势需满足以下两个条件:
(1)临界半径 rc内,每一伪波函数的范数需与其所对应的全电子波函数相同,即
其中 和 分别表示原子上有效势的全电子参照态和伪参照态。
(2)全电子波函数和伪波函数在临界半径外需要完全一致。
超软(Ultrasoft)有效势为了进一步缩小必须的基组集合,松弛(relax)了范数守恒有效势中的限制条件,引入了一个广义的本征值问题。若范数间的差别非零,则可以定义:
因此赝哈密顿量的归一化本征态满足推广后的方程:
其中,算符 被定义为:
是在截断频率内通过赝参照态(pseudo reference state)形成对偶空间的投影(projector),在截断频率外取的值为零:
投影缀加平面波方法(PAW)与此相关。
费米有效势是恩里科·费米为了描述自由中子受原子核的散射而引入的。散射被假设为s波散射,因此具有球对称性,是一个与半径 r相关的函数:
其中,为约化普朗克常数, m为质量,是狄拉克δ函数,b是束缚相干(bound coherent)中子散射长度。对此δ函数进行傅里叶变换将得到为常数的中子构型因素。