Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog.

An atomic-level understanding of radiation-induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out-gassing products are computed.

Predicted Fracture Tendency of Naturally Occurring Aluminum Surface Coatings under Tensile Loading.

Naturally occurring coatings on aluminum metal, such as its oxide or hydroxide, serve to protect the material from corrosion. Understanding the conditions under which these coatings mechanically fail is therefore expected to be an important aspect of predictive models for aluminum component lifetimes. To this end, we develop and apply a molecular dynamics (MD) modeling framework for conducting tension tests that is capable of isolating factors governing the mechanical strength as a function of coating chemistry, defect morphology, and variables associated with the loading path.

Maximum Entropy Theory of Multiscale Coarse-Graining via Matching Thermodynamic Forces: Application to a Molecular Crystal (TATB).

The MSCG/FM (multiscale coarse-graining via force-matching) approach is an efficient supervised machine learning method to develop microscopically informed coarse-grained (CG) models. We present a theory based on the principle of maximum entropy (PME) enveloping the existing MSCG/FM approaches. This theory views the MSCG/FM method as a special case of matching the thermodynamic forces from the extended ensemble described by the set of thermodynamic (relevant) system coordinates.

Model for Humidity-Mediated Diffusion on Aluminum Surfaces and Its Role in Accelerating Atmospheric Aluminum Corrosion.

Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity.

United atom and coarse grained models for crosslinked polydimethylsiloxane with applications to the rheology of silicone fluids.

Siloxane systems consisting primarily of polydimethylsiloxane (PDMS) are versatile, multifaceted materials that play a key role in diverse applications. However, open questions exist regarding the correlation between their varied atomic-level properties and observed macroscale features. To this effect, we have created a systematic workflow to determine coarse-grained simulation models for crosslinked PDMS in order to further elucidate the effects of network changes on the system's rheological properties below the gel point.

Multiscale Strategy for Predicting Radiation Chemistry in Polymers.

A primary mode for radiation damage in polymers arises from ballistic electrons that induce electronic excitations, yet subsequent chemical mechanisms are poorly understood. We develop a multiscale strategy to predict this chemistry starting from subatomic scattering calculations. Nonadiabatic molecular dynamics simulations sample initial bond-breaking events following the most likely excitations, which feed into semiempirical simulations that approach chemical equilibrium.

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy.

Regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive.

Polymer degradation through chemical change: a quantum-based test of inferred reactions in irradiated polydimethylsiloxane.

Chemical reaction schemes are key conceptual tools for interpreting the results of experiments and simulations, but often carry implicit assumptions that remain largely unverified for complicated systems. Established schemes for chemical damage through crosslinking in irradiated silicone polymers comprised of polydimethylsiloxane (PDMS) date to the 1950's and correlate small-molecule off-gassing with specific crosslink features. In this regard, we use a somewhat reductionist model to develop a general conditional probability and correlation analysis approach that tests these types of causal connections between proposed experimental observables to reexamine this chemistry through quantum-based molecular dynamics (QMD) simulations.

Mechanochemical synthesis of glycine oligomers in a virtual rotational diamond anvil cell.

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.

A Hotspot's Better Half: Non-Equilibrium Intra-Molecular Strain in Shock Physics.

Shockwave interactions with a material's microstructure localizes energy into hotspots, which act as nucleation sites for complex processes such as phase transformations and chemical reactions. To date, hotspots have been described via their temperature fields. Nonreactive, all-atom molecular dynamics simulations of shock-induced pore collapse in a molecular crystal show that more energy is localized as potential energy (PE) than can be inferred from the temperature field and that PE localization persists beyond thermal diffusion.

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105).

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.

A Quantum-Based Approach to Predict Primary Radiation Damage in Polymeric Networks.

Initial atomistic-level radiation damage in chemically reactive materials is thought to induce reaction cascades that can result in undesirable degradation of macroscale properties. Ensembles of quantum-based molecular dynamics (QMD) simulations can accurately predict these cascades, but extracting chemical insights from the many underlying trajectories is a labor-intensive process that can require substantial intuition. We develop here a general and automated graph-based approach to extract all chemically distinct structures sampled in QMD simulations and apply our approach to predict primary radiation damage of polydimethylsiloxane (PDMS), the main constituent of silicones.

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity.

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 2/ phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic 1̅ phase.

High Explosive Ignition through Chemically Activated Nanoscale Shear Bands.

Shock initiation and detonation of high explosives is considered to be controlled through hot spots, which are local regions of elevated temperature that accelerate chemical reactions. Using classical molecular dynamics, we predict the formation of nanoscale shear bands through plastic failure in shocked 1,3,5-triamino-2,4,6-trinitrobenzene high explosive crystal. By scale bridging with quantum-based molecular dynamics, we show that shear bands exhibit lower reaction barriers.

An Isosymmetric High-Pressure Phase Transition in α-Glycylglycine: A Combined Experimental and Theoretical Study.

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.

Anisotropic Hydrolysis Susceptibility in Deformed Polydimethylsiloxanes.

Chemical reactions involving the polydimethylsiloxane (PDMS) backbone can induce significant network rearrangements and ultimately degrade macro-scale mechanical properties of silicone components. Using two levels of quantum chemical theory, we identify a possible electronic driver for chemical susceptibility in strained PDMS chains and explore the complicated interplay between hydrolytic chain scissioning reactions, mechanical deformations of the backbone, water attack vector, and chain mobility. Density functional theory (DFT) calculations reveal that susceptibility to hydrolysis varies significantly with the vector for water attacks on silicon backbone atoms, which matches strain-induced anisotropic changes in the backbone electronic structure.

Synthesis of functionalized nitrogen-containing polycyclic aromatic hydrocarbons and other prebiotic compounds in impacting glycine solutions.

Proteinogenic amino acids can be produced on or delivered to a planet impacting abiotic sources and consequently were likely present before the emergence of life on Earth. However, the role that these materials played in prebiotic scenarios remains an open question, in part because little is known about the survivability and reactivity of astrophysical organic compounds upon impact with a planetary surface. To this end, we use a force-matched semi-empirical quantum simulation method to study impacts of aqueous proteinogenic amino acids at conditions reaching 48 GPa and 3000 K.

Chemical Degradation Pathways in Siloxane Polymers Following Phenyl Excitations.

We use ensembles of quantum-based molecular dynamics simulations to predict the chemical reactions that follow radiation-induced excitations of phenyl groups in a model copolymer of polydimethylsiloxane and polydiphenylsiloxane. Our simulations span a wide range of highly porous and condensed phase densities and include both wet and dry conditions. We observe that in the absence of water, excited phenyl groups tend to abstract hydrogen from other methyl or phenyl side groups to produce benzene, with the under-hydrogenated group initiating subsequent intrachain cyclization reactions.

Role of filler and its heterostructure on moisture sorption mechanisms in polyimide films.

Moisture sorption and diffusion exacerbate hygrothermal aging and can significantly alter the chemical and mechanical properties of polymeric-based components over time. In this study, we employ a multi-pronged multi-scale approach to model and understand moisture diffusion and sorption processes in polyimide polymers. A reactive transport model with triple-mode sorption (i.

Effects of pressure on the structure and lattice dynamics of ammonium perchlorate: A combined experimental and theoretical study.

Ammonium perchlorate NHClO (AP) was studied using synchrotron angle-dispersive X-ray powder diffraction (XRPD) and Raman spectroscopy. A diamond-anvil cell was used to compress AP up to 50 GPa at room temperature (RT). Density functional theory (DFT) calculations were performed to provide further insight and comparison to the experimental data.

Generating Converged Accurate Free Energy Surfaces for Chemical Reactions with a Force-Matched Semiempirical Model.

We demonstrate the capability of creating robust density functional tight binding (DFTB) models for chemical reactivity in prebiotic mixtures through force matching to short time scale quantum free energy estimates. Molecular dynamics using density functional theory (DFT) is a highly accurate approach to generate free energy surfaces for chemical reactions, but the extreme computational cost often limits the time scales and range of thermodynamic states that can feasibly be studied. In contrast, DFTB is a semiempirical quantum method that affords up to a thousandfold reduction in cost and can recover DFT-level accuracy.

Characteristics of energy exchange between inter- and intramolecular degrees of freedom in crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) with implications for coarse-grained simulations of shock waves in polyatomic molecular crystals.

In this report, we characterize the kinetics and dynamics of energy exchange between intramolecular and intermolecular degrees of freedom (DoF) in crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). All-atom molecular dynamics (MD) simulations are used to obtain predictions for relaxation from certain limiting initial distributions of energy between the intra- and intermolecular DoF. The results are used to parameterize a coarse-grained Dissipative Particle Dynamics at constant Energy (DPDE) model for TATB.

Theoretical determination of anisotropic thermal conductivity for initially defect-free and defective TATB single crystals.

The anisotropic thermal conductivity was determined for initially defect-free and defective crystals of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a material that exhibits a graphitic-like packing structure with stacked single-molecule-thick layers, using the reverse non-equilibrium molecular dynamics method and an established TATB molecular dynamics force field. Thermal conduction in TATB is predicted to be substantially higher and more anisotropic than in other related organic molecular explosives, with conduction along directions nominally in the plane of the molecular layers at least 68% greater than conduction along the direction exactly perpendicular to the layers. Finite-size effects along the conduction directions were assessed.

Theoretical determination of anisotropic thermal conductivity for crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB).

Bond stretching and three-center angle bending potentials have been developed to extend an existing rigid-bond 1,3,5-triamino-2,4,6-trinitrobenzene molecular dynamics force field [D. Bedrov, O. Borodin, G.