中科大-Materials-Studio-弹性常数计算.ppt
- 文档编号:18866584
- 上传时间:2024-02-03
- 格式:PPT
- 页数:25
- 大小:759KB
中科大-Materials-Studio-弹性常数计算.ppt
《中科大-Materials-Studio-弹性常数计算.ppt》由会员分享,可在线阅读,更多相关《中科大-Materials-Studio-弹性常数计算.ppt(25页珍藏版)》请在冰点文库上搜索。
背景背景:
当前,可应用于大周期性体系的密度泛函理论(当前,可应用于大周期性体系的密度泛函理论(DFTDFT)取得了显著)取得了显著的进展,已经成为解决材料设计、加工中难题的有效方法。
人们依据这个的进展,已经成为解决材料设计、加工中难题的有效方法。
人们依据这个理论可以使解释实验数据,预测新晶体的结构、结合能和表面活性等基本理论可以使解释实验数据,预测新晶体的结构、结合能和表面活性等基本性质。
这些工具可以用来指导设计新材料,允许研究人员理解基本的化学性质。
这些工具可以用来指导设计新材料,允许研究人员理解基本的化学和物理过程。
和物理过程。
绪论绪论:
在本教程中,将学习如何使用在本教程中,将学习如何使用CASTEPCASTEP来计算弹性常数和其他的力来计算弹性常数和其他的力学性能。
首先我们要优化学性能。
首先我们要优化BNBN立方晶体的结构,然后计算它的弹性常数。
立方晶体的结构,然后计算它的弹性常数。
本指南主要包括以下内容:
本指南主要包括以下内容:
11优化优化BNBN立方晶体的结构立方晶体的结构22计算计算BNBN的弹性常数的弹性常数33弹性常数文件的描述弹性常数文件的描述使用使用CASTEPCASTEP计算计算BNBN的弹性常数的弹性常数目的:
使用目的:
使用CASTEPCASTEP计算弹性常数计算弹性常数模块:
模块:
MaterialsVisualizer,CASTEPMaterialsVisualizer,CASTEP前提:
已使用前提:
已使用firstprinciplesfirstprinciples预测了预测了AlAsAlAs的晶格常数的晶格常数1.1.优化优化BNBN立方晶体的结构立方晶体的结构在计算弹性常数之前并不一定要进行几何优化,可以由实验观测到的在计算弹性常数之前并不一定要进行几何优化,可以由实验观测到的结构计算出结构计算出CCijij数据。
尽管如此,如果我们完成晶胞的几何优化,可以获得数据。
尽管如此,如果我们完成晶胞的几何优化,可以获得更多相容的结果,进而计算与理论基态对应的弹性常数。
更多相容的结果,进而计算与理论基态对应的弹性常数。
弹性常数的精确度,尤其是切变常数的精确度,主要取决于弹性常数的精确度,尤其是切变常数的精确度,主要取决于SCFSCF计算计算的品质,特别是布里渊区取样和波函数收敛程度的品质。
所以我们设置的品质,特别是布里渊区取样和波函数收敛程度的品质。
所以我们设置SCFSCF、kk点取样和点取样和FFTFFT格子的精度为格子的精度为FineFine。
首先导入首先导入BNBN结构结构在菜单栏中选择在菜单栏中选择File/File/ImportImport,从,从structures/semiconductorsstructures/semiconductors中中选中选中BN.msiBN.msi,按ImportImport按钮,输入BNBN的晶体结构,见右图的晶体结构,见右图。
为了节省计算时间,由为了节省计算时间,由Build/Symmetry/Build/Symmetry/PrimitiveCellPrimitiveCell将此将此conventionalrepresentationconventionalrepresentation转化为转化为primitiveprimitiverepresentation.representation.现在设置几何优化现在设置几何优化从工具栏中选择从工具栏中选择CASTEPCASTEP工具工具,然后从下拉列中选择,然后从下拉列中选择CalculationCalculation(或(或从菜单栏中选择从菜单栏中选择Modules/CASTEP/Modules/CASTEP/CalculationCalculation)。
)。
CASTEPCalculationCASTEPCalculation对话框见右图:
对话框见右图:
在在SetupSetup标签中,把标签中,把TaskTask设置为设置为GeometryOptimizationGeometryOptimization,把,把QualityQuality设置为设置为FineFine,并且把,并且把FunctionalFunctional设置为设置为GGAandPW91GGAandPW91。
按下按下moremore按钮,选中按钮,选中OptimizecellOptimizecell。
关闭关闭CASTEPGeometryOptimizationCASTEPGeometryOptimization对对话框。
话框。
选择选择ElectronicElectronic标签,按下标签,按下More.More.按钮以得到按钮以得到CASTEPElectronicCASTEPElectronicOptionsOptions对话框。
把对话框。
把DerivedgridDerivedgrid的设置从的设置从StandardStandard改为改为FineFine。
关闭。
关闭CASTEPElectronicOptionsCASTEPElectronicOptions对话框。
对话框。
选择选择JobControlJobControl标签,设定本地机运算标签,设定本地机运算。
按下按下CASTEPCalculationCASTEPCalculation对话框中的对话框中的RunRun按钮按钮。
优化之后,此结构的晶胞参数应为优化之后,此结构的晶胞参数应为a=b=c=2.574a=b=c=2.574。
现在我们可以。
现在我们可以继续计算优化结构的弹性常数。
继续计算优化结构的弹性常数。
或按右键显示或按右键显示2.2.计算计算BNBN的弹性常数的弹性常数BNCASTEPGeomOpt/BN.xsdBNCASTEPGeomOpt/BN.xsd处处于激活状态。
于激活状态。
选择选择CASTEPCASTEPCalculationCalculation对话框中的对话框中的SetupSetup标签标签,从,从TaskTask的下拉清单中选择的下拉清单中选择ElasticConstantsElasticConstants。
按下按下More.More.按钮,按钮,CASTEPCASTEPElasticConstantsElasticConstants对话框见右图。
对话框见右图。
将将NumberofstepsforeachNumberofstepsforeachstrainstrain由由44增加为增加为66,按,按RunRun运行运行。
CASTEPCASTEP的弹性常数计算任务的的弹性常数计算任务的结果以一批结果以一批.castep.castep输出文件的形输出文件的形式给出。
这些文件中的每一个文件式给出。
这些文件中的每一个文件都代表确定的晶胞在假设的应变模都代表确定的晶胞在假设的应变模式和应变振幅下的几何优化运行结式和应变振幅下的几何优化运行结果。
这些文件的命名约定为:
果。
这些文件的命名约定为:
seedname_cij_m_nseedname_cij_m_n。
对于给定的。
对于给定的模式来说,模式来说,mm代表当前的应变模式代表当前的应变模式,nn代表当前的应变振幅。
代表当前的应变振幅。
6仅取一种应仅取一种应变模式变模式从属性清单中选择从属性清单中选择ElasticconstantsElasticconstants,从,从BNBN的弹性常数计算工作中得到的结果文件的弹性常数计算工作中得到的结果文件BN.castepBN.castep应自动显示在应自动显示在ResultsfileResultsfile选框中。
按选框中。
按下下CalculateCalculate按钮。
计算结束后产生一个新的文档按钮。
计算结束后产生一个新的文档BNElasticConstants.txtBNElasticConstants.txt。
此文档中的信息包括此文档中的信息包括:
*输入的应变和计算出的应力的总结输入的应变和计算出的应力的总结*每一种应变模式线性拟合和拟合质量的计算结果每一种应变模式线性拟合和拟合质量的计算结果*给定对称性下计算出的应力与弹性常数之间的对给定对称性下计算出的应力与弹性常数之间的对应应*弹性常数弹性常数CijCij和弹性柔量和弹性柔量SijSij的表格的表格*导出量:
体积模量和其倒数、压缩系数、杨氏模导出量:
体积模量和其倒数、压缩系数、杨氏模量、量、PoissonPoisson比、比、LameLame常数常数(用于模拟各向同性用于模拟各向同性介质介质)CASTEPCASTEP可以使用这些结果来分析每一个运行计算出来的压力张量,产可以使用这些结果来分析每一个运行计算出来的压力张量,产生一个有关弹性性质的文件。
生一个有关弹性性质的文件。
从工具栏中选择从工具栏中选择CASTEPCASTEP工具,然后选择工具,然后选择AnalysisAnalysis或者从菜单或者从菜单栏中选择栏中选择Modules|CASTEP|AnalysisModules|CASTEP|Analysis。
33弹性常数文件的描述弹性常数文件的描述对于这种点阵类型,需要考虑两种应变模式对于这种点阵类型,需要考虑两种应变模式(本教程只计算了一种本教程只计算了一种)。
对于每一种应变模式,都有一个计算出的应力的总结对于每一种应变模式,都有一个计算出的应力的总结(由各自的由各自的.castep.castep文文件得到件得到)。
=ElasticconstantsfromMaterialsStudio:
CASTEP=Summaryofthecalculatedstresses*Strainpattern:
1=Currentamplitude:
1Transformedstresstensor(GPa):
-4.9905780.0000000.0000000.000000-6.9071590.9536580.0000000.953658-6.908215Currentamplitude:
2Transformedstresstensor(GPa):
-5.9490420.0000000.0000000.000000-7.0936250.5713070.0000000.571307-7.094263提供了应力,应变的组成和弹性常数张量之间联系的所有信息。
在这一阶段,每一个弹性常数均有一个简洁的指数代表而不是由一对ij指数代表。
稍后会在文件夹中给出压缩符和常规的指数标定之间的对应。
和弹性系数相对应的应力(压缩符):
177400asinducedbythestraincomponents:
111400在下面的表格中给出了每一种应力组成的应力-应变线性适配关系:
StressCijvalueofvalueofindexindexstressstrain11-4.990578-0.00300011-5.949042-0.00180011-6.891618-0.00060011-7.8385970.00060011-8.7849590.00180011-9.7265620.003000C(gradient):
788.920238ErroronC:
0.945626Correlationcoeff:
0.999997Stressintercept:
-7.363559此梯度提供了弹性常数的数值(或弹性常数的线性组合),适配的质量,由相关系数表示,提供了另人满意的弹性常数的不确定度。
在进一步的分析中没有使用压力的切点值,它很简单的指示出收敛的基态离最初的结构有多远。
所有应变模式的结果总结如下:
=Summaryofelasticconstants=idijCij(GPa)111788.92024+/-0.946444447.55108+/-0.749712148.70983+/-0.754Theerrorsareonlyprovidedwhenmorethantwovaluesforthestrainamplitudewereused,sincethereisnostatisticaluncertaintyassociatedwithfittingastraightlinetoonlytwopoints.弹性常数以常规的6x6张量的形式显示出,随后弹性柔量(compliances)以相似的6x6形式显示出:
=ElasticStiffnessConstantsCij(GPa)=ElasticComplianceConstantsSij(1/GPa)=0.0014282-0.0002075-0.00020750.00000000.00000000.0000000-0.00020750.0014282-0.00020750.00000000.00000000.0000000-0.0002075-0.00020750.00142820.00000000.00000000.00000000.00000000.00000000.00000000.00235330.00000000.00000000.00000000.00000000.00000000.00000000.00235330.00000000.00000000.00000000.00000000.00000000.00000000.0023533文件的最后部分包含推出的属性:
736.57379125.20883125.208830.000000.000000.00000125.20883736.57379125.208830.000000.000000.00000125.20883125.20883736.573790.000000.000000.000000.000000.000000.00000424.939740.000000.000000.000000.000000.000000.00000424.939740.000000.000000.000000.000000.000000.00000424.93974Bulkmodulus=362.11330+/-0.593(GPa)Compressibility=0.00276(1/GPa)AxisYoungModulusPoissonRatios(GPa)X741.74894Exy=0.1586Exz=0.1586Y741.74894Eyx=0.1586Eyz=0.1586Z741.74894Ezx=0.1586Ezy=0.1586Lameconstantsforisotropicmaterial(GPa)Lambda=-106.1819,Mu=447.5511END预测锗的热力学属性背景线性响应或密度功能混乱理论是点阵动力学从头开始计算中最受欢迎的方法之一,尽管如此,这种方法的应用已经扩充到对振动属性的研究。
线性响应提供了一种分析方法用于计算给定混乱的二级派生的整体能量。
可以计算出许多属性,主要依赖于混乱的种类。
在离子位置的混乱可以引起动力矩阵和声子;在磁场中引起NMR效应;在单位晶格矢量中产生弹性常数;在电场中引起非传导性效应等。
在本指南中,我们将要学习为了计算声子散射和能态密度以及预测热力学属性如焓和自由能,如何使用CASTEP来完成线性响应计算。
本指南主要包含以下内容:
1优化锗单胞的结构2计算声子散射和能态密度3显示声子散射和能态密度4显示热力学属性1优化锗单胞的结构首先我们要导入锗的结构,它包含在MaterialsStudio所提供的结构库中。
在菜单栏中选择File|Import。
遵循下列路径structures/metals/puremetals选中Ge.xsd。
把它转换为原胞结构后,对它的计算会更快。
从菜单栏中选择Build|Symmetry|PrimitiveCell。
锗的原胞结构如右所示:
现在开始使用CASTEP来优化锗的几何结构。
从工具栏中选择CASTEP工具然后选中Calculation或从菜单栏中选择Modules|CASTEP|Calculation。
CASTEPCalculation的对话框如下:
几何优化的默认值不包括对单胞的优化。
在Setup标签上,把Task从Energy改为GeometryOptimization,把Functional改为LDA。
在CASTEPGeometryOptimization对话框中,按下More.按钮,勾选上OptimizeCell。
选中Electronic标签,把Energycutoff设置为Ultra-fine,把SCFtolerance设置为Ultra-fine,把k-pointset设置为Coarse以及把Pseudo-potentials设置为Norm-conserving。
选中JobControl标签。
选择你想要在其上运行工作的Gatewaylocation。
把Runtimeoptimization设置为Speed。
按下Run按钮开始运行。
工作递交后开始运行。
它大概需要2分钟时间,这主要取决于你的电脑的速度。
结果放在一个被称为GeCASTEPGeomOpt的新文件夹中。
2计算声子散射和能态密度为了计算声子散射和声子的能态密度,在从CASTEPCalculation对话框的Properties标签选定适当的属性后,我们必须完成一个单点能量计算。
确定GeCASTEPGeomOpt文件夹中的Ge.xsd文件示激活文档。
选中CASTEPCalculation对话框中的Setup标签,把Task设置为Energy。
选择Electronic标签,按下More按钮,显示出CASTEPElectronicOptions对话框。
选择SCF标签,勾选上Fixoccupancy。
关闭CASTEPElectronicOptions对话框。
在CASTEPCalculation对话框中选定Properties标签。
选择Phonondispersion并且把q-vectorset设置为medium。
选择Phonondensityofstates并且把q-vectorset设置为medium。
选中JobControl标签。
选择你想要在其上运行工作的Gatewaylocation。
按下Run按钮。
在GeCASTEPGeomOpt文件夹中创建了一个名为GeCASTEPEnergy的新文件。
当能量计算完成后,两个新文件Ge_PhononDisp.castep和Ge_PhononDOS.castep放在此文件夹中。
3显示声子散射和能态密度声子散射曲线显示出声子能量沿着布里渊区高对称性方向如何依赖于q向量。
此信息可以从单晶的中子散射实验中获得。
只有为数不多的物质可以获得这样的信息,所以用来确定建模方法是否正确的理论偏差曲线对于论证在卷首的从头开始计算方法的预测性能力是非常有用的。
在一定情形下,它可能测量态密度而不是声子散射。
而且与声子的态密度有直接关系的电子声子交感作用可以通过隧道实验直接测量。
所以能够从第一原理计算出声子的态密度是非常重要的。
MaterialsStudio可以从任何.phononCASTEP输出文件中产生声子散射图和态密度图。
这些文件隐藏在ProjectExplorer里,但是一个.phonon文件会和每一个拥有PhonDisp或PhonDOS后缀的.castep文件一起产生。
现在,我们使用先前的计算结果来创建声子散射图。
从MaterialsStudio的菜单栏中选择Modules|CASTEP|Analysis。
从属性列表中选择Phonondispersion。
确定Resultsfile选择框中显示的是Ge_PhononDisp.castep。
从Units下拉列表中选择cm-1。
从Graphstyle下拉列表中选择Line。
按下View按钮。
在结果文档中创建了一个新的图形文档GePhononDispersion.xcd。
它应和下图相似:
声子散射的实验图如下所示:
预测的频率可从Ge_PhononDisp.castep文件中得到。
在ProjectExplorer中双击Ge_PhononDisp.castep。
按下CRTL+F键,搜索VibrationalFrequencies。
结果文件中的部分内容如下所示:
=+VibrationalFrequencies+-+BranchnumberFrequency(cm-1)+=+-+q-pt=1(0.5000000.2500000.750000)0.076923+Effectivecut-off=252.3267EV+-+1142.467048+2154.267057+3226.736315+4228.208353+5285.654435+6292.950626+-+q-pt=2(0.5000000.3000000.700000)0.038462+Effectivecut-off=246.7136EV+-+每一个q点和每一个分支(纵向光波或声波(LO/LA),横向光波或声波(TO/TA)的频率以cm-1给出,同时也给出了q点在倒易空间中的位置。
高对称性点,L和X在倒易空间中的位置各自为(000),(0.50.50.5)和(0.500.5)。
这些点和q点12,6以及19相对应。
预测的频率(cm-1)和实验的频率(cm-1)如下:
总体来说,计算的精度是可以接受的。
在Gamma点错误的声学波频率,3cm-1而不是0cm-1,使我们感觉这只是一般的精确度。
通过对更好的SCFk点格子的计算,我们可以获得更加另人满意的实验结果。
现在创建声子态密度图:
从MaterialsStudio菜单栏中选择Modules|CASTEP|Analysis。
从属性列表中选择Phonondensityofstates。
确定Resultsfile选择框中显示的是Ge_PhononDOS.castep。
把DOSdisplay设置为Full。
按下View按钮。
创建了一个新的图形文档GeDOS.xcd。
它应当与下图相似:
4显示热力学属性在CASTEP中的声子计算可以用来评价近似准谐波晶体的焓,熵,自由能,格子的热容对于温度的依赖性。
可以用这些结果和实验数据(如,热容的测量)相比较以预测不同的结构经过修正后的相稳定性和相转变。
所有与能量相关的属性均画在同一种曲线图中,并且0点能量的计算值也包括在内。
热容被独自画在图表的右侧。
现在使用声子计算的结果创建热力学属性图表。
从MaterialsStudio菜单栏中选择Modules|CASTEP|An
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 中科大 Materials Studio 弹性 常数 计算