Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Geometrically frustrated interactions drive structural complexity in amorphous calcium carbonate

Amorphous calcium carbonate is an important precursor for biomineralization in marine organisms. Key outstanding problems include understanding the structure of amorphous calcium carbonate and rationalizing its metastability as an amorphous phase. Here we report high-quality atomistic models of amorphous calcium carbonate generated using state-of-the-art interatomic potentials to help guide fits to X-ray total scattering data. Exploiting a recently developed inversion approach, we extract from these models the effective Ca∙∙∙Ca interaction potential governing the structure. This potential contains minima at two competing distances, corresponding to the two different ways that carbonate ions bridge Ca 2+ -ion pairs. We reveal an unexpected mapping to the Lennard-Jones–Gauss model normally studied in the context of computational soft matter. The empirical model parameters for amorphous calcium carbonate take values known to promote structural complexity. We thus show that both the complex structure and its resilience to crystallization are actually encoded in the geometrically frustrated effective interactions between Ca 2+ ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The quest for safer nuclear fuels

Here, to increase the safety and performance of nuclear reactors, researchers are working on developing ‘accident tolerant’ fuels. Ideally, these fuel pellets will be less susceptible to swelling and be more resistant to extreme heat. This led chemists to develop uranium dioxide (UO 2 ) doped with divalent elements, E 2+ . These dopants can introduce vacancies into the lattice that lead to larger grain sizes, which create longer pathways for fission gasses to migrate out of the grains. In undoped fuel, fission gas build-up at grain boundaries causes the fuel pellet to swell. Therefore, having larger grains that slow fission gas migration into grain boundaries is highly desired. Divalent Cr 2+ ions are a commonly explored dopant for UO 2 fuels. These materials have larger grain sizes while maintaining the dissolution kinetics of undoped UO 2 , making their long-term storage and disposal paths similarly safe. Atomistic models have suggested that Mn 2+ is a conceivable alternative doping material, with the potential to increase the sizes of grains even more than Cr 2+ dopants. As such, experimental studies are required to understand the physical chemistry and reactivity of these potential fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Origins of the change in mechanical strength of silicon/gold nanocomposites during irradiation

Abstract Silicon-based layered nanocomposites, comprised of covalent-metal interfaces, have demonstrated elevated resistance to radiation. The amorphization of the crystalline silicon sublayer during irradiation and/or heating can provide an additional mechanism for accommodating irradiation-induced defects. In this study, we investigated the mechanical strength of irradiated Si-based nanocomposites using atomistic modeling. We first examined dose effects on the defect evolution mechanisms near silicon-gold crystalline and amorphous interfaces. Our simulations reveal the growth of an emergent amorphous interfacial layer with increasing dose, a dominant factor mitigating radiation damage. We then examined the effect of radiation on the mechanical strength of silicon-gold multilayers by constructing yield surfaces. These results demonstrate a rapid onset strength loss with dose. Nearly identical behavior is observed in bulk gold, a phenomenon that can be rooted to the formation of radiation-induced stacking fault tetrahedra which dominate the dislocation emission mechanism during mechanical loading. Taken together, these results advance our understanding of the interaction between radiation-induced point defects and metal-covalent interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structural water in amorphous carbonate minerals: ab initio molecular dynamics simulations of X-ray pair distribution experiments

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO 3 ·nH 2 O systems, where M 2+ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Finally, water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Size-dependent shape distributions of platinum nanoparticles

Transmission electron microscopy revealed size-dependent shape distributions in platinum nanoparticles, which were consistent with trends observed by applying Boltzmann statistics to the energy computed with atomistic models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tuning martensitic transformations via coherent second phases in nanolaminates using free energy landscape engineering

In this work we explore the possibilities and limitations of using a coherent second phase to engineer the thermo-mechanical properties of a martensitic alloy by modifying the underlying free energy landscape that controls the transformation. We use molecular dynamics simulations of a model atomistic system where the properties of a coherent, nanoscale second phase can be varied systematically. With a base martensitic material that undergoes a temperature-induced transformation from a cubic austenite to a monoclinic martensite, simulations show significant ability to engineer the transformation temperatures, from a ~50% reduction to a ~200% increase, with 50 at.% of the cubic second phase. We establish correlations between the properties of the second phase, the transformation characteristics, and the microstructure via the free energy landscape of the twophase systems. Coherency stresses have a strong influence on the martensitic variants observed and can even cause the non-martensitic second phase to undergo a transformation. Reducing the stiffness of the second phase increases the transformation strain and modifies the martensitic microstructure, increasing the volume fraction of the transformed material. This increase in transformation strain is accompanied by a significant increase in Af and thermal hysteresis, while Ms remains unaltered. Our findings on the tunability of martensitic transformations can be used for informed searches of second phases to achieve desired material properties, such as achieving room temperature, lightweight shape memory alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulations of nonradiative processes in semiconductor nanocrystals

The description of carrier dynamics in spatially confined semiconductor nanocrystals (NCs), which have enhanced electron–hole and exciton–phonon interactions, is a great challenge for modern computational science. These NCs typically contain thousands of atoms and tens of thousands of valence electrons with discrete spectra at low excitation energies, similar to atoms and molecules, that converge to the continuum bulk limit at higher energies. Computational methods developed for molecules are limited to very small nanoclusters, and methods for bulk systems with periodic boundary conditions are not suitable due to the lack of translational symmetry in NCs. This perspective focuses on our recent efforts in developing a unified atomistic model based on the semiempirical pseudopotential approach, which is parameterized by first-principle calculations and validated against experimental measurements, to describe two of the main nonradiative relaxation processes of quantum confined excitons: exciton cooling and Auger recombination. Here, we focus on the description of both electron–hole and exciton–phonon interactions in our approach and discuss the role of size, shape, and interfacing on the electronic properties and dynamics for II–VI and III–V semiconductor NCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Chemical order transitions within extended interfacial segregation zones in NbMoTaW

Interfacial segregation and chemical short-range ordering influence the behavior of grain boundaries in complex concentrated alloys. In this study, we use atomistic modeling of a NbMoTaW refractory complex concentrated alloy to provide insight into the interplay between these two phenomena. Hybrid Monte Carlo and molecular dynamics simulations are performed on columnar grain models to identify equilibrium grain boundary structures. Here our results reveal extended near-boundary segregation zones that are much larger than traditional segregation regions, which also exhibit chemical patterning that bridges the interfacial and grain interior regions. Furthermore, structural transitions pertaining to an A2-to-B2 transformation are observed within these extended segregation zones. Both grain size and temperature are found to significantly alter the widths of these regions. An analysis of chemical short-range order indicates that not all pairwise elemental interactions are affected by the presence of a grain boundary equally, as only a subset of elemental clustering types are more likely to reside near certain boundaries. The results emphasize the increased chemical complexity that is associated with near-boundary segregation zones and demonstrate the unique nature of interfacial segregation in complex concentrated alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Inherent-state melting and the onset of glassy dynamics in two-dimensional supercooled liquids

Below the onset temperature T o , the equilibrium relaxation time of most glass-forming liquids exhibits glassy dynamics characterized by a super-Arrhenius temperature dependence. In this supercooled regime, the relaxation dynamics also proceeds through localized elastic excitations corresponding to hopping events between inherent states, i.e., potential-energy-minimizing configurations of the liquid. Despite its importance in distinguishing the supercooled regime from the high-temperature regime, the microscopic origin of T o is not yet known. Here, we construct a theory for the onset temperature in two dimensions and find that an inherent-state melting transition, described by the binding–unbinding transition of dipolar elastic excitations, delineates the supercooled regime from the high-temperature regime. The corresponding melting transition temperature is in good agreement with the onset temperature found in various two-dimensional (2D) atomistic models of glass formers and an experimental binary colloidal system confined to a water–air interface. Additionally, we find the predictions for the renormalized elastic moduli to agree with the experimentally observed values for the latter 2D colloidal system. We further discuss the predictions of our theory on the displacement and density correlations at supercooled conditions, which are consistent with observations of the Mermin–Wagner fluctuations in experiments and molecular simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Tuning the surface reactivity of oxides by peroxide species

The Mars–van Krevelen mechanism is the foundation for oxide-catalyzed oxidation reactions and relies on spatiotemporally separated redox steps. Herein, we demonstrate the tunability of this separation with peroxide species formed by excessively adsorbed oxygen, thereby modifying the catalytic activity and selectivity of the oxide. Using CuO as an example, we show that a surface layer of peroxide species acts as a promotor to significantly enhance CuO reducibility in favor of H 2 oxidation but conversely as an inhibitor to suppress CuO reduction against CO oxidation. Together with atomistic modeling, we identify that this opposite effect of the peroxide on the two oxidation reactions stems from its modification on coordinately unsaturated sites of the oxide surface. By differentiating the chemical functionality between lattice oxygen and peroxide, these results are closely relevant to a wide range of catalytic oxidation reactions using excessively adsorbed oxygen to activate lattice oxygen and tune the activity and selectivity of redox sites.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗
DOE OSTI2021

Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.

58 GEOSCIENCES↗
DOE OSTI2026

Momentum-Resolved Direct Observation of Chiral Phonons in Elemental Tellurium

Chiral phonons carry finite phonon angular momentum, yet their momentum-resolved behavior in acoustic branches remains largely unexplored. In this work, we report the direct detection of chiral acoustic phonons in elemental Te. Through a combination of high-energy-resolution inelastic x-ray scattering experiments and atomistic modeling, we confirm that longitudinal acoustic modes develop pronounced circular atomic motion by mixing with transverse modes. Our results establish a direct link between unexpected scattering intensity in nominally forbidden geometries and circular phonons, offering a framework for characterizing momentum-resolved phonon chirality and highlighting the emergence of intrinsic phonon angular momentum in acoustic branches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Effect of surface steps on chemical ordering in the subsurface of Cu(Au) solid solutions

Atomic steps are typically assumed to only influence surface phenomena of solids. In contrast, in this paper we show, using in-situ atomic-resolution electron microscopy observations and atomistic modeling, the pronounced effect from surface stepsin inducing compositional and structural evolution in the subsurface of a Cu(Au) solid solution. We find that Au surface segregation results in a stacked sequence of Cu-Au ordered phases in the subsurface. The presence of a monatomic step at the surface induces the formation of an anti-phase boundary that extends from the surface step to deeper layers by maintaining the same composition profile associated with each terrace. The bunching of surface steps induces chemical disordering of the Cu-Au ordered phases in the subsurface region of the bunched steps. These results demonstrate the instant propagation of surface dynamics of atomic steps into the subsurface region and can find broader applicability in utilizing surface defects to tune the composition and structure in the subsurface of alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pressure-induced polyamorphic transition in CaA l 2 O 4 glass

In situ high-pressure ultrasonic velocity measurements of CaAl 2 O 4 glass reveal abrupt irreversible discontinuities in the elastic wave velocities at ~8–10 GPa. Total structure factor and pair distribution functions measured by synchrotron x-ray diffraction show a rapid change in the intermediate range structure attributed to a rearrangement of calcium ions over this narrow pressure region. Atomistic models obtained from molecular dynamics simulations reveal that this intermediate range structure is explained by a transition of Ca–O void radius distribution from a bimodal distribution with peaks at ~2.1 and ~2.4 Å to a single distribution centered at ~2.1 Å. Further, these abrupt structural changes involving the rapid increase in elastic wave velocity are markedly different to the continuous transformations observed in conventional network-forming glasses, such as SiO 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗