赝势

✍ dations ◷ 2025-07-11 19:26:59 #计算物理学,原子物理学

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

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

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

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

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

在现代的平面波电子结构数值计算(英语: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年首先提出的。 最初的HSC范数守恒赝势的形式如下:

V ^ ps ( r ) = l m | Y l m V l m ( r ) Y l m | {\displaystyle {\hat {V}}_{\textit {ps}}(r)=\sum _{l}\sum _{m}|Y_{lm}\rangle V_{lm}(r)\langle Y_{lm}|}

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

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

1. 临界半径 r c {\displaystyle r_{c}} 内,每一伪波函数的范数需与其所对应的全电子波函数相同,即

2. 全电子波函数和伪波函数在临界半径 r c {\displaystyle r_{c}} 外需要完全一致。

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

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

其中,算符 S ^ {\displaystyle {\hat {S}}} 被定义为:

p R , i {\displaystyle p_{\mathbf {R} ,i}} 是在截断频率内通过赝参照态(pseudo reference state)形成对偶空间的投影(projector),在截断频率外取的值为零:

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

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

V ( r ) = 4 π 2 m b δ ( r ) {\displaystyle V(r)={\frac {4\pi \hbar ^{2}}{m}}b\,\delta (r)} ,

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

相关

  • VIAbr /16固体、 液体、 气体氧族元素是指元素周期表上第16族(ⅥA族)的元素,位于氮族元素和卤素之间。氧族元素包含氧(O)、硫(S)、硒(Se)、碲(Te)、钋(Po)、钅立(Lv),其中氧、硫、硒为非金属,碲为类金
  • 测绘学形貌测绘学研究测定和推算地面几何位置、地球形状及地球重力场,据此测量地球表面自然物体和人工设施的几何分布,编制各种比例尺地图的理论和技术的学科。测绘学的研究对象是地球的
  • 南丹麦大区南丹麦大区 (Region Syddanmark)是丹麦五大区之一,2007年1月1日由原来的菲英郡、里伯郡、南日德兰郡和瓦埃勒郡的十个市合并而成。下分22市。面积12,191平方公里(半岛部分8705k
  • 平方百米公顷(hectare)为面积的公制单位(国际单位),为国际单位制并用单位之一。国际上,缩写为ha;在中国大陆的法定缩写为hm²。一块面积一公顷的土地,大约与标准足球场近似,此外,公顷也称为(HA)
  • DDD2,2-双(对氯苯基)-1,1-二氯乙烷(DDD)是DDT主要代谢产物之一。1943年首先由Geigy公司合成。后来由Rohm Haas公司开发为商品。可由二氯乙醛与氯苯缩合而得。工业品DDT根据来源一
  • 宗室载萼宗室载萼(1825年-1900年),载萼,字兰舫,号瓣香,爱新觉罗氏,清朝政治人物、进士出身。道光二十九年,乡试中举。后授宗人府理事官。光绪三年,任杀虎口监督。光绪二十年,任内阁侍读学士,后改
  • 反应机制化学中,反应机理用来描述某一化学变化所经由的全部基元反应。虽然整个化学变化所发生的物质转变可能很明显,但为了探明这一过程的反应机理,常常需要实验来验证。机理详细描述了
  • 托莱多 (西班牙)托莱多(西班牙语:Toledo,西班牙语:.mw-parser-output .IPA{font-family:"Charis SIL","Doulos SIL","Linux Libertine","Segoe UI","Lucida Sans Unicode","Code2000","Gentium"
  • 中井猛之进中井猛之进(1882年11月9日-1952年),日本植物学家,东京大学教授、小石川植物园园长、国立科学博物馆馆长。作家中井英夫是其儿子。中井猛之进命名的生物分类
  • 革囊袋鼬属革囊袋鼬属(红腹袋鼬),哺乳纲的一属,而与革囊袋鼬属(红腹袋鼬)同科的动物尚有侏袋鼬属(侏袋鼬)、宽足带鼩属(红耳宽足带鼩)、帚尾袋鼩属(红帚尾袋鼩)、斑袋鼬属(斑袋鼬)等之数种哺乳动物。