Search

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.

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. Here we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.

Interlayer interactions in layered transition metal dichalcogenides are known to be important for describing their electronic properties. Here, we demonstrate that the absence of interlayer coupling in monolayer VTe 2 also causes their structural modification from a distorted 1T' structure in bulk and multilayer samples to a hexagonal 1T structure in the monolayer. X-ray photoemission spectroscopy indicates that this structural transition is associated with electron transfer from the vanadium d bands to the tellurium atoms for the monolayer. This charge transfer may reduce the in-plane d orbital hybridization and thus favor the undistorted 1T structure. Phonon-dispersion calculations show that, in contrast to the 1T' structure, the 1T structure exhibits imaginary phonon modes that lead to a charge density wave (CDW) instability, which is also observed by low-temperature scanning tunneling microscopy as a 4 × 4 periodic lattice distortion. Thus, this work demonstrates a novel CDW material, whose properties are tuned by interlayer interactions.