Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic simulation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Machine Learning Reveals Memory of the Parent Phases in Ferroelectric Relaxors Ba(Ti$_{1-x}$,Zr x )O 3

Machine learning has been establishing its potential in multiple areas of condensed matter physics and materials science. Here, in this work, an unsupervised machine learning workflow is developed and used within a framework of first-principles-based atomistic simulations to investigate phases, phase transitions, and their structural origins in ferroelectric relaxors, Ba(Ti 1-x ,Zr x )O 3 . The applicability of the workflow is first demonstrated to identify phases and phase transitions in the parent compound, a prototypical ferroelectric BaTiO 3 . Then the workflow is applied for Ba(Ti 1-x ,Zrx)O 3 with x ≤ 0.25 to reveal i) that some of the compounds bear a subtle memory of BaTiO 3 phases beyond the point of the pinched phase transition, which could contribute to their enhanced electromechanical response; ii) the existence of peculiar phases with delocalized precursors of nanodomains—likely candidates for the controversial polar nanoregions; and iii) nanodomain phases for the largest concentrations of x.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Direct Observation of Elusive (DTBM‐SEGPHOS)CuH Monomer Enables Mechanistic Insights Into Hydrocupration, Aggregation, and Dynamics of Alkene Functionalization Catalysis

The bulky diphosphine DTBM-SEGPHOS is widely employed in CuH-catalyzed transformations as it provides remarkably active catalyst systems. The transient (DTBM-SEGPHOS)CuH monomer (LCuH) is the often-invoked active species. However, its instability has prevented spectroscopic characterization and mechanistic elucidation, hindering mechanistic understanding. We report low-temperature NMR spectroscopic characterization of LCuH, enabling quantitative kinetic analysis of the stoichiometric hydrocupration and catalytic hydroboration of cyclopentene, as well as the structural identification of two CuH clusters. LCuH inserts cyclopentene at −43°C, reaffirming its high reactivity toward olefins. LCuH deactivates to form L 2 Cu 3 H 3 and L 2 Cu 4 H 4 clusters, in which LCuH dimerization initiates aggregation. Kinetic analysis of reactions of unactivated alkenes indicates that competing on-cycle alkene hydrocupration and LCuH dimerization impact performance, as catalyst deactivation and turnover occur on comparable timescales. Structure–activity analysis using atomistic simulations shows that the steric profile of DTBM-SEGPHOS increases the CuH dimerization barrier by ∼7.7 kcal mol−1 compared to that of SEGPHOS, rationalizing the unique ability of DTBM-SEGPHOS to stabilize a reactive monomer for hydrocupration of broader alkene substrates. These findings illustrate the fundamental design principle that steric control of aggregation governs CuH catalyst performance, explaining both the exceptional activity of (DTBM-SEGPHOS)CuH and the limitations imposed by competing deactivation.

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

Active Learning A Neural Network Model For Gold Clusters & Bulk From Sparse First Principles Training Data

Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Given the high computational cost of first-principles calculations, molecular modeling and atomistic simulations such as molecular dynamics (MD) has proven to be an important complementary tool to aid this understanding. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities (similar to 10 4 or greater) of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. A highly desirable goal is minimize the number of training data, especially if the underlying reference model is first-principles based and hence expensive. In this work, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (similar to 1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with similar to 500 total reference calculations. Using an extensive DFT test set of similar to 1100 configurations, we show that our AL-NN is able to accurately predict both the DFT energies and the forces for clusters of a myriad of different sizes. Our NN predictions are within 30 meV/atom and 40 meV/angstrom of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments. We finally show that our AL-NN model also captures bulk properties reasonably well, even though they were not included in the training data.

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

Unlocking High Capacity and Reversible Alkaline Iron Redox Using Silicate-Sodium Hydroxide Hybrid Electrolytes

Alkaline iron (Fe) batteries are attractive due to the high abundance, low cost, and multiple valent states of Fe but show limited columbic efficiency and storage capacity when forming electrochemically inert Fe 3 O 4 on discharging and parasitic H 2 on charging. Herein, sodium silicate is found to promote Fe(OH) 2 /FeOOH against Fe(OH) 2 /Fe 3 O 4 conversions. Electrochemical experiments, operando X-ray characterization, and atomistic simulations reveal that improved Fe(OH) 2 /FeOOH conversion originates from (i) strong interaction between sodium silicate and iron oxide and (ii) silicate-induced strengthening of hydrogen-bond networks in electrolytes that inhibits water transport. Furthermore, the silicate additive suppresses hydrogen evolution by impairing energetics of water dissociation and hydroxyl de-sorption on iron surfaces. In conclusion, this new silicate-assisted redox chemistry mitigates H 2 and Fe 3 O 4 formation, improving storage capacity (199 mAh g -1 in half-cells) and coulombic efficiency (94 % after 400 full-cell cycles), paving a path to realizing green battery systems built from earth-abundant materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine learning for accuracy in density functional approximations

Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. In this paper, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.

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

Tailoring the Angular Mismatch in MoS 2 Homobilayers through Deformation Fields

Ultrathin MoS 2 has shown remarkable characteristics at the atomic scale with an immutable disorder to weak external stimuli. Ion beam modification unlocks the potential to selectively tune the size, concentration, and morphology of defects produced at the site of impact in 2D materials. Combining experiments, first-principles calculations, atomistic simulations, and transfer learning, it is shown that irradiation-induced defects can induce a rotation-dependent moiré pattern in vertically stacked homobilayers of MoS 2 by deforming the atomically thin material and exciting surface acoustic waves (SAWs). Furthermore, the direct correlation between stress and lattice disorder by probing the intrinsic defects and atomic environments are demonstrated. The method introduced in this paper sheds light on how engineering defects in the lattice can be used to tailor the angular mismatch in van der Waals (vdW) solids.

2D materials↗
DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗
DOE OSTI2021

A Review: Microstructural and Phase Evolution in Alloys during Extended Plastic Deformation

Materials in service and during processing are often subjected to plastic deformation. For multi-phase metallic alloys, simple geometric models and atomistic simulations suggest that two distinctive regimes in these materials’ evolution during deformation exist. At low strains, evolutions are often dominated by the kinetic roughening of interfaces, which results from the superdiffusive transport of matter in sheared crystals. Furthermore, at high strains and temperatures where thermal diffusion is sluggish, on the other hand, shearing- induced forced atomic mixing dominates these evolutions, resulting in significantly enhanced solubility. Distinguishing these two regimes is shown to provide a convenient framework for rationalizing and analyzing recent experiments and simulations on wear of layered structures in the low strain-regime and nonequilibrium phase co-existence and self-organization in highly immiscible or reactive alloy systems in the high-strain regime.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Enhanced Radiation Damage Tolerance of Amorphous Interphase and Grain Boundary Complexions in Cu-Ta

Amorphous interfacial complexions are particularly resistant to radiation damage and have been primarily studied in alloys with good glass-forming ability, yet recent reports suggest that these features can form even in immiscible alloys such as Cu-Ta under irradiation. In this work, the mechanisms of damage production and annihilation due to primary knock-on atom collisions are investigated for amorphous interphase and grain boundaries in a Cu-Ta alloy using atomistic simulations. Amorphous complexions, in particular amorphous interphase complexions that separate Cu and Ta grains, result in less residual defect damage than their ordered counterparts. Stemming from the nanophase chemical separation in this alloy, the amorphous complexions exhibit a highly heterogeneous distribution of atomic excess volume, as compared to a good glass former like Cu-Zr. Complexion thickness, a tunable structural descriptor, plays a vital role in damage resistance. Thicker interfacial films are more damage-tolerant because they alter the defect production rate due to differences in intrinsic displacement threshold energies during the collision cascade. Overall, the findings of this work highlight the importance of interfacial engineering in enhancing the properties of materials operating in radiation-prone environments and the promise of amorphous complexions as particularly radiation damage-tolerant microstructural features.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Carbon nanotube (CNT) metal composites exhibit greatly reduced radiation damage

Radiation damage of structural materials leads to mechanical property degradation, eventually inducing failure. Secondary-phase dispersoids or other radiation defect sinks are often added to materials to boost their radiation resistance. We demonstrate that a metal composite made by adding 1D carbon nanotubes (CNTs) to aluminum (Al) exhibits superior radiation resistance. In situ ion irradiation with transmission electron microscopy (TEM) and atomistic simulations together reveal the mechanisms of rapid defect migration to CNTs, facilitating defect recombination and enhancing radiation tolerance. The origin of this effect is an evolving stress gradient in the Al matrix resulting from CNT transformation under irradiation, and the stability of resulting carbides. Extreme value statistics of large defect behavior in our simulations highlight the role of CNTs in reducing accumulated damage. Furthermore, this approach to controlling defect migration represents a promising opportunity to enhance the radiation resistance of nuclear materials without detrimental effects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Three-dimensional atomic scale characterization of {11$\overline{2}$2} twin boundaries in titanium

The {11$\overline{2}$2}<11$\overline{23}$> compression twin can accommodate a considerable amount of strain under c-axis compression in Ti. However, unlike the tensile twin, the structure of the {10$\overline{1}$2}$\langle$$\overline{1}$011$\rangle$ compressive twins has not been completely characterized. In this study, we apply a combined technique of HR-TEM characterization, topological analysis, and atomistic simulations to explore the facets that bound the {11$\overline{2}$2}<11$\overline{23}$> twin in Ti. In addition to the currently known facets (CTB and B-Py), six new facets are observed and categorized for the first time from atomic-scale TEM observations along five crystallographic directions. The six new facets are (11$\overline{2}$0)//(11$\overline{2}$6), PrPr1, PyPy1, (2$\overline{11}$1)//($\overline{1}$2$\overline{1}$2), (1$\overline{1}$04)//(01$\overline{11}$), and (01$\overline{1}$0)//(2$\overline{11}$4). Results from the topological and computational analysis are in reasonable agreement with and support the HRTEM observations. Specifically, (1) the observed facets align with low-index interfaces in both twin and matrix domains, (2) the facets with lower surface energies are found to form extended interfaces, and (3) high-surface-energy facets are found at the twin tip region and explained by the fact that the energy of the combined facet and facet junction configuration is energetically preferred in the twin tip region. These results not only provide a comprehensive understanding of the 3D structure of the {11$\overline{2}$2}<11$\overline{23}$> compressive twins in Ti, but also validate the MD procedure and the Ti interatomic potential employed. This is extremely important for future study of the {11$\overline{2}$2} twin mobility and interactions with other defects, both features that remain extremely challenging to capture in experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Grain boundary softening from stress assisted helium cavity coalescence in ultrafine-grained tungsten

The formation of helium cavities in coarse-grained materials produces hardening proportional to the number density and size of the cavities and due to the interaction of dislocations with intragranular helium defects. In nanostructured metals containing a high density of interfacial sinks, preferential cavity formation in the grain boundaries instead produces softening that is often attributed to enhanced interfacial plasticity. Here, employing two grades of ultrafine-grained tungsten, we explore this effect using targeted implantation studies to map cavity evolution as a function of the irradiation conditions and quantify its impact on the mechanical response through nanoindentation. Softening is reported at implantation temperatures above the threshold for preferential grain boundary cavity formation but at a sufficiently low fluence prior to the growth of intragranular cavities. Collective changes in the mean cavity size, density, and morphology beneath a residual impression on an implanted surface indicate that cavity coalescence accompanied the reduction in hardness. Complementary atomistic simulations demonstrate that, in tungsten grain structures exhibiting softening, grain boundary bubble coalescence is driven by stress concentrations that further act to localize strain in the grain boundaries through cooperative deformation processes involving local atomic shuffling and sliding, dislocation emission, and even the nucleation of unstable twinning events.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Effect of Grain Boundary Facet Junctions on Segregation and Embrittlement

Junctions are discontinuities in flat grain boundaries that arise in all polycrystalline materials and are thought to play important roles in the response of a grain boundary network to thermal and mechanical loads. A key open question concerns the mechanisms by which solute segregation to junctions impacts properties of the grain boundary. Here, in this work, we investigate the influence of grain boundary facet junctions on solute embrittlement, and we present an analytical model that uses the hydrostatic stress field contributed by dislocations at multiple junctions to describe these effects. Specifically, we study junctions between {112} facets of various lengths in Au $\langle111\rangle$ Σ3 tilt grain boundaries. Copper and silver solutes are employed to determine if the effect of junctions on solute segregation and embrittlement is dependent on size relative to the host. Combined, atomistic simulation data and the analytical model show that Cu and Ag have opposite segregation responses to junctions due to the sign of the hydrostatic stress field induced by junctions. However, a positive shift in the embrittling potency is computed near junctions regardless of solute type or the stress state of the segregation site. Hence, for the conditions studied, junctions consistently shift the energetic landscape towards embrittlement.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ultimate compressive strength and severe plastic deformation of equilibrated single-crystalline copper nanoparticles

Mechanical properties and deformation mechanisms of defect-free copper nanoparticles are investigated by combining experiments with atomistic simulations. The compressive strength of the particles increases with decreasing size and tends to saturate near the theoretical strength in the small-size limit. In this limit, the intrinsic size dependence of the strength is governed by the stochastic nature of dislocation nucleation near the particle surface. The particle deformation process evolves from the initial strain softening to strain hardening as the particle accumulates residual damage. The normalized strength-size relation for Cu is compared with those for Au, Ni, and Pt. The lack of universal behavior among the four FCC metals is discussed. Heavily deformed Cu nanoparticles develop polycrystalline structures and change the lattice orientation from [111] to [110]. The experiments and simulations reveal the twinning mechanism of the lattice rotation leading to the new grain formation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Formation of carbon homonuclear bonds in β-SiC under neutron irradiation at various temperatures and neutron doses

To elucidate radiation defect processes in SiC, Raman spectroscopy was systematically applied to high-purity, polycrystalline β-SiC that was neutron irradiated at a range of temperature and dose conditions. The analysis specifically focused on formation of carbon homonuclear bonds by irradiation; these bonds were indicated by D and G bands and amorphous carbon peaks. Intensity of the carbon peaks relative to SiC peaks significantly decreased in the case of high temperature and/or high neutron dose of 500 °C to 29 displacements per atom (dpa) and about 800 °C to 1.38 and 29 dpa. The absence of carbon bond peaks under those conditions was explained by growth of stoichiometric defect clusters, consistent with previous atomistic simulations on SiC defect stability. The lack of Raman bands associated with carbon clusters under high-temperature and high-dose radiation conditions accounts for the resistance of SiC to phase separation under irradiation. The findings further suggest that material compositions and chemical properties that are inherently resistant to chemical disordering under high-dose radiation conditions are indicative of the long-term durability of ceramic compounds in radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Porosity in nuclear graphite and its impact on nuclear reactor science and criticality safety applications

Porosity in nuclear-grade graphite significantly influences its low-energy neutron scattering, yet its effect on underlying phonon properties remains debated. Here, this work integrates inelastic and small-angle neutron scattering (INS/SANS) experiments, advanced atomistic simulations with a novel machine-learned potential (DeepMD), total cross-section measurements, and neutronics calculations (SCALE, MCNP, OpenMC) to investigate porosity’s impact on neutron thermalization. INS measurements on diverse graphite grades reveal no discernible porosity effect on phonon spectra, which align with crystalline graphite. Conversely, total cross-section data below ≈10 meV show increased scattering attributable to SANS. Our DeepMD simulations demonstrate that realistic micropores do not distort phonon spectra, challenging the assumptions in current ENDF/B-VIII.1 porosity thermal scattering laws (TSLs). These TSLs, based on random atom removal, produce unphysical phonon spectra and inflate inelastic cross-sections. Augmenting a crystalline TSL with an SANS component accurately captures experimental total cross-sections. Neutronics benchmarks (ICSBEP/IRPhE) show ENDF porosity TSLs unphysically increase neutron multiplication factor, keff. Crucially, incorporating SANS physics (NCrystal/OpenMC) indicates accurately modeled porosity negligibly affects keff, reactor physics, or criticality safety.

Critical benchmarks↗
DOE OSTI2020

Polymer-cement composites with adhesion and re-adhesion (healing) to casing capability for geothermal wellbore applications

Deterioration of cement/casing adhesion in wellbore scenarios can result in unwanted and potentially harmful leakage with the potential of serious repair costs. In this work, we explore the use of self-healing polymers added to conventional wellbore cements as a way to bring about self-healing and readhering (to steel casing) properties to the composite material. The polymers are pH resistant and seem to improve the cement integrity after exposure to typical chemical and thermal stresses encountered under geothermal wellbore conditions. We find that addition of about 10-15 wt% of polymer to the cement visually increases its resistance to fracturing from exposure to geothermal conditions, while the adhesive strength of cement/stainless steel increases with curing time for a period of about 10 days. Self-healing capability was demonstrated by permeability analysis showing that polymer-cement composites reduce flow by 50-70% at cement bulk and at the cement/steel interface. Use of atomistic simulations imply that these polymers have good wetting properties on the steel surfaces. Analysis of the interactions between steel/polymer and cement/polymer show that they are complementary, resulting in a wider range of bonding patterns. Cracks are likely to expose under-coordinated sites that result in more bonding interactions, which agrees well with the permeability measurements showing high degree of healed cracks and healed (cement-steel) interfacial gaps together with an overall increased in structural integrity of these advanced polymer-cement composite materials.

Rod, Kenton A.↗
DOE OSTI2021

New insights into creep characteristics of calcium silicate hydrates at molecular level

The fundamental mechanisms under concrete creep are far from being fully understood, especially at the molecular level. Hereby, a calcium-silicate-hydrate (C-S-H) molecular model is developed to explain, for the first time, the creep characteristics at various stress states, temperature levels and water contents, which are not accessible experimentally. Rather tensile and compressive loadings, C-S-H only creeps under shear loadings originating from the sliding of the calcium silicate layers over each other as the interlayer component (water and ions) acts as a lubricator. A heterogeneous creep characteristic is observed. Elevated temperature reduces the interlayer lubricator viscosity and weakens the interfacial adhesion between the layers and the interlayer lubricator, which accelerates C-S-H creep. The removal of interlayer water enhances the creep resistance, resulting from the reduced interlayer space and enhanced interfacial adhesion. The atomic-level mechanisms explain the inter-CSH-particle behaviours at the microscale, which bridges the gap between atomistic simulation and microcosmic phenomenon.

36 MATERIALS SCIENCE↗