第一性原理计算判断材料稳定性的几种方法
材料科学中第一性原理计算方法研究
材料科学中第一性原理计算方法研究近年来,材料科学领域的研究取得了许多重大突破,其中第一性原理计算方法成为材料设计和研究的重要工具之一。
这种方法通过基本的物理原理和数学方程来研究材料的性质和行为,为材料设计和性能优化提供了新的途径。
第一性原理计算方法是基于量子力学的一种计算方法,从第一性原理出发,通过求解薛定谔方程以及其他相关方程来研究材料的性质。
它不依赖于任何经验参数或假设,能够提供对材料的精确描述和准确预测。
第一性原理计算方法的核心是密度泛函理论(Density Functional Theory,简称DFT),它将体系的物理性质与体系中电子的密度联系起来。
根据Kohn-Sham方程,DFT通过对电子的运动方程进行求解,得到体系的基态电子密度。
通过计算得到的电子密度,可以进一步计算出材料的能带结构、电子态密度、态密度、声子谱、磁性及其它性质。
与传统的实验方法相比,第一性原理计算方法具有独特的优势。
首先,它能够提供物理性质的原子尺度描述,可以捕捉到材料内部微观原子结构的信息。
其次,该方法能够计算和预测材料的多种性质,如电子能带结构、晶格常数、弹性性能、热力学性质等,为材料设计和开发提供了重要参考。
此外,第一性原理计算方法可以帮助解释材料性能背后的基本物理机制,揭示材料特性的微观本质。
近年来,随着计算机性能的不断提升和计算方法的进步,第一性原理计算方法在材料科学中的应用得到了广泛拓展。
例如,它在材料的合成、器件的设计和材料的特性优化等方面发挥了重要作用。
通过预测和优化材料的能带结构和电子态密度,可以筛选出具有优异性能的新材料,为新能源、环境友好材料、传感器和光电器件的研发提供重要支持。
此外,第一性原理计算方法还可以帮助优化材料的力学、热力学和电磁性能,提高材料的功能性能。
尽管第一性原理计算方法为材料科学提供了强大的工具和理论基础,但也面临一些挑战。
首先,该方法对计算所需的资源要求较高,需要大量计算时间和计算机内存。
材料科学中的相变与相稳定性研究
材料科学中的相变与相稳定性研究相变和相稳定性是材料科学中重要的研究方向,它们对于理解材料的性质和应用具有关键作用。
相变是指物质在一定条件下从一种相态转变为另一种相态的过程,而相稳定性则是指物质在一定条件下保持某种相态的能力。
本文将从理论和实验两个方面介绍材料科学中相变与相稳定性的研究进展。
一、理论研究材料的相变和相稳定性研究常常需要依靠理论模型和计算方法。
在相变研究中,一种常见的理论模型是自由能模型。
自由能是描述物质系统能量的函数,通过计算不同相态的自由能变化,可以确定相变发生的条件和过程。
此外,晶格动力学模型也被广泛应用于相变研究中,它考虑了晶格振动对物质结构的影响,并能够模拟相变的动力学行为。
对于相稳定性的研究,理论计算方法在材料科学中发挥着重要的作用。
例如,通过密度泛函理论等第一性原理计算方法可以计算材料的能带结构、原子结构和电子态密度等物理性质,从而预测材料的相稳定性和性能。
此外,相场模型和微观动力学模拟等方法也可用于研究材料的相稳定性。
这些理论模型和计算方法为相变和相稳定性的研究提供了有力的理论支持。
二、实验研究材料科学中的相变和相稳定性也需要通过实验手段进行研究。
实验方法主要包括热分析、X射线衍射、电子显微镜等。
热分析是一种常用的方法,可以通过对材料加热或降温过程中的热量变化进行测量,来研究相变的发生温度、吸放热等性质。
X射线衍射可以确定材料的晶体结构和晶格常数,进而研究相变过程中晶体结构的变化。
电子显微镜则可以观察材料微观结构的演变过程,揭示相变的原位特性。
在实验研究中,科学家们还开展了很多有关相变和相稳定性的前沿研究。
例如,自发晶化相变研究表明,一些非晶态材料在一定条件下会发生自发结晶,从而形成具有特殊性质的晶态材料。
此外,通过调控材料中的缺陷结构,科学家们也成功实现了一些新型相态的实验构建,为相变和相稳定性的研究提供了新的思路和方法。
三、在材料应用中的意义相变和相稳定性的研究对于材料科学的应用具有重要意义。
Y_掺杂MgZn2_稳定性、电子结构和力学性能的第一性原理计算
表面技术第52卷第8期Y掺杂MgZn2稳定性、电子结构和力学性能的第一性原理计算陈敬昶1,刘建国2,陆绍敏1,钟丽辉1,陈丽萍1,王远1(1.西南林业大学,昆明 650224;2.云南建投第七建设有限公司,昆明 650100)摘要:目的稀土元素Y掺杂是改善7xxx系铝合金断裂韧性的重要途径,然而因其掺杂量极低,通过实验很难测定微量Y对7xxx系铝合金析出相及强韧机制产生的作用,限制了7xxx系铝合金的进一步发展。
采用第一性原理计算方法探究Y掺杂对7xxx系铝合金中重要析出相MgZn2的影响机理,为7xxx系铝合金的微合金化强韧机理研究提供理论依据。
方法构建适于第一性原理计算、Mg/Zn的原子数分数比为1∶2的晶体模型,Y原子通过替换Mg或Zn原子的方式进行掺杂,通过能量计算、电子计算和弹性常数计算等分析Y掺杂对MgZn2能量稳定性、电子结构和力学性能的影响机理。
结果经Y掺杂后,形成3种固溶体Mg3Zn8Y、Mg4Zn7Y-1和Mg4Zn7Y-2,它们的形成热均小于0,即它们均可自发形成且稳定存在。
通过结合能计算发现,3种固溶体的结合能都小于MgZn2的结合能,说明Y掺杂促进了MgZn2的稳定性。
通过电子结构分析发现,Y掺杂后与Mg、Zn原子形成强的共价键,增强了体系的稳定性,Mg-Zn原子间形成了强离子键,MgZn2中Zn-Zn原子间的共价键变为强离子键。
力学性能计算结果表明,经Y掺杂后MgZn2的硬度降低、韧性上升,即Y掺杂增强了7xxx系铝合金中重要弥散析出相MgZn2的韧性,从而提升了7xxx 系铝合金的断裂韧性和抗疲劳能力。
结论基于计算结果分析得出,Y掺杂提升了MgZn2的稳定性、键合强度和断裂韧性,相关计算分析为微量Y掺杂增强7xxx系铝合金断裂韧性的实验分析提供了指导。
关键词:MgZn2;Y掺杂;电子结构;力学性能;第一性原理计算中图分类号:TG135 文献标识码:A 文章编号:1001-3660(2023)08-0444-07DOI:10.16490/ki.issn.1001-3660.2023.08.040First-principles Calculations of Stability, Electronic Structure andMechanical Properties of Y-doped MgZn2CHEN Jing-chang1, LIU Jian-guo2, LU Shao-min1, ZHONG Li-hui1, CHEN Li-ping1, WANG Yuan1(1. Southwest Forestry University, Kunming 650224, China; 2. YYCIN No.7 Construction Co., Ltd., Kunming 650100, China)ABSTRACT: Micro-alloying of rare earth element Y is an important way to strengthen the fracture toughness of 7xxx series收稿日期:2022-07-20;修订日期:2023-02-16Received:2022-07-20;Revised:2023-02-16基金项目:国家自然科学基金(51301144);云南省教育厅科学研究基金(2022Y574)Fund:National Natural Science Foundation of China (51301144); Scientific Research Foundation of Education Department of Yunnan Province (2022Y574)作者简介:陈敬昶(1996—),男,硕士,主要研究方向为材料表面改性。
材料热稳定性评估方法总结
材料热稳定性评估方法总结材料的热稳定性是指材料在高温或长时间暴露下的保持稳定性能和不发生明显物理或化学变化的能力。
热稳定性评估方法的选择对于材料的开发、制备和应用至关重要。
本文将综述几种常见的材料热稳定性评估方法,包括热重分析法、差示扫描量热法、动态热机械分析法、厨师自燃法和氧指数测定法。
热重分析法(Thermogravimetric Analysis, TGA)是一种广泛应用于材料热稳定性评估的常用方法。
该方法通过在恒定升温速率下测量样品的质量变化,来研究材料在不同温度下的热分解、挥发、燃烧等行为。
热重分析法可以定量得到材料的热分解温度、热分解速率、残渣含量等参数,进而评估材料的热稳定性。
这种方法具有操作简便、测量精度高的优点,适用于各种材料的热稳定性评估。
差示扫描量热法(Differential Scanning Calorimetry, DSC)是一种常见的用于研究材料热性质的方法,也可用于热稳定性评估。
该方法通过测量样品与参比物之间的温度差异和吸热/放热效应来分析材料的热分解、熔融等行为。
差示扫描量热法可以得到材料的熔点、熔融焓、热分解焓等参数,进而评估材料的热稳定性。
这种方法具有灵敏度高、分辨率好的优点,适用于大多数材料的热稳定性评估。
动态热机械分析法(Dynamic Mechanical Analysis, DMA)是一种通过在恒定频率或恒定应变下测量材料的动态力学性能来评估材料热稳定性的方法。
该方法可以测定材料的弹性模量、损耗因子、玻璃化转变温度等参数,以及材料在不同温度下的力学性能变化。
动态热机械分析法可以评估材料的粘弹性行为和蠕变行为,进而判断材料的热稳定性。
这种方法具有测试频率范围广、测试结果可靠的优点,适用于研究材料的热稳定性。
厨师自燃法(Cook's Self-ignition Test)是一种常见的用于评估材料热稳定性的方法。
该方法将样品置于恒定温度条件下,观察样品的自燃或燃烧表现。
第一性原理计算方法的应用
第一性原理计算方法的应用1. 简介第一性原理计算方法是一种基于量子力学原理的计算方法,用于模拟材料的性质和相互作用。
该方法基于物理学的基本原理,通过求解薛定谔方程来计算材料的电子结构和材料的相互作用。
2. 第一性原理计算方法的原理第一性原理计算方法的关键是通过求解薛定谔方程来模拟材料的电子结构。
薛定谔方程描述了系统的整体波函数和能量。
通过求解薛定谔方程可以得到材料的电子布居和能带结构,从而了解其性质和相互作用。
常用的第一性原理计算方法包括密度泛函理论 (Density Functional Theory, DFT)、平面波赝势方法 (Plane Wave Pseudopotential Method, PWPP) 等。
3. 第一性原理计算方法的应用领域第一性原理计算方法在材料科学、物理化学和生物化学等领域有广泛的应用。
下面列举了一些常见的应用领域:•材料设计:通过第一性原理计算方法可以预测材料的性质,并设计出具有特定性能的材料。
例如,可以通过计算得到材料的晶格常数、弹性常数、能带结构等,从而实现对材料性质的精确控制。
•催化剂研究:第一性原理计算方法可以用于研究催化剂的活性和选择性。
通过计算能够揭示催化剂的电子结构和表面吸附特性,进而提供理论指导和预测。
•表面和界面科学:表面和界面对材料的性能和反应具有重要影响。
第一性原理计算方法可以揭示材料表面和界面的结构、能量和化学反应动力学,有助于理解材料表面的吸附、扩散和反应机理。
•电池材料:第一性原理计算方法可以用于研究电池材料的离子扩散、电荷转移和储能机理。
通过计算,可以揭示材料的析出机制、电化学性能和循环稳定性,为电池材料设计提供理论指导。
•生物分子模拟:第一性原理计算方法可以用于模拟生物分子的结构和相互作用。
通过计算能够揭示生物分子的构象和能量,从而理解生物分子的功能和生物过程。
4. 第一性原理计算方法的优势和挑战第一性原理计算方法具有以下优势:•精确性:第一性原理计算方法基于基本物理原理,能够提供精确的计算结果,有助于深入理解材料的性质和行为。
第一性原理计算方法在材料科学中的应用
第一性原理计算方法在材料科学中的应用1.引言第一性原理计算方法(First Principles Calculation)是近年来发展的新型计算方法,用于准确计算分子和固体物质的能量、结构和物理性质。
它的优势在于不依赖于实验数据,可以直接从基本原理推导出体系的特性。
在材料科学领域,第一性原理计算方法已经成为研究材料的重要工具,可以为合成新材料和设计功能材料提供理论依据,并指导实验研究。
2.第一性原理计算方法的基本原理第一性原理计算方法的基本原理是量子力学中的密度泛函理论,它的基本假设是所有粒子的运动都可以描述为波函数的运动。
根据波函数理论,一个由N个电子和原子核组成的体系的波函数可以用N个单电子波函数表示。
通过求解薛定谔方程,可以确定体系的基态能量和电子的密度,从而得到体系的性质。
3.第一性原理计算方法在材料科学中的应用(1)材料合成第一性原理计算方法可以模拟材料的结构、动力学和化学反应,为材料合成提供理论指导。
例如,使用第一性原理计算方法可以预测材料的稳定性、生长机制和晶体缺陷,从而为材料的设计和制备提供指导。
(2)材料性能第一性原理计算方法可以计算材料的电子结构、热力学性质、光电性质和磁学性质等,从而为材料的性能研究提供理论基础。
例如,通过计算材料的电子结构,可以预测材料的导电性、热导率和热电性能等,为相关应用提供指导。
(3)材料改性第一性原理计算方法可以模拟材料的界面和表面结构,研究材料的改性效果。
例如,可以通过计算材料与其他材料的界面能量来评估材料的附着性和界面稳定性,从而指导材料的改性设计。
(4)功能材料设计借助第一性原理计算方法,可以针对具体的应用需求,设计出具有特定功能的材料。
例如,通过计算材料的光电性质、催化活性和磁学性质等,可以指导材料的功能设计,为实现特定的应用提供理论指导。
4.发展趋势随着材料科学和计算科学的发展,第一性原理计算方法的应用前景越来越广阔。
未来,第一性原理计算方法将会与机器学习和高通量计算等技术结合,为材料科学的研究提供更多的可能性。
稳定性分析结构的稳定性判断与计算方法
稳定性分析结构的稳定性判断与计算方法稳定性分析在结构工程中具有重要的意义,它用于评估结构在受力情况下的稳定性和可靠性。
本文将讨论结构的稳定性判断和计算方法,并介绍一些常用的工程实践。
一、稳定性判断方法1. 静力刚度法静力刚度法是最简单且常用的稳定性判断方法之一。
该方法基于结构在稳定状态下,受力平衡和变形满足静力学方程的假设。
根据结构的初始几何形状和受力情况,可以得到结构的初始刚度矩阵。
通过判断结构的刚度矩阵的特征值是否为正,可以确定结构的稳定性。
2. 弹性屈曲分析法弹性屈曲分析法是一种精确的稳定性判断方法,适用于具有复杂几何形状和较大位移的结构。
该方法基于弹性力学原理,通过对结构的弹性刚度矩阵进行特征值分析,得到结构的屈曲荷载和屈曲模式。
如果结构在设计荷载下的实际荷载小于屈曲荷载,那么结构就是稳定的。
3. 极限平衡法极限平衡法是一种基于能量平衡原理的稳定性分析方法。
该方法通过建立稳定状态下结构的能量平衡方程,利用极限状态下的能量变化来判断结构的稳定性。
当结构受到外力作用时,如果能量平衡方程能够满足,那么结构就是稳定的。
否则,结构将失去稳定性。
二、稳定性计算方法1. 弯曲稳定性计算在结构设计中,弯曲稳定性是最常见的稳定性问题之一。
弯曲稳定性计算可以通过欧拉公式进行。
欧拉公式是计算压杆稳定性的经典方法,它可以用来计算弯曲后的截面失稳荷载。
根据欧拉公式,弯曲稳定性计算可以通过截面惯性矩、截面形状和截面材料的参数来进行。
2. 局部稳定性计算除了弯曲稳定性,局部稳定性也是一个重要的考虑因素。
局部稳定性通常涉及到薄弱的结构构件,如薄壁构件和薄板。
局部稳定性计算可以通过截面失稳计算、临界载荷计算和局部屈曲分析来进行。
这些方法可以帮助设计人员确定结构是否足够抵抗局部失稳的力量。
三、工程实践1. 结构稳定性设计在结构设计中,稳定性是一个基本的要求。
设计人员需要根据结构的空间几何形状、荷载情况和材料特性,综合考虑弯曲稳定性和局部稳定性。
castep动力学计算
castep动力学计算CASTEP是一种基于密度泛函理论(DFT)的第一性原理计算软件,用于模拟和预测材料的物理和化学性质。
它是由剑桥大学的理论化学家Richard Needs开发的,能够提供高精度和可靠的计算结果。
CASTEP可以用于研究固体材料、表面、分子等多种系统,并有广泛的应用领域,如能源材料、催化剂、纳米材料等。
动力学计算是CASTEP的一个重要功能之一。
动力学计算能够模拟计算体系在一定时间范围内的演化过程,从而得到材料的结构、能量、振动、热力学性质等信息。
通过动力学计算,我们可以了解材料的物理和化学性质如何随时间变化,推断出材料的稳定性、相变、催化反应等信息。
在CASTEP中进行动力学计算通常需要以下步骤:首先,需要确定计算模型和参数。
这一步骤包括选择所研究的系统,确定原子种类、初态、边界条件等,并设定计算的时间尺度和步长。
第二步是进行能量最小化计算。
这一步骤通过优化系统的几何结构来找到材料的平衡构型。
在能量最小化计算中,CASTEP使用了Broyden-Fletcher-Goldfarb-Shanno(BFGS)方法来寻找能量最小值点。
完成能量最小化计算后,就可以开始进行动力学计算了。
CASTEP使用Verlet算法来求解牛顿运动方程。
Verlet算法是一种基于位置和速度的运算方法,通过迭代计算来得到系统在给定时间内的演化过程。
在动力学计算中,CASTEP可以计算各种物理量的演化情况,包括原子的运动轨迹、动能、势能、温度、压强等。
这些物理量的计算可以帮助我们了解材料的热力学性质、力学性能等。
动力学计算的结果需要通过对数据进行统计和分析来得出结论。
CASTEP提供了强大的后处理分析工具,如计算系统的结构动力学、参考能量和力学性质等。
总的来说,CASTEP的动力学计算功能使科学家能够模拟和研究材料的演化过程,为材料科学的研究和应用提供了重要的工具。
通过动力学计算,我们可以更好地了解材料的物理和化学性质,预测材料的性能、相变等信息,为设计和改进新材料提供指导和帮助。
第一性原理计算CaNi2和CaMgNi4合金的相稳定性
2 计 算 模 型 和 方 法
计算 由 Maeil tdo 3 1中 的 C se tr u i . aS atp软 件 包 完 成 , atp软件 是一 种 基 于 密度 泛 函 方法 的从 头算 C se 量 子力 学程序 , 它利 用 总能量平 面波 赝势 方法 , 离子 将
的收敛标 准设 为 1 o e 作 用 于每 个 原 子 上 的力 ×1 V,
低 于 0 0 e n 内 应 力 低 于 0 0 G a 位 移 < . 3 V/ m, . 3P, 00 0n . 0 1 m。开始计 算时 , 根据 实验数 据口。 立计算 初 。 建
数 ) ln [ 。 i g3 使 用 机 械 合 金 技 术 制 备 了 C , a 等 a一
势用 赝 势替代 , 电子 波 函数通过 平 面波基 组展 开 , 电子
一
电 子 相 互 作 用 的交 换 和 相 关 势 由 局 域 密 度 近 似
关 键 词 : C M g 相 稳 定 性 ; 性 ; 一 性 原 理 a Ni; 弹 第 中 图 分 类 号 : TG1 9 7 T 4 . 文 献 标 识 码 : 3 . ; G1 6 4 A
定 。弹性 系数计 算结 果 显 示, a Ni C Mg 4的体 模 量远 大
于 C Ni。利 用 弹 性 系数 、 模 量 的 经 验 公 式 计 算 了 a 。 体
合 金 的 熔 点 , 果 显 示 , Mg 的 熔 点 ( 3 9 6 或 结 Ca Ni 17. K 1 5 . K) 于 C Ni 熔 点 ( 0 5 4 或 9 9 3 K) 4 84 大 a 2的 1 7.K 7. 7 。
第一性原理的应用ms
第一性原理的应用什么是第一性原理第一性原理指的是通过独立思考,从基本的事实和原理出发,而非依赖于经验和传统的做法来解决问题的方法。
它是一种基于基本的定律和原则,从头开始推导出结论的方法。
在科学领域,第一性原理方法被广泛应用于物理学、化学、材料科学等领域,因为它可以提供最准确和可靠的结果。
然而,近年来,第一性原理也被越来越多地应用于其他领域,如工程设计、数据分析、经济预测等。
第一性原理在材料科学中的应用•原子结构分析:通过第一性原理计算,可以精确地计算材料的原子结构和晶格参数,从而了解材料的性质和行为。
•能量计算:第一性原理方法可以准确地计算材料的能量,包括电子能量、晶格能、声子能等,从而评估材料的稳定性和反应性。
•界面模拟:利用第一性原理方法,可以模拟和预测材料之间的界面和界面上的反应,从而指导材料的设计和合成。
•总能量优化:通过第一性原理计算,可以确定材料的最低能量构型和稳定结构,从而优化材料的性能和性质。
第一性原理在工程设计中的应用•材料选择:通过第一性原理计算,可以评估不同材料的性能和特性,从而选择最适合的材料用于特定的工程设计。
•设备优化:利用第一性原理方法,可以优化工程设备的结构和性能,提高工艺效率和产品质量。
•反应预测:通过第一性原理计算,可以预测化学反应的速率和产物,指导实验操作和工艺设计。
•材料性能预测:基于第一性原理方法,可以预测材料的电子结构、光学性质、磁性等特性,为材料的设计和开发提供重要的参考。
第一性原理在数据分析中的应用•数据处理:第一性原理方法可以用于处理和分析大规模数据,从中提取出有用的信息和模式。
•特征选择:通过第一性原理计算,可以确定数据中最重要和相关性最强的特征,减少数据维度和复杂度。
•模型训练:基于第一性原理的方法,可以利用数据生成和训练模型,用于预测和分类任务。
•数据挖掘:第一性原理方法可以用于挖掘隐藏在数据背后的规律和趋势,为决策和规划提供支持。
第一性原理在经济预测中的应用•宏观经济模型:通过第一性原理建立宏观经济模型,对经济的发展趋势和影响因素进行分析和预测。
第一性原理计算方法讲义
第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。
而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。
第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。
量子力学是反映微观粒子运动规律的理论。
量子力学的出现,使得人们对于物质微观结构的认识日益深入。
原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。
量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。
以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。
目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。
但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。
绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。
Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。
但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。
1964年,Hohenberg和Kohn提出了严格的密度泛函理论(Density Functional Theory, DFT)。
它建立在非均匀电子气理论基础之上, 作为基本变量。
计算机在材料科学与工程中的应用题集
计算机在材料科学与工程中的应用题集一、选择题1. 在材料科学中,正交实验设计主要用于:A. 预测材料性能B. 单一变量研究C. 复杂系统建模D. 多因素交互作用分析答案: D. 多因素交互作用分析2. 正交表在正交实验设计中的主要作用是:A. 确定实验次数B. 预测实验结果C. 安排实验因素水平D. 评估实验误差答案: C. 安排实验因素水平3. 下列哪项不是正交实验设计的优点?A. 减少实验次数B. 提高实验效率C. 精确控制单一变量D. 便于分析多因素交互效应答案: C. 精确控制单一变量(正交实验设计主要关注多因素交互,而非单一变量)4. 在使用正交实验设计软件时,首先需要:A. 设定实验结果预期B. 确定实验因素及水平C. 分析历史数据D. 选择合适的正交表答案: B. 确定实验因素及水平5. 正交实验设计中的“水平”指的是:A. 实验的重复次数B. 实验因素的不同取值C. 实验结果的精确度D. 实验设备的型号答案: B. 实验因素的不同取值6. 下列哪项是分析正交实验结果时常用的方法?A. 方差分析B. 回归分析C. 因子分析D. 以上都是答案: D. 以上都是7. 正交实验设计能够有效减少实验次数,主要是因为:A. 减少了实验误差B. 简化了实验流程C. 利用正交表合理安排实验D. 提高了实验结果的可靠性答案: C. 利用正交表合理安排实验8. 正交实验设计中,如果某个因素有3个水平,另一个因素有2个水平,那么至少需要进行的实验次数是:A. 2B. 3C. 5D. 6答案: D. 6(正交表的选择会基于各因素的最大水平数,通常需要选择能容纳所有因素及其水平的最小正交表)9. 在分析正交实验结果时,如果某个因素的极差(Range)最大,说明:A. 该因素对实验结果影响最小B. 该因素对实验结果影响最大C. 该因素与其他因素无交互作用D. 实验结果完全由该因素决定答案: B. 该因素对实验结果影响最大10. 正交实验设计的结果分析不包括以下哪项内容?A. 计算各因素的极差B. 绘制因素与指标的关系图C. 预测未来实验结果D. 分析多因素间的交互作用答案: C. 预测未来实验结果(正交实验设计主要用于分析现有实验数据,而非预测未来结果)11. 在选择正交表时,以下哪个因素不是主要考虑的?A. 实验因素的数量B. 各因素的水平数C. 实验结果的精确度要求D. 实验人员的个人喜好答案: D. 实验人员的个人喜好(正交表的选择应基于实验设计的科学性和统计要求)12. 正交实验设计在材料改性研究中,常用于:A. 确定最佳改性条件B. 评估改性材料的市场价值C. 预测改性材料的寿命D. 分析改性材料的微观结构答案: A. 确定最佳改性条件13. 当正交实验设计的结果显示两个因素之间存在显著的交互作用时,应进一步:A. 忽略该交互作用B. 仅分析其中一个因素C. 进行交互作用图的绘制与分析D. 增加实验次数以确认结果答案: C. 进行交互作用图的绘制与分析14. 在材料科学研究中,采用正交实验设计后,通常需要进一步进行:A. 重复实验验证B. 单一变量补充实验C. 数据拟合分析D. 以上都是答案: D. 以上都是15. 在材料力学分析中,用于模拟复杂应力状态下材料行为的常用数值方法是?A. 有限元法B. 有限差分法C. 蒙特卡洛模拟D. 分子动力学模拟答案: A. 有限元法16. 哪种软件常用于材料科学中的三维结构模拟与优化?A. MATLABB. SolidWorksC. COMSOL MultiphysicsD. Autodesk Inventor答案: C. COMSOL Multiphysics17. 在进行材料断裂力学分析时,主要关注的是哪个参数?A. 弹性模量B. 应力强度因子C. 泊松比D. 屈服强度答案: B. 应力强度因子18. 下列哪项技术可用于预测材料在极端条件下的性能变化?A. 密度泛函理论B. X射线衍射分析C. 分子动力学模拟D. 扫描电子显微镜答案: C. 分子动力学模拟19. 在材料结构分析中,哪种方法能够直接观察材料的内部微观结构?A. 有限元分析B. 透射电子显微镜C. 数值模拟D. 红外光谱分析答案: B. 透射电子显微镜20. 材料科学中,用于模拟材料在加载条件下应力分布的软件通常是?A. AutoCADB. ANSYSC. OriginD. SPSS答案: B. ANSYS21. 哪种分析技术可以评估复合材料中各组分间的界面结合强度?A. 纳米压痕测试B. 扫描隧道显微镜C. 微观力学模型D. 拉曼光谱分析答案: C. 微观力学模型22. 在材料疲劳分析中,为了评估材料的寿命,常用的方法是?A. 应力-应变曲线分析B. 疲劳裂纹扩展速率测试C. 硬度测试D. 热重分析答案: B. 疲劳裂纹扩展速率测试23. 下列哪个软件常用于材料的相图计算和热力学模拟?A. MATLABB. Thermo-CalcC. SolidWorksD. Gaussian答案: B. Thermo-Calc24. 在进行材料的热应力分析时,主要考虑的是材料的哪项性质?A. 导热系数B. 弹性模量C. 密度D. 熔点答案: A. 导热系数25. 材料科学中,模拟材料在腐蚀环境下的行为常用哪种方法?A. 蒙特卡洛模拟B. 有限元腐蚀模拟C. 电化学测试D. 密度泛函理论答案: B. 有限元腐蚀模拟(注意:虽然实际中“有限元腐蚀模拟”不是标准术语,但这里为了题目设计而使用,意在表达使用有限元方法进行腐蚀行为的模拟)26. 在材料科学研究中,为了分析材料的微观缺陷,常用的技术是?A. 超声波检测B. 透射电子显微镜C. 红外热成像D. 激光粒度分析答案: B. 透射电子显微镜27. 在进行材料的力学性能测试时,用于评估材料韧性的主要指标是?A. 硬度B. 弹性模量C. 冲击韧性D. 屈服强度答案: C. 冲击韧性28. 在模拟材料的高温蠕变行为时,关键考虑的是材料的哪个性质?A. 弹性模量B. 蠕变极限C. 断裂韧性D. 熔点答案: B. 蠕变极限29. 下列哪种技术用于分析材料在加载过程中的变形和应力分布?A. 扫描电子显微镜B. 数字图像相关法(DIC)C. 红外光谱D. 能量色散X射线光谱(EDS)答案: B. 数字图像相关法(DIC)30. 在材料科学中,为了模拟材料的热传导过程,常用的软件是?A. ANSYSB. GaussianC. MATLABD. AutoCAD答案: A. ANSYS31. 第一性原理计算中,用于描述材料电子结构的理论框架主要是?A. 密度泛函理论B. 经典力学C. 统计力学D. 量子力学答案: D32. 在第一性原理计算中,哪个参数对于描述固体的能带结构至关重要?A. 晶胞大小B. 原子间距C. 离子电荷D. 布里渊区答案: D33. 下列哪个软件常用于第一性原理计算的密度泛函理论模拟?A. MATLABB. AutoCADC. GaussianD. VASP答案: D34. 第一性原理计算中,用于近似处理多电子体系中电子间相互作用的常见方法是?A. 玻尔模型B. 哈特里-福克方法C. 密度矩阵方法D. 局域密度近似(LDA)或广义梯度近似(GGA)答案: D35. 哪种类型的计算在第一性原理模拟中通常用于研究材料的表面和界面性质?A. 分子动力学模拟B. 蒙特卡洛模拟C. 量子蒙特卡洛模拟D. 平板模型计算答案: D36. 在第一性原理计算中,为了获得更准确的电子结构信息,经常需要优化哪个参数?A. 原子质量B. 离子半径C. 晶格常数D. 截断能答案: D37. 下列哪个物理量在第一性原理计算中通常用于描述材料的磁性?A. 电导率B. 介电常数C. 热导率D. 磁矩答案: D38. 在进行第一性原理计算时,哪个步骤通常涉及求解Kohn-Sham方程?A. 初始化参数B. 能量最小化C. 自治场迭代D. 能带结构分析(尽管不直接求解K-S方程,但C项更接近实际求解过程)答案: C(注意:D项不直接相关,但C项是求解K-S方程的核心步骤)39. 第一性原理计算中,用于描述材料光学性质的物理量通常是什么?A. 折射率B. 电阻率C. 介电函数D. 复介电常数答案: D40. 下列哪个算法在第一性原理计算中常用于处理材料的声子谱?A. 密度泛函微扰理论(DFPT)B. 分子动力学C. 蒙特卡洛方法D. 超胞法结合有限位移法答案: A(但D项也是实际操作中可能用到的方法之一,但A项更直接相关)(注意:第40题的D选项虽然在实际中可能用于计算声子谱,但A 选项的DFPT是第一性原理计算中更直接用于声子谱计算的方法。
材料科学中的第一性原理计算
材料科学中的第一性原理计算第一性原理计算是材料科学研究中一种重要的计算方法。
它是基于量子力学理论和电子结构理论的计算模型,通过求解薛定谔方程,从基本粒子(原子、离子、电子)的特性出发,利用数学方法预测和描述材料的结构、能量、性质等基本信息。
本文将对第一性原理计算的原理、方法和应用进行详细介绍。
第一性原理计算的核心是量子力学。
量子力学是描述微观粒子行为的理论,它认为微观粒子的运动和相互作用需要用波函数描述,而波函数可以通过薛定谔方程求解。
在材料科学中,我们关心的是材料中电子的结构和性质。
通过解薛定谔方程,可以得到材料中电子的轨道分布、能带结构和电子密度等信息,进而预测和研究材料的各种性质。
第一性原理计算分为两个主要步骤:构建模型和求解薛定谔方程。
首先,需要确定材料的晶胞结构,即原子的排列方式和间距。
其次,需要选择合适的计算方法,如密度泛函理论(DFT)等。
DFT是一种基于电子密度的近似方法,它将材料中的电子相互作用简化为一个电子密度函数。
然后,需要选取计算所需的参数,包括平面波基组、能量截断和k点网格等。
最后,通过求解薛定谔方程,可以得到材料中电子的波函数和能量等信息。
第一性原理计算在材料科学中有广泛的应用。
首先,它可以用于材料的结构预测和优化。
通过计算不同原子和离子的结合能、晶格参数和局域构型能等信息,可以预测新材料的结构和稳定性,为材料设计和合成提供指导。
其次,第一性原理计算可以用于研究材料的电子性质。
通过计算材料的能带结构、禁带宽度和电子态密度等信息,可以预测材料的导电性和光学性质。
此外,第一性原理计算还可以用于模拟材料的机械性质、热学性质和磁学性质等。
尽管第一性原理计算有广泛的应用,但其存在一些限制。
首先,求解薛定谔方程是一项复杂且计算量大的任务,需要高性能计算机和大量的计算时间。
其次,第一性原理计算通常采用一些近似方法,如DFT等,会带来一定的误差。
此外,由于计算的复杂性,第一性原理计算通常只能研究小尺寸的体系,难以模拟大尺寸和复杂的材料。
第一性原理计算方法讲义
第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。
而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。
第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。
量子力学是反映微观粒子运动规律的理论。
量子力学的出现,使得人们对于物质微观结构的认识日益深入。
原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。
量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。
以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。
目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。
但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。
绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。
Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。
但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。
1964年,Hohenberg和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT)。
它建立在非均匀作为基本变量。
1965年,Kohn和Sham 电子气理论基础之上,以粒子数密度()r提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。
针对固体可靠而准确的稳定性第一性原理的详细概述
针对固体可靠而准确的稳定性第一性原理的详细概述可靠而准确的稳定性第一性原理计算对材料的合成、反应性和性质的研究至关重要,对于探索新的化学空间和难以观察的相结构来说也很重要。
稳定性涉及化学稳定能力(不同化学环境中化学计量组成保持不变的能力)和结构选择(一定化学计量组成的晶体结构取向)两个方面,对于已知材料,可用实验测量生成焓来预测。
然而对于新材料和化学组成相同而空间结构不同的材料来说,则不可能或难以用实测生成焓来预测。
来自美国的华裔科学家Haowei Peng 和Jianwei Sun,应用密度泛函理论结合最近开发的强约束-适当规范(strongly constrained and appropriately normed, SCAN)泛函,对材料稳定问题的两个方面(化学稳定能力">Efficient first-principles prediction of solid stability: Towards chemical accuracyYubo Zhang, Daniil A. Kitchaev, Julia Yang, Tina Chen, Stephen T. Dacek, Rafael A. Sarmiento-Pérez, Maguel A. L. Marques, Haowei Peng, Gerbrand Ceder, John P. Perdew ">Abstract The question of material stability is of fundamental importance to any analysis of system properties in condensed matter physics and materials science. The ability to evaluate chemical stability, i.e., whether a stoichiometry will persist in some chemical environment, and structure selection, i.e. what crystal structure a stoichiometry will adopt, is critical to the prediction of materials synthesis, reactivity and properties. Here, we demonstrate that density functional theory, with the recently developed strongly constrained and appropriately normed (SCAN) functional, has advanced to a point where both facets of the stability problem can be reliably and efficiently predicted for main group compounds, while transition metal compounds are improved but remain a challenge. SCAN therefore offers a robust model for a significant portion of the periodic table, presenting an opportunity for the development of novel materials and the study of fine phase transformations even in largely unexplored systems with little to no experimental data.。
第一性原理计算Mgn(n≤11)团簇的基态结构及稳定性
第一性原理计算Mgn(n≤11)团簇的基态结构及稳定性摘要:本文应用第一性原理平面波模守恒赝势方法研究了中性和一价阴离子(n=2~11)团簇的几何结构,稳定性和电子性质。
结合能计算表明,中性团簇尺寸依赖的稳定性与电子壳层闭合和几何结构“闭合”密切相关。
关键词:镁团簇第一性原理几何结构1 引言近年来人们发现随着团簇尺寸的变化,团簇结合键的性质发生了变化。
含原子数目少的金属团簇往往不具备金属的一些特性,随着原子数目的增加,金属性才逐步显示出来。
镁团簇非常适合于研究随原子数目增加团簇金属性发生变化这一性质。
近年来许多文献用实验和理论模拟的方法对镁团簇进行了研究[1~5]。
但是对小团簇的最稳定结构仍然存在许多争议。
主要是因为,团簇对于光谱技术来说它显得太大;对于散射技术来说它又显得过小,因此目前还没有一种能直接观察在自由状态下小团簇结构的方法。
本文应用基于密度泛函理论的第一性原理方法对2~11个原子尺度的中性和负一价团簇的几何学结构及稳定性进行了研究。
2 计算方法采用基于密度泛函理论的CASTEP软件包对各种初始结构进行充分的几何优化。
用模守恒赝势描述电子-离子之间相互作用,用PBE 形式的广义梯度近似计算交换关联能,用BFGS方法优化模型的几何结构。
平面波截断能量为330eV;布里渊区积分在倒异空间Γ点。
当最大的能量公差、应力公差、最大原子位移公差的值分别小于5.0×10-6eV/Aatom,0.01eV/A,5.0×10-4时认为计算收敛。
计算温度为0K。
计算采用周期性边界条件,将一个原子团簇放入具有周期性的超晶胞内。
超晶胞的晶格常数为20A,这样团簇与团簇之间的距离足够大以消除相邻团簇之间的电子云相互影响所带来的误差。
(注:1A=10-10m=0.1nm)。
3 结果和讨论3.1 镁团簇的几何结构参考文献[1]Th.Diederich,T.Dppner,J. Braune,J.Tiggesbumker,K.-H. Meiwes-Broer,Phys.Rev.Lett.86(2001) 4807[2]O.C.Thomas,W.Zheng,S.Xu,K.H.Bowen,Phys.Rev.Lett.89(2002) 213403.[3]K.P.Huber,G.Herzberg,Molecular Spectra and Molecular Structure.IV.Constants of Diatomic Molecules(Van Nostrand Reinhold,New York,1979).[4]V.Kumar,R.Car,Phys.Rev.B 44 (1991)8243.[5]J.Akola,K.RytKonen,M.Manninen,Eur.Phys.J.D 16(2001)21.[6]A.Lyalin,I.A.Solov’yov,A.V.Solov’yov,W.Greiner,Phys.Rev.A 67(2003) 063203.[7]C.Kittel,Introduction to Solid State Physics,7th edn.(John Wiley Sons, New York,1996).[8]S.Srinivas,J.Jellinek,J.Chem.Phys.121(2004)7243.[9]M.Kabir,A.Mookerjee,D.G. Kanhere,Phys.Rev.B 73(2006)224439.[10]M.Kabir,D.G.Kanhere,A.Mookerjee,Phys.Rev.B75(2007)214433.。
材料第一性原理操作和试验方法
CASTAP性质任务
CASTAP性质任务允许在完成能量,几何优化或动力学运行之后求出电 子和结构性质。可以产生的性质如下:
* 态密度(DOS):利用原始模拟中产生的电荷密度和势能,非自恰计算价 带和导带的精细Monkhorst-Pack 网格上的电子本征值。
* 带结构:利用原始模拟中产生的电荷密度和势能,非自恰计算价带和导 带的布里渊区高对称性方向电子本征值。
CASTAP中能量的缺损单位是电子伏特(eV),各种能量单位的换 算关系见Mohr.P.J(2000).
1 eV=0.036749308 Ha=23.0605 kcal/mole=96.4853 kJ/mole
CASTAP几何优化任务
CASTAP几何优化任务允许改善结构的几何,获得稳定结构 或多晶型物。通过一个迭代过程来完成这项任务,迭代过程中 调整原子坐标和晶胞参数使结构的总能量最小化。
2 第一原理所能解决的问题
Mo2C 是合金钢中的 二次硬化相,它均匀、 细小、弥散的分布于 基体中,具有极高的 硬度和良好的高温稳 定性。
Hexagoal Mo2C
①合金钢中合金元素种类繁多,是否存在其他的碳化物相,如Cr2C、 V2C、Ni2C来取代其作用。 ② Mo2C 是否具有其他更加稳定的结构 ③ 空位或其他原子置换对其稳定性和性能有何影响
Metal Crystal
atomic Meaured Calculated
structure configuration (B/atom) (DFT) (B)
Fe
BCC
3d64s2
2.22
2.25
Co
HCP
3d74s2
1.72
1.63
Ni
CCP
第一性原理计算判断材料稳定性的几种方法
第一性原理计算判断材料稳定性的几种方法当我们通过一些方法,如:人工设计、机器学习和结构搜索等,设计出一种新材料的时候,首先需要做的一件事情就是去判断这个材料是否稳定。
如果这个材料不稳定,那么后续的性能分析就犹如空中楼阁。
因此,判断材料是否稳定是材料设计领域中非常关键的一个环节。
接下来,我们介绍几种通过第一性原理计算判断材料是否稳定的方法。
011.结合能结合能是指原子由自由状态形成化合物所释放的能量,一般默认算出来能量越低越稳定。
对于简单的二元化合物A m B n(A,B为该化合物中包含的两种元素,m,n为相应原子在化学式中的数目),其结合能可表示为:其中E(A m B n)为化学式A m B n的能量,E(A)和E(B)分别为自由原子A和B的能量,E b越低,越稳定。
022.形成能形成能是指由相应单质合成化合物所释放的能量。
同样,对于二元化合物A m B n,其形成能可表示为:其中E(A)和E(B)分别为对应单质A和B归一化后的能量。
用能量判断某一材料稳定性的时候,选择形成能可能更符合实际。
因为实验合成某一材料的时候,我们一般使用其组成单质进行合成。
如果想进一步判断该材料是处于稳态还是亚稳态,那么需要用凸包图(convex hull)进行。
如图1所示,计算已知稳态A x B y的形成能,构成凸包图(红色虚线),其横轴为B在化学式中所占比例,纵轴为形成能。
通过比较考察化合物与红色虚线的相对位置,如果在红色虚线上方则其可能分解(如:图1 插图中的D,将分解为A和B)或处于亚稳态(D的声子谱没有虚频);如果在红色虚线下方(如:图 1 插图中的C),则该化合物稳定。
图1:凸包图用于判断亚稳态和稳态[[1]]033.声子谱声子谱是表示组成材料原子的集体振动模式。
如果材料的原胞包含n个原子,那么声子谱总共有3n支,其中有3条声学支,3n-3条光学支。
声学支表示原胞的整体振动,光学支表示原胞内原子间的相对振动。
高考化学物质稳定性判断四招
高考化学物质稳定性判断四招(总1
页)
--本页仅作为文档封面,使用时请直接删除即可--
--内页可以根据需求调整合适字体及大小--
高中阶段的话,只需要掌握这几点就可以了.
1.单质,考察原子半径,半径小,自然结合力就大,稳定性就高.
2.氢氧化物,金属性越强,碱的热稳定性越强(碱性越强,热稳定性越强).?
3.含氧酸,高中好像不会考的很细致,因为涉及到大学的反极化作用.只要记住几个常见的,比如硝酸不稳定,硫酸很稳定等等就可以了.
4.气态氢化物,元素的非金属性越强,形成的气态氢化物就越稳定.同主族的非金属元素,从上到下,随核电荷数的增加,非金属性渐弱,气态氢化物的稳定性渐弱;同周期的非金属元素,从左到右,随核电荷数的增加,非金属性渐强,气态氢化物的稳定性渐强.
2。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第一性原理计算判断材料稳定性的几种方法
当我们通过一些方法,如:人工设计、机器学习和结构搜索等,设计出一种新材料的时候,首先需要做的一件事情就是去判断这个材料是否稳定。
如果这个材料不稳定,那么后续的性能分析就犹如空中楼阁。
因此,判断材料是否稳定是材料设计领域中非常关键的一个环节。
接下来,我们介绍几种通过第一性原理计算判断材料是否稳定的方法。
1.结合能
结合能是指原子由自由状态形成化合物所释放的能量,一般默认算出来能量越低越稳定。
对于简单的二元化合物A m B n(A,B为该化合物中包含的两种元素,m,n为相应原子在化学式中的数目),其结合能可表示为:
其中E(A m B n)为化学式A m B n的能量,E(A)和E(B)分别为自由原子A和B的能量,E b越低,越稳定。
2.形成能
形成能是指由相应单质合成化合物所释放的能量。
同样,对于二元化合物A m B n,其形成能可表示为:
其中E(A)和E(B)分别为对应单质A和B归一化后的能量。
用能量判断某一材料稳定性的时候,选择形成能可能更符合实际。
因为实验合成某一材料的时候,我们一般使用其组成单质进行合成。
如果想进一步判断该材料是处于稳态还是亚稳态,那
么需要用凸包图(convex hull)进行。
如图1所示,计算已知稳态A x B y的形成能,构成凸包图(红色虚线),其横轴为B在化学式中所占比例,纵轴为形成能。
通过比较考察化合物与红色虚线的相对位置,如果在红色虚线上方则其可能分解(如:图1 插图中的D,将分解为A和B)或处于亚稳态(D的声子谱没有虚频);如果在红色虚线下方(如:图1 插图中的C),则该化合物稳定。
图 1:凸包图用于判断亚稳态和稳态[[1]]
3.声子谱
声子谱是表示组成材料原子的集体振动模式。
如果材料的原胞包含n个原子,那么声子谱总共有3n支,其中有3条声学支,3n-3条光学支。
声学支表示原胞的整体振动,光学支表示原胞内原子间的相对振动。
计算出的声子谱有虚频,往往表示该材料不稳定。
因为
其中ω为振动频率,β可理解为弹性常量,E(x)表示原子间相互作用能,x表示原子偏离平衡位置的位移,m为原子质量。
由上式可以看出,当ω为虚频时,
,也就是表示原子平衡位置位于能量的“山顶”(类似抛物线顶点)。
很明显,处于该平衡位置的原子是不稳定的。
图2 单层2H-NbSe2的声子谱[[2]]
有些情况下,我们可以利用虚频信息使不稳定的材料变得稳定。
如图2所示,单层2H-NbSe2声子谱的一条声学支存在虚频,主要位于Γ点和M点1/2处(对应倒格矢的1/4位置)。
倒格矢的1/4,对应晶格长度的4倍。
我们可能需要将原胞沿上述倒格矢方向扩大四倍,进一步优化原子位置,才可能得到比较稳定的晶胞。
4.分子动力学和吉布斯自由能
通过能量和声子谱判断材料比较稳定之后,便可通过分析动力学或吉布斯自由能来进一步判断材料在一定温度下的稳定性。
分子动力学方法:首先构建超胞,然后施加一定温度,运行一段时间之后观察原胞结构是否遭到破坏来判断该材料能否在该温度下稳定存在。
吉布斯自由能可以用来比较不同构型材料在不同温度下的稳定性,如图3所示。
图3 几种碳的同素异形体在不同温度下的吉布斯自由能[[3]]
5.波恩稳定性判据
材料的弹性势能可以表示为,
其中V0为材料晶胞不受外力时的体积,C ij为弹性常量矩阵元,εi为应力。
如果一个材料的是稳定的,得到的弹性能E一定大于0。
这样就可以得到材料的弹性稳定性条件:矩阵C是正定的;矩阵C的所有本征值是正的;矩阵C的所有顺序主子式是正的;矩阵C的任意子式都是正的。
因此,不同晶系材料的弹性常量矩阵元需要满足不同的条件,具体可查看文献“Necessary and s ufficient elastic stability conditions in various crystal systems”。