计算材料学概述.docVIP

计算材料学概述 计算材料学是基于物理建模与数值计算方法，通过理论计算主动对材料-器件-微系统的本征特性、结构与组分、使用性能以及合成与制备工艺进行综合设计，达到对材料结构与功能的调控，并提供优化设计和协同制造技术的一门交叉边缘学科。 1 密度泛函理论 密度泛函理论是一种研究多电子体系电子结构的量子力学方法。密度泛函理论在物理和化学上都有广泛的应用，特别是用来研究分子和凝聚态的性质，是凝聚态物理和计算化学领域最常用的方法之一。 在通常的多体问题电子结构的计算中，原子核可以看作静止不动的（波恩-奥本海默近似），这样电子可看作在原子核产生的静电势 中运动。电子的定态可由满足多体薛定谔方程的波函数 描述： 其中 为电子数目， 为电子间的相互作用势。算符 和 称为普适算符，它们在所有系统中都相同，而算符则依赖于系统，为非普适的。可以看出，单粒子问题和比较复杂的多粒子问题的区别在于交换作用项 。目前有很多成熟的方法来解多体薛定谔方程，例如：物理学里使用的图形微扰理论和量子化学里使用的基于斯莱特行列式中波函数系统展开的组态相互作用（CI）方法。然而，这些方法的问题在于较大的计算量，很难用于大规模复杂系统的计算。 相比之下，密度函理论将含 的多体问题转化为不含 的单体问题上，成为解决此类问题的一个有效方法。在密度泛函理论中，最关键的变量为粒子密度 ，它由下式给出 霍恩伯格和沃尔特·科恩在1964年提出 [1]，上面的关系可以反过来，即给出基态电子密度 ，原则上可以计算出对应的基态波函数 。也就是说， 是 的唯一泛函，即 对应地，所有其它基态可观测量 均为 的泛函 进而可以得出，基态能量也是 的泛函 , 其中外势场的贡献 可以用密度表示成 泛函 和 称为普适泛函，而 显然不是普适的，它取决于所考虑的系统。对于确定的系统，即 已知，需要将泛函 对于 求极小值。这里假定能够得出 和 的表达式。对能量泛函求极值可以得到基态能量 ，进而求得所有基态可观测量。 对能量泛函 求变分极值可以用不定算子的拉格朗日方法，这由科恩和沈吕久在1965年完成 [2]。这里我们使用如下结论：上面方程中的泛函可以写成一个无相互作用的体系的密度泛函 其中 为无相互作用的动能， 为粒子运动感受到的外势场。显然， ，若 取为 这样，可以解这个辅助的无相互作用体系的科恩-沈吕久方程 可以得到一系列的电子轨域 ，并由此求得原来的多体体系的电子密度 等效的单粒子势 可以表示成 其中第二项为描述电子间库仑斥力的哈特里项，最后一项 叫做交换关联势，包含所有多粒子的相互作用。由于哈特里项和交换关联项 都依赖于 , 又依赖于 , 而 又依赖于 , 科恩-沈吕九方程的求解需要用自洽方法。通常首先假设一个初始的 , 然后计算对应的 并求解科恩-沈吕九方程中的 。进而可以计算出新的密度分布，并开始新一轮计算。此过程不断重复，直到计算结果收敛。 2 分子动力学 分子动力学是一套HYPERLINK /wiki/%B7%D6%D7%D3%C4%A3%94M分子模拟方法，该方法主要是依靠HYPERLINK /wiki/%C5%A3%B6%D9%C1%A6%D1%A7牛顿力学来模拟分子体系的运动，以在由分子体系的不同状态构成的HYPERLINK /wiki/%CF%B5%D7%DB系综中抽取样本，从而计算体系的HYPERLINK /w/index.php?title=%B9%B9%D0%CD%BB%FD%B7%D6action=editredlink=1构型积分，并以构型积分的结果为基础进一步计算体系的热力学量和其他宏观性质。 基本步骤为： (1)确定起始构型 进行分子动力学模拟的第一步是确定起始HYPERLINK /w/index.php?title=%B9%B9%D0%CDaction=editredlink=1构型， 一个能量较低的起始构型是进行分子模拟的基础 ，一般分子的起始构型主要来自实验数据或HYPERLINK /wiki/%C1%BF%D7%D3%BB%AF%D1%A7量子化学计算。 在确定起始构型之后要赋予构成分子的各个原子速度，这一速度是根据HYPERLINK /wiki/%C2%F3%BF%CB%CB%B9%CE%A4-%B2%A3%B6%FB%D7%C8%C2%FC%CD%B3%BC%C6玻尔兹曼分布随机生成的，由于速度的分布符合玻尔兹曼统计，因此在这个阶段，体系的HYPERLINK /wiki/%CE%C2%B6%C8温度是恒定的。另外，在随机生成各个原子的运动速度之后须进行调整，使得体系总体在各个方向上的动量之和为零，即保证体系没有HYPERLINK /w/index.php?title=%C6%BD%