Search

α-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.

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. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out‐gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro‐scale properties that can aid in the interpretation of future radiation damage experiments on plastic materials. [ABSTRACT FROM AUTHOR]

The response of rapidly compressed highly oriented pyrolytic graphite (HOPG) normal to its basal plane was investigated at a pressure of ∼80 GPa. Ultrafast x-ray diffraction using ∼100 fs pulses at the Materials Under Extreme Conditions sector of the Linac Coherent Light Source was used to probe the changes in crystal structure resulting from picosecond timescale compression at laser drive energies ranging from 2.5 to 250 mJ. A phase transformation from HOPG to a highly textured hexagonal diamond structure is observed at the highest energy, followed by relaxation to a still highly oriented, but distorted graphite structure following release. We observe the formation of a highly oriented lonsdaleite within 20 ps, subsequent to compression. This suggests that a diffusionless martensitic mechanism may play a fundamental role in phase transition, as speculated in an early work on this system, and more recent static studies of diamonds formed in impact events. [ABSTRACT FROM AUTHOR]

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. Our approach leverages a fine-grained united atom model for linear PDMS, which we extend to include crosslinking terms, and applies iterative Boltzmann inversion to obtain a coarse-grain "bead-spring-type" model. We then perform extensive molecular dynamics simulations to explore the effect of crosslinking on the rheology of silicone fluids, where we compute systematic increases in both density and shear viscosity that compare favorably to experiments that we conduct here. The kinematic viscosity of partially crosslinked fluids follows an empirical linear relationship that is surprisingly consistent with Rouse theory, which was originally derived for systems comprised of a uniform distribution of linear chains. The models developed here serve to enable quantitative bottom-up predictions for curing- and age-induced effects on macroscale rheological properties, allowing for accurate prediction of material properties based on fundamental chemical data. [ABSTRACT FROM AUTHOR]

There have been several numerical studies on the collapse of internal voids in energetic grains but fewer investigations have probed grain‐grain interface effects. In this study, we examine the effects of grain surface morphology and binder conditions on hot spot mechanisms during shock loading using the multi‐physics hydrocode, ALE3D, coupled with the thermochemical code, Cheetah. To improve the accuracy of our grain‐scale simulations, the HMX material models have been updated from previous studies to incorporate new property predictions from molecular dynamics simulations. In our simulations of the interface between two neighboring energetic grains, the upstream grain surface is described by a single sinusoid with amplitude (0.5–2 μm) and wavelength (0.5–2 μm). The effects of grain surface coating, i. e., no binder, partially coated, and fully coated with binder, is also explored. A range of shock loading pressures (10–30 GPa) is considered. Binder coating reduces hot spots at grain‐grain interfaces by limiting plastic deformation and heating of upstream grain surfaces as well as preventing transmission of localized deformation modes (e. g., jetting) to the downstream grain. More plastic work and heating is observed with new HMX material models versus the previous models because the strength does not thermally soften as much as a result of the higher melt temperature; this also reduces vortical flow and jetting mechanisms with the new HMX material models. A surrogate model based on neural networks is developed for the early stage reaction rate. Physics‐based input features improve the performance of the neural networks over basic input features alone with a root mean squared error=1.05 μs−1 and R2=0.98. Accurate and fast‐running surrogate models can effectively serve as structure‐property‐sensitivity relationships. [ABSTRACT FROM AUTHOR]

Ammonium perchlorate NH4ClO4 (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. A high-pressure barite-type structure (Phase II) forms at ≈4 GPa and appears stable up to 40 GPa. Refined atomic coordinates for Phase II are provided, and details for the Phase I → II transition mechanics are outlined. Pressure-dependent enthalpies computed for DFT-optimized crystal structures confirm the Phase I → II transition sequence, and the interpolated transition pressure is in excellent agreement with the experiment. Evidence for additional (underlying) structural modifications include a marked decrease in the Phase II b′-axis compressibility starting at 15 GPa and an unambiguous stress relaxation in the normalized stress-strain response at 36 GPa. Above 47 GPa, XRD Bragg peaks begin to decrease in amplitude and broaden. The apparent loss of crystalline long-range order likely signals the onset of amorphization. Three isostructural modifications were discovered within Phase II via Raman spectroscopy. A revised RT isothermal phase diagram is discussed based on the findings of this study. [ABSTRACT FROM AUTHOR]

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stressâ€"strain response of TATB, with a focus on understanding flow behavior. The dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stressâ€"strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure. [ABSTRACT FROM AUTHOR]

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. Each TATB molecule in the DPDE model is represented as an all-atom, rigid-molecule mesoparticle, with explicit external (molecular translational and rotational) DoF and coarse-grained implicit internal (vibrational) DoF. In addition to conserving linear and angular momentum, the DPDE equations of motion conserve the total system energy provided that particles can exchange energy between their external and internal DoF. The internal temperature of a TATB molecule is calculated using an internal equation of state, which we develop here, and the temperatures of the external and internal DoF are coupled using a fluctuation-dissipation relation. The DPDE force expression requires specification of the input parameter σ that determines the rate at which energy is exchanged between external and internal DoF. We adjusted σ based on the predictions for relaxation processes obtained from MD simulations. The parameterized DPDE model was employed in large-scale simulations of shock compression of TATB.We show that the rate of energy exchange governed by s can significantly influence the transient behavior of the system behind the shock. [ABSTRACT FROM AUTHOR]