Shock-Driven Decomposition of Polymers and Polymeric Foams.

Polymers and foams are pervasive in everyday life, as well as in specialized contexts such as space exploration, industry, and defense. They are frequently subject to shock loading in the latter cases, and will chemically decompose to small molecule gases and carbon (soot) under loads of sufficient strength. We review a body of work-most of it performed at Los Alamos National Laboratory-on polymers and foams under extreme conditions.

Shockwave compression and dissociation of ammonia gas.

We performed a series of plate impact experiments on NH gas initially at room temperature and at a pressure of ∼100 psi. Shocked states were determined by optical velocimetry and the temperatures by optical pyrometry, yielding compression ratios of ∼5-10 and second shock temperatures in excess of 7500 K. A first-principles statistical mechanical (thermochemical) approach that included chemical dissociation yielded reasonable agreement with experimental results on the principal Hugoniot, even with interparticle interactions neglected.

Reactive Monte Carlo sampling with an ab initio potential.

We present the first application of reactive Monte Carlo in a first-principles context. The algorithm samples in a modified NVT ensemble in which the volume, temperature, and total number of atoms of a given type are held fixed, but molecular composition is allowed to evolve through stochastic variation of chemical connectivity. We discuss general features of the method, as well as techniques needed to enhance the efficiency of Boltzmann sampling.

An efficient approach to ab initio Monte Carlo simulation.

We present a Nested Markov chain Monte Carlo (NMC) scheme for building equilibrium averages based on accurate potentials such as density functional theory. Metropolis sampling of a reference system, defined by an inexpensive but approximate potential, was used to substantially decorrelate configurations at which the potential of interest was evaluated, thereby dramatically reducing the number needed to build ensemble averages at a given level of precision. The efficiency of this procedure was maximized on-the-fly through variation of the reference system thermodynamic state (characterized here by its inverse temperature β(0)), which was otherwise unconstrained.

Intermolecular stabilization of 3,3'-diamino-4,4'-azoxyfurazan (DAAF) compressed to 20 GPa.

The room temperature stability of 3,3'-diamino-4,4'-azoxyfurazan (DAAF) has been investigated using synchrotron far-infrared, mid-infrared, Raman spectroscopy, and synchrotron X-ray diffraction (XRD) up to 20 GPa. The as-loaded DAAF samples exhibited subtle pressure-induced ordering phenomena (associated with positional disorder of the azoxy "O" atom) resulting in doubling of the a-axis, to form a superlattice similar to the low-temperature polymorph. Neither high pressure synchrotron XRD, nor high pressure infrared or Raman spectroscopies indicated the presence of structural phase transitions up to 20 GPa.

Nested Markov chain Monte Carlo sampling of a density functional theory potential: equilibrium thermodynamics of dense fluid nitrogen.

An optimized variant of the nested Markov chain Monte Carlo [n(MC)(2)] method [J. Chem. Phys.

Optimal sampling efficiency in Monte Carlo simulation with an approximate potential.

Building on the work of Iftimie et al. [J. Chem.

An "optimal" spawning algorithm for adaptive basis set expansion in nonadiabatic dynamics.

The full multiple spawning (FMS) method has been developed to simulate quantum dynamics in the multistate electronic problem. In FMS, the nuclear wave function is represented in a basis of coupled, frozen Gaussians, and a "spawning" procedure prescribes a means of adaptively increasing the size of this basis in order to capture population transfer between electronic states. Herein we detail a new algorithm for specifying the initial conditions of newly spawned basis functions that minimizes the number of spawned basis functions needed for convergence.

On the extent and connectivity of conical intersection seams and the effects of three-state intersections.

We discuss the connectivity of intersection spaces and the role of minimal energy points within these intersection spaces (minimal energy conical intersections or MECIs) in promoting nonadiabatic transitions. We focus on malonaldeyde as a specific example, where there is a low-lying three-state conical intersection. This three-state intersection is the global minimum on the bright excited electronic state, but it plays a limited role in population transfer in our ab initio multiple spawning (AIMS) simulations because the molecule must traverse a series of two-state conical intersections to reach the three-state intersection.

Optimizing conical intersections without derivative coupling vectors: application to multistate multireference second-order perturbation theory (MS-CASPT2).

We introduce a new method for optimizing minimal energy conical intersections (MECIs), based on a sequential penalty constrained optimization in conjunction with a smoothing function. The method is applied to optimize MECI geometries using the multistate formulation of second-order multireference perturbation theory (MS-CASPT2). Resulting geometries and energetics for conjugated molecules including ethylene, butadiene, stilbene, and the green fluorescent protein chromophore are compared with state-averaged complete active space self-consistent field (SA-CASSCF) and, where possible, benchmark multireference single- and double-excitation configuration interaction (MRSDCI) optimizations.

Ab initio molecular dynamics of excited-state intramolecular proton transfer using multireference perturbation theory.

We present the first calculations of excited-state dynamics using ab initio molecular dynamics with a multireference perturbation theory description of the electronic structure. The new AIMS-CASPT2 method is applied to a paradigmatic excited-state intramolecular proton-transfer reaction in methyl salicylate, and the results are compared with previous ultrafast spectroscopic experiments. Agreement of AIMS-CASPT2 and experimental results is quantitative.

Ab initio molecular dynamics of excited-state intramolecular proton transfer around a three-state conical intersection in malonaldehyde.

Excited-state potential energy surface (PES) characterization is carried out at the CASSCF and MRSDCI levels, followed by ab initio dynamics simulation of excited-state intramolecular proton transfer (ESIPT) on the S2(pipi*) state in malonaldehyde. The proton-transfer transition state lies close to an S2/S1 conical intersection, leading to substantial coupling of proton transfer with electronic relaxation. Proton exchange proceeds freely on S2, but its duration is limited by competition with twisting out of the molecular plane.

Competitive decay at two- and three-state conical intersections in excited-state intramolecular proton transfer.

We demonstrate the existence of a simultaneous degeneracy (not required by symmetry) of three electronic states in malonaldehyde. This is one of the first reports of such a triple degeneracy involving S0, S1, and S2 in a molecule with a closed-shell ground state. We further report on a two-state S2/S1 conical intersection which is higher in energy than the three-state intersection, but closer to the Franck-Condon point.