Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probing degradation at solid-state battery interfaces using machine-learning interatomic potential

Solid-state batteries featuring fast ion-conducting solid electrolytes are promising next-generation energy storage technologies, yet challenges remain for practical deployment due to electro-chemo-mechanical instabilities at solid-solid interfaces. These interfaces, which include homogeneous/internal interfaces such as grain boundaries (GBs) and heterogeneous/external interfaces between solid-electrolyte and electrode materials, can impede Li-ion transport, deteriorate performance, and eventually lead to cell failure. Here, in this study, we leverage large-scale molecular simulations, enabled by validated machine-learning interatomic potentials, to directly probe the onset of interfacial degradation at the garnet Li 7 La 3 Zr 2 O 12 (LLZO) solid-electrolyte/LiCoO 2 (LCO) cathode interface. By surveying different interfacial geometries and compositions, it is found that Li-deficient interfaces can lead to severe interfacial disordering with cation mixing and Co interdiffusion from LCO into LLZO. By contrast, Li-sufficient interfaces are less disordered, although elemental segregation with local ordering is observed. As a consequence of Co interdiffusion, Co-rich regions are formed at the GBs of LLZO due to cation segregation and trapping effects. This behavior is independent of the GB tilting axis, degree of disorder at the GBs, and Co concentration, which implies Co clustering at GBs is a general phenomenon in polycrystalline LLZO and can dictate its overall transport and mechanical properties. Our findings elucidate the underlying fundamental mechanisms that give rise to experimentally observed physicochemical properties and provide guidelines for interface design that can mitigate interfacial degradation and improve cycling performance.

25 ENERGY STORAGE↗
DOE OSTI2025

Hydrogen diffusion induced dislocation transformations in a nickel superalloy

The diffusion of hydrogen in metals and alloys induces embrittlement that can adversely affect the structural properties. We examine the adsorption and diffusion of hydrogen in Inconel-718 (IN-718), and scrutinize the ensuing effects on the dislocation behavior in the alloy to elucidate the fundamental mechanisms of hydrogen-microstructure interactions from classical molecular simulations. Hydrogen adsorption increases with time until the surface saturates, while hydrogen diffusion exhibits strong temperature dependence, with diffusion coefficients converging above 1300 K regardless of the initial hydrogen concentration in the alloy. The diffusion in IN-718 is significantly sluggish than in pure Ni, Fe, or Cr, and is strongly impacted by hydrogen concentrations, resulting in an order of magnitude higher diffusion coefficient for hydrogen (10-14 m2/s relative to 10-15 m2/s) at high concentrations, especially below 600 K. Hydrogen diffusion coefficient varies from 10-12 to 10-15 m2/s in IN-718 depending on temperature (500–1400 K). More critically, our results reveal that increasing hydrogen concentration induces microstructural changes in the alloy, transforming perfect dislocations into stair-rods and Shockley partials, with higher temperatures favoring the latter. The results are significant for hydrogen fuel applications to gain insights into the materials chemistry for designing safer and more efficient propulsion systems, particularly in high-performance environments related to controlled hydrogen combustion applications.

Banerjee, Tanumoy↗
DOE OSTI2025

LuGo: An enhanced quantum phase estimation implementation

Quantum Phase Estimation (QPE) is a cardinal algorithm in quantum computing that plays a crucial role in various applications, including cryptography, molecular simulation, and solving systems of linear equations. However, the standard implementation of QPE faces challenges related to time complexity and circuit depth, which limit its practicality for large-scale computations. We introduce LuGo, a novel framework designed to enhance the performance of QPE by reducing circuit duplication, as well as using parallelization techniques to achieve faster generation of the QPE circuit and gate reduction. We validate the effectiveness of our framework by generating quantum linear solver circuits, which require both QPE and inverse QPE, to solve linear systems of equations. LuGo achieves significant improvements in both computational efficiency and hardware requirements without compromising on accuracy. Compared to a standard QPE implementation, LuGo reduces time consumption to generate a circuit that solves a 2 6 × 2 6 system matrix by a factor of 50.68 and over 31× reduction of quantum gates and circuit depth, with no fidelity loss on an ideal quantum simulator. Furthermore, we demonstrated the versatility and scalability of LuGo enabled HHL algorithm by simulating a canonical Hele-Shaw fluid problem using a quantum simulator. With these advantages, LuGo paves the way for more efficient implementations of QPE, enabling broader applications across several quantum computing domains.

Quantum algorithm↗
DOE OSTI2025

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Distinguishing Surface and Bulk Reactivity: Concentration-Dependent Kinetics of Iodide Oxidation by Ozone in Microdroplets

Iodine oxidation reactions play an important role in environmental, biological, and industrial contexts. The multiphase reaction between aqueous iodide and ozone is of particular interest due to its prevalence in the marine atmosphere and unique reactivity at the air–water interface. Here, we explore the concentration dependence of the I – + O 3 reaction in levitated microdroplets under both acidic and basic conditions. To interpret the experimental kinetics, molecular simulations are used to benchmark a kinetic model, which enables insight into the reactivity of the interface, the nanometer-scale subsurface region, and the bulk interior of the droplet. For all experiments, a kinetic description of gas- and liquid-phase diffusion is critical to interpreting the results. We find that the surface dominates the iodide oxidation kinetics under concentrated and acidic conditions, with the reactive uptake coefficient approaching an upper limit of 10 –2 at pH 3. In contrast, reactions in the subsurface dominate under more dilute and alkaline conditions, with inhibition of the surface reaction at pH 12 and an uptake coefficient that is 10× smaller. In conclusion, the origin of a changing surface mechanism with pH is explored and compared to previous ozone-dependent measurements.

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

Interfacial Ice Density Fluctuations Inform Surface Ice-Philicity

The propensity of a surface to nucleate ice or bind to ice is governed by its ice-philicity─its relative preference for ice over liquid water. However, the relationship between the features of a surface and its ice-philicity is not well understood, and for surfaces with chemical or topographical heterogeneity, such as proteins, their ice-philicity is not even well-defined. In the analogous problem of surface hydrophobicity, it has been shown that hydrophobic surfaces display enhanced low water-density (vapor-like) fluctuations in their vicinity. To interrogate whether enhanced ice-like fluctuations are similarly observed near ice-philic surfaces, here we use molecular simulations and enhanced sampling techniques. Using a family of model surfaces for which the wetting coefficient, k , has previously been characterized, we show that the free energy of observing rare interfacial ice-density fluctuations decreases monotonically with increasing k . By utilizing this connection, we investigate a set of fcc systems and find that the (110) surface is more ice-philic than the (111) or (100) surfaces. By additionally analyzing the structure of interfacial ice, we find that all surfaces prefer to bind to the basal plane of ice, and the topographical complementarity of the (110) surface to the basal plane explains its higher ice-philicity. Using enhanced interfacial ice-like fluctuations as a measure of surface ice-philicity, we then characterize the ice-philicity of chemically heterogeneous and topologically complex systems. In particular, we study the spruce budworm antifreeze protein (sbwAFP), which binds to ice using a known ice-binding site (IBS) and resists engulfment using nonbinding sites of the protein (NBSs). We find that the IBS displays enhanced interfacial ice-density fluctuations and is therefore more ice-philic than the two NBSs studied. We also find the two NBSs are similarly ice-phobic. By establishing a connection between interfacial ice-like fluctuations and surface ice-philicity, our findings thus provide a way to characterize the ice-philicity of heterogeneous surfaces.

crystals↗
DOE OSTI2024

Degradation of Anion Exchange Membranes by Cation Elimination: Impact on Water Uptake, Nanostructure, and Ionic Mobility

Anion exchange membranes (AEMs) are an attractive platform for fuel cell and electrolysis technologies since they enable the use of cheaper, nonplatinum group metal electrodes and components. However, the widespread adoption of AEM-based devices is limited by chemical degradation of the AEM in the highly alkaline medium of operation. Experimental studies report three major pathways of degradation, including substitution (S N 2), elimination (E2), and backbone cleavage. The decline in membrane performance is likely due to a cumulative effect of several pathways, which has proven to be difficult to disentangle through experiments. Here we use coarse-grained molecular simulations to isolate the impact of E2 degradation, where cationic sites are replaced by an alkene group, on AEM water uptake, nanoscale morphology, structure of the ion-conducting channels, water mobility, and ionic conductivity. Our studies focus on the well-studied model AEM polyphenylene oxide with tetraalkylammonium cationic sites (PPO-TMA). We find that about half of the hydrophobic groups resulting from the degradation retract into the hydrophobic domains of the membrane. The reduction in ion exchange capacity (IEC) and the hydrophobicity of the alkene groups decrease both the equilibrium water uptake (%WU) and the number of water molecules per cation (λ) of the membrane. Interestingly, the evolution of λ with the IEC seems to follow that of undegraded PPO-TMA membranes. Additionally, we ascribe this to the similarity between the structure of the degraded monomer and the neutral one in the polymer. We conclude that most of the performance losses originate in the lower hydration of the degraded membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting Partial Atomic Charges in Metal–Organic Frameworks: An Extension to Ionic MOFs

Molecular simulation is an invaluable tool to predict and understand the usage of metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges, commonly obtained from density functional theory (DFT) calculations, are often required to model the electrostatic interactions between the MOF and adsorbates, especially when the adsorbates have dipole or quadrupole moments, such as water and CO 2 . Machine learning (ML) models have been previously employed to predict partial charges and avoid the computational cost associated with DFT calculations. However, previous ML models suffer from small training data sets, which limit their scope of application. In this work, we introduce two novel machine learning models, PACMOF2-neutral and PACMOF2-ionic, aimed at predicting the density-derived electrostatic and chemical (DDEC6) partial atomic charges for both neutral and ionic MOFs. These models not only yield DFT-level accuracy at a fraction of the computational cost but also demonstrate a remarkable improvement in prediction of adsorption, as validated with grand canonical Monte Carlo simulations. Furthermore, the robustness and fast computational time of the PACMOF2 models, along with their transferability to other porous materials such as covalent organic frameworks and zeolites, underscores their potential in high-throughput screening of MOFs for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nuclear Quantum Effects on the Electronic Structure of Water and Ice

The electronic properties and optical response of ice and water are intricately shaped by their molecular structure, including the quantum mechanical nature of the hydrogen atoms. Despite numerous previous studies, a comprehensive understanding of the nuclear quantum effects (NQEs) on the electronic structure of water and ice at finite temperatures remains elusive. Here, we utilize molecular simulations that harness efficient machine-learning potentials and many-body perturbation theory to assess how NQEs impact the electronic bands of water and hexagonal ice. By comparing path-integral and classical simulations, we find that NQEs lead to a larger renormalization of the fundamental gap of ice, compared to that of water, ultimately yielding similar bandgaps in the two systems, consistent with experimental estimates. Our calculations suggest that the increased quantum mechanical delocalization of protons in ice, relative to water, is a key factor leading to the enhancement of NQEs on the electronic structure of ice.

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

Elastomer Mechanics of Cross-Linked Linear-Ring Polymer Blends

Cross-linking a blend of linear and ring polymers creates a new topology-based dual-network elastomer in which the two components differ significantly in their topology. We use molecular simulations and topological analysis to examine key mechanical properties as functions of ring polymer volume fraction Φ R . For Φ R < Φ R *, where the rings begin to overlap, the network shear modulus G and the maximum stretch ratio λ p are weakly dependent on Φ R . For Φ R > Φ R *, entanglements trapped in the network are diluted as the rings overlap, leading to a significant decrease in G and an increase in λ p with increasing Φ R . Here, the peak tensile stress, σ p , exhibits a maximum around Φ R *, indicating an enhancement of network strength due to the stronger cohesion from the entanglements between linear and ring polymers.

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

The Effects of Morphology and Hydration on Anion Transport in Self-Assembled Nanoporous Membranes

Ordered nanoporous polymer membranes offer opportunities for systematically probing the mechanisms of ion transport under confinement and for realizing useful materials for electrochemical devices. Here, we examine the impact of morphology and ion hydration on the transport of hydroxide and bromide anions in nanostructured polymer membranes with 1 nm scale pores. We use aqueous lyotropic self-assembly of an amphiphilic monomer, with a polymerizable surfactant to create direct hexagonal (H I ) and gyroid mesophases. UV-induced cross-linking leads to the formation of nanoporous polymers with water continuous channels. The membranes are mechanically robust and chemically durable, resisting degradation during extended exposure to 1 M NaOH solutions. We use a combination of electrochemical impedance spectroscopy, pulsed-field gradient NMR spectroscopy, and molecular simulations to elucidate anion and water transport. The as-prepared hexagonal systems display higher conductivity and lower activation energies for both anions relative to the gyroid system. When compared at equivalent hydration, however, gyroid and hexagonal membranes show similar activation energies, with nearly identical conductivities at ambient temperatures. Both ionic conductivity and water diffusivity increase with increasing hydration. Here, the water uptake as a function of relative humidity for the hexagonal and gyroid mesophases ultimately dictates the water diffusion and magnitude of the ionic conductivity, with the hexagonal system showing overall higher capacity for hydration and thus faster ion transport. The durability of these materials under aggressive alkaline conditions and their relatively high hydroxide ion conductivity suggest that these nanostructured polymers could be of interest as membranes for alkaline fuel cells.

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

Molecular Insights into Geochemical Reactions of Iron-Bearing Minerals: Implications for Hydrogen Geo-Storage

Here, this study investigates the reaction of hydrogen (H 2 ) with pyrite (FeS 2 ), focusing on how temperature and the presence of water influence the reaction pathways and kinetics. Utilizing computational molecular simulations and kinetic analyses, we explore the impact of these factors on the formation of hydrogen sulfide (H 2 S) and related species. First, grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations reveal that physical H 2 adsorption occurs in distinct layers on the pyrite surface. In addition, increased temperatures reduce the absolute adsorption capacities. Reactive MD simulations demonstrate that H 2 interacts differently with pyrite under varying conditions. At 298 K, H 2 reacts with pyrite to form HS – , leading to the formation of HS – through covalent bonding with sulfur of pyrite. However, no H 2 S is produced at this temperature, suggesting that a kinetic barrier (i.e., activation energy) may prevent this reaction. At higher temperatures, H 2 S production significantly increases. The presence of water introduces additional complexity to the reaction mechanism. Unlike dry conditions, water enhances H 2 S generation, even at low temperatures. Water also facilitates the formation of additional products, such as SOH, indicating a more intricate chemical environment on the pyrite surface. Our findings identify the association of HS – ions to form H 2 S as the rate-limiting step, with temperature influencing this process. This finding suggests that while the presence of water can create a more dynamic reaction environment, the overall mechanisms leading to H 2 S formation remain consistent. These outcomes suggest the need for developing targeted strategies to manage and control H 2 S emissions within the context of underground hydrogen storage.

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

Interzeolite Transformation through Cross-Nucleation: A Molecular Mechanism for Seed-Assisted Synthesis

Polymorph selection and efficient crystallization are central goals in zeolite synthesis. Crystalline seeds are used for both purposes. While it has been proposed that zeolite seeds induce interzeolite transformation by dissolving into structural units that promote nucleation of the daughter crystal, the seed’s structural elements do not always match those of the target zeolite. This discrepancy raises the question of how the seed promotes the daughter phase. Here, we present the first molecularly resolved investigation of seed-assisted zeolite synthesis. Using molecular simulations, we reproduce the experimental finding that a parent zeolite can promote the nucleation of a daughter zeolite even when it lacks common composite building units (CBUs) or crystal planes. Modeling the seed-assisted synthesis of an AFI-type zeolite using zeolite CHA, our simulations indicate that stand-alone CBUs from the parent seed do not facilitate daughter crystal formation. However, introducing the intact seed significantly reduces the synthesis time, supporting that seed integrity is key to increased efficiency. This reduction arises from the cross-nucleation of the AFI-type zeolite on the CHA (001) face. We find that parent and daughter zeolites are connected by an interfacial transition layer with an order distinct from that of both zeolites. Simulations reveal that cross-nucleation occurs over a broad range of synthesis conditions. We argue that cross-nucleation would be most favorable for zeolite pairs that share crystalline planes such as those forming intergrowths. In conclusion, our findings suggest that the prevalence of intergrowths with a common lattice plane in zeolite synthesis is likely a kinetic effect of accelerated cross-nucleation.

Crystallization↗
DOE OSTI2026

Mechanistic Insights into Defect-Mediated Crystallization Revealed by Lattice Strain Evolution

Structural defects and lattice strain are intrinsic to many crystalline materials, yet their roles in controlling chemical reaction mechanisms and directing crystallization pathways remain poorly understood. Here, in this study, we revealed the three-dimensional evolution of strain and dislocation defects at the nanoscale during the growth of heterogeneously nucleated barite (BaSO 4 ) and calcite (CaCO 3 ) crystals by using coherent X-ray scattering, electron microscopy, and molecular simulations. Unlike barite, which formed with minimal internal strain, calcite developed dislocation defects and exhibited spatially varying strain that increased during growth. During growth in Sr-rich solutions, calcite likely incorporates Sr 2+ into the defects, which further modulates the local lattice structure and increases both the compressive and tensile strain. These findings suggest that calcite crystallization was likely dominated by attachment of precursor phases, which gave rise to defect-enriched domain structures not predicted by classical growth models. By linking defect formation to ion incorporation and growth dynamics, this work provides fundamental insight into how lattice-level strain heterogeneity governs the chemical reactivity of ionic crystals.

Bragg coherent diffractive imaging↗
DOE OSTI2025

Molecular model of TFIIH recruitment to the transcription-coupled repair machinery

Transcription-coupled repair (TCR) is a vital nucleotide excision repair sub-pathway that removes DNA lesions from actively transcribed DNA strands. Binding of CSB to lesion-stalled RNA Polymerase II (Pol II) initiates TCR by triggering the recruitment of downstream repair factors. Yet it remains unknown how transcription factor IIH (TFIIH) is recruited to the intact TCR complex. Combining existing structural data with AlphaFold predictions, we build an integrative model of the initial TFIIH-bound TCR complex. We show how TFIIH can be first recruited in an open repair-inhibited conformation, which requires subsequent CAK module removal and conformational closure to process damaged DNA. In our model, CSB, CSA, UVSSA, elongation factor 1 (ELOF1), and specific Pol II and UVSSA-bound ubiquitin moieties come together to provide interaction interfaces needed for TFIIH recruitment. STK19 acts as a linchpin of the assembly, orienting the incoming TFIIH and bridging Pol II to core TCR factors and DNA. Molecular simulations of the TCR-associated CRL4CSA ubiquitin ligase complex unveil the interplay of segmental DDB1 flexibility, continuous Cullin4A flexibility, and the key role of ELOF1 for Pol II ubiquitination that enables TCR. Collectively, these findings elucidate the coordinated assembly of repair proteins in early TCR.

Paul, Tanmoy↗
DOE OSTI2025

Metastability and Ostwald step rule in the crystallisation of diamond and graphite from molten carbon

Abstract Experimental challenges in determining the phase diagram of carbon at temperatures and pressures near the graphite-diamond-liquid triple point are often related to the persistence of metastable crystalline or glassy phases, superheated crystals, or supercooled liquids. A deeper understanding of the crystallisation kinetics of diamond and graphite is crucial for effectively interpreting the outcomes of these experiments. Here, we reveal the microscopic mechanisms of diamond and graphite nucleation from liquid carbon through molecular simulations with first-principles machine learning potentials. Our simulations accurately reproduce the experimental phase diagram of carbon near the triple point and show that liquid carbon crystallises spontaneously upon cooling. Metastable graphite crystallises in the domain of diamond thermodynamic stability at pressures above the triple point. Furthermore, whereas diamond crystallises through a classical nucleation pathway, graphite follows a two-step process in which low-density fluctuations forego ordering. Calculations of the nucleation rates of the two competing phases confirm this result and reveal a manifestation of Ostwald’s step rule, where the strong metastability of graphite hinders the transformation to the stable diamond phase. Our results provide a key to interpreting melting and recrystallisation experiments and shed light on nucleation kinetics in polymorphic materials with deep metastable states.

Science & Technology - Other Topics↗
DOE OSTI2025

Phosphorylation toggles the SARS-CoV-2 nucleocapsid protein between two membrane-associated condensate states

Abstract The Nucleocapsid protein (N) of SARS-CoV-2 plays a critical role in the viral lifecycle by regulating RNA replication and by packaging the viral genome. N and RNA phase separate to form condensates that may be important for these functions. Both functions occur at membrane surfaces, but how N toggles between these two membrane-associated functional states is unclear. Here, we reveal that phosphorylation switches how N condensates interact with membranes, in part by modulating condensate material properties. Our studies also show that phosphorylation alters N’s interaction with viral membrane proteins. We gain mechanistic insight through structural analysis and molecular simulations, which suggest phosphorylation induces a conformational change in N that softens condensate material properties. Together, our findings identify membrane association as a key feature of N condensates and provide mechanistic insights into the regulatory role of phosphorylation. Understanding this mechanism suggests potential therapeutic targets for COVID infection.

Science & Technology - Other Topics↗
DOE OSTI2025

Translocation mechanism of xeroderma pigmentosum group D protein on single-stranded DNA and genetic disease etiology

Abstract XPD is a key nucleotide excision repair (NER) protein whose function is vital for genome integrity. During NER, XPD serves as a 5′−3′ single-strand DNA translocase that enables lesion scanning and verification in genomic DNA. Yet, its translocation mechanism is incompletely understood. Here we use molecular simulations and chain-of-replicas path optimization methods to model the ATP-driven translocation mechanisms of XPD and its bacterial homolog DinG, revealing all on-path metastable intermediates and corresponding kinetic rates. We identify the XPD(DinG) global domain motions that modulate the strength of DNA association at the opposing ends of the DNA-binding groove. During the ATP hydrolysis cycle, alternating weak and strong interactions at two defined groove constrictions enable DNA reptation and forward displacement of the ATPase. Moreover, we show that DNA- or ATP-binding residues directly involved in translocation are hotspots for genetic disease mutations. Thus, our findings shed light on the etiology of XPD-associated genetic syndromes.

Paul, Tanmoy↗