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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Elucidating the temperature and density dependence of silver chloride hydration numbers in high-temperature water vapor: A first-principles molecular simulation study

Hydration numbers of metal complexes in low-density aqueous solutions are required for developing geochemical models for ore-forming metals and for designing supercritical desalination processes. In this study, we investigate the temperature and density dependence of the hydration numbers of silver chloride at 623 K, 673 K, and 713 K and densities of 10–100 kg/m 3 . Experimental estimates of the hydration number in the literature for AgCl at these conditions are inconclusive and possibly contradictory as to the temperature and density dependence. First-principles molecular simulation presents an attractive alternative to experimental measurements. Specifically, recent work shows that machine-learning-accelerated nested Monte Carlo simulations provide reliable estimates for the hydration numbers of CuCl at 623 K from 10 to 100 kg/m 3 . Using the same technique, we find a monotonic temperature dependence, with the hydration number decreasing slightly with increasing temperature. In addition, the simulation-predicted hydration numbers steadily increase with increasing density. These temperature and density trends are in agreement with certain experimental data sets. Therefore, this work demonstrates how first-principles Monte Carlo simulations assist in resolving discrepancies between experimental data sets. Our simulation results also correctly predict that the hydration number, and thus also the solubility, of AgCl is lower than that of CuCl under the same conditions. Furthermore, the bond length and angle formed between the water complex and AgCl differ from those for CuCl, consistent with the lower solubility of AgCl.

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

Molecular Simulation Models of Carbon Dioxide Intercalation in Hydrated Sodium Montmorillonite

In this study, classical molecular dynamics simulations and density functional theory (DFT)-based molecular dynamics are used to elucidate the process of CO 2 intercalation into hydrated Na-montmorillonite at P-T conditions relevant to geological formations suitable for CO 2 storage. Of particular interest are the structural and transport properties of interlayer species after CO 2 intercalation. The conducted simulations allowed the research team to quantify expansion/contraction of smectite as a function of CO 2 and H 2 O compositions. The resulting swelling curves can be used to gauge the amount of stored CO 2 , compare it to the experiment, and estimate changes in geomechanical properties of the storage formation. The obtained results showed that the infrared signal of the asymmetric stretch vibration of CO 2 molecule is extremely sensitive to the solvent environment. The extent of the frequency shift relative to the gas-phase value can be used to probe hydration level in the interlayer with intercalated CO 2 . Interaction of supercritical CO 2 with brine in deep geological formations promotes an increase of hydrophobicity of clay surfaces. As a result of wettability alteration, estimated diffusion constants of CO 2 and H 2 O increase with the increased CO 2 load; this can contribute to faster migration of CO 2 throughout the formation.

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

Transferable Kinetic Monte Carlo Models with Thousands of Reactions Learned from Molecular Dynamics Simulations

Molecular dynamics (MD) simulation of complex chemistry typically involves thousands of atoms propagating over millions of time steps, generating a wealth of data. Traditionally these data are used to calculate some aggregate properties of the system and then discarded, but we propose that these data can be reused to study related chemical systems. Using approximate chemical kinetic models and methods from statistical learning, we study hydrocarbon chemistries under extreme thermodynamic conditions. We discover that a single MD simulation can contain sufficient information about reactions and rates to predict the dynamics of related yet different chemical systems using kinetic Monte Carlo (KMC) simulation. Our learned KMC models identify thousands of reactions and run four orders of magnitude faster than MD. In conclusion, the transferability of these models suggests that we can viably reuse data from existing MD simulations to accelerate future simulation studies and reduce the number of new MD simulations required.

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

Long-Range Dispersion Governs the Behavior of Near-Critical Fluids: Universal Scaling and Implications for Accurate Molecular Simulation

Computer simulations of near-critical and supercritical fluids often deviate from experimental results, a discrepancy commonly attributed to force-field resolutions and inaccuracies. We demonstrate that these errors primarily stem from using finite cutoff schemes for dispersion interactions, which become effectively long-range as the correlation length grows near the critical point. By employing the smooth particle-mesh Ewald method to account for full long-range dispersion, we show that coarse-grained models can also achieve high quantitative accuracy for n-undecane and propylbenzene. This approach enables the precise determination of density fluctuations and correlation lengths, which exhibit universal critical scaling and depend on the molecular size and shape. Our findings resolve a longstanding debate in molecular simulations of near-critical fluids, establishing that capturing long-wavelength fluctuations is essential for bridging the gap between the microscopic behavior and macroscopic critical phenomena.

Chen, Guang↗
DOE OSTI2015

Application of bi-Helmholtz nonlocal elasticity and molecular simulations to the dynamical response of carbon nanotubes

The nonlocal theory of elasticity is employed for the study of the free vibrations of carbon nanotubes (CNT). For the first time, a bi-Helmholtz operator has been used instead of the standard Helmholtz operator in a nonlocal beam model. Alongside the continuum formulation and its numerical solution, atomistic Molecular Dynamics (MD) simulations have been conducted in order to directly evaluate the eigenfrequencies of vibrating CNTs with a minimum of adjustable parameters. Our results show that the bi-Helmholtz operator is the most appropriate one to fit MD simulation results. However, the estimation of vibration eigenfrequencies from molecular simulations still remains an open (albeit well-posed) problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Molecular Simulations of Heterogeneous Ice Nucleation. II. Peeling back the Layers

Coarse grained molecular dynamics simulations are presented in which the sensitivity of the ice nucleation rate to the hydrophilicity of a graphene nanoflake is investigated. We find that an optimal interaction strength for promoting ice nucleation exists, which coincides with that found previously for a face centered cubic (111) surface. We further investigate the role that the layering of interfacial water plays in heterogeneous ice nucleation and demonstrate that the extent of layering is not a good indicator of ice nucleating ability for all surfaces. Our results suggest that to be an efficient ice nucleating agent, a surface should not bind water too strongly if it is able to accommodate high coverages of water.

Cox, Stephen J.↗
DOE OSTI2022

Molecular Simulations of CH 4 and CO 2 Diffusion in Rigid Nanoporous Amorphous Materials

Molecular diffusion in nanoporous materials is important in determining the rate of equilibration of various adsorption processes and plays a pivotal role in kinetic separations and membrane-based separations. Because generating realistic structures of amorphous nanoporous materials is difficult, far less is known about diffusion in amorphous nanoporous materials than in their crystalline counterparts. Here, we present molecular dynamics simulations assessing the room-temperature self-diffusion of CH 4 and CO 2 in a wide range of rigid amorphous nanoporous materials, including porous carbons, kerogens, polymers of intrinsic microporosity, and hyper-cross-linked polymers. Our results are the largest collection of molecular diffusivities in amorphous nanoporous materials to date. In each material, the diffusivity increases with the adsorbate concentration at low and moderate adsorbate concentrations, reaching a maximum before decreasing due to steric effects at higher concentrations. The observed diffusivities are much slower than that would be expected based on standard descriptions of Knudsen diffusivity. Here we show that the observed diffusivities are not correlated in a simple way with scalar descriptors of the pore structures such as the pore limiting diameter.

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

Xe Recovery from Nuclear Power Plants Off-Gas Streams: Molecular Simulations of Gas Permeation through DD3R Zeolite Membrane

Recent experimental work has shown zeolite membrane-based separation as a promising potential technology for Kr/Xe gas mixtures due to its much lower energy requirements in comparison to cryogenic distillation, the conventional separation method for such mixtures. Such a separation is also economically rewarding because Xe is in high demand, as a valuable product for many applications/processes. In this work, we have used Molecular Dynamics (MD) simulations to study the effects of different conditions, i.e., temperature, pressure, and gas feed composition, on Kr/Xe separation performance via DD3R zeolite membranes. We provide a comprehensive study of the permeation of the different gas species, density profiles, and diffusion coefficients. Molecular simulations show that if the feed is changed from pure Kr/Xe to an equimolar mixture, the Kr/Xe separation factor increases, which agrees with experiments. In addition, when Ar is introduced as a sweep gas, the adsorption of both Kr and Xe increases, while the permeation of pure Kr increases. A similar behavior is observed with equimolar mixtures of Kr/Xe with Ar as the sweep gas. High-separation Kr/Xe selectivity is observed at 50 atm and 425 K but with low total permeation rates. Changing pressure and temperature are found to have profound effects on optimizing the separation selectivity and the permeation throughput.

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

Molecular Simulations of the Hydration Behavior of a Zwitterion Brush Array and Its Antifouling Property in an Aqueous Environment

We carried out umbrella sampling and molecular dynamics (MD) simulations to investigate molecular interactions between sulfobetaine zwitterions or between sulfobetaine brushes in different media. Simulation results show that it is more energetically favorable for the two sulfobetaine zwitterions or brushes to be fully hydrated in aqueous solutions than in vacuum where strong ion pairs are formed. Structural properties of the hydrated sulfobetaine brush array and its antifouling behavior against a foulant gel are subsequently studied through steered MD simulations. We find that sulfobetaine brush arrays with different grafting densities have different structures and antifouling mechanisms. At a comparably higher grafting density, the sulfobetaine brush array exhibits a more organized structure which can hold a tightly bound hydration water layer at the interface. Compression of this hydration layer results in a strong repulsive force. However, at a comparably lower grafting density, the brush array exhibits a randomly oriented structure in which the antifouling of the brush array is through the deformation of the sulfobetaine branches.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Molecular simulations of competitive adsorption of carbon dioxide – methane mixture on illitic clay surfaces

Molecular dynamics and Monte Carlo simulation studies were carried out here to investigate adsorption behavior of carbon dioxide and methane mixtures on illitic clay surfaces under dry conditions. Various compositions of the mixtures and distributions of isomorphic substitutions in clay layers were chosen to explore competitive adsorption depending on component concentration and charge localization. The simulations show that carbon dioxide is preferentially sorbed on the illitic surface and is capable to promote methane desorption. Density distributions of the molecular species in pore space reveal formation of multilayers on the clay surfaces at elevated pressures. Finally, mixed adsorption isotherms were compared with adsorption isotherms of pure compounds and thermodynamic quantities were reported to characterize the interaction of the carbon dioxide and methane with the clay surface.

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

Molecular Simulation Study of Piezoelectric Potential Distribution in a ZnO Nanowire under Mechanical Bending

ABSTRACT Molecular dynamics (MD) simulations are carried out to investigate the interfacial sliding dynamics of a ZnO piezoelectric nanogenerator in an atomic force microscope (AFM) setting. The molecular system includes a vertically aligned ZnO nanowire along the [0001] direction and a Pt (111) metal tip sliding over it. We calculate the piezoelectric potential distributions based on the equilibrium molecular configurations using classical ionic charges. Simulation results reveal the very detailed evolution changes of the piezopotential within the nanowire, which are largely contributed from the different internal tensile or compressive strains induced by mechanical bending of the ZnO nanowire. Variations of the normal contact and lateral frictional forces versus the sliding distance of the Pt metal tip are also presented.

Tan, Dan↗
DOE OSTI2025

Molecular Simulation of Functionalized Covalent Organic Framework Membranes for Inorganic Salt Separation

Covalent organic frameworks (COFs) enable molecular-level design of nanochannels for selective separation in pressure-driven membrane processes. Through variation of building blocks and, subsequently, the pore structure and chemistry, membrane performance can be tailored. This study employs nonequilibrium molecular dynamics simulations to theoretically demonstrate the tunable selectivity of COF membranes through a bottom-up functionalization approach. Water and salt transport are evaluated for six β-ketoenamine-linked COFs with varying multilayer thicknesses. For the thinnest multilayer (0.64–0.72 nm), all COFs exhibit low sodium sulfate (Na 2 SO 4 ) rejection (55–66%). However, 20 stacked sheets (6.4–7.2 nm) provide 67–98% rejection, with the sulfonated COF providing the highest Na 2 SO 4 rejection. Analysis of time-resolved ion density profiles reveals that solute rejection is primarily governed by interfacial exclusion arising from pore size and functional group chemistry. Although increasing salt rejection compromises water permeance, the permeance of all COF membranes is at least two orders of magnitude greater than that of a commercially available nanofiltration membrane. Overall, this work guides the rational design of COF membranes for aqueous salt separation.

Nanofiltration↗
DOE OSTI2019

Assessing the Quality of Molecular Simulations for Vapor–Liquid Equilibria: An Analysis of the TraPPE Database

As molecular modeling and simulation techniques become increasingly important sources of thermophysical property and phase equilibrium data, the ability to assess the robustness of that data becomes more critical. Recently, the use of the compressibility factor (Z) has been suggested as a metric for testing the quality of simulation data for vapor–liquid equilibria (VLE). Int his paper, we analyze predicted VLE data from the transferable potentials for phase equilibria (TraPPE) database and show that, apart from data entry or typographical errors, Z will always be well-behaved in Gibbs ensemble Monte Carlo (GEMC) simulations even when the simulations are not sufficiently equilibrated. However, this is not true for grand canonical Monte Carlo simulations. When the pressure is calculated from the internal forces, then pressure and density are strongly correlated for the vapor phase and, for GEMC simulations, it is recommended to treat Z as an instantaneous mechanical property. From analysis of the TraPPE VLE data, we propose a complementary metric based on the predicted vapor pressures at three neighboring temperatures and their deviation from a local Clausius–Clapeyron fit.

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

Molecular Simulations in Astrobiology

One of the main goals of astrobiology is to understand the origin of cellular life. The most direct approach to this problem is to construct laboratory models of protocells. Such efforts, currently underway in the NASA Astrobiology Program, are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing energy and responding to changes in the environment, are carried out by proteins bound to membranes. We use computer simulations to address the following questions about these proteins: (1) How do small proteins self-organize into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides form membrane-spanning structures (e.g. channels)? (3) By what mechanisms do such structures perform their functions? The simulations are performed using the molecular dynamics method. In this method, Newton's equations of motion for each atom in the system are solved iteratively. At each time step, the forces exerted on each atom by the remaining atoms are evaluated by dividing them into two parts. Short-range forces are calculated in real space while long-range forces are evaluated in reciprocal space, using a particle-mesh algorithm which is of order O(NInN). With a time step of 2 femtoseconds, problems occurring on multi-nanosecond time scales (10(exp 6)-10(exp 8) time steps) are accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude, which requires algorithmic improvements and codes scalable to a large number of processors. Work in this direction is in progress. Two series of simulations are discussed. In one series, it is shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and fold into helical structures (see Figure). Once in the membrane, the peptides exhibit orientational flexibility with changing conditions, which may have provided a mechanism of transmitting signals between the protocell and its environment. In another series of simulations, the mechanism by which a simple protein channel efficiently mediates proton transport across membranes was investigated. This process is a key step in cellular bioenergetics. In the channel under study, proton transport is gated by four histidines that occlude the channel pore. The simulations identify the mechanisms by which protons move through the gate.

Pohorille, Andrew↗
NASA NTRS2001

Molecular Simulations in Astrobiology

One of the main goals of astrobiology is to understand the origin of cellular life. In the absence of any record of the earliest ancestors of contemporary cells, protocells, the most direct way to test our understanding of their characteristics is to construct laboratory models of protocells. Such efforts, currently underway in the NASA Astrobiology Program, are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures and developing designs of molecules that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing and storing energy, and responding to changes in the environment, are carried out by proteins bound to membranes. We use computer simulations to address the following, questions about these proteins: (1) How do small proteins (peptides) organize themselves into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides aggregate to form membrane-spannin(y structures (e.g., channels)? (3) By what mechanisms do such aggregates perform their functions? The simulations are performed using the molecular dynamics (MD) method. In this method, Newton's equations of motion for each atom in the system are solved iteratively. At each time step, the forces exerted on each atom by the remaining atoms are evaluated by dividing them into two parts. Short-range forces are calculated directly in real space while long-range forces are evaluated in reciprocal space, usually using a particle-mesh algorithm which is of order O(NlnN). Currently, a time step of 2 femtoseconds is typically used, thereby making studies of problems occurring on multi-nanosecond time scales (10(exp 6) - 10(exp 8) time steps) accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude. Such an extension requires both algorithmic improvements and codes scalable to a large number of parallel processors. Work in this direction is in progress. Two specific series of simulations that demonstrate how peptides self-organize and function in membranes are discussed. In one series of simulations, it was shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and, simultaneously, fold into two different helical structures, which remain in equilibrium. Once in the membrane, the peptides can readily change their orientation, especially in response to local electric fields. This structural and orientational flexibility of peptides with changing conditions may have provided a mechanism of transmitting signals between the environment and the interior of the protocell. In another series of simulations, the mechanism by which a simple protein channel efficiently mediates proton transport across membranes was investigated. This process is a key step in cellular bioenergetics. In the channel under study, proton transport is gated by four histidines that occlude the channel pore. The simulations demonstrate that protons move through the gate by a "shuttle" mechanism, wherein one histidine is protonated on the extracellular side and, subsequently, the proton bound on the opposite side is released.

Pohorille, Andrew↗
DOE OSTI2019

Antimicrobial Zn-Based “TSOL” for Citrus Greening Management: Insights from Spectroscopy and Molecular Simulation

Huanglongbing (HLB), also known as citrus greening, is a bacterial disease that poses a devastating threat to the citrus industry worldwide. To manage this disease efficiently, we developed and characterized a ternary aqueous solution (TSOL) that contains zinc nitrate, urea, and hydrogen peroxide. In this paper, we report that TSOL exhibits better antimicrobial activity than commercial bactericides for growers. X-ray fluorescence analysis demonstrates that zinc is delivered to citrus leaves, where the bacteria reside. FTIR and Raman spectroscopy, molecular dynamics simulations, and density functional theory calculations elucidate the solution structure of TSOL and reveal a water-mediated interaction between Zn 2+ and H 2 O 2 , which may facilitate the generation of highly reactive hydroxyl radicals contributing to superior antimicrobial activity of TSOL. Lastly, our results not only suggest TSOL as a potent antimicrobial agent to suppress bacterial growth in HLB-infected trees, but also provide a structure–property relationship that explains the superior performance of TSOL.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance

For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λ max value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λ max values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λ max of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λ max using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λ max values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λ max of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λ max values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

The alloying element dependence of the local lattice deformation and the elastic properties of Ni{sub 3}Al: A molecular dynamics simulation

Molecular dynamics (MD) together with the modified analytical embedded atom method (MAEAM) is employed to study the alloying elements (Re, Ru, Co, and Ta) dependence of the elastic properties of L1{sub 2}-Ni{sub 3}Al. The investigations indicate that the calculated elastic properties of Ni{sub 3}Al are in reasonable agreement with the previous results. The substituting formation energies of the alloying elements in Ni{sub 3}Al are calculated to determine the site preference. It is found out that Re, Ru, and Ta atoms prefer to occupy the Al sites, and the Co atom prefers to occupy the Ni site. Based on Re, Ru, and Ta substituting the 1st, 2nd, 3rd, and 4th nearest-neighbor atoms, we ascertain that the substituting manners of these alloying elements have a decisive effect on the bulk modulus and the local crystal lattice of Ni{sub 3}Al. Moreover, for Re, the bonding interaction plays a predominant role in the improvement in the bulk modulus of Ni{sub 3}Al, whereas the size effect of Ru and Ta on the improvement in the bulk modulus is more obvious. Finally, the potential reasons of which the alloying elements enhance the bulk modulus are discussed in details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗