《第一性原理计算判断材料稳定性的几种方法》由会员分享，可在线阅读，更多相关《第一性原理计算判断材料稳定性的几种方法（4页珍藏版）》请在金锄头文库上搜索。
1、第一性原理计算判断材料稳定性的几种方法当我们通过一些方法，如：人工设计、机器学习和结构搜索等,设计出一种新材料的 时候，首先需要做的一件事情就是去判断这个材料是否稳走。如果这个材料不稳 定，那么后续的性能分析就犹如空中楼阁。因此，判断材料是否稳走是材料设计领 域中非常关键的一个坏节。接下来，我们介绍几种通过第一性原理计算判断材料是 否稳定的方法。011结合能结合能是指原子由自由状态形成化合物所释放的能量，一般默认算出来能量越低越稳 定。对于简单的二元化合物AmBn( A.B为该化合物中包含的两种元素，为相应原子在 化学式中的数目)，其结合能可表示为：其中E(AmBn)为化学式AmBn的能量,E
2、(A)和E(B)分别为自由原子A和B的能量，Eb 越低，越稳定。022.形成能形成能是指由相应单质合成化合物所弟睾放的能量。同样，对于二元化合物AmBn ,其形成能可表示为: 其中E （A）和E （B）分别为对应单质A和B归一化后的能量。用能量判断某一材料稳定性的时候，选择形成能可能更符合实际。因为实验合成某一 材料的时候，我们一般使用其组成单质进行合成。如果想进一步判断该材料是处于稳 态还是亚稳态，那么需要用凸包图（convex hull）进行。如图1所示，计算已知稳态 AxBy的形成能，构成凸包图（红色虚线），其横轴为B在化学式中所占比例，纵轴 为形成能。通过比较考察化合物与红色虚线的相对
3、位置，如果在红色虚线上方则其可 能分解（如：图1插图中的D，将分解为A和B）或处于亚稳态（D的声子谱没有虚 频）；如果在红色虚线下方（如：图1插图中的C），则该化合物稳定。图1 :凸包图用于判断亚稳态和稳态03M声子谱声子谱是表示组成材料原子的集体振动模式。如果材料的原胞包含n个原子，那么 声子谱总共有3n支，其中有3条声学支，3n-3条光学支。声学支表示原胞的整体 振动，光学支表示原胞内原子间的相对振动。计算出的声子谱有虚频，往往表示该材料不稳定。因为其中U）为振动频率，P可理解为弹性常量，E （x）表示原子间相互作用能，x表示原子偏离平衡位置的位移，m为原子质量。由上式可以看出，当co为虚
4、频0a?，也就是表示原子平衡位置位于能量的山顶（类似抛物线顶点）。很明显，处于该平衡位置的原子是不稳定的。图2单层2H-NbSe2的声子谱2有些情况下，我们可以利用虚频信息使不稳走的材料变得稳定。如图2所示，单层 2H-NbSe2声子谱的一条声学支存在虚频，主要位于点和M点1/2处（对应倒格矢 的1/4位置）。倒格矢的1/4,对应晶格长度的4倍。我们可能需要将原胞沿上述倒 格矢方向扩大四倍，进一步优化原子位置，才可能得到比较稳定的晶胞。044.分子动力学和吉布斯自由能通过能量和声子谱判断材料比较稳定之后，便可通过分析动力学或吉布斯自由能来 进一步判断材料在一定温度下的稳定性。分子动力学方法：首
5、先构建超胞，然后施 加一定温度，运行一段时间之后观察原胞结构是否遭到破坏来判断该材料能否在该 温度下稳定存在。吉布斯自由能可以用来比较不同构型材料在不同温度下的稳定 性，如图3所示。图3几种碳的同素异形体在不同温度下的吉布斯自由能3055.波恩稳定性判据材料的弹性势能可以表示为,E二瓦+飙乞q苗+提勺.L沪其中Vo为材料晶胞不受外力时的体积，5为弹性常量矩阵元，Ei为应力。如果一 个材料的是稳定的，得到的弹性能E一定大于0。这样就可以得到材料的弹性稳定 性条件：矩阵C是正定的；矩阵C的所有本征值是正的；矩阵C的所有顺序主子 式是正的；矩阵C的任意子式都是正的。因此，不同晶系材料的弹性常量矩阵元需 要满足不同的条件，具体可查看文献“Necessary and sufficient elastic stabilityconditions in various crystal systems,4o
