Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Fast and Accurate Machine Learning Strategy for Calculating Partial Atomic Charges in Metal–Organic Frameworks

Computational high-throughput screening using molecular simulations is a powerful tool for identifying top-performing metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges are often required to model the electrostatic interactions between the MOF and the adsorbate, especially when the adsorption involves molecules with dipole or quadrupole moments such as water and CO 2 . Although ab initio methods can be used to calculate accurate partial atomic charges, these methods are impractical for screening large material databases because of the high computational cost. We developed a random forest machine learning model to predict the partial atomic charges in MOFs using a small yet meaningful set of features that represent both the elemental properties and the local environment of each atom. The model was trained and tested on a collection of about 320 000 density-derived electrostatic and chemical (DDEC) atomic charges calculated on a subset of the Computation-Ready Experimental Metal–Organic Framework (CoRE MOF-2019) database and separately on charge model 5 (CM5) charges. The model predicts accurate atomic charges for MOFs at a fraction of the computational cost of periodic density functional theory (DFT) and is found to be transferable to other porous molecular crystals and zeolites. In conclusion, a strong correlation is observed between the partial atomic charge and the average electronegativity difference between the central atom and its bonded neighbors.

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

Quantification of the REE 3+ aqua ions and chloride species in aqueous fluids by in situ Raman spectroscopy using perturbations of the water band

Acidic NaCl-rich aqueous fluids play a crucial role in forming hydrothermal rare earth elements (REE) mineral deposits. Aqueous REE mobility is mostly controlled by the stabilities of REE 3+ and REE chloride species. Our current knowledge of REE speciation is based on solubility data, thermodynamic models and in situ spectroscopic measurements, sometimes coupled with molecular simulations. Here, in this study, we investigate Nd and Yb speciation in pH2 Cl-bearing solutions at 25 °C and 0.1 MPa with variable Cl/REE ratios using Raman Spectroscopy in solutions with 0.1 to 0.6 mol/kg NdCl 3 or YbCl 3 and 0.2 to 3.2 mol/kg NaCl. Due to the challenges in resolving the REE-Cl band, we developed a new method using the water vibrational mode and multivariate curve resolution (MCR) analysis. The Raman spectra for the vibrational band of water (2700 to 3900 cm –1 ) were collected at 25 °C and fitted by three Gaussian sub-peaks, then quantified using MCR analysis to de-convolute the water band into bulk H 2 O and the perturbations caused by of Cl – , REE 3+ , and REE chloride species. REE speciation based on the perturbations of the water band indicates that REE 3+ aqua ions dominate acidic solutions at 25 °C, but up to ~20 mol% YbCl 2+ forms at high YbCl 3 concentrations. The new method is promising for quantifying in situ speciation of the REE 3+ aqua ions and REE chloride species in aqueous fluids while providing information on the hydration of ions. This method improves our molecular level understanding of REE aqueous species stability and their role in REE mobilization during fluid-rock interaction.

58 GEOSCIENCES↗
DOE OSTI2024

Role of Dynamic Polarization Interactions in the Electrical Double Layer at Calcite (104) Interfaces with Aqueous Solutions

Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.

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

Protein–Protein Complex Stability Controls Substrate Scope in a β‐Ketoacyl‐ACP Reductase Specific for Medium Chains

Assembly-line enzymes carry out multistep synthesis of important metabolites by using acyl carrier proteins (ACPs) to shuttle intermediates along defined sequences of active sites. Despite longstanding interest in reprogramming these systems for metabolic engineering and biosynthetic chemistry, the mechanisms underlying their reaction order remain poorly understood and difficult to control. In this work, we describe a β-ketoacyl-ACP reductase from Pseudomonas putida (PpFabG4) with an unusual selectivity for medium chains and use it to explore the molecular basis of substrate specificity in enzymes that pull intermediates from fatty acid synthesis, a common route to specialized products. X-ray crystallography shows no obvious barriers to short-chain binding. Molecular simulations and supporting mutational analyses indicate that substrate preference arises instead from a weak enzyme–ACP interaction that is stabilized by medium acyl chains but not by short chains. Indeed, mutations that strengthen this interaction for PpFabG4 or weaken it for EcFabG, an Escherichia coli β-ketoacyl-ACP reductase with a broad substrate specificity, can enhance or reduce activity on short-chain substrates by over 100-fold. Our findings show how the stability of enzyme-ACP interactions can control substrate scope in promiscuous enzymes and guide the exchange of intermediates between (and within) assembly-line systems.

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

DeePKS + ABACUS as a Bridge between Expensive Quantum Mechanical Models and Machine Learning Potentials

Recently, the development of machine learning (ML) potentials has made it possible to perform large-scale and long-time molecular simulations with the accuracy of quantum mechanical (QM) models. However, for different levels of QM methods, such as density functional theory (DFT) at the meta-GGA level and/or with exact exchange, quantum Monte Carlo, etc., generating a sufficient amount of data for training an ML potential has remained computationally challenging due to their high cost. In this work, we demonstrate that this issue can be largely alleviated with Deep Kohn–Sham (DeePKS), an ML-based DFT model. DeePKS employs a computationally efficient neural network-based functional model to construct a correction term added upon a cheap DFT model. Upon training, DeePKS offers closely matched energies and forces compared with high-level QM method, but the number of training data required is orders of magnitude less than that required for training a reliable ML potential. As such, DeePKS can serve as a bridge between expensive QM models and ML potentials: one can generate a decent amount of high-accuracy QM data to train a DeePKS model and then use the DeePKS model to label a much larger amount of configurations to train an ML potential. Further, this scheme for periodic systems is implemented in a DFT package ABACUS, which is open source and ready for use in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Sub-5 Ångstrom Porosity Tuning in Calixarene-Derived Porous Liquids via Supramolecular Complexation Construction

Sub-Ångstrom-level porosity engineering, which is appealing in gas separations, has been demonstrated in solid carbon, polymer, and framework materials but rarely achieved in the liquid phase. In this work, a gas molecular sieving effect in the liquid phase at sub-5 Ångstrom scale is created via sophisticated porosity tuning in calixarene-derived porous liquids (PLs). Type II PLs are constructed via supramolecular complexation between the sodium salts of calixarene derivatives and crown ether solvents. The chemical structure variation and assembly behavior of the porous host upon PL construction are monitored by spectroscopy-, X-ray-, and neutron-scattering techniques. The presence of permanent porosity in calixarene-derived PLs is verified by pressure swing gas uptake, altered CO 2 physisorption behavior, and molecular simulations. Sub-5 Ångstrom porosity tuning within the PL phase is achieved by introducing bulky substituted groups on the benzene ring of the calixarene host, which then greatly affects the dynamic motion and transport behavior of CO 2 molecules and the Xe uptake performance. Further, the approach being demonstrated in this work represents a promising pathway to tune and leverage the porosity effect for enhanced gas uptake capacity and selectivity in liquid sorbents.

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

Novel Zwitterionic Polyurethane-in-Salt Electrolytes with High Ion Conductivity, Elasticity, and Adhesion for High-Performance Solid-State Lithium Metal Batteries

This study presents a novel polymer-in-salt (PIS) zwitterionic polyurethane-based solid polymer electrolyte (zPU-SPE) that offers high ionic conductivity, strong interaction with electrodes, and excellent mechanical and electrochemical stabilities, making it promising for high-performance all solid-state lithium batteries (ASSLBs). The zPU-SPE exhibits remarkable lithium-ion (Li+) conductivity (3.7 × 10⁻⁴ S cm−1 at 25 °C), enabled by exceptionally high salt loading of up to 90 wt.% (12.6 molar ratio of Li salt to polymer unit) without phase separation. It addresses the limitations of conventional SPEs by combining high ionic conductivity with a Li+ transference number of 0.44, achieved through the incorporation of zwitterionic groups that enhance ion dissociation and transport. The high surface energy (338.4 J m−2) and elasticity ensure excellent adhesion to Li anodes, reducing interfacial resistance and ensuring uniform Li+ flux. When tested in Li||zPU||LiFePO₄ and Li||zPU||S/C cells, the zPU-SPE demonstrated remarkable cycling stability, retaining 76% capacity after 2000 cycles with the LiFePO4 cathode, and achieving 84% capacity retention after 300 cycles with the S/C cathode. Molecular simulations and a range of experimental characterizations confirm the superior structural organization of the zPU matrix, contributing to its outstanding electrochemical performance. The findings strongly suggest that zPU-SPE is a promising candidate for next-generation ASSLBs.

Wang, Kun↗
DOE OSTI2025

Computational screening of fly ash zeolite sorbents for boric acid removal

In the United States, many impoundments at coal-fired power plants contain elevated contaminants like arsenic, boron, barium, and selenium. Zeolites synthesized from fly ash show promise as sorbents for these contaminants. However, optimizing sorption capacity is challenging due to numerous possible topologies, silicon to aluminum (Si/Al) ratios, and cation types. In this study, molecular simulations are used to design cationic zeolites for boric acid adsorption. Force field models based on quantum mechanical calculations (PBE + D2) for Na-, Ca-, Mn-, and Fe-exchanged chabazite and LTA are presented. The new D2FF force fields reproduce DFT energies with about half the error of UFF. Zeolite performance depends on Si/Al ratio and cation type, with low Si/Al ratio chabazite (CHA) and phillipsite (PHI) zeolite frameworks exchanged with Ca 2+ or Na + /Ca 2+ mixtures showing the highest adsorption. In conclusion, these findings suggest tailored fly ash-derived zeolites could provide effective boron removal from leachate ponds.

CCR impoundment↗
DOE OSTI2021

Porous Covalent Organic Polymers for Efficient Fluorocarbon‐Based Adsorption Cooling

Abstract Adsorption‐based cooling is an energy‐efficient renewable‐energy technology that can be driven using low‐grade industrial waste heat and/or solar heat. Here, we report the first exploration of fluorocarbon adsorption using porous covalent organic polymers (COPs) for this cooling application. High fluorocarbon R134a equilibrium capacities and unique overall linear‐shaped isotherms are revealed for the materials, namely COP‐2 and COP‐3. The key role of mesoporous defects on this unusual adsorption behavior was demonstrated by molecular simulations based on atomistic defect‐containing models built for both porous COPs. Analysis of simulated R134a adsorption isotherms for various defect‐containing atomistic models of the COPs shows a direct correlation between higher fluorocarbon adsorption capacities and increasing pore volumes induced by defects. Combined with their high porosities, excellent reversibility, fast kinetics, and large operating window, these defect‐containing porous COPs are promising for adsorption‐based cooling applications.

Zheng, Jian↗
DOE OSTI2022

What Is the Smallest Zeolite That Could Be Synthesized?**

Abstract Zeolites with a few unit cells are promising as catalyst and adsorbents. The quest to synthesize the smallest zeolites has recently resulted in 4 to 8 nm nanozeolites, about 2 to 4 unit cells. These findings pose the question of what is the smallest zeolite that could be obtained by hydrothermal synthesis. Here we address this question using molecular simulations and thermodynamic analysis. The simulations predict that amorphous precursors as small as 4 nm can crystallize zeolites, in agreement with the experiments. We find that interfacial forces dominate the structure of smaller particles, resulting in size‐dependent compact isomers that have ring and pore distributions different from open framework zeolites. The instability of zeolites smaller than 3±0.5 nm precludes a classical mechanism of nucleation from solution or through assembly of small nanoslabs.

Dhabal, Debdas↗
DOE OSTI2025

Sub‐5 Ångstrom Porosity Tuning in Calixarene‐Derived Porous Liquids via Supramolecular Complexation Construction

Abstract Sub‐Ångstrom‐level porosity engineering, which is appealing in gas separations, has been demonstrated in solid carbon, polymer, and framework materials but rarely achieved in the liquid phase. In this work, a gas molecular sieving effect in the liquid phase at sub‐5 Ångstrom scale is created via sophisticated porosity tuning in calixarene‐derived porous liquids (PLs). Type II PLs are constructed via supramolecular complexation between the sodium salts of calixarene derivatives and crown ether solvents. The chemical structure variation and assembly behavior of the porous host upon PL construction are monitored by spectroscopy‐, X‐ray‐, and neutron‐scattering techniques. The presence of permanent porosity in calixarene‐derived PLs is verified by pressure swing gas uptake, altered CO 2 physisorption behavior, and molecular simulations. Sub‐5 Ångstrom porosity tuning within the PL phase is achieved by introducing bulky substituted groups on the benzene ring of the calixarene host, which then greatly affects the dynamic motion and transport behavior of CO 2 molecules and the Xe uptake performance. The approach being demonstrated in this work represents a promising pathway to tune and leverage the porosity effect for enhanced gas uptake capacity and selectivity in liquid sorbents.

Li, Errui [Department of Chemistry University of T↗
DOE OSTI2021

Porous Covalent Organic Polymers for efficient Fluorocarbon-based Adsorption Cooling

Adsorption-based cooling is an energy-efficient renewable-energy technology that can be driven using low-grade industrial waste heat and/or solar heat. Here, we report the first exploration of fluorocarbon adsorption using porous covalent organic polymers (COPs) for this cooling application. High fluorocarbon R134a equilibrium capacities and unique overall linear-shaped isotherms are revealed for the materials, namely COP-2 and COP-3. The key role of mesoporous defects on this unusual adsorption behavior was demonstrated by molecular simulations based on atomistic defect-containing models built for both porous COPs. Analysis of simulated R134a adsorption isotherms for various defect-containing atomistic models of the COPs shows a direct correlation between higher fluorocarbon adsorption capacities and increasing pore volumes induced by defects. Combined with their high porosities, excellent reversibility, fast kinetics, and large operating window, these defect-containing porous COPs are promising for adsorption-based cooling applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Peptoid-Directed Formation of Five-Fold Twinned Au Nanostars through Particle Attachment and Facet Stabilization

While bio-inspired synthesis offers great potential for controlling nucleation and growth of inorganic particles, precisely tuning biomolecule-particle interactions is a long-standing challenge. Herein, we used variations in peptoid sequence to manipulate peptoid-Au interaction, leading to synthesis of concave five-fold twinned, five-pointed Au nanostars via a process of repeated particle attachment and facet stabilization. Ex situ and liquid-phase TEM observations show that a balance between particle attachment biased to occur near the star points, preferential growth along the [100] direction, and stabilization of (111) facets is critical to forming star-shaped particles. Molecular simulations predict that interaction strengths between peptoids and distinct Au facets differ significantly and thus can alter attachment kinetics and surface energies to form the stars. This work provides new insights into how sequence-defined ligands affect particle growth to regulate crystal morphology.

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

What Is the Smallest Zeolite That Could Be Synthesized?**

Zeolites with a few unit cells are promising as catalyst and adsorbents. The quest to synthesize the smallest zeolites has recently resulted in 4 to 8 nm nanozeolites, about 2 to 4 unit cells. These findings pose the question of what is the smallest zeolite that could be obtained by hydrothermal synthesis. We address this question using molecular simulations and thermodynamic analysis. The simulations predict that amorphous precursors as small as 4 nm can crystallize zeolites, in agreement with the experiments. We find that interfacial forces dominate the structure of smaller particles, resulting in size-dependent compact isomers that have ring and pore distributions different from open framework zeolites. The instability of zeolites smaller than 3±0.5 nm precludes a classical mechanism of nucleation from solution or through assembly of small nanoslabs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗