======================================== Elastic Compliance Constants Sij (1/GPa) ========================================
0.0014282 -0.0002075 -0.0002075 -0.0002075 0.0014282 -0.0002075 -0.0002075 -0.0002075 0.0014282 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0023533 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0023533 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0023533
选择Job Control标签，设定本地机运算。 按下CASTEP Calculation对话框中的Run按钮。
优化之后，此结构的晶胞参数应为a=b=c=2.574Å。在我们可以继 续计算优化结构的弹性常数。
或按右键显示
2. 计算BN的弹性常数
BN CASTEP GeomOpt/BN.xsd处 于激活状态。选择CASTEP Calculation对话框中的Setup标签， 从Task的下拉清单中选择Elastic Constants。 按下More...按钮，CASTEP Elastic Constants对话框见右图。 将Number of steps for each strain由4增加为6，按Run运行。 CASTEP 的弹性常数计算任务的 结果以一批 .castep 输出文件的形式 给出。这些文件中的每一个文件都 代表确定的晶胞在假设的应变模式 和应变振幅下的几何优化运行结果。 这 些 文 件 的 命 名 约 定 为 ： seedname_cij__m__n 。 对 于 给 定 的 模式来说， m代表当前的应变模式， n代表当前的应变振幅。
在进一步的分析中没有使用压力的切点值 , 它很简单的指示出收敛 的基态离最初的结构有多远。 所有应变模式的结果总结如下： ============================ Summary of elastic constants ============================ id i j Cij (GPa) 1 1 1 788.92024 +/- 0.946 4 4 4 447.55108 +/- 0.749 7 1 2 148.70983 +/0.754 The errors are only provided when more than two values for the strain amplitude were used, since there is no statistical uncertainty associated with fitting a straight line to only two points. 弹性常数以常规的6x6张量的形式显示出,随后弹性柔量 （compliances）以相似的6x6形式显示出：
现在设置几何优化 从工具栏中选择CASTEP工具 ， 然后从下拉列中选择Calculation（或从 菜单栏中选择Modules / CASTEP / Calculation）。 CASTEP Calculation对话框见右图：
在Setup标签中，把Task设置为 Geometry Optimization，把Quality 设 置为Fine，并且把Functional设置为GGA and PW91。
预测锗的热力学属性
的对应。 和弹性系数相对应的应力 (压缩符)： 1 7 7 4 0 0 as induced by the strain components: 1 1 1 4 0 0 在下面的表格中给出了每一种应力组成的应力-应变线性适配关系: Stress Cij value of value of index index stress strain 1 1 -4.990578 0.003000 1 1 -5.949042 0.001800 1 1 -6.891618 -0.000600 1 1 -7.838597 0.000600 1 1 -8.784959 0.001800 1 1 -9.726562 C (gradient) : 788.920238 0.003000 Error on C : Correlation 0.945626 coeff: Stress intercept : 0.999997 7.363559 此梯度提供了弹性常数的数值（或弹性常数的线性组合），适配的 质量，由相关系数表示，提供了另人满意的弹性常数的不确定度。
3 弹性常数文件的描述 对于这种点阵类型，需要考虑两种应变模式(本教程只计算了一种)。对 于每一种应变模式，都有一个计算出的应力的总结(由各自的.castep文件得 到) 。
=============================================== Elastic constants from Materials Studio: CASTEP =============================================== Summary of the calculated stresses ********************************** Strain pattern: 1 ====================== Current amplitude: 1 Transformed stress tensor (GPa) : -4.990578 0.000000 0.000000 0.000000 -6.907159 0.000000 0.953658 0.953658 Current amplitude: 2 6.908215 Transformed stress tensor (GPa) : -5.949042 0.000000 0.000000 -7.093625 0.000000 0.000000 0.571307 0.571307 提供了应力，应变的组成和弹性常数张量之间联系的所有信息。在 7.094263 这一阶段，每一个弹性常数均有一个简洁的指数代表而不是由一对 ij指数代表。稍后会在文件夹中给出压缩符和常规的指数标定之间
文件的最后部分包含推出的属性:
Bulk modulus = 362.11330 +/- 0.593 (GPa) Compressibility = 0.00276 (1/GPa) Axis Young Modulus Poisson Ratios (GPa) X 741.74894 Exy= 0.1586 Exz= 0.1586 Y 741.74894 Eyx= 0.1586 Eyz= 0.1586 Z 741.74894 Ezx= 0.1586 Ezy= 0.1586 Lame constants for isotropic material (GPa) Lambda = -106.1819, Mu = 447.5511 END
按下more按钮，选中Optimize cell。 关闭CASTEP Geometry Optimization对话 框。
选 择 Electronic 标 签 ， 按 下 More... 按 钮 以 得 到 CASTEP Electronic Options 对 话 框 。 把 Derived grid 的 设 置 从 Standard 改 为 Fine 。 关 闭 CASTEP Electronic Options对话框。
===================================== Elastic Stiffness Constants Cij (GPa) =====================================
736.57379 125.20883 125.20883 0.00000 0.00000 0.00000 125.20883 736.57379 125.20883 0.00000 0.00000 0.00000 125.20883 0.00000 0.00000 125.20883 0.00000 0.00000 736.57379 0.00000 0.00000 0.00000 424.93974 0.00000 0.00000 0.00000 424.93974 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 424.93974
6
仅取一种应 变模式
CASTEP可以使用这些结果来分析每一个运行计算出来的压力张量，产生 一个有关弹性性质的文件。 从工具栏中选择 CASTEP 工具，然后选择 Analysis 或者从菜单栏 中选择Modules | CASTEP | Analysis。 从属性清单中选择Elastic constants，从BN 的弹性常数计算工作中得到的结果文件BN.castep应 自动显示在Results file选框中。按下Calculate按 钮。计算结束后产生一个新的文档BN Elastic Constants.txt。 此文档中的信息包括： *输入的应变和计算出的应力的总结 *每一种应变模式线性拟合和拟合质量的计算结果 *给定对称性下计算出的应力与弹性常数之间的对应 *弹性常数Cij和弹性柔量Sij的表格 *导出量：体积模量和其倒数、压缩系数、杨氏模量、 Poisson比、 Lame 常数(用于模拟各向同性介质)
首先导入BN结构 在 菜 单 栏 中 选 择 File/ Import ， 从 structures/semiconductors 中 选 中 BN.msi ， 按 Import 按 钮 ， 输入BN的晶体结构，见右图。 为了节省计算时间，由 Build / Symmetry / Primitive Cell将此 conventional representation 转化为primitive representation.