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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Force field for calcium sulfate minerals to predict structural, hydration, and interfacial properties

Calcium sulfates such as anhydrite, hemihydrate, and gypsum find widespread use in building materials, implants, and tissue healing. We introduce a simple and compatible atomistic force field for all calcium sulfate phases that reproduces a wide range of experimental data including lattice parameters, surface, hydration, mechanical, and thermal properties in 1% to 5% accuracy relative to experiments. The performance is several times better than prior force fields and DFT methods, which lead to errors in structures and energies up to 100%. We explain (hkl) cleavage energies, the dynamics of (hkl) water interfaces, and new insights into molecular origins of crystal-facet specific hydration and solubility. Impressive agreement of computed and experimentally measured hydration energies is shown. The models add to the Interface force field (IFF) and are compatible with multiple force fields (CHARMM, AMBER, GROMOS, CVFF, PCFF, OPLS-AA) for property predictions of sulfate-containing materials from atoms to the large nanometer scale.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

A General Model for the Electrochemical Double Layer in Solids

It has long been clear that electrochemical double layer / space charge zone in the vicinity of interfaces and other extended defects greatly influences transport and reactivity. Over the past several years, multiple microscopic measurements of defect concentrations in the vicinity of solid-solid interfaces made clear that dilute-case thermodynamic theories to describe the behavior of charged species in space charge zones are untenable. Inhomogeneous thermodynamics – i.e., phase field models – have delivered the most success in describing both microscopic and macroscopic experiments. We still need a generalized approach that quantitatively reproduces the microscopic and macroscopic evidence, while remaining accessible to a wide range of researchers and engineers. In my talk I will present such a framework, which builds on the previously introduced Poisson-Cahn theory but utilizes data science methods which come into the thermodynamic framework in a well-defined way. I will review recent applications to both microscopic and macroscopic datasets while previewing our plans for connecting the theory to atomistic / first principles calculations and experiment simultaneously. I will also review plans for releasing open source code that runs on desktop-level computational resources.

David Mebane↗
DOE OSTI2023

Mending cracks atom-by-atom in rutile TiO2 with electron beam radiolysis

Abstract Rich electron-matter interactions fundamentally enable electron probe studies of materials such as scanning transmission electron microscopy (STEM). Inelastic interactions often result in structural modifications of the material, ultimately limiting the quality of electron probe measurements. However, atomistic mechanisms of inelastic-scattering-driven transformations are difficult to characterize. Here, we report direct visualization of radiolysis-driven restructuring of rutile TiO 2 under electron beam irradiation. Using annular dark field imaging and electron energy-loss spectroscopy signals, STEM probes revealed the progressive filling of atomically sharp nanometer-wide cracks with striking atomic resolution detail. STEM probes of varying beam energy and precisely controlled electron dose were found to constructively restructure rutile TiO 2 according to a quantified radiolytic mechanism. Based on direct experimental observation, a “two-step rolling” model of mobile octahedral building blocks enabling radiolysis-driven atomic migration is introduced. Such controlled electron beam-induced radiolytic restructuring can be used to engineer novel nanostructures atom-by-atom.

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

Computing Thermodynamic Properties of Fluids Augmented by Nanoconfinement: Application to Pressurized Methane

Nanoconfined fluids exhibit remarkably different thermodynamic behavior compared to the bulk phase. These confinement effects render predictions of thermodynamic quantities of nanoconfined fluids challenging. In particular, confinement creates a spatially varying density profile near the wall that is primarily responsible for adsorption and capillary condensation behavior. Significant fluctuations in thermodynamic quantities, inherent in such nanoscale systems, coupled to strong fluid–wall interactions give rise to this near-wall density profile. Empirical models have been proposed to explain and model these effects, yet no first-principles based formulation has been developed. We present a statistical mechanics framework that embeds such a coupling to describe the effect of the fluid–wall interaction in amplifying the near-wall density behavior for compressible gases at elevated pressures such as pressurized methane in confinement. We show that the proposed theory predicts accurately the adsorbed layer thickness as obtained with small-angle neutron scattering measurements. Furthermore, the predictions of density under confinement from the proposed theory are shown to be in excellent agreement with available experimental and atomistic simulations data for a range of temperatures for nanoconfined methane. While the framework is presented for evaluating the near-wall density, owing to its rigorous foundation in statistical mechanics, the proposed theory can also be generalized for predicting phase-transition and nonequilibrium transport of nanoconfined fluids.

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

Assessment of thermally driven local structural phase changes in 1⁢𝑇′−MoTe 2

The role of layer disorder is important in establishing the topological phases of MoTe 2 . A rich tapestry of atomic ordering influences the structural phase transitions (SPTs), but there is little understanding of the mechanistic details of the phase transition. An atomistic level study was conducted to investigate the local structure of the 1⁢𝑇′ and 𝑇 𝑑 phases of MoTe 2 by using the pair distribution function (PDF) technique. While the average structure exhibits an SPT and coexistence of phases as a function of temperature, the local structure showed the suppression of SPT. The sample retained its monoclinic structure at all temperatures in short-range order. A sharp PDF peak observed at short distances indicated a strong atom-atom correlation between the Mo and Te atoms within the Mo octahedra. In addition, a large-box modeling of the PDF data indicated a preferential motion of Te atoms towards 𝑐 axis at all temperatures. Structural defects, such as stacking faults, likely result in the coexistence of phases in the average structure and suppress the local SPT of MoTe 2 . These results are stepping stones to understand the long-debated origins of structural, vibrational, and electronic properties of MoTe 2 and similar transition metal dichalcogenides.

2-dimensional systems↗
DOE OSTI2021

Phase-field modeling of the interactions between an edge dislocation and an array of obstacles

Obstacles, such as voids and precipitates, are prevalent in crystalline materials. They strengthen crystals by serving as barriers to dislocation glide. Here in this work, we develop a phase-field dislocation dynamics (PFDD) technique for investigating the interactions between dislocations and second-phase obstacles, which can be either voids or precipitates. The PFDD technique is constructed to account for elastic heterogeneity, elastic anisotropy, dissociation of the dislocation, and dislocation transmission across bicrystalline interfaces. Within the framework, we present a model for “pseudo-voids”, which are voids shearable by dislocations, in contrast to unphysical, unshearable voids in conventional phase-field dislocation formulations. We employ the PFDD technique to investigate the in-plane interactions between an edge dislocation and an array of nano-scale obstacles with different spacings. In this application, the interactions take place in glide planes of either a face-centered cubic (FCC) Cu or a body-centered cubic (BCC) Nb matrix, while the precipitates have a Cu 1-x Nb x composition, with x varying from 0.1 to 0.9. Our atomistic simulations find that the alloy precipitates can have an FCC, an amorphous, or a BCC phase, depending on the compositional ratio between Cu and Nb, i.e., value of x. Among all types of obstacles, the critical stresses for dislocation bypass are the highest for unshearable amorphous precipitates, followed by shearable crystalline precipitates, and then the pseudo-voids.

42 ENGINEERING↗
DOE OSTI2025

Diffusion and Deformation Mechanism Maps in Li Metal from Atomistic Simulations

The rate of Li transport in the Li metal anode is important for the operation of Li metal-solid state batteries. Here, transport rates due to diffusion and creep are predicted in Li using atomistic simulations. First, molecular dynamics is used to estimate the rate of Li diffusion along dislocations and in grain boundary triple junctions. By combining this data with that from a prior study of grain boundary diffusion the dominant mechanisms and rates of self-diffusion in Li polycrystals are predicted as a function of grain size, grain shape, dislocation density, and temperature. Second, the dominant creep mechanisms are predicted and used to estimate critical current densities and void annihilation times. Grain boundary sliding and coble creep are the dominant mechanisms for micron-sized grains. Lastly, a continuum model for interfacial contact loss reveals that high dislocation densities of ∼10 12 /cm 2 enable achieving battery performance targets for Li grain sizes of ∼10 μm.

Batteries↗
DOE OSTI2021

Dynamical origins of weakly coupled relaxor behavior in Sn-doped (Ba,Ca)TiO 3 -BiScO 3

The peculiar characteristics of relaxors, viz., a frequency-dependent dielectric permittivity peak and good functional properties (dielectric, electromechanical, electrocaloric, etc.), are attributed to nanoscale regions with correlated dipoles, or polar nanoregions (PNRs). However, the exact nature of PNRs and their contribution to relaxor behavior remains debatable. In recent years, solid solutions of BaTiO3-BiMeO3 (where Me is a metal), have emerged as an interesting system with characteristics in between that of relaxors and dipole glasses. Here, we have examined the atomistic origins of weakly coupled relaxor behavior, specifically with regard to formation of PNRs, in Sn-doped (1-x)(Ba,Ca)TiO 3-x BiScO 3 using macroscopic polarization and neutron dynamic pair distribution function measurements. We show that the short-range atomic correlations observed within the PNRs dynamically fluctuate with frequencies of the order of THz. Furthermore the composition-dependent dielectric and polarization behaviors are critically influenced by the relative stability of the atomic correlations near ~1 THz, while the instantaneous atomic correlations are largely independent of x. The current results are discussed based on a model of intrinsic local modes distributed in a dielectrically soft matrix.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Perspective—Lithium Metal Nucleation and Growth on Conductive Substrates: A Multi-Scale Understanding from Atomistic, Nano-, Meso-, to Micro-Scales

Reversible lithium (Li) deposition and stripping on conductive substrates like copper is vital for long-life high-energy rechargeable Li batteries. The reversibility is believed hindered by solid electrolyte interphase (SEI) formation and Li dendritic growth. Via in situ, operando cryo-microscopy and molecular dynamics simulation, we discovered amorphous Li before a disorder-order phase transition to bcc phase. The kinetics plays a significant role in Li nucleation and morphology evolution. This perspective on Li nucleation and growth from atomistic to nano- (<20 nm), meso- (20–100 nm), and micro-scales (>100 nm) provides a practical guidance on regulating dense Li deposits reversibly for long life performance.

25 ENERGY STORAGE↗
DOE OSTI2022

Temperature Effects in Flexible Adsorption Processes for Amorphous Microporous Polymers

A collection of atomistic molecular simulations is reported that illustrate the impact of adsorption temperature on species uptake and adsorbate-induced structural rearrangement for amorphous polymers of intrinsic microporosity. Temperature-sensitive structural rearrangement is evaluated by contrasting two methods: standard grand canonical Monte Carlo simulations using a rigid framework approximation and a combined Monte Carlo/molecular dynamics approach that fully incorporates framework flexibility. We report single-component gas phase adsorption isotherms for CH 4 , C 2 H 4 , C 2 H 6 , C 3 H 6 , C 3 H 8 , and CO 2 across a temperature range of 250–400 K for models of an archetypal polymer of intrinsic microporosity, PIM-1. A quadratic model is presented that captures two main mechanisms of temperature-dependent adsorption-induced deformation of PIM-1 up to a relative swelling of 1.15: thermal expansion and an increased propensity to swell as a function of species uptake. Two case studies are reported that highlight the critical role of operating temperature in industrial storage and separation applications. The first study focuses on methane storage and delivery applications using a pressure–temperature swing adsorption application (PTSA). We demonstrate that larger working capacities are accompanied by increased volumetric strain between adsorption–desorption steps. The second case study considers PIM-1 as an adsorbent to separate an exemplar ternary syngas mixture at operating temperatures ranging 300–550 K. Here, a temperature threshold of ~400 K is identified, beyond which adsorption-induced PIM-1 swelling is negligible and the solubility selectivity-loading curve transitions to exhibiting a nearly linear relationship.

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

Machine learning the metastable phase diagram of covalently bonded carbon

Abstract Conventional phase diagram generation involves experimentation to provide an initial estimate of the set of thermodynamically accessible phases and their boundaries, followed by use of phenomenological models to interpolate between the available experimental data points and extrapolate to experimentally inaccessible regions. Such an approach, combined with high throughput first-principles calculations and data-mining techniques, has led to exhaustive thermodynamic databases (e.g. compatible with the CALPHAD method), albeit focused on the reduced set of phases observed at distinct thermodynamic equilibria. In contrast, materials during their synthesis, operation, or processing, may not reach their thermodynamic equilibrium state but, instead, remain trapped in a local (metastable) free energy minimum, which may exhibit desirable properties. Here, we introduce an automated workflow that integrates first-principles physics and atomistic simulations with machine learning (ML), and high-performance computing to allow rapid exploration of the metastable phases to construct “metastable” phase diagrams for materials far-from-equilibrium. Using carbon as a prototypical system, we demonstrate automated metastable phase diagram construction to map hundreds of metastable states ranging from near equilibrium to far-from-equilibrium (400 meV/atom). We incorporate the free energy calculations into a neural-network-based learning of the equations of state that allows for efficient construction of metastable phase diagrams. We use the metastable phase diagram and identify domains of relative stability and synthesizability of metastable materials. High temperature high pressure experiments using a diamond anvil cell on graphite sample coupled with high-resolution transmission electron microscopy (HRTEM) confirm our metastable phase predictions. In particular, we identify the previously ambiguous structure of n -diamond as a cubic-analog of diaphite-like lonsdaelite phase.

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

Squeezing Out Interfacial Solvation: The Role of Hydrogen-Bonding in the Structural and Orientational Freedom of Molecular Self-Assembly

We report bioinspired membrane molecules with improved physical properties and enhanced stability can serve as functional models for conventional lipid or amphiphilic species. Importantly, these molecules can also provide new insights into emergent phenomena that manifest during self-assembly at interfaces. Here, we elucidate the structural response and mechanistic steps underlying the self-assembly of the amphiphilic, charged oligodimethylsiloxane imidazolium cation (ODMS-MIM+) at the air–aqueous interface using Langmuir trough methods with coincident surface-specific vibrational sum-frequency generation (SFG) spectroscopy. We find evidence for a new compression-induced desolvation step that precedes commonly known disordered-to-ordered phase transitions to form nanoscopic assemblies. The experimental data was supported by atomistic molecular dynamics (MD) simulations to provide a detailed mechanistic picture underlying the assembly and the role of water in these phase transitions. The sensitivity of the hydrophobic ODMS tail conformations to compression-owing to distinct water–ODMS interactions and tail–tail solvation properties-offers new strategies for the design of interfaces that can be further used to develop soft-matter electronics and low-dimensional materials using physical and chemical controls.

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

Topology induced crossover between Langevin, subdiffusion, and Brownian diffusion regimes in supercooled water

Despite extensive studies of supercooled water, it remains challenging to understand its peculiar dynamic anomalous properties. In this work, we integrated full atomistic simulations of supercooled water over the temperature range of room temperature to 200 K using a quantum-mechanics-based polarizable force field with the dressed dynamics method that couples fast collision events and slow reorganization dynamics of hydrogen-bond networks. Our analysis unveils the salient multiscale features in the transient relaxation dynamics of supercooled water. A classical Langevin behavior dominates at fast timescales, while long-time relaxations unveil two different activation barriers in two temperature regions: below and above 230 K. The modulation of the entropy spectrum by temperature is elucidated in terms of a three-state model underlined by the complexity of the water dynamics associated with a topological transition of a strong hydrogen-bond network. This state-dependent network topology is quantitatively characterized by power-law exponents of inverse network connectivity from 200 to 298 K. Finally, this work provides valuable guidance for further studies on the transient relaxation dynamics of supercooled water.

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

Impact of Irradiation-induced Lattice Defects and Dopants on Thermal Transport in Oxide Nuclear Fuels

Advanced nuclear reactor concepts aim to use fuels that must withstand unprecedented temperature and radiation extremes. In these fuels, thermal energy transport under irradiation is directly related to fuel longevity, reactor safety, and is arguably one of the most important performance metrics. Here we provide a comprehensive, first-principles-informed treatment of phonon mediated thermal transport in a defect-bearing actinide oxide with direct comparison to experimental measurements. Pristine and proton irradiated thorium dioxide was chosen as a model system to treat the complexity of thermal transport in the presence of lattice defects. A thermal transport model is implemented using the linearized Boltzmann transport equation (LBTE) with input from first principles calculations and defect evolution models. The output of the LBTE is compared directly to mesoscopic measurements of thermal conductivity on length scales commensurate with defect accumulation. Parametric measurements of conductivity with irradiation dose and temperature reveal a monotonic decrease in conductivity with irradiation dose, in the range of 0.001 dpa to 0.1 dpa. Besides radiation-induced defects, doping the thorium dioxide crystal with small quantities of uranium atoms also results in a drastic reduction in thermal conductivity. Temperature-dependent measurements of thermal conductivity in uranium-doped thorium dioxide single crystals suggests a resonant scattering mechanism is responsible for the observed reduction. This comprehensive, atomistic- to meso-scale treatment provides the necessary basis to investigate thermal transport under irradiation in more complex systems that exhibit strong electron correlation. This is a poster presentation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Atomistic Insight into Cation Effects on Binding Energies in Cu-Catalyzed Carbon Dioxide Reduction

Electrolyte properties in general and electrolyte ions in particular have been shown to have significant effects on the kinetics of electrochemical CO 2 and CO reduction at metal electrodes, but these effects have not yet been fully understood. Here, we investigate the effects of cations at the aqueous electrolyte–Cu(211) interface on adsorbate binding energies and the electrolyte structure using density functional theory (DFT). Charging the interface via explicit Na + has systematic effects on adsorbate–electrolyte interactions and conformations. We describe specific local adsorbate–ion interactions, including direct alkali ion–adsorbate coordination and hydrogen bonding via ion-coordinated water molecules. The relative importance of these specific interactions and purely electrostatic field–adsorbate interactions is investigated by comparing the DFT-calculated ion effects to those predicted by purely electrostatic models of the interface and the adsorbates. We find that the trend in ion effects among different adsorbates at constant surface charge density is well explained by a purely electrostatic interaction model. The binding energy of OCCHO is found to depend strongly on the surface charge density as well as the spatial distribution of ions at constant surface charge density. These effects are also explained by a purely electrostatic local electrostatic field–adsorbate model. This indicates that alkali ion effects can be mainly attributed to purely electrostatic field interactions and that the local field at ion-stabilized active sites can depend significantly on both the overall charge density and on the spatial distribution of ions at constant charge density. This work provides new insight on alkali ion effects on a variety of adsorbates relevant to the CO 2 and CO reduction reactions by describing specific local ion–adsorbate interactions, systematic changes in adsorbate–electrolyte interactions induced by explicit surface charge, and comparisons of the calculated effects on binding energies to simple electrostatic field–adsorbate models to explain trends among various adsorbates and at varying surface charge densities and supercell configurations.

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

Data-driven electron-diffraction approach reveals local short-range ordering in CrCoNi with ordering effects

Abstract The exceptional mechanical strength of medium/high-entropy alloys has been attributed to hardening in random solid solutions. Here, we evidence non-random chemical mixing in a CrCoNi alloy, resulting from short-range ordering. A data-mining approach of electron nanodiffraction enabled the study, which is assisted by neutron scattering, atom probe tomography, and diffraction simulation using first-principles theory models. Two samples, one homogenized and one heat-treated, are observed. In both samples, results reveal two types of short-range-order inside nanoclusters that minimize the Cr–Cr nearest neighbors (L1 2 ) or segregate Cr on alternating close-packed planes (L1 1 ). The L1 1 is predominant in the homogenized sample, while the L1 2 formation is promoted by heat-treatment, with the latter being accompanied by a dramatic change in dislocation-slip behavior. These findings uncover short-range order and the resulted chemical heterogeneities behind the mechanical strength in CrCoNi, providing general opportunities for atomistic-structure study in concentrated alloys for the design of strong and ductile materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning enabled discovery of superhard and ultrahard carbon polymorphs

The demand for multifunctional materials has motivated the move from near-equilibrium materials to metastable i.e. out-of-equilibrium phases that can meet several desired target properties. The search for such metastable phases with exotic properties is non-trivial and often serendipitous. Inverse design approaches based on evolutionary search have been powerful tools, but such traditional searches have focused on identifying primarily stable and metastable materials with the lowest enthalpy. The inverse design of materials, with a focus on a desired property such as, for example, hardness is a challenging task because of the expensive computational cost involved in sampling multiple structures. The recent advances in machine learning have brought new powerful AI techniques to the forefront which can potentially revolutionize the inverse design and discovery of materials, especially metastable phases capable of meeting multifunctionality. Here, in this work, we develop and apply an automated reinforcement learning workflow for inverse design that integrates first principles physics and atomistic simulations with machine learning (ML), and high-performance computing to allow rapid exploration of the superhard and ultrahard metastable phases of Carbon. We demonstrate an automatic machine learning based inverse design workflow to map new undiscovered metastable states ranging from near equilibrium to those far-from-equilibrium that satisfy multiple property objectives, specifically bulk moduli, shear moduli and hardness. We create a comprehensive library of carbon stable and metastable phases with varying hardness and subsequently shortlist 10 top performing candidate carbon structures, including two newly reported phases, based on their hardness and characterize their temperature dependent mechanical properties. A neural network model is built using featurization of allotropes of carbon to predict the quasi-harmonic Gibbs free energies. The Gibbs free energies of the top performing phases are analyzed to get an estimate of the experimental synthesizability of these superhard and ultrahard carbon phases. In general, we show using machine learning based inverse design approaches how hitherto inaccessible metastable states can be identified and potentially synthesized to meet the demand for multifunctional materials.

Balasubramanian, Karthik [Univ. of Illinois, Chica↗
DOE OSTI2023

Identifying crack tip position and stress intensity factors from displacement data

Fracture prognosis and characterization efforts require knowledge of crack tip position and the Stress Intensity Factors (SIFs) acting in the vicinity of the crack. Here, in this work, we present an efficient numerical approach to infer both of these characteristics under a consistent theoretical framework from noisy, unstructured displacement data. The novel approach utilizes the separability of the asymptotic linear elastic fracture mechanics fields to expedite the search for crack tip position and is particularly useful for noisy displacement data. The manuscript begins with an assessment of the importance of accurately locating crack tip position when quantifying the SIFs from displacement data. Next, the proposed separability approach for quickly inferring crack tip position is introduced. Comparing to the widely used displacement correlation approach, the performance of the separability approach is assessed. Cases involving both noisy data and systematic deviation from the asymptotic linear elastic fracture mechanics model are considered, e.g. inelastic material behavior and finite geometries. An open source python implementation of the proposed approach is available for use by those doing field and laboratory work involving digital image correlation and simulations, e.g. finite element, discrete element, molecular dynamics and peridynamics, where the crack tip position is not explicitly defined.

36 MATERIALS SCIENCE↗