Nanosecond formation of diamond and lonsdaleite by shock compression of graphite

Accuracy, transferability, and computational efficiency of interatomic potentials for simulations of carbon under extreme conditions

Large-scale atomistic molecular dynamics (MD) simulations provide an exceptional opportunity to advance the fundamental understanding of carbon under extreme conditions of high pressures and temperatures. However, the fidelity of these simulations depends heavily on the accuracy of classical interatomic potentials governing the dynamics of many-atom systems. This study critically assesses several popular empirical potentials for carbon, as well as machine learning interatomic potentials (MLIPs), in their ability to simulate a range of physical properties at high pressures and temperatures, including the diamond equation of state, its melting line, shock Hugoniot, uniaxial compressions, and the structure of liquid carbon. Empirical potentials fail to accurately predict the behavior of carbon under high pressure-temperature conditions. In contrast, MLIPs demonstrate quantum accuracy, with Spectral Neighbor Analysis Potential (SNAP) and atomic cluster expansion (ACE) being the most accurate in reproducing the density functional theory results. ACE displays remarkable transferability despite not being specifically trained for extreme conditions. Furthermore, ACE and SNAP exhibit superior computational performance on graphics processing unit-based systems in billion atom MD simulations, with SNAP emerging as the fastest. In addition to offering practical guidance in selecting an interatomic potential with a fine balance of accuracy, transferability, and computational efficiency, this work also highlights transformative opportunities for groundbreaking scientific discoveries facilitated by quantum-accurate MD simulations with MLIPs on emerging exascale supercomputers.

X-ray Thomson scattering absolute intensity from the f-sum rule in the imaginary-time domain

We present a formally exact and simulation-free approach for the normalization of X-ray Thomson scattering (XRTS) spectra based on the f-sum rule of the imaginary-time correlation function (ITCF). Our method works for any degree of collectivity, over a broad range of temperatures, and is applicable even in nonequilibrium situations. In addition to giving us model-free access to electronic correlations, this new approach opens up the intriguing possibility to extract a plethora of physical properties from the ITCF based on XRTS experiments.

Ultrafast Heating-Induced Suppression of d-Band Dominance in the Electronic Excitation Spectrum of Cuprum

The combination of isochoric heating of solids by free-electron lasers (FELs) and in situ diagnostics by X-ray Thomson scattering (XRTS) allows for measurements of material properties at warm dense matter (WDM) conditions relevant for astrophysics, inertial confinement fusion, and materials science. In the case of metals, the FEL beam pumps energy directly into electrons with the lattice structure of ions being nearly unaffected. This leads to a unique transient state that gives rise to a set of interesting physical effects, which can serve as a reliable testing platform for WDM theories. In this work, we present extensive linear-response time-dependent density functional theory (TDDFT) results for the electronic dynamic structure factor of isochorically heated copper with a face-centered cubic lattice. At ambient conditions, the plasmon is heavily damped due to the presence of d-band excitations, and its position is independent of the wavenumber. In contrast, the plasmon feature starts to dominate the excitation spectrum and has a Bohm-Gross-type plasmon dispersion for temperatures T ≥ 4 eV, where the quasi-free electrons in the interstitial region are in the WDM regime. In addition, we analyze the thermal changes in the d-band excitations and outline the possibility to use future XRTS measurements of isochorically heated copper as a controlled testbed for WDM theories.

Release dynamics of nanodiamonds created by laser-driven shock-compression of polyethylene terephthalate

Laser-driven dynamic compression experiments of plastic materials have found surprisingly fast formation of nanodiamonds (ND) via X-ray probing. This mechanism is relevant for planetary models, but could also open efficient synthesis routes for tailored NDs. We investigate the release mechanics of compressed NDs by molecular dynamics simulation of the isotropic expansion of finite size diamond from different P-T states. Analysing the structural integrity along different release paths via molecular dynamic simulations, we found substantial disintegration rates upon shock release, increasing with the on-Hugnoiot shock temperature. We also find that recrystallization can occur after the expansion and hence during the release, depending on subsequent cooling mechanisms. Our study suggests higher ND recovery rates from off-Hugoniot states, e.g., via double-shocks, due to faster cooling. Laser-driven shock compression experiments of polyethylene terephthalate (PET) samples with in situ X-ray probing at the simulated conditions found diamond signal that persists up to 11 ns after breakout. In the diffraction pattern, we observed peak shifts, which we attribute to thermal expansion of the NDs and thus a total release of pressure, which indicates the stability of the released NDs.

A comparative review of time-resolved x-ray and electron scattering to probe structural dynamics

The structure of molecules, particularly the dynamic changes in structure, plays an essential role in understanding physical and chemical phenomena. Time-resolved (TR) scattering techniques serve as crucial experimental tools for studying structural dynamics, offering direct sensitivity to molecular structures through scattering signals. Over the past decade, the advent of x-ray free-electron lasers (XFELs) and mega-electron-volt ultrafast electron diffraction (MeV-UED) facilities has ushered TR scattering experiments into a new era, garnering significant attention. In this review, we delve into the basic principles of TR scattering experiments, especially focusing on those that employ x-rays and electrons. We highlight the variations in experimental conditions when employing x-rays vs electrons and discuss their complementarity. Additionally, cutting-edge XFELs and MeV-UED facilities for TR x-ray and electron scattering experiments and the experiments performed at those facilities are reviewed. As new facilities are constructed and existing ones undergo upgrades, the landscape for TR x-ray and electron scattering experiments is poised for further expansion. Through this review, we aim to facilitate the effective utilization of these emerging opportunities, assisting researchers in delving deeper into the intricate dynamics of molecular structures.

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.

Virial coefficients of the uniform electron gas from path-integral Monte Carlo simulations

The properties of plasmas in the low-density limit are described by virial expansions. Analytical expressions are known from Green's function approaches only for the first three virial coefficients. Accurate path-integral Monte Carlo (PIMC) simulations have recently been performed for the uniform electron gas, allowing the virial expansions to be analyzed and interpolation formulas to be derived. The exact expression for the second virial coefficient is used to test the accuracy of the PIMC simulations and the range of validity of the interpolation formula of Groth et al. [Phys. Rev. Lett. 119, 135001 (2017)0031-900710.1103/PhysRevLett.119.135001], and we discuss the fourth virial coefficient, which is not exactly known yet. Combining PIMC simulations with benchmarks from exact virial expansion results would allow us to obtain more accurate representations of the equation of state for parameter ranges of conditions which are of interest, e.g., for helioseismology.

Interplay of Mechanochemistry and Material Processes in the Graphite to Diamond Phase Transformation

The manifestation of intramolecular strains in covalent systems is widely known to accelerate chemical reactions and open alternative reaction paths. This process is moderately well understood for isolated molecules and unimolecular processes. However, in condensed matter processes such as phase transformations, material properties and structure may influence typical mechanochemical effects. Therefore, we utilize steered molecular dynamics to induce out of plane strains in graphite and compress the system under a constant strain rate to induce phase transformation. We show that the out of plane strain allows phase transformations to initiate at small amounts of compressive strain. However, in contrast to typical mechanochemical results, the sum of compressive and out of plane work needed to form a diamond has a local minimum due to altered defect formation processes during phase transformation. Additionally, these altered processes slow the kinetics of the phase transformation, taking longer from initiation to total material transformation.

Analysing the dynamic structure of warm dense matter in the imaginary-time domain: theoretical models and simulations

Rigorous diagnostics of experiments with warm dense matter are notoriously difficult. A key method is X-ray Thomson scattering (XRTS), but the interpretation of XRTS measurements is usually based on theoretical models that entail various approximations. Recently, Dornheim et al. [Nat. Commun. 13, 7911 (2022)] introduced a new framework for temperature diagnostics of XRTS experiments that is based on imaginary-time correlation functions. On the one hand, switching from the frequency to the imaginary-time domain gives one direct access to a number of physical properties, which facilitates the extraction of the temperature of arbitrarily complex materials without relying on any models or approximations. On the other hand, the bulk of theoretical work in dynamic quantum many-body theory is devoted to the frequency domain, and, to the best of our knowledge, the manifestation of physics properties within the imaginary-time density-density correlation function (ITCF) remains poorly understood. In the present work, we aim to fill this gap by introducing a simple, semi-analytical model for the imaginary-time dependence of two-body correlations within the framework of imaginary-time path integrals. As a practical example, we compare our new model to extensive ab initio path integral Monte Carlo results for the ITCF of a uniform electron gas, and find excellent agreement over a broad range of wavenumbers, densities and temperatures. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Cu3As: Uncommon Crystallographic Features, Low-Temperature Phase Transitions, Thermodynamic and Physical Properties

The formation and crystal structure of the binary Cu₃As phase have been re-investigated. Some physical properties were then measured on both single crystal and polycrystalline bulk. Cu₃As melts congruently at 835 °C. At room temperature (RT), this compound has been found to crystallize in the hexagonal Cu₃P prototype (hP24, P63cm) with lattice parameters: a = 7.1393(1) Å and c = 7.3113(1) Å, rather than in the anti HoH₃-type (hP24, P-3c1) as indicated in literature. A small compositional range of 74.0-75.5 at.% Cu (26.0-24.5 at.% As) was found for samples synthesized at 300 and 400 °C; a corresponding slight understoichiometry is found in one out of the four Cu atomic sites, leading to the final refined composition Cu_2.882(1)As. The present results disprove a change in the crystal structure above RT actually reported in the phase diagram (from γ' to γ on heating). Instead, below RT, at T = 243 K (-30 °C), a first-order structural transition to a trigonal low-temperature superstructure, LT-Cu_3-xAs (hP72, P-3c1) has been found. The LT polymorph is metrically related to the RT one, having the c lattice parameter three times larger: a = 7.110(2) Å and c = 21.879(4) Å. Both the high- and low-temperature polymorphs are characterized by the presence of a tridimensional (3D) uncommon and rigid Cu sublattice of the lonsdaleite type (Cu atoms tetrahedrally bonded), which remains almost unaffected by the structural change(s), and characteristic layers of triangular 'Cu₃As'-units (each hosting one As atom at the center, interconnected each other by sharing the three vertices). The first-order transition is then followed by an additional structural change when lowering the temperature, which induces doubling of also the lattice parameter a. Differential scanning calorimetry nicely detects the first low-temperature structural change occurring at T = 243 K, with an associated enthalpy difference, ΔH_(TR), of approximately 2 J/g (0.53 kJ/mol). Low-temperature electrical resistivity shows a typical metallic behavior; clear anomalies are detected in correspondence to the solid-state transformations. The Seebeck coefficient, measured as a function of temperature, highlights a conduction of n-type. The temperature dependence of the magnetic susceptibility displays an overall constant diamagnetic response.

Ab Initio Static Exchange-Correlation Kernel across Jacob's Ladder without Functional Derivatives

The electronic exchange─correlation (XC) kernel constitutes a fundamental input for the estimation of a gamut of properties such as the dielectric characteristics, the thermal and electrical conductivity, or the response to an external perturbation. In this work, we present a formally exact methodology for the computation of the system specific static XC kernel exclusively within the framework of density functional theory (DFT) and without employing functional derivatives─no external input apart from the usual XC-functional is required. We compare our new results with exact quantum Monte Carlo (QMC) data for the archetypical uniform electron gas model under both ambient and warm dense matter conditions. This gives us unprecedented insights into the performance of different XC functionals, and it has important implications for the development of new functionals that are designed for the application at extreme temperatures. In addition, we obtain new DFT results for the XC kernel of warm dense hydrogen as it occurs in fusion applications and astrophysical objects. The observed excellent agreement to the QMC reference data demonstrates that presented framework is capable to capture nontrivial effects such as XC-induced isotropy breaking in the density response of hydrogen at large wave numbers.

Production of diamond using intense heavy ion beams at the FAIR facility and application to planetary physics

Diamonds are supposedly abundantly present in different objects in the Universe including meteorites, carbon-rich stars as well as carbon-rich extrasolar planets. Moreover, the prediction that in deep layers of Uranus and Neptune, methane may undergo a process of phase separation into diamond and hydrogen, has been experimentally verified. In particular, high power lasers have been used to study this problem. It is therefore important from the point of view of astrophysics and planetary physics, to further study the production processes of diamond in the laboratory. In the present paper, we present numerical simulations of implosion of a solid carbon sample using an intense uranium beam that is to be delivered by the heavy ion synchrotron, SIS100, that is under construction at the Facility for Antiprotons and Ion Research (FAIR), at Darmstadt. These calculations show that using our proposed experimental scheme, one can generate the extreme pressure and temperature conditions, necessary to produce diamonds of mm3 dimensions.

Accurate temperature diagnostics for matter under extreme conditions

The experimental investigation of matter under extreme densities and temperatures, as in astrophysical objects and nuclear fusion applications, constitutes one of the most active frontiers at the interface of material science, plasma physics, and engineering. The central obstacle is given by the rigorous interpretation of the experimental results, as even the diagnosis of basic parameters like the temperature T is rendered difficult at these extreme conditions. Here, we present a simple, approximation-free method to extract the temperature of arbitrarily complex materials in thermal equilibrium from X-ray Thomson scattering experiments, without the need for any simulations or an explicit deconvolution. Our paradigm can be readily implemented at modern facilities and corresponding experiments will have a profound impact on our understanding of warm dense matter and beyond, and open up a variety of appealing possibilities in the context of thermonuclear fusion, laboratory astrophysics, and related disciplines.

Sequential Lonsdaleite to Diamond Formation in Ureilite Meteorites via In Situ Chemical Fluid/Vapor Deposition

Ureilite meteorites are arguably our only large suite of samples from the mantle of a dwarf planet and typically contain greater abundances of diamond than any known rock. Some also contain lonsdaleite, which may be harder than diamond. Here, we use electron microscopy to map the relative distribution of coexisting lonsdaleite, diamond, and graphite in ureilites. These maps show that lonsdaleite tends to occur as polycrystalline grains, sometimes with distinctive fold morphologies, partially replaced by diamond + graphite in rims and cross-cutting veins. These observations provide strong evidence for how the carbon phases formed in ureilites, which, despite much conjecture and seemingly conflicting observations, has not been resolved. We suggest that lonsdaleite formed by pseudomorphic replacement of primary graphite shapes, facilitated by a supercritical C-H-O-S fluid during rapid decompression and cooling. Diamond + graphite formed after lonsdaleite via ongoing reaction with C-H-O-S gas. This graphite > lonsdaleite > diamond + graphite formation process is akin to industrial chemical vapor deposition but operates at higher pressure (∼1-100 bar) and provides a pathway toward manufacture of shaped lonsdaleite for industrial application. It also provides a unique model for ureilites that can reconcile all conflicting observations relating to diamond formation.

Crystal structure, electronic structure, and optical properties of the novel Li4CdGe2S7, a wide-bandgap quaternary sulfide with a polar structure derived from lonsdaleite

The novel quaternary thiogermanate Li4CdGe2S7 (tetralithium cadmium digermanium heptasulfide) was discovered from a solid-state reaction at 750 °C. Single-crystal X-ray diffraction data were collected and used to solve and refine the structure. Li4CdGe2S7 is a member of the small, but growing, class of I4–II–IV2–VI7 diamond-like materials. The compound adopts the Cu5Si2S7 structure type, which is a derivative of lonsdaleite. Crystallizing in the polar space group Cc, Li4CdGe2S7 contains 14 crystallographically unique ions, all residing on general positions. Like all diamond-like structures, the compound is built of corner-sharing tetrahedral units that create a relatively dense three-dimensional assembly. The title compound is the major phase of the reaction product, as evidenced by powder X-ray diffraction and optical diffuse reflectance spectroscopy. While the compound exhibits a second-harmonic generation (SHG) response comparable to that of the AgGaS2 (AGS) reference material in the IR region, its laser-induced damage threshold (LIDT) is over an order of magnitude greater than AGS for λ = 1.064 µm and τ = 30 ps. Bond valence sums, global instability index, minimum bounding ellipsoid (MBE) analysis, and electronic structure calculations using density functional theory (DFT) were used to further evaluate the crystal structure and electronic structure of the compound and provide a comparison with the analogous I2–II–IV–VI4 diamond-like compound Li2CdGeS4. Li4CdGe2S7 appears to be a better IR nonlinear optical (NLO) candidate than Li2CdGeS4 and one of the most promising contenders to date. The exceptional LIDT is likely due, at least in part, to the wider optical bandgap of ∼3.6 eV.

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

Regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

Nonequilibrium band occupation and optical response of gold after ultrafast XUV excitation

Free electron lasers offer unique properties to study matter in states far from equilibrium as they combine short pulses with a large range of photon energies. In particular, the possibility to excite core states drives new relaxation pathways that, in turn, also change the properties of the optically and chemically active electrons. Here, we present a theoretical model for the dynamics of the nonequilibrium occupation of the different energy bands in solid gold driven by exciting deep core states. The resulting optical response is in excellent agreement with recent measurements and, combined with our model, provides a quantitative benchmark for the description of electron-phonon coupling in strongly driven gold. Focusing on sub-picosecond time scales, we find essential differences between the dynamics induced by XUV and visible light.

A Review of Binderless Polycrystalline Diamonds: Focus on the High-Pressure-High-Temperature Sintering Process

Nowadays, synthetic diamonds are easy to fabricate industrially, and a wide range of methods were developed during the last century. Among them, the high-pressure-high-temperature (HP-HT) process is the most used to prepare diamond compacts for cutting or drilling applications. However, these diamond compacts contain binder, limiting their mechanical and optical properties and their substantial uses. Binderless diamond compacts were synthesized more recently, and important developments were made to optimize the P-T conditions of sintering. Resulting sintered compacts had mechanical and optical properties at least equivalent to that of natural single crystal and higher than that of binder-containing sintered compacts, offering a huge potential market. However, pressure-temperature (P-T) conditions to sinter such bodies remain too high for an industrial transfer, making this the next challenge to be accomplished. This review gives an overview of natural diamond formation and the main experimental techniques that are used to synthesize and/or sinter diamond powders and compact objects. The focus of this review is the HP-HT process, especially for the synthesis and sintering of binderless diamonds. P-T conditions of the formation and exceptional properties of such objects are discussed and compared with classic binder-diamonds objects and with natural single-crystal diamonds. Finally, the question of an industrial transfer is asked and outlooks related to this are proposed.

Occurrence of SiC and Diamond Polytypes, Chromite and Uranophane in Breccia from Nickel Laterites (New Caledonia): Combined Analyses

Different techniques have been combined to identify the structure and the chemical composition of siliceous breccia from a drill core of nickel laterites in New Caledonia (Tiebaghi mine). XRD analyses show quartz as a major phase. Micro-Raman spectroscopy confirmed the presence of reddish microcrystalline quartz as a major phase with inclusion of microparticles of iron oxides and oxyhydroxide. Lithoclasts present in breccia are composed of lizardite, chrysotile, forsterite, hedenbergite and saponite. The veins cutting through the breccia are filled with Ni-bearing talc. Furthermore, for the first time, we discovered the presence of diamond microcrystals accompanied by moissanite polytypes (SiC), chromite (FeCr2O4) and uranophane crystals (Ca(UO2)2(SiO3OH)2.5(H2O)) and lonsdaleite (2H-[C-C]) in the porosities of the breccia. The origin of SiC and diamond polytypes are attributed to ultrahigh-pressure crystallization in the lower mantle. The SiC and diamond polytypes are inherited from serpentinized peridotites having experienced interaction with a boninitic melt. Serpentinization, then weathering of the peridotites into saprolite, did not affect the resistant SiC polytypes, diamond and lonsdaleite. During karstification and brecciation, silica rich aqueous solutions partly digested the saprolite. Again, the SiC polymorph represent stable relicts from this dissolution process being deposited in breccia pores. Uranophane is a neoformed phase having crystallized from the silica rich aqueous solutions. Our study highlights the need of combining chemical and mineralogical analytical technologies to acquire the most comprehensive information on samples, as well as the value of Raman spectroscopy in characterizing structural properties of porous materials.

Momentum distribution of the uniform electron gas at finite temperature: Effects of spin polarization

We carry out extensive direct path integral Monte Carlo (PIMC) simulations of the uniform electron gas (UEG) at finite temperature for different values of the spin-polarization ξ. This allows us to unambiguously quantify the impact of spin effects on the momentum distribution function n(k) and related properties. We find that interesting physical effects like the interaction-induced increase in the occupation of the zero-momentum state n(0) substantially depend on ξ. Our results further advance the current understanding of the UEG as a fundamental model system, and are of practical relevance for the description of transport properties of warm dense matter in an external magnetic field. All PIMC results are freely available online and can be used as a benchmark for the development of methods and applications.

Fabricating Graphene and Nanodiamonds from Lignin by Femtosecond Laser Irradiation

This study demonstrates a new transformation path from lignin to graphene and nanodiamonds (NDs) by femtosecond laser writing in air at ambient temperature and pressure. Graphene nanoribbon rolls were generated at lower laser power. When the laser power was high, NDs could be obtained apart from graphene and onion-like carbon intermediates. These structures were confirmed by scanning electron microscopy, X-ray diffraction, Raman spectroscopy, X-ray photoelectron spectroscopy, and high-resolution transmission electron microscopy. The effects of laser power and laser writing speed on the structure of laser-induced patterns were investigated. The results show that the laser power was more important than the writing speed for the synthesis of carbon nanoparticles, and high laser power contributed to enhanced electrically conductive performance. Therefore, the direct laser irradiation technique leads a simple, low-cost, and sustainable way to synthesize graphene and NDs and is promising for the fabrication of sensors and electric devices.

Cu4MnGe2S7 and Cu2MnGeS4: two polar thiogermanates exhibiting second harmonic generation in the infrared and structures derived from hexagonal diamond

The new, quaternary diamond-like semiconductor (DLS) Cu₄MnGe₂S₇ was prepared at high-temperature from a stoichiometric reaction of the elements under vacuum. Single crystal X-ray diffraction data were used to solve and refine the structure in the polar space group Cc. Cu₄MnGe₂S₇ features [Ge₂S₇]^6- units and adopts the Cu₅Si₂S₇ structure type that can be considered a derivative of the hexagonal diamond structure. The DLS Cu₂MnGeS₄ with the wurtz-stannite structure was similarly prepared at a lower temperature. The achievement of relatively phase-pure samples, confirmed by X-ray powder diffraction data, was nontrival as differential thermal analysis shows an incongruent melting behaviour for both compounds at relatively high temperature. The dark red Cu₂MnGeS₄ and Cu₄MnGe₂S₇ compounds exhibit direct optical bandgaps of 2.21 and 1.98 eV, respectively. The infrared (IR) spectra indicate potentially wide windows of optical transparency up to 25 μm for both materials. Using the Kurtz-Perry powder method, the second-order nonlinear optical susceptibility, χ⁽²⁾, values for Cu₂MnGeS₄ and Cu₄MnGe₂S₇ were estimated to be 16.9 ± 2.0 pm V^-1 and 2.33 ± 0.86 pm V^-1, respectively, by comparing with an optical-quality standard reference material, AgGaSe₂ (AGSe). Cu₂MnGeS₄ was found to be phase matchable at λ = 3100 nm, whereas Cu₄MnGe₂S₇ was determined to be non-phase matchable at λ = 1600 nm. The weak SHG response of Cu₄MnGe₂S₇ precluded phase-matching studies at longer wavelengths. The laser-induced damage threshold (LIDT) for Cu₂MnGeS₄ was estimated to be ∼0.1 GW cm^-2 at λ = 1064 nm (pulse width: τ = 30 ps), while the LIDT for Cu₄MnGe₂S₇ could not be ascertained due to its weak response. The significant variance in NLO properties can be reasoned using the results from electronic structure calculations.

Carbon clusters formed from shocked benzene

Benzene (C₆H₆), while stable under ambient conditions, can become chemically reactive at high pressures and temperatures, such as under shock loading conditions. Here, we report in situ x-ray diffraction and small angle x-ray scattering measurements of liquid benzene shocked to 55 GPa, capturing the morphology and crystalline structure of the shock-driven reaction products at nanosecond timescales. The shock-driven chemical reactions in benzene observed using coherent XFEL x-rays were a complex mixture of products composed of carbon and hydrocarbon allotropes. In contrast to the conventional description of diamond, methane and hydrogen formation, our present results indicate that benzene’s shock-driven reaction products consist of layered sheet-like hydrocarbon structures and nanosized carbon clusters with mixed sp²-sp³ hybridized bonding. Implications of these findings range from guiding shock synthesis of novel compounds to the fundamentals of carbon transport in planetary physics.

Shock-wave driven reactions of organic molecules may have played a key role in prebiotic chemistry, but their mechanisms are difficult to investigate. The authors, using time-resolved x-ray diffraction and small-angle x-ray scattering experiments, observe the transformation of liquid benzene during a shock, identifying carbon and hydrocarbon solid products.

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.

Graphite to diamond transition induced by photoelectric absorption of ultraviolet photons

The phase transition from graphite to diamond is an appealing object of study because of many fundamental and also, practical reasons. The out-of-plane distortions required for the transition are a good tool to understand the collective behaviour of layered materials (graphene, graphite) and the van der Waals forces. As today, two basic processes have been successfully tested to drive this transition: strong shocks and high energy femtolaser excitation. They induce it by increasing either pressure or temperature on graphite. In this work, we report a third method consisting in the irradiation of graphite with ultraviolet photons of energies above 4.4 eV. We show high resolution electron microscopy images of pyrolytic carbon evidencing the dislocation of the superficial graphitic layers after irradiation and the formation of crystallite islands within them. Electron energy loss spectroscopy of the islands show that the sp² to sp³ hybridation transition is a surface effect. High sensitivity X-ray diffraction experiments and Raman spectroscopy confirm the formation of diamond within the islands.

X-ray powder diffraction in reflection geometry on multi-beam kJ-type laser facilities

An ultrafast x-ray powder diffraction setup for laser-driven dynamic compression has been developed at the LULI2000 laser facility. X-ray diffraction is performed in reflection geometry from a quasi-monochromatic laser-generated plasma x-ray source. In comparison to a transmission geometry setup, this configuration allows us to probe only a small portion of the compressed sample, as well as to shield the detectors against the x-rays generated by the laser-plasma interaction on the front side of the target. Thus, this new platform facilitates probing of spatially and temporarily uniform thermodynamic conditions and enables us to study samples of a large range of atomic numbers, thicknesses, and compression dynamics. As a proof-of-concept, we report direct structural measurements of the bcc-hcp transition both in shock and ramp-compressed polycrystalline iron with diffraction signals recorded between 2θ ∼ 30° and ∼150°. In parallel, the pressure and temperature history of probed samples is measured by rear-side visible diagnostics (velocimetry and pyrometry).

Finite-size effects in the reconstruction of dynamic properties from ab initio path integral Monte Carlo simulations

We systematically investigate finite-size effects in the dynamic structure factor S(q,ω) of the uniform electron gas obtained via the analytic continuation of ab initio path integral Monte Carlo data for the imaginary-time density-density correlation function F(q,τ). Using the recent scheme by Dornheim et al. [Phys. Rev. Lett. 121, 255001 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.255001], we find that the reconstructed spectra are not afflicted with any finite-size effects for as few as N=14 electrons both at warm dense matter (WDM) conditions and at the margins of the strongly correlated electron liquid regime. Our results further corroborate the high quality of our current description of the dynamic density response of correlated electrons, which is of high importance for many applications in WDM theory and beyond.

Investigation of Room Temperature Formation of the Ultra-Hard Nanocarbons Diamond and Lonsdaleite

Diamond is an attractive material due to its extreme hardness, high thermal conductivity, quantum optical, and biomedical applications. There is still much that is not understood about how diamonds form, particularly at room temperature and without catalysts. In this work, a new route for the formation of nanocrystalline diamond and the diamond-like phase lonsdaleite is presented. Both diamond phases are found to form together within bands with a core-shell structure following the high pressure treatment of a glassy carbon precursor at room temperature. The crystallographic arrangements of the diamond phases revealed that shear is the driving force for their formation and growth. This study gives new understanding of how shear can lead to crystallization in materials and helps elucidate how diamonds can form on Earth, in meteorite impacts and on other planets. Finally, the new shear induced formation mechanism works at room temperature, a key finding that may enable diamond and other technically important nanomaterials to be synthesized more readily.

An approach for the measurement of the bulk temperature of single crystal diamond using an X-ray free electron laser

We present a method to determine the bulk temperature of a single crystal diamond sample at an X-Ray free electron laser using inelastic X-ray scattering. The experiment was performed at the high energy density instrument at the European XFEL GmbH, Germany. The technique, based on inelastic X-ray scattering and the principle of detailed balance, was demonstrated to give accurate temperature measurements, within [Formula: see text] for both room temperature diamond and heated diamond to 500 K. Here, the temperature was increased in a controlled way using a resistive heater to test theoretical predictions of the scaling of the signal with temperature. The method was tested by validating the energy of the phonon modes with previous measurements made at room temperature using inelastic X-ray scattering and neutron scattering techniques. This technique could be used to determine the bulk temperature in transient systems with a temporal resolution of 50 fs and for which accurate measurements of thermodynamic properties are vital to build accurate equation of state and transport models.

Investigating off-Hugoniot states using multi-layer ring-up targets

Laser compression has long been used as a method to study solids at high pressure. This is commonly achieved by sandwiching a sample between two diamond anvils and using a ramped laser pulse to slowly compress the sample, while keeping it cool enough to stay below the melt curve. We demonstrate a different approach, using a multilayer ‘ring-up’ target whereby laser-ablation pressure compresses Pb up to 150 GPa while keeping it solid, over two times as high in pressure than where it would shock melt on the Hugoniot. We find that the efficiency of this approach compares favourably with the commonly used diamond sandwich technique and could be important for new facilities located at XFELs and synchrotrons which often have higher repetition rate, lower energy lasers which limits the achievable pressures that can be reached.

Femtosecond laser produced periodic plasma in a colloidal crystal probed by XFEL radiation

With the rapid development of short-pulse intense laser sources, studies of matter under extreme irradiation conditions enter further unexplored regimes. In addition, an application of X-ray Free-Electron Lasers (XFELs) delivering intense femtosecond X-ray pulses, allows to investigate sample evolution in IR pump - X-ray probe experiments with an unprecedented time resolution. Here we present a detailed study of the periodic plasma created from the colloidal crystal. Both experimental data and theory modeling show that the periodicity in the sample survives to a large extent the extreme excitation and shock wave propagation inside the colloidal crystal. This feature enables probing the excited crystal, using the powerful Bragg peak analysis, in contrast to the conventional studies of dense plasma created from bulk samples for which probing with Bragg diffraction technique is not possible. X-ray diffraction measurements of excited colloidal crystals may then lead towards a better understanding of matter phase transitions under extreme irradiation conditions.

Demonstration of X-ray Thomson scattering as diagnostics for miscibility in warm dense matter

The gas and ice giants in our solar system can be seen as a natural laboratory for the physics of highly compressed matter at temperatures up to thousands of kelvins. In turn, our understanding of their structure and evolution depends critically on our ability to model such matter. One key aspect is the miscibility of the elements in their interiors. Here, we demonstrate the feasibility of X-ray Thomson scattering to quantify the degree of species separation in a 1:1 carbon-hydrogen mixture at a pressure of ~150 GPa and a temperature of ~5000 K. Our measurements provide absolute values of the structure factor that encodes the microscopic arrangement of the particles. From these data, we find a lower limit of [Formula: see text]% of the carbon atoms forming isolated carbon clusters. In principle, this procedure can be employed for investigating the miscibility behaviour of any binary mixture at the high-pressure environment of planetary interiors, in particular, for non-crystalline samples where it is difficult to obtain conclusive results from X-ray diffraction. Moreover, this method will enable unprecedented measurements of mixing/demixing kinetics in dense plasma environments, e.g., induced by chemistry or hydrodynamic instabilities.

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.

A spherical crystal diffraction imager for Sandia's Z Pulsed Power Facility

Sandia's Z Pulsed Power Facility is able to dynamically compress matter to extreme states with exceptional uniformity, duration, and size, which are ideal for investigating fundamental material properties of high energy density conditions. X-ray diffraction (XRD) is a key atomic scale probe since it provides direct observation of the compression and strain of the crystal lattice and is used to detect, identify, and quantify phase transitions. Because of the destructive nature of Z-Dynamic Material Property (DMP) experiments and low signal vs background emission levels of XRD, it is very challenging to detect a diffraction signal close to the Z-DMP load and to recover the data. We have developed a new Spherical Crystal Diffraction Imager (SCDI) diagnostic to relay and image the diffracted x-ray pattern away from the load debris field. The SCDI diagnostic utilizes the Z-Beamlet laser to generate 6.2-keV Mn-He_α x rays to probe a shock-compressed material on the Z-DMP load. A spherically bent crystal composed of highly oriented pyrolytic graphite is used to collect and focus the diffracted x rays into a 1-in. thick tungsten housing, where an image plate is used to record the data.

Ab initio dielectric response function of diamond and other relevant high pressure phases of carbon

The electronic structure and dielectric properties of the diamond, body centered cubic diamond (bc8), and hexagonal diamond (lonsdaleite) phases of carbon are computed using density functional theory and many-body perturbation theory with the emphasis on the excitonic picture of the solid phases relevant in the regimes of high-pressure physics and warm dense matter. We also discuss the capabilities of reproducing the inelastic x-ray scattering spectra in comparison with the existing models in light of recent x-ray scattering experiments on carbon and carbon bearing materials in the Megabar range.

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.

Nanoleite: a new semiconducting carbon allotrope predicted by density functional theory

A new carbon allotrope with an indirect bandgap of 2.06 eV has been predicted by density functional theory, which has a high absorption coefficient in the visible spectral range that is suitable for solar cell application.

Quantifying hexagonal stacking in diamond

Diamond is a material of immense technological importance and an ancient signifier for wealth and societal status. In geology, diamond forms as part of the deep carbon cycle and typically displays a highly ordered cubic crystal structure. Impact diamonds, however, often exhibit structural disorder in the form of complex combinations of cubic and hexagonal stacking motifs. The structural characterization of such diamonds remains a challenge. Here, impact diamonds from the Popigai crater were characterized with a range of techniques. Using the MCDIFFaX approach for analysing X-ray diffraction data, hexagonality indices up to 40% were found. The effects of increasing amounts of hexagonal stacking on the Raman spectra of diamond were investigated computationally and found to be in excellent agreement with trends in the experimental spectra. Electron microscopy revealed nanoscale twinning within the cubic diamond structure. Our analyses lead us to propose a systematic protocol for assigning specific hexagonality attributes to the mineral designated as lonsdaleite among natural and synthetic samples.

Stabilizing the metastable superhard material wurtzite boron nitride by three-dimensional networks of planar defects

Wurtzite boron nitride (w-BN) is a metastable superhard material that is a high-pressure polymorph of BN. Clarifying how the metastable high-pressure material can be stabilized at atmospheric pressure is a challenging issue of fundamental scientific importance and promising technological value. Here, we fabricate millimeter-size w-BN bulk crystals via the hexagonal-to-wurtzite phase transformation at high pressure and high temperature. By combining transmission electron microscopy and ab initio molecular dynamics simulations, we reveal a stabilization mechanism for w-BN, i.e., the metastable high-pressure phase can be stabilized by 3D networks of planar defects which are constructed by a high density of intersecting (0001) stacking faults and {10[Formula: see text]0} inversion domain boundaries. The 3D networks of planar defects segment the w-BN bulk crystal into numerous nanometer-size prismatic domains with the reverse crystallographic polarities. Our findings unambiguously demonstrate the retarding effect of crystal defects on the phase transformations of metastable materials, which is in contrast to the common knowledge that the crystal defects in materials will facilitate the occurrence of phase transformations.

Evidence for Crystalline Structure in Dynamically-Compressed Polyethylene up to 200 GPa

We investigated the high-pressure behavior of polyethylene (CH₂) by probing dynamically-compressed samples with X-ray diffraction. At pressures up to 200 GPa, comparable to those present inside icy giant planets (Uranus, Neptune), shock-compressed polyethylene retains a polymer crystal structure, from which we infer the presence of significant covalent bonding. The A2/m structure which we observe has previously been seen at significantly lower pressures, and the equation of state measured agrees with our findings. This result appears to contrast with recent data from shock-compressed polystyrene (CH) at higher temperatures, which demonstrated demixing and recrystallization into a diamond lattice, implying the breaking of the original chemical bonds. As such chemical processes have significant implications for the structure and energy transfer within ice giants, our results highlight the need for a deeper understanding of the chemistry of high pressure hydrocarbons, and the importance of better constraining planetary temperature profiles.

An experimental platform using high-power, high-intensity optical lasers with the hard X-ray free-electron laser at SACLA

An experimental platform using X-ray free-electron laser (XFEL) pulses with high-intensity optical laser pulses is open for early users' experiments at the SACLA XFEL facility after completion of the commissioning. The combination of the hard XFEL and the high-intensity laser provides capabilities to open new frontiers of laser-based high-energy-density science. During the commissioning phase, characterization of the XFEL and the laser at the platform has been carried out for the combinative utilization as well as the development of instruments and basic diagnostics for user experiments. An overview of the commissioning and the current capabilities of the experimental platform is presented.

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.

Liquid Structure of Shock-Compressed Hydrocarbons at Megabar Pressures

We present results for the ionic structure in hydrocarbons (polystyrene, polyethylene) that were shock compressed to pressures of up to 190 GPa, inducing rapid melting of the samples. The structure of the resulting liquid is then probed using in situ diffraction by an x-ray free electron laser beam, demonstrating the capability to obtain reliable diffraction data in a single shot, even for low-Z samples without long range order. The data agree well with ab initio simulations, validating the ability of such approaches to model mixed samples in states where complex interparticle bonds remain, and showing that simpler models are not necessarily valid. While the results clearly exclude the possibility of complete carbon-hydrogen demixing at the conditions probed, they also, in contrast to previous predictions, indicate that diffraction is not always a sufficient diagnostic for this phenomenon.

Femtosecond diffraction studies of solid and liquid phase changes in shock-compressed bismuth

Bismuth has long been a prototypical system for investigating phase transformations and melting at high pressure. Despite decades of experimental study, however, the lattice-level response of Bi to rapid (shock) compression and the relationship between structures occurring dynamically and those observed during slow (static) compression, are still not clearly understood. We have determined the structural response of shock-compressed Bi to 68 GPa using femtosecond X-ray diffraction, thereby revealing the phase transition sequence and equation-of-state in unprecedented detail for the first time. We show that shocked-Bi exhibits a marked departure from equilibrium behavior - the incommensurate Bi-III phase is not observed, but rather a new metastable phase, and the Bi-V phase is formed at significantly lower pressures compared to static compression studies. We also directly measure structural changes in a shocked liquid for the first time. These observations reveal new behaviour in the solid and liquid phases of a shocked material and give important insights into the validity of comparing static and dynamic datasets.

Analysis of stacking disorder in ice I using pair distribution functions

Stacking-disordered materials display crystalline order in two dimensions but are disordered along the direction in which layered structural motifs are stacked. Countless examples of stacking disorder exist, ranging from close-packed metals, ice I and diamond to open-framework materials and small-molecule pharmaceuticals. In general, the presence of stacking disorder can have profound consequences for the physical and chemical properties of a material. Traditional analyses of powder diffraction data are often complicated by the presence of memory effects in the stacking sequences. Here it is shown that experimental pair distribution functions of stacking-disordered ice I can be used to determine local information on the fractions of cubic and hexagonal stacking. Ice is a particularly challenging material in this respect, since both the stacking disorder and the orientational disorder of the water molecules need to be described. Memory effects are found to contribute very little to the pair distribution functions, and consequently, the analysis of pair distribution functions is the method of choice for characterizing stacking-disordered samples with complicated and high-order memory effects. In the context of this work, the limitations of current structure-reconstruction approaches are also discussed.

Simulations of X-ray diffraction of shock-compressed single-crystal tantalum with synchrotron undulator sources

Polychromatic synchrotron undulator X-ray sources are useful for ultrafast single-crystal diffraction under shock compression. Here, simulations of X-ray diffraction of shock-compressed single-crystal tantalum with realistic undulator sources are reported, based on large-scale molecular dynamics simulations. Purely elastic deformation, elastic-plastic two-wave structure, and severe plastic deformation under different impact velocities are explored, as well as an edge release case. Transmission-mode diffraction simulations consider crystallographic orientation, loading direction, incident beam direction, X-ray spectrum bandwidth and realistic detector size. Diffraction patterns and reciprocal space nodes are obtained from atomic configurations for different loading (elastic and plastic) and detection conditions, and interpretation of the diffraction patterns is discussed.

Comment on "Incipient plasticity of diamond during nanoindentation" by C. Xu, C. Liu and H. Wang, RSC Advances, 2017, , 36093

A recent molecular dynamics simulation study on nanoindentation of diamond carried out by Xu et al. 1 has reported observation of the presence of a controversial hexagonal lonsdaleite phase of carbon in the indentation area. In this comment, we question the reported observation and attribute this anomaly to shortcomings of the long range bond order potential (LCBOP) employed in the nanoindentation study.