呋咱醚类化合物是俄罗斯科学院最新报道的一类化合物，具有能量高、塑性强、感度适中，热安定性及水稳定性较好的特点[1-8]。2012年西安近代化学研究所首次报道了含能呋咱环醚三呋咱并氧杂环庚三烯(TFO)[9]，其熔点为76.5~77.0 ℃，密度为1.935 g·cm-3，计算爆速为8646 m·s-1，撞击感度为24%，摩擦感度为4%，H50为72.4 cm。TFO具有熔点低、能量密度高、感度低、热稳定性好的特点，综合性能较好，有望用于熔铸炸药中的液相载体炸药组份，TFO与高能炸药1, 3, 5-三硝基-1, 3, 5-三氮杂环已烷(RDX)等进行混合浇铸，可有效提高混合炸药的密度，降低感度，增加热稳定性，也可用作固体推进剂中的含能增塑剂、氧化剂使用[9-13]。目前，报道TFO的合成以3, 4-双(3x-氨基呋咱-4x-基)氧化呋咱(DATF)为原料[9-13]，经还原、氧化及醚化3步反应合成目标化合物TFO。该合成方法存在如下问题：(1)总收率较低，仅为32.5%；(2)中间体3, 4-双(3′-氨基呋咱-4′-基)呋咱(BATF)、3, 4-双(3′-硝基呋咱-4′-基)呋咱(BNTF)仅限于小试合成，不易得；(3)氧化、醚化步骤需要进行萃取，操作复杂，难于工业化。针对文献[9-13]报道的上述TFO合成存在问题，本文自主设计了TFO新的合成路线及工艺条件，以已实现工业化的3, 4-双(3x-硝基呋咱-4x-基)氧化呋咱(DNTF)为原料，经醚化、还原2步反应合成目标化合物TFO，总收率为71.1%。该方法总收率高、原料易得、操作简便，具有潜在工业化前景。研究遴选最佳还原剂，探讨温度、料比及反应时间等关键影响还原反应因素，确定最佳工艺条件；采用B3LYP方法在6-31G(d, p)基组水平上对TFO的结构进行计算，得到其稳定的几何构型和键级。在振动分析的基础上求得体系在不同温度下的热力学性质，得到温度对热力学性能影响的关系式，结果表明TFO稳定性较好。
2 实验(材料与方法) 2.1 实验材料与仪器3, 4-双(3x-硝基呋咱-4x-基)氧化呋咱(DNTF)自制[14]；无水碳酸钠、亚磷酸三乙酯、氯化亚锡、盐酸、醋酸、乙腈等均为分析纯。
NEXUS 870型傅里叶变换红外光谱仪，美国Nicolet公司；AV 500型(500 MHz)超导核磁共振仪，瑞士Bruker公司；Vario EL Ⅲ型自动微量有机元素分析仪，德国Elementar公司；SYNAPT型UPLC –Q - TOFMS液质联用仪，美国Waters公司；ZF-Ⅱ型三用紫外仪，上海市安亭电子仪器厂；DSC-204差示扫描量热仪，德国Netzsch公司；GC-2010型高效液相色谱仪，日本岛津公司。
2.2 合成路线以3, 4-双(3x-硝基呋咱-4x-基)氧化呋咱(DNTF)为原料，经醚化、还原2步反应合成了TFO(见图 1)。
|
图 1 TFO合成路线 Fig.1 A novel synthesis route of TFO |
25 ℃下搅拌，将1 g (3.2 mmol) DNTF加入到3 mL无水乙腈中，再加入1.53 g(14.43 mmol)无水碳酸钠，加料完毕后，将反应液加热至80 ℃反应4 h后，过滤除掉固体，将反应液减压浓缩，加入6 mL质量百分比为70%的醋酸溶液，冰箱中冷冻一夜，过滤，干燥后得0.639 g BFFO，收率84.6%，纯度为99.1%。红外光谱：IR (KBr) ν：1 655, 1 623, 1 562, 1 543, 1 470, 1 384, 1 151 cm-1；13C NMR (DMSO-d6，125 MHz)，δ：160.50, 159.98, 144.39, 137.78, 135.52, 105.03；MS (EI) m/z (%):236(M+，48)；元素分析：(C6N6O5，%)计算值：C 30.51, N 35.59；实测值：C 30.87, N 35.99。
2.4 三呋咱并氧杂环庚三烯(TFO)的合成25 ℃下搅拌，将1 g (4.23 mmol) BFFO加入到5 mL无水乙腈中，溶液澄清后再加入1.4 g (8.4 mmol)亚磷酸三乙酯，加料完毕后，升温至80 ℃，保温反应18 h，反应完毕后，冷却至室温，加入反应助剂反应1 h后，倒入水中，加入10 g质量百分比为36.5%的盐酸，静置析出固体，过滤洗涤得0.779 g白色TFO固体，收率83.7%，纯度为99.4% HPLC法；见图 2。IR (KBr) ν：1 569, 1 481, 1 165, 1 087, 1 003, 985 cm-1；13C NMR (DMSO-d6，125 MHz)，δ：160.79, 144.01, 137.28；MS (EI) m/z (%): 220(M+, 45), 30(100)；元素分析：(C6N6O4，%)计算值：C 32.73, N 38.18;实测值：C 32.55，N 37.99。
|
图 2 TFO高效液相色谱图 Fig.2 HPLC spectrum of TFO |
差示扫描量热仪(DSC)操作条件：试样量0.800 mg，试样皿为铝盘，气氛为流动氮气，流量为50.0 mL·min-1，压力为2 MPa，采用10.0 ℃·min-1的升温速率从25升至400 ℃。
2.6 量子化学计算理论计算部分采用Gaussian 09软件包[a]中的密度泛函(DFT)理论方法。由Becke的3参数混合泛函(B3)和Lee-Yang-Parr参数作为相关泛函(LYP)组成的B3LYP方法，配合6-31G(d, p)基组，在含能材料领域已有广泛应用[17-22]。
3 实验结果与讨论 3.1 还原反应因素的探讨 3.1.1 还原剂采用4种还原剂对BFFO进行了还原反应研究，考察了还原剂对合成TFO收率的影响，结果见表 1。
|
表 1 还原剂对TFO收率的影响 Table 1 Effect of reductants on yield of TFO |
由表 1可看出：还原体系采用锡与盐酸、锌与醋酸均可得到产物TFO，但收率低；采用二氯化锡与冰醋酸作为还原体系，从液相色谱以及TLC监测显示反应副产物较少，反应率较高，但由于需要萃取，减压蒸馏等一系列复杂的操作，损失较大，以致收率较低，收率仅为56.5%；采用亚磷酸三乙酯的还原体系收率最高，可达到83.7%，因此较佳还原体系为亚磷酸三乙酯。
3.1.2 还原反应料比的影响还原温度为80 ℃、还原时间为18 h，考察了料比对还原反应收率的影响，结果见表 2。
|
|
表 2 料比对TFO收率的影响 Table 2 Effect of reductant ratios on yield of TFO |
由表 2可看出，还原体系的料比对反应收率影响较大，当还原剂与BFFO料比为1.5:1时，TLC监测BFFO反应不完全；当还原剂与BFFO料比为2:1，收率最高达到了83.7%，纯度为99.4%，继续增加还原剂的量，收率未见提高，但是增加的亚磷酸三乙酯会带来后处理的麻烦。
3.1.3 反应温度的影响亚磷酸三乙酯与BFFO摩尔比为2:1，还原时间为18 h，研究了还原温度对TFO收率的影响，结果见表 3。
|
|
表 3 反应温度对TFO收率的影响 Table 3 Effect of reaction temperature on yield of TFO |
由表 3可看出，反应温度对反应收率的影响较大，温度越高，收率越高，当还原温度为反应溶剂-乙腈的回流温度80 ℃时，反应收率最高达到83.7%，纯度达到了99.4%，因此最适宜还原温度应为80 ℃。
3.1.4 反应时间的影响亚磷酸三乙酯与BFFO摩尔比为2:1，还原温度为80 ℃，考察了还原反应时间对TFO收率的影响，结果见表 4。
|
|
表 4 反应时间对TFO收率的影响 Table 4 Effect of reaction time on yield of TFO |
从表 4可以看出，当还原时间为18 h时，TFO收率达到最高(83.7%)，继续延长还原反应时间，收率未见进一步提高，因此还原最佳反应时间为18 h。
3.2 热性能研究采用DSC法，研究了TFO的热性能，结果如图 3所示：
|
图 3 TFO的DSC曲线 Fig.3 DSC curves of TFO |
高压(2MPa)下TFO的DSC曲线(升温速率β =10 ℃·min-1)如图 3所示，呈现一个吸热峰和一个放热峰。吸热峰峰型尖锐，峰值温度为79.3 ℃，与TFO熔点为76.5~77.5 ℃基本相符，表明样品吸热后由固体熔融相变为液态；放热峰峰型较宽，放热峰程为343~406 ℃，峰程跨度较大，说明样品熔融为液体后相变为气态分解；峰值温度为377.4 ℃，说明样品热稳定性较好。
3.3 几何构型和键级优化后的TFO几何构型及原子编号见图 4，通过对其优化结构的振动频率计算，所得频率均为正值，表明所得构型为相对稳定构型。该构型属于Cs点群，唯一对称面为C(2)-C(3)键与O(10)和O(1)连线所形成的平面，而平面C2v构型有一个向图 5构型方向振动的虚频。
|
图 4 B3LYP/6-31G**优化后的TFO结构 Fig.4 The optimized geometry of TFO at B3LYP/6-31G** |
|
图 5 TFO的LUMO轨道示意图 Fig.5 LUMO of TFO |
一般而言，键级的大小决定了分子中键的强弱。TFO经B3LYP/6-31G(d, p)优化后的自然轨道(NBO)分析所得键级见表 5，可以看出：TFO结构中，N─O键及C─C键的键级比C─N键级小，相对C─N更容易断裂。
|
|
表 5 B3LYP/6-31G(d, p)优化后的TFO键级 Table 5 The bond order of TFO at B3LYP/6-31G(d, p) |
TFO分子中，各原子的Mulliken电荷分布见表 6，总体呈现出对称分布的趋势，呋咱环上与醚键O(10)原子相连的C(8)及C(9)原子带有较多的正电荷，这是由于O(10)原子电负性较强，且由于O(10)原子与呋咱环共轭，负电荷被分散，表明分子结构的稳定性较好。
|
|
表 6 TFO的原子电荷 Table 6 The atomic charges of TFO |
分子轨道理论认为，分子轨道对化合物活性影响巨大，因此其研究结果可以为确定分子活性部位以及探讨作用机制等提供重要信息。TFO的最低空轨道(LUMO) (图 5)和最高占有轨道(HOMO) (图 6)的能量分别为-3.256 7和-8.330 6 eV；前线轨道的能级差为5.073 9 eV。C(7)在LUMO中的比重较大，尽管其所带的正电荷比C(9)少，说明C(7)更容易发挥电子受体的作用；O(10)原子带的负电荷最多，在HOMO轨道中占的比重并不大，因此不容易发挥供电子作用，表明醚键不易断裂。
|
图 6 TFO的HOMO轨道示意图 Fig.6 HOMO of TFO |
经B3LYP/6-31G(d, p)几何优化后，求得TFO在273~1 000 K标准热力学函数与温度关系趋势图，见图 7。
|
图 7 TFO热力学性质与温度关系 Fig.7 The relationships between TFO thermodynamic properties and temperature |
TFO在273~1 000 K标准热力学函数与温度关系式如式(1)~(3)。其中Et, m为热能，Cp, m为热容，Sp, m为熵。
| $ E_{{\rm{t,m}}}^{\rm{ \mathsf{ θ} }} = 0.000\;1{T^2} + 0.118\;8T + 194.694$ | (1) |
| $ C_{{\text{p}},{\text{m}}}^{\rm{ \mathsf{ θ} }} = - 0.000\;3{T^2} + 0.622\;2T + 11.097 $ | (2) |
| $ S_{{\text{p}},{\text{m}}}^{\rm{ \mathsf{ θ} }} = - 0.000\;1{T^2} + 0.671\;9T + 236.62 $ | (3) |
可以看出，所有热力学函数值均随温度的升高而增加，这是由于在高温时振动加剧所致。此外，对Cp, m做T的二阶导数，结果为一负常数，表明Cp, m随温度的变化越来越慢。
4 结论(1) 设计了三呋咱并氧杂环庚三烯(TFO)便捷新合成方法，以易得的3, 4-双(3x-硝基呋咱-4x-基)氧化呋咱(DNTF)为原料，采用醚化、亚磷酸三乙酯还原的方法合成了TFO，总收率为71.1%，并采用红外光谱、核磁共振光谱、质谱及元素分析等确定了结构。
(2) 探讨了还原反应的关键影响因素，确定了最佳的反应条件为：还原剂是亚磷酸三乙酯，料比n(亚磷酸三乙酯):n(BFFO)为2:1，还原反应温度为80 ℃，还原反应时间为18 h。
(3) 通过量子化学计算表明：分子结构中O(10)原子带的负电荷最多，在HOMO轨道中占的比重并不大，因此不容易发挥供电子作用，醚键较为稳定；此外，O(10)原子具有较强电负性，加之O(10)原子与呋咱环的共轭作用，负电荷得到分散，表明TFO分子结构的稳定性较好。
| [1] |
SHEREMETEV A B. 3, 3'-Bis(1-fluoro-1, 1-dinitromethyl)difurazanyl ether[C]//Proceeding of the 29th International ICT Conference. Karlsruhe: 1998: 1-13.
|
| [2] |
SHEREMETEV A B, ALEKSANDROVA N S, MELNIKOVA T M, et al. Synthesis of difurazanyl ethers from 4, 4'-dinitroazoxy furazan[J]. Heteroatom Chemistry, 2000, 11(1): 48-56. DOI:10.1002/(SICI)1098-1071(2000)11:1<48::AID-HC8>3.0.CO;2-T |
| [3] |
范艳洁, 王伯周, 来蔚鹏, 等. 3, 3'-二氰基二呋咱基醚(FOF-2)的合成、表征及量子化学研究[J]. 有机化学, 2009, 29(4): 614-620. FAN Y J, WANG B Z, LAI W P, et al. Synthesis, characterization and quantum chemistry study on 3, 3'-dicyanodifurazanyl ether(FOF-2)[J]. Chinese Journal of Organic Chemistry, 2009, 29(4): 614-620. |
| [4] |
SHEREMETEV A B, KHARITONOVA O V, MELNIKOVA T M, et al. Synthesis of symmetrical difurazanyl ethers[J]. Mendeleev Communication, 1996, 4(6): 141-143. |
| [5] |
SHEREMETEV A B, KULAGINA V O, BATOG L V, et al. Furazan derivatives: High energetic materials from diaminofurazan[C]//Proceeding of the 22nd International Pyseminar. Colorado: [s.n.], 1996: 377-388.
|
| [6] |
李辉, 王伯周, 于倩倩, 等. 3, 3'-二(四唑-5-基)二呋咱基醚的合成与表征[J]. 含能材料, 2012, 20(1): 18-21. LI H, WANG B Z, YU Q Q, et al. Synthesis and characterization of 3, 3'-bis(tetrazol-5-y1)difurazanyl ether[J]. Chinese Journal of Energetic Materials, 2012, 20(1): 18-21. DOI:10.3969/j.issn.1006-9941.2012.01.005 |
| [7] |
翟连杰, 王伯周, 李亚南, 等. 3, 4-双(3-氰基呋咱-4-氧基)呋咱的晶体结构及热行为[J]. 火炸药学报, 2013, 36(5): 63-67. ZHAI L J, WANG B Z, LI Y N, et al. Crystal structure and thermal behavior of 3, 4-bis(3-cyanofurazan -4-oxy)furazan (FOF-12)[J]. Chinese Journal of Explosives &Propellants, 2013, 36(5): 63-67. DOI:10.3969/j.issn.1007-7812.2013.05.011 |
| [8] |
张俊林, 肖川, 翟连杰, 等. 多硝基氮杂稠环化合物合成及性能[J]. 有机化学, 2016, 36(6): 1197-1207. ZHANG J L, XIAO C, ZHAI L J, et al. Synthesis and properties of the fused aza-polynitrocyclic compounds[J]. Chinese Journal of Organic Chemistry, 2016, 36(6): 1197-1207. |
| [9] |
王锡杰, 王伯周, 贾思媛, 等. 三呋咱并氧杂环庚三烯(TFO)合成研究[J]. 含能材料, 2012, 20(6): 258-259. WANG X J, WANG B Z, JIA S Y, et al. Synthesis study on trifurazanooxacyclohetpatriene[J]. Chinese Journal of Energetic Materials, 2012, 20(6): 258-259. |
| [10] |
王彬, 周彦水, 吴敏杰, 等. 3-氨基-4-偕氨肟基呋咱及其含能衍生物合成研究进展[J]. 火炸药学报, 2018, 41(3): 213-222. WANG B, ZHOU Y S, WU M J, et al. Research progress on synthesis of 3-amino-4-aminoximiofurazan and its energetic derivatives[J]. Chinese Journal of Explosives & Propellants, 2018, 41(3): 213-222. |
| [11] |
WANG X J, XU K Z, SUN Q, et al. A novel insensitive energetic material trifurazano oxacyclo heptatriene(TFO):Synthesis and detonation properties[J]. Propellants, Explosives, Pyrotechnics, 2015, 40(1): 9-12. DOI:10.1002/prep.201400148 |
| [12] |
周彦水, 李建康, 黄新萍, 等. 3, 4-双(4'-氨基呋咱-3')氧化呋咱的合成及性能[J]. 火炸药学报, 2007, 30(1): 54-56. ZHOU Y S, LI J K, HUANG X P, et al. Synthesis and properties of 3, 4-bis (4'-aminofurazano -3')furoxan[J]. Chinese Journal of Explosives & Propellants, 2007, 30(1): 54-56. DOI:10.3969/j.issn.1007-7812.2007.01.015 |
| [13] |
贾思媛, 张海昊, 王伯周, 等. 3, 4-双(3'-氨基呋咱-4'-基)呋咱的合成与表征[J]. 含能材料, 2013, 21(3): 289-293. JIA S Y, ZHANG H H, WANG B Z, et al. Synthesis and characterization of 3, 4-bis-(3'-aminofurazal-4'-y1)-furazan(BATF)[J]. Chinese Journal of Energetic Materials, 2013, 21(3): 289-293. DOI:10.3969/j.issn.1006-9941.2013.03.002 |
| [14] |
周彦水, 王伯周, 李建康, 等. 3, 4-双(4'-硝基呋咱基-3')氧化呋咱合成、表征与性能研究[J]. 化学学报, 2011, 69(14): 1673-1680. ZHOU Y S, WANG B Z, LI J K, et al. Study on synthesis, characterization and properties of 3, 4-bis (4'-nitrofurazano-3'-yl)furoxan[J]. Acta Chimca Sinica, 2011, 69(14): 1673-1680. |
| [15] |
周彦水, 王伯周, 王锡杰, 等. 双呋咱并[3", 4"-d]氧杂环庚三烯的合成[J]. 含能材料, 2012, 20(1): 137-138. ZHOU Y S, WANG B Z, WANG X J, et al. Synthesis of bifurazano[3, 4-b:3', 4'-f] furoxano[3", 4"-d]oxacyclohetpatriene[J]. Chinese Journal of Energetic Materials, 2012, 20(1): 137-138. DOI:10.3969/j.issn.1006-9941.2012.01.032 |
| [16] |
ZHOU Y S, XU K Z, WANG B Z, et al. Synthesis, structure and thermalproperties of bifurazano[3, 4-b:3', 4'-f]furoxano[3", 4"-d] -oxacyclohetpatriene (BFFO)[J]. Bulletin of the Korean Chemical Society, 2012, 33(10): 3317-3320. DOI:10.5012/bkcs.2012.33.10.3317 |
| [17] |
POLITZER P J, MURRAY S, GRICE M E, et al. Calculation of heats of sublimation and solid phase heats of formation[J]. Molecular Physics, 1997, 91(5): 923-928. |
| [18] |
QING H Z, JIA H Z, XIU J Q, et al. Molecular design and property prediction of high density polynitro[3.3.3] -propellane -derivatized frameworks as potential high explosive[J]. The Journal of Physical Chemistry A, 2014, 118(45): 10857-10865. DOI:10.1021/jp509549q |
| [19] |
马海霞, 肖鹤鸣, 宋纪蓉, 等. 1, 2, 4-三唑-5-酮结构和性质的理论研究[J]. 含能材料, 2005, 13(3): 166-168. MA H X, XIAO H M, SONG J R, et al. Theoretical study on the structure and properties of 1, 2, 4-triazol-5-one[J]. Chinese Journal of Energetic Materials, 2005, 13(3): 166-168. DOI:10.3969/j.issn.1006-9941.2005.03.009 |
| [20] |
梁晓琴. 四唑衍生物结构及性质的理论研究[J]. 四川师范大学学报(自然科学版), 2008, 31(2): 219-223. LIANG X Q. Theoretical study on geometries and properties of tetrazole derivatives[J]. Journal of Sichuan Normal University (Natural Science), 2008, 31(2): 219-223. DOI:10.3969/j.issn.1001-8395.2008.02.024 |
| [21] |
来蔚鹏, 廉鹏, 王伯周, 等. 二硝基吡唑并吡唑性能的量子化学研究[J]. 计算机与应用化学, 2007, 24(8): 1025-1028. LAI W P, LIAN P, WANG B Z, et al. Study of the properties of 3, 6-dinitropyrazolo[4, 3-c]pyrazoles by quantum chemistry[J]. Computers and Applied Chemistry, 2007, 24(8): 1025-1028. DOI:10.3969/j.issn.1001-4160.2007.08.004 |
| [22] |
林智辉, 高莉, 李敏霞, 等. 几种呋咱类含能化合物的合成、热行为及理论爆轰性能预估[J]. 火炸药学报, 2014, 37(3): 6-10. LIN Z H, GAO L, LI M X, et al. Synthesis, crystal, thermal behaviour and prediction of theoretical detonation performance for some energetic compounds derived from furazan[J]. Chinese Journal of Explosives & Propellants, 2014, 37(3): 6-10. DOI:10.3969/j.issn.1007-7812.2014.03.002 |