第一性原理计算及相关理论方法
- 格式:doc
- 大小:282.34 KB
- 文档页数:8
材料科学中第一性原理计算方法研究近年来,材料科学领域的研究取得了许多重大突破,其中第一性原理计算方法成为材料设计和研究的重要工具之一。
这种方法通过基本的物理原理和数学方程来研究材料的性质和行为,为材料设计和性能优化提供了新的途径。
第一性原理计算方法是基于量子力学的一种计算方法,从第一性原理出发,通过求解薛定谔方程以及其他相关方程来研究材料的性质。
它不依赖于任何经验参数或假设,能够提供对材料的精确描述和准确预测。
第一性原理计算方法的核心是密度泛函理论(Density Functional Theory,简称DFT),它将体系的物理性质与体系中电子的密度联系起来。
根据Kohn-Sham方程,DFT通过对电子的运动方程进行求解,得到体系的基态电子密度。
通过计算得到的电子密度,可以进一步计算出材料的能带结构、电子态密度、态密度、声子谱、磁性及其它性质。
与传统的实验方法相比,第一性原理计算方法具有独特的优势。
首先,它能够提供物理性质的原子尺度描述,可以捕捉到材料内部微观原子结构的信息。
其次,该方法能够计算和预测材料的多种性质,如电子能带结构、晶格常数、弹性性能、热力学性质等,为材料设计和开发提供了重要参考。
此外,第一性原理计算方法可以帮助解释材料性能背后的基本物理机制,揭示材料特性的微观本质。
近年来,随着计算机性能的不断提升和计算方法的进步,第一性原理计算方法在材料科学中的应用得到了广泛拓展。
例如,它在材料的合成、器件的设计和材料的特性优化等方面发挥了重要作用。
通过预测和优化材料的能带结构和电子态密度,可以筛选出具有优异性能的新材料,为新能源、环境友好材料、传感器和光电器件的研发提供重要支持。
此外,第一性原理计算方法还可以帮助优化材料的力学、热力学和电磁性能,提高材料的功能性能。
尽管第一性原理计算方法为材料科学提供了强大的工具和理论基础,但也面临一些挑战。
首先,该方法对计算所需的资源要求较高,需要大量计算时间和计算机内存。
第一性原理理论介绍第一性原理理论是一种基于量子力学的理论,用于解释材料和分子的性质和行为。
它是通过求解薛定谔方程来描述原子核和电子之间的相互作用,从而得出系统的总能量。
第一性原理理论被广泛应用于材料科学、物理化学和计算物理等领域,为设计新材料、预测化学反应和模拟材料性质提供了重要的工具。
第一性原理理论的核心是薛定谔方程,它描述了系统的波函数随时间的演化。
薛定谔方程包括了系统的势能和动能项,其中势能项描述了原子核和电子之间的相互作用,动能项则描述了电子的运动。
解薛定谔方程可以得到系统的波函数,进而可以计算系统的总能量。
在求解薛定谔方程时,第一性原理理论通常采用密度泛函理论(DFT)作为基础。
DFT是一种将电子系统的性质与电子密度之间建立关联的方法。
根据Kohn-Sham方程,系统的能量可以表示为电子密度的泛函形式。
为了将电子相互作用考虑在内,通常使用电子交换关联泛函来近似描述系统的能量。
第一性原理理论已经成为材料科学和计算物理的重要工具。
它可以用于预测材料的结构和稳定性,计算材料的力学性质和电子结构,模拟化学反应和催化过程,设计新的材料和催化剂等。
特别是在材料发现和设计中,第一性原理理论具有重要的意义,可以指导实验研究,加速材料研发过程。
总之,第一性原理理论是一种基于量子力学的理论,通过求解薛定谔方程来描述原子核和电子之间的相互作用。
它是预测和解释材料和分子性质的重要工具,广泛应用于材料科学、物理化学和计算物理等领域。
通过第一性原理计算,我们可以更好地理解和控制材料的性质,促进科学研究和技术创新的发展。
新能源材料研究中的第一性原理计算近年来,随着节能减排和环保意识的逐步加强,新能源的开发和利用已成为世界各国共同关注的焦点。
而为了更有效地提高新能源的利用效率和降低成本,科学家们开始转向新能源材料的研究和开发。
在这一过程中,第一性原理计算发挥着越来越重要的作用。
第一性原理计算是指基于量子力学理论和数学方法对材料的电子结构和性质进行计算和模拟。
这种计算方法的好处在于既能提供高精度的计算结果,又能对材料的微观结构和电子能带等性质进行深入分析,为新材料的设计和开发提供有力的支持。
在新能源材料研究中,第一性原理计算可以帮助科学家们确定材料的电子结构、晶格结构、热力学性质、光电特性等重要参数。
以太阳能电池材料为例,研究者可以通过第一性原理计算预测材料的光吸收性能、载流子输运特性和光电转换效率等重要指标,从而优化材料的能带结构和界面特性,提高太阳能电池的转化效率。
除了太阳能电池材料之外,第一性原理计算在其他新能源领域的研究中也发挥着重要作用。
比如,在固态氢储存材料的研究中,第一性原理计算可以用来预测材料的结晶形态、氢吸附能力和释放能力等关键性质,为研发更高效、更安全的氢储存材料提供支持。
在燃料电池材料的研究中,第一性原理计算可以预测氧化还原反应的能垒、电子传输特性和催化活性等参数,为提高燃料电池的效率和寿命提供重要帮助。
需要指出的是,尽管第一性原理计算具有高计算精度和深入分析的优点,但该方法也存在一些挑战和限制。
其中,计算复杂度是最主要的问题之一。
由于第一性原理计算需要对大量的原子和电子进行计算,因此计算量非常大,需要使用高性能计算机进行处理。
而由于计算复杂度高,一些材料的性质无法通过第一性原理计算来预测,需要通过实验来验证。
另一方面,第一性原理计算还需要与实验相结合,以验证计算结果的准确性和可靠性。
特别是在新能源材料研究中,第一性原理计算和实验之间的结合非常重要。
通过实验,科学家们可以验证计算结果,并不断优化计算模型,提高计算精度和可靠性。
材料力学性质的第一性原理计算研究第一性原理计算是现代材料科学研究中重要的工具之一。
通过基于量子力学的第一性原理计算方法,可以预测材料的各种性质。
在材料科学的发展中,探索并理解材料的力学性质一直是一个关键问题。
本文将以第一性原理计算为基础,深入研究材料力学性质的预测和分析。
首先,我们来了解一下第一性原理计算的基本原理。
第一性原理计算方法是基于薛定谔方程和密度泛函理论的计算方法,能够计算出材料的电子结构、能带结构、晶格常数等基本性质。
在这个基础上,可以进一步计算材料的弹性性质、断裂强度、塑性形变等力学性质。
通过计算材料的原子结构和晶体缺陷,可以预测力学性能的变化规律,为材料设计和优化提供理论指导。
接下来,我们将通过具体的案例来说明第一性原理计算在材料力学性质研究中的应用。
以金属材料为例,我们可以通过第一性原理计算方法研究材料的弹性性质。
首先,需要计算金属的晶体结构和原子排列。
然后,通过计算材料的局部应变和应力,可以得到金属的弹性常数。
这些弹性常数包括剪切模量、杨氏模量等,可以描述金属在外力作用下的变形特性。
另外,我们还可以通过第一性原理计算来研究材料的断裂性质。
断裂是材料损伤和失效的重要原因之一。
通过计算材料的断裂表面能量和断裂韧性等参数,可以预测材料的断裂强度和断裂模式。
这对于材料的设计和改良具有重要的指导意义。
例如,在薄膜材料领域,第一性原理计算可以用来研究薄膜的断裂行为,为薄膜的应用和制备提供理论依据。
除了弹性性质和断裂性质,第一性原理计算还可以用于研究材料的塑性形变和变形机制。
材料的塑性形变是材料经历应力后产生永久形变的过程。
通过计算材料中的晶体缺陷如位错和空位,可以模拟材料的塑性变形。
通过分析位错的运动和材料的能量变化,可以揭示材料塑性形变和变形机制的微观本质。
这对于提高材料的塑性变形能力和改善材料的塑性加工性能具有重要意义。
最后,我们可以看到第一性原理计算在研究材料力学性质中的潜力和应用前景。
第一性原理,英文First Principle,是一个计算物理或计算化学专业名词,广义的第一性原理计算指的是一切基于量子力学原理的计算。
我们知道物质由分子组成,分子由原子组成,原子由原子核和电子组成。
量子力学计算就是根据原子核和电子的相互作用原理去计算分子结构和分子能量(或离子),然后就能计算物质的各种性质。
从头算(ab initio)是狭义的第一性原理计算,它是指不使用经验参数,只用电子质量,光速,质子中子质量等少数实验数据去做量子计算。
但是这个计算很慢,所以就加入一些经验参数,可以大大加快计算速度,当然也会不可避免的牺牲计算结果精度。
量子力学第一性原理:仅需五个物理基本常数——电子质量、电子电量、普郎克常数、光速和玻耳兹曼常数,通过求薛定谔方程得到材料的电子结构,而不依赖于任何经验常数即可以预测微观体系的状态和性质,预测材料的组分、结构、性能之间的关系,进一步设计具有特定性能的新材料。
作为评价事物的依据,第一性原理和经验参数是两个极端。
第一性原理是某些硬性规定或推演得出的结论,而经验参数则是通过大量实例得出的规律性的数据,这些数据可以来自第一性原理(称为理论统计数据),也可以来自实验(称为实验统计数据)。
如果某些原理或数据来源于第一性原理,但推演过程中加入了一些假设(这些假设当然是很有说服力的),那么这些原理或数据就称为“半经验的”。
量子化学的第一性原理是指多电子体系的Schrödinger方程,但是光有这个方程是无法解决任何问题的,量子力学能够准确的解决的问题很少很少,绝大多数都是有各种各样的近似,为此计算量子力学提出一个称为“从头计算”的原理作为第一性原理,除了Schrödinger 方程外还允许使用下列参数和原理:(1) 物理常数,包括光速c、Planck常数h、电子电量e、电子质量me以及原子的各种同位素的质量,尽管这些常数也是通过实验获得的。
(在国际单位值中,光速是定义值,Planck 常数是测量值,在原子单位制中则相反。
性原理计算2023-11-08CATALOGUE目录•引言•声声耦合理论•电声耦合理论•声声耦合的第一性原理计算•电声耦合的第一性原理计算•声声耦合与电声耦合的模拟方法与软件01引言声声耦合和电声耦合在物理领域具有重要应用价值理解声声耦合和电声耦合的机理对于新材料设计、器件优化等方面具有指导意义研究背景与意义研究内容研究声声耦合和电声耦合的物理机理,探索影响耦合的因素,并预测新材料和器件的性能。
研究方法采用第一性原理计算方法,结合量子力学和经典力学理论,对声声耦合和电声耦合进行模拟和计算。
研究内容与方法02声声耦合理论声波传播原理声速与媒质特性声波的传播速度取决于媒质的密度和弹性模量。
边界条件与反射在媒质的边界上,声波的传播方向会发生改变,这导致反射和折射现象。
波动现象的基本方程描述声波传播的波动方程是偏微分方程,它描述了声源引起的媒质密度和压力的变化。
03散射和吸收的微观机制了解微观机制有助于理解散射和吸收现象的本质。
声波的散射与吸收01散射现象当声波遇到比声波波长大的物体时,会发生散射现象,声波的传播方向会发生改变。
02吸收现象当声波通过媒质传播时,由于媒质的内摩擦和热传导,声波的能量会逐渐被吸收。
当两个或多个声波相遇时,它们会相互叠加,产生干涉现象。
干涉现象衍射现象干涉和衍射的应用当声波遇到障碍物或狭缝时,会发生衍射现象,声波会绕过障碍物或穿过狭缝。
了解干涉和衍射现象在声学中的应用有助于更好地理解声波的传播和控制。
03声波的干涉与衍射020103电声耦合理论描述电荷之间的相互作用力。
库伦定律研究静止电荷产生的电场和电势。
静电学描述电流的流动和电阻对电流的阻碍。
电流和电阻电学基础电极化与电导电极化描述电介质在电场作用下的极化现象。
电导描述电介质的电导率,即电场作用下电介质的导电性能。
电-声相互作用电致伸缩描述电场作用下电介质的形变。
声波传播描述声波在电介质的传播特性。
04声声耦合的第一性原理计算密度泛函理论是研究固体材料的电子结构和物理性质的理论框架,可以准确地描述材料的电子状态、电荷密度分布以及能量等。
第一性原理计算方法在材料科学中的应用1.引言第一性原理计算方法(First Principles Calculation)是近年来发展的新型计算方法,用于准确计算分子和固体物质的能量、结构和物理性质。
它的优势在于不依赖于实验数据,可以直接从基本原理推导出体系的特性。
在材料科学领域,第一性原理计算方法已经成为研究材料的重要工具,可以为合成新材料和设计功能材料提供理论依据,并指导实验研究。
2.第一性原理计算方法的基本原理第一性原理计算方法的基本原理是量子力学中的密度泛函理论,它的基本假设是所有粒子的运动都可以描述为波函数的运动。
根据波函数理论,一个由N个电子和原子核组成的体系的波函数可以用N个单电子波函数表示。
通过求解薛定谔方程,可以确定体系的基态能量和电子的密度,从而得到体系的性质。
3.第一性原理计算方法在材料科学中的应用(1)材料合成第一性原理计算方法可以模拟材料的结构、动力学和化学反应,为材料合成提供理论指导。
例如,使用第一性原理计算方法可以预测材料的稳定性、生长机制和晶体缺陷,从而为材料的设计和制备提供指导。
(2)材料性能第一性原理计算方法可以计算材料的电子结构、热力学性质、光电性质和磁学性质等,从而为材料的性能研究提供理论基础。
例如,通过计算材料的电子结构,可以预测材料的导电性、热导率和热电性能等,为相关应用提供指导。
(3)材料改性第一性原理计算方法可以模拟材料的界面和表面结构,研究材料的改性效果。
例如,可以通过计算材料与其他材料的界面能量来评估材料的附着性和界面稳定性,从而指导材料的改性设计。
(4)功能材料设计借助第一性原理计算方法,可以针对具体的应用需求,设计出具有特定功能的材料。
例如,通过计算材料的光电性质、催化活性和磁学性质等,可以指导材料的功能设计,为实现特定的应用提供理论指导。
4.发展趋势随着材料科学和计算科学的发展,第一性原理计算方法的应用前景越来越广阔。
未来,第一性原理计算方法将会与机器学习和高通量计算等技术结合,为材料科学的研究提供更多的可能性。
化学反应动力学的第一性原理计算方法化学反应是物质变化的一种形式,通常是指原子或分子之间的化学键被打破或形成,从而形成新的化合物。
化学反应动力学研究的是化学反应速率的研究,也就是反应物转变为产物的速率。
动力学的研究对于理解化学反应机理和制定化学反应工艺有着重要的意义。
在现代化学研究中,化学反应动力学的第一性原理计算方法已经成为重要的工具。
化学反应动力学的第一性原理计算方法指的是运用量子力学原理和分子动力学模拟技术对化学反应动力学过程进行精确的计算和模拟。
这种方法无需依靠实验数据,而是直接从微观层面分析分子之间的相互作用。
通过对分子结构和动力学过程的分析,可以计算得到反应动力学的速率常数、反应机理、反应能垒等其它重要参数,从而能够深入理解化学反应的本质。
化学反应动力学的第一性原理计算方法主要应用于分子动力学模拟和量子化学计算两个方面。
其中,分子动力学模拟方法主要是基于原子力场,通过数值积分求解牛顿方程,模拟反应过程。
它可以计算物质的结构、能量以及动力学过程。
量子化学计算方法则是基于量子力学理论,通过求解薛定谔方程,计算分子间的相互作用和反应机理。
这种方法可以计算各种化学反应的能垒、活化能、反应速率以及反应机理。
对于化学反应动力学的第一性原理计算方法,其中一个比较重要的问题就是如何评估理论计算的准确性。
实际上,在计算化学的过程中,化学反应动力学的第一性原理计算方法也不能完全避免计算误差。
因此,如何评估计算误差以及如何优化理论计算模型是这个领域研究者一直在关注的问题。
面对以上问题,化学反应动力学的第一性原理计算方法的研究者们借鉴了机器学习的思想,开发出了一种基于数据库和机器学习的化学反应动力学数据驱动模型。
该模型基于已有的关于反应动力学的实验数据和理论计算数据,通过机器学习方法对反应动力学模型进行训练、验证和优化。
这种模型可以有效地降低计算误差,提高计算准确性,并能够提高计算速度。
总之,化学反应动力学的第一性原理计算方法是指在量子力学和分子动力学的基础上,通过计算和模拟分子间的相互作用和反应过程来研究化学反应动力学的的方法。
晶胞计算技巧晶胞计算是固态物理和材料科学中的一项重要技术,它可以用来研究晶体结构、物性和材料性能。
在进行晶胞计算时,有一些技巧和方法可以提高计算效率和准确度。
本文将介绍一些晶胞计算的相关参考内容。
1. 第一性原理计算方法第一性原理计算方法是基于量子力学原理和密度泛函理论的计算方法。
它可以用来计算晶体的电子结构和能带结构。
在进行第一性原理计算时,需要选择适当的计算方法和基组,其中一些常用的方法包括密度泛函理论方法、赝势方法和平面波基组方法等。
2. 晶体结构优化晶体结构优化是指对晶格参数和原子坐标进行优化,以得到最稳定的晶体结构。
在晶胞计算中,晶体结构优化是一个重要的步骤,可以通过调整晶格参数和原子位置来获得准确的晶体结构。
常用的结构优化方法包括经典力学方法、分子动力学方法和共轭梯度方法等。
3. 原子间相互作用模型原子间相互作用模型是描述晶体中原子间相互作用力的模型。
晶体中的原子间相互作用力可以通过势能函数来描述,其中一些常用的势能函数包括Lennard-Jones势能函数、Morse势能函数和Stillinger-Weber势能函数等。
选择合适的原子间相互作用模型可以提高晶胞计算的准确度。
4. 基态能量计算基态能量是指晶体在零温下的稳定能量。
在晶胞计算中,基态能量的计算是一个重要的指标,可以用来评估晶体的稳定性。
基态能量的计算可以通过第一性原理计算方法来进行,其中一些常用的方法包括密度泛函理论方法和Hartree-Fock方法等。
5. 增加k点数目对于周期性较强的晶体体系,使用较少数量的k点可能会导致计算结果的欠收敛。
增加k点数目可以提高计算精度,得到更准确的结果。
在进行晶胞计算时,可以通过网格搜索方法来确定合适的k点数目,以获得最优的计算结果。
6. 减少计算误差在进行晶胞计算时,会存在一些误差,如数值误差、模拟误差和近似误差等。
减少计算误差可以提高计算结果的精度和准确度。
一些常用的减少误差的方法包括使用更高的基组、增加计算网格密度和使用更精确的数值算法等。
金属材料晶体缺陷的第一性原理计算与研究金属材料是现代工业中最为重要的材料之一,广泛应用于各个领域中。
其中,金属材料的力学性能是需要被考虑的一个极为关键的因素。
在金属材料中,晶体缺陷是导致力学性能变化的重要因素。
因此,晶体缺陷的研究和计算已经成为了金属材料研究领域的一个重要分支。
晶体缺陷可以分为点缺陷、线缺陷和面缺陷三种类型。
点缺陷通常指晶体的点位上缺失了一个原子或者晶格中插入了一个附加原子。
而线缺陷描述了一系列连续的点缺陷,如晶体中的位错线或嵌错。
面缺陷则描述了晶体表面或内部晶界(晶格方向或晶格步长不连续处)的变形。
这些缺陷会导致金属材料的性质明显变化,如塑性、疲劳、蠕变、断裂等力学性能。
第一性原理计算方法指的是利用从一组基本物理原理出发的计算方法进行计算,不需要任何实验数据和经验参数。
第一性原理方法被广泛应用于各种领域,特别是在材料科学和物理化学领域中,其应用范围涵盖了材料结构、动力学、热力学、电学和磁学性质等。
首先,第一性原理方法计算的结果通常具有很高的可靠性,与实验结果吻合度很高,比起修正参数和拟合本质上更精确。
其次,该方法的可应用范围非常广泛,从凝聚态物理到表面和界面等微观尺度都有涉及。
第一性原理方法可以用于模拟材料中的各种晶体缺陷,例如金属材料中常见的位错、孪生、晶界、空隙和间隙原子等。
通过计算这些缺陷的性质,可以更全面地了解金属材料的力学性质。
对于金属材料的不同性质,第一性原理方法也可通过调整缺陷的类型、密度和形状实现定量预测。
举例而言,位错是晶体内最简单和基本的线缺陷,是晶体中原子排列存在缺陷的部分。
假设晶体属于立方晶系,一个最简单的位错被称为edge位错,即在某一平面结束的晶胞中有一个平移,在平面另一侧开始另一个平移,造成缺陷的部分为一个原子。
由于材料中存在的缺陷引起了原子排列内的扰动,这会导致位错的能量、弛豫和移动特性与材料力学性能的变化。
计算位错的能量、弛豫和移动需要采用第一性原理方法中的密度泛函理论和分子动力学等方法。
第一性原理计算简述第一性原理,英文Firs t Principle,是一个计算物理或计算化学专业名词,广义的第一性原理计算指的是一切基于量子力学原理的计算。
我们知道物质由分子组成,分子由原子组成,原子由原子核和电子组成。
量子力学计算就是根据原子核和电子的相互作用原理去计算分子结构和分子能量(或离子),然后就能计算物质的各种性质。
从头算(ab initio)是狭义的第一性原理计算,它是指不使用经验参数,只用电子质量,光速,质子中子质量等少数实验数据去做量子计算。
但是这个计算很慢,所以就加入一些经验参数,可以大大加快计算速度,当然也会不可避免的牺牲计算结果精度。
那为什么使用“第一性原理”这个字眼呢?据说这是来源于“第一推动力”这个宗教词汇。
第一推动力是牛顿创立的,因为牛顿第一定律说明了物质在不受外力的作用下保持静止或匀速直线运动。
如果宇宙诞生之初万事万物应该是静止的,后来却都在运动,是怎么动起来的呢?牛顿相信这是由于上帝推了一把,并且牛顿晚年致力于神学研究。
现代科学认为宇宙起源于大爆炸,那么大爆炸也是有原因的吧。
所有这些说不清的东西,都归结为宇宙“第一推动力”问题。
科学不相信上帝,我们不清楚“第一推动力”问题只是因为我们科学知识不完善。
第一推动一定由某种原理决定。
这个可以成为“第一原理”。
爱因斯坦晚年致力与“大统一场理论”研究,也是希望找到统概一切物理定律的“第一原理”,可惜,这是当时科学水平所不能及的。
现在也远没有答案。
但是为什么称量子力学计算为第一性原理计算?大概是因为这种计算能够从根本上计算出来分子结构和物质的性质,这样的理论很接近于反映宇宙本质的原理,就称为第一原理了。
广义的第一原理包括两大类,以Hartr ee-Fork自洽场计算为基础的abinitio从头算,和密度泛函理论(DFT)计算。
分析能带图能带结构是目前采用第一性原理(从头abinitio)计算所得到的常用信息,可用来结合解释金属、半导体和绝缘体的区别.能带可分为价带、禁带和导带三部分,倒带和价带之间的空隙称为能隙,基本概念如图所示:如何能隙很小或为0 ,则固体为金属材料,在室温下电子很容易获得能量而跳跃至传倒带而导电;而绝缘材料则因为能隙很大(通常大于9电子伏特),电子很难跳跃至传导带,所以无法导电.一般半导体材料的能隙约为1至3电子伏特,介于导体和绝缘体之间。
因此只要给予适当条件的能量激发,或是改变其能隙之间距,此材料距能导电。
能带用来定性地阐明了晶体中电子运动的普遍特点。
价带(valence band),或称价电带,通常指绝对零度时,固体材料里电子的最高能量.在导带(conduction band)中,电子的能量范围高于价带,而所有在传导带中的电子均可经由外在的电场加速而形成电流。
对与半导体以及绝缘体而言,价带的上方有一个能隙(band gap),能隙上方的能带则是传导带,电子进入传导带后才能在固体材料内自由移动,形成电流.对金属而言,则没有能隙介于价带与传导带之间,因此价带是特指半导体与绝缘体的状况。
费米能级(fermi level)是绝对零度下的最高能级。
根据泡利不相容原理,一个量子态不能容纳两个或两个以上的费米子(电子),所以在绝度零度下,电子将从低到高依次填充各能级,除最高能级外均被填满,形成电子态的“费米海”。
“费米海”中每个电子的平均能量为(绝对零度下)为费米能级的3/5。
海平面即是费米能级.一般来说,费米能级对应态密度为0的地方,但对于绝缘体而言,费米能级就位于价带顶。
成为优良电子导体的先决条件是费米能级与一个或更多的能带相交。
能量色散(dispersion of energy)。
同一个能带内之所以会有不同能量的量子态,原因是能带的电子具有不同波向量(wave vector),或是k-向量。
在量子力学中,k—向量即为粒子的动量,不同的材料会有不同的能量-动量关系(E—K relationship).能量色散决定了半导体材料的能隙是直接能隙还是间接能隙.如导带最低点与价带最高点的K值相同,则为直接能隙,否则为间接能隙。
第二章第一性原理计算及相关理论方法固体能带理论是凝聚态物理中最成功的理论之一,是固体电子论的支柱。
原则上固体能带理论能够阐明和解释固体的许多基本物理性质,如力学,电学,光学及磁学等性质。
固体能带理论的主要任务是确定固体电子能级,也就是能带。
而要确定固体电子能带,其出发点便是求解组成固体的多粒子系统的定态薛定愕方程。
但是对固体这样有1029/m3数量级个原子核和电子的复杂的多粒子系统而言,其薛定愕方程是无法精确求解的,必须作近似简化求解:通过绝热近似将核运动和电子运动分开;通过Hartree-Fock 自洽场方法或更严格、更精确的密度泛函理论(DFT) 将多电子问题简化为单电子问题;通过将固体抽象为平移对称性的理想晶体,将多体问题归结为单电子在周期性势场中的运动。
而上述的这些近似简化最终把求解薛定愕方程的全部复杂性都归入了所谓的交换关联泛函,可见交换关联泛函在密度泛函理论中占有重要地位。
目前密度泛函理论已经成为探索具有重要应用背景的功能材料和结构材料的重要理论方法。
基于密度泛函理论,根据基函数选取的不同有多种具体的计算方法,通常都称为第一性原理计算(ab initio calculation)。
所谓第一性原理,即从最基本的物理规律出发,求解体系的薛定愕方程以获取材料性能的信息,从而理解材料中出现的一些现象,预测材料的性能。
第一性原理计算方法[79, 80]有着半经验(HF)方法不可比拟的优势,因为它只需要知道构成微观体系各元素的原子序数,而不需要任何可调(经验或拟合)参数,就可以基于量子力学来处理体系中的电子运动,来计算出该微观体系电子的波函数和对应的本征能量,从而求得系统的总能量、电子结构以及成键、弹性、稳定性等性质。
它们被广泛应用于原子,分子,固体,固体表面,界面,超晶格材料,低维材料的电子结构和物理性质的计算,并取得了惊人的成功。
随着最近几十年计算机技术的飞速发展,第一性原理的规模和效率都有了极大的提高。
目前计算物理已经从解析理论物理和实验物理中完全独立了出来,成为物理学中不可缺少的一个独立分支– 计算物理学[81]。
能带理论的周期性表象为布洛赫(Bloch)表象,同时还存在着局域化的瓦尼(Wannier)表象,Marzari N.等人[82, 83]发展了一套构造最局域Wannier 函数的数值方法,使得Wannier 轨道与第一性原理很好的结合,在传统的化学键分析和最近比较热的拓扑体系的非平庸的表面态计算上都有很大的应用。
传统的半导体模型分析k · p 方法[84],对于分析有效质量和光学性质有很大的帮助,同时它也一直是很好的能带近似手段,能准确的抓住物理实质,便于解析理论分析。
本章的主要任务是对密度泛函理论[85, 86],以及相应的数值的Wannier 方法[82, 83]和解析的k · p 能带[84]分析方法作简要介绍。
2.1 密度泛函理论(DFT)对于简单的晶体模型,描述多粒子体系的多体薛定谔方程可以写成如下形式: (),(,)H r R E r R ψ=ψ (2.1)其中r,R 分别表示所有电子和原子核的坐标。
不考虑外部势场时,其哈密顿量可以写为:2222222,1122||22||||I J I i Ii i j I I J i I e i j I I J i I Z Z e Z e e H m r r M R R r R ≠≠=-∇+-∇+----∑∑∑∑∑ (2.2) 右边第一项为电子动能部分,第二项为电子电子相互作用,第三项为原子核的 动能部分,第四项为原子原子相互作用,第五项为电子和原子核的相互作用。
这是一个复杂的多体哈密顿量,其希尔伯特空间(态空间)的维度与体系中的原子核和电子数目是指数相关的。
在实际体系中,总粒子数往往是巨大的。
因此尽管我们知道体系的各种性质是由薛定谔方程确定的,但是当总粒子数较大时,我们仍然不可能求解该方程,这样对于体系的性质我们仍然一无所知。
通过直接求解薛定谔方程来决定体系的各种性质, 仅限于用在理想模型或者简单的分子原子体系中,如氢分子H2,氦原子He 等。
对于实际的复杂体系而言,我们只能采取单电子近似。
下面介绍如何由晶格模型的多体哈密顿量到单电子方程。
2.1.1 从多体哈密顿量到单电子方程原则上,凝聚态物理的所有性质都已包括在(2.2)式中。
但是,要从对应的 薛定谔方程中得出所有的信息是不可能的。
既然直接求解是行不通的,人们提出了一些近似和模型来简化这个问题。
我们考虑一个具有N 个全同粒子的固体的简单模型,其中一定数量的近核电子可认为被紧紧的束缚在核的周围。
这些电子是定域的,因而对固体性质的贡献很小。
然而,外层轨道的电子可以是离域的,对体系的物理性质影响较大。
实际上,当原子结合在一起形成固体时,这些价电子(电子)的组态变化很大,而近核电子却比较固定,因此固体的电学和光学性质主要由价电子决定,这样我们就引入了价电子近似。
每个核和它的束缚电子看作一个整体,称为离子实(离子),固体可以看成是由价电子和离子实组合而成。
与此同时,对于电子系统,总体的电子波函数取决于作热振动的离子的瞬间位置。
电子运动的很快,以至于绝热的来适应缓慢运动的离子的影响。
另一方面,就离子的缓慢运动而言,快速运动的电子可以看成一个均匀的背景,这样我们就把运动的电子和振动的离子看成两个独立的子系统,这就是著名的Born-Oppenheimer 近似,或所谓的绝热近似。
这样,多粒子系统的Schrodinger 方程(2.1)的解就可以写成如下分离的形式: (,)()(,)r R R r R χψψ= (2.3)其中,()R χ描述的是离子实的运动;(,)r R ψ描述的是电子的运动。
在讨论后面电子近似时,我们主要关注电子的运动,离子的仅提供一个周期的晶格势场。
方便起见,在后面的电子波函数略去了对R 的依赖。
采取了绝热近似之后,我们的薛定谔方程仍然是多电子的,是一个多体问题,具体如下:2221(())()()22||i i i i j i e i j e r r E r m r r νψψ≠-∇++=-∑∑∑ (2.4)为了解决多体问题的需要,人们想到各种方法来把多电子方程变成单电子方程。
直到密度泛函理论的建立,才真正严谨地给出了单电子方程,为单电子方法建立了坚实物理理论基础。
这里我们简要的描述一下各种单电子近似方法的相继承关系。
最简单的近似就是自由电子气近似,这种近似认为金属中价电子是完全自由的,没有任何相互作用,尽管看上去这是很不切实际的想法,但是这种模型确实可以对简单金属的物理性质给以理论解释,得到了很好的结果。
最重要的是给出了金属的能谱,引出了重要的费米面的概念。
当然,由于这种模型的过于简单化,人们一直对这个结论心存怀疑,但是从后面提到的赝势理论可以看出,这种简化在一定条件下是完全合理的,只是这里得到的波函数应该认为是赝波函数,而在某种程度上赝波函数可以和真实波函数同样来描述这个系统。
哈特利(Hartree )近似[87]比自由电子气体近似多考虑了电子之间的库仑相互作用,把复杂的库仑作用看成是一种平均势。
哈特利多体波函数写成每个电子波函数()i i r ϕ连成的形式:112233()()()()()N N r r r r r ψϕϕϕϕ= (2.5)我们假定每个电子波函数i ϕ满足正交归一化条件,即是|i j ij ϕϕδ=。
然后把哈特利波函数的能量期待值对每个电子波函数i ϕ做变分,在电子波函数的正交归一条件下求能量的极小值。
这样就得到了单电子的哈特利方程:2222'|(')|(()')()()2|'|i i i i i i r r e dr r E r m r r ϕνϕϕ≠-∇++=-∑⎰ (2.6) 其中,()r ν表示单电子受到的晶格势,22'|(')|'|'|i i i r e dr r r ϕ≠-∑⎰表示单电子受到其它电子的平均库伦势。
所以他描写了r 处单个电子在晶格势()r ν和其他所有电子的平均势场中的运动。
这里单电子能量Ei 其实是变分的拉格朗日乘子。
我们可以看到平均库伦势明显地包含了单电子波函数()i r ϕ,所以解此类方程时必须自洽求解,这是该方法的特点。
首先,先假设电子处于一组特殊的状态,计算有效势,然后代入方程重新计算本征态,反复进行此过程,直到得到自洽结果。
哈特利近似明显没有考虑电子间交换相互作用,这个缺陷是提出哈特利-福克近似[88, 89]的出发点。
明显的修正体现在哈特利-福克多体波函数波函数上,他取单电子波函数φi 的Slater 行列式形式:1122212212()()()()()()1()!()()()i N i N N N N N q q q q q q q N q q q ϕϕϕϕϕϕψϕϕϕ= (2.7)i q 表示第i 个电子的坐标,包括位置i r 和自旋指标。
同样利用单电子波函数i ϕ为约束条件,对哈特利-福克波函数的能量期待值变分求极小值,得到如下单电子的哈特利-福克方程:2*2222'''',|||(')|(')(')(()')()'()()2|'||'|i i i i i i i i i i i r r r r e dr r e dr r E r m r r r r ϕϕϕνϕϕϕ≠≠-∇++-=--∑∑⎰⎰(2.8)很显然,哈特利-福克单电子方程(2.8)只计及电子间的交换作用,完全忽略了自旋反平行电子之间的相关能,这是这一近似的本质性缺陷。
虽然它在能带理论中长期被应用,但是由于这一本质性的缺陷,哈特利-福克近似不是单电子近似的严格理论。
后来提出的密度泛函理论突破了这个障碍,是严格的单电子理论。
2.1.2 密度泛函理论(DFT)对比传统的单电子自洽方程,密度泛函理论的不同之处是将单电子的密度函数作为基本变量,体系总能E 可以表示成电子密度函数的泛函。
20世纪六十年代,P. Hohenberg 和W. Kohn 在研究非均匀电子气理论时提出了著名的Hohenberg-Kohn(HK)定理[85],成功的建立起了密度泛函理论:定理一:不记自旋的全同费米子系统的基态能量是粒子数密度函数()r ρ的唯一泛函。
定理二:能量泛函在粒子数不变条件下在且仅在基态粒子数密度函数()o r ρ处取极小值,并等于基态能量。
此两定理的证明,可以在原始文献[85]中找到,这里不作赘述。
考虑N 个电子的相互作用系统,其哈密顿量可写成如下的两部分:int ext H H V =+(2.9) 其中:222int 122||ee i i i j e i j e H T V m r r ≠=+=-∇+-∑∑ (2.10)表示电子的动能项和电子之间的库伦相互作用项之和。