Anisotropic structure deformation of ß-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine under high pressure: vibration spectra calculation and resolution based on AIMD simulation.

The phase transition of the ß-HMX crystal has been widely studied under high pressure, but the microscopic transition mechanism is not sufficiently understood. In this article, we perform a series of ab initio molecular dynamics simulations focusing on structure deformation and the corresponding vibration spectra resolution of ß-HMX at 0-40 GPa. Several typical pressure-induced phase transition processes are confirmed by analyzing the chemical bond, dihedral angle, charge transfer, and IR and Raman spectra. The corresponding relationship between molecular structure and spectral signal is constructed through the partial spectra calculations of special functional groups within the HMX molecule. The anisotropic effects of different groups on the initial structural phase transition are uncovered. The equatorial C-N and axial N-N bonds have the largest compression ratio as pressure increases, which is the intrinsic factor for the initiation of structure transformation. The C-N molecular ring plays an important role in the entire phase transition process. In addition, the phase transition of ß â†’ ζ is also closely related to the deformation of NO2, while that of ζ → ε is induced by the axial N-NO2 group. Regarding the higher-pressure phase transition, the synergetic effect of N-NO2, CH2 groups, and molecular rings becomes more considerable.

The Structure Properties of Carbon Materials Formed in 2,4,6-Triamino-1,3,5-Trinitrobenzene Detonation: A Theoretical Insight for Nucleation of Diamond-like Carbon.

The structure and properties of nano-carbon materials formed in explosives detonation are always a challenge, not only for the designing and manufacturing of these materials but also for clearly understanding the detonation performance of explosives. Herein, we study the dynamic evolution process of condensed-phase carbon involved in 2,4,6-Triamino-1,3,5-trinitrobenzene (TATB) detonation using the quantum-based molecular dynamics method. Various carbon structures such as, graphene-like, diamond-like, and "diaphite", are obtained under different pressures. The transition from a C sp2- to a sp3-hybrid, driven by the conversion of a hexatomic to a non-hexatomic ring, is detected under high pressure. A tightly bound nucleation mechanism for diamond-like carbon dominated by a graphene-like carbon layer is uncovered. The graphene-like layer is readily constructed at the early stage, which would connect with surrounding carbon atoms or fragments to form the tetrahedral structure, with a high fraction of sp3-hybridized carbon. After that, the deformed carbon layers further coalesce with each other by bonding between carbon atoms within the five-member ring, to form the diamond-like nucleus. The complex "diaphite" configuration is detected during the diamond-like carbon nucleation, which illustrates that the nucleation and growth of detonation nano-diamond would accompany the intergrowth of graphene-like layers.

Modifying a D-A-π-A-D HTM system for higher hole mobility by the meta-substitution strategy to weaken the electron-donating ability of the donor unit: a DFT study.

Tuning the electron-donating ability (EDA) of the donor units of hole transporting materials (HTMs) is an efficient strategy to modulate the optoelectronic properties of HTMs. Based on this strategy, we first theoretically investigated the effects of the EDA of donor units on D-A-π-A-D architectural HTMs. The results show that the enhanced EDA of the donor unit leads to larger hole reorganization energy and poorer molecular stability of HTMs. In contrast, meta-substitution of side groups is an effective strategy to reduce the EDA of the donor unit. We found that the application of the meta-substitution strategy in the D-A-π-A-D system not only successfully improves the molecular stability, but also achieves higher hole mobility by promoting the electronic coupling between the molecular dimers and decreasing the hole reorganization energies simultaneously. Studies on interfacial properties indicate that intermolecular coupling also synergistically enhances the interfacial charge extraction performance and reduces carrier recombination. In conclusion, by utilizing the meta-substitution strategy to reduce the EDA of donor units on D-A-π-A-D architectural HTMs, we successfully designed four superior performance HTMs mD1, mD2, mD3, and mD4.

Molecular Dynamics Study on the Reaction of RDX Molecule with Si Substrate.

RDX is widely used in various explosion situations, and there are many studies on its detonation performance, safety, preparation, etc. Research on preparation of ß-RDX is mainly conducted by experiments. In recent years, part of the research points to the use of substrate as a medium to produce ß-RDX faster. Based on this guidance, our work aims to theoretically solve the physical and chemical processes that RDX may experience in the production process through numerical simulation. In this work, molecular dynamics simulation is set up for the interaction between RDX and a Si clean surface and a Si hydroxyl saturated surface separately, and a higher precision simulation is set up to verify the reliability of the results. NCI analysis is also used to guess the possible phase transition mechanism in the simulation results. In the simulation process, a 7 × 7 Si clean surface, a 3 × 3 Si clean surface, and a 7 × 7 Si-OH surface are set, and each surface adsorbs one α-RDX. The semiempirical Gfn1-xtb method is used for the 7 × 7 surface, and the DFT method is used for the 3 × 3 surface. The calculation results confirmed by high-precision results show that RDX molecules will react with the dangling bonds on the Si surface. Three conformations of RDX were found on the hydroxyl saturated surface of Si. The isosurface generated by the NCI method is used to analyze the reasons for the formation of these conformations.

Enhanced Effect of an External Electric Field on NH3BH3 Dehydrogenation: an AIMD Study for Thermolysis.

How to improve the dehydrogenation properties of ammonia borane (AB, NH3BH3) is always a challenge for its practical application in hydrogen storage. In this study, we reveal the enhanced effect of an external electric field (E ext) on AB dehydrogenation by means of the ab initio molecular dynamics method. The molecular rotation induced by an electrostatic force can facilitate the formation of the H-N···B-H framework, which would aggregate into poly-BN species and further suppress the generation of the volatile byproducts. Meanwhile, the dihydrogen bond (N-Hδ+···Î´-H-B) is favorably formed under E ext, and the interaction between relevant H atoms is enhanced, leading to a faster H2 liberation. Correspondingly, the apparent activation energy for AB dissociation is greatly reduced from 18.42 to around 15 kcal·mol-1 with the application of an electric field, while that for H2 formation decreases from 20.4 to about 16 kcal·mol-1. In the whole process, the cleavage of the B-H bond is more favorable than that of the N-H bond, no matter whether the application of E ext. Our results give a deep insight into a positive effect of an electric field on AB dehydrogenation, which would provide an important inspiration for hydrogen storage in industry applications.

Anisotropic Reaction Properties for Different HMX/HTPB Composites: A Theoretical Study of Shock Decomposition.

Plastic-bonded explosives (PBXs) consisting of explosive grains and a polymer binder are commonly synthesized to improve mechanical properties and reduce sensitivity, but their intrinsic chemical behaviors while subjected to stress are not sufficiently understood yet. Here, we construct three composites of ß-HMX bonded with the HTPB binder to investigate the reaction characteristics under shock loading using the quantum-based molecular dynamics method. Six typical interactions between HMX and HTPB molecules are detected when the system is subjected to pressure. Although the initial electron structure is modified by the impurity states from HTPB, the metallization process for HMX does not significantly change. The shock decompositions of HMX/HTPB along the (100) and (010) surface are initiated by molecular ring dissociation and hydrogen transfer. The initial oxidations of C and H within HTPB possess advantages. As for the (001) surface, the dissociation is started with alkyl dehydrogenation oxidation, and a stronger hydrogen transfer from HTPB to HMX is detected during the following process. Furthermore, considerable fragment aggregation is observed, which mainly derives from the formation of new C-C and C-N bonds under high pressure. The effect of cluster evolution on the progression of the following reaction is further studied by analyzing the bonded structure and displacement rate.

Intermolecular Vibration Energy Transfer Process in Two CL-20-Based Cocrystals Theoretically Revealed by Two-Dimensional Infrared Spectra.

Inspired by the recent cocrystallization and theory of energetic materials, we theoretically investigated the intermolecular vibrational energy transfer process and the non-covalent intermolecular interactions between explosive compounds. The intermolecular interactions between 2,4,6-trinitrotoluene (TNT) and 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) and between 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) and CL-20 were studied using calculated two-dimensional infrared (2D IR) spectra and the independent gradient model based on the Hirshfeld partition (IGMH) method, respectively. Based on the comparison of the theoretical infrared spectra and optimized geometries with experimental results, the theoretical models can effectively reproduce the experimental geometries. By analyzing cross-peaks in the 2D IR spectra of TNT/CL-20, the intermolecular vibrational energy transfer process between TNT and CL-20 was calculated, and the conclusion was made that the vibrational energy transfer process between CL-20 and TNTII (TNTIII) is relatively slower than between CL-20 and TNTI. As the vibration energy transfer is the bridge of the intermolecular interactions, the weak intermolecular interactions were visualized using the IGMH method, and the results demonstrate that the intermolecular non-covalent interactions of TNT/CL-20 include van der Waals (vdW) interactions and hydrogen bonds, while the intermolecular non-covalent interactions of HMX/CL-20 are mainly comprised of vdW interactions. Further, we determined that the intermolecular interaction can stabilize the trigger bond in TNT/CL-20 and HMX/CL-20 based on Mayer bond order density, and stronger intermolecular interactions generally indicate lower impact sensitivity of energetic materials. We believe that the results obtained in this work are important for a better understanding of the cocrystal mechanism and its application in the field of energetic materials.

Revealing the Relationship between Electric Fields and the Conformation of Oxytocin Using Quasi-Static Amide-I Two-Dimensional Infrared Spectra.

It is reported that the cis/trans conformation change of the peptide hormone oxytocin plays an important role in its receptors and activation and the cis conformation does not lead to antagonistic activity. Motivated by recent experiments and theories, the quasi-static amide-I 2D IR spectra of oxytocin are investigated using DFT/B3LYP (D3)/6-31G (d, p) in combination with the isotope labeling method under different electric fields. The theoretical amide-I IR spectra and bond length of the disulfide bond are consistent with the experimental values, which indicates that the theoretical modes are reasonable. Our theoretical results demonstrate that the oxytocin conformation is transformed from the cis conformation to the trans conformation with the change of the direction of the electric field, which is confirmed by the distance of the backbone carbonyl oxygen of Cys6 and Pro7, the Ramachandran plot of Cys6 and Pro7, the dihedral angle of Cß-S-S-Cß, and the rmsd of the oxytocin backbone. Moreover, the trans conformation as the result of the turn in the vicinity of Pro7 has a tighter secondary spatial structure than the cis conformation, including stronger hydrogen bonds, longer γ-turn geometry involving five amino acids, and a more stable disulfide bond. Our work provides new insights into the relationship between the conformation, the activation of the peptide hormone oxytocin, and the electric fields.

Quasi-Static Two-Dimensional Infrared Spectra of the Carboxyhemoglobin Subsystem under Electric Fields: A Theoretical Study.

There is no doubt that electric fields of a specific frequency and intensity could excite certain vibrational modes of a macromolecule, which alters its mode coupling and conformation. Motivated by recent experiments and theories, we study the mode coupling between the Fe-CO mode and CO-stretch mode and vibration energy transfer among the active site and proteins in carboxyhemoglobin (HbCO) under different electric fields using the quasi-static two-dimensional infrared spectra. This study uses iron-porphyrin-imidazole-CO and two distal histidines in HbCO as the subsystem. The potential energy and dipole moment surfaces of the subsystem are calculated using an all-electron ab initio (B3LYP-D3(BJ)) method with the basis set Lanl2dz for the Fe atom and 6-31G(d,p) for C, H, O, and N atoms. Although the subsystem is reduced dimensionally, the anharmonic frequency and anharmonicity of the CO-stretch mode show excellent agreement with experimental values. We use the revealing noncovalent interaction method to confirm the hydrogen bond between the Hε atom of the His63 and the CO molecule. Our study confirms that the mode coupling between the Fe-CO mode and CO-stretch mode does not exist when the subsystem is free of electric field perturbation, which is coupled when the electric field is -0.5142 V/nm. In addition, with the increases of distance between the active site and the His92, there is no vibrational energy transfer between them when the electric field is 1.028 V/nm. We believe that our work could provide new ideas for increasing the dissociation efficiency of the Fe-CO bond and theoretical references for experimental research.

Towards the low-sensitive and high-energetic co-crystal explosive CL-20/TNT: from intermolecular interactions to structures and properties.

Employing molecular dynamic (MD) simulations and solid-state density functional theory (DFT), we carried out thorough studies to understand the interaction-structure-property interrelationship of the co-crystal explosive 1 : 1 CL-20 : TNT. Our results revealed that the co-crystallization of CL-20 and TNT molecules enhances the intermolecular binding forces, where the main driving force for the formation of the co-crystal CL-20/TNT comes from HO and CO interactions, while OO contributes to the co-crystal stabilization. Furthermore, we also used the concept of atoms in molecule (AIM) and the reduced density gradient (RDG) to describe the spatial arrangements and interactions of co-crystal compositions, which showed that although the adjoining TNT molecules possess two symmetry groups and the adjoining CL-20 molecules possess the same symmetry group, their interactions are not identical. These spatial arrangements provide a good reference to the formation of other co-crystals. When the obtained structural and detonation properties of the three crystals were compared, it was observed that the CL-20/TNT co-crystal achieved the desirable properties of explosives, i.e., low-sensitivity and high-energy, possessing the advantages of both CL-20 and TNT explosives.

Shock response of condensed-phase RDX: molecular dynamics simulations in conjunction with the MSST method.

We have performed molecular dynamics simulations in conjunction with the multiscale shock technique (MSST) to study the initial chemical processes of condensed-phase RDX under various shock velocities (8 km s-1, 10 km s-1 and 11 km s-1). A self-consistent charge density functional tight-binding (SCC-DFTB) method was used. We find that the N-NO2 bond dissociation is the primary pathway for RDX with the NO2 groups facing (group 1) the shock, whereas the C-N bond scission is the dominant primary channel for RDX with the NO2 groups facing away from (group 2) the shock. In addition, our results present that the NO2 groups facing away from the shock are rather inert to shock loading. Moreover, the reaction pathways of a single RDX molecule under the 11 km s-1 shock velocity have been mapped out in detail, NO2, NO, N2O, CO and N2 were the main products.

Detoxification of aflatoxins on prospective approach: effect on structural, mechanical, and optical properties under pressures.

Aflatoxins are sequential of derivatives of coumarin and dihydrofuran with similar chemical structures and well-known carcinogenic agent. Many studies performed to detoxify aflatoxins, but the result is not ideal. Therefore, we studied structural, infrared spectrum, mechanical, and optical properties of these compounds in the aim of perspective physics. Mulliken charge distributions and infrared spectral analysis performed to understand the structural difference between the basic types of aflatoxins. In addition, the effect of pressure, different polarized, and incident directions on their structural changes was determined. It is found that AFB1 is most stable structure among four basic types aflatoxins (AFB1, AFB2, AFG1, and AFG2), and IR spectra are analyzed to exhibit the difference on structures of them. The mechanical properties of AFB1 indicate that the structure of this toxin can be easily changed by pressure. The real [Formula: see text] and imaginary [Formula: see text] parts of the dielectric function, and the absorption coefficient [Formula: see text] and energy loss spectrum [Formula: see text] were also obtained under different polarized and incident directions. Furthermore, biological experiments needed to support the toxic level of AFB1 using optical technologies.

IR Spectra of Different O2-Content Hemoglobin from Computational Study: Promising Detector of Hemoglobin Variant in Medical Diagnosis.

IR spectra of heme and different O2-content hemoglobin were studied by the quantum computation method at the molecule level. IR spectra of heme and different O2-content hemoglobin were quantificationally characterized from 0 to 100 THz. The IR spectra of oxy-heme and de-oxy-heme are obviously different at the frequency regions of 9.08-9.48, 38.38-39.78, 50.46-50.82, and 89.04-91.00 THz. At 24.72 THz, there exists the absorption peak for oxy-heme, whereas there is not the absorption peak for de-oxy-heme. Whether the heme contains Fe-O-O bond or not has the great influence on its IR spectra and vibration intensities of functional groups in the mid-infrared area. The IR adsorption peak shape changes hardly for different O2-content hemoglobin. However, there exist three frequency regions corresponding to the large change of IR adsorption intensities for containing-O2 hemoglobin in comparison with de-oxy-hemoglobin, which are 11.08-15.93, 44.70-50.22, and 88.00-96.68 THz regions, respectively. The most differential values with IR intensity of different O2-content hemoglobin all exceed 1.0 × 104 L mol-1 cm-1. With the increase of oxygen content, the absorption peak appears in the high-frequency region for the containing-O2 hemoglobin in comparison with de-oxy-hemoglobin. The more the O2-content is, the greater the absorption peak is at the high-frequency region. The IR spectra of different O2-content hemoglobin are so obviously different in the mid-infrared region that it is very easy to distinguish the hemoglobin variant by means of IR spectra detector. IR spectra of hemoglobin from quantum computation can provide scientific basis and specific identification of hemoglobin variant resulting from different O2 contents in medical diagnosis.

Dynamic Responses and Initial Decomposition under Shock Loading: A DFTB Calculation Combined with MSST Method for ß-HMX with Molecular Vacancy.

Despite extensive efforts on studying the decomposition mechanism of HMX under extreme condition, an intrinsic understanding of mechanical and chemical response processes, inducing the initial chemical reaction, is not yet achieved. In this work, the microscopic dynamic response and initial decomposition of ß-HMX with (1 0 0) surface and molecular vacancy under shock condition, were explored by means of the self-consistent-charge density-functional tight-binding method (SCC-DFTB) in conjunction with multiscale shock technique (MSST). The evolutions of various bond lengths and charge transfers were analyzed to explore and understand the initial reaction mechanism of HMX. Our results discovered that the C-N bond close to major axes had less compression sensitivity and higher stretch activity. The charge was transferred mainly from the N-NO2 group along the minor axes and H atom to C atom during the early compression process. The first reaction of HMX primarily initiated with the fission of the molecular ring at the site of the C-N bond close to major axes. Further breaking of the molecular ring enhanced intermolecular interactions and promoted the cleavage of C-H and N-NO2 bonds. More significantly, the dynamic response behavior clearly depended on the angle between chemical bond and shock direction.

Thermal decomposition of solid phase nitromethane under various heating rates and target temperatures based on ab initio molecular dynamics simulations.

The Car-Parrinello molecular dynamics simulation was applied to study the thermal decomposition of solid phase nitromethane under gradual heating and fast annealing conditions. In gradual heating simulations, we found that, rather than C-N bond cleavage, intermolecular proton transfer is more likely to be the first reaction in the decomposition process. At high temperature, the first reaction in fast annealing simulation is intermolecular proton transfer leading to CH3NOOH and CH2NO2, whereas the initial chemical event at low temperature tends to be a unimolecular C-N bond cleavage, producing CH3 and NO2 fragments. It is the first time to date that the direct rupture of a C-N bond has been reported as the first reaction in solid phase nitromethane. In addition, the fast annealing simulations on a supercell at different temperatures are conducted to validate the effect of simulation cell size on initial reaction mechanisms. The results are in qualitative agreement with the simulations on a unit cell. By analyzing the time evolution of some molecules, we also found that the time of first water molecule formation is clearly sensitive to heating rates and target temperatures when the first reaction is an intermolecular proton transfer.

Pressure-induced metallization of condensed phase ß-HMX under shock loadings via molecular dynamics simulations in conjunction with multi-scale shock technique.

The electronic structure and initial decomposition in high explosive HMX under conditions of shock loading are examined. The simulation is performed using quantum molecular dynamics in conjunction with multi-scale shock technique (MSST). A self-consistent charge density-functional tight-binding (SCC-DFTB) method is adapted. The results show that the N-N-C angle has a drastic change under shock wave compression along lattice vector b at shock velocity 11 km/s, which is the main reason that leads to an insulator-to-metal transition for the HMX system. The metallization pressure (about 130 GPa) of condensed-phase HMX is predicted firstly. We also detect the formation of several key products of condensed-phase HMX decomposition, such as NO2, NO, N2, N2O, H2O, CO, and CO2, and all of them have been observed in previous experimental studies. Moreover, the initial decomposition products include H2 due to the C-H bond breaking as a primary reaction pathway at extreme condition, which presents a new insight into the initial decomposition mechanism of HMX under shock loading at the atomistic level.

Anisotropic responses and initial decomposition of condensed-phase ß-HMX under shock loadings via molecular dynamics simulations in conjunction with multiscale shock technique.

Molecular dynamics simulations in conjunction with multiscale shock technique (MSST) are performed to study the initial chemical processes and the anisotropy of shock sensitivity of the condensed-phase HMX under shock loadings applied along the a, b, and c lattice vectors. A self-consistent charge density-functional tight-binding (SCC-DFTB) method was employed. Our results show that there is a difference between lattice vector a (or c) and lattice vector b in the response to a shock wave velocity of 11 km/s, which is investigated through reaction temperature and relative sliding rate between adjacent slipping planes. The response along lattice vectors a and c are similar to each other, whose reaction temperature is up to 7000 K, but quite different along lattice vector b, whose reaction temperature is only up to 4000 K. When compared with shock wave propagation along the lattice vectors a (18 Å/ps) and c (21 Å/ps), the relative sliding rate between adjacent slipping planes along lattice vector b is only 0.2 Å/ps. Thus, the small relative sliding rate between adjacent slipping planes results in the temperature and energy under shock loading increasing at a slower rate, which is the main reason leading to less sensitivity under shock wave compression along lattice vector b. In addition, the C-H bond dissociation is the primary pathway for HMX decomposition in early stages under high shock loading from various directions. Compared with the observation for shock velocities V(imp) = 10 and 11 km/s, the homolytic cleavage of N-NO2 bond was obviously suppressed with increasing pressure.

Molecular dynamics study on the correlation between structure and sensitivity for defective RDX crystals and their PBXs.

Molecular dynamics simulation was applied to investigate the sensitivities of perfect and defective RDX (cyclotrimethylene trinitramine) crystals, as well as their PBXs (polymer-bonded explosives) with the polymeric binder F(2311), in the NPT (constant number of particles, constant pressure, constant temperature) ensemble using the COMPASS force field. Five kinds of defects-two dislocations, one vacancy, and two types of doping-were considered separately. The bond length distribution and the maximum (L (max)) and average (L (ave)) bond lengths of the N-NO(2) trigger bonds in RDX were obtained and their relationships to the sensitivities of RDX and PBXs are discussed. L (max) was found to be an important structural parameter for judging the relative sensitivity, and defects were observed to have little effect on the sensitivities of PBXs, due to the strong desensitizing effect of the polymer F(2311).

Initial decomposition of the condensed-phase ß-HMX under shock waves: molecular dynamics simulations.

We have performed quantum-based multiscale simulations to study the initial chemical processes of condensed-phase octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) under shock wave loading. A self-consistent charge density-functional tight-binding (SCC-DFTB) method was employed. The results show that the initial decomposition of shocked HMX is triggered by the N-NO(2) bond breaking under the low velocity impact (8 km/s). As the shock velocity increases (11 km/s), the homolytic cleavage of the N-NO(2) bond is suppressed under high pressure, the C-H bond dissociation becomes the primary pathway for HMX decomposition in its early stages. It is accompanied by a five-membered ring formation and hydrogen transfer from the CH(2) group to the -NO(2) group. Our simulations suggest that the initial chemical processes of shocked HMX are dependent on the impact velocity, which gain new insights into the initial decomposition mechanism of HMX upon shock loading at the atomistic level, and have important implications for understanding and development of energetic materials.

First-principles study of high-pressure behavior of solid beta-HMX.

A first-principles plane-wave method with an ultrasoft pseudopotential scheme in the framework of the generalized gradient approximation (GGA) was used to calculate the lattice parameters, bulk modulus and its pressure derivative, energy band structures, density of states, phonon density of states, thermodynamic properties, and absorption spectra of solid beta-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (beta-HMX). The current study is focused on the thermodynamics and electronic properties that were not reported previously. The bulk modulus and its pressure derivative are also consistent with experimental data and other theoretical results. From the results for the band gaps and density of states, it was found that beta-HMX has the tendency to become a semiconductor with increasing pressure. As the temperature increases, the heat capacity, enthalpy, product of temperature and entropy, and Debye temperature all increase, whereas the free energy decreases. The optical absorption coefficients shift to higher frequencies/energies with increasing pressure. The present study leads to a better understanding of how energetic materials respond to compression.