高压下G<sub>2</sub>ZT晶体结构、电子结构和光学性质的第一性原理研究

李佐; 刘云; 廖大麟; 成丽红

doi:10.11858/gywlxb.20220514

高压下G₂ZT晶体结构、电子结构和光学性质的第一性原理研究

doi: 10.11858/gywlxb.20220514

贵州工程应用技术学院理学院, 贵州毕节　551700

基金项目: 贵州省科技计划项目（[2020]1Y023）；贵州省教育厅青年科技人才成长项目（KY[2018]398）

详细信息

作者简介:
李佐（1983－），男，硕士，讲师，主要从事高压下含能材料物性研究. E-mail：lizuo212@163.com

中图分类号: O521.2
计量
- 文章访问数: 400
- HTML全文浏览量: 327
- PDF下载量: 53
出版历程
- 收稿日期: 2022-02-17
- 修回日期: 2022-03-06
- 录用日期: 2022-03-06
- 网络出版日期: 2022-06-22
- 刊出日期: 2022-07-28

First-Principles Study on Structural, Electronic and Optical Properties of G₂ZT Crystal under High Pressure

School of Science, Guizhou University of Engineering Science, Bijie 551700, Guizhou, China

摘要

摘要: 基于密度泛函理论的第一性原理研究了高压下富氮含能材料（双3, 4, 5-三氨基-1, 2, 4-三唑）-5, 5′-偶氮四唑（G₂ZT）的几何结构、电子结构和光学性质。结果表明，在考虑范德瓦尔斯色散修正和密度泛函色散修正的情况下, 分子晶体结构数据与实验结果的相对误差均在3%以内。Hirshfeld表面分析结果表明，随着压强增大，分子间氢键的相互作用减弱。G₂ZT晶体在零压下的能带带隙为2.03 eV，是一种p型半导体。随着压强增大，带隙变窄，吸收系数可达到3.0×10⁶ cm⁻¹。研究结果为进一步分析高压下G₂ZT晶体的特征提供了理论参考。
- G₂ZT晶体 /
- 高压 /
- 电子结构 /
- 光学性质
Abstract: Geometric structure, electronic structure and optical properties of nitrogen-rich energetic materials (bis 3, 4, 5-triamino-1, 2, 4-triazole)-5, 5′-azotetrazole (G₂ZT) at high pressures are investigated using first-principles based on density functional theory. The calculated results obtained by using vdW-DF2 and PBE-D2 methods show that the crystal structure data fit well with the experimental results, and the error rates are all within 3%. Hirshfeld surface analyses indicate that interactions of the inter-molecular hydrogen bond are weaken with the increasing pressure. G₂ZT crystal possesses a band gap of 2.03 eV at zero pressure, and it is a p-type semiconductor. As the pressure increases, the band gap becomes narrower and the absorption coefficient can approach 3.0×10⁶ cm⁻¹. These results provide a theoretical reference for further analysis of G₂ZT crystal’s characteristics under high pressure.
- G₂ZT crystal /
- high pressure /
- electronic structure /
- optical properties

HTML全文

破片式战斗部爆炸后，会形成破片和冲击波两种毁伤元。在战斗部爆炸初始时刻，炸药爆炸产生的高温高压气体急剧膨胀，使壳体碎裂，爆轰产物溢出，破片在爆轰产物的作用下由静止状态逐渐加速到最大初速。在爆炸后一段时间内，冲击波运动在破片之前，随着运动距离的增加，冲击波速度衰减，破片运动在冲击波之前。因此，在破片和冲击波运动过程中，存在一个相遇位置。而在相遇位置前后，破片和冲击波作用于目标的顺序不同，对目标的毁伤效果也不相同。因此，对破片和冲击波相遇位置的研究具有重要的意义。针对破片和冲击波在空气中的传播特性，相关学者进行了大量研究^[1-5]。刘刚^[6]采用ANSYS/LS-DYNA非线性软件对破片式战斗部的爆炸过程进行了数值计算，得到破片和冲击波在12.6倍装药直径处相遇；梁为民等^[7]对破片式战斗部爆炸后破片和冲击波的运动过程进行了理论推导，分析了不同比例距离和装填系数下，破片与冲击波的运动规律，并通过试验得到了破片与冲击波的相遇位置，约为装药半径的19.3倍；Grisaro^[8]、Nyström^[9]等用理论计算的方法对破片和冲击波相遇位置进行了研究。目前，对破片式战斗部爆炸后破片和冲击波相遇位置的研究主要是通过理论计算和数值仿真的方法，通过试验的方法对破片和冲击波相遇位置的研究较少。

针对当前研究的不足，设计了一个破片式战斗部，对破片和冲击波相遇位置进行了数值计算，并设计了一种测量破片和冲击波相遇位置的试验方法，在此基础上分析了战斗部的炸药装填系数、破片质量、爆热和爆速等因素对破片和冲击波相遇位置的影响规律。

1. 破片和冲击波相遇位置数值计算

1.1 战斗部结构

为研究破片和冲击波的相遇位置，设计了一个破片式战斗部，其结构如图 1所示。其外部是环形的预制破片层，中间为铝合金壳体，内部装填1 kg的8701炸药；战斗部长度L为160 mm，装药半径R为31.3 mm，破片层厚度t为5.5 mm，单枚破片高h为8 mm，破片材料为钢，单枚破片质量为3.3 g，破片数目为480枚。

图 1 战斗部结构示意图

Figure 1. Schematic structure of the warhead

下载: 全尺寸图片幻灯片

1.2 战斗部数值计算模型

采用ANSYS/LS-DYNA对战斗部爆炸后毁伤元形成过程进行了数值计算。建立了如图 2所示的1/2数值计算模型，该模型由炸药、壳体、破片和空气域组成。空气域整体尺寸为300 cm×300 cm×1200 cm，网格大小为4 mm^[1]。破片、壳体采用Lagrange网格, 炸药和空气采用Euler网格，使用多物质单元ALE算法。采用关键字*CONSTRAINED_LAGRANGE_IN_SOLID定义Lagrange网格和空气网格间的耦合算法，采用中心点火的方式起爆战斗部，在模型的对称面设置对称边界条件，空气域边界采用透射边界。网格单元类型均为八节点六面体SOLID164实体单元，采用g-cm-μs单位制。

图 2 战斗部三维计算模型(局部)

Figure 2. Three-dimensional finite element model of the warhead (part)

下载: 全尺寸图片幻灯片

战斗部壳体材料采用Johnson_Cook模型和Grüneisen状态方程描述，破片材料采用弹塑性Plastic_Kinematic本构模型，材料主要参数见表 1^[10]。用High_Explosive_Burn本构模型和JWL状态方程描述炸药，空气材料采用Null材料模型和Linear_Polynomial状态方程描述，材料主要参数分别见表 2和表 3^[11]。

表 1 壳体和破片的主要材料参数表

Table 1. Material parameters of shell and fragments

Material	ρ/(kg·m^-3)	E/GPa	μ	A/MPa	B/MPa	n	m	c/(m·s^-1)	S₁	γ
Alloy	2797	69.63	0.33	265	426	0.34	1	5280	1.4	2
Steel	7850	210	0.28	335

下载: 导出CSV

| 显示表格

表 2 炸药材料参数及状态方程参数

Table 2. Material and equation of state (EOS) parameters for explosive

ρ/(kg·m^-3)	D/(m·s^-1)	A/GPa	B/GPa	R₁	R₂	ω	p_CJ/GPa	V₀
1780	8390	581.4	6.8	4.1	1	0.35	34	1

下载: 导出CSV

| 显示表格

表 3 空气材料参数及状态方程参数

Table 3. Material and equation of state (EOS) parameters for air

ρ/(kg·m^-3)	C₀/MPa	C₁	C₂	C₃	C₄	C₅	C₆	V	e₀/MPa
1.25	-0.1	0	0	0	0.4	0.4	0	1	0.25

下载: 导出CSV

| 显示表格

1.3 数值计算结果

战斗部爆炸后破片和冲击波形成过程如图 3所示，战斗部装药爆炸后，壳体及破片在高温高压的爆轰产物下膨胀，最后破裂，泄露的气体产物压缩周围空气形成冲击波。由图 3可知，柱形装药爆炸后产生的冲击波传播一段距离后，冲击波波阵面近似以球面波的形式向四周传播。在爆炸初始时刻，冲击波速度高于破片速度，所以在425 μs之前冲击波运动在破片之前；随着运动距离的增加，冲击波的强度衰减很快，相应的传播速度也逐渐降低，而破片的速度衰减相对冲击波较慢，所以运动一段时间后，在425 μs时刻破片追赶上冲击波和其相遇；在425 μs之后，破片超过冲击波，运动在冲击波之前。

图 3 不同时刻破片和冲击波相对位置

Figure 3. Relative locations of fragments and blast wave for different time instants

下载: 全尺寸图片幻灯片

1.4 破片和冲击波相遇位置

图 4为起爆点对应的中间层破片和冲击波波阵面前沿运动距离R_C随时间的变化曲线。由图 4可知：在所计算的时间范围内破片的运动距离随时间近似线性增加，即速度变化不大；而冲击波运动距离的增量随时间的增加逐渐变缓，表明冲击波的运动速度逐渐降低。在425 μs时，破片和冲击波相遇，相遇点距起爆点0.72 m处(约为11.5倍装药直径)，在此位置之前，破片运动在冲击波之后；在此位置之后，破片运动在冲击波之前，随后两者之间的距离差逐渐增大。

图 4 传播距离随时间的变化(数值计算)

Figure 4. Variation of propagation distances with time (Numerical computation)

下载: 全尺寸图片幻灯片

2. 破片和冲击波相遇位置试验

2.1 试验设计

为了研究破片式战斗部爆炸后破片和冲击波的相遇位置，设计了一种测量相遇位置的试验方法。试验现场布置及试验原理如图 5所示，试验中使用图 1所示的战斗部，在平行于战斗部轴线不同位置处放置5块硬质木块(24 mm×54 mm×240 mm)，战斗部起爆后，通过高速录像记录木块的倾斜状态，当冲击波到达木块时，木块会发生倾斜，以此时刻作为冲击波到达对应位置处的时间, 每发试验木块距战斗部的距离如表 4所示。在木块后面布置一块白色背景布，便于高速摄像机捕捉目标位置；在距战斗部一定距离处设置测速靶，用于测量战斗部爆炸后破片的速度；木块、测速靶和战斗部的中心保持在同一平面上。

图 5 试验装置及布置示意图

Figure 5. Photograph of experimental setup and the arrangement of test devices

下载: 全尺寸图片幻灯片

表 4 木块距战斗部距离

Table 4. Distance between wood blocks and warhead

Shot	Distance/m
Shot	Wood block 1	Wood block 2	Wood block 3	Wood block 4	Wood block 5
1	0.4	0.6	1.2	1.6	2.0
2	1.6	2.0	2.4	2.8	3.2
3	1.6	2.0	2.8	3.2	3.6

下载: 导出CSV

| 显示表格

距离战斗部起爆点分别为2.5、4.3、6.2、8.0和10.0 m处的地面上分别布置kistler 211B4压电传感器，用于测量不同距离处冲击波超压，现场布置和试验原理如图 6所示。

图 6 冲击波超压测试现场布置图

Figure 6. Schematic diagram of experimental setup for blast wave overpressure measurement

下载: 全尺寸图片幻灯片

采用中心起爆的方式起爆战斗部，在战斗部壳体上缠绕细铜丝，当战斗部爆炸壳体碎裂时，铜丝断裂，作为高速录像和压电传感器的启动零时刻。高速录像的帧数为24000帧每秒。

2.2 试验过程

共进行了3发试验，图 7为高速录像记录的第2发试验中破片和冲击波典型时刻的运动过程。如图 7(c)所示，第1块木块开始倾斜，以1.791 ms作为冲击波达到第1块木块的时间；如图 7(d)所示，第2块木块发生倾斜，以2.334 ms作为冲击波到达第2块木块的时间。试验之后的木块如图 8所示，由图 8可知，木块表面没有被破片击中的痕迹，表明木块是在冲击波的作用下发生了移动，排除了木块是因为破片击中而发生移动的可能性。

图 7 破片和冲击波运动过程图

Figure 7. Relative locations of fragments and blast wave

下载: 全尺寸图片幻灯片

图 8 试验后木块照片

Figure 8. Wood blocks photo after test

下载: 全尺寸图片幻灯片

2.3 试验结果及分析

2.3.1 破片和冲击波相遇位置

3发试验中，测速靶测得破片的平均速度为1600 m/s。冲击波到达各个测点的时间如表 5所示，由于战斗部爆炸后火光太强，部分试验数据无法获得。将每发试验距战斗部相同位置处测得的试验数据取平均值，得到距战斗部不同位置处冲击波到达时间，如表 6所示。

表 5 冲击波到达时间

Table 5. Arrival time of blast wave

Shot	Time/ms
Shot	Wood block 1	Wood block 2	Wood block 3	Wood block 4	Wood block 5
1			0.752	1.832	2.586
2	1.791	2.334	3.009	3.825	4.771
3		2.876	3.762	4.689

下载: 导出CSV

| 显示表格

表 6 冲击波到达平均时间

Table 6. Average arrival time of blast wave

Distance/m	Time/ms
1.2	0.752
1.6	1.812
2.0	2.599
2.4	3.009
2.8	3.794
3.2	4.730

下载: 导出CSV

| 显示表格

根据冲击波到达不同距离处的时间，使用幂函数对冲击波传播距离随时间的变化进行拟合^[11]，由于破片在短时间内速度衰减不明显, 所以破片速度近似取1600 m/s，破片和冲击波运动距离随时间的变化曲线如图 9所示，由图可知破片和冲击波在465 μs时相遇，相遇点距战斗部爆炸中心0.76 m(约为12.14倍装药直径)，数值计算得到的相遇距离与试验测得的相遇距离相比，误差为5.26%。

图 9 传播距离随时间的变化(试验)

Figure 9. Variation of propagation distances with time (Test)

下载: 全尺寸图片幻灯片

2.3.2 冲击波超压值

在1.3节中对毁伤元形成过程进行了数值计算，图 10是数值计算中对应的测点处冲击波超压随时间的变化曲线。由图 10可知，随着冲击波传播距离的增加，冲击波超压峰值急剧降低，超压作用时间增加。表 7显示地面5个测点的冲击波超压峰值的3次试验平均值和数值计算的对比结果，误差均在9.05%以下，说明该数值计算具有一定的合理性。在此基础上，下文以相同的数值计算方法，计算研究战斗部装填系数、破片质量、爆热和爆速对破片和冲击波相遇位置的影响。

图 10 不同点处超压-时间曲线

Figure 10. Shock wave overpressure versus time at different locations

下载: 全尺寸图片幻灯片

表 7 计算超压值与试验超压值对比

Table 7. Comparison of the numerical and experimental overpressure results

Distance/m	Over pressure/MPa		Error/%
Distance/m	Test	Numerical	Error/%
2.5	0.356	0.3557	8.43
4.3	0.086	0.0912	6.05
6.2	0.04565	0.0489	7.12
8.0	0.0256	0.0241	5.86
10.0	0.0199	0.0217	9.05

下载: 导出CSV

| 显示表格

3. 破片和冲击波相遇位置影响因素分析

为了研究战斗部结构参数对破片和冲击波相遇位置的影响规律，改变图 1中战斗部的结构参数，战斗部结构参数和计算方案设计如表 8所示。方案1、方案2和方案3分别研究装填系数k、破片质量m、炸药类型对破片和冲击波相遇位置的影响。

表 8 战斗部及破片参数(L=260 mm，R=96 mm)

Table 8. Parameters of warhead and fragment (L=260 mm, R=96 mm)

Arrangement	k/%	m/g	t/mm	h/mm	Explosive
	37.31	6.00	9	10	8701
1	40.33	6.00	8	10	8701
	43.73	6.00	7	10	8701
	52.34	2.12	5	4	8701
	52.34	4.24	5	6	8701
2	52.34	6.00	5	8.2	8701
	52.34	8.48	5	9.6	8701
	52.34	12.00	5	12	8701
3		6.00	5	8.2	8701/TNT/B/HMX

下载: 导出CSV

| 显示表格

3.1 装填系数对破片和冲击波相遇位置的影响

不同装填系数下战斗部爆炸后形成的破片和冲击波相遇位置的关系曲线如图 11所示。由曲线可以看出，随着炸药装填系数的增加，冲击波和破片的相遇位置距爆炸中心的距离越小。这主要是因为装填系数增大时，爆炸后破片的初速也提高，虽然冲击波的速度也随着装填系数的增加而增加，但是在空气中破片速度的衰减幅度比冲击波的衰减幅度小，所以破片在距爆炸中心更近的距离内赶上了冲击波波阵面。

图 11 装填系数k对相遇位置的影响

Figure 11. Influence of k on meeting location

下载: 全尺寸图片幻灯片

3.2 单枚破片质量对破片和冲击波相遇位置的影响

破片质量的变化对破片和冲击波相遇位置的影响如图 12所示。由图可知，随着单枚破片质量的增大，破片和冲击波的相遇位置距爆炸中心的距离减小；因为单枚破片质量增大，其抗速度衰减的能力增强，追上冲击波的距离缩短，所以破片和冲击波的相遇位置距爆炸中心的距离随单枚破片质量的增大而减小。

图 12 单枚破片质量对相遇位置的影响

Figure 12. Influence of fragment mass on meeting location

下载: 全尺寸图片幻灯片

但随着破片质量的提高，相遇位置距爆炸中心的距离减小幅度很小，破片质量平均增加1倍，相遇位置距爆炸中心的距离减小2.4%，表明破片质量对相遇位置的影响程度很小。这主要是因为在战斗部爆炸后，在较短时间内破片在空气中的速度衰减幅度较低，而质量差别不大的破片在空气中的速度衰减幅度差别不大，所以破片质量对破片和冲击波相遇位置的影响不大。

3.3 装药类型对破片和冲击波相遇位置的影响

方案3中，战斗部结构参数不变，只改变装药类型，战斗部装药类型的差异主要体现在装药爆热及爆速上。为方便数值计算，将不同类型的炸药装药均转换为TNT当量质量W_T，W_T的计算公式为

${W_{\rm{T}}} = {W_{\rm{i}}}\frac{{{Q_{{\rm{Ti}}}}}}{{{Q_{{\rm{TTNT}}}}}}$

(1)

式中:W_i为该炸药质量，单位kg；Q_Ti为该炸药爆热，单位J/kg，Q_TTNT为TNT的爆热。不同类型炸药的相关参数与TNT当量如表 9所示^[4]。

表 9 不同类型炸药的相关参数

Table 9. Parameters of different explosives

Explosive	ρ/(kg·m^-3)	Q_V/(m·s^-1)	Q_T/(J·g^-1)	Explosive mass/kg	TNT equivalent	k/%
TNT	1.58	6856	4225	11.00	1.00	44.51
B	1.71	7680	4690	12.37	1.11	50.05
8701	1.72	7980	5300	11.95	1.25	52.34
HMX	1.89	9100	6188	13.13	1.46	58.30

下载: 导出CSV

| 显示表格

根据TNT当量质量计算了装填系数对破片和冲击波的相遇位置的影响结果，如图 13所示。炸药的爆热、爆速增大时，破片和冲击波的相遇位置距爆炸中心的距离减小。主要原因为：一方面，战斗部装药爆热的增加，其装药等效的TNT当量也增加，进而炸药的装填系数也增加，所以，装药爆热对破片和冲击波相遇位置的影响规律和装填系数对相遇位置的影响相同；另外，由计算破片初速的格林公式可知，格林常数与炸药的爆速有关，且近似呈线性变化。炸药的爆速越高，格林常数越大，破片的初始速度越高，破片追上冲击波波阵面的距离缩短。

图 13 不同炸药类型对相遇位置的影响

Figure 13. Influence of different types of explosives on meeting location

下载: 全尺寸图片幻灯片

4. 结论

通过数值计算和试验的方法，对破片和冲击波的相遇位置进行了研究，分析了破片和冲击波相遇位置的影响因素及其影响规律。

(1) 采用木块测量冲击波传播距离和测速靶测破片初速的试验方法，可以有效地测量破片和冲击波的相遇位置。

(2) 破片和冲击波相遇位置随炸药装填系数的增加而减小，装填系数增加31%，相遇点距爆炸中心的距离减小11.5%。

(3) 装填系数相同的情况下，单枚破片的质量越大，破片和冲击波的相遇位置距爆炸中心的距离减小；破片质量增加1倍，相遇距离减小2.4%，可见破片质量的变化对相遇位置的影响较小。

(4) 炸药的爆热、爆速越大，破片和冲击波的相遇位置距爆炸中心的距离越小。

图 1 G₂ZT晶体的晶胞结构和分子模型

Figure 1. Cell structure and molecular model of G₂ZT crystal

下载: 全尺寸图片幻灯片

图 2 G₂ZT晶体的晶格参数和体积比随压强的变化

Figure 2. Lattice parameters and V/V₀ of G₂ZT crystal versus pressure

下载: 全尺寸图片幻灯片

图 3 C₂N₁₀分子在0 GPa下的Hirshfeld面和二维指纹图

Figure 3. Hirshfeld surface and two-dimensional finger patterns of C₂N₁₀ molecule at 0 GPa

下载: 全尺寸图片幻灯片

图 4 C₂N₁₀分子在40 GPa下的Hirshfeld面和二维指纹图

Figure 4. Hirshfeld surface and two-dimensional finger patterns of C₂N₁₀ molecule at 40 GPa

下载: 全尺寸图片幻灯片

图 5 G₂ZT晶体在不同压强下的能带色散关系

Figure 5. Energy band dispersion of G₂ZT crystal at different pressures

下载: 全尺寸图片幻灯片

图 6 不同压强下G₂ZT晶体的态密度分布

Figure 6. Distributions of DOS of G₂ZT crystal at different pressures

下载: 全尺寸图片幻灯片

图 7 不同压强下G₂ZT晶体的介电函数随入射能量的变化

Figure 7. Dielectric functions of G₂ZT crystal versus incident energy at different pressures

下载: 全尺寸图片幻灯片

图 8 G₂ZT晶体在0和40 GPa压强下的反射率随入射能量的变化

Figure 8. Reflection index of G₂ZT crystal versus incident energy at 0 and 40 GPa

下载: 全尺寸图片幻灯片

图 9 G₂ZT晶体在0和40 GPa压强下的折射率随入射能量的变化

Figure 9. Refraction index of G₂ZT crystal versus incident energy at 0 and 40 GPa

下载: 全尺寸图片幻灯片

图 10 G₂ZT晶体在0和40 GPa压强下的吸收系数随入射能量的变化

Figure 10. Absorption coefficient of G₂ZT crystal versus incident energy at 0 and 40 GPa

下载: 全尺寸图片幻灯片

表 1 G₂ZT晶格参数的计算结果和实验结果比较^[2]

Table 1. Calculated crystal lattice parameters of G₂ZT compared with experimental data^[2]

Method	a/Å	b/Å	c/Å	$\alpha$ /(°)	$\,\beta$ /(°)	$\gamma$ /(°)	V/Å³
vdW-DF2	5.3733	6.6044	11.9931	102.36	91.14	109.42	390.1694
PBE-D2	5.3540	6.7920	11.7190	100.23	91.80	110.55	390.6893
Expt.^[2]	5.2619	6.6980	11.8840	102.05	90.80	109.96	383.3499
$\delta$ _vdW-DF2/%	2.10	−1.40	0.92	0.30	0.37	−0.49	1.80
$\delta$ _PBE-D2/%	1.70	1.40	−1.40	−1.80	1.10	0.53	1.90

下载: 导出CSV

表 2 三阶BM和Vinet物态方程拟合得到的G₂ZT晶体的体弹模量及其一阶导数

Table 2. Bulk moduli and their pressure-derivatives of G₂ZT crystal determined by third-order BM and Vinet equations of state

Method	Third-order BM		Vinet
Method	B₀/GPa	$B_0'$	B₀/GPa	$B_0'$
PBE-D2	19.69±0.80	5.50±0.22	18.84±0.63	6.06±0.16
vdW-DF2	27.58±0.76	4.85±0.14	26.48±0.66	5.39±0.13

下载: 导出CSV

表 3 不同压强下第一布里渊区高对称 $k$ 点在价带顶 $E_{\rm{v}}$ 和导带底 $E_{\rm{c}}$ 的特征能量

Table 3. Characteristic energy values for top of valence band $E_{\rm{v}}$ and bottom of conduction band $E_{\rm{c}}$ in highly symmetric $k$ points of the first Brillouin region at different pressures

Pressure/GPa	E_c/eV				E_v/eV
Pressure/GPa	R	Γ	X	M	R	Γ	X	M
0	2.03	2.25	2.03	2.03	0.00	−0.05	0.00	0.00
3	1.73	2.01	1.73	1.73	0.00	−0.04	0.00	0.00
10	1.33	1.66	1.31	1.33	0.00	−0.09	−0.03	0.00
40	0.93	1.46	0.81	0.91	−0.04	−0.25	−0.18	0.00

下载: 导出CSV

参考文献(28)

[1]	KLAPÖTKE T M, SABATÉ C M. Bistetrazoles: nitrogen-rich, high-performing, insensitive energetic compounds [J]. Chemistry of Materials, 2008, 20(11): 3629–3637. doi: 10.1021/cm703657k
[2]	KLAPÖTKE T M, SABATÉ 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.
[3]	KARAGHIOSOFF K, KLAPÖTKE T M, SABATÉ 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. doi: 10.1002/ejic.200800939
[4]	EBESPÄCHER M, KLAPÖTKE T M, SABATÉ C M. Nitrogen-rich alkali metal 5, 5′-hydrazinebistetrazolate salts: environmentally friendly compounds in pyrotechnic mixtures [J]. New Journal of Chemistry, 2009, 33(3): 517–527. doi: 10.1039/B818927G
[5]	DE LUCIA F C JR, 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
[6]	DREGER Z A, STASH A I, YU Z G, et al. High-pressure structural response of an insensitive energetic crystal: dihydroxylammonium 5, 5′-bistetrazole-1, 1′-diolate (TKX-50) [J]. The Journal of Physical Chemistry C, 2017, 121(10): 5761–5767. doi: 10.1021/acs.jpcc.7b00867
[7]	FANG X, STONE M, STENNETT C. Pulsed laser irradiation of a nanoparticles sensitized RDX crystal [J]. Combustion and Flame, 2020, 214: 387–393. doi: 10.1016/j.combustflame.2020.01.009
[8]	ISBELL R A, BREWSTER M Q. Optical properties of energetic materials: RDX, HMX, AP, NC/NG, and HTPB [J]. Propellants, Explosives, Pyrotechnics, 1998, 23(4): 218–224. doi: 10.1002/(SICI)1521-4087(199808)23:4<218::AID-PREP218>3.0.CO;2-A
[9]	ORDEJÓN P, ARTACHO E, SOLER J M. Self-consistent order-N density-functional calculations for very large systems [J]. Physical Review B, 1996, 53(16): R10441–R10444. doi: 10.1103/PhysRevB.53.R10441
[10]	SÁNCHEZ-PORTAL D, ORDEJÓN P, ARTACHO E, et al. Density-functional method for very large systems with LCAO basis sets [J]. International Journal of Quantum Chemistry, 1997, 65(5): 453–461. doi: 10.1002/(SICI)1097-461X(1997)65:5<453::AID-QUA9>3.0.CO;2-V
[11]	SOLER J M, ARTACHO E, GALE J D, et al. The SIESTA method for ab initio order-N materials simulation [J]. Journal of Physics: Condensed Matter, 2002, 14(11): 2745–2779. doi: 10.1088/0953-8984/14/11/302
[12]	DION M, RYDBERG H, SCHRÖDER E, et al. Van der Waals density functional for general geometries [J]. Physical Review Letters, 2004, 92(24): 246401. doi: 10.1103/PhysRevLett.92.246401
[13]	HAMADA I. Van der Waals density functional made accurate [J]. Physical Review B, 2014, 89(12): 121103. doi: 10.1103/PhysRevB.89.121103
[14]	GRIMME S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction [J]. Journal of Computational Chemistry, 2006, 27(15): 1787–1799. doi: 10.1002/jcc.20495
[15]	GRIMME S, EHRLICH S, GOERIGK L. Effect of the damping function in dispersion corrected density functional theory [J]. Journal of Computational Chemistry, 2011, 32(7): 1456–1465. doi: 10.1002/jcc.21759
[16]	TROULLIER N, MARTINS J L. Efficient pseudopotentials for plane-wave calculations [J]. Physical Review B, 1991, 43(3): 1993–2006. doi: 10.1103/PhysRevB.43.1993
[17]	SANKEY O F, NIKLEWSKI D J. Ab initio multicenter tight-binding model for molecular dynamics simulations and other applications in covalent systems [J]. Physical Review B, 1989, 40(6): 3979–3995. doi: 10.1103/PhysRevB.40.3979
[18]	JUNQUERA J, PAZ Ó, SÁNCHEZ-PORTAL D, et al. Numerical atomic orbitals for linear-scaling calculations [J]. Physical Review B, 2001, 64(23): 235111. doi: 10.1103/PhysRevB.64.235111
[19]	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
[20]	TRAN F, BLAHA P. Accurate band gaps of semiconductors and insulators with a semilocal exchange-correlation potential [J]. Physical Review Letters, 2009, 102(22): 226401. doi: 10.1103/PhysRevLett.102.226401
[21]	AMBROSCH-DRAXL C, SOFO J O. Linear optical properties of solids within the full-potential linearized augmented planewave method [J]. Computer Physics Communications, 2006, 175(1): 1–14. doi: 10.1016/j.cpc.2006.03.005
[22]	TOLL J S. Causality and the dispersion relation: logical foundations [J]. Physical Review, 1956, 104(6): 1760–1770. doi: 10.1103/PhysRev.104.1760
[23]	RIVAS-SILVA J F, BLAS M A, HOAT D M. Theoretical study of electronic and optical properties of antiferromagnetic β-MnS using the modified becke johnson (mBJ) potential [J]. Journal of Physics and Chemistry of Solids, 2017, 128: 310–315.
[24]	BIRCH F. Finite elastic strain of cubic crystals [J]. Physical Review, 1947, 71(11): 809–824. doi: 10.1103/PhysRev.71.809
[25]	VINET P, FERRANTE J, SMITH J R, et al. A universal equation of state for solids [J]. Journal of Physics C: Solid State Physics, 1986, 19(20): L467–L473. doi: 10.1088/0022-3719/19/20/001
[26]	MCKINNON J J, JAYATILAKA D, SPACKMAN M A. Towards quantitative analysis of intermolecular interactions with Hirshfeld surfaces [J]. Chemistry Communications, 2007(37): 3814–3816. doi: 10.1039/b704980c
[27]	SPACKMAN M A, JAYATILAKA D. Hirshfeld surface analysis [J]. CrystEngComm, 2009, 11(1): 19–32. doi: 10.1039/B818330A
[28]	SPACKMAN P R, TURNER M J, MCKINNON J J, et al. CrystalExplorer: a program for Hirshfeld surface analysis, visualization and quantitative analysis of molecular crystals [J]. Journal of Applied Crystallography, 2021, 54(3): 1006–1011. doi: 10.1107/S1600576721002910