Ultrafast transformation of graphite to diamond: An ab initio study of graphite under shock compression

Crystal-Structure Matches in Solid-Solid Phase Transitions

The exploration of solid-solid phase transition suffers from the uncertainty of how atoms in two crystal structures match. We devised a theoretical framework to describe and classify crystal-structure matches (CSM). Such description fully exploits the translational and rotational symmetries and is independent of the choice of supercells. This is enabled by the use of the Hermite normal form, an analog of reduced echelon form for integer matrices. With its help, exhausting all CSMs is made possible, which goes beyond the conventional optimization schemes. In an example study of the martensitic transformation of steel, our enumeration algorithm finds many candidate CSMs with lower strains than known mechanisms. Two long-sought CSMs accounting for the most commonly observed Kurdjumov-Sachs orientation relationship and the Nishiyama-Wassermann orientation relationship are unveiled. Given the comprehensiveness and efficiency, our enumeration scheme provide a promising strategy for solid-solid phase transition mechanism research.

Extreme Metastability of Diamond and its Transformation to the BC8 Post-Diamond Phase of Carbon

Diamond possesses exceptional physical properties due to its remarkably strong carbon-carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum-accurate multimillion atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure-temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through premelting followed by BC8 nucleation and growth in the metastable carbon liquid. We propose a double-shock compression pathway for BC8 synthesis, which is currently being explored in experiments at the National Ignition Facility.

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. This suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

Carbonization of a Molybdenum Substrate Surface and Nanoparticles by a One-Step Method of Femtosecond Laser Ablation in a Hexane Solution

Molybdenum carbides (MoC and Mo₂C) are being reported for various applications, for example, catalysts for sustainable energies, nonlinear materials for laser applications, protective coatings for improving tribological performance, and so on. A one-step method for simultaneously fabricating molybdenum monocarbide (MoC) nanoparticles (NPs) and MoC surfaces with a laser-induced periodic surface structure (LIPSS) was developed by using pulsed laser ablation of a molybdenum (Mo) substrate in hexane. Spherical NPs with an average diameter of 61 nm were observed by scanning electron microscopy. The X-ray diffraction pattern and electron diffraction (ED) pattern results indicate that a face-centered cubic MoC was successfully synthesized for the NPs and on the laser-irradiated area. Notably, the ED pattern suggests that the observed NPs are nanosized single crystals, and a carbon shell was observed on the surface of MoC NPs. The X-ray diffraction pattern of both MoC NPs and LIPSS surface indicates the formation of FCC MoC, agreeing with the results of ED. The results of X-ray photoelectron spectroscopy also showed the bonding energy attributed to Mo-C, and the sp²-sp³ transition was confirmed on the LIPSS surface. The results of Raman spectroscopy have also supported the formation of MoC and amorphous carbon structures. This simple synthesis method for MoC may provide new possibilities for preparing Mo _x C-based devices and nanomaterials, which may contribute to the development of catalytic, photonic, and tribological fields.

Study on the Polycrystalline Mechanism of Polycrystalline Diamond Synthesized from Graphite by Direct Detonation Method

In this paper, a polycrystalline diamond was synthesized by the direct detonation method using graphite as the carbon source. By comparing the numbers of the obtained diamond particles and the original graphite particles, it was found that when the graphite phase transformed into the polycrystalline diamond during the detonation process, a single graphite particle would form multiple diamond nuclei, and the nuclei would grow simultaneously to form polycrystals. Accordingly, a validation experiment was designed, which added different ratios of inert additives while keeping the ratio of graphite to hexogen (RDX) unchanged. It was found that increasing the ratio of inert additives within a certain range could increase the grain size of a polycrystalline diamond, which is consistent with the obtained polycrystalline mechanism.

Phase transition mechanism of hexagonal graphite to hexagonal and cubic diamond:ab initiosimulation

Ab initiomolecular dynamics simulations are used to elucidate the mechanism of the phase transition in shock experiments from hexagonal graphite (HG) to hexagonal diamond (HD) or to cubic diamond (CD). The transition from HG to HD is found to occur swiftly in very small time of 0.2 ps, with large cooperative displacements of all the atoms. We observe that alternate layers of atoms in HG slide in opposite directions by 1/6 along the ±[2, 1, 0], which is about 0.7 Å, while simultaneously puckering by about ±0.25 Å perpendicular to thea-bplane. The transition from HG to CD occurred with more complex cooperative displacements. In this case, six successive HG layers slide in pairs by 1/3 along [0, 1, 0], [-1, -1, 0] and [1, 0, 0], respectively along with the puckering as above. We have also performed calculations of the phonon spectrum in HG at high pressure, which reveal soft phonon modes that may facilitate the phase transition involving the sliding and puckering of the HG layers. We have further calculated the Gibbs free energy, including the vibrational energy and entropy, and derived the phase diagram between HG and CD phases.

Molecular Insight into the Deformation of Single Crystal Copper Loaded by High-Speed Shock Wave

Molecular dynamics simulations were performed to study the evolution of single crystal copper with and without a nanovoid (located at the middle of crystal with a diameter of ~2.9 nm) when loaded with shock waves of different velocities. The simulation results show that the average particle velocity of single crystal copper linearly relates to the velocity of the loaded shock wave for both the systems (crystal with and without a nanovoid). When loaded by the shock wave, the equilibrated temperature and pressure of the system with a nanovoid are found to be slightly larger than those of the system without the nanovoid, while the volume of the system with the nanovoid is found to be lower than that of the void-free system. The single crystal copper undergoes a phase transition from face-centered cubic (FCC) to hexagonal-close packed (HCP) and a dislocation structure forms around the nanovoid. The existence of a nanovoid can induce the rearrangement and deformation of the crystalline structure and eventually lead to the plastic deformation of the system. This work provides molecular-level insight into the effect of nanovoids on the shock plasticity of metals, which can aid in the ultimate application of the control of material structure damage in shock-wave propagation.

Application of First-Principles-Based Artificial Neural Network Potentials to Multiscale-Shock Dynamics Simulations on Solid Materials

The use of artificial neural network (ANN) potentials trained with first-principles calculations has emerged as a promising approach for molecular dynamics (MD) simulations encompassing large space and time scales while retaining first-principles accuracy. To date, however, the application of ANN-MD has been limited to near-equilibrium processes. Here we combine first-principles-trained ANN-MD with multiscale shock theory (MSST) to successfully describe far-from-equilibrium shock phenomena. Our ANN-MSST-MD approach describes shock-wave propagation in solids with first-principles accuracy but a 5000 times shorter computing time. Accordingly, ANN-MD-MSST was able to resolve fine, long-time elastic deformation at low shock speed, which was impossible with first-principles MD because of the high computational cost. This work thus lays a foundation of ANN-MD simulation to study a wide range of far-from-equilibrium processes.

Diamond-Graphene Composite Nanostructures

The search for new nanostructural topologies composed of elemental carbon is driven by technological opportunities as well as the need to understand the structure and evolution of carbon materials formed by planetary shock impact events and in laboratory syntheses. We describe two new families of diamond-graphene (diaphite) phases constructed from layered and bonded sp³ and sp² nanostructural units and provide a framework for classifying the members of this new class of materials. The nanocomposite structures are identified within both natural impact diamonds and laboratory-shocked samples and possess diffraction features that have previously been assigned to lonsdaleite and postgraphite phases. The diaphite nanocomposites represent a new class of high-performance carbon materials that are predicted to combine the superhard qualities of diamond with high fracture toughness and ductility enabled by the graphitic units and the atomically defined interfaces between the sp³- and sp²-bonded nanodomains.

Decompression-Induced Diamond Formation from Graphite Sheared under Pressure

Graphite is known to transform into diamond under dynamic compression or under combined high pressure and high temperature, either by a concerted mechanism or by a nucleation mechanism. However, these mechanisms fail to explain the recently reported discovery of diamond formation during ambient temperature compression combined with shear stress. Here we report a new transition pathway for graphite to diamond under compression combined with shear, based on results from both theoretical simulations and advanced experiments. In contrast to the known model for thermally activated diamond formation under pressure, the shear-induced diamond formation takes place during the decompression process via structural transitions. At a high pressure with large shear, graphite transforms into ultrastrong sp^{3} phases whose structures depend on the degree of shear stress. These metastable sp^{3} phases transform into either diamond or graphite upon decompression. Our results explain several recent experimental observations of low-temperature diamond formation. They also emphasize the importance of shear stress for diamond formation, providing new insight into the graphite-diamond transformation mechanism.

Structural studies of carbons by neutron and x-ray scattering

Carbon can have many different forms and the characterisation of structural features on a length scale of 1 Å to 10 μm is important in defining its physical and chemical properties for the various forms. The use of either electro-magnetic (x-ray) or particle (neutron) beams plays an important role in determining these characteristics. In this paper, we review the various techniques that are used to determine the structural features by experimental means and how the data are processed to give the required information in a suitable form for detailed analysis by computer simulation. Diffraction methods are used for studies of the atomic arrangement and small-angle scattering techniques are used for studies of microporosity in the sample materials. The experimental data obtained from a wide range of different carbon materials are considered and how these results can be used as a basis for modelling the structures in a quantitative manner is also considered. This information underpins their use as active components in a wide range of functional materials.

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⁻¹, 10 km s⁻¹ and 11 km s⁻¹).

Transformation of shock-compressed graphite to hexagonal diamond in nanoseconds

The graphite-to-diamond transformation under shock compression has been of broad scientific interest since 1961. The formation of hexagonal diamond (HD) is of particular interest because it is expected to be harder than cubic diamond and due to its use in terrestrial sciences as a marker at meteorite impact sites. However, the formation of diamond having a fully hexagonal structure continues to be questioned and remains unresolved. Using real-time (nanosecond), in situ x-ray diffraction measurements, we show unequivocally that highly oriented pyrolytic graphite, shock-compressed along the c axis to 50 GPa, transforms to highly oriented elastically strained HD with the (100)HD plane parallel to the graphite basal plane. These findings contradict recent molecular dynamics simulation results for the shock-induced graphite-to-diamond transformation and provide a benchmark for future theoretical simulations. Additionally, our results show that an earlier report of HD forming only above 170 GPa for shocked pyrolytic graphite may lead to incorrect interpretations of meteorite impact events.

A metastable phase of shocked bulk single crystal copper: an atomistic simulation study

Structural phase transformation in bulk single crystal Cu in different orientation under shock loading of different intensities has been investigated in this article. Atomistic simulations, such as, classical molecular dynamics using embedded atom method (EAM) interatomic potential and ab-initio based molecular dynamics simulations, have been carried out to demonstrate FCC-to-BCT phase transformation under shock loading of 〈100〉 oriented bulk single crystal copper. Simulated x-ray diffraction patterns have been utilized to confirm the structural phase transformation before shock-induced melting in Cu(100).

Nanosecond formation of diamond and lonsdaleite by shock compression of graphite

The shock-induced transition from graphite to diamond has been of great scientific and technological interest since the discovery of microscopic diamonds in remnants of explosively driven graphite. Furthermore, shock synthesis of diamond and lonsdaleite, a speculative hexagonal carbon polymorph with unique hardness, is expected to happen during violent meteor impacts. Here, we show unprecedented in situ X-ray diffraction measurements of diamond formation on nanosecond timescales by shock compression of pyrolytic as well as polycrystalline graphite to pressures from 19 GPa up to 228 GPa. While we observe the transition to diamond starting at 50 GPa for both pyrolytic and polycrystalline graphite, we also record the direct formation of lonsdaleite above 170 GPa for pyrolytic samples only. Our experiment provides new insights into the processes of the shock-induced transition from graphite to diamond and uniquely resolves the dynamics that explain the main natural occurrence of the lonsdaleite crystal structure being close to meteor impact sites.

Shock melting method to determine melting curve by molecular dynamics: Cu, Pd, and Al

A melting simulation method, the shock melting (SM) method, is proposed and proved to be able to determine the melting curves of materials accurately and efficiently. The SM method, which is based on the multi-scale shock technique, determines melting curves by preheating and/or prepressurizing materials before shock. This strategy was extensively verified using both classical and ab initio molecular dynamics (MD). First, the SM method yielded the same satisfactory melting curve of Cu with only 360 atoms using classical MD, compared to the results from the Z-method and the two-phase coexistence method. Then, it also produced a satisfactory melting curve of Pd with only 756 atoms. Finally, the SM method combined with ab initio MD cheaply achieved a good melting curve of Al with only 180 atoms, which agrees well with the experimental data and the calculated results from other methods. It turned out that the SM method is an alternative efficient method for calculating the melting curves of materials.

Synthesis of diamond-like phase from graphite by ultrafast laser driven dynamical compression

Rapid variations of the environmental energy caused by ultrashort laser pulses have induced phase transitions in carbon allotropes, therefore bringing the promise of revealing new carbon phases. Here, by exposing polycrystalline graphite to 25 fs laser pulses at 4 J/cm² fluence under standard air atmosphere, we demonstrated the synthesis of translucent micrometer-sized structures carrying diamond-like and onion-like carbon phases. Texturized domains of the diamond phase were also identified. Concerning different synthesized carbon forms, pulse superposition and singularities of the thermodynamical process, we pinpoint the synthesis mechanism by the laser-induced subsequent products energetically evolving to attain the diamond-like phase.

Two-step nucleation mechanism in solid-solid phase transitions

The microscopic kinetics of ubiquitous solid-solid phase transitions remain poorly understood. Here, by using single-particle-resolution video microscopy of colloidal films of diameter-tunable microspheres, we show that transitions between square and triangular lattices occur via a two-step diffusive nucleation pathway involving liquid nuclei. The nucleation pathway is favoured over the direct one-step nucleation because the energy of the solid/liquid interface is lower than that between solid phases. We also observed that nucleation precursors are particle-swapping loops rather than newly generated structural defects, and that coherent and incoherent facets of the evolving nuclei exhibit different energies and growth rates that can markedly alter the nucleation kinetics. Our findings suggest that an intermediate liquid should exist in the nucleation processes of solid-solid transitions of most metals and alloys, and provide guidance for better control of the kinetics of the transition and for future refinements of solid-solid transition theory.

Graphene and nano-diamond synthesis in expansions of molten liquid carbon

Despite their widespread use in high-pressure experiments, little is known about the physical and chemical properties of carbon-containing materials as they expand and cool to ambient conditions. As a result, interpretation of experiments can rely on use of unconstrained models with poor accuracy for the ensuing equation of state properties and final chemical products. To this end, we use quantum simulations to study the free expansion and cooling of carbon from metallic liquid states achieved during shock compression. Expansions from three different sets of shock conditions yielded of a variety of chain and ring structures. We then quantify the relative amounts of graphite-like and diamond-like particles formed during cooling and equilibration. We observe that for all cases, graphene sheets are the majority product formed with more extreme initial conditions producing increasingly larger amounts of diamond particles. Our results can address key needs for future meso-scale models of experiments, where knowledge of material properties and chemical end products can have a pronounced effect on interpreting experimental observables.

Mechanism for direct graphite-to-diamond phase transition

Using classical molecular dynamics with a more reliable reactive LCBOPII potential, we have performed a detailed study on the direct graphite-to-diamond phase transition. Our results reveal a new so-called "wave-like buckling and slipping" mechanism, which controls the transformation from hexagonal graphite to cubic diamond. Based on this mechanism, we have explained how polycrystalline cubic diamond is converted from hexagonal graphite, and demonstrated that the initial interlayer distance of compressed hexagonal graphite play a key role to determine the grain size of cubic diamond. These results can broaden our understanding of the high pressure graphite-to-diamond phase transition.

A density functional tight binding model with an extended basis set and three-body repulsion for hydrogen under extreme thermodynamic conditions

We present a new DFTB-p3b density functional tight binding model for hydrogen at extremely high pressures and temperatures, which includes a polarizable basis set (p) and a three-body environmentally dependent repulsive potential (3b). We find that use of an extended basis set is necessary under dissociated liquid conditions to account for the substantial p-orbital character of the electronic states around the Fermi energy. The repulsive energy is determined through comparison to cold curve pressures computed from density functional theory (DFT) for the hexagonal close-packed solid, as well as pressures from thermally equilibrated DFT-MD simulations of the liquid phase. In particular, we observe improved agreement in our DFTB-p3b model with previous theoretical and experimental results for the shock Hugoniot of hydrogen up to 100 GPa and 25000 K, compared to a standard DFTB model using pairwise interactions and an s-orbital basis set, only. The DFTB-p3b approach discussed here provides a general method to extend the DFTB method for a wide variety of materials over a significantly larger range of thermodynamic conditions than previously possible.

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.

Prebiotic chemistry within a simple impacting icy mixture

We present results of prebiotic organic synthesis in shock compressed mixtures of simple ices from quantum molecular dynamics (MD) simulations extended to close to equilibrium time scales. Given the likelihood of an inhospitable prebiotic atmosphere on early Earth, it is possible that impact processes of comets or other icy bodies were a source of prebiotic chemical compounds on the primitive planet. We observe that moderate shock pressures and temperatures within a CO2-rich icy mixture (36 GPa and 2800 K) produce a number of nitrogen containing heterocycles, which dissociate to form functionalized aromatic hydrocarbons upon expansion and cooling to ambient conditions. In contrast, higher shock conditions (48-60 GPa, 3700-4800 K) resulted in the synthesis of long carbon-chain molecules, CH4, and formaldehyde. All shock compression simulations at these conditions have produced significant quantities of simple C-N bonded compounds such as HCN, HNC, and HNCO upon expansion and cooling to ambient conditions. Our results elucidate a mechanism for impact synthesis of prebiotic molecules at realistic impact conditions that is independent of external constraints such as the presence of a catalyst, illuminating UV radiation, or pre-existing conditions on a planet.

Optimal calculation of the pair correlation function for an orthorhombic system

We present a computational method to calculate arbitrary pair correlation functions of an orthorhombic system in the most efficient way. The algorithm is demonstrated by the calculation of the radial distribution function of shock compressed liquid hydrogen.

Nucleation mechanism for the direct graphite-to-diamond phase transition

Graphite and diamond have comparable free energies, yet forming diamond from graphite in the absence of a catalyst requires pressures that are significantly higher than those at equilibrium coexistence. At lower temperatures, the formation of the metastable hexagonal polymorph of diamond is favoured instead of the more stable cubic diamond. These phenomena cannot be explained by the concerted mechanism suggested in previous theoretical studies. Using an ab initio quality neural-network potential, we carried out a large-scale study of the graphite-to-diamond transition assuming that it occurs through nucleation. The nucleation mechanism accounts for the observed phenomenology and reveals its microscopic origins. We demonstrate that the large lattice distortions that accompany the formation of diamond nuclei inhibit the phase transition at low pressure, and direct it towards the hexagonal diamond phase at higher pressure. The proposed nucleation mechanism should improve our understanding of structural transformations in a wide range of carbon-based materials.

Anomalous sound propagation and slow kinetics in dynamically compressed amorphous carbon

We have performed molecular-dynamics simulations of dynamic compression waves propagating through amorphous carbon using the Tersoff potential and find that a variety of dynamic compression features appear for two different initial densities. These features include steady elastic shocks, steady chemically reactive shocks, unsteady elastic waves, and unsteady chemically reactive waves. We show how these features can be distinguished by analyzing time-dependent propagation speeds, time-dependent sound speeds, and comparison to multiscale shock technique (MSST) simulations. Understanding such features is a key challenge in quasi-isentropic experiments involving phase transformations. In addition to direct simulations of dynamic compression, we employ the MSST and find agreement with the direct method for this system for the shocks observed. We show how the MSST can be extended to include explicit material viscosity and demonstrate on an amorphous Lennard-Jones system.