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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

J-OCTA Applications to in Space Manufacturing and Nanofluidic Technologies

We consider examples of J-OCTA software application to physics-based analysis of in-space manufacturing technologies. Firstly, we consider an atomistic model of thermal welding at the polymer-polymer interface for applications to 3D manufacturing by following diffusion of polymer chains. We show that each component of the polymer blend has its own characteristic time of diffusion and that both strain–stress and shear viscosity curves agree with experimental data. Next, nonlinear shrinkage of the metal part during manufacturing by bound metal deposition under microgravity, is considered using a multi-scale physics-based approach that spans from atomistic dynamics to full-part shrinkage. At the smallest scale we use atomistic dynamics to simulated sintering of a few Ti6Al4V nanoparticles, study sintering mechanisms and estimate important parameters including grain boundary width and diffusion rates of the atoms at the surface that are used at the other scales of modelling. Finally, we analyse selective ionic conduction through an artificial nanopore in a graphene sheet motivated by applications to desalination and energy harvesting. We build distributions of the number and orientation of water molecules in the ion’s hydration shells providing novel insight into the ion-water-carbon interaction at the nanopore.

molecular dynamics↗
DOE OSTI2022

LANL Contributions to pyiron: An integrated development environment for computational materials science

To support and accelerate the development of simulation protocols in atomistic modelling, we introduce an integrated development environment (IDE) for computational materials science called pyiron (http://pyiron.org). The pyiron IDE combines a web based source code editor, a job management system for build automation, and a hierarchical data management solution. The core components of the pyiron IDE are pyiron objects based on an abstract class, which links application structures such as atomistic structures, projects, jobs, simulation protocols and computing resources with persistent storage and an interactive user environment. The simulation protocols within the pyiron IDE are constructed using the Python programming language. To highlight key concepts of this tool as well as to demonstrate its ability to simplify the implementation and testing of simulation protocols we discuss two applications. In these examples we show how pyiron supports the whole life cycle of a typical simulation, seamlessly combines ab initio with empirical potential calculations, and how complex feedback loops can be implemented. While originally developed with focus on ab initio thermodynamics simulations, the concepts and implementation of pyiron are general thus allowing to employ it for a wide range of simulation topics.

Janssen, Jan↗
NASA NTRS2004

adwTools Developed: New Bulk Alloy and Surface Analysis Software for the Alloy Design Workbench

A suite of atomistic modeling software, called the Alloy Design Workbench, has been developed by the Computational Materials Group at the NASA Glenn Research Center and the Ohio Aerospace Institute (OAI). The main goal of this software is to guide and augment experimental materials research and development efforts by creating powerful, yet intuitive, software that combines a graphical user interface with an operating code suitable for real-time atomistic simulations of multicomponent alloy systems. Targeted for experimentalists, the interface is straightforward and requires minimum knowledge of the underlying theory, allowing researchers to focus on the scientific aspects of the work. The centerpiece of the Alloy Design Workbench suite is the adwTools module, which concentrates on the atomistic analysis of surfaces and bulk alloys containing an arbitrary number of elements. An additional module, adwParams, handles ab initio input for the parameterization used in adwTools. Future modules planned for the suite include adwSeg, which will provide numerical predictions for segregation profiles to alloy surfaces and interfaces, and adwReport, which will serve as a window into the database, providing public access to the parameterization data and a repository where users can submit their own findings from the rest of the suite. The entire suite is designed to run on desktop-scale computers. The adwTools module incorporates a custom OAI/Glenn-developed Fortran code based on the BFS (Bozzolo- Ferrante-Smith) method for alloys, ref. 1). The heart of the suite, this code is used to calculate the energetics of different compositions and configurations of atoms.

Bozzolo, Guillermo↗
DOE OSTI2021

Probing the electronic properties of the electrified silicon/water interface by combining simulations and experiments.

Silicon (Si) is broadly used in electrochemical and photoelectrochemical devices, where the capacitive and Faradaic reactions at the Si/water interfaces are critical for signal transduction or noise generation. However, probing the electrified Si/water interface at the microscopic level remains a challenging task. Here we focus on hydrogenated Si surfaces in contact with water, relevant to transient electronics and photoelectrochemical modulation of biological cells and tissues. We show that by carrying out first principles molecular dynamics simulations of the Si(100)/water interface in the presence of an electric field we can realistically correlate the computed flat-band potential and tunneling current images at the interface with experimentally measured capacitive and Faradaic currents. Specifically, we validate our simulations in the presence of bias by performing pulsed chronoamperometry measurements on Si wafers in solution. Consistent with prior experiments, our measurements and simulations indicate the presence of voltage-dependent capacitive currents at the interface. We also find that Faradaic currents are weakly dependent on the applied bias, which we relate to surface defects present in newly prepared samples.

Ye, Zifan↗
NASA NTRS2003

Alloy Interface Interdiffusion Modeled

With renewed interest in developing nuclear-powered deep space probes, attention will return to improving the metallurgical processing of potential nuclear fuels so that they remain dimensionally stable over the years required for a successful mission. Previous work on fuel alloys at the NASA Glenn Research Center was primarily empirical, with virtually no continuing research. Even when empirical studies are exacting, they often fail to provide enough insight to guide future research efforts. In addition, from a fundamental theoretical standpoint, the actinide metals (which include materials used for nuclear fuels) pose a severe challenge to modern electronic-structure theory. Recent advances in quantum approximate atomistic modeling, coupled with first-principles derivation of needed input parameters, can help researchers develop new alloys for nuclear propulsion.

Bozzolo, Guillermo H.↗
DOE OSTI2025

Nonadiabatic Force Matching for Alchemical Free-Energy Estimation

We propose a method to compute free-energy differences from nonadiabatic alchemical transformations by using flow-based generative models. The method, nonadiabatic force matching, hinges on estimating the dissipation along an alchemical switching process in terms of a nonadiabatic force field that can be learned through stochastic flow matching. The learned field can be used in conjunction with short-time trajectory data to evaluate upper and lower bounds on the alchemical free energy that variationally converge to the exact value if the field is optimal. Applying the method to evaluate the alchemical free energy of atomistic models shows that it can substantially reduce the simulation cost of a free-energy estimate at a negligible loss of accuracy when compared with thermodynamic integration.

Computational chemistry↗
DOE OSTI2024

Predicting RNA structure and dynamics with deep learning and solution scattering

Advanced deep learning and statistical methods can predict structural models for RNA molecules. However, RNAs are flexible, and it remains difficult to describe their macromolecular conformations in solutions where varying conditions can induce conformational changes. Small-angle x-ray scattering (SAXS) in solution is an efficient technique to validate structural predictions by comparing the experimental SAXS profile with those calculated from predicted structures. There are two main challenges in comparing SAXS profiles to RNA structures: the absence of cations essential for stability and charge neutralization in predicted structures and the inadequacy of a single structure to represent RNA’s conformational plasticity. We introduce a solution conformation predictor for RNA (SCOPER) to address these challenges. This pipeline integrates kinematics-based conformational sampling with the innovative deep learning model, IonNet, designed for predicting Mg 2+ ion binding sites. Validated through benchmarking against 14 experimental data sets, SCOPER significantly improved the quality of SAXS profile fits by including Mg 2+ ions and sampling of conformational plasticity. We observe that an increased content of monovalent and bivalent ions leads to decreased RNA plasticity. Therefore, carefully adjusting the plasticity and ion density is crucial to avoid overfitting experimental SAXS data. SCOPER is an efficient tool for accurately validating the solution state of RNAs given an initial, sufficiently accurate structure and provides the corrected atomistic model, including ions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Interfacial Defect of Lithium Metal in Solid‐State Batteries

Abstract All‐solid‐state battery with Li metal anode is a promising rechargeable battery technology with high energy density and improved safety. Currently, the application of Li metal anode is plagued by the failure at the interfaces between lithium metal and solid electrolyte (SE). However, little is known about the defects at Li–SE interfaces and their effects on Li cycling, impeding further improvement of Li metal anodes. Herein, by performing large‐scale atomistic modeling of Li metal interfaces with common SEs, we discover that lithium metal forms an interfacial defect layer of nanometer‐thin disordered lithium at the Li–SE interfaces. This interfacial defect Li layer is highly detrimental, leading to interfacial failure such as pore formation and contact loss during Li stripping. By systematically studying and comparing incoherent, coherent, and semi‐coherent Li–SE interfaces, we find that the interface with good lattice coherence has reduced Li defects at the interface and has suppressed interfacial failure during Li cycling. Our finding discovered the critical roles of atomistic lithium defects at interfaces for the interfacial failure of Li metal anode, and motivates future atomistic‐level interfacial engineering for Li metal anode in solid‐state batteries.

Yang, Menghao↗
DOE OSTI2021

Interfacial Defect of Lithium Metal in Solid-State Batteries

All-solid-state battery with Li metal anode is a promising rechargeable battery technology with high energy density and improved safety. Currently, the application of Li metal anode is plagued by the failure at the interfaces between lithium metal and solid electrolyte (SE). However, little is known about the defects at Li–SE interfaces and their effects on Li cycling, impeding further improvement of Li metal anodes. Herein, by performing large-scale atomistic modeling of Li metal interfaces with common SEs, we discover that lithium metal forms an interfacial defect layer of nanometer-thin disordered lithium at the Li–SE interfaces. Further, this interfacial defect Li layer is highly detrimental, leading to interfacial failure such as pore formation and contact loss during Li stripping. By systematically studying and comparing incoherent, coherent, and semi-coherent Li–SE interfaces, we find that the interface with good lattice coherence has reduced Li defects at the interface and has suppressed interfacial failure during Li cycling. Our finding discovered the critical roles of atomistic lithium defects at interfaces for the interfacial failure of Li metal anode, and motivates future atomistic-level interfacial engineering for Li metal anode in solid-state batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Dissociated prismatic loop punching by bubble growth in FCC metals

Abstract Materials performance can be significantly degraded due to bubble generation. In this work, the bubble growth process is elaborated in Cu by atomistic modeling to bridge the gap of experimental observations. Upon continuous He implantation, bubble growth is accommodated first by nucleation of dislocation network from bubble surface, then formation of dissociated prismatic dislocation loop (DPDL), and final DPDL emission in $$\langle 110\rangle$$ ⟨ 110 ⟩ directions. As the DPDL is found capable of collecting He atoms, this process is likely to assist the formation of self-organized bubble superlattice, which has been reported from experiments. Moreover, the pressurized bubble in solid state manifests the shape of an imperfect octahedron, built by Cu $$\{111\}$$ { 111 } surfaces, consistent with experiments. These atomistic details integrating experimental work fill the gap of mechanistic understanding of athermal bubble growth in Cu. Importantly, by associating with nanoindentation testings, DPDL punching by bubble growth arguably applies to various FCC (face-centered cubic) metals such as Au, Ag, Ni, and Al.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A theory of localized excitations in supercooled liquids

Here, a new connection between the structure and dynamics in glass-forming liquids is presented. We show how the origin of spatially localized excitations, as defined by the dynamical facilitation (DF) theory, can be understood from a structure-based framework. This framework is constructed by associating excitation events in the DF theory to hopping events between energy minima in the potential energy landscape (PEL). By reducing the PEL to an equal energy well picture and applying a harmonic approximation, we develop a field theory to describe elastic fluctuations about inherent states, which are energy minimizing configurations of the PEL. We model an excitation as a shear transformation zone (STZ) inducing a localized pure shear deformation onto an inherent state. We connect STZs to T1 transition events that break the elastic bonds holding the local structure of an inherent state. A formula for the excitation energy barrier, denoted as J σ , is obtained as a function of inherent-state elastic moduli and the radial distribution function. The energy barrier from the current theory is compared to the one predicted by the DF theory where good agreement is found in various two-dimensional continuous poly-disperse atomistic models of glass formers. These results strengthen the role of structure and elasticity in driving glassy dynamics through the creation and relaxation of localized excitations.

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

Oscillatory redox behavior in oxides: Cyclic surface reconstruction and reactivity modulation via the Mars–van Krevelen mechanism

The breaking of translational symmetry at oxide surfaces gives rise to coordinatively unsaturated cations/anions and surface restructuring—key factors that govern surface reactivity. Using direct in situ environmental transmission electron microscopy (TEM) observations along with atomistic modeling, we report oscillatory redox behavior in CuO under H 2 , where cyclic surface reconstruction and reactivity modulation occur via the Mars–van Krevelen (MvK) mechanism. We observe self-switching between oxygen-rich and oxygen-deficient surface reconstructions, alternately activating and deactivating the surface for H 2 O formation. During periods of chemical inactivity, the oxygen-deficient surface undergoes slow reoxidation via lattice oxygen diffusing from subsurface and bulk reservoirs, restoring the active oxygen-rich surface termination. The inherent disparity in chemical activity among undercoordinated surface ions, along with sluggish subsurface-to-surface oxygen replenishment, drives this oscillatory redox cycle, modulating H2-induced loss of lattice oxygen at the surface and its delayed replenishment from the subsurface. This creates spatiotemporally separated redox steps at the oxide surface. The phenomena and atomistic insights presented here have significant implications for manipulating the surface reactivity of oxides by tuning the separation of these redox steps.

36 MATERIALS SCIENCE↗
NASA NTRS2005

Computational Materials Program for Alloy Design

The research program sponsored by this grant, "Computational Materials Program for Alloy Design", covers a period of time of enormous change in the emerging field of computational materials science. The computational materials program started with the development of the BFS method for alloys, a quantum approximate method for atomistic analysis of alloys specifically tailored to effectively deal with the current challenges in the area of atomistic modeling and to support modern experimental programs. During the grant period, the program benefited from steady growth which, as detailed below, far exceeds its original set of goals and objectives. Not surprisingly, by the end of this grant, the methodology and the computational materials program became an established force in the materials communitiy, with substantial impact in several areas. Major achievements during the duration of the grant include the completion of a Level 1 Milestone for the HITEMP program at NASA Glenn, consisting of the planning, development and organization of an international conference held at the Ohio Aerospace Institute in August of 2002, finalizing a period of rapid insertion of the methodology in the research community worlwide. The conference, attended by citizens of 17 countries representing various fields of the research community, resulted in a special issue of the leading journal in the area of applied surface science. Another element of the Level 1 Milestone was the presentation of the first version of the Alloy Design Workbench software package, currently known as "adwTools". This software package constitutes the first PC-based piece of software for atomistic simulations for both solid alloys and surfaces in the market.Dissemination of results and insertion in the materials community worldwide was a primary focus during this period. As a result, the P.I. was responsible for presenting 37 contributed talks, 19 invited talks, and publishing 71 articles in peer-reviewed journals, as detailed later in this Report.

Bozzolo, Guillermo↗
DOE OSTI2022

Tackling Disorder in γ‐Ga 2 O 3

Abstract Ga 2 O 3 and its polymorphs are attracting increasing attention. The rich structural space of polymorphic oxide systems such as Ga 2 O 3 offers potential for electronic structure engineering, which is of particular interest for a range of applications, such as power electronics. γ‐Ga 2 O 3 presents a particular challenge across synthesis, characterization, and theory due to its inherent disorder and resulting complex structure–electronic‐structure relationship. Here, density functional theory is used in combination with a machine‐learning approach to screen nearly one million potential structures, thereby developing a robust atomistic model of the γ‐phase. Theoretical results are compared with surface and bulk sensitive soft and hard X‐ray photoelectron spectroscopy, X‐ray absorption spectroscopy, spectroscopic ellipsometry, and photoluminescence excitation spectroscopy experiments representative of the occupied and unoccupied states of γ‐Ga 2 O 3 . The first onset of strong absorption at room temperature is found at 5.1 eV from spectroscopic ellipsometry, which agrees well with the excitation maximum at 5.17 eV obtained by photoluminescence excitation spectroscopy, where the latter shifts to 5.33 eV at 5 K. This work presents a leap forward in the treatment of complex, disordered oxides and is a crucial step toward exploring how their electronic structure can be understood in terms of local coordination and overall structure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamic Interfacial Architectures: Cruciferin‐Stabilized Oil/Water Interfaces for Sustainable Emulsions

Stabilizing oil-water interfaces in emulsions by plant-based proteins provides sustainable and tunable ways for designing emulsions with specific properties, for food, healthcare, and pharmaceuticals. Cruciferin, a protein from rapeseed, has great potential as green emulsifier, but details about its structure and mobility at oil-water interfaces are largely unknown. Here, these properties are studied with small angle neutron and x-ray scattering, and neutron spin echo spectroscopy, analyzed by atomistic modelling of scattering curves and coarse-grained modelling, to gain insight into interface coverage, and molecular conformation and mobility at the interface. Cruciferin assumes trimeric conformations at the interface, as in solution, but with its protrusions from the central core of the subunits (“arms”) more compressed. Interfacial mobility is only marginally lower than in solution, indicating the arms still transiently extend and preserve a network, for the first time revealing the mechanism how cruciferin forms highly elastic 2d gel-like oil-water interfaces, as observed in macroscopic rheology. The high interfacial mobility may help in self-repairing non-stabilized interfacial fractions, reducing coalescence. These findings provide a deeper molecular level understanding of proteins at oil-water interfaces, which can stimulate development of new plant-based emulsion products, and contribute to the global protein transition.

cruciferin↗
DOE OSTI2024

Molecular‐Level Insight into the Chlorofluorocarbons Adsorption by Defective Covalent Organic Polymers

Abstract Halocarbons have important industrial applications, but because of their contribution to global warming and the fact that they can cause ozone depletion, they are considered highly toxic. Hence, the techniques that can capture and recover the used halocarbons with energy‐efficient methods have been recently received greater attention. In this contribution, we report the capture of dichlorodifluoromethane (R12), which has high global warming and ozone depletion potential, using covalent organic polymers (COPs). The defect‐engineered COPs were synthesized and demonstrated outstanding sorption capacities, ~226 wt % of R12 combined with linear‐shaped adsorption isotherms. We further identified the plausible microscopic adsorption mechanism of the investigated COPs via grand canonical Monte Carlo simulations applied to non‐defective and a collection of atomistic models of the defective COPs. The modeling work suggests that significant R12 adsorption performance is attributed to a gradual increment of porosities due to isolated/interconnected micro‐/meso‐pore channels and the change of the long‐range ordering of both COPs. The successive hierarchical‐pore‐filling mechanism promotes R12 molecular adsorption via moderate van der Waals adsorbate‐adsorbent interactions in the micropores of both COPs at low pressure followed by adsorbate‐adsorbate interactions in the extra‐voids created at moderate to high pressure ranges. This continuous pore‐filling mechanism makes defective COPs as promising sorbents for halocarbon adsorption.

Shen, Jian↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

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

Computational analysis of the tryptophan cation radical energetics in peroxidase Compound $\mathrm{I}$

Three well-characterized heme peroxidases (cytochrome c peroxidase = CCP, ascorbate peroxidase = APX, and Leishmania major peroxidase = LMP) all have a Trp residue tucked under the heme stacked against the proximal His heme ligand. The reaction of peroxidases with H 2 O 2 to give Compound I results in the oxidation of this Trp to a cationic radical in CCP and LMP but not in APX. Considerable experimental data indicate that the local electrostatic environment controls whether this Trp or the porphyrin is oxidized in Compound I. Attempts have been made to place the differences between these peroxidases on a quantitative basis using computational methods. These efforts have been somewhat limited by the approximations required owing to the computational cost of using fully solvated atomistic models with well-developed forcefields. This now has changed with available GPU computing power and the associated development of software. Here we employ thermodynamic integration and multistate Bennett acceptance ratio methods to help fine-tune our understanding on the energetic differences in Trp radical stabilization in all three peroxidases. These results indicate that the local solvent structure near the redox active Trp plays a significant role in stabilization of the cationic Trp radical.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗