Search

The energetics of the regioselective mononitration of 9,10-BN-naphthalene with acetyl nitrate (H

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. We adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

The high-pressure equation of state (EOS) of energetic materials (EMs) is important for continuum and mesoscale models of detonation performance and initiation safety. Obtaining a high-fidelity EOS of the insensitive EM 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has proven to be difficult because of challenges in experimental characterization at high pressures (HPs). In this work, powder X-ray diffraction patterns were fitted using the recently discovered monoclinic I2/a phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic P1 phase. First-principles calculations were performed with Perdew-Burke-Ernzerhof (PBE) and PBE0 functionals including thermal effects using the P1 phase. PBE0 improves the description of hydrogen bonding and thus predicts accurate planar a and b lattice parameters under ambient conditions. However, discrepancies in the predicted lattice parameters above 4-10 GPa compared with experimental measurements indexed with P1 are further evidence of a structural modification at high pressure. Layer sliding defects are formed during molecular dynamics simulations, which induces an anharmonic effect on the thermal expansion of the c lattice parameter. In short, the results provide several insights into determining high-fidelity EOS parameters for TATB and other molecular crystals.

Mechanochemistry of glycine under compression and shear at room temperature is predicted using quantum-based molecular dynamics (QMD) and a simulation design based on rotational diamond anvil cell (RDAC) experiments. Ensembles of high throughput semiempirical density functional tight binding (DFTB) simulations are used to identify chemical trends and bounds for glycine chemistry during rapid shear under compressive loads of up to 15.6 GPa. Significant chemistry is found to occur during compressive shear above 10 GPa. Recovered products consist of small molecules such as water, structural analogs to glycine, heterocyclic molecules, large oligomers, and polypeptides including the simplest polypeptide glycylglycine at up to 4% mass fraction. The population and size of oligomers generally increases with pressure. A number of oligomeric polypeptide precursors and intermediates are also identified that consist of two or three glycine monomers linked together through C–C, C–N, and/or C–O bridges. Even larger oligomers also form that contain peptide C–N bonds and exhibit branched structures. Many of the product molecules exhibit one or more chiral centers. Our simulations demonstrate that athermal mechanical compressive shearing of glycine is a plausible prebiotic route to forming polypeptides., Compressive shearing forces can induce mechanochemical oligomerization reactions in glycine.

α-Glycine is studied up to 50 GPa using synchrotron angle-dispersive X-ray powder diffraction (XRD), Raman spectroscopy, and quantum chemistry calculations performed at multiples levels of theory. Results from both XRD and Raman experiments reveal an extended pressure stability of the α phase up to 50 GPa and the room temperature (RT) equation of state (EOS) was determined up to this pressure. This extended stability is corroborated by density functional theory (DFT) based calculations using the USPEX evolutionary structural search algorithm. Two calculated EOSs, as determined by DFT at T = 0 K and semiempirical density functional tight-binding (DFTB) at RT, and the calculated Raman modes frequencies show a good agreement with the corresponding experimental results. Our work provides a definitive phase diagram and EOS for α-glycine up to 50 GPa, which informs prebiotic synthesis scenarios that can involve pressures well in excess of 10 GPa.

We investigated the effects of hydrostatic pressure on α-glycylglycine (α-digly) using a combined experimental and theoretical approach. The results of powder X-ray diffraction show a change in compressibility of the axes above 6.7 GPa, but also indicate that the structure remains in the same monoclinic space group, suggesting an isosymmetric phase transition. A noticeable change in the Raman spectra between 6 and 7.5 GPa further supports the observed phase transition. First-principles-based calculations combined with the crystal structure prediction code USPEX predict a number of possible polymorphs at high pressure. An orthorhombic structure with a bent peptide backbone is the lowest enthalpy polymorph above 6.4 GPa; however, it is not consistent with experimental observations. A second monoclinic structure isosymmetric to α-digly, α'-digly, is predicted to become more stable above 11.4 GPa. The partial atomic charges in α'-digly differ from α-digly, and the molecule is bent, possibly indicating different reactivity of α'-digly. The similarity in the lattice parameters predicted from calculations and the axial changes observed experimentally support that the α'-digly phase is likely observed at high pressure. A possible explanation for the isosymmetric phase transition is discussed in terms of relaxing strained hydrogen bonding interactions. Such combined experimental and modeling efforts provide atomic-level insight into how pressure-driven conformational changes alter hydrogen-bonding networks in complicated molecular crystals.

Productivity from day one on supercomputers that leverage new technologies requires significant preparation. An institution that procures a novel system architecture often lacks sufficient institutional knowledge and skills to prepare for it. Thus, the "Center of Excellence" (CoE) concept has emerged to prepare for systems such as Summit and Sierra, currently the top two systems in the Top 500. This paper documents CoE experiences that prepared a workload of diverse applications and math libraries for a heterogeneous system. We describe our approach to this preparation, including our management and execution strategies, and detail our experiences with and reasons for using different programming approaches. Our early science and performance results show that the project enabled significant early seismic science with up to a l4X throughput increase over Cori. In addition to our successes, we discuss our challenges and failures so others may benefit from our experience.

We report gas-phase enthalpies of formation for the set of energetic molecules NTO, DADE, LLM-105, TNT, RDX, TATB, HMX, and PETN using the G2, G3, G4, and ccCA-PS3 quantum composite methods. Calculations for HMX and PETN hitherto represent the largest molecules attempted with these methods. G3 and G4 calculations are typically close to one another, with a larger difference found between these methods and ccCA-PS3. Although there is significant uncertainty in experimental values, the mean absolute deviation between the average experimental value and calculations are 12, 6, 7, and 3 kcal/mol for G2, G3, G4, and ccCA-PS3, respectively.

The propensities of the water self-ions, H3O(+) and OH(-), for the air-water interface have implications for interfacial acid-base chemistry. Despite numerous experimental and computational studies, no consensus has been reached on the question of whether or not H3O(+) and/or OH(-) prefer to be at the water surface or in the bulk. Here we report a molecular dynamics simulation study of the bulk vs interfacial behavior of H3O(+) and OH(-) that employs forces derived from density functional theory with a generalized gradient approximation exchange-correlation functional (specifically, BLYP) and empirical dispersion corrections. We computed the potential of mean force (PMF) for H3O(+) as a function of the position of the ion in the vicinity of an air-water interface. The PMF suggests that H3O(+) has equal propensity for the interface and the bulk. We compare the PMF for H3O(+) to our previously computed PMF for OH(-) adsorption, which contains a shallow minimum at the interface, and we explore how differences in solvation of each ion at the interface vs in the bulk are connected with interfacial propensity. We find that the solvation shell of H3O(+) is only slightly dependent on its position in the water slab, while OH(-) partially desolvates as it approaches the interface, and we examine how this difference in solvation behavior is manifested in the electronic structure and chemistry of the two ions.

Determining the unreacted equation of state of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) is challenging because it exhibits low crystal symmetry and low X-ray scattering strength. Here, we present the first high-pressure single-crystal X-ray diffraction (SXD) study of this material. Our SXD results reveal a previously unknown transition to a monoclinic phase above 4 GPa. No abrupt change of the volume occurs but the compressibility changes. Concomitant first principles evolutionary crystal structure prediction USPEX calculations confirm this transition and show that it involves a pressure-induced in-plane shift of the layers of TATB molecules with respect to the ambient-pressure phase.

While it is well established that the dominant degradation mechanism for the chemical warfare agent, sarin, is hydrolysis via an SN2 mechanism in aqueous environments, the same cannot be said when hydrophilic or hydrophobic surfaces are present. Utilizing first-principles molecular dynamics simulations based on density functional theory, potential degradation pathways for sarin were investigated in the vicinity of both hydrophilic and hydrophobic surfaces. In an idealized situation where sarin is within hydrogen-bonding distance to a hydrophilic surface, the most probable degradation occurred via an SN2 mechanism, but with a lowering of the reaction barrier (ΔΔG⧧) of ∼15 kcal/mol as compared with bulk aqueous SN2 decomposition. However, in the presence of an idealized hydrophobic surface, it was found that the degradation mechanism for sarin proceeded via an SN1-type mechanism where the phosphorus–fluorine bond breaks first, followed by a second step in which water acts as a nucleophile. The reaction barr...

A review is presented on recent progress of the application of molecular dynamics simulation methods with the inclusion of polarizability for the understanding of aqueous interfaces. Comparisons among a variety of models, including those based on density functional theory of the neat air-water interface, are given. These results are used to describe the effect of polarizability on modeling the microscopic structure of the neat air-water interface, including comparisons with recent spectroscopic studies. Also, the understanding of the contribution of polarization to the electrostatic potential across the air-water interface is elucidated. Finally, the importance of polarizability for understanding anion transfer across an organic-water interface is shown.

A series of 30 ps first principles molecular dynamics simulations in the microcanonical ensemble were carried out to investigate transport and vibrational properties of liquid water. To allow for sufficient sampling, the thermodynamic constraints were set to an elevated temperature of around 423 K and a density of 0.71 g cm(-)(3) corresponding to the saturated liquid density for the Becke-Lee-Yang-Parr (BLYP) representation of water. Four simulations using the Car-Parrinello molecular dynamics (CPMD) technique with varying values of the fictitious electronic mass (μ) and two simulations using the Born-Oppenheimer molecular dynamics (BOMD) technique are analyzed to yield structural and dynamical information. At the selected state point, the simulations are found to exhibit nonglassy dynamics and yield consistent results for the liquid structure and the self-diffusion coefficient, although the statistical uncertainties in the latter quantity are quite large. Consequently, it can be said that the CPMD and BOMD methods produce equivalent results for these properties on the time scales reported here. However, it was found that the choice of μ affects the frequency spectrum of the intramolecular modes, shifting them slightly to regions of lower frequency. Using a value of μ = 400 au results in a significant drift in the electronic kinetic energy of the system over the course of 30 ps and a downward drift in the ionic temperature. Therefore, for long trajectories at elevated temperatures, lower values of this parameter are recommended for CPMD simulations of water.

Whether aqueous interfaces are acidic or basic has implications for interfacial chemistry, but the question remains open. We employ first-principles molecular dynamics simulations to determine the intrinsic propensity of OH− for the air–water interface and find that OH− is stabilized by roughly kBT at the interface vs. the bulk. We predict, therefore, that the surface population OH− is slightly enhanced. Our simulations suggest that the solvation of OH− at the interface is similar to that observed in small water clusters, and they reveal changes in the orientation of solvating water molecules that are consistent with surface-sensitive vibrational spectra.

Efficient Monte Carlo algorithms are combined with the Quickstep energy routines of CP2K to develop a program that allows for Monte Carlo simulations in the canonical, isobaric-isothermal, and Gibbs ensembles using a first principles description of the physical system. Configurational-bias Monte Carlo techniques and pre-biasing using an inexpensive approximate potential are employed to increase the sampling efficiency and to reduce the frequency of expensive ab initio energy evaluations. The new Monte Carlo program has been validated through extensive comparison with molecular dynamics simulations using the programs CPMD and CP2K . Preliminary results for the vapor–liquid coexistence properties ( T = 473 K ) of water using the Becke–Lee–Yang–Parr exchange and correlation energy functionals, a triple-zeta valence basis set augmented with two sets of d-type or p-type polarization functions, and Goedecker–Teter–Hutter pseudopotentials are presented. The preliminary results indicate that this description of water leads to an underestimation of the saturated liquid density and heat of vaporization and, correspondingly, an overestimation of the saturated vapor pressure.

A series of first principles molecular dynamics and Monte Carlo simulations were carried out for liquid water to investigate the reproducibility of different sampling approaches. These simulations include Car−Parrinello molecular dynamics simulations using the program cpmd with different values of the fictitious electron mass in the microcanonical and canonical ensembles, Born−Oppenheimer molecular dynamics using the programs cpmd and cp2k in the microcanonical ensemble, and Metropolis Monte Carlo using cp2k in the canonical ensemble. With the exception of one simulation for 128 water molecules, all other simulations were carried out for systems consisting of 64 molecules. Although the simulations yield somewhat fortuitous agreement in structural properties, analysis of other properties demonstrate that one should exercise caution when assuming the reproducibility of Car−Parrinello and Born−Oppenheimer molecular dynamics simulations for small system sizes in the microcanonical ensemble. In contrast, the m...

In this work we investigate the fast anomalous diffusion of hydrogen molecules in water using Car–Parrinello molecular dynamics simulations. We employ Voronoi polyhedra analysis to distinguish between void diffusion and void hopping. Our results indicate that a combination of both mechanism is sufficient to explain anomalous diffusion. Furthermore, we investigate the geometrical and electronical structure of the first solvation shell.

We report dispersion-corrected density functional theoretical calculations of the unreacted equation of state (EOS) of crystal 2,6-diamino-3, 5-dinitropyrazine-1-oxide (LLM-105) under hydrostatic compression of up to 45 GPa. Convergence tests for k-points sampling in the Brillouin zone show that a 3 × 1 × 2 mesh is required to reproduce the X-ray crystal structure at ambient conditions, and we confirm our finding with a separate supercell calculation. Our high-pressure EOS yields a bulk modulus of 19.2 GPa, and indicates a tendency towards anisotropic compression along the b lattice vector due to molecular orientations within the lattice. We find that the electronic energy band gap decreases from a semiconductor type of 1.3 eV at 0 GPa to quasi-metallic type of 0.6 eV at 45 GPa. The extensive intermolecular hydrogen bonds involving the oxide (–NO) and dioxide (–NO2) interactions with the amine (–NH2) group showed enhanced interactions with increasing pressure that should be discernible in the mid IR spectral region. We do not find evidence for structural phase transitions or chemically induced transformations within the pressure range of our study. The gas phase heat of formation is calculated at the G4 level of theory to be 22.48 kcal/mol, while we obtain 25.92 kcal/mol using the ccCA-PS3 method. Density functional theory calculations of the crystal and the gas phases provided an estimate for the heat of sublimation of 32.4 kcal/mol. We thus determine the room-temperature solid heat of formation of LLM-105 to be −9.9 or −6.5 kcal/mol based on the G4 or ccCA-PS3 methods, respectively.

We present an ab initio molecular dynamics simulation of the aqueous liquid-vapor interface. Having successfully stabilized a region of bulk water in the center of a water slab, we were able to reproduce and further quantify the experimentally observed abundance of surface “acceptor-only”(19%) and “single-donor”(66%) moieties as well as substantial surface relaxation approaching the liquid-vapor interface. Examination of the orientational dynamics points to a faster relaxation in the interfacial region. Furthermore, the average value of the dipole decreases and the average value of the highest occupied molecular orbital for each water molecule increases approaching the liquid-vapor interface. Our results support the idea that the surface contains, on average, far more reactive states than the bulk.

Sliding clamps are toroidal proteins that encircle DNA and act as mobile platforms for DNA replication and repair machinery. To be loaded onto DNA, the eukaryotic sliding clamp Proliferating Cell Nuclear Antigen (PCNA) must be splayed open at one of the subunit-subunit interfaces by the ATP-dependent clamp loader, Replication Factor C, whose clamp-interacting sites form a right-handed spiral. Earlier molecular dynamics (MD) studies suggested that when PCNA opens, it preferentially adopts a right-handed spiral to match the spiral of the clamp loader. Here, analysis of considerably longer MD simulations shows that although the opened form of PCNA can achieve conformations matching the helical pitch of Replication Factor C, it is not biased toward a right-handed spiral structure. A coarse-grained elastic model was also built; its strong correspondence to the all-atom MD simulations of PCNA suggests that the behavior of the open clamp is primarily due to elastic deformation governed by the topology of the clamp domains. The elastic model was further used to construct the energy landscape of the opened PCNA clamp, including conformations that would allow binding to the clamp loader and loading onto double-stranded DNA. A picture of PCNA emerges of a rather flexible protein that, once opened, is mechanically compliant in the clamp opening process.

We utilize QM/MM methodology with different size QM regions to elucidate the spectral fingerprint of the pentagonal water cluster in Bacteriorhodopsin. This work was supported by the U.S. Department of Energy's (DOE) Office of Basic Energy Sciences, Chemical Sciences program. The Pacific Northwest National Laboratory is operated by Battelle for DOE.

We apply a quantum-mechanics/molecular-mechanics (QM/MM) scheme to simulate the migration of charged defects in alpha-quartz. A QM/MM framework together with a self-interaction correction (SIC) scheme was used to investigate the stability of the charged dimer defect and the E'(1) center in alpha-quartz. The use of a SIC scheme makes the charged dimer center unstable in alpha-quartz, in agreement with the experimental failure of identifying this defect, while the effects of the correction of the self-interaction error on the E-1(') center are negligible. Moreover, we conjecture that by overcoming a low energy barrier, the puckering mechanism (ideally heading from the charged dimer defect to the E'(1) center in SIC-free calculations) can be reiterated allowing the drift of the positive charge localized on an overcoordinated oxygen atom. This process can be regarded as an important channel of structural reorganization of oxygen deficient silica in the presence of strong polarizing electric fields.

Recent theoretical studies of 2,6-diamino-3,5-dinitropyrazine-1-oxide (C4H4N6O5 Lawrence Livermore Molecule No. 105, LLM-105) report unreacted high pressure equations of state that include several structural phase transitions, between 8 and 50 GPa, while one published experimental study reports equation of state (EOS) data up to a pressure of 6 GPa with no observed transition. Here we report the results of a synchrotron-based X-ray diffraction study and also ambient temperature isobaric-isothermal atomistic molecular dynamics simulations of LLM-105 up to 20 GPa. We find that the ambient pressure phase remains stable up to 20 GPa; there is no indication of a pressure induced phase transition. We do find a prominent decrease in b-axis compressibility starting at approximately 13 GPa and attribute the stiffening to a critical length where inter-sheet distance becomes similar to the intermolecular distance within individual sheets. The ambient temperature isothermal equation of state was determined through refinements of measured X-ray diffraction patterns. The pressure-volume data were fit using various EOS models to yield bulk moduli with corresponding pressure derivatives. We find very good agreement between the experimental and theoretically derived EOS.

The structure, dynamical, and electronic properties of liquid water utilizing different hybrid density functionals were tested within the plane wave framework of first-principles molecular dynamics simulations. The computational approach, which employs modified functionals with short-ranged Hartree-Fock exchange, was first tested in calculations of the structural and bonding properties of the water dimer and cyclic water trimer. Liquid water simulations were performed at the state point of 350 K at the experimental density. Simulations included three different hybrid functionals, a meta-functional, four gradient-corrected functionals, and the local density and Hartree-Fock approximations. It is found that hybrid functionals are superior in reproducing the experimental structure and dynamical properties as measured by the radial distribution function and self-diffusion constant when compared to the pure density functionals. The local density and Hartree-Fock approximations show strongly over- and understructured liquids, respectively. Hydrogen bond analysis shows that the hybrid functionals give slightly smaller average numbers of hydrogen bonds than pure density functionals but similar hydrogen bond populations. The average molecular dipole moments in the liquid from the three hybrid functionals are lower than those of the corresponding pure density functionals.

Efficient Monte Carlo algorithms and a mixed-basis set electronic structure program were used to compute from first principles the vapor−liquid coexistence curve of water. A water representation based on the Becke−Lee−Yang−Parr exchange and correlation functionals yields a saturated liquid density of 900 kg/m3 at 323 K and normal boiling and critical temperatures of 350 and 550 K, respectively. An analysis of the structural and electronic properties of the saturated liquid phase shows an increase of the asymmetry of the local hydrogen-bonded structure despite the persistence of a 4-fold coordination and decreases of the molecular dipole moment and of the spread of the lowest unoccupied molecular orbital with increasing temperature.

The predicted high pressure superionic phases of water and HF are investigated via ab initio molecular dynamics. These phases could potentially be achieved through either static compression with heating or through shock compression. We study water at densities of 2.0–3.0 g/cc (34 –115 GPa) along the 2000 K isotherm. We find that extremely rapid (superionic) diffusion of protons occurs in a fluid phase at pressures between 34 and 58 GPa. A transition to a stable body‐centered cubic (bcc) O lattice with superionic proton conductivity is observed between 70 and 75 GPa, a much higher pressure than suggested in prior work. We find that all molecular species at pressures greater than 75 GPa are too short lived to be classified as bound states. Up to 95 GPa, we find a solid superionic phase characterized by covalent O‐H bonding. Above 95 GPa, a transient network phase is found characterized by symmetric O‐H hydrogen bonding with nearly 50% covalent character. Ab initio molecular dynamics simulations of HF were conducted at densities of 1.8 – 4.0 g/cc along the 900 K isotherm. According to our simulations, a unique form of (symmetric) hydrogen bonding could play a significant role in superionic conduction. Our work shows that superionic phases could be more prevalent in hydrogen bonded systems than previously thought, such as HCl and HBr.

A series of first-principles Monte Carlo simulations in the isobaric-isothermal ensemble were carried out for liquid water at ambient conditions (T = 298 K and p = 1 atm). The Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and norm-conserving Goedecker-Teter-Hutter (GTH) pseudopotentials were employed with the CP2K simulation package to examine systems consisting of 64 water molecules. The fluctuations in the system volume encountered in simulations in the isobaric-isothermal ensemble require a reconsideration of the suitability of the typical charge-density cutoff and the regular grid-generation method previously used for the computation of the electrostatic energy in first-principles simulations in the microcanonical or canonical ensembles. In particular, it is noted that a much higher cutoff is needed and that the most computationally efficient method of creating grids con result in poor simulations. Analysis of the simulation trajectories using a very large charge-density cutoff at 1200 Ry and four different grid-generation methods point to a significantly underestimated liquid density of about 0.8 g cm(-3) resulting in a somewhat understructured liquid (with a value of about 2.7 for the height of the first peak in the oxygen-oxygen radial distribution function) for BLYP-GTH water at ambient conditions. In addition, a simulation using a charge-density cutoff at 280 Ry yields a higher density of 0.9 g cm(-3) showing the sensitivity of the simulation outcome to this parameter.

The predicted superionic phase of water is investigated via ab initio molecular dynamics at densities of 2.0--3.0 g/cc (34-115 GPa) along the 2000 K isotherm. We find that extremely rapid (superionic) diffusion of protons occurs in a fluid phase at pressures between 34 and 58 GPa. A transition to a stable body-centered cubic O lattice with superionic proton conductivity is observed between 70 and 75 GPa, a much higher pressure than suggested in prior work. We find that all molecular species at pressures greater than 75 GPa are too short lived to be classified as bound states. Up to 95 GPa, we find a solid superionic phase characterized by covalent O-H bonding. Above 95 GPa, a transient network phase is found characterized by symmetric O-H hydrogen bonding with nearly 50% covalent character. In addition, we describe a metastable superionic phase with quenched O disorder.

Raman spectroscopy in a laser heated diamond anvil cell and first principles molecular dynamics simulations have been used to study water in the temperature range 300 to 1500 K and at pressures to 56 GPa. We find a substantial decrease in the intensity of the O-H stretch mode in the liquid phase with pressure, and a change in slope of the melting line at 47 GPa and 1000 K. Consistent with these observations, theoretical calculations show that water beyond 50 GPa is ``dynamically ionized'' in that it consists of very short-lived ($l10\text{ }\text{ }\mathrm{fs}$) ${\mathrm{H}}_{2}\mathrm{O}$, ${\mathrm{H}}_{3}{\mathrm{O}}^{+}$, and ${\mathrm{OH}}^{\ensuremath{-}}$ species, and also that the mobility of the oxygen ions decreases abruptly with pressure, while hydrogen ions remain very mobile. We suggest that this regime corresponds to a superionic state.

We report an ab initio molecular dynamics simulation study of a complex composed of the superoxide anion and four water molecules. We probe the symmetry and stability of the superoxide tetrahydrate anion cluster at a low temperature (25 K) where there are vibrational spectra and a higher temperature (111 K) where experimental data are not available. From the ab initio simulation trajectory, we found the superoxide tetrahydrate complex to be highly symmetric, exhibiting a collective breathing motion on the picosecond time scale at low temperature that subsequently disappears at higher temperatures. Furthermore, we establish that superoxide tetrahydrate is a highly stable species at the higher temperature with little or no distortion to its geometry. The vibrational spectrum calculated from our ab initio trajectory is in reasonable agreement with experimental spectra, and we were able to assign individual peaks in the spectra to specific modes in the superoxide tetrahydrate anion cluster.

Computational studies of condensed-phase chemical reactions are challenging in part because of complexities in understanding the effects of the solvent environment on the reacting chemical species. Such studies are further complicated due to the demanding computational resources required to implement high-level ab initio quantum chemical methods when considering the solvent explicitly. Here, we use first-principles molecular dynamics simulations to examine condensed-phase decontamination reactions of V-type nerve agents in an explicit aqueous solvent. Our results include a detailed study of hydrolysis, base-hydrolysis, and nucleophilic oxidation of both VX and R-VX, as well as their protonated counterparts (i.e., VXH(+) and R-VXH(+)). The decontamination mechanisms and chemical reaction energy barriers, as determined from our simulations, are found to be in good agreement with experiment. The results demonstrate the applicability of using such simulations to assist in understanding new decontamination technologies or other applications that require computational screening of condensed-phase chemical reaction mechanisms.