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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Nonadditive CO 2 Uptake of Type II Porous Liquids Based on Imine Cages

Type II porous liquids can potentially exploit the fluidity of liquids and sorption properties of porous sorbents, yet CO 2 uptake in porous liquids is still poorly understood. Molecular simulations and experiments are used to examine CO 2 uptake by a prototypical porous liquid composed of porous organic cages (CC13) in 2′-hydroxyacetophenone (2′-HAP). The simulations are in reasonable agreement with experimental measurements of CO 2 solubility and provide unambiguous information on the partitioning of CO 2 within microenvironments in the liquid. Analysis of CO 2 dynamics is performed using these simulations, including assessing the self-diffusivity of CO 2 in both the neat solvent and porous liquid. This offers insights into the kinetics of CO 2 uptake and transport in type II porous liquids based on imine cages. Experiments with type II porous liquids formed by dissolving CC13 in three different size-excluded solvents show nonadditive CO 2 absorption relative to predictions based on ideal volume additivity. This nonadditive absorption behavior is also observed in simulations. Nonadditive CO 2 uptake is also demonstrated in type II porous liquids based on another imine-based porous cage, CC19.

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

Characterization of promoter elements of isoprene‐responsive genes and the ability of isoprene to bind START domain transcription factors

Abstract Isoprene has recently been proposed to be a signaling molecule that can enhance tolerance of both biotic and abiotic stress. Not all plants make isoprene, but all plants tested to date respond to isoprene. We hypothesized that isoprene interacts with existing signaling pathways rather than requiring novel mechanisms for its effect on plants. We analyzed the cis‐regulatory elements (CREs) in promoters of isoprene‐responsive genes and the corresponding transcription factors binding these promoter elements to obtain clues about the transcription factors and other proteins involved in isoprene signaling. Promoter regions of isoprene‐responsive genes were characterized using the Arabidopsis cis‐regulatory element database. CREs bind ARR1, Dof, DPBF, bHLH112, GATA factors, GT‐1, MYB, and WRKY transcription factors, and light‐responsive elements were overrepresented in promoters of isoprene‐responsive genes; CBF‐, HSF‐, WUS‐binding motifs were underrepresented. Transcription factors corresponding to CREs overrepresented in promoters of isoprene‐responsive genes were mainly those important for stress responses: drought‐, salt/osmotic‐, oxidative‐, herbivory/wounding and pathogen‐stress. More than half of the isoprene‐responsive genes contained at least one binding site for TFs of the class IV (homeodomain leucine zipper) HD‐ZIP family, such as GL2, ATML1, PDF2, HDG11, ATHB17. While the HD‐zipper‐loop‐zipper (ZLZ) domain binds to the L1 box of the promoter region, a special domain called the steroidogenic acute regulatory protein‐related lipid transfer, or START domain, can bind ligands such as fatty acids (e.g., linolenic and linoleic acid). We tested whether isoprene might bind in such a START domain. Molecular simulations and modeling to test interactions between isoprene and a class IV HD‐ZIP family START‐domain‐containing protein were carried out. Without membrane penetration by the HDG11 START domain, isoprene within the lipid bilayer was inaccessible to this domain, preventing protein interactions with membrane bound isoprene. The cross‐talk between isoprene‐mediated signaling and other growth regulator and stress signaling pathways, in terms of common CREs and transcription factors could enhance the stability of the isoprene emission trait when it evolves in a plant but so far it has not been possible to say what how isoprene is sensed to initiate signaling responses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Effects of Galactic Irradiation on Thermal and Electronic Transport in Tungsten

The impact of irradiation on the thermal and electronic properties of materials is a persistent puzzle, particularly defect formation at the atomic and nanoscales. This work examines the nanoscale effects of low-energy irradiation on tungsten (W), focusing on defect-induced modifications to thermal and electronic transport. Using the site-projected thermal conductivity (SPTC) method [A. Gautam et al. PSS-RRL, 2400306, 2024], bulk and twin-grain boundary W with vacancy defects based on the Norgett–Robinson–Torrens displacements per atom model is analyzed. SPTC provides a detailed prediction of postcascade spatial thermal conductivity distribution. Electronic conductivity activity is estimated using the “N 2 method” [K. Nepal et al. Carbon, 119711, 2025] to explore the consequences of vacancies and grain boundaries, highlighting the defect-dependent nature of charge transport behavior. These findings offer high-resolution insights into irradiation-driven transport phenomena, with implications for space-exposed materials and nanoscale thermal/electronic management.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Broad range material-to-system screening of metal–organic frameworks for hydrogen storage using machine learning

Hydrogen is pivotal in the transition to sustainable energy systems, playing major roles in power generation and industrial applications. Metal–organic frameworks (MOFs) have emerged as promising mediums for efficient hydrogen storage. However, identifying potential candidates for deployment is challenging due to the vast number of currently available synthesized MOFs. This study integrates molecular simulations, machine learning, and techno-economic analysis to evaluate the performance of MOFs across broad operation conditions for hydrogen storage applications. While previous screenings of MOF databases have predominantly emphasized high hydrogen capacities under cryogenic conditions, this study reveals that optimal temperatures and pressures for cost minimization depend on the raw price of the MOF. Specifically, when MOFs are priced at $15/kg, among the 9720 MOFs tested, 9692 MOFs achieve the lowest cost at temperatures between 170 K and 250 K and a pressure of 150 bar. Under these optimal conditions, 362 MOFs deliver a lower levelized cost of storage than 350 bar compressed gas hydrogen storage. Furthermore, this study reveals key material properties that result in low system cost, such as high surface areas (>3000 m2/g), large void fractions (>0.78), and large pore volumes (>1.1 cm3/g).

Hydrogen storage↗
DOE OSTI2024

3D printed lignin/polymer composite with enhanced mechanical and anti-thermal-aging performance

Lignin is the most abundant natural aromatic polymer globally but is still underutilized as a renewable material, even with its versatile properties. One approach to using lignin is to incorporate it in polymer composites, but this application is often limited by the poor mechanical performance of the resultant composite due to the poor interfacial adhesion. Following the structure–property relationship, stronger interactions were designed to be enhanced by alternating the functional groups of lignin. Here, this study applied a demethylation method to hardwood lignin (Ori-Lig) to introduce more phenolic hydroxyl groups. Research studies with rheology behavior and molecular simulations demonstrated that an increased phenolic hydroxyl content could improve the adhesion between the modified lignin (OH-Lig) and the polymer matrix at the interface. The tensile performance of the samples from the fused depositional modeling (FDM) 3D printing technique was also improved due to the improved interfacial adhesion. Specifically, by adding 10 wt% of the modified lignin, the tensile strength of the 3D printed samples could reach 46.1 MPa (40.2 MPa of Polyamide 12, PA12), and Young’s Modulus could be improved to 1.73 GPa (1.48 GPa of PA12). OH-Lig also improved anti-aging performance so that the tensile strength of the OH-Lig composite after thermo-aging (100 h at 140 °C) remained ~48MPa. The enhanced mechanical performance with anti-aging properties indicated that the OH-Lig composite may replace PA12 to reduce the use of petrol-based materials. This work showed the method of purposefully designing and modifying lignin structures to fulfill the interaction requirements between lignin and polymer matrix. The as-prepared lignin could be used to prepare functional composites to achieve improved mechanical performance and targeted functions at the same time.

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

Membrane-based solvent extraction for the recovery of rare earths from phosphate mining process streams

This study reports on the capture of rare earth elements (REEs) from phosphate industry process streams, including phosphoric acid (PA) sludge and phosphogypsum (PG), using a membrane solvent extraction (MSX) process. While MSX has been proven effective for a relatively concentrated feed, its effectiveness for dilute REEs solutions remains unexplored. Investigated PA-sludge and PG particles contain total REEs concentrations of ∼1100 and ∼320 ppm, respectively. Acid leaching, implemented to dissolve the REEs, significantly dilutes the REEs concentration to ∼210 ppm for PA-sludge leachate and ∼60 ppm for PG leachate. These low concentrations, compounded by the higher levels of non-REE ions and radioactive species, uranium (U) and thorium (Th), poses challenges to the MSX process. Here, we demonstrated that N,N,N′,N′-tetraoctyl-diglycolamide (TODGA) selectively binds REEs from a >3 M nitric-acid leachate while effectively rejecting U and Th. Concentrations of light REEs in strip solution were doubled compared to the feed, while heavy REEs were preferentially extracted. Furthermore, >99% purity gypsum, free of U and Th, was precipitated during the acid leaching process, aiding separation by removing significant amounts of non-REEs species (e.g., calcium) prior to the MSX process. Molecular simulations support the experimental data, suggesting preferential separation of heavy over light REEs. Based on these results, a cost-effective integrated process including pretreatment, acid leaching, MSX, and wastewater treatment is proposed for the co-recovery of REEs, phosphoric acid, gypsum, and U. This study shows MSX as a technically and economically feasible process for the recovery of REEs from low-concentration process streams, offering advantages over conventional solvent extraction.

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

A physics-informed operator regression framework for extracting data-driven continuum models

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

42 ENGINEERING↗
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 OSTI2024

Energetics of water expulsion from intervening space between two particles during aggregation

Solvent expulsion away from an intervening region between two approaching particles plays important roles in particle aggregation yet remains poorly understood. Here, in this work, we use metadynamics molecular simulations to study the free energy landscape of removing water molecules from gibbsite and pyrophyllite slit pores representing the confined spaces between two approaching particles. For gibbsite, removing water from the intervening region is both entropically and enthalpically unfavorable. The closer the particles approach each other, the harder it is to expel water molecules. For pyrophyllite, water expulsion is spontaneous, which is different from the gibbsite system. A smaller pore makes the water removal more favorable. When water is being drained from the intervening region, single chains of water molecules are observed in gibbsite pore, while in pyrophyllite pore water cluster is usually observed. Water-gibbsite hydrogen bonds help stabilize water chains, while water forms clusters in pyrophyllite pore to maximize the number of hydrogen bonds among themselves. This work provides the first assessment into the energetics and structure of water being drained from the intervening region between two approaching particles during oriented attachment and aggregation.

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

Environmental application of chlorine-doped graphitic carbon nitride: Continuous solar-driven photocatalytic production of hydrogen peroxide

Solar-driven photocatalytic generation of H 2 O 2 over metal-free catalysts is a sustainable approach for value-added chemical production. Here, we synthesized chlorine-doped graphitic carbon nitride (Cl-doped g-C 3 N 4 ) through a solvothermal method to effectively produce H 2 O 2 with a rate of 1.19 ± 0.06 µM min –1 under visible light irradiation, which was improved by 104 times compared to pristine g-C 3 N 4 . Continuous net production of H 2 O 2 was realized at a rate of 2.78 ± 0.10 µM min –1 up to 54 h with isopropanol as the hole scavenger, whereas H 2 O 2 production was only sustained for ~ 6 h without scavengers. Both molecular simulations and advanced spectroscopic characterizations elucidated that the Cl dopant increased the charge transfer rate, decreased the bandgap, and reduced the activation energy of the rate-limiting step of O 2 reduction, all of which favored H 2 O 2 production. Furthermore, this work implemented a novel metal-free photocatalyst for sustainable H 2 O 2 production and elucidated the mechanism for promoting H 2 O 2 production that can guide future photoreactive nanomaterial design.

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

Critical role of solvation on CC13 porous organic cages for design of porous liquids

Efficient carbon capture requires the design of new materials with high CO 2 selectivity and gas adsorption capacity that can be incorporated into existing industrial processes. Porous liquids (PLs) are promising candidate materials that consist of a nanoporous host and a solvent forming a liquid with permanent porosity based on exclusion of the solvent from the interior of the nanoporous host. Stable PLs are based on solvent-nanoporous host interactions, which can be evaluated through molecular simulations. Here, time- and temperature-dependent density functional theory simulations were performed between four solvents, 2-bromophenol, 4-methylphenol, 2,4-dimethylphenol, and cyclohexanone and the CC13 porous organic cage (POC) as a prototypical PL composition. Overall, minimal reactions occurred in the PL including no changes in the POC structure. Additionally, POC-solvent coordination occurred through interactions of neighboring functional groups such as methyl/bromide and hydroxyl on the solvent molecules with the POC surface. Therefore, the location rather than the number of functional groups on the solvent molecule controls the POC-solvent interactions. Additionally, the POC pore window contracted or expanded up to 8% during solvation, which correlates with the experimental solubility and static solvent-POC binding, where solvents that caused less contraction of the POC pore window increased POC solubility. Finally, these results allow for the design of optimized POC-based PL compositions based on solvent-nanoporous host binding and variation in the pore window during solvation.

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

Ultra-efficient polymer binder for silicon anode in high-capacity lithium-ion batteries

As a highly promising anode material for high-capacity lithium-ion batteries (LIBs), the low electronic conductivity and large volume variation of silicon (Si) make the slurry-coating Si based electrode requiring high content of “inert” materials and suffering rapid capacity fading. In this work, a polyimine, synthesized via one-step condensation reaction, has been demonstrated as an ultra-efficient polymer binder that can resolve the above issues. The polyimine binder containing Si electrode delivers superior electrochemical performance: a delithiation specific capacity of 804.4 mAh g –1 with capacity retention of 82.4% after 1000 cycles at the current density of 2 A g –1 . The high efficiency of polyimine binder for Si electrode has also been demonstrated with ultrahigh weight ratio of “active” material to “inert” material (R A/I ). The electrode with 95 wt% of Si (95Si/Polyimine, R A/I = 19) reveals a reversible delithiation capacity of 2114 mAh g –1 (capacity retention ~ 80.4%) over 200 cycles at the current density of 400 mA g –1 . Even at the high current density of 2 A g –1 , a delithiation capacity of 1087.8 mAh g –1 after 500 cycles can be obtained. Molecular simulations and atomic force microscopy (AFM) indentation are utilized to investigate the ultra-efficiency of polyimine binder. With simple manufacturing process and ultra-efficient binder performance, the designed polyimine binder will be definitely meaningful in achieving low-cost and high-capacity LIBs with prolonged cycle life.

25 ENERGY STORAGE↗
DOE OSTI2022

Tuning Enzyme Thermostability via Computationally Guided Covalent Stapling and Structural Basis of Enhanced Stabilization

Enhancing the thermostability of enzymes without impacting their catalytic function represents an important yet challenging goal in protein engineering and biocatalysis. We recently introduced a novel method for enzyme thermostabilization that relies on the computationally guided installation of genetically encoded thioether “staples” into a protein via cysteine alkylation with the noncanonical amino acid O-2-bromoethyl tyrosine (O2beY). Here, we demonstrate the functionality of an expanded set of electrophilic amino acids featuring chloroacetamido, acrylamido, and vinylsulfonamido side-chain groups for protein stapling using this strategy. Using a myoglobin-based cyclopropanase as a model enzyme, our studies show that covalent stapling with pchloroacetamido-phenylalanine (pCaaF) provides higher stapling efficiency and enhanced stability (thermodynamic and kinetic) compared to the other stapled variants and the parent protein. Interestingly, molecular simulations of conformational flexibility of the cross-links show that the pCaaF staple allows fewer energetically feasible conformers than the other staples, and this property may be a broader indicator of stability enhancement. Using this strategy, pCaaF-stapled variants with significantly enhanced stability against thermal denaturation (ΔT m ' = +27 °C) and temperature-induced heme loss (ΔT 50 = +30 °C) were obtained while maintaining high levels of catalytic activity and stereoselectivity. Crystallographic analyses of singly and doubly stapled variants provide key insights into the structural basis for stabilization, which includes both direct interactions of the staples with protein residues and indirect interactions through adjacent residues involved in heme binding. This work expands the toolbox of protein stapling strategies available for protein stabilization.

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

Multistage Condensation Pathway Minimizes Hysteresis in Water Harvesting with Large-Pore Metal–Organic Frameworks

Metal-organic frameworks (MOFs) have emerged as promising materials for atmospheric water harvesting (AWH). Large-pore MOFs provide high water capacity, but their significant hysteresis between sorption and desorption makes them unsuitable for AWH. Co2Cl2(BTDD) is a noteworthy exception. This MOF has large, 2.2 nm diameter one-dimensional pores, and combines both record-high water capacity and minimal hysteresis, making it an excellent material for water capture in arid areas. Sorption reversibility in Co2Cl2(BTDD) has been attributed to continuous water uptake. However, the sharp adsorption/desorption in the isotherms supports a discontinuous first-order transition. Here we use molecular simulations to compute the water adsorption and desorption pathways and isotherms in a Co2Cl2(BTDD) model, to elucidate how does this MOF achieve reversibility despite its large pore size. The simulations reveal a multi-stage mechanism of discontinuous water uptake facilitated by spatial segregation of rows of hydrophilic metal sites bridged by ~1 nm hydrophobic ligands. The multi-stage mechanism breaks the barrier of capillary condensation into smaller, easier to surmount ones, resulting in a facile process despite the sharp density discontinuity between confined liquid and vapor. Our results explain why exchanging Co2+ for Ni2+ or Clfor F- in the MOF have minimal impact on the condensation and desorption pressures. On the other hand, we predict that a decrease in hydrophilicity of the MOF vertices would strongly increase the hysteresis. We expect that the relationships between spatial distribution of hydrophilic sites and hysteresis unraveled in this study assist the design of water harvesting materials with maximal capacity and reversibility.

Zaragoza, Alberto↗
DOE OSTI2023

Identifying High-Performance Metal–Organic Frameworks for Low-Temperature Oxygen Recovery from Helium by Computational Screening

Metal–organic frameworks (MOFs) are an important class of porous crystalline materials for applications ranging from gas adsorption and separation to catalysis. There are thousands of potential MOFs available for separation applications. Here, we developed a computational approach to screen MOFs for the separation of oxygen–helium mixtures at low temperatures (100–200 K), conditions that were motivated by issues associated with propulsion in space-based settings. We used detailed molecular simulations for a small number of MOFs to develop screening methods that were then used to estimate the optimum temperatures for separations using pressure swing adsorption for 2932 MOFs from the CoRE MOF database and the swing capacity and oxygen–helium selectivity at these temperatures. We used the stability of the best-performing structures in the presence of moisture as a means to provide a short list of high-performance materials. In addition to identifying specific materials for oxygen–helium separations, this approach could prove useful for selecting adsorbents for other gas separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗