Publications

Scripta Materialia

Ion implantation emerges as a promising technique to address the persistent challenge of lithium (Li) filament growth in solid-state electrolytes as it can induce compressive stresses inhibiting crack growth and deflect dendrites, de facto mitigating early electrolyte failure. In this study, we examine the potential paradox of ion implantation: while aiming to enhance electrolyte performance, the radiation damage associated with implantation might inadvertently compromise both the ionic conductivity and the intrinsic fracture toughness of the material, rendering the material unsuitable for battery applications. Specifically, we employed molecular dynamics simulations to examine the scope of the downsides of ion implantation, specifically: (i) reduced ionic conductivity (due to radiation-induced amorphization) and (ii) mechanical stability (due to radiation-induced embrittlement) in ion-implanted Li₇La₃Zr₂O₁₂ (LLZO) solid-state electrolytes. We explore how radiation damage impacts LLZO’s crystalline structure, Li-ion diffusion, and fracture properties at various temperatures and radiation damage levels. The study aims to provide insights into the competing effects of ion implantation and suggest potential engineering strategies for developing more robust solid-state electrolytes with improved conductivity and dendrite resistance.

The creep behavior of 316 L stainless steel at room temperature was evaluated as a function of time and applied stress using a new high-Throughput approach. Several common creep models were evaluated against the observations, leading to deeper analysis of a stress-dependent modified logarithmic creep model. Within this model, multiple sources of uncertainty were compared. Aleatoric stochastic variation between samples under nominally identical conditions was identified as the primary contributor to uncertainty in creep response. Under any particular set of conditions, the sample-To-sample variability in creep strain was as high as a factor of two, highlighting the engineering importance of characterizing large statistical datasets. The model's extrapolation capabilities were assessed by comparing predictions derived from calibration on partial, shorter-duration subsets of the data. These findings underscore the importance of accounting for stochastic effects in predictive modeling of aging phenomena.

Sputter-deposited, nanocrystalline Cu-Ag thin films produced across a broad compositional and deposition-parameter space were evaluated to unravel the process-structure-property relationships important for creating hard, conductive electrical contacts and coatings. Combinatorial deposition involving pulsed direct current magnetron sputtering of elemental targets enabled swift examination of nearly the full range of alloy compositions and a relevant portion of deposition atomistics. Several high-throughput characterization modalities were employed to evaluate the chemistry, structure, and properties of the films. The resultant hardness, modulus, film density, crystal texture, and resistivity were analyzed in terms of key deposition characteristics (incident atom kinetic energy and incidence angle) predicted by binary-collision, kinematic Monte Carlo simulations. The study revealed improved hardness, parabolic resistivity dependence on composition, and compositional and process dependencies of film tarnishing. The results are discussed in the context of variations in microstructure and film density. Transmission electron microscopy and X-ray diffraction demonstrate several forms of compositional variation including solute segregation to grain boundaries as well as periodic, intragranular compositional modulations. Annealing of a Cu-rich alloy film exhibiting grain boundary segregation showed that this as-deposited, compositional variation is not stable above 100 °C. Finally, the Cu-Ag system is shown to have potential for hard, conductive, tarnish-resistant and room temperature-stable nanocrystalline thin films across the composition space.

We introduce physics-informed multimodal autoencoders (PIMA)-a variational inference framework for discovering shared information in multimodal datasets. Individual modalities are embedded into a shared latent space and fused through a product-of-experts formulation, enabling a Gaussian mixture prior to identify shared features. Sampling from clusters allows cross-modal generative modeling, with a mixture-of-experts decoder that imposes inductive biases from prior scientific knowledge and thereby imparts structured disentanglement of the latent space. This approach enables cross-modal inference and the discovery of features in high-dimensional heterogeneous datasets. Consequently, this approach provides a means to discover fingerprints in multimodal scientific datasets and to avoid traditional bottlenecks related to high-fidelity measurement and characterization of scientific datasets.

CrMnFeCoNi, also called the Cantor alloy, is a well-known high-entropy alloy whose magnetic properties have recently become a focus of attention. Here, we present a detailed muon spin relaxation study of the influence of chemical composition and sample processing protocols on the magnetic phase transitions and spin dynamics of several different Cantor alloy samples. Specific samples studied include a pristine equiatomic sample, samples with deficient and excess Mn content, and equiatomic samples magnetized in a field of 9 T or plastically deformed in pressures up to 0.5 GPa. The results confirm the sensitive dependence of the transition temperature on composition and demonstrate that post-synthesis pressure treatments cause the transition to become significantly less homogeneous throughout the sample volume. In addition, we observe critical spin dynamics in the vicinity of the transition in all samples, reminiscent of canonical spin glasses and magnetic materials with ideal continuous phase transitions. Application of an external magnetic field suppresses the critical dynamics in the Mn-deficient sample, while the equiatomic and Mn-rich samples show more robust critical dynamics. The spin-flip thermal activation energy in the paramagnetic phase increases with Mn content, ranging from 3.1⁢(3) × 10^-21 J for 0% Mn to 1.2⁢(2) × 10^-20 J for 30% Mn content. These results shed light on critical magnetic behavior in environments of extreme chemical disorder and demonstrate the tunability of spin dynamics in the Cantor alloy via chemical composition and sample processing.

Abstract not provided.

Abstract not provided.

BeyondFingerprinting was a 2021-2024 Sandia Grand Challenge LDRD exploring the potential to develop new resilient materials and manufacturing processes by taking an artificial-intelligence (AI)-guided approach that integrates human-subject-matter expertise with algorithms enriched with physics-based constraints to unearth process-structure-property correlations. Such algorithms, trained on high-throughput experiments and simulations, are shown to serve as surrogate models that efficiently detect key “fingerprints” in materials data, prognose material performance, and guide effective process improvements. To accelerate broader adoption across mission areas, this AI-guided approach was demonstrated with three complex process-centric exemplars: electroplating, physical vapor deposition, and laser powder bed fusion. Together, these exemplars impact nearly every hardware component relevant to DOE and NNSA national security missions.

Sputter-deposited Pt-Au thin films have been reported to develop a hard, stable, nanocrystalline structure, yet little is known about how these characteristics vary with Pt_xAu_1-x composition and process conditions. Toward this end, this document describes an extensive, combinatorial Pt-Au thin film library including characterized film compositions, structure, and properties. Complemented by kinematic Monte Carlo simulations of codeposition, a broad range of Pt_xAu_1-x compositions (from x ~ 0.02 to 0.93) was first established by sputtering with varied magnetron powers and gun tilt angles. Further, the produced films were subsequently interrogated using automated nanoindentation, x-ray reflectivity, x-ray diffraction, atomic force microscopy, surface profilometry, four-point probe sheet resistance techniques, and wavelength dispersive spectroscopy in order to determine how hardness, modulus, density, surface roughness, structure, and resistivity vary with film stoichiometry and process parameters. Combinatorial films displayed an assortment of properties with the hardness of some films exceeding values reported previously for this material system. High hardness, high modulus, and low resistivity were generally attained when using increased deposition energy and reduced angle-of-incidence processes. Overall, the research identified promising, new Pt_xAu_1-x compositions for future study and pinpointed strategies for improved deposition.

Here, this work focuses on the representation of model-form uncertainties in molecular dynamics simulations in various statistical ensembles. In prior contributions, the modeling of such uncertainties was formalized and applied to quantify the impact of, and the error generated by, pair-potential selection in the microcanonical ensemble (NVE). In this work, we extend this formulation and present a linear-subspace reduced-order model for the canonical (NVT) and isobaric (NPT) ensembles. The symplectic reduced-order basis is randomized on the tangent space of the Stiefel manifold to provide topological relationships and capture model-form uncertainty. Using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), we assess the relevance of these stochastic reduced-order atomistic models on canonical problems involving a Lennard-Jones fluid and an argon crystal melt.

Abstract not provided.

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

We present a comprehensive benchmarking framework for evaluating machine-learning approaches applied to phase-field problems. This framework focuses on four key analysis areas crucial for assessing the performance of such approaches in a systematic and structured way. Firstly, interpolation tasks are examined to identify trends in prediction accuracy and accumulation of error over simulation time. Secondly, extrapolation tasks are also evaluated according to the same metrics. Thirdly, the relationship between model performance and data requirements is investigated to understand the impact on predictions and robustness of these approaches. Finally, systematic errors are analyzed to identify specific events or inadvertent rare events triggering high errors. Quantitative metrics evaluating the local and global description of the microstructure evolution, along with other scalar metrics representative of phase-field problems, are used across these four analysis areas. This benchmarking framework provides a path to evaluate the effectiveness and limitations of machine-learning strategies applied to phase-field problems, ultimately facilitating their practical application.

Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

This dataset is comprised of a library of atomistic structure files and corresponding X-ray diffraction (XRD) profiles and vibrational density of states (VDoS) profiles for bulk single crystal silicon (Si), gold (Au), magnesium (Mg), and iron (Fe) with and without disorder introduced into the atomic structure and with and without mechanical loading. Included with the atomistic structure files are descriptor files that measure the stress state, phase fractions, and dislocation content of the microstructures. All data was generated via molecular dynamics or molecular statics simulations using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. This dataset can inform the understanding of how local or global changes to a materials microstructure can alter their spectroscopic and diffraction behavior across a variety of initial structure types (cubic diamond, face-centered cubic (FCC), hexagonal close-packed (HCP), and body-centered cubic (BCC) for Si, Au, Mg, and Fe, respectively) and overlapping changes to the microstructure (i.e., both disorder insertion and mechanical loading).

Abstract not provided.

This is a presentation on using multi-objective bayesian optimization and machine methods to guide the optimization of metal alloy compositions.

Abstract not provided.

Understanding the mobility of twinning dislocations is important for multiscale modeling of crystal plasticity, especially at high strain rates, where such dislocations may reach transonic or supersonic speeds. Here, we used molecular dynamics simulations to investigate the relationship between dislocation velocity and the applied resolved shear stress of an edge twinning dislocation in copper up to supersonic speeds. The twinning dislocation mobility relation is composed of two branches separated by a band of forbidden velocities. The lower velocity branch is limited by the first transverse sound speed ~2000 m/s while the upper branch stretches from ~3500 m/s in the transonic regime to supersonic velocities. Twinning dislocations cannot undergo uniform steady-state motion at velocities within the forbidden band. Our simulation results also reveal that edge twinning dislocation motion in copper is kink-mediated. We discuss the implications of our findings for the motion of twins, twinning dislocations, and twinning dislocation kinks in copper.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.