Detonation Performance of Insensitive Nitrogen-Rich Nitroenamine Energetic Materials Predicted from First-Principles Reactive Molecular Dynamics Simulations.

Because of the excellent combination of high detonation and low sensitivity properties of the 1,1-diamino-2,2-dinitroethylene (FOX-7) energetic material (EM), it is useful to explore new energetic derivatives that start with the FOX-7 structure. However, most such derivatives are highly sensitive, making them unsuitable for EM applications. An exception is the new nitroenamine EM, 1,1-diamino-2-tetrazole-2-nitroethene (FOX-7-T) (synthesized by replacing a nitro group with a tetrazole ring), which exhibits good stability. Unfortunately, FOX-7-T shows an unexpected much lower detonation performance than FOX-7, despite its higher nitrogen content. To achieve an atomistic understanding of the insensitivity and detonation performance of FOX-7 and FOX-7-T, we carried out reactive molecular dynamics (RxMD) using the ReaxFF reactive force field and combined quantum mechanics MD (QM-MD). We found that the functional group plays a significant role in the initial decomposition reaction. For FOX-7, the initial decomposition involves only simple hydrogen transfer reactions at high temperature, whereas for FOX-7-T, the initial reaction begins at much lower temperature with a tetrazole ring breaking to form N2, followed by many subsequent reactions. Our first-principles-based simulations predicted that FOX-7-T has 34% lower CJ pressure, 15% lower detonation velocity, and 45% lower CJ temperature than FOX-7. This is partly because a larger portion of the FOX-7-T mass gets trapped into condensed phase carbon clusters at the CJ point, suppressing generation of gaseous CO2 and N2 final products, leading to reduced energy delivery. Our findings suggest that the oxygen balance is an important factor to be considered in the design of the next generation of high-nitrogen-containing EMs.

Synthesis of a Porous [14]Annulene Graphene Nanoribbon and a Porous [30]Annulene Graphene Nanosheet on Metal Surfaces.

The electrical and mechanical properties of graphene-based materials can be tuned by the introduction of nanopores, which are sensitively related to the size, morphology, density, and location of nanopores. The synthesis of low-dimensional graphene nanostructures containing well-defined nonplanar nanopores has been challenging due to the intrinsic steric hindrance. Herein, we report the selective synthesis of one-dimensional (1D) graphene nanoribbons (GNRs) containing periodic nonplanar [14]annulene pores on Ag(111) and two-dimensional (2D) porous graphene nanosheet containing periodic nonplanar [30]annulene pores on Au(111), starting from a same precursor. The formation of distinct products on the two substrates originates from the different thermodynamics and kinetics of coupling reactions. The reaction mechanisms were confirmed by a series of control experiments, and the appropriate thermodynamic and kinetic parameters for optimizing the reaction pathways were proposed. In addition, the combined scanning tunneling spectroscopy (STS) and density functional theory (DFT) calculations revealed the electronic structures of porous graphene structures, demonstrating the impact of nonplanar pores on the π-conjugation of molecules.

Increasing Oxygen Balance Leads to Enhanced Performance in Environmentally Acceptable High-Energy Density Materials: Predictions from First-Principles Molecular Dynamics Simulations.

Environmental concerns have stimulated the development of green alternatives to environmentally pollutive nitramine compounds used for high-energy density materials (HEDMs). The excellent energetic properties of CL20 make it a promising candidate, but its negative oxygen balance limits its efficiency for industrial and military applications. We predict here that CL20-EO formed by introducing ether links into the CC bonds of the original CL20 structure to attain balanced CO2 and H2O production leads to improved performance while minimizing the formation of carbonaceous clusters and toxic gases. To test this concept, we predicted the detonation properties at the Chapman-Jouguet (CJ) state using reactive molecular dynamics simulations with the ReaxFF force field combined with quantum mechanics based moleculear dynamics. We predict that CL20-EO enhances energetic performance compared to CL20 with a 6.0% increase in the CJ pressure and a 1.1% increase in the detonation velocity, which we attribute to achieving the correct oxygen balance to produce fully oxidized gaseous products. After expansion to normal conditions from the CJ state, CL20-EO leads only to nontoxic fully oxidized gases instead of forming the carbonaceous clusters and toxic gases found with CL-20. Thus, CL20-EO is predicted to be environmentally green. These results indicate that oxygen balance plays an important role in both energy availability and end-product toxicity and that balanced CO2 and H2O production systems provide promising candidates for the next generation of environmentally acceptable alternatives to toxic HEDMs while also enhancing the detonation performance.

Chemisorption-Induced Formation of Biphenylene Dimer on Ag(111).

We report an example that demonstrates the clear interdependence between surface-supported reactions and molecular-adsorption configurations. Two biphenyl-based molecules with two and four bromine substituents, i.e., 2,2'-dibromobiphenyl (DBBP) and 2,2',6,6'-tetrabromo-1,1'-biphenyl (TBBP), show completely different reaction pathways on a Ag(111) surface, leading to the selective formation of dibenzo[e,l]pyrene and biphenylene dimer, respectively. By combining low-temperature scanning tunneling microscopy, synchrotron radiation photoemission spectroscopy, and density functional theory calculations, we unravel the underlying reaction mechanism. After debromination, a biradical biphenyl can be stabilized by surface Ag adatoms, while a four-radical biphenyl undergoes spontaneous intramolecular annulation due to its extreme instability on Ag(111). Such different chemisorption-induced precursor states between DBBP and TBBP consequently lead to different reaction pathways after further annealing. In addition, using bond-resolving scanning tunneling microscopy and scanning tunneling spectroscopy, we determine with atomic precision the bond-length alternation of the biphenylene dimer product, which contains 4-, 6-, and 8-membered rings. The 4-membered ring units turn out to be radialene structures.

Elucidating the Cooperative Roles of Water and Lewis Acid-Base Pairs in Cascade C-C Coupling and Self-Deoxygenation Reactions.

Water plays pivotal roles in tailoring reaction pathways in many important reactions, including cascade C-C bond formation and oxygen elimination. Herein, a kinetic study combined with complementary analyses (DRIFTS, isotopic study, 1H solid-state magic angle spinning nuclear magnetic resonance) and density functional theory (DFT) calculations are performed to elucidate the roles of water in cascade acetone-to-isobutene reactions on a Zn x Zr y O z mixed metal oxide with balanced Lewis acid-base pairs. Our results reveal that the reaction follows the acetone-diacetone alcohol-isobutene pathway. Isobutene is produced through an intramolecular rearrangement of the eight-membered ring intermediate formed via the adsorption of diacetone alcohol on the Lewis acid-base pairs in the presence of cofed water. OH adspecies, formed by the dissociative adsorption of water on the catalyst surface, were found to distort diacetone alcohol's hydroxyl functional group toward its carbonyl functional group and facilitate the intramolecular rearrangement of diacetone alcohol to form isobutene. In the absence of water, diacetone alcohol binds strongly to the Lewis acid site, e.g., at a Zr4+ site, via its carbonyl functional group, leading to its dramatic structural distortion and further dehydration reaction to form mesityl oxide as well as subsequent polymerization reactions and the formation of coke. The present results provide insights into the cooperative roles of water and Lewis acid-base pairs in catalytic upgrading of biomass to fuels and chemicals.

Dislocation-mediated shear amorphization in boron carbide.

The failure of superhard materials is often associated with stress-induced amorphization. However, the underlying mechanisms of the structural evolution remain largely unknown. Here, we report the experimental measurements of the onset of shear amorphization in single-crystal boron carbide by nanoindentation and transmission electron microscopy. We verified that rate-dependent loading discontinuity, i.e., pop-in, in nanoindentation load-displacement curves results from the formation of nanosized amorphous bands via shear amorphization. Stochastic analysis of the pop-in events reveals an exceptionally small activation volume, slow nucleation rate, and lower activation energy of the shear amorphization, suggesting that the high-pressure structural transition is activated and initiated by dislocation nucleation. This dislocation-mediated amorphization has important implications in understanding the failure mechanisms of superhard materials at stresses far below their theoretical strengths.

Thermal decomposition and diffusion of methane in clathrate hydrates from quantum mechanics simulations.

Clathrate hydrates are ice-like crystalline substances in which small gas molecules are trapped inside the polyhedral cavities of water molecules. They are of great importance in both scientific research and the petroleum industry because of their applications in modern energy and environmental technologies. To achieve an atomistic-level understanding of the diffusion and decomposition of trapped molecules in clathrate hydrate, we used methane hydrates (MHs) as the prototype system and examined the methane diffusion and decomposition mechanism by employing quantum mechanics (QM) and quantum mechanics molecular dynamics (QMD) simulations. Our QMD simulations illustrated that the initial decomposition reaction in MHs initiates from hydrogen transfer among water molecules and attacks by fragments of O and OH on CH4 molecules are responsible for the destruction of the methane molecules. Next, our QM simulations revealed that the methane molecule prefers to escape from the ice cage through the hexagonal face at low temperature. To suppress the methane diffusion, we demonstrated that the diffusion barrier is significantly enhanced by adding electron or hole carriers. This is because the extra electrons and holes enhance the electrostatic interaction between methane and water molecules, leading to an increased diffusion barrier. Thus, the clathrate hydrates could be stabilized by adding extra free electron or hole carriers.

Thermal Stability and Detonation Properties of Potassium 4,4'-Bis(dinitromethyl)-3,3'-azofurazanate, an Environmentally Friendly Energetic Three-Dimensional Metal-Organic Framework.

Environmentally acceptable alternatives to toxic lead-based primary explosives have become increasingly demanding for energetic materials (EMs) because of environmental concerns. Recent experiments suggested that energetic three-dimensional (3D) metal-organic frameworks (MOFs) are promising candidates for the next generation of environmentally friendly primary explosives. A new energetic 3D MOF, denoted as potassium 4,4'-bis(dinitromethyl)-3,3'-azofurazanate, was synthesized and suggested as an excellent candidate for green primary explosives. To achieve an atomistic-level understanding of the thermal stability and detonation properties of this material, we carried out quantum mechanics simulations to examine its initial decomposition mechanism and the Chapman-Jouguet state for sustainable detonation. We find that the initial decomposition reaction of potassium 4,4'-bis(dinitromethyl)-3,3'-azofurazanate is to break the C2N2O five-member ring in which K+ ions play a significant role in stabilizing the molecule structure, leading to an excellent thermal stability. Furthermore, this MOF system has a higher detonation velocity than that of lead azide, a comparable detonation pressure and temperature, and no toxic gases are produced at detonation. The combination of these detonation properties makes it a promising candidate for green EMs. Our results suggest that synthesizing 3D MOFs is an effective approach to develop environmentally acceptable alternatives to toxic EMs with enhanced thermal stability.

Grain Boundary Sliding and Amorphization are Responsible for the Reverse Hall-Petch Relation in Superhard Nanocrystalline Boron Carbide.

The recent observation of the reverse Hall-Petch relation in nanocrystalline ceramics offers a possible pathway to achieve enhanced ductility for traditional brittle ceramics via the nanosize effect, just as nanocrystalline metals and alloys. However, the underlying deformation mechanisms of nanocrystalline ceramics have not been well established. Here we combine reactive molecular dynamics (RMD) simulations and experimental transmission electron microscopy to determine the atomic level deformation mechanisms of nanocrystalline boron carbide (B_{4}C). We performed large-scale (up to ∼3 700 000 atoms) ReaxFF RMD simulations on finite shear deformation of three models of grain boundaries with grain sizes from 4.84 (135 050 atoms) to 14.64 nm (3 702 861 atoms). We found a reverse Hall-Petch relationship in nanocrystalline B_{4}C in which the deformation mechanism is dominated by the grain boundary (GB) sliding. This GB sliding leads to the amorphous band formation at predistorted icosahedral GB regions with initiation of cavitation within the amorphous bands. Our simulation results are validated by the experimental observations of an intergranular amorphous GB phase due to GBs sliding under indentation experiments. These theoretical and experimental results provide an atomistic explanation for the influence of GBs on the deformation behavior of nanocrystalline ceramics, explaining the reverse Hall-Petch relation.

Prediction of the Chapman-Jouguet chemical equilibrium state in a detonation wave from first principles based reactive molecular dynamics.

The combustion or detonation of reacting materials at high temperature and pressure can be characterized by the Chapman-Jouguet (CJ) state that describes the chemical equilibrium of the products at the end of the reaction zone of the detonation wave for sustained detonation. This provides the critical properties and product kinetics for input to macroscale continuum simulations of energetic materials. We propose the ReaxFF Reactive Dynamics to CJ point protocol (Rx2CJ) for predicting the CJ state parameters, providing the means to predict the performance of new materials prior to synthesis and characterization, allowing the simulation based design to be done in silico. Our Rx2CJ method is based on atomistic reactive molecular dynamics (RMD) using the QM-derived ReaxFF force field. We validate this method here by predicting the CJ point and detonation products for three typical energetic materials. We find good agreement between the predicted and experimental detonation velocities, indicating that this method can reliably predict the CJ state using modest levels of computation.