Ying Liu‡
ab,
Huifang Du‡c,
Leiming Fang‡d,
Fei Sune,
Haipeng Sua,
Zhongxue Ge*a,
Wei Guo*c and
Jinlong Zhu*be
aXi'an Modern Chemistry Research Institute, Xi'an, 710065, China. E-mail: gzx204@sina.com
bDepartment of Physics, Shenzhen Engineering Research Center for Frontier Materials Synthesis at High Pressures, Southern University of Science and Technology, Shenzhen 518055, China. E-mail: zhujl@sustech.edu.cn
cBeijing Key Laboratory of Nanophotonics and Ultrafine Optoelectronic Systems, School of Physics, Beijing Institute of Technology, Beijing 100081, People's Republic of China. E-mail: weiguo7@bit.edu.cn
dKey Laboratory for Neutron Physics, Institute of Nuclear Physics and Chemistry, China Academy of Engineering Physics, Mianyang 621900, China
eCenter for High Pressure Science and Technology Advanced Research (HPSTAR), Beijing 100094, China
First published on 17th June 2021
High-energy-density materials (HEDMs) require new design rules collected from experimental and theoretical results and a proposed mechanism. One of the targeted systems is the nitrogen-rich compounds as precursors for possible polymeric nitrogen or its counterpart in a reasonable pressure range. 1H-tetrazole (CH2N4) with hydrogen bonds was studied under pressure by both diffraction and spectroscopy techniques. The observed crystal structure phase transition and hydrogen bond-assisted electronic structure anomaly were confirmed by first-principles calculation. The rearrangement of the hydrogen bonds under pressure elucidates the bonding interactions of the nitrogen-rich system in local 3D chemical environments, allowing the discovery and design of a feasible materials system to make new-generation high-energy materials.
In this study, we investigated the phase transition of 1H-tetrazole under high pressure at room temperature using in situ Raman, IR spectra, synchrotron X-ray diffraction (XRD) and neutron diffraction techniques. Combined with first-principles calculations, we propose a possible hydrogen-bond-assisted phase transition in the high-pressure counterpart, which has not been reported in previous studies. We uncover the novel electronic and stereo-chemical evolution of 1H-tetrazole (CH2N4), providing a fundamental understanding of the nitrogen-rich system and a potential pathway to find a targeted system for new-generation high-energy materials.
Fig. 1 Ambient and high-pressure crystal structures for 1H-tetrazole. The details are listed in Tables S1 & S2.† Type-O has N1–H1 along 4 o'clock, C5–H5 along 6 o'clock; type-A has N1–H1 along 7 o'clock, C5–H5 along 9 o'clock; type-B has N1–H1 along 1 o'clock, C5–H5 along 3 o'clock, respectively. The molecular orientation of type-A and type-B on the top is the up-view of the crystal structure. The orange dashed lines represent the hydrogen bonds and the translucent green dashed lines indicates that newly proposed hydrogen bonds will form at ∼10 GPa. |
As the pressure increases, the intralayer hydrogen-bonds of the ambient phase gradually change to the interlayer hydrogen close contacts of the high-pressure phase, where new hydrogen-bonds C5–H5⋯N2 and C5–H5⋯N3 form between the two layers. Thus, a pair of type-As and type-Bs molecules were coupled by hydrogen bonding interaction presented by orange dashed lines. Based on this configuration, the behavior of directional hydrogen bonding interaction and repulsion of the lone pair electrons on nitrogen in the type-A and type-B molecules under external compression leads to gradually decreasing of layer spacing, as well as sliding between interlayers, as demonstrated in Fig. 1. The sub-lattice of type-A heterocyclic ring shifts to the left relatively to sub-lattice of type-B as the increase of pressure, resulting in type-A-type-B molecules more overlapped, and an anticipated H–π interaction denoted by translucent green dashed lines in Fig. 1. Simultaneously, the entangled angle changes of P1 structure were plotted in Fig. 2. The angle γ is mainly related to the layered distance of the coupled type-As and type-Bs molecular planes. The increase of angle α and decrease of angle β both related to the relative position sliding in the direction perpendicular to the molecular planes. As the pressure increasing, these will further wrinkle the layers and increase the energy of crystal as the distortion increasing in Fig. 1, the schematic show of 8 GPa structure. Two runs of synchrotron X-ray diffraction (XRD) of higher pressure were conducted to further investigate the structure evolution, as plotted in Fig. S1 & S2.† By employing the rigid body model in the GSAS II software package, we conducted a sequential refinement with pressures of 4 GPa, 6 GPa and 8 GPa of the neutron data and the refined volume was plotted in Fig. S2† with X-ray results. The abnormal volume increase at ∼10 GPa as shown in Fig. S2† could be from the shear stresses in the DAC, which is exaggerated by the layered character of the crystal structure. A proposed new H–π interaction between type-As and type-Bs layers forms at pressure of ∼10 GPa, where the anomaly of crystal angle shows up in Fig. 2b and c. This new hydrogen bonding interaction, confirmed further by the spectroscopic data, will release part of the energy for staying in the same structure at high pressure. The detailed, refined parameters are listed in Tables S1 & S2.†
Fig. 2 (a) Neutron diffractions patterns of 1H-tetrazole as a function of pressure (the diffraction patterns of X-ray under pressure are plotted in Fig. S1†). The excluded region is the diffraction from the Al gasket; (b) and (c) angle parameter evolutions of the high-pressure 1H-tetrazole phase (the dot lines are guiding for the eyes). The anomaly at around 9 GPa is proposed to be related to the newly-formed hydrogen bonding interaction. |
The Raman and IR transmittance spectra collected at near ambient pressure are shown in Fig. S3† of the ESI and are consistent with previous reports.26,27 The detailed assignment of different crystal vibrational modes is provided in Table S3.†9,26,28 We also observed three additional Raman modes at 107, 130, and 159 cm−1 in the lattice-based vibrational region, which correspond to the intermolecular vibrations. Briefly, the readers should pay attention to the C–H vibrational mode at wavenumbers higher than 3000 cm−1 and N–N vibrational modes located between 1000 cm−1 and 1100 cm−1 for both Raman and IR, which show an apparent anomaly at ∼9 GPa.
In situ high-pressure Raman and IR spectra are plotted in Fig. 3a and b, and the Raman spectra in a particular spectral region are presented in Fig. S4 in the ESI.† As shown, upon compression to 4.0 GPa, new Raman peaks emerge at 180 cm−1, 956 cm−1, 1083 cm−1, and 1138 cm−1, respectively, which is more likely to be indicative of the appearance of Davydov splitting such as stated in Cui's work.26,29 However the result of non-hydrostatic compression can not be excluded. This type of splitting may be caused by enhanced interactions between non-equivalent molecules, as shown in the high-pressure crystal structure in Fig. 1. Similarly, the IR mode at 951.5 cm−1 observed at ambient conditions and corresponding to the in-plane bending mode δ(NC–N) disappears. Simultaneously, a new IR mode at 1258 cm−1 caused by the splitting of the ν(C5–N1) mode emerges, further indicating the relevance of the first phase transition to the reduced crystal structure symmetry. The dependencies of both the Raman and IR shift on pressure are plotted in Fig. 4.
As expected, most of the vibrational modes exhibit a substantial blue shift under increasing pressure, and only a few modes display a red shift. First, the Raman vibrational mode of δ(NC–N) at 947 cm−1 and the IR mode of δ(NN–N) at 937 cm−1 exhibit a red shift under pressure, which may indicate the disappearance of the original hydrogen bond accompanied by the gradual formation of the new interlayer hydrogen-bonds of C5–H5⋯N2 and C5–H5⋯N3 with increasing pressure, eventually resulting in the phase transition. From a static crystal structure framework perspective, the refined NC–N and NN–N bond length and the bend angle of the newly formed hydrogen bond support such an unusual red shift. This also implies that the newly formed hydrogen bond causes a decrease of the vibration energy, which accordingly stabilizes the molecule in the new high-pressure phase. Furthermore, the new Raman peak at 956 cm−1 that splits from the δ(NC–N) mode also shows a pressure-dependent red shift before vanishing. Additionally, the γ(C5–H5) mode at 906 cm−1 shows a sudden drop at 4 GPa, which might be connected with the change of spatial configuration and the conversion from the original hydrogen-bond C5–H5⋯N2 and C5–H5⋯N3 between neighboring molecules in the same layer, which are connected between the adjacent layers during the structural phase transition to the high-pressure phase.
Both Raman and IR spectroscopic data indicated that, upon a continuous increase of pressure, there is a transition at ∼9–10 GPa, while the crystal structure remains unchanged. The 9–10 GPa anomaly has some signature changes; both the Raman and IR spectra show that at wave number 3211 cm−1, the C–H stretching vibration splits into two modes along with the original ν(C5–H5) mode showing a temporary significant red shift. At the same time, the out-of-plane C5–H5 bending IR mode at 907 cm−1 and the C5–H5 bending Raman mode at 1085 cm−1 display a sudden drop. Meanwhile, the N4C5 vibrational mode at 1550 cm−1 disappears at a pressure above 10 GPa; the N4–N3 mode at 1075 cm−1 softens simultaneously; the IR mode of δ(NN–N) at 937 cm−1 exhibits a V-shape via pressure, and the N1–N2 vibration mode remains unchanged. In addition, the peaks at 1083 cm−1 and 1138 cm−1 appear when the phase transition displays a slight red shift. Several small peaks or shoulders appear near the C5–N1 stretching vibrations, and also between C5N4 and N2N3 stretching vibrations.
This comprehensive evolution indicates that the new H–π interaction formed and its surrounding electronic configuration changed with a further increase of pressure. Considering the interlayer slip via pressure, we speculate that at around 9–10 GPa the N1–H1⋯π hydrogen bond formed between the new type-As and type-Bs molecule pairs in the adjacent layer as forecasted in Fig. 1. Based on this speculation, the tetrazole heterocyclic ring of one type-B(A) molecule can be donors of two different type-A(B) molecules in the adjacent layer. With the new N1–H1⋯π bond formed, the charge centre of the donor deviates to the new linked adjacent molecule. Therefore, the bonds of C5–H5 and the ring vibration change accordingly, resulting in the observed anomalies in the Raman and IR spectra. According to the refined high-pressure angle parameters, the angle γ turned to increasing from decreasing at ∼9–10 GPa, which is caused by the change of interaction between the coupled molecular planes and consistent with the trends of several vibration modes, such as the stretching of N3–N4 and bending mode of C5–H5. Thus, the new hydrogen bonding interaction between the interlayer molecules is most likely due to the occurrence of an electronic phase transition at around 9–10 GPa.
Using DFT calculations to optimize the lattice obtained from our experiments, we found that the lattice parameters and bond lengths were abnormal near 9 GPa, consistent with the observed anomaly in our experimental results. We have calculated phonon dispersion relations at 8 GPa, shown as Fig. S5.† No imaginary phonon frequency is found in the whole Brillouin zone, which reveals that the structure of 1H-tetrazole may be stable under this pressure in theory. To observe the abnormality of its electronic structure, we have calculated its electronic band structure and density of states (DOS) at different pressures. Fig. S6† shows the electronic band at different pressures and the change of band gap. By observing the change in band gap, we found that an anomaly occurred near 9 GPa. Fig. 5 summarizes the variation of the DOS at 8 GPa, 9 GPa, and 10 GPa. Fig. 5a reveals that both the conduction band and the valence band reduce in energy from 8 GPa to 9 GPa, while they shift higher energy from 9 GPa to 10 GPa. We quantitatively calculated the band centers of the conduction band and valence band around 9 GPa, confirming that these energy shifts caused the structural anomaly. For the specific values, see Table S4.† The PDOS of C 2p N 2p and H 1s is shown in Fig. 5b. We also calculated the band centers of these atoms, as shown in Table S5.† The changes in the unit-cell, bond length, and electronic structures are all consistent, indicating that they are nonlinear with respect to the increasing pressure from 8 to 10 GPa.
To understand the structural anomaly near 9 GPa, we calculated the phonon density of states (PHDOS) of 1H-tetrazole, as shown in Fig. 6. The frequencies of 1H-tetrazole at Gamma at 8 GPa, 9 GPa, 10 GPa are listed in Table S6.† We observed that the most obvious change in PHDOS at pressures from 8 GPa to 10 GPa is the intensity of the peak at 860 cm−1 and the peak position at 2400 cm−1. Fig. 6b–d show that the vibrational modes from the N and H atoms are the main contribution from the low frequency at 860 cm−1 and the high frequency at 2400 cm−1, respectively. The vibrational intensity of the N atoms at 860 cm−1 increases from 8 GPa to 9 GPa, but decreases from 9 GPa to 10 GPa. At 2400 cm−1, the vibration of the H atom moves toward high frequency from 8 GPa to 9 GPa and then moves toward low frequency from 9 GPa to 10 GPa. This nonlinear change in the vibrational modes of the N and H atoms near 9 GPa may be due to the formation of new hydrogen bonds between the interlayer molecules, which is consistent with the experimental results.
Fig. 6 (a) The total and (b)–(d) projected phonon density of states (PHDOS) of 1H-tetrazole at 8 GPa, 9 GPa, and 10 GPa. |
The high-pressure Raman experiments were conducted through a HORIBA HR EVOLUTION Raman microscope with a liquid-nitrogen-cooled CCD detector and spectral resolution of about 1 cm−1. The exciting laser had a wavelength of 532 nm. High-pressure in situ IR spectra were collected in a Bruker VERTEX 80V FTIR spectrometer and HYPERION 2000 IR microscope, which was operated under a vacuum of <3 hPa. The IR beam was set to a spot size of 20 × 30 μm by a pair of knife-edge diaphragms. The IR spectra were recorded using a nitrogen-cooled broad-band mercury cadmium telluride detector. The spectral resolution was 4 cm−1. The Raman spectrum of the sample quenched from 50 GPa returned to the original state, indicating that the local chemical environment was still preserved, and the amorphization reported in ref. 26 could be due to different pressurization routes and history. For instance, the pressurization rate/step and hydrostatic degree could play important roles in the amorphization of small molecular crystals. However, this is beyond the scope of this work and will not be discussed here.
High-pressure angle-dispersive X-ray was used for monitoring the structural change of the 1H-tetrazole sample and performed at Cornell High Energy Synchrotron Source (CHESS) at room temperature.34 The X-ray wavelength was λ = 0.4859 Å, and the diffracted X-rays were detected by an imaging plate (Mar345). To elaborate the refined results, the second experimental run was conducted at the 4W2 High Pressure Station of Beijing Synchrotron Radiation Facility (BSRF), where the X-ray wavelength was λ = 0.6199 Å.
High-pressure neutron diffractions was collected at HPND beamline of CMRR. PE-cell with polycrystalline nano-diamond anvil and Al gasket was used for high pressure experiments at room temperature. The wavelength was λ = 1.5925 Å.
XRD patterns were integrated from the images with FIT2D and GSAS II software,35,36 and then refined with GSAS I and GSAS II37,38 for additional analysis. The lattice parameters were obtained by LeBail fitting of all the X-ray data. Since a P1 structure has 6 refinable parameters, plus one for intensity, and 3 more per 1H-tetrazole (using a rigid body model,37 we only refined 3 parameters (x,y,z) for the molecule's center of mass position). Additionally, peak shape, thermal displacements were also refined. The background was fitted in advance for reducing the refine parameters. In total, 15 parameters were refined. Considering the compressibility of C–C, as in graphite only 0.01–0.02 A changes in 10 GPa, which is within the error bar and resolution range, so the bond length in the rigid body was fixed as ambient ones and not refined for the three pressure data points.38
A DFT calculation based on our refined crystal structure was performed. The atomic positions were optimized at each fixed pressure. The geometric optimization and electronic property calculation were performed using the Vienna Ab Initio simulation package (VASP) program39 based on the density functional theory (DFT).40 The exchange–correlation function was described by the Perdew–Burke–Ernzerh (PBE) of the generalized gradient approximation (GGA).41 The C 2s22p2, N 2s22p3, and H 1s1 orbitals were explicitly included as valence electrons. The intermolecular van der Waals interactions42,43 (DFT-D3) were implemented in the structural optimization. A plane-wave basis set cutoff of 550 eV and a k-mesh grid of 5 × 7 × 7 for Brillouin zone sampling44 were chosen. The convergence criterion of the energy and force in the ionic relaxation was set as 10−5 eV and 0.02 eV Å−1 in the simulation. The phonon density of state was calculated using the supercell method with the PHONOPY code.45
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ra00522g |
‡ Y. L., H. D. and L. F. contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2021 |