赝势

能带结构理论概念

赝势(pseudopotential),或有效势(effective potential),是指在对能带结构进行数值计算时所引入的一个虚拟的。引入赝势有助于实现一个复杂的系统的近似计算。事实上,赝势近似法是正交平面波方法(Orthogonalized Plane Wave method,OPW method)的延伸,其应用范围包括原子物理学中子散射英语Neutron scattering。“赝势”这个概念是由汉斯·赫尔曼英语Hans Hellmann于1934年首先发表的。[1]

库伦势以及库伦势下的波函数(蓝色)与赝势以及赝势下的波函数(红色)的比较。赝势 以及通过赝势计算出的波函数 与两者所对应的真实值在 大于某个临界半径 之上均是相符的。

概述

在赝势近似中,将原子的核电子(即價電子)以及原子核共同产生的一个复杂的势置换成一个“有效势”(赝势)之后,薛定谔方程中的库仑势能项会变成一个有利于进行下一步计算的有效势能项。所构造的赝势通过替代原子中所有电子共同产生的势,简化了原子中心部分的,从而可用包含较少节点的赝波函数来描述价电子。较少的节点意味着可以用较少的傅里叶级数项写出波函数,这也使平面波基组的计算变得实用。通常的计算中只考虑那些有化学反应活性的价电子;核电子则被看作和原子核“冻结”在一起,形成了一个刚性的不可极化的“粒子核”。根据所在的化学环境,自洽地更新赝势是一种修正上述“冻结的核心”的方法;但此做法较少见。

第一性原理的赝势是通过原子参照态(atomic reference state)推导出来的。这要求赝电子价本征态和全电子价本征态(pseudo- and all-electron valence eigenstates)在某个临界半径   之外有相同的能量和振幅。

临界半径较大的赝势被称作“软”赝势,具有更快的收敛速度,同时也更难模拟出现实系统的特征。

早期的赝势基于对原子光谱的拟合,并没有取得较大的成功。赝势在如今能获得广泛应用,很大一部分应归功于沃尔特·哈里森(Walter Harrison)在1958年对近自由电子费米面,以及詹姆斯·C·菲利普斯英语James Charles Phillips于同年对的共价能隙的成功拟合。后来,菲利普斯及其同事将此工作推广到其他的半导体中,并称其为“半经验赝势”(semiempirical pseudopotential)。[2]

范数守恒赝势和超软赝势

在现代的平面波电子结构数值计算英语List of quantum chemistry and solid-state physics software中,范数守恒(Norm-conserving)和超软(Ultrasoft)赝势是两种最常用的赝势。这两种赝势使基组可用较低的截断频率(即傅里叶展开项中的最高频率)来描述电子的波函数,从而在有限的计算资源下达到一定的数值收敛。这些方法的一个变种是线性缀加平面波方法英语Muffin-tin approximation(Linear Augmented Plane Wave,LAPW),即在原子核周围加上一些原子函数作为基组。

范数守恒赝势

范数守恒(Norm-conserving)赝势是由 Hamann,Schlüter 和 Chiang(HSC)于1979年首先提出的。[3] 最初的HSC范数守恒赝势的形式如下:

 

其中   将某一单粒子波函数,如科恩-沈吕九轨道,映射至由   标记的角动量。  是作用在被映射部分的赝势。不同的角动量态会受到不同的赝势作用,也就是说HSC范数守恒赝势是非局域性的;这一点与作用在整个单粒子波函数上的局域性赝势是不同的。

构造的范数守恒赝势需满足以下两个条件:

1. 临界半径   内,每一伪波函数的范数需与其所对应的全电子波函数相同,即[4]

 
其中    分别表示原子   上赝势的全电子参照态和伪参照态。

2. 全电子波函数和伪波函数在临界半径   外需要完全一致。

超软赝势

超软(Ultrasoft)赝势为了进一步缩小必须的基组集合,松弛(relax)了范数守恒赝势中的限制条件,引入了一个广义的本征值问题[5]。若范数间的差别非零,则可以定义:

 ,

因此赝哈密顿量的归一化本征态满足推广后的方程:

 ,

其中,算符   被定义为:

 ,

  是在截断频率内通过赝参照态(pseudo reference state)形成对偶空间的投影(projector),在截断频率外取的值为零:

 .

投影缀加平面波方法英语projector augmented wave method(PAW)与此相关[6]

费米赝势

费米赝势是恩里科·费米为了描述自由中子原子核的散射而引入的[7]。散射被假设为s波英语partial wave analysis散射,因此具有球对称性,是一个与半径   相关的函数:

 ,

其中, 约化普朗克常数 质量 狄拉克δ函数  是束缚相干(bound coherent)中子散射长度[8]。对此δ函数进行傅里叶变换将得到为常数的中子构型因素英语Atomic form factor

参见

相关文章

外部链接

参考文献

  1. ^ Schwerdtfeger, P., The Pseudopotential Approximation in Electronic Structure Theory, ChemPhysChem, August 2011, doi:10.1002/cphc.201100387 
  2. ^ M. L. Cohen, J. R. Chelikowsky, "Electronic Structure and Optical Spectra of Semiconductors", (Springer Verlag, Berlin 1988)
  3. ^ Hamann, D. R.; Schlüter, M.; Chiang, C. Norm-Conserving Pseudopotentials. Physical Review Letters. 1979-11-12, 43 (20): 1494–1497. doi:10.1103/PhysRevLett.43.1494. 
  4. ^ Bachelet, G. B.; Hamann, D. R.; Schlüter, M., Pseudopotentials that work: From H to Pu, Physical Review B 26 (8) (American Physical Society), October 1982, 26 (8): 4199–4228, Bibcode:1982PhRvB..26.4199B, doi:10.1103/PhysRevB.26.4199 
  5. ^ Vanderbilt, David, Soft self-consistent pseudopotentials in a generalized eigenvalue formalism, Physical Review B 41 (11) (American Physical Society), April 1990, 41 (11): 7892–7895, Bibcode:1990PhRvB..41.7892V, doi:10.1103/PhysRevB.41.7892 
  6. ^ Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. 1999. Bibcode:1999PhRvB..59.1758K. doi:10.1103/PhysRevB.59.1758. 
  7. ^ E. Fermi, Motion of neutrons in hydrogenous substances, Ricerca Scientifica, July 1936, 7: 13–52 
  8. ^ Squires, Introduction to the Theory of Thermal Neutron Scattering, Dover Publications (1996) ISBN 0-486-69447-X