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

Relationship between microporous structure and light gas transport through glassy polymeric membranes revealed by molecular simulations

Microporous glassy polymers are attractive materials for gas separation membranes, due to their high permeability and tailorable selectivity, provided physical aging can be delayed. The archetypal microporous glassy polymer, PTMSP, can be blended with a hyper-crosslinked isatin–triptycene porous polymer network (PPN) to delay physical aging. However, while PPN is effective at reducing physical aging, it also affects the permeability of light gases through PTMSP. Molecular dynamics simulations were used here to shed fundamental light on the mechanisms responsible for these effects. Atomistic models are developed that satisfactorily reproduce experimental observations such as matrix density and cavity size distributions for neat PTMSP as well as for PTMSP–PPN blends. Analysis of the simulation results suggests that physical aging is delayed because the PPN inclusions slow down PTMSP relaxation while reducing the connectivity between free volume pockets. To understand how PPN inclusions affect light gas permeability, the atomistic models developed are used to probe CO 2 and CH 4 diffusion and sorption. These simulations are conducted for PTMSP matrices exhibiting varying density, towards reproducing experimental permeability data. Interrogating the simulation trajectories, it is discovered that while CH 4 travels preferentially through the free volume cavities, CO 2 preferentially interacts with the available surfaces, especially in the presence of PPN. These differences help interpret experimental observations. The transport of C 2 H 6 and H 2 S through PTMSP matrices was also investigated. These gases, and in particular H 2 S, were found to absorb within PTMSP, yielding very low diffusion coefficients. The ability to predict differences in diffusion pathways for various natural gas components unveils the possibility of engineering membranes to control permeability and selectivity towards large-scale applications.

Bao Le, Tran Thi [Univ. of Oklahoma, Norman, OK (U↗
NASA NTRSDate not supplied

Validation of Accurate Cryogenic Fluid Vapor-Liquid Boundary Conditions Via Molecular Simulations

Vapor-liquid interfaces drive many important phenomena in cryogenic fluid management, including heat transfer, evaporation, and capillary flow. Design of cryogenic fluid systems, such as propellant storage, requires accurate predictions of fluid behavior, including evaporation rates. Many models have been proposed for heat and mass transfer at vapor liquid interfaces, but the accuracy of these models in the context of cryogenic fluids has not been performed. We use molecular dynamics simulations, which allow for nanometer scale resolution of fluid phenomena, to evaluate the accuracy of a variety of vapor-liquid boundary conditions at evaporating and condensing interfaces. We find that an anisotropic temperature distribution is a critical ingredient for accurate prediction of intensive evaporation and condensation.

Daniel Vigil↗
DOE OSTI2022

MISPR : an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Antiviral discovery using sparse datasets by integrating experiments, molecular simulations, and machine learning

Computational methods have demonstrated success in identifying virucidal agents, effectively contributing to the discovery of novel virucidal molecules. In this study, we developed a machine learning (ML) model, trained on a small dataset, to predict inhibitors of human enterovirus 71 (EV71), a pathological agent that causes severe disease in children and immunocompromised adults. Despite the dataset’s limitation, comprising of only 36 compounds tested, our ML framework demonstrated significant predictive capability. Notably, experimental validation revealed that five out of the eight compounds predicted by our model from the Chinese cosmetic material list exhibited virucidal activity. The inhibitor effects displayed by the main active compounds were further confirmed by molecular dynamics simulation. This underscores the potential of our AI-driven approach to bypass data constraints in identifying active molecules against viral pathogens.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Effect of Collective Dynamics and Anharmonicity on Entropy in Heterogenous Catalysis: Building the Case for Advanced Molecular Simulations

In this paper, we present a perspective on the computational determination of entropy and its effects and consequences on heterogeneous catalysis. Special attention is paid to the role of anharmonicity and the deviations from the standard harmonic oscillator approximations which can fail to provide a reliable assessment of entropy. To address these challenges, advanced methodologies are needed that can explicitly account for these thermodynamic drivers through the appropriate statistical sampling of reactive free energy surfaces. We discuss where anharmonicity should be expected, where it has been observed from a theoretical perspective, and the methods currently employed to address it. We concentrate on three types of systems where we have observed major, nonnegligible anharmonic effects: (1) supported nanoparticles, where the migration of metal atoms, complexes, and entire clusters exhibit anharmonic behavior in their dynamic motion; (2) porous solids where confinement effects distort potential energy surfaces and hinder molecular motions resulting in large entropic terms; and (3) solid/liquid interfaces, where interactions between solvent molecules and adsorbed species can result in large solvent organization free energy and unique reactivity.

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

Unravelling the dynamics of the maturation protein in MS2 bacteriophage via molecular simulations

The MS2 bacteriophage capsid serves as a model system for studying viral structure and function. Mature MS2 virus consists of 178 capsid proteins and a single maturation protein (MP), which is essential for host receptor binding and infection initiation. Despite its critical role, the dynamic behavior of the capsid with the MP remains poorly understood. To address this, we conducted 0.5 µs all-atom molecular dynamics (MD) simulations of the MS2 capsid with and without the MP, revealing key insights into its structural dynamics. Our simulations showed that MP exhibits high flexibility, particularly in the “tip” and “side-loop” regions, which undergo significant motions that likely enhance its ability to engage with the F-pilus receptor. Detailed analysis of MP conformational states revealed that loop rearrangements around H357 enable transient switching between “semi-closed” and “open” conformations, suggesting a conformational selection mechanism for pilus binding. Additionally, ion interaction analyses revealed distinct sodium and chloride binding patterns, where sodium ions were mostly found at the outer capsid shell, while chloride ions interacted with the basic residues on the RNA-facing side. We also found that the presence of the MP enhances salt-bridge interactions, contributing to increased capsid stability, yet it does not significantly alter the pore sizes of pentameric and hexameric units. Together, these findings provide new insights into the functional role of the MP, highlighting its contribution to capsid stability and host receptor engagement. This study offers a foundation for understanding capsid dynamics relevant to viral infectivity and may guide future rational strategies aimed at disrupting host-virus interactions.

Capsid stability↗
DOE OSTI2024

Enhanced Interfacial Electron Transfer in Photocatalyst-Natural Enzyme Coupled Artificial Photosynthesis System: Tuning Strategies and Molecular Simulations

Laccase is capable of catalyzing a vast array of reactions, but its low redox potential limits its potential applications. The use of photocatalytic materials offers a solution to this problem by converting absorbed visible light into electrons to facilitate enzyme catalysis. Herein, MIL-53(Fe) and NH 2 -MIL-53(Fe) serve as both light absorbers and enzyme immobilization carriers, and laccase is employed for solar-driven chemical conversion. Electron spin resonance spectroscopy results confirm that visible light irradiation causes rapid transfer of photogenerated electrons from MOF excitation to T1 Cu(II) of laccase, significantly increasing the degradation rate constant of tetracycline (TC) from 0.0062 to 0.0127 min -1 . Conversely, there is only minimal or no electron transfer between MOF and laccase in the physical mixture state. Theoretical calculations demonstrate that the immobilization of laccase's active site and its covalent binding to the metal-organic framework surface augment the coupled system's activity, reducing the active site accessible from 27.8 to 18.1 Å. In conclusion, the constructed photo-enzyme coupled system successfully combines enzyme catalysis’ selectivity with photocatalysis's high reactivity, providing a promising solution for solar energy use.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Data-Centric Development of Lignin Structure–Solubility Relationships in Deep Eutectic Solvents Using Molecular Simulations

Lignin is a natural source of aromatic chemicals with significant potential as an abundant, renewable feedstock for value-added products. Deep eutectic solvents (DES)–solvents composed of a hydrogen bond donor (HBD) and acceptor (HBA) in varying ratios–have emerged as a highly tunable class of solvents for lignin solubilization. However, the variety of possible DES compositions and limited molecular-scale understanding of lignin solubility makes solvent selection a challenge without laborious trial-and-error experimentation. To address these challenges, we use classical molecular dynamics (MD) simulations to study the interactions of lignin model compounds with various DES–water systems. Quantitative parameters (descriptors) were calculated by postprocessing the MD results and used to train a regression model that predicts experimentally determined solubilities of lignin model compounds. This approach revealed that the most important descriptors of solubility are the system temperature, solute hydrophilicity, and metrics quantifying hydrogen bonding. Maximizing the interactions between solute–HBD (hydrophobic group), water–HBD (hydrophilic group), and water–HBA molecules led to the highest model compound solubility. Our results support a hydrotropic mechanism in which extensive DES–water hydrogen bonding and favorable HBD interactions with the solute promote high solubility. We applied the regression model derived using model compounds to predict the solubility of representative lignin oligomers. The model predicted lignin oligomers’ solubilities in good agreement with experiments, indicating that the simulations of model compounds can be extended to predict the solubility of larger lignin compounds across a range of solvent compositions and temperatures. Furthermore, these findings provide new molecular-scale insight into lignin solubilization mechanisms and a new method for computationally screening potential solvent systems for lignin valorization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Interfacial Structure and Capacitance of Li-doped Ionic Liquid Electrolytes from Molecular Simulations

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries, supercapacitors, and hybrid energy storage devices. Though Li-salt is often present in these systems, its influence on interfacial properties is largely uncharacterized. We, thereby, present an extensive computational analysis, supported by experimental comparisons, of the properties of a representative set of these electrolytesat an ideal carbon interface as a function of Li-salt doping and voltage. We have performed polarizable molecular (MD) dynamics simulations, using the APPLEP force field, to evaluate electric double layer (EDL) capacitance and distribution of Li+ in the EDL. Differential capacitance exhibits the characteristic camel profile and is insensitive to Li-doping. Li+ localizes in the second molecular layer of the EDL, which is a result of confinement from free energy barriers associated with ion layering. Joint MDelectronic structure computations show the electrochemical window of the electrolytes to be a weak function of Li-doping. Estimates of supercapacitor specific energy are made using the computed window and capacitance. The magnitude and trends in specific energy are in good agreement with experiment.

batteries↗
DOE OSTI2025

Molecular simulation and artificial intelligence for the circular economy of bioenergy and bioproducts

The concept of the circular bioeconomy is a carbon neutral, sustainable system with zero waste. One vision for such an economy is based upon lignocellulosic biomass. This lignocellulosic circular bioeconomy requires CO 2 absorption from biomass growth and the efficient deconstruction of recalcitrant biomass into solubilized and fractionated biopolymers, which are then used as precursors for the sustainable production of high-quality liquid fuels, chemical bioproducts, and bio-based materials. Here, in this study, we summarize the roles that molecular dynamics (MD) simulations and machine learning (ML) are playing in overcoming several fundamental challenges hindering the adoption of a circular bioeconomy. Specifically, we discuss the role of MD and ML/AI in overcoming lignocellulose recalcitrance by designing biomass pretreatment methods to efficiently produce solubilized cellulose/lignin/hemicellulose and of that in improving energy-intensive manufacturing of biomass-based materials and their structural and mechanical properties. Quantum mechanical methods and MD simulations, in addition to offering a mechanistic understanding of biomass deconstruction and biomaterials design, can provide meaningful structural, energetics, and physiochemical properties as inputs to train AI/ML models. The ML models can guide the experimental prioritization of materials/solvents and process parameters that significantly accelerate the development of biofuel and biomaterial components of the circular bioeconomy.

Smith, Jeremy C. [Oak Ridge National Laboratory (O↗
DOE OSTI2026

Quantifying Membrane Structure and Dynamics during Bioproduct Production in Zymomonas mobilis by Molecular Simulation

The conversion of lignocellulosic biomass into biofuels and bioproducts by microbial biorefineries is central to a sustainable chemical industry. Zymomonas mobilis is one such biorefinery chassis and is resistant to ethanol stress, leading to its use in biomass conversion to biofuels and bioproducts. However, Z. mobilis growth is often inhibited by organic acids, aldehydes, alcohols, ketones, and amides found in biomass hydrolysate. The resulting slow growth inhibits production and as a result drives up the price for the resulting products. One hypothesis is that these molecules interact with or disrupt the bacterial membrane, triggering stress responses and hindering growth. To test this hypothesis at the molecular level, we employ all-atom molecular dynamics (MD) simulations to investigate lignocellulose-derived small molecules and their impact on a biologically relevant Z. mobilis membrane model. Simulations were conducted across a range of inhibitor concentrations from 0 to 2.5 mol %, analyzing key membrane properties such as area per lipid (APL), membrane thickness, lipid-order parameter (−S CH ), lateral diffusion coefficient (D xy ), and permeability coefficient (Pm). From simulation, we observed altered membrane structure and dynamics at these modest small molecule concentrations commonly found in hydrolysates. Generally, the membranes become thinner, with a higher area per lipid and lower-order parameter as the small molecule concentration increases. These trends are stronger for more hydrophobic molecules with greater hydrophobic bulk, as isobutanol, propanol, and propanoic acid showed greater membrane perturbations as the concentration increased compared to other small molecules. Tracking small molecule distributions directly in our equilibrium simulations allows us to determine concentration-dependent free energy profiles for these molecules. While the trends are noisy, generally the barriers to crossing the membrane decrease as the concentration increases, indicating that the membranes become leakier as small molecule concentrations rise. Comparing between native Z. mobilis membranes with hopanoids and membranes sharing the same phospholipid composition but without hopanoids, hopanoids stabilize and order the membrane for smaller molecules to maintain membrane structure but appear insufficient for larger hydrophobic molecules like isobutanol. These findings provide a mechanistic understanding of how small molecules found in biomass degradation streams interact with the Z. mobilis membrane, offering valuable insights for future strain engineering efforts to optimize biofuel and bioproduct synthesis from biomass feedstocks by highlighting limits to small molecule tolerance. This knowledge can guide the modification of membrane composition to develop more robust microbes, thereby improving microbial survival and yields in industrial contexts.

Singh, Nitin Kumar [Michigan State Univ., East Lan↗
DOE OSTI2021

Machine-learning accelerated geometry optimization in molecular simulation

Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2017

Small Ion Channel Linking Molecular Simulations and Electrophysiology

Ion channels are pore-forming protein assemblies that mediate the transport of small ions across cell membranes. Otherwise, membrane bilayers would be almost impermeable to ions incapable to traverse the low dielectric constant, hydrophobic membrane core. Ion channels are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses, cardiac functions, muscle contraction and apoptosis. On the other extreme of biological complexity, viral ion channels (viroporins) influence many stages of the virus infection cycle either through regulating virus replication, such as entry, assembly and release or modulating the electrochemical balance in the subcellular compartments of host cells. Ion channels were crucial components of protocells. Their emergence facilitated adaptation of nascent life to different environmental conditions. The earliest ion channels must have been much simpler than most of their modern ancestors. Viral channels are among only a few naturally occurring models to study the structure, function and evolution of primordial channels. Experimental studies of these properties are difficult and often unreliable. In principle, computational methods, and molecular dynamics (MD) simulations in particular, can aid in providing information about both the structure and the function of ion channels. However, MD suffers from its own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. It is, therefore, essential to determine the reliability of MD simulations. We propose to do so on the basis of two criteria. One is channel stability on time scales that extend for several microseconds or longer. The other is the ability to reproduce the measured ionic conductance as a function of applied voltage. If both the stability and the calculated ionic conductance are satisfactory it will greatly increase our confidence that the structure and the function of a channel are described sufficiently accurately. To our knowledge, long time scale stability (approx.10 micro-sec) and the correct electrophysiology have been shown so far for only one channel - the synthetic LS3 hexamer). In this presentation, this approach will be discussed in application to two viral channels - Vpu, encoded by the HIV-1 genome and p7 of hepatitis C.

Pohorille, Andrzej↗
DOE OSTI2023

Multiscale molecular simulations for the solvation of lignin in ionic liquids

Lignin, the second most abundant biopolymer found in nature, has emerged as a potential source of sustainable fuels, chemicals, and materials. Finding suitable solvents, as well as technologies for efficient and affordable lignin dissolution and depolymerization, are major obstacles in the conversion of lignin to value-added products. Certain ionic liquids (ILs) are capable of dissolving and depolymerizing lignin but designing and developing an effective IL for lignin dissolution remains quite challenging. To address this issue, the COnductor-like Screening MOdel for Real Solvents (COSMO-RS) model was used to screen 5670 ILs by computing logarithmic activity coefficients (ln(γ)) and excess enthalpies (H E ) of lignin, respectively. Based on the COSMO-RS computed thermodynamic properties (ln(γ) and H E ) of lignin, anions such as acetate, methyl carbonate, octanoate, glycinate, alaninate, and lysinate in combination with cations like tetraalkylammonium, tetraalkylphosphonium, and pyridinium are predicted to be suitable solvents for lignin dissolution. The dissolution properties such as interaction energy between anion and cation, viscosity, Hansen solubility parameters, dissociation constants, and Kamlet–Taft parameters of selected ILs were evaluated to assess their propensity for lignin dissolution. Furthermore, molecular dynamics (MD) simulations were performed to understand the structural and dynamic properties of tetrabutylammonium [TBA] + -based ILs and lignin mixtures and to shed light on the mechanisms involved in lignin dissolution. MD simulation results suggested [TBA] + -based ILs have the potential to dissolve lignin because of their higher contact probability and interaction energies with lignin when compared to cholinium lysinate.

09 BIOMASS FUELS↗
DOE OSTI2022

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗