Publications by authors named "Guang-Fu Ji"

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 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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

How to improve the dehydrogenation properties of ammonia borane (AB, NHBH) is always a challenge for its practical application in hydrogen storage. In this study, we reveal the enhanced effect of an external electric field ( ) 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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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, 10 km s and 11 km s). A self-consistent charge density functional tight-binding (SCC-DFTB) method was used. We find that the N-NO bond dissociation is the primary pathway for RDX with the NO groups facing (group 1) the shock, whereas the C-N bond scission is the dominant primary channel for RDX with the NO groups facing away from (group 2) the shock.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

IR spectra of heme and different O-content hemoglobin were studied by the quantum computation method at the molecule level. IR spectra of heme and different O-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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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).

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

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.

View Article and Find Full Text PDF

The structural parameters, elastic constants and electronic structure of NbB(2) under pressure are investigated by using first-principles plane-wave pseudopotential density functional theory within the generalized gradient approximation (GGA). The obtained results are in agreement with the available theoretical data. It is found that the elastic constants and the Debye temperature of NbB(2) increase monotonically and the anisotropies weaken with pressure.

View Article and Find Full Text PDF