Optical phonons and elasticity of diamond at megabar stresses

Metallization of diamond

Experimental discovery of ultralarge elastic deformation in nanoscale diamond and machine learning of its electronic and phonon structures have created opportunities to address new scientific questions. Can diamond, with an ultrawide bandgap of 5.6 eV, be completely metallized, solely under mechanical strain without phonon instability, so that its electronic bandgap fully vanishes? Through first-principles calculations, finite-element simulations validated by experiments, and neural network learning, we show here that metallization/demetallization as well as indirect-to-direct bandgap transitions can be achieved reversibly in diamond below threshold strain levels for phonon instability. We identify the pathway to metallization within six-dimensional strain space for different sample geometries. We also explore phonon-instability conditions that promote phase transition to graphite. These findings offer opportunities for tailoring properties of diamond via strain engineering for electronic, photonic, and quantum applications.

First-Principles Calculation of Third-Order Elastic Constants via Numerical Differentiation of the Second Piola-Kirchhoff Stress Tensor

A general method is presented to calculate from first principles the full set of third-order elastic constants of a material of arbitrary symmetry. The method here illustrated relies on a plane-wave density functional theory scheme to calculate the Cauchy stress and the numerical differentiation of the second Piola-Kirchhoff stress tensor to evaluate the elastic constants. It is shown that finite difference formulas lead to a cancellation of the finite basis set errors, whereas simple solutions are proposed to eliminate numerical errors arising from the use of Fourier interpolation techniques. Applications to diamond, silicon, aluminum, magnesium, graphene, and a graphane conformer give results in excellent agreement with both experiments and previous calculations based on fitting energy density curves, demonstrating both the accuracy and generality of our new methodology to investigate nonlinear elastic behaviors of materials.

Transformation of shock-compressed graphite to hexagonal diamond in nanoseconds

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

Pressure Self-focusing Effect and Novel Methods for Increasing the Maximum Pressure in Traditional and Rotational Diamond Anvil Cells

The main principles of producing a region near the center of a sample, compressed in a diamond anvil cell (DAC), with a very high pressure gradient and, consequently, with high pressure are predicted theoretically. The revealed phenomenon of generating extremely high pressure gradient is called the pressure self-focusing effect. Initial analytical predictions utilized generalization of a simplified equilibrium equation. Then, the results are refined using our recent advanced model for elastoplastic material under high pressures in finite element method (FEM) simulations. The main points in producing the pressure self-focusing effect are to use beveled anvils and reach a very thin sample thickness at the center. We find that the superposition of torsion in a rotational DAC (RDAC) offers drastic enhancement of the pressure self-focusing effect and allows one to reach the same pressure under a much lower force and deformation of anvils.

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.

Thermo-mechanical stability and strength of peptide nanostructures from molecular dynamics: self-assembled cyclic peptide nanotubes

Peptide nanostructures present a wide range of opportunities for applications in biomedicine and bionanotechnology; hence experimental and theoretical studies aiming at determination of thermo-mechanical stability of peptide-based nanostructures are critical for the design and development of their technological applications. Here, we present a homogeneous deformation method combined with the finite elasticity theory and molecular dynamics simulations (MD) for the calculation of second-order anisotropic elastic constants for a membrane model made up of self-assembled cyclic peptide nanotubes. We have computed the values of all anisotropic elastic constants at 300 K. The value of the engineering Young's modulus (in the z direction) is 19.6 GPa. We observed a yield behavior in the z direction for a strain value of 6%. Furthermore, we also report calculated heat capacity, thermal expansion coefficient and isothermal compressibility of the system under study.

Nonlinear elasticity of monolayer graphene

By combining continuum elasticity theory and tight-binding atomistic simulations, we work out the constitutive nonlinear stress-strain relation for graphene stretching elasticity and we calculate all the corresponding nonlinear elastic moduli. Present results represent a robust picture on elastic behavior and provide the proper interpretation of recent experiments. In particular, we discuss the physical meaning of the effective nonlinear elastic modulus there introduced and we predict its value in good agreement with available data. Finally, a hyperelastic softening behavior is observed and discussed, so determining the failure properties of graphene.

Shock compressing diamond to a conducting fluid

Laser generated shock reflectance data show that diamond undergoes a continuous transition from optically absorbing to reflecting between Hugoniot pressures 600<P(H)<1000 GPa. The data are consistent with diamond having a thermal population of carriers at P(H) approximately 600 GPa, undergoing band overlap metallization at P(H) approximately 1000 GPa and melting at 800<P(H)<1000 GPa. The results agree well with an equation of state model that predicts that elemental carbon remains solid throughout the interior of Neptune.

Properties of diamond under hydrostatic pressures up to 140 GPa

Diamond is the archetypal covalent material. Each atom in an sp(3) configuration is bonded to four nearest neighbours. Because of its remarkable properties, diamond has been extensively studied. And yet our knowledge of the properties of diamond under very high pressure is still incomplete. Although diamond is known to be the preferred allotrope of carbon at high pressure, the possibility of producing under pressure high-density polymorphs of diamond, including metallic forms, has been discussed. Structural changes have already been reported in diamond under non-hydrostatic pressures around 150 GPa and large deformation. However, measurements of the properties of diamond under hydrostatic pressure have been limited to below 40 GPa. Here, we report accurate measurements of the volume and of the optical phonon frequency of diamond under hydrostatic pressure up to 140 GPa. We show that diamond is more compressible than currently expected. By combining the volume and the frequency pressure shifts, we deduce that diamond remains very stable under pressure: it is a Gruneisen solid up to at least 140 GPa, and the covalent bond is even slightly strengthened under pressure. Finally, the optical phonon frequency versus pressure is calibrated here to be used as a pressure gauge for diamond anvil cell studies in the multi-megabar range.

Spectroscopic studies of the vibrational and electronic properties of solid hydrogen to 285 GPa

We report Raman scattering and visible to near-infrared absorption spectra of solid hydrogen under static pressure up to 285 GPa between 20 and 140 K. We obtain pressure dependences of vibron and phonon modes consistent with results previously determined to lower pressures. The results indicate the stability of the ordered molecular phase III to the highest pressure reached and provide constraints on the insulator-to-metal transition pressure.

Instabilities in diamond under high shear stress

We investigate, through first-principles calculations, lattice instabilities induced in diamond by the application of high shear stresses. For shear stresses as low as 95 GPa a lattice instability will occur, leading to graphitelike layered structures. This effect is highly anisotropic. The reversal of the direction of the applied shear forces may cause a change of 80 GPa in the shear stress value at which the instability develops. The same reversal also causes different bonds to be broken, resulting in a drastic change in the orientation of the resulting graphitelike structures. We also find that an additional compressive stress of 50 GPa along the (111) direction does not eliminate the shear-induced instability.