• 815.50 KB
  • 2022-04-29 14:48:23 发布

最新中科大-Materials-Studio-培训教程-16(包你学会!)请将这一系列全看完-一定有收获。概要课件PPT.ppt

  • 35页
  • 当前文档由用户上传发布,收益归属用户
  1. 1、本文档共5页,可阅读全部内容。
  2. 2、本文档内容版权归属内容提供方,所产生的收益全部归内容提供方所有。如果您对本文有版权争议,可选择认领,认领后既往收益都归您。
  3. 3、本文档由用户上传,本站不保证质量和数量令人满意,可能有诸多瑕疵,付费之前,请仔细先通过免费阅读内容等途径辨别内容交易风险。如存在严重挂羊头卖狗肉之情形,可联系本站下载客服投诉处理。
  4. 文档侵权举报电话:19940600175。
'进入夏天,少不了一个热字当头,电扇空调陆续登场,每逢此时,总会想起那一把蒲扇。蒲扇,是记忆中的农村,夏季经常用的一件物品。  记忆中的故乡,每逢进入夏天,集市上最常见的便是蒲扇、凉席,不论男女老少,个个手持一把,忽闪忽闪个不停,嘴里叨叨着“怎么这么热”,于是三五成群,聚在大树下,或站着,或随即坐在石头上,手持那把扇子,边唠嗑边乘凉。孩子们却在周围跑跑跳跳,热得满头大汗,不时听到“强子,别跑了,快来我给你扇扇”。孩子们才不听这一套,跑个没完,直到累气喘吁吁,这才一跑一踮地围过了,这时母亲总是,好似生气的样子,边扇边训,“你看热的,跑什么?”此时这把蒲扇,是那么凉快,那么的温馨幸福,有母亲的味道!  蒲扇是中国传统工艺品,在我国已有三千年多年的历史。取材于棕榈树,制作简单,方便携带,且蒲扇的表面光滑,因而,古人常会在上面作画。古有棕扇、葵扇、蒲扇、蕉扇诸名,实即今日的蒲扇,江浙称之为芭蕉扇。六七十年代,人们最常用的就是这种,似圆非圆,轻巧又便宜的蒲扇。  蒲扇流传至今,我的记忆中,它跨越了半个世纪,也走过了我们的半个人生的轨迹,携带着特有的念想,一年年,一天天,流向长长的时间隧道,袅中科大-Materials-Studio-培训教程-16(包你学会!)请将这一系列全看完-一定有收获。概要 1.优化BN立方晶体的结构在计算弹性常数之前并不一定要进行几何优化,可以由实验观测到的结构计算出Cij数据。尽管如此,如果我们完成晶胞的几何优化,可以获得更多相容的结果,进而计算与理论基态对应的弹性常数。弹性常数的精确度,尤其是切变常数的精确度,主要取决于SCF计算的品质,特别是布里渊区取样和波函数收敛程度的品质。所以我们设置SCF、k点取样和FFT格子的精度为Fine。首先导入BN结构在菜单栏中选择File/Import,从structures/semiconductors中选中BN.msi,按Import按钮,输入BN的晶体结构,见右图。为了节省计算时间,由Build/Symmetry/PrimitiveCell将此conventionalrepresentation转化为primitiverepresentation. 现在设置几何优化从工具栏中选择CASTEP工具,然后从下拉列中选择Calculation(或从菜单栏中选择Modules/CASTEP/Calculation)。CASTEPCalculation对话框见右图:在Setup标签中,把Task设置为GeometryOptimization,把Quality设置为Fine,并且把Functional设置为GGAandPW91。按下more按钮,选中Optimizecell。关闭CASTEPGeometryOptimization对话框。 从属性清单中选择Elasticconstants,从BN的弹性常数计算工作中得到的结果文件BN.castep应自动显示在Resultsfile选框中。按下Calculate按钮。计算结束后产生一个新的文档BNElasticConstants.txt。此文档中的信息包括:*输入的应变和计算出的应力的总结*每一种应变模式线性拟合和拟合质量的计算结果*给定对称性下计算出的应力与弹性常数之间的对应*弹性常数Cij和弹性柔量Sij的表格*导出量:体积模量和其倒数、压缩系数、杨氏模量、Poisson比、Lame常数(用于模拟各向同性介质)CASTEP可以使用这些结果来分析每一个运行计算出来的压力张量,产生一个有关弹性性质的文件。从工具栏中选择CASTEP工具,然后选择Analysis或者从菜单栏中选择Modules|CASTEP|Analysis。 3弹性常数文件的描述对于这种点阵类型,需要考虑两种应变模式(本教程只计算了一种)。对于每一种应变模式,都有一个计算出的应力的总结(由各自的.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.93974 Bulkmodulus=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点在倒易空间中的位置。高对称性点G,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|Analysis。从属性列表中选择Thermodynamicproperties。确定Resultsfile选择框中显示的是Ge_PhononDOS.castep。勾选上Debyetemperature图,按下View按钮。在结果文件夹中创建了两个新的图形文档GeThermodynamicProperties.xcd和GeDebyeTemperature.xcd。所示图形如下: 没有非谐性的实验结果表明在高温极限的Debye温度是395(3)K。模拟的Debye温度是396K,很好的与实验值相符。考虑到所完成的计算的级别(giventhelevelofcalculationperformed),这个结果还是可以接受的。使用更好的k点格子可以提高精确度。通过实验,推导出Debye温度在低温极限时的数值为374K,而CASTEP的预测值为460K。在低温极限时的错误可通过下面的事实来解释,在G点的频率并不是严格为0。我们也应该注意到类似的谐波近似是声子计算的基础,但是当温度高于Debye温度的三分之一时,类似的谐波近似是无效的。结果,在温度高于Debye温度时,非谐性影响出现,不能够仅仅依靠量的改变来解释。尽管如此,总的来说,实验图表和由CASTEP产生的图表在质量上是非常相似的。当Debye最低温度达到255K时,在实验图表上会骤然降低25K;CASTEP显示出温度骤降出现在相同的位置并且预测出Debye温度的最小值大约为270K。END 怎样滚得远 教学目标1、让学生参与探索斜面与地面成怎样的角度时物体滚得最远的实验活动,进行收集数据、求平均数、角的测量等数学知识和方法的综合运用,体会数学的应用价值。2、让学生感受做实验是研究事物的一种方法,培养学生实事求是的科学态度,进一步激发学生应用数学知识解决实际问题的兴趣。 教学准备:平滑的木板,圆柱形物体,卷尺,量角器,实验报告单。 (一)      创设情境,揭示课题这里利用木板斜坡和山的斜坡来运东西,即省钱又省力。 你能举一些像上面的例子吗?你知道斜坡和地面成什么角度时,物体滚动得远一些吗? 实验⑴搭一搭用30-50厘米的木板搭出一个与地面成30度的斜坡 .⑵滚一滚圆形物体在轻轻放在斜坡顶上,让它自动往下滑落.⑶量一量从木板底部量出物体在地面上滚动的长度. 填写下表:第一次第二次第三次平均数物体滚动的长度是不是每次滚动的长度都一样长? 学生每8人一组,按刚才师生合作的方法和步骤进行,边实验边填写实验报告单上的表格。木板与地面的角度:45度60度第一次第二次第三次平均数物体滚动的长度 阅读“你知道吗”通过阅读你知道了什么? 这节课你学到了什么?'