Unsupervised Learning-Based Multiscale Model of Thermochemistry in 1,3,5-Trinitro-1,3,5-triazinane (RDX).

The response of high-energy-density materials to thermal or mechanical insults involves coupled thermal, mechanical, and chemical processes with disparate temporal and spatial scales that no single model can capture. Therefore, we developed a multiscale model for 1,3,5-trinitro-1,3,5-triazinane, RDX, where a continuum description is informed by reactive and nonreactive molecular dynamics (MD) simulations to describe chemical reactions and thermal transport. Reactive MD simulations under homogeneous isothermal and adiabatic conditions are used to develop a reduced-order chemical kinetics model. Coarse graining is done using unsupervised learning via non-negative matrix factorization. Importantly, the components resulting from the analysis can be interpreted as reactants, intermediates, and products, which allows us to write kinetics equations for their evolution. The kinetics parameters are obtained from isothermal MD simulations over a wide temperature range, 1200-3000 K, and the heat evolved is calibrated from adiabatic simulations. We validate the continuum model against MD simulations by comparing the evolution of a cylindrical hotspot 10 nm in diameter. We find excellent agreement in the time evolution of the hotspot temperature fields both in cases where quenching is observed and at higher temperatures for which the hotspot transitions into a deflagration wave. The validated continuum model is then used to assess the criticality of hotspots involving scales beyond the reach of atomistic simulations that are relevant to detonation initiation.

Engineering curvature in graphene ribbons using ultrathin polymer films.

We propose a method to induce curvature in graphene nanoribbons in a controlled manner using an ultrathin thermoset polymer in a bimaterial strip setup and test it via molecular dynamics (MD) simulations. Continuum mechanics shows that curvature develops to release the residual stress caused by the chemical and thermal shrinkage of the polymer during processing and that this curvature increases with decreasing film thickness; however, significant deformation is only achieved for ultrathin polymer films. Quite surprisingly, explicit MD simulations of the curing and annealing processes show that the predicted trend not just continues down to film thicknesses of 1-2 nm but that the curvature development is enhanced significantly in such ultrathin films due to surface tension effects. This combination of effects leads to very large curvatures of over 0.14 nm(-1) that can be tuned via film thickness. This provides a new avenue to engineer curvature and, thus, electromagnetic properties of graphene.

Avalanches and scaling in plastic deformation.

Plastic deformation of crystalline materials is a complex nonhomogeneous process characterized by avalanches in the motion of dislocations. We study the evolution of dislocation loops using an analytically solvable phase-field model of dislocations for ductile single crystals during monotonic loading. The distribution of dislocation loop sizes is given by P(A) approximately A-sigma, with sigma=1.8+/-0.1. The exponent is in agreement with those found in acoustic emission experiments. This model also predicts a range of macroscopic behaviors in agreement with observation, including hardening with monotonic loading, and a maximum in the acoustic emission signal at the onset of yielding.

Dislocation structures and the deformation of materials.

We present results from phase-field simulations of a two-dimensional model of dislocation microstructure development under increasing strain that incorporates the effects of the full, three-dimensional, microstructure in an approximate way. Despite its simplicity, the model yields quantitative predictions of both the deformation properties of face-centered cubic metals as well as key descriptors of the evolving microstructure over a wide range of stress and strain. The present results have important implications for how we interpret and describe the deformation properties of fcc materials.