-
含能材料的弹性性质包括弹性常数、体弹模量、杨氏模量和剪切模量等。它反映材料在常温、静荷载作用下的宏观力学性能,不但决定材料对施加应力的响应方式,还能反映含能材料的稳定性等[1]。目前,含能材料的弹性性质得到了人们的广泛关注,并对RDX[2]、TATB[3]、HMX[4]、FOX-7[5]、CL-20[6]等含能材料开展了研究。5, 5′-双四唑肼(5, 5′-Hydrazinebistetrazole,HBT)晶体是一种著名的单质富氮、高能、钝感含能材料[7]。HBT含有83.7%的氮元素,爆炸燃烧后释放的气体主要为氮气、二氧化碳和水,因此被称为“绿色”含能材料[8-11]。本研究采用基于第一性原理的准谐振近似方法,计算固态HBT晶体在0~120 GPa高压下的弹性性质,同时分析其弹性的各向异性特征。
-
采用基于第一性原理的Quantum ESPRESSO软件包中的PWSCF模块[12]计算电子结构的总能量。广义梯度近似(GGA)下交换相关势采用Perdew-Burke-Ernzerhof (PBE)形式[13-14],原子赝势平面波基组采用超软赝势加核修正的C.pbe_v1.2.uspp.F.UPF,N.pbe_v1.2.uspp.F.UPF,H.pbe_v1.2.uspp.F.UPF形式[15],几何优化采用BFGS算法[16]。为增加计算精度,平面波基函数的截断能取为544 eV,布里渊空间采用Monkhorst-Pack方法[17],k点方案为6 × 6 × 3,总能自洽,离子位移最大值和应力最大值都能达到设置的收敛精度。在高压下几何结构优化的基础上,利用准谐振近似方法计算HBT晶体的弹性性质,详细的计算过程在软件Thermo_pw中体现[18]。
-
固态物质弹性理论计算中引入应变εj(j = 1, 2, ···, 6),晶体的一个应力用σi表示。当应力较小时,σi和εj满足[19]
${\sigma _i} = \sum\limits_{i,j = 1}^6 {{C_{ij}}} {\varepsilon _j}$ ${\varepsilon _i} = \sum\limits_{i,j = 1}^6 {{S_{ij}}} {\sigma _j}$ 式中:系数Cij为弹性常数(或者弹性劲度),矩阵S为矩阵C的逆矩阵,Sij为弹性柔顺系数。初始构型采用Klapötke等[7]获得的HBT晶体的实验结构,它属于单斜结构的分子晶体,点群为C2/c空间群,其原胞晶系属于三斜晶系,对称性最低,包含最多的矩阵元。弹性常数矩阵C是一个对称矩阵,共有21个独立矩阵元。计算中采用原胞结构,弹性常数矩阵表达式为
${{ C}} = \left( {\begin{array}{*{20}{c}} {{C_{11}}}&{{C_{12}}}&{{C_{13}}}&{{C_{14}}}&{{C_{15}}}&{{C_{16}}} \\ {{C_{12}}}&{{C_{22}}}&{{C_{23}}}&{{C_{24}}}&{{C_{25}}}&{{C_{26}}} \\ {{C_{13}}}&{{C_{23}}}&{{C_{33}}}&{{C_{34}}}&{{C_{35}}}&{{C_{36}}} \\ {{C_{14}}}&{{C_{24}}}&{{C_{34}}}&{{C_{44}}}&{{C_{45}}}&{{C_{46}}} \\ {{C_{15}}}&{{C_{25}}}&{{C_{35}}}&{{C_{45}}}&{{C_{55}}}&{{C_{56}}} \\ {{C_{16}}}&{{C_{26}}}&{{C_{36}}}&{{C_{46}}}&{{C_{56}}}&{{C_{66}}} \end{array}} \right)$ 考虑到二阶弹性常数与能量满足如下关系[20]
${C_{ij}} = \frac{1}{{{V_0}}}\left( {\frac{{{\partial ^2}E_{\rm c}}}{{\partial {\varepsilon _i}\partial {\varepsilon _j}}}} \right)$ 因此,只要知道稳定构型下的晶体体积V0和能量Ec,就可以获得弹性常数。
-
基于密度泛函理论的第一性原理,得到零温零压下HBT晶体的晶格参数,如表1所示,晶体几何结构如图1所示。从表1可以看出,晶格参数b、c和V的计算值均高于实验结果,而a和β却低于实验结果。在目前的第一性原理方法中,晶格参数的计算精度取决于原子赝势和计算方法的选择,因此,无法做到与实验结果完全符合。此外,采用5 GPa间隔优化HBT晶体在0~120 GPa压强范围内的晶体结构。优化后获得了稳定的晶体构型,压强-体积比(p-V/V0)曲线与实验结果的比较如图2所示。从图2可以看出,在0~26 GPa压强范围内,理论计算结果与Ciezak-Jenkins等[21]的实验结果具有相同的变化趋势,并且在1和22 GPa处曲线相交,说明本理论计算结果与实验结果有一定的可比性。同时,给出了HBT晶体在压强高于25 GPa时的体积变化情况。理论计算显示,随着压强的增大,晶体体积变化受压强的影响越来越明显,60 GPa以上时压强与体积比呈线性关系。
表 1 HBT晶体晶格参数的计算值和实验值
Table 1. Calculated lattice parameters of HBT crystal along with experimental data
-
通过高压下结构优化,得到不同压力下稳定的HBT晶体结构,在此基础上计算弹性系数和弹性模量。三斜结构的弹性常数矩阵包含21个独立矩阵元。Born[22]d给出了晶体材料力学稳定性的充分必要条件为弹性常数矩阵C正定。按照线性代数知识可知,对称矩阵C正定的充要条件为C的各阶顺序主子式都为正,即
${C_{11}} > 0,\;\;\left| {\begin{array}{*{20}{c}} {{C_{11}}}&{{C_{12}}} \\ {{C_{12}}}&{{C_{22}}} \end{array}} \right| = {C_{11}}{C_{22}} - C_{12}^2 > 0$ $\left| {\begin{array}{*{20}{c}} {{C_{11}}}&{{C_{12}}}&{{C_{13}}} \\ {{C_{12}}}&{{C_{22}}}&{{C_{23}}} \\ {{C_{13}}}&{{C_{23}}}&{{C_{33}}} \end{array}} \right| = ({C_{11}}{C_{22}}{C_{33}} - {C_{11}}C_{23}^2 - C_{12}^2{C_{33}} - C_{13}^2{C_{22}} + 2{C_{12}}{C_{13}}{C_{23}}) > 0$ $\left| {\begin{array}{*{20}{c}} {{C_{11}}}&{{C_{12}}}&{{C_{13}}}&{{C_{14}}} \\ {{C_{12}}}&{{C_{22}}}&{{C_{23}}}&{{C_{24}}} \\ {{C_{13}}}&{{C_{23}}}&{{C_{33}}}&{{C_{34}}} \\ {{C_{14}}}&{{C_{24}}}&{{C_{34}}}&{{C_{44}}} \end{array}} \right| > 0,\;\;\left| {\begin{array}{*{20}{c}} {{C_{11}}}&{{C_{12}}}&{{C_{13}}}&{{C_{14}}}&{{C_{15}}} \\ {{C_{12}}}&{{C_{22}}}&{{C_{23}}}&{{C_{24}}}&{{C_{25}}} \\ {{C_{13}}}&{{C_{23}}}&{{C_{33}}}&{{C_{34}}}&{{C_{35}}} \\ {{C_{14}}}&{{C_{24}}}&{{C_{34}}}&{{C_{44}}}&{{C_{45}}} \\ {{C_{15}}}&{{C_{25}}}&{{C_{35}}}&{{C_{45}}}&{{C_{55}}} \end{array}} \right| > 0$ $\left| {\begin{array}{*{20}{c}} {{C_{11}}}&{{C_{12}}}&{{C_{13}}}&{{C_{14}}}&{{C_{15}}}&{{C_{16}}} \\ {{C_{12}}}&{{C_{22}}}&{{C_{23}}}&{{C_{24}}}&{{C_{25}}}&{{C_{26}}} \\ {{C_{13}}}&{{C_{23}}}&{{C_{33}}}&{{C_{34}}}&{{C_{35}}}&{{C_{36}}} \\ {{C_{14}}}&{{C_{24}}}&{{C_{34}}}&{{C_{44}}}&{{C_{45}}}&{{C_{46}}} \\ {{C_{15}}}&{{C_{25}}}&{{C_{35}}}&{{C_{45}}}&{{C_{55}}}&{{C_{56}}} \\ {{C_{16}}}&{{C_{26}}}&{{C_{36}}}&{{C_{46}}}&{{C_{56}}}&{{C_{66}}} \end{array}} \right| > 0$ 式(5)~式(8)给出了三斜晶体稳定性的充要条件。四阶以上行列式的表达式比较复杂,一般采用直接带入行列式元素来计算验证。图3(a)、 图3(b)、图3(c)、图3(d)给出了21个弹性常数随压强的变化关系。从图3可以看出,C15、C56、C24、C46随压强的变化比较剧烈,其中C15、C56和C24先减小后增大,而C46则先增大后减小随后又增大。将零温零压下晶体的弹性常数代入式(5)~式(8),验证后发现符合稳定性条件。对高压下晶体进行验证,发现也符合稳定性条件,说明高压下HBT晶体是稳定的。
计算体弹模量、剪切模量的普适公式为[23]
$\left\{ \begin{aligned} & 9{B^{\rm V}} = {C_{11}} + {C_{22}} + {C_{33}} + 2{C_{12}} + 2{C_{13}} + 2{C_{23}} \\ & 1/{B^{\rm R}} = {S_{11}} + {S_{22}} + {S_{33}} + 2{S_{12}} + 2{S_{13}} + 2{S_{23}} \\ & 15{G^{\rm V}} = {C_{11}} + {C_{22}} + {C_{33}} - {C_{12}} - {C_{13}} - {C_{23}} + 3{C_{44}} + 3{C_{55}} + 3{C_{66}} \\ & 15/{G^{\rm R}} = 4({S _{11}} + {S_{22}} + {S_{33}} - {S_{12}} - {S_{13}} - {S_{23}} + 3{S_{44}} + 3{S_{55}} + 3{S_{66}}) \end{aligned} \right.$ 根据Voigt-Reuss-Hill均值方法,有
$B = \frac{{{B^{\rm V}} + {B^{\rm R}}}}{2},\;\;G = \frac{{{G^{\rm V}} + {G^{\rm R}}}}{2}$ 再根据B、G可以获得体弹模量E和泊松比
$\mu$ $E = \frac{{9BG}}{{3B + G}},\;\;\mu = \frac{{3\dfrac{B}{G} - 2}}{{6\dfrac{B}{G} + 2}}$ 依据式(9)~式(11),可以获得体弹模量、剪切模量和杨氏模量随压强的变化关系,如图4所示。从图4可以看出,随着压强的增大,3种模量的数值也增大,其中体弹模量在15 GPa以上时表现出与压强的线性关系。而3种模量中,杨氏模量的变化最显著,高压下的数值也最大。此外,根据Pugh[24]的理论,B/G可作为鉴别晶体韧脆性的标准,即当B/G<1.74时,晶体材料表现为脆性;反之,晶体材料表现为韧性。因此,根据式(11)中泊松比
$\mu$ 与B/G的关系来计算B/G的数值,并作出B/G与压强的关系曲线,如图5所示。随着压强增大,B/G逐渐增大,与图4的变化趋势一致。当压强低于85 GPa时,B/G<1.74,材料表现出脆性;而当压强高于85 GPa时,HBT晶体表现出韧性。类似的层状分子晶体TATB在0~50 GPa时表现为韧性[25],而同为层状分子晶体的HBT在0~85 GPa下却表现为脆性,说明HBT晶体对压力的响应较为敏感。图 4 体弹、剪切和杨氏模量随压强的变化关系
Figure 4. Variations of bulk modulus B, shear modulus G and Young’s modulus E with pressure
同时,根据以下公式可以得到HBT晶体材料的压缩波声速
$v_{\rm p}$ 、剪切波声速${v_{\rm{s}}}$ 及平均声速${v _{\rm{m}}}$ ${v\rm_ p} = \sqrt {\frac{{3B + 4G}}{{3\rho }}} $ ${v\rm_s} = \sqrt {{G}/{\rho }} $ ${v\rm_m} = {\left[ {\frac{1}{3}\left( {\frac{1}{{v \rm_p^3}} + \frac{2}{{v \rm_s^3}}} \right)} \right]^{ - 1/3}}$ 式中:ρ为HBT晶体密度。平均声速νm与德拜温度
${\varTheta} $ 满足如下关系[26]${\varTheta}= \frac{h}{{{k\rm_B}}}{\left[ {\frac{{3n}}{{4\text{π} }}\left( {\frac{{{N\rm_A}\rho }}{M}} \right)} \right]^{1/3}}{v \rm_m}$ 式中:h为普朗克常数,kB为玻尔兹曼常数,n为原子数,NA为阿伏加德罗常数,M为相对分子质量。图6给出了压缩波声速
${v _{\rm{p}}}$ 、剪切波声速${v _{\rm{s}}}$ 、平均声速${v _{\rm{m}}}$ 和德拜温度${\varTheta} $ 随压强的变化关系。可见,声速随压强的增大而增大,其中压缩波声速${v _{\rm{p}}}$ 最大,且随压强变化最明显。在0 GPa下得到的德拜温度为975 K,而120 GPa时则为1 820 K,压强对德拜温度的影响明显。 -
弹性性质的各向异性是材料力学性能通常受微裂纹和晶格畸变影响的重要条件,因此,研究弹性性质的各向异性对于提高材料的力学性能有着重要的意义[27]。目前有3种处理方案。第1种方案是Chung等[28]提出的体弹模量、杨氏模量和剪切模量分数比各向异性的概念,定义为
${A_B} = \frac{{{B^{\rm V}} - {B^{\rm R}}}}{{{B^{\rm V}} + {B^{\rm R}}}}$ ${A_E} = \frac{{{E^{\rm V}} - {E^{\rm R}}}}{{{E^{\rm V}} + {E^{\rm R}}}}$ ${A_G} = \frac{{{G^{\rm V}} - {G^{\rm R}}}}{{{G^{\rm V}} + {G^{\rm R}}}}$ 式(16)~式(18)分别用于表征材料的体弹模量、杨氏模量和剪切模量的各向异性程度。据此可以计算出模量各向异性分数比与压强的关系,如图7所示。对于弹性各向同性材料,AB、AE和AG都为零,而AB、AG为1表示最大可能的弹性各向异性。从图7可以看出,AB、AE和AG均大于零而小于1,说明在0~120 GPa的压强范围内HBT晶体的弹性模量表现出各向异性。此外,体弹模量分数比只有很小的波动,杨氏模量和剪切模量则表现出相同的变化趋势,即先增大后逐渐减小,说明HBT晶体在高压下的各向异性程度减弱。
第2种方案是Ranganathan等[29]提出的全局弹性各向异性系数(AU)的概念,用来描述材料弹性模量的各向异性,具体表达式如下
${A^{\rm U}} = 5\frac{{{G^{\rm V}}}}{{{G^{\rm R}}}} + \frac{{{B^{\rm V}}}}{{{B^{\rm R}}}} - 6$ ${A^{\rm C}} = \frac{{{G^{\rm V}} - {G^{\rm R}}}}{{2G}}$ 图8的纵轴表示0~120 GPa压强范围内的全局弹性各向异性系数AU的数值,其中,最高点表示0 GPa下的数值。从结果看都大于零,因此,HBT晶体在高压下表现出各向异性。
此外,式(20)的AC与式(18)的AG是相等的,都是由Chung等于1967年提出,均用于表征剪切模量的各向异性。为了比较两种方案的关系,作AU-AC曲线。可见,两种方法都很好地描述了HBT晶体在0~120 GPa压强范围弹性模量的各向异性特征。随着压强的增大,AU、AC和AU/AC的数值都在减少,说明弹性模量的各向异性程度在减弱,两种方案得到的结论一致。
第3种方案是Toher总结和发展的Ranganathan关于材料弹性各向异性的理论表述[30]
${A^{\rm L}} = \sqrt 5 \ln \left(1 + \frac{{{A^{\rm U}}}}{5}\right)$ AL随压强的变化关系如图9所示。可以发现,曲线的变化趋势与第1种方案(弹性模量分数比理论)得到的结果类似,即AL先增大后减小,5 GPa处出现转折。整体上看,高压下HBT晶体弹性模量的各向异性程度减弱。由此可知,AC、AU和AL三者在描述HBT晶体各向异性程度上是等价的。
从图7、图9可以明显看出,5、25和55 GPa下曲线有较为明显的变化。为了进一步考察这3个压强下HBT晶体模量的各向异性程度,利用ELATE开源软件[31],将这3个压强下的弹性系数6×6矩阵元36个数值全部代入,绘制出剪切模量的二维分布,如图10所示。因单位换算,图10中坐标轴均被放大了10倍。图10(a)、图10(b)、图10(c)分别显示了5、25和55 GPa下剪切模量的二维分布。每个压强下的3幅图分别代表剪切模量三维分布在xy、xz、yz 3个面上的投影。蓝线和绿线分别代表
$\overrightarrow a $ 、$\overrightarrow b $ 两个方向剪切模量的矢量投影。从图10可以看出:压强越大,每个面上投影矢量的各向异性越强;在确定的压强下,每个面上的投影也表现出差异性。综上所述,HBT晶体剪切模量的各向异性受压强的影响显著。 -
采用基于密度泛函理论的第一性原理方法计算了HBT晶体在零温零压及高压下的结构,研究了高压下HBT的力学稳定性和弹性性质,具体分析了HBT晶体的弹性常数Cij、弹性模量(B、G、E)、声速(
${v _{\rm{p}}}$ 、${v _{\rm{s}}}$ 、${v _{\rm{m}}}$ )和德拜温度Θ随压强的变化关系。在此基础上,基于3种不同理论研究了弹性模量的各向异性。结果表明,弹性常数和模量都随压强的增加而增加。根据B/G估测了HBT晶体在高压下的韧脆性。弹性模量各向异性参量随压强表现出不同的变化规律,总体而言,高压下HBT晶体的弹性模量各向异性程度减弱。
高压下HBT晶体的弹性性质
Elastic Properties of HBT Crystal under High Pressure
-
摘要: 采用基于密度泛函理论的第一性原理方法,研究了HBT晶体在常压和高压下的晶格常数和弹性性质及各向异性,利用3种不同理论模型研究了高压下HBT晶体的各向异性性质。研究结果表明:高压下HBT晶体的弹性常数和弹性模量显著增加,晶体表现出高压韧性;同时,高压下HBT晶体具有较大的弹性模量和力学各向异性,随着压强的增大,HBT晶体的各向异性程度减弱。此外,热力学性质计算结果表明,HBT晶体具有较高的德拜温度,并且德拜温度随着压强的增大而升高。Abstract: The lattice constants, elastic and anisotropy properties of HBT crystals at normal and high pressure were investigated by using the first-principle method based on density functional theory. The anisotropic properties of HBT crystal under high pressure were studied by using three different theoretical models. The results show that the elastic constant and elastic modulus of HBT crystal increase significantly under high pressure, and HBT crystal shows high-pressure toughness. Simultaneously, HBT crystal has large elastic modulus and mechanical anisotropy under high pressure. With increasing pressure, the extent of anisotropy of HBT crystal decreases. In addition, the thermodynamic properties show that HBT crystal has a higher Debye temperature, which increases with increasing pressure.
-
Key words:
- HBT crystal /
- high pressure /
- elastic properties /
- anisotropic
-
-
[1] 范开敏, 杨莉, 孙庆强, 等. 六角相ErA x (A = H, He)体系弹性性质的第一性原理研究 [J]. 物理学报, 2013, 62(11): 116201. doi: 10.7498/aps.62.116201
FAN K M, YANG L, SUN Q Q, et al. First-principles study on elastic properties of hexagonal phase ErAx(A = H, He) [J]. Acta Physica Sinica, 2013, 62(11): 116201. doi: 10.7498/aps.62.116201[2] ANIRBAN P, CATALIN R P. Contribution of molecular flexibility to the elastic–plastic properties of molecular crystal α-RDX [J]. Modelling and Simulation in Materials Science and Engineering, 2017, 25(1): 015006. doi: 10.1088/1361-651X/25/1/015006 [3] RYKOUNOV A A. Investigation of the pressure dependent thermodynamic and elastic properties of 1,3,5-triamino-2,4,6-trinitrobenzene using dispersion corrected density functional theory [J]. Journal of Applied Physics, 2015, 117(21): 215901. doi: 10.1063/1.4921815 [4] LIAN D, LU L Y, WEI D Q. High-pressure behavior of β-HMX crystal studied by DFT-LDA [J]. Chinese Physics Letters, 2008, 25(3): 899–902. doi: 10.1088/0256-307X/25/3/026 [5] WU Q, ZHU W H, XIAO H M. DFT study on crystalline 1,1-diamino-2,2-dintroethylene under high pressures [J]. Journal of Molecule Modeling, 2013, 19(9): 4039–4047. doi: 10.1007/s00894-013-1931-8 [6] XU X J, ZHU W H, XIAO H M. DFT studies on the four polymorphs of crystalline CL-20 and the influences of hydrostatic pressure on epsilon-CL-20 crystal [J]. The Journal of Physical Chemistry B, 2007, 111(8): 2090–2097. doi: 10.1021/jp066833e [7] KLAPÖTKE T M, SABATE C M. Bistetrazoles: nitrogen-rich, high-performing, insensitive energetic compounds [J]. Chemistry of Materials, 2008, 20(11): 3629–3637. doi: 10.1021/cm703657k [8] KLAPÖTKE T M, SABATE C M. 5, 5′-Hydrazinebistetrazole: an oxidation-stable nitrogen-rich compound and starting material for the synthesis of 5, 5′-Azobistetrazolates [J]. Zeitschrift für Anorganische und Allgemeine Chemie, 2007, 633(15): 2671–2677. doi: 10.1002/zaac.200700271 [9] KARAGHIOSOFF K, KLAPOTKE T M, SABATE C M. Nitrogen-rich compounds in pyrotechnics: alkaline earth metal salts of 5, 5′-Hydrazine-1, 2-diylbis(1H-tetrazole) [J]. European Journal of Inorganic Chemistry, 2009, 2: 238–250. [10] EBESPACHER M, KLAPOTKE T M, SABATE C M. Nitrogen-rich alkali metal 5, 5′-Hydrazinebistetrazolate salts: environmentally friendly compounds in pyrotechnic mixtures [J]. New Journal of Chemistry, 2009, 33: 517–527. doi: 10.1039/B818927G [11] DE LUCIA F C, GOTTFRIED J L. Characterization of a series of nitrogen-rich molecules using laser induced breakdown spectroscopy [J]. Propellants Explosives Pyrotechnics, 2010, 35(3): 268–277. doi: 10.1002/prep.201000009 [12] GIANNOZZI P, ANDREUSSI O, BRUMME T, et al. Advanced capabilities for materials modelling with Quantum ESPRESSO [J]. Journal of Physics: Condensed Matter, 2017, 29(46): 465901. doi: 10.1088/1361-648X/aa8f79 [13] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Physical Review Letter, 1996, 77(18): 3865–3868. doi: 10.1103/PhysRevLett.77.3865 [14] PERDEW J P, BURKE K, ERNZERHOF M. Errata: generalized gradient approximation made simple [J]. Physical Review Letter, 1997, 78(7): 1396. [15] HAMANN D R. Optimized norm-conserving vanderbilt pseudopotentials [J]. Physical Review B, 2013, 88(8): 085117. doi: 10.1103/PhysRevB.88.085117 [16] FLETCHER R. Practical methods of optimization [M]. New York: Wiley, 1980: 126−127. [17] MONKHORST H J, PACK J D. Special points for brillouin-zone integrations [J]. Physical Review B, 1976, 13(12): 5188–5192. doi: 10.1103/PhysRevB.13.5188 [18] CORSO D A. Pseudopotentials periodic table: from H to Pu [J]. Computational Materials Science, 2014, 95: 337–350. doi: 10.1016/j.commatsci.2014.07.043 [19] CORSO D A. Elastic constants of beryllium: a first-principles investigation [J]. Journal of Physics: Condensed Matter, 2016, 28(7): 075401. doi: 10.1088/0953-8984/28/7/075401 [20] MOUHAT F, COUDERT F X. Necessary and sufficient elastic stability conditions in various crystal systems [J]. Physical Review B, 2014, 90(22): 224104. doi: 10.1103/PhysRevB.90.224104 [21] CIEZAK-JENKINS J A, JENKINS T A. Shear induced weakening of the hydrogen bonding lattice of the energetic material 5,5′-Hydrazinebistetrazole at high-pressure [J]. Journal of Molecular Structure, 2017, 1129: 313–318. doi: 10.1016/j.molstruc.2016.09.084 [22] BORN M. On the stability of crystal lattices. I [J]. Mathematical Proceedings of the Cambridge Philosophical Society, 1940, 36(2): 160–172. doi: 10.1017/S0305004100017138 [23] HILL R. The elastic behaviour of a crystalline aggregate [J]. Proceedings of the Physical Society (Section A), 1952, 65(5): 349–354. doi: 10.1088/0370-1298/65/5/307 [24] PUGH S F. XCII Relations between the elastic moduli and the plastic properties of polycrystalline pure metals [J]. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1954, 45(367): 823–843. doi: 10.1080/14786440808520496 [25] 范航, 聂福德, 龙瑶, 等. 钝感高能炸药三氨基三硝基苯高温高压下热力学性质的分子动力学模拟研究 [J]. 物理学报, 2016, 65(6): 066201. doi: 10.7498/aps.65.066201
FAN H, NIE F D, LONG Y, et al. A molecular dynamics simulation of thermodynamic properties of 1, 3, 5-triamino-2, 4, 6-trinitrobenzene under high pressure and high temperature [J]. Acta Physica Sinica, 2016, 65(6): 066201. doi: 10.7498/aps.65.066201[26] ANDERSON O L. A simplified method for calculating the Debye temperature from elastic constants [J]. Journal of Physics and Chemistry of Solids, 1963, 24(7): 909–917. doi: 10.1016/0022-3697(63)90067-2 [27] 陈志谦, 林兴香, 李春梅. 硅化锆(ZrSi)弹性及其各向异性研究 [J]. 西南大学学报(自然科学版), 2016, 38(3): 1–8.
CHEN Z Q, LIN X X, LI C M. Investigations of elasticity anisotropy of ZrSi [J]. Journal of Southwest University (Natural Science Edition), 2016, 38(3): 1–8.[28] CHUNG D H, BUESSEM W R. The elastic anisotropy of crystals [J]. Journal of Applied Physics, 1967, 38(5): 2010–2012. doi: 10.1063/1.1709819 [29] RANGANATHAN S I, OSTOJA-STARZEWSKI M. Universal elastic anisotropy index [J]. Physical Review Letter, 2008, 101(5): 055504. doi: 10.1103/PhysRevLett.101.055504 [30] TOHER C, PLATA J J, LEVY O, et al. High-throughput computational screening of thermal conductivity, Debye temperature, and Grüneisen parameter using a quasiharmonic Debye model [J]. Physical Review B, 2014, 90(17): 174107. doi: 10.1103/PhysRevB.90.174107 [31] GAILLAC R, PULLUMBI P, COUDERT F X. ELATE: an open-source online application for analysis and visualization of elastic tensors [J]. Journal of Physics: Condensed Matter, 2016, 28(27): 275201. doi: 10.1088/0953-8984/28/27/275201 -