Diamond at 800 GPa

Extended X-ray absorption fine structure of dynamically-compressed copper up to 1 terapascal

Large laser facilities have recently enabled material characterization at the pressures of Earth and Super-Earth cores. However, the temperature of the compressed materials has been largely unknown, or solely relied on models and simulations, due to lack of diagnostics under these challenging conditions. Here, we report on temperature, density, pressure, and local structure of copper determined from extended x-ray absorption fine structure and velocimetry up to 1 Terapascal. These results nearly double the highest pressure at which extended x-ray absorption fine structure has been reported in any material. In this work, the copper temperature is unexpectedly found to be much higher than predicted when adjacent to diamond layer(s), demonstrating the important influence of the sample environment on the thermal state of materials; this effect may introduce additional temperature uncertainties in some previous experiments using diamond and provides new guidance for future experimental design.

Double-Shock Compression Pathways from Diamond to BC8 Carbon

Carbon is one of the most important elements for both industrial applications and fundamental research, including life, physics, chemistry, materials, and even planetary science. Although theoretical predictions on the transition from diamond to the BC8 (Ia3[over ¯]) carbon were made more than thirty years ago, after tremendous experimental efforts, direct evidence for the existence of BC8 carbon is still lacking. In this study, a machine learning potential was developed for high-pressure carbon fitted from first-principles calculations, which exhibited great capabilities in modeling the melting and Hugoniot line. Using the molecular dynamics based on this machine learning potential, we designed a thermodynamic pathway that is achievable for the double shock compression experiment to obtain the elusive BC8 carbon. Diamond was compressed up to 584 GPa after the first shock at 20.5  km/s. Subsequently, in the second shock compression at 24.8 or 25.0  km/s, diamond was compressed to a supercooled liquid and then solidified to BC8 in around 1 ns. Furthermore, the critical nucleus size and nucleation rate of BC8 were calculated, which are crucial for nano-second x-ray diffraction measurements to observe BC8 carbon during shock compressions. The key to obtaining BC8 carbon lies in the formation of liquid at a sufficient supercooling. Our work provides a feasible pathway by which the long-sought BC8 phase of carbon can be reached in experiments.

Imaging of the electronic bonding of diamond at pressures up to 2 million atmospheres

Diamond shows unprecedented hardness. Because hardness is a measure of resistance of chemical bonds in a material to external indentation, the electronic bonding nature of diamond beyond several million atmospheres is key to understanding the origin of hardness. However, probing the electronic structures of diamond at such extreme pressure has not been experimentally possible. The measurements on the inelastic x-ray scattering spectra for diamond up to 2 million atmospheres provide data on the evolution of its electronic structures under compression. The mapping of the observed electronic density of states allows us to obtain a two-dimensional image of the bonding transitions of diamond undergoing deformation. The spectral change near edge onset is minor beyond a million atmospheres, while its electronic structure displays marked pressure-induced electron delocalization. Such electronic responses indicate that diamond's external rigidity is supported by its ability to reconcile internal stress, providing insights into the origins of hardness in materials.

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/cm^{3} to 16 g/cm^{3} and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. We find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

Imaging velocity interferometer system for any reflector (VISAR) diagnostics for high energy density sciences

Two variants of optical imaging velocimetry, specifically the one-dimensional streaked line-imaging and the two-dimensional time-resolved area-imaging versions of the Velocity Interferometer System for Any Reflector (VISAR), have become important diagnostics in high energy density sciences, including inertial confinement fusion and dynamic compression of condensed matter. Here, we give a brief review of the historical development of these techniques, then describe the current implementations at major high energy density (HED) facilities worldwide, including the OMEGA Laser Facility and the National Ignition Facility. We illustrate the versatility and power of these techniques by reviewing diverse applications of imaging VISARs for gas-gun and laser-driven dynamic compression experiments for materials science, shock physics, condensed matter physics, chemical physics, plasma physics, planetary science and astronomy, as well as a broad range of HED experiments and laser-driven inertial confinement fusion research.

Exceptionally Enhanced Thermal Conductivity of Aluminum Driven by Extreme Pressures: A First-Principles Study

Extreme pressure conditions reveal fundamental insights into the physical properties of elemental metals that are otherwise not evident under ambient conditions. Herein, we use the density functional perturbation theory to demonstrate that the change in thermal conductivity as a result of large hydrostatic pressures at room temperature for aluminum is the largest of any known material. More specifically, in comparison to ambient conditions, we find that the change in thermal conductivity for aluminum is greater than the relative changes in thermal conductivities of diamond and cubic boron nitride combined, which are two of the most thermally conductive bulk materials known to date. We attribute this to the relatively larger increase in mean free paths and lifetimes of electrons in aluminum as a result of weaker electron-phonon coupling at higher pressures. Our work reveals direct insights into the exceptional electronic transport properties of pressurized aluminum and advances a broad paradigm for understanding thermal transport in metals under extreme pressure.

Structure and density of silicon carbide to 1.5 TPa and implications for extrasolar planets

There has been considerable recent interest in the high-pressure behavior of silicon carbide, a potential major constituent of carbon-rich exoplanets. In this work, the atomic-level structure of SiC was determined through in situ X-ray diffraction under laser-driven ramp compression up to 1.5 TPa; stresses more than seven times greater than previous static and shock data. Here we show that the B1-type structure persists over this stress range and we have constrained its equation of state (EOS). Using this data we have determined the first experimentally based mass-radius curves for a hypothetical pure SiC planet. Interior structure models are constructed for planets consisting of a SiC-rich mantle and iron-rich core. Carbide planets are found to be ~10% less dense than corresponding terrestrial planets.

Metastability of diamond ramp-compressed to 2 terapascals

Carbon is the fourth-most prevalent element in the Universe and essential for all known life. In the elemental form it is found in multiple allotropes, including graphite, diamond and fullerenes, and it has long been predicted that even more structures can exist at pressures greater than those at Earth's core^1-3. Several phases have been predicted to exist in the multi-terapascal regime, which is important for accurate modelling of the interiors of carbon-rich exoplanets^4,5. By compressing solid carbon to 2 terapascals (20 million atmospheres; more than five times the pressure at Earth's core) using ramp-shaped laser pulses and simultaneously measuring nanosecond-duration time-resolved X-ray diffraction, we found that solid carbon retains the diamond structure far beyond its regime of predicted stability. The results confirm predictions that the strength of the tetrahedral molecular orbital bonds in diamond persists under enormous pressure, resulting in large energy barriers that hinder conversion to more-stable high-pressure allotropes^1,2, just as graphite formation from metastable diamond is kinetically hindered at atmospheric pressure. This work nearly doubles the highest pressure at which X-ray diffraction has been recorded on any material.

Thinking about the Development of High-Pressure Experimental Chemistry

High-pressure chemistry is an interdisciplinary science which uses high-pressure experiments and theories to study the interactions, reactions, and transformations among atoms or molecules. It has been extensively studied thus far and achieved rapid development over the past decades. However, what is next for high-pressure chemistry? In this Perspective, we mainly focus on the development of high-pressure experimental chemistry from our own viewpoint. An overview of the series of topics is as follows: (I) high pressure used as an effective tool to help resolve scientific disputes regarding phenomena observed under ambient conditions; (II) high-pressure reactions of interest to synthetic chemists; (III) utilizing chemical methods to quench the high-pressure phase; (IV) using high pressure to achieve what chemists want to do but could not do; (V) potential applications of in situ properties under high pressure. This Perspective is expected to offer future research opportunities for researchers to develop high-pressure chemistry and to inspire new endeavors in this area to promote the field of compression chemistry science.

X-ray diffraction at the National Ignition Facility

We report details of an experimental platform implemented at the National Ignition Facility to obtain in situ powder diffraction data from solids dynamically compressed to extreme pressures. Thin samples are sandwiched between tamper layers and ramp compressed using a gradual increase in the drive-laser irradiance. Pressure history in the sample is determined using high-precision velocimetry measurements. Up to two independently timed pulses of x rays are produced at or near the time of peak pressure by laser illumination of thin metal foils. The quasi-monochromatic x-ray pulses have a mean wavelength selectable between 0.6 Å and 1.9 Å depending on the foil material. The diffracted signal is recorded on image plates with a typical 2θ x-ray scattering angle uncertainty of about 0.2° and resolution of about 1°. Analytic expressions are reported for systematic corrections to 2θ due to finite pinhole size and sample offset. A new variant of a nonlinear background subtraction algorithm is described, which has been used to observe diffraction lines at signal-to-background ratios as low as a few percent. Variations in system response over the detector area are compensated in order to obtain accurate line intensities; this system response calculation includes a new analytic approximation for image-plate sensitivity as a function of photon energy and incident angle. This experimental platform has been used up to 2 TPa (20 Mbar) to determine the crystal structure, measure the density, and evaluate the strain-induced texturing of a variety of compressed samples spanning periods 2-7 on the periodic table.

Theoretical development and experimental validation on the measurement of temperature by extended X-ray absorption fine structure

A systematic investigation on the theoretical framework of the ultra-fast measurement of temperature by extended X-ray absorption fine structure (EXAFS) applied in laser-driven-compression experiments has been carried out and a new temperature measurement scheme based on the EXAFS cumulant expansion analysis and anharmonic correlated Debye model has been advanced. By considering the anharmonic effect of thermal vibration and avoiding the employment of the empirical model as well as parameters which have large inherent uncertainties in the temperature determination, this new scheme is theoretically more accurate than traditional ones. Then the performance of the new measurement scheme and traditional methods were validated on a synchrotron radiation platform by temperature-dependent EXAFS (TDEXAFS) experiments on Au, Fe, V and Ti; the results showed that the new scheme could provide the most accurate measured temperatures with much lower uncertainties. This accurate scheme gives a firmer physical ground to the EXAFS temperature measurement technique and can expect to be applied in laser-driven compression experiments and promote the development of matter state research at extreme conditions.

Heat-release equation of state and thermal conductivity of warm dense carbon by proton differential heating

Warm dense carbon is generated at 0.3-2.0 g/cc and 1-7 eV by proton heating. The release equation of state (EOS) after heating and thermal conductivity of warm dense carbon are studied experimentally in this regime using a Au/C dual-layer target to initiate a temperature gradient and two picosecond time-resolved diagnostics to probe the surface expansion and heat flow. Comparison between the data and simulations using various EOSs and thermal conductivity models is quantified with a statistical χ^{2} analysis. Out of seven EOS tables and five thermal conductivity models, only L9061 with the Lee-More model provides a probability above 50% to match all data.

Non-iterative characteristics analysis for high-pressure ramp loading

In the canonical ramp compression experiment, a smoothly increasing load is applied to the surface of the sample, and the particle velocity history is measured at two or more different distances into the sample, at interfaces where the surface of the sample can be probed. The velocity histories are used to deduce a stress-density relation, usually using iterative Lagrangian analysis to account for the perturbing effect of the impedance mismatch at the interface. In that technique, a stress-density relation is assumed in order to correct for the perturbation and is adjusted until it becomes consistent with the deduced stress-density relation. This process is subject to the usual difficulties of nonlinear optimization, such as the existence of local minima (sensitivity to the initial guess), possible failure to converge, and relatively large computational effort. We show that, by considering the interaction of successive characteristics reaching a free surface, the stress-density relation can be deduced directly by recursion rather than iteration. This calculation is orders of magnitude faster than iterative analysis and does not require an initial guess. Direct recursion may be less suitable for very noisy data, but it was robust when applied to trial data. The stress-density relation deduced was identical to the result from iterative Lagrangian analysis.

Deducing density and strength of nanocrystalline Ta and diamond under extreme conditions from X-ray diffraction

In situ X-ray diffraction with advanced X-ray sources offers unique opportunities for investigating materials properties under extreme conditions such as shock-wave loading. Here, Singh's theory for deducing high-pressure density and strength from two-dimensional (2D) diffraction patterns is rigorously examined with large-scale molecular dynamics simulations of isothermal compression and shock-wave compression. Two representative solids are explored: nanocrystalline Ta and diamond. Analysis of simulated 2D X-ray diffraction patterns is compared against direct molecular dynamics simulation results. Singh's method is highly accurate for density measurement (within 1%) and reasonable for strength measurement (within 10%), and can be used for such measurements on nanocrystalline and polycrystalline solids under extreme conditions (e.g. in the megabar regime).

Design optimization through thermomechanical finite-element analysis of a hybrid piston-clamped anvil cell for nuclear magnetic resonance experiments

The investigation of materials under extreme pressure conditions requires high-performance cells whose design invariably involves trade-offs between the maximum achievable pressure, the allowed sample volume, and the possibility of real-time pressure monitoring. With a newly conceived hybrid piston-clamped anvil cell, we offer a relatively simple and versatile system, suitable for nuclear magnetic resonance experiments up to 4.4 GPa. Finite-element models, taking into account mechanical and thermal conditions, were used to optimize and validate the design prior to the realization of the device. Cell body and gaskets were made of beryllium-copper alloy and the pistons and pusher were made of tungsten carbide, while the anvils consist of zirconium dioxide. The low-temperature pressure cell performance was tested by monitoring in situ the pressure-dependent ⁶³Cu nuclear-quadrupole-resonance signal of Cu₂O.

Crystal structure and equation of state of Fe-Si alloys at super-Earth core conditions

The high-pressure behavior of Fe alloys governs the interior structure and dynamics of super-Earths, rocky extrasolar planets that could be as much as 10 times more massive than Earth. In experiments reaching up to 1300 GPa, we combine laser-driven dynamic ramp compression with in situ x-ray diffraction to study the effect of composition on the crystal structure and density of Fe-Si alloys, a potential constituent of super-Earth cores. We find that Fe-Si alloy with 7 weight % (wt %) Si adopts the hexagonal close-packed structure over the measured pressure range, whereas Fe-15wt%Si is observed in a body-centered cubic structure. This study represents the first experimental determination of the density and crystal structure of Fe-Si alloys at pressures corresponding to the center of a ~3-Earth mass terrestrial planet. Our results allow for direct determination of the effects of light elements on core radius, density, and pressures for these planets.

Free energy model for solid high-pressure phases of carbon

Analytic free energy models for three solid high-pressure phases--diamond, body centered cubic phase with eight atoms in the unit cell (BC8), and simple cubic (SC)--are developed using density functional theory. We explicitly include anharmonic effects by performing molecular dynamics simulations and investigate their density and temperature dependence in detail. Anharmonicity in the nuclear motion shifts the phase transitions significantly compared to the harmonic approximation. Furthermore, we apply a thermodynamically constrained correction that brings the equation of state in accordance with diamond anvil cell experiments. The performance of our thermodynamic functions is validated against Hugoniot experiments.

X-Ray Diffraction of Solid Tin to 1.2 TPa

We report direct in situ measurements of the crystal structure of tin between 0.12 and 1.2 TPa, the highest stress at which a crystal structure has ever been observed. Using angle-dispersive powder x-ray diffraction, we find that dynamically compressed Sn transforms to the body-centered-cubic (bcc) structure previously identified by ambient-temperature quasistatic-compression studies and by zero-kelvin density-functional theory predictions between 0.06 and 0.16 TPa. However, we observe no evidence for the hexagonal close-packed (hcp) phase found by those studies to be stable above 0.16 TPa. Instead, our results are consistent with bcc up to 1.2 TPa. We conjecture that at high temperature bcc is stabilized relative to hcp due to differences in vibrational free energy.

Perspective: crystal structure prediction at high pressures

Crystal structure prediction at high pressures unbiased by any prior known structure information has recently become a topic of considerable interest. We here present a short overview of recently developed structure prediction methods and propose current challenges for crystal structure prediction. We focus on first-principles crystal structure prediction at high pressures, paying particular attention to novel high pressure structures uncovered by efficient structure prediction methods. Finally, a brief perspective on the outstanding issues that remain to be solved and some directions for future structure prediction researches at high pressure are presented and discussed.

Single Hit Energy-resolved Laue Diffraction

In situ white light Laue diffraction has been successfully used to interrogate the structure of single crystal materials undergoing rapid (nanosecond) dynamic compression up to megabar pressures. However, information on strain state accessible via this technique is limited, reducing its applicability for a range of applications. We present an extension to the existing Laue diffraction platform in which we record the photon energy of a subset of diffraction peaks. This allows for a measurement of the longitudinal and transverse strains in situ during compression. Consequently, we demonstrate measurement of volumetric compression of the unit cell, in addition to the limited aspect ratio information accessible in conventional white light Laue. We present preliminary results for silicon, where only an elastic strain is observed. VISAR measurements show the presence of a two wave structure and measurements show that material downstream of the second wave does not contribute to the observed diffraction peaks, supporting the idea that this material may be highly disordered, or has undergone large scale rotation.

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.

Ramp compression of diamond to five terapascals

The recent discovery of more than a thousand planets outside our Solar System, together with the significant push to achieve inertially confined fusion in the laboratory, has prompted a renewed interest in how dense matter behaves at millions to billions of atmospheres of pressure. The theoretical description of such electron-degenerate matter has matured since the early quantum statistical model of Thomas and Fermi, and now suggests that new complexities can emerge at pressures where core electrons (not only valence electrons) influence the structure and bonding of matter. Recent developments in shock-free dynamic (ramp) compression now allow laboratory access to this dense matter regime. Here we describe ramp-compression measurements for diamond, achieving 3.7-fold compression at a peak pressure of 5 terapascals (equivalent to 50 million atmospheres). These equation-of-state data can now be compared to first-principles density functional calculations and theories long used to describe matter present in the interiors of giant planets, in stars, and in inertial-confinement fusion experiments. Our data also provide new constraints on mass-radius relationships for carbon-rich planets.

Ramp compression of diamond to five terapascals

The recent discovery of more than a thousand planets outside our Solar System, together with the significant push to achieve inertially confined fusion in the laboratory, has prompted a renewed interest in how dense matter behaves at millions to billions of atmospheres of pressure. The theoretical description of such electron-degenerate matter has matured since the early quantum statistical model of Thomas and Fermi, and now suggests that new complexities can emerge at pressures where core electrons (not only valence electrons) influence the structure and bonding of matter. Recent developments in shock-free dynamic (ramp) compression now allow laboratory access to this dense matter regime. Here we describe ramp-compression measurements for diamond, achieving 3.7-fold compression at a peak pressure of 5 terapascals (equivalent to 50 million atmospheres). These equation-of-state data can now be compared to first-principles density functional calculations and theories long used to describe matter present in the interiors of giant planets, in stars, and in inertial-confinement fusion experiments. Our data also provide new constraints on mass-radius relationships for carbon-rich planets.

Single photon energy dispersive x-ray diffraction

With the pressure range accessible to laser driven compression experiments on solid material rising rapidly, new challenges in the diagnosis of samples in harsh laser environments are emerging. When driving to TPa pressures (conditions highly relevant to planetary interiors), traditional x-ray diffraction techniques are plagued by increased sources of background and noise, as well as a potential reduction in signal. In this paper we present a new diffraction diagnostic designed to record x-ray diffraction in low signal-to-noise environments. By utilising single photon counting techniques we demonstrate the ability to record diffraction patterns on nanosecond timescales, and subsequently separate, photon-by-photon, signal from background. In doing this, we mitigate many of the issues surrounding the use of high intensity lasers to drive samples to extremes of pressure, allowing for structural information to be obtained in a regime which is currently largely unexplored.

Stable all-nitrogen metallic salt at terapascal pressures

The phase diagram and equation of state of dense nitrogen are of interest in understanding the fundamental physics and chemistry under extreme conditions, including planetary processes, and in discovering new materials. We predict several stable phases of nitrogen at multi-TPa pressures, including a P4/nbm structure consisting of partially charged N(2)(δ+) pairs and N(5)(δ-) tetrahedra, which is stable in the range 2.5-6.8 TPa. This is followed by a modulated layered structure between 6.8 and 12.6 TPa, which also exhibits a significant charge transfer. The P4/nbm metallic nitrogen salt and the modulated structure are stable at high pressures and temperatures, and they exhibit strongly ionic features and charge density distortions, which is unexpected in an element under such extreme conditions and could represent a new class of nitrogen materials. The P-T phase diagram of nitrogen at TPa pressures is investigated using quasiharmonic phonon calculations and ab initio molecular dynamics simulations.

Solid iron compressed up to 560 GPa

Dynamic compression by multiple shocks is used to compress iron up to 560 GPa (5.6 Mbar), the highest solid-state pressure yet attained for iron in the laboratory. Extended x-ray absorption fine structure (EXAFS) spectroscopy offers simultaneous density, temperature, and local-structure measurements for the compressed iron. The data show that the close-packed structure of iron is stable up to 560 GPa, the temperature at peak compression is significantly higher than expected from pure compressive work, and the dynamic strength of iron is many times greater than the static strength based on lower pressure data. The results provide the first constraint on the melting line of iron above 400 GPa.

Physical and chemical transformations of highly compressed carbon dioxide at bond energies

Carbon dioxide exhibits a richness of high-pressure polymorphs with a great diversity in intermolecular interaction, chemical bonding, and crystal structures. It ranges from typical molecular solids to fully extended covalent solids with crystal structures similar to those of SiO2. These extended solids of carbon dioxide are fundamentally new materials exhibiting interesting optical nonlinearity, low compressibility and high energy density. Furthermore, the large disparity in chemical bonding between the extended network and molecular structures results in a broad metastability domain for these phases to room temperature and almost to ambient pressure and thereby offers enhanced opportunities for novel materials developments. Broadly speaking, these molecular-to-non-molecular transitions occur due to electron delocalization manifested as a rapid increase in electron kinetic energy at high density. The detailed mechanisms, however, are more complex with phase metastabilities, path-dependent phases and phase boundaries, and large lattice strains and structural distortions - all of which are controlled by well beyond thermodynamic constraints to chemical kinetics associated with the governing phases and transitions. As a result, the equilibrium phase boundary is difficult to locate precisely (experimentally or theoretically) and is often obscured by the presence of metastable phases (ordered or disordered). This paper will review the pressure-induced transformations observed in highly compressed carbon dioxide and present chemistry perspectives on those molecular-to-non-molecular transformations that can be applied to other low-Z molecular solids at Mbar pressures where the compression energy rivals the chemical bond energies.

The liquid<−>amorphous transition and the high pressure phase diagram of carbon

The phase diagram of carbon is mapped to high pressure using a computationally-tractable potential model. The use of a relatively simple (Tersoff-II) potential model allows a large range of phase space to be explored. The coexistence (melting) curve for the diamond crystal/liquid dyad is mapped directly by modelling the solid/liquid interfaces. The melting curve is found to be re-entrant and belongs to a conformal class of diamond/liquid coexistence curves. On supercooling the liquid a phase transition to a tetrahedral amorphous form (ta-C) is observed. The liquid <−> amorphous coexistence curve is mapped onto the pT plane and is found to also be re-entrant. The entropy changes for both melting and the amorphous −> liquid transitions are obtained from the respective coexistence curves and the associated changes in molar volume. The structural change on amorphization is analysed at different points on the coexistence curve including for transitions that are both isochoric and isocoordinate (no change in nearest-neighbour coordination number). The conformal nature of the melting curve is highlighted with respect to the known behaviour of Si. The relationship of the observed liquid/amorphous coexistence curve to the Si low- and high-density amorphous (LDA/HDA) transition is discussed.

Thermodynamically stable phases of carbon at multiterapascal pressures

Phases of carbon are studied up to pressures of 1 petapascal (PPa) using first-principles density-functional-theory methods and a structure searching algorithm. Our extensive search over the potential energy surface supports the sequence of transitions diamond → BC8 → simple cubic under increasing pressure found in previous theoretical studies. At higher pressures we predict a soft-phonon driven transition to a simple hexagonal structure at 6.4 terapascals (TPa), and further transitions to the face centered cubic electride structure at 21 TPa, a double hexagonal close packed structure at 270 TPa, and the body centered cubic structure at 650 TPa.

Persistence and eventual demise of oxygen molecules at terapascal pressures

Computational searches for structures of solid oxygen under high pressures in the multi-TPa range are carried out using density-functional-theory methods. We find that molecular oxygen persists to about 1.9 TPa at which it transforms into a semiconducting square-spiral-like polymeric structure (I4(1)/acd) with a band gap of ~3.0 eV. Solid oxygen forms a metallic zigzag chainlike structure (Cmcm) at about 3.0 TPa, but the chains in each layer gradually merge as the pressure is increased and a structure of Fmmm symmetry forms at about 9.3 TPa in which each atom has four nearest neighbors. The superconducting properties of molecular oxygen do not vary much with compression, although the structure becomes more symmetric. The electronic properties of oxygen have a complex evolution with pressure, swapping between insulating, semiconducting, and metallic.

High-pressure phases of hydrogen cyanide: formation of hydrogenated carbon nitride polymers and layers and their electronic properties

We have performed a set of first-principles simulations to consider the possible phase transitions in molecular crystals of HCN under high pressure. Our calculations reveal several transition paths from the orthorhombic phase to tetragonal and then to triclinic phases. The transitions from the orthorhombic to the tetragonal phases are of the second order, whereas those from the tetragonal to the triclinic phases turn out to be of the first-order type and characterized by an abrupt decrease in volume. Our calculations show that, by adjustment of the temperature and pressure of the HCN molecular crystal, novel layered and polymeric crystals with insulating, semiconducting or metallic properties can be found. Based on our simulation results, two different crystal formation mechanisms are deduced. The stabilities of the predicted structures at ambient pressure are further assessed by performing phonon or MD simulations. In addition, the electron transport properties of the predicted polymers are obtained using the non-equilibrium Green's function technique combined with density functional theory. The results show that the polymers have metallic-like I-V characteristics.

Evidence of superdense aluminium synthesized by ultrafast microexplosion

At extreme pressures and temperatures, such as those inside planets and stars, common materials form new dense phases with compacted atomic arrangements and unusual physical properties. The synthesis and study of new phases of matter at pressures above 100 GPa and temperatures above 10(4) K--warm dense matter--may reveal the functional details of planet and star interiors, and may lead to materials with extraordinary properties. Many phases have been predicted theoretically that may be realized once appropriate formation conditions are found. Here we report the synthesis of a superdense stable phase of body-centred-cubic aluminium, predicted by first-principles theories to exist at pressures above 380 GPa. The superdense Al phase was synthesized in the non-equilibrium conditions of an ultrafast laser-induced microexplosion confined inside sapphire (α-Al(2)O(3)). Confined microexplosions offer a strategy to create and recover high-density polymorphs, and a simple method for tabletop study of warm dense matter.

Experimental determination of third-order elastic constants of diamond

To determine the nonlinear elastic response of diamond, single crystals were shock compressed along the [100], [110], and [111] orientations to 120 GPa peak elastic stresses. Particle velocity histories and elastic wave velocities were measured by using laser interferometry. The measured elastic wave profiles were used, in combination with published acoustic measurements, to determine the complete set of third-order elastic constants. These constants represent the first experimental determination, and several differ significantly from those calculated by using theoretical models.

Plasmons in strongly coupled shock-compressed matter

We present the first measurements of the plasmon dispersion and damping in laser shock-compressed solid matter. Petawatt laser produced K-α radiation scatters on boron targets compressed by a 10 ns-long 400 J laser pulse. In the vicinity of the Fermi momentum, the scattering spectra show dispersionless, collisionally damped plasmons, indicating a strongly coupled electron liquid. These observations agree with calculations that include the Born-Mermin approximation to account for electron-ion collisional damping and local field corrections reflecting electron-electron correlations.

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.

High-Rate Plastic Deformation of Nanocrystalline Tantalum to Large Strains: Molecular Dynamics Simulation

Recent advances in the ability to generate extremes of pressure and temperature in dynamic experiments and to probe the response of materials has motivated the need for special materials optimized for those conditions as well as a need for a much deeper understanding of the behavior of materials subjected to high pressure and/or temperature. Of particular importance is the understanding of rate effects at the extremely high rates encountered in those experiments, especially with the next generation of laser drives such as at the National Ignition Facility. Here we use large-scale molecular dynamics (MD) simulations of the high-rate deformation of nanocrystalline tantalum to investigate the processes associated with plastic deformation for strains up to 100%. We use initial atomic configurations that were produced through simulations of solidification in the work of Streitz et al [Phys. Rev. Lett. 96, (2006) 225701]. These 3D polycrystalline systems have typical grain sizes of 10-20 nm. We also study a rapidly quenched liquid (amorphous solid) tantalum. We apply a constant volume (isochoric), constant temperature (isothermal) shear deformation over a range of strain rates, and compute the resulting stress-strain curves to large strains for both uniaxial and biaxial compression. We study the rate dependence and identify plastic deformation mechanisms. The identification of the mechanisms is facilitated through a novel technique that computes the local grain orientation, returning it as a quaternion for each atom. This analysis technique is robust and fast, and has been used to compute the orientations on the fly during our parallel MD simulations on supercomputers. We find both dislocation and twinning processes are important, and they interact in the weak strain hardening in these extremely fine-grained microstructures.