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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Molecular Simulations Probing the Adsorption and Diffusion of Ammonia, Nitrogen, Hydrogen, and Their Mixtures in Bulk MFI Zeolite and MFI Nanosheets at High Temperature and Pressure

Recent advances in the synthesis of MFI zeolite nanosheets have led to highly selective membranes that are promising candidates for small-scale ammonia separation from ammonia/nitrogen/hydrogen mixtures in distributed green ammonia production plants. Using force-field-based molecular simulations, we evaluate the performance of bulk all-silica MFI zeolite and a 3 nm thick all-silica MFI nanosheet for ammonia/nitrogen/hydrogen separations at moderate temperature and pressure conditions (T = 373 K, p = 5 bar) as well as conditions relevant for a membrane-based reactor–separator process (T = 523 or 623 K, p = 80 bar). Isobaric–isothermal Gibbs ensemble Monte Carlo (GEMC) simulations were carried out to understand the selective adsorption behavior for these mixtures. Molecular dynamics simulations in the canonical ensemble were performed to examine the transport behavior of ammonia in a mixture with nitrogen or hydrogen at loadings obtained from the GEMC simulations. Our results show that both bulk MFI and the nanosheet with explicit surface silanols are highly selective toward ammonia adsorption, but the adsorption selectivity decreases by factors of 4 and 10 for ammonia/nitrogen and ammonia/hydrogen mixtures as the temperature is increased from 373 to 623 K. Conversely, the diffusion selectivities toward ammonia are more favorable at process-relevant temperatures and for the MFI nanosheet. Here, at T = 523 K and p = 80 bar, our simulations predict overall separation factors of 3.8 ± 0.2 and 4.5 ± 0.3 for ammonia/nitrogen and ammonia/hydrogen mixtures, respectively, in the bulk MFI zeolite, and separation factors of 8.6 ± 0.4 and 12.5 ± 0.8, respectively, for the MFI nanosheet.

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

Molecular Simulation of Lithium Carbonate Reactive Vapor–Liquid Equilibria Using a Deep Potential Model

We developed a first-principles machine learning model for the reactive vapor–liquid phase behavior of molten Li 2 CO 3 . The model was trained on ab initio electronic density functional theory data using the Deep Potential (DP) methodology, and its accuracy was evaluated by comparing model predictions of density and viscosity to experimental measurements. Direct coexistence simulations with the DP model over time scales of tens of nanoseconds were used to observe equilibrium dissociation of Li 2 CO 3 into CO 2 residing primarily in the vapor phase and Li 2 O which remains dissolved in the liquid. The simulations covered a range of temperatures, overall system sizes, and vapor-to-liquid volume ratios. Results were analyzed in terms of the observed chemical composition of the liquid and vapor phases, product structure, and CO 2 partial pressures. In addition, we calculated equilibrium constants for the dissociation reaction by assuming ideal-solution behavior for the liquid. As expected on the basis of thermodynamic arguments and prior experiments for this system, the observed partial pressure of CO 2 in the gas phase depends on both the temperature and the ratio of vapor to liquid volumes, while the calculated equilibrium constants only depend on temperature. DP model predictions for the equilibrium constant of the reaction are generally consistent with the available experimental measurements. Furthermore, the present study establishes the validity of the DP methodology for the description of reactive, multiphase equilibria from first principles, with possible applications to many other systems of scientific and technological interest even in the absence of relevant experimental measurements.

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

Nonane and Hexanol Adsorption in the Lamellar Phase of a Nonionic Surfactant: Molecular Simulations and Comparison to Ideal Adsorbed Solution Theory

Adsorption of n-nonane/1-hexanol (C9/C6OH) mixtures into the lamellar phase formed by a 50/50 w/w triethylene glycol mono-n-decyl ether (C10E3)/water system was studied using configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The interactions were described by the Shinoda–Devane–Klein coarse-grained force field. Prior simulations probing single-component adsorption indicated that C9 molecules preferentially load near the center of the bilayer increasing the bilayer thickness, whereas C6OH molecules are more likely to be found near the interface of the polar and non-polar moieties swelling the bilayer in the lateral dimension. Here, we extend this work to binary C9/C6OH adsorption to probe whether the difference in the spatial preferences may lead to a synergistic effect and enhanced loadings for the mixture. Comparing loading trends and the thermodynamics of binary adsorption to unary adsorption reveals that C9–C9 interactions lead to the largest enhancement, whereas C9–C6OH and C6OH–C6OH interactions are less favorable for this bilayer system. As a result, ideal adsorbed solution theory yields satisfactory predictions of the binary loading.

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

Ideal conductor/dielectric model (ICDM): A generalized technique to correct for finite-size effects in molecular simulations of hindered ion transport

Molecular simulations serve as indispensable tools for investigating the kinetics and elucidating the mechanism of hindered ion transport across nanoporous membranes. In particular, recent advancements in advanced sampling techniques have made it possible to access translocation timescales spanning several orders of magnitude. In our prior study [Shoemaker et al., J. Chem. Theory Comput. 18, 7142 (2022)], we identified significant finite size artifacts in simulations of pressure-driven hindered ion transport through nanoporous graphitic membranes. We introduced the ideal conductor model, which effectively corrects for such artifacts by assuming the feed to be an ideal conductor. In the present work, we introduce the ideal conductor dielectric model (Icdm), a generalization of our earlier model, which accounts for the dielectric properties of both the membrane and the filtrate. Using the Icdm model substantially enhances the agreement among corrected free energy profiles obtained from systems of varying sizes, with notable improvements observed in regions proximate to the pore exit. Moreover, the model has the capability to consider secondary ion passage events, including the transport of a co-ion subsequent to the traversal of a counter-ion, a feature that is absent in our original model. We also investigate the sensitivity of the new model to various implementation details. The Icdm model offers a universally applicable framework for addressing finite size artifacts in molecular simulations of ion transport. It stands as a significant advancement in our quest to use molecular simulations to comprehensively understand and manipulate ion transport processes through nanoporous membranes.

Chemistry↗
DOE OSTI2022

HPC Molecular Simulation Tries Out a New GPU: Experiences on Early AMD Test Systems for the Frontier Supercomputer

Molecular simulation is an important tool for nu- merous efforts in physics, chemistry, and the biological sciences. Simulating molecular dynamics requires extremely rapid cal- culations to enable sufficient sampling of simulated temporal molecular processes. The Hewlett Packard Enterprise (HPE) Cray EX Frontier supercomputer installed at the Oak Ridge Leadership Computing Facility (OLCF) will provide an exascale resource for open science, and will feature graphics processing units (GPUs) from Advanced Micro Devices (AMD). The future LUMI supercomputer in Finland will be based on an HPE Cray EX platform as well. Here we test the ports of several widely used molecular dynamics packages that have each made substantial use of acceleration with NVIDIA GPUs, on Spock, the early Cray pre-Frontier testbed system at the OLCF which employs AMD GPUs. These programs are used extensively in industry for pharmaceutical and materials research, as well as academia, and are also frequently deployed on high-performance computing (HPC) systems, including national leadership HPC resources. We find that in general, performance is competitive and installation is straightforward, even at these early stages in a new GPU ecosystem. Our experiences point to an expanding arena for GPU vendors in HPC for molecular simulation.

Sedova, Ada↗
DOE OSTI2026

Transfer Learning Meets Embedded Correlated Wavefunction Theory for Chemically Accurate Molecular Simulations: Application to Calcium Carbonate Ion Pairing

Achieving chemical accuracy for molecular simulations remains a central challenge in computational chemistry. Here, we present an embedded correlated wavefunction transfer learning (ECW-TL) framework for accurately simulating molecular dynamics in the condensed phase. ECW-TL incorporates high-level electron exchange and correlation effects in ECW theory while preserving the training and computational efficiency of machine-learned interatomic potentials. We demonstrate the framework on Ca 2+ –CO 3 2– ion pairing in aqueous solution, a key process underlying CO 2 mineralization in seawater. As proof of principle, we first show that fine-tuning a DFT-revPBE-D3(BJ) baseline model with embedded-DFT-SCAN data reproduces the DFT-SCAN free-energy surface within 1 kcal/mol across all solvation states. Extending the framework to embedded MP2 and localized natural-orbital CCSD(T) further refines the free-energy profile, revealing the crucial role of exact electron exchange and correlation in determining ion-pair stability and structure. The computed ion-pair association free energy is in quantitative agreement with experimental measurements, further validating the accuracy of the ECW-TL framework. ECW-TL thus provides a general, data-efficient route for transferring CW accuracy to efficient simulations of complex aqueous and interfacial chemical processes.

cluster chemistry↗
DOE OSTI2021

Improving scalability of electronic structure code for molecular simulations in the presence of environment

A scalable density functional electronic code with Gaussian basis set, called UTEP-NRLMOL, is developed to perform simulations of molecular systems in the presence of the environment with particular attention to the memory requirements. In the electronic structure calculations, the memory and computation time are proportional to the number of atoms. Memory requirements for density functional calculations scale as N*N, where N is the number of atoms. While the recent advances in HPC offer platforms with large numbers of cores, the limited amount of memory available on a given node and poor scalability of the electronic structure codes hinder their efficient usage of these platforms. We have introduced new scaling and parallelization paradigms using MPI-3 shared-memory functionality combined with usage of sparse algebra and storage of matrices in sparse format. This extends the range of applicability of the UTEP-NRLMOL code to large systems over 10,000 atoms, or using up to 67,000 basis functions, and making use of HPC architectures using over 6,000 processors utilizing all available cores. We have also interfaced code with effective fragment potential and polarizable continuum model libraries. The code was used in simulations of several applications which are published in reputed scientific journals.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2017

Dissociation and Internal Excitation of Molecular Nitrogen Due to N + N2 Collisions Using Direct Molecular Simulation

In this work we present a molecular level study of N2+N collisions, focusing on excitation of internal energy modes and non-equilibrium dissociation. The computation technique used here is the direct molecular simulation (DMS) method and the molecular interactions have been modeled using an ab−initio potential energy surface (PES) developed at NASA's Ames Research Center. We carried out vibrational excitation calculations between 5000K and 30000K and found that the characteristic vibrational excitation time for the N + N2 process was an order of magnitude lower than that predicted by the Millikan and White correlation. It is observed that during vibrational excitation the high energy tail of the vibrational energy distribution gets over populated first and the lower energy levels get populated as the system evolves. It is found that the non-equilibrium dissociation rate coefficients for the N + N2 process are larger than those for the N2 + N2 process. This is attributed to the non-equilibrium vibrational energy distributions for the N + N2 process being less depleted than that for the N2 +N2 process. For an isothermal simulation we find that the probability of dissociation goes as 1/T(sub tr) for molecules with internal energy (epsilon(sub int)) less than approximately 9.9eV, while for molecules with epsilon (sub int) greater than 9.9eV the dissociation probability was weakly dependent on translational temperature of the system. We compared non-equilibrium dissociation rate coefficients and characteristic vibrational excitation times obtained by using the ab-initio PES developed at NASA's Ames Research Center to those obtained by using an ab-initio PES developed at the University of Minnesota. Good agreement was found between the macroscopic properties and molecular level description of the system obtained by using the two PESs.

Grover, Maninder S.↗
DOE OSTI2022

Combining Physics and Machine Learning for the Next Generation of Molecular Simulation

Simulating molecules and atomic systems at quantum accuracy is a grand challenge for science in the 21 st century. Quantum-accurate simulations would enable the design of new medicines and the discovery of new materials. The defining problem in this challenge is that quantum calculations on large molecules, like proteins or DNA, are fundamentally impossible with current algorithms. In this work, we explore a range of different methods that aim to make large, quantum-accurate simulations possible. We show that using advanced classical models, we can accurately simulate ion channels, an important biomolecular system. We show how advanced classical models can be implemented in an exascale-ready software package. Lastly, we show how machine learning can learn the laws of quantum mechanics from data and enable quantum electronic structure calculations on thousands of atoms, a feat that is impossible for current algorithms. Altogether, this work shows that combining advances in physics models, computing, and machine learning, we are moving closer to the reality of accurately simulating our molecular world.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Sampling Rare Events in Aqueous Systems Using Molecular Simulations

Birth of a new distinct phase is a phenomenon encountered in a myriad of processes, and has wide ranging consequences in material processing, biological self-assembly, separations and several other processes. Several phase transitions are nucleation driven. The nucleation events occur over nanosecond timescales and involve hundreds to thousands of molecules. These length and timescales are difficult to access in experiments, thereby making experimental studies of nucleation challenging. On the other hand, molecular simulations sample the nanosecond and nanometer scales making them ideal to study nucleation. However, nucleation is a rare event, meaning that the waiting time to observe one nucleation event is significant. This makes simulation studies of rare events challenging. The project focused on a multi-pronged approach to address such challenges to develop the next generation rare event sampling methods for molecular simulations. The key outcomes of our work include developing more effective methods for sampling rare events, utilizing machine learning to better elucidate nucleation mechanisms, development of software for easy implementation of the methodologies, and applications of the methods to realistic systems to push the method applicability beyond model systems. Overall, this work has enabled pushing the frontiers of molecular simulations to study rare events with a focus on nucleation in aqueous solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

BETO 2021 Peer Review - Inverse Bioproduct Design Through Machine Learning and Molecular Simulation

This work aims to identify performance advantaged bioproducts (PABPs) through property prediction, which will guide experimental synthesis. The impact of this work will be faster market adoption of bioproducts with greater performance relative to incumbent products. We have identified >106 bioproduct candidates, but only some will have superior performance to create a market pull. High-throughput property prediction, enabled by machine learning, and elucidation of structure-function relationships, enabled by molecular simulation, provide a hypothesis driven approach for down selection of candidate biomolecules to pursue experimentally. To enable machine learning and molecular simulation for bioproduct discovery, automated structure generation and embedding must capture relevant features for prediction, databases must cover domains applicable to biobased products, and best practices for simulation of polymer systems must be developed. To address these challenges, we have established bioproduct relevant datasets, developed high-throughput polymer structure generation, and built end-to-end neural networks that have predicted 8 properties for >1.4 x 106 biopolymers. A molecular simulation pipeline for building, running, and analyzing polymers and polymer additives is being used to predict performance and develop design principles of biobased products. In collaboration with the PABP synthesis project, these computational tools are guiding synthesis and informing design of PABPs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Using molecular simulation to understand the skin barrier

Skin's effectiveness as a barrier to permeation of water and other chemicals rests almost entirely in the outermost layer of the epidermis, the stratum corneum (SC), which consists of layers of corneocytes surrounded by highly organized lipid lamellae. As the only continuous path through the SC, transdermal permeation necessarily involves diffusion through these lipid layers. The role of the SC as a protective barrier is supported by its exceptional lipid composition consisting of ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs) and the complete absence of phospholipids, which are present in most biological membranes. Molecular simulation, which provides molecular level detail of lipid configurations that can be connected with barrier function, has become a popular tool for studying SC lipid systems. We review this ever-increasing body of literature with the goals of (1) enabling the experimental skin community to understand, interpret and use the information generated from the simulations, (2) providing simulation experts with a solid background in the chemistry of SC lipids including the composition, structure and organization, and barrier function, and (3) presenting a state of the art picture of the field of SC lipid simulations, highlighting the difficulties and best practices for studying these systems, to encourage the generation of robust reproducible studies in the future. This review describes molecular simulation methodology and then critically examines results derived from simulations using atomistic and then coarse-grained models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Combining Molecular Simulations and Machine Learning to Comprehensively Explore Adsorption Space (Final Report)

This is the final report for the project "Combining Molecular Simulations and Machine Learning to Comprehensively Explore Adsorption Space". The central aim of this project was to develop modeling tools that vastly increase the diversity of molecules for which reliable adsorption isotherms in nanoporous materials are available. Within the project, molecular simulation data for a diverse set of examples was combined with machine-learning (ML) methods to rapidly predict adsorption equilibria in crystalline nanoporous materials (primarily MOFs).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of Molecular Simulations in the Design of Metal–Organic Frameworks for Gas-Phase Thermocatalysis: A Perspective

Here, metal organic frameworks (MOFs) are highly tunable porous crystalline solids with spatially and electronically isolated catalytically active sites that have been demonstrated for a variety of thermo-, redox-, and photocatalytic reactions. Their tunable natures and relatively well-defined active sites make them advantageous for catalyst design, and molecular simulations have proven highly valuable in this endeavor. However, complexities in the MOF structure require advanced simulation strategies that accurately capture quantum chemistry at the strongly correlated transition metal cation active sites, compositional and structural changes caused by finite temperature and pressure reaction conditions, and transport effects in variously sized pore environments. Luckily, several groups have started using such simulation strategies, paving the way for rich opportunities to design MOF catalysts. Herein, we highlight such examples and provide a perspective on the needs for molecular simulations in the design of MOF catalysts moving forward.

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

Nuclear magnetic resonance and molecular simulation study of H 2 and CH 4 adsorption onto shale and sandstone for hydrogen geological storage

Understanding pure H 2 and H 2 /CH 4 adsorption and diffusion in earth materials is one vital step toward a successful and safe H 2 storage in depleted gas reservoirs. Despite recent research efforts such understanding is far from complete. In this work we first use Nuclear Magnetic Resonance (NMR) experiments to study the NMR response of injected H 2 into Duvernay shale and Berea sandstone samples, representing materials in confining and storage zones. Then we use molecular simulations to investigate H 2 /CH 4 competitive adsorption and diffusion in kerogen, a common component of shale. Our results indicate that in shale there are two H 2 populations, i.e., free H 2 and adsorbed H 2 , that yield very distinct NMR responses. However, only free gas presents in sandstone that yields a H 2 NMR response similar to that of bulk H 2 . About 10 % of injected H 2 can be lost due to adsorption/desorption hysteresis in shale, and no H 2 loss (no hysteresis) is observed in sandstone. Here, our molecular simulation results support our NMR results that there are two H 2 populations in nanoporous materials (kerogen). The simulation results also indicate that CH 4 outcompetes H 2 in adsorption onto kerogen, due to stronger CH 4 -kerogen interactions than H 2 -kerogen interactions. Nevertheless, in a depleted gas reservoir with low CH 4 gas pressure, about ~30 % of residual CH 4 can be desorbed upon H 2 injection. The simulation results also predict that H 2 diffusion in porous kerogen is about one order of magnitude higher than that of CH 4 and CO 2 . This work provides an understanding of H 2 /CH 4 behaviors in deleted gas reservoirs upon H 2 injection and predictions of H 2 loss and CH 4 desorption in H 2 storage.

08 HYDROGEN↗
DOE OSTI2020

Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein

Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BP SS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BP SS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BP SS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BP SS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R 0 ) and fluorophore orientation factor (κ 2 ) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Understanding Gas Solubility of Pure Component and Binary Mixtures within Multivalent Ionic Liquids from Molecular Simulations

Understanding the molecular-level solubility of CO 2 and its mixtures is essential to the progress of gas treating technologies. Herein, we use grand canonical Monte Carlo simulations to study the single component gas absorption of SO 2 , N 2 , CH 4 , and H 2 , and binary mixtures of CO 2 /SO 2 , CO 2 /N 2 , CO 2 /CH 4 , and CO 2 /H 2 of varying mole fractions within multivalent ionic liquids (ILs). Our results highlight the importance of the free volume effect and the anion effect when interpreting the absorption behavior of these mixtures, similar to the behavior of CO 2 found in our previous study. The deviation of the gas solubility between the pure component absorption versus the binary absorption, as well as the solubility selectivity, highlight the importance of the relative affinity of gas species within a mixture to the different anions. Further, the absorption selectivity within a specific IL system can be predicted based on the relative gas affinity to the anion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗