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

Kinetic Growth of Multicomponent Microcompartment Shells

An important goal of systems and synthetic biology is to produce high value chemical species in large quantities. Microcompartments, which are protein nanoshells encapsulating catalytic enzyme cargo, could potentially function as tunable nanobioreactors inside and outside cells to generate these high value species. Modifying the morphology of microcompartments through genetic engineering of shell proteins is one viable strategy to tune cofactor and metabolite access to encapsulated enzymes. However, this is a difficult task without understanding how changing interactions between the many different types of shell proteins and enzymes affect microcompartment assembly and shape. Here, we use multiscale molecular dynamics and experimental data to describe assembly pathways available to microcompartments composed of multiple types of shell proteins with varied interactions. As the average interaction between the enzyme cargo and the multiple types of shell proteins is weakened, the shell assembly pathway transitions from (i) nucleating on the enzyme cargo to (ii) nucleating in the bulk and then binding the cargo as it grows to (iii) an empty shell. Atomistic simulations and experiments using the 1,2-propanediol utilization microcompartment system demonstrate that shell protein interactions are highly varied and consistent with our multicomponent, coarse-grained model. Furthermore, our results suggest that intrinsic bending angles control the size of these microcompartments. Altogether, our simulations and experiments provide guidance to control microcomparmtent size and assembly by modulating the interactions between shell proteins.

defects↗
DOE OSTI2026

Extreme Nanoconfinement Dramatically Enhances Small Molecule Solubility in Nonpolar Polymers

Elucidating gas solubility in confined polymer systems addresses a fundamental gap in polymer physics and has important implications for gas barrier and separation technologies and polymer upcycling reactions. In this study, we examine the solubility of methanol and n-hexane in polystyrene and low-density polyethylene confined within the interstitial pores of disordered silica nanoparticle packings. Using capillary rise infiltration, these polymers are infiltrated into the nanoparticle packings. Gas solubility in these confined polymers is measured using a quartz crystal microbalance. Remarkably, confinement leads to a 10- to ∼100-fold increase in gas solubility. Systematic experiments reveal that pore size plays a dominant role in increasing solubility. In contrast, the molecular weight of the polymer and the surface wetting characteristics of the nanoparticles, achieved by hydrophobically modifying the nanoparticles, have minimal effects. By systematically varying penetrant polarity, polymer crystallinity, pore size, and nanopore surface chemistry, this study isolates geometric nanoconfinement as the dominant factor governing solubility enhancement in confined polymers. Atomistic simulations revealed that confinement and surface−polymer interactions both contribute to enhanced solubility, with polymer packing playing a significant role in modulating gas uptake. Here, these results suggest that the changes in the molecular arrangements of polymer segments underlie the observed trends. This study highlights the potential of confined polymers in engineering the separation performance of membranes and heterogeneously catalyzed polymer upcycling reactions.

alcohols↗
DOE OSTI2025

Understanding (La,Sr)(Co,Fe)O$_{3-δ}$ Phase Instability within SOECs Using a Combined Experimental and Atomistic Modeling Approach

Understanding the onset of degradation in the air electrode within solid oxide electrolysis cells (SOECs), and the subsequent impact on cell performance, is a critical step in mitigating the performance losses and stability issues of SOECs. In an effort to identify early onset degradation phenomena, SOECs were characterized as fabricated and after testing potentiostatically at 1.3 V for 1000 h at 750 °C. SOEC air electrodes composed of a 1:1 composite of La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3–δ (6428-LSCF) and Gd 0.1 Ce 0.9 O 1.95 (GDC) were studied using synchrotron X-ray diffraction (XRD), scanning transmission electron microscopy coupled with energy dispersive X-ray spectroscopy (STEM-EDS), and X-ray absorption near-edge spectroscopy (XANES) to evaluate the changes in the air electrode structurally and chemically. These techniques show the migration of Sr species from the air electrode through pores in the GDC barrier layer, progressing to the electrolyte boundary, where it accumulates and reacts with (Zr 0.84 Y 0.16 )O 2–δ (YSZ) to form SrZrO 3 . Microscopy results are paired with atomistic simulations to better understand the relationship between the thermodynamic instability of 6428-LSCF and cell fabrication/testing conditions. First-principles calculations reveal that LSCF-6428 is not stable during cell manufacturing and testing conditions, which supports the experimental identification of secondary phases in both as-fabricated and tested cells. Together, these results demonstrate that the challenging environments encountered by SOECs during cell manufacturing and operation lead to instabilities of the target 6428-LSCF anode material and underscore the need for more durable, high-performing SOEC components.

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

Electrochemical Reduction Pathways from Goethite to Green Iron in Alkaline Solution with Silicate Additive

Energy-efficient and low-temperature iron electrolysis in alkaline solutions is a low-cost and sustainable ironmaking process with zero-carbon emissions when renewable electrical sources are involved. However, its implementation is hindered by electrochemically inert Fe 3 O 4 and parasitic H 2 gas formation during the electrochemical reduction process, resulting in the low energy efficiency of iron electrolysis. Here, we further explore the potential of electrochemical reduction of goethite (FeOOH) by employing a low concentration of silicate additive in an alkaline solution to mitigate Fe 3 O 4 accumulation and H 2 generation. Electrochemical measurements coupled with operando X-ray diffraction and X-ray absorption spectroscopy suggested FeOOH → Fe 3 O 4 → Fe(OH) 2 → Fe reduction pathways. Interestingly, a poorly crystalline or amorphous Fe(OH) 2 phase formed in the NaOH/silicate mixed electrolyte, possibly due to the inhibitive effect of silicate on water and ion transport, which eventually contributed to the improved reduction of Fe 3 O 4 , also supported by atomistic simulations. This work demonstrates the potential for silicate as a low-cost and effective electrolyte additive to improve room-temperature green iron formation via electrolysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhanced Sampling Aided Design of Molecular Photoswitches

Advances in the evolving field of atomistic simulations promise important insights for the design and fundamental understanding of novel molecular photoswitches. Here, in this study, we use state-of-the-art enhanced simulation techniques to unravel the complex, multistep chemistry of donor–acceptor Stenhouse adducts (DASAs). Our reaction discovery workflow consists of enhanced sampling for efficient chemical space exploration, refinement of newly observed pathways with more accurate ab initio electronic structure calculations, and structural modifications to introduce design principles within future generations of DASAs. We showcase our discovery workflow by not only recovering the full photoswitching mechanism of DASA but also predicting a plethora of new plausible thermal pathways and suggesting a way for their experimental validation. Furthermore, we illustrate the tunability of these newly discovered reactions, leading to a potential avenue for controlling DASA dynamics through multiple external stimuli. Overall, these insights could offer alternative routes to increase the efficiency and control of DASA’s photoswitching mechanism, providing new elements to design more complex light-responsive materials.

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

Atomic-Scale Dynamics of Five-Fold Twin Mediated Coalescence: Pathway-Dependent and Defect-Governed Nonclassical Growth Mechanisms

Defective crystals with distinct properties have been discovered in many systems. However, the growth mechanism of defective crystals is still poorly understood. Here, in this work, using a 5-fold twinned gold (Au) nanocrystal (NC) as a model system, three new coalescence pathways involving detwinning or twinning have been uncovered through atomic-scale dynamic observations in an aberration-corrected transmission electron microscope coupled with atomistic simulations. This demonstrates that beyond crystal size, coalescence growth dynamics involving 5-fold twins (5-FTs) are highly dependent on crystal defect density and the approach pathways of the crystals. When a 5-FT encounters a smaller 5-FT or a smaller NC in a face-to-face way, a new, larger 5-FT is produced at a relatively fast coalescence growth rate; while in a corner-to-corner way, the coalescence dynamics are more retarded and sluggish, which is conducive to the formation of complex multitwined structures rather than 5-FTs. This highlights that the planar defect density and crystal approach pathway influence the coalescence dynamics of crystals containing 5-FT. Moreover, a column-by-column grain boundary (GB) migration mechanism, which results in bent GBs, was also unveiled during the crystal coalescence process. These results enrich the general understanding of the crystal growth theory and provide new insights into the controllable fabrication of 5-FTs by crystal coalescence mechanisms.

crystallization↗
DOE OSTI2022

Molecular-Level Control over Ionic Conduction and Ionic Current Direction by Designing Macrocycle-Based Ionomers

Poor ionic conductivity of the catalyst-binding, submicrometer- thick ionomer layers in energy conversion and storage devices is a huge challenge. However, ionomers are rarely designed keeping in mind the specific issues associated with nanoconfinement. Here, we designed nature-inspired ionomers (calix-2) having hollow, macrocyclic, calix[4]arene-based repeat units with precise, sub-nanometer diameter. In ≤100 nm-thick films, the in-plane proton conductivity of calix-2 was up to 8 times higher than the current benchmark ionomer Nafion at 85% relative humidity (RH), while it was 1-2 orders of magnitude higher than Nafion at 20-25% RH. Confocal laser scanning microscopy and other synthetic techniques allowed us to demonstrate the role of macrocyclic cavities in boosting the proton conductivity. The systematic self-assembly of calix-2 chains into ellipsoids in thin films was evidenced from atomic force microscopy and grazing incidence small-angle X-ray scattering measurements. Moreover, the likelihood of alignment and stacking of macrocyclic units, the presence of one-dimensional water wires across this macrocycle stacks, and thus the formation of long-range proton conduction pathways were suggested by atomistic simulations. We not only did see an unprecedented improvement in thin-film proton conductivity but also saw an improvement in proton conductivity of bulk membranes when calix-2 was added to the Nafion matrices. Nafion-calix-2 composite membranes also took advantage of the asymmetric charge distribution across calix[4]arene repeat units collectively and exhibited voltage-gating behavior. The inclusion of molecular macrocyclic cavities into the ionomer chemical structure can thus emerge as a promising design concept for highly efficient ion-conducting and ion-permselective materials for sustainable energy applications.

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

A Universal Model of Cation Effects in Electrocatalysis

Electrolyte cations are conventionally viewed as inert spectators in electrocatalysis. However, a wealth of observations show that catalytic rates are often highly sensitive to cation identity. Despite their prevalence, these cation effects have resisted a unified mechanistic explanation, with different physical phenomena implicated across reaction chemistries, catalyst compositions, and choice of solvent. In this perspective, we describe a general framework for understanding cation effects in electrocatalysis based on electrostatics. We argue that cations influence reaction rates by modifying the strength of the electric field present at the catalyst surface, which alters the energetics of adsorbed intermediates and transition states according to their dipole moments and polarizabilities. The magnitude of this field depends on how cations arrange at the electrode surface, controlled by their size, shape, solvation, and packing efficiency. Cations that can arrange more densely result in a steeper potential drop at the electrode surface and consequently a stronger electric field. Our model further identifies two criteria for observing cation effects: (1) the operating potential must be negative of the electrode’s potential of zero total charge, ensuring that cations accumulate at the interface, and (2) the energetics of the kinetically relevant elementary step must be field sensitive. This framework reconciles previously inconsistent trends, including why cation effects appear only for some catalysts, why reaction selectivity is sensitive to cation identity, and why activity can increase with cation size on certain metals but decrease on others. Supported by kinetic measurements, spectroscopy, and atomistic simulations, the model provides both conceptual value for building intuition about catalysis at charged interfaces and predictive value for anticipating trends for new reactions, catalysts, and electrolytes. We conclude by highlighting the importance of electric fields across electrochemical, thermochemical, and biological catalysis and propose that considering the electrostatic environment around active sites offers new opportunities for improving activity and selectivity.

cation effects↗
DOE OSTI2022

Antisite defect qubits in monolayer transition metal dichalcogenides

Abstract Being atomically thin and amenable to external controls, two-dimensional (2D) materials offer a new paradigm for the realization of patterned qubit fabrication and operation at room temperature for quantum information sciences applications. Here we show that the antisite defect in 2D transition metal dichalcogenides (TMDs) can provide a controllable solid-state spin qubit system. Using high-throughput atomistic simulations, we identify several neutral antisite defects in TMDs that lie deep in the bulk band gap and host a paramagnetic triplet ground state. Our in-depth analysis reveals the presence of optical transitions and triplet-singlet intersystem crossing processes for fingerprinting these defect qubits. As an illustrative example, we discuss the initialization and readout principles of an antisite qubit in WS 2 , which is expected to be stable against interlayer interactions in a multilayer structure for qubit isolation and protection in future qubit-based devices. Our study opens a new pathway for creating scalable, room-temperature spin qubits in 2D TMDs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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

Dual phase patterning during a congruent grain boundary phase transition in elemental copper

The phase behavior of grain boundaries can have a strong influence on interfacial properties. Little is known about the emergence of grain boundary phases in elemental metal systems and how they transform. Here, we observe the nanoscale patterning of a grain boundary by two alternating grain boundary phases with distinct atomic structures in elemental copper by atomic resolution imaging. The same grain boundary phases are found by computational grain boundary structure search indicating a first-order transformation. Finite temperature atomistic simulations reveal a congruent, diffusionless transition between these phases under ambient pressure. The patterning of the grain boundary at room temperature is dominated by the grain boundary phase junctions separating the phase segments. Our analysis suggests that the reduced mobility of the phase junctions at low temperatures kinetically limits the transformation, but repulsive elastic interactions between them and disconnections could additionally stabilize the pattern formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ubiquitous short-range order in multi-principal element alloys

Recent research in multi-principal element alloys (MPEAs) has increasingly focused on the role of short-range order (SRO) on material performance. However, the mechanisms of SRO formation and its precise control remain elusive, limiting the progress of SRO engineering. Here, leveraging advanced additive manufacturing techniques that produce samples with a wide range of cooling rates (up to 10 7 K s –1 ) and an enhanced semi-quantitative electron microscopy method, we characterize SRO in three CoCrNi-based face-centered-cubic (FCC) MPEAs. Surprisingly, irrespective of the processing and thermal treatment history, all samples exhibit similar levels of SRO. Atomistic simulations reveal that during solidification, prevalent local chemical order arises in the liquid-solid interface (solidification front) even under the extreme cooling rate of 10 11 K s –1 . This phenomenon stems from the swift atomic diffusion in the supercooled liquid, which matches or even surpasses the rate of solidification. Therefore, SRO is an inherent characteristic of most FCC MPEAs, insensitive to variations in cooling rates and even annealing treatments typically available in experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Intermittent cluster dynamics and temporal fractional diffusion in a bulk metallic glass

Glassy solids evolve towards lower-energy structural states by physical aging. This can be characterized by structural relaxation times, the assessment of which is essential for understanding the glass’ time-dependent property changes. Conducted over short times, a continuous increase of relaxation times with time is seen, suggesting a time-dependent dissipative transport mechanism. By focusing on micro-structural rearrangements at the atomic-scale, we demonstrate the emergence of sub-diffusive anomalous transport and therefore temporal fractional diffusion in a metallic glass, which we track via coherent x-ray scattering conducted over more than 300,000 s. At the longest probed decorrelation times, a transition from classical stretched exponential to a power-law behavior occurs, which in concert with atomistic simulations reveals collective and intermittent atomic motion. Our observations give a physical basis for classical stretched exponential relaxation behavior, uncover a new power-law governed collective transport regime for metallic glasses at long and practically relevant time-scales, and demonstrate a rich and highly non-monotonous aging response in a glassy solid, thereby challenging the common framework of homogeneous aging and atomic scale diffusion.

42 ENGINEERING↗
DOE OSTI2025

Terahertz-field activation of polar skyrons

Unraveling collective modes arising from coupled degrees of freedom is crucial for understanding complex interactions in solids and developing new functionalities. Unique collective behaviors emerge when two degrees of freedom, ordered on distinct length scales, interact. Polar skyrmions, three-dimensional electric polarization textures in ferroelectric superlattices, disrupt the lattice continuity at the nanometer scale with nontrivial topology, leading to previously unexplored collective modes. Here, using terahertz-field excitation and femtosecond x-ray diffraction, we discover subterahertz collective modes, dubbed “skyrons”, which appear as swirling patterns of atomic displacements functioning as atomic-scale gearsets. The key to activating skyrons is the use of the THz field that couples primarily to skyrmion domain walls. Momentum-resolved time-domain measurements of diffuse scattering reveal an avoided crossing in the dispersion relation of skyrons. Atomistic simulations and dynamical phase-field modeling provide microscopic insights into the three-dimensional crystallographic and polarization dynamics. The amplitude and dispersion of skyrons are demonstrated to be controlled by sample temperature and electric-field bias. The discovery of skyrons and their coupling with terahertz fields opens avenues for ultrafast control of topological polar structures.

ferroelectrics↗
DOE OSTI2025

Structural basis of undecaprenyl phosphate glycosylation leading to polymyxin resistance in Gram-negative bacteria

In Gram-negative bacteria, the enzymatic modification of Lipid A with aminoarabinose (L-Ara4N) leads to resistance against polymyxin antibiotics and cationic antimicrobial peptides. ArnC, an integral membrane glycosyltransferase, attaches a formylated form of aminoarabinose to the lipid undecaprenyl phosphate, enabling its association with the bacterial inner membrane. Here, we present cryo-electron microscopy structures of ArnC from S. enterica in apo and nucleotide-bound conformations. These structures reveal a conformational transition that takes place upon binding of the partial donor substrate. Using coarse-grained and atomistic simulations, we provide insights into substrate coordination before and during catalysis, and we propose a catalytic mechanism that may operate on all similar metal-dependent polyprenyl phosphate glycosyltransferases. The reported structures provide a new target for drug design aiming to combat polymyxin resistance.

Ashraf, Khuram U. [Rutgers University, Newark, NJ ↗
DOE OSTI2025

Quantification and prediction of solidification textures under additive manufacturing conditions

Crystallographic textures are a major determinant of the macroscale anisotropic properties of polycrystalline metallic alloys produced in a wide range of additive manufacturing (AM) processes. Here, we introduce a statistical method that can accurately quantify the degree of orientational order of textures despite the large random fluctuations in the orientation of individual grains inherent in AM processes. The method, demonstrated for laser and resolidification of AlSi thin films, extends Z-scoring to a dynamical regime to assess the statistical significance of observed textures compared to randomly generated ones at different stages of solidification. We further show that, combined with phase-field modeling, this method can be used to infer fundamental anisotropic properties of the solid-liquid interface that are essential for texture prediction, and are compared here to the results of atomistic simulations. In addition, phase-field modeling reveals that, even at rapid AM solidification rates, the observed 〈110〉-dominated textures in the AlSi thin films are controlled predominantly by the anisotropy of the interface free-energy and sheds light on the physical mechanism of grain competition. These results significantly enhance both the existing tools for the quantification and prediction of AM crystallographic textures and our basic understanding of their formation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multiscale computational understanding and growth of 2D materials: a review

Abstract The successful discovery and isolation of graphene in 2004, and the subsequent synthesis of layered semiconductors and heterostructures beyond graphene have led to the exploding field of two-dimensional (2D) materials that explore their growth, new atomic-scale physics, and potential device applications. This review aims to provide an overview of theoretical, computational, and machine learning methods and tools at multiple length and time scales, and discuss how they can be utilized to assist/guide the design and synthesis of 2D materials beyond graphene. We focus on three methods at different length and time scales as follows: (i) nanoscale atomistic simulations including density functional theory (DFT) calculations and molecular dynamics simulations employing empirical and reactive interatomic potentials; (ii) mesoscale methods such as phase-field method; and (iii) macroscale continuum approaches by coupling thermal and chemical transport equations. We discuss how machine learning can be combined with computation and experiments to understand the correlations between structures and properties of 2D materials, and to guide the discovery of new 2D materials. We will also provide an outlook for the applications of computational approaches to 2D materials synthesis and growth in general.

Momeni, Kasra↗
DOE OSTI2021

Decoding defect statistics from diffractograms via machine learning

Abstract Diffraction techniques can powerfully and nondestructively probe materials while maintaining high resolution in both space and time. Unfortunately, these characterizations have been limited and sometimes even erroneous due to the difficulty of decoding the desired material information from features of the diffractograms. Currently, these features are identified non-comprehensively via human intuition, so the resulting models can only predict a subset of the available structural information. In the present work we show (i) how to compute machine-identified features that fully summarize a diffractogram and (ii) how to employ machine learning to reliably connect these features to an expanded set of structural statistics. To exemplify this framework, we assessed virtual electron diffractograms generated from atomistic simulations of irradiated copper. When based on machine-identified features rather than human-identified features, our machine-learning model not only predicted one-point statistics (i.e. density) but also a two-point statistic (i.e. spatial distribution) of the defect population. Hence, this work demonstrates that machine-learning models that input machine-identified features significantly advance the state of the art for accurately and robustly decoding diffractograms.

36 MATERIALS SCIENCE↗