Optical Characterization of Shock-Induced Chemistry in the Explosive Nitromethane using DFT and Time-Dependent DFT

  • PDF / 6,305,565 Bytes
  • 6 Pages / 612 x 792 pts (letter) Page_size
  • 61 Downloads / 143 Views

DOWNLOAD

REPORT


Optical Characterization of Shock-Induced Chemistry in the Explosive Nitromethane using DFT and Time-Dependent DFT Lenson A. Pellouchoud, Evan J. Reed Stanford University Materials Science & Engineering 496 Lomita Mall Stanford, CA 94024, U.S.A. ABSTRACT A thorough understanding of chemistry in extreme environments is a major challenge in experimental as well as theoretical work. With continual improvements in ultrafast optical measurements and new methods for simulations of shock-induced chemistry for timescales approaching a nanosecond, the opportunity is beginning to exist to connect experiments with simulations on the same timescale. In the present work, we compute the optical properties of the energetic material nitromethane (CH3NO2) for the first 100 picoseconds behind the detonation shock front in a molecular dynamics simulation. We compute optical spectra using the KuboGreenwood approach with DFT Kohn-Sham electronic states and compare with spectra computed by linear-response time-dependent density functional theory (TDDFT). The latter typically yields more accurate spectra for molecular systems. At optical wavelengths, the TDDFT method offers a correction of up to 25% in the real part of conductivity relative to the Kubo-Greenwood calculation. We also study the effects of thermal electronic excitations on the calculated spectra, and find no discernible change at optical wavelengths. In all of our calculations, we observe a non-monotonic change over time in the entire spectrum of optical properties as decomposition products evolve. The most optically relevant decomposition products are found to be NO, CNO, CNOH, water, and larger transient molecules. In particular, the disappearance of transient NO and CNO molecules (about 90 picoseconds behind the shock front) is coincident with a substantial decrease in conductivity across the optical spectrum. INTRODUCTION Shock-induced chemistry is an active and accelerating area of research in simulations as well as experiments1. In this work, we compute measurable optical properties of the initial stages of shock-induced chemistry in liquid nitromethane using density-functional theory-based approaches. Our calculations utilize the results of tight-binding-based multi-scale shock technique (MSST)2-5 molecular dynamics (performed in earlier work) for the first 100 picoseconds behind the shock front. The MSST simulation used 16 nitromethane molecules (112 atoms) with forces from self-consistent-charge density-functional tight binding (SCC-DFTB)6. To our knowledge, this is the first ab-initio investigation of optical properties of an energetic material shocked to conditions near those of detonation. In this work, we identify major chemical products over time and employ DFT-based Kubo-Greenwood perturbation theory and TDDFT to compute time-dependent optical properties of the mixture. The Kubo-Greenwood approach with Kohn-Sham states is often employed for

the study of matter at high temperatures7-12, but TDDFT approaches are typically found to yield more accurate optical s