Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular simulations”

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 37 records · Page 2

[Using Molecular Simulations to Understand Complex Nanoscale Dynamic Phenomena in Polymer Solutions]

The first half of the project concentrated on molecular simulation studies of the translocation of model molecules for single-stranded DNA through a nanosized pore. This has resulted in the publication, Translocation of a polymer chain across a nanopore: A Brownian dynamics simulation study, by Pu Tian and Grant D. Smith, JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 21 1 DECEMBER 2003, which is attached to this report. In this work we carried out Brownian dynamics simulation studies of the translocation of single polymer chains across a nanosized pore under the driving of an applied field (chemical potential gradient) designed to mimic an electrostatic field. The translocation process can be either dominated by the entropic barrier resulted from restricted motion of flexible polymer chains or by applied forces (or chemical gradient). We focused on the latter case in our studies. Calculation of radius of gyration of the translocating chain at the two opposite sides of the wall shows that the polymer chains are not in equilibrium during the translocation process. Despite this fact, our results show that the one-dimensional diffusion and the nucleation model provide an excellent description of the dependence of average translocation time on the chemical potential gradients, the polymer chain length and the solvent viscosity. In good agreement with experimental results and theoretical predictions, the translocation time distribution of our simple model shows strong non-Gaussian characteristics. It is observed that even for this simple tube-like pore geometry, more than one peak of translocation time distribution can be generated for proper pore diameter and applied field strengths. Both repulsive Weeks-Chandler-Anderson and attractive Lennard-Jones polymer-nanopore interaction were studied. Attraction facilitates the translocation process by shortening the total translocation time and dramatically improve the capturing of polymer chain. The width of the translocation time distribution was found to decrease with increasing temperature, increasing field strength, and decreasing pore diameter.

Smith, Grant↗

A comparison of probabilistic generative frameworks for molecular simulations

Generative artificial intelligence is now a widely used tool in molecular science. Despite the popularity of probabilistic generative models, numerical experiments benchmarking their performance on molecular data are lacking. Here, in this work, we introduce and explain several classes of generative models, broadly sorted into two categories: flow-based models and diffusion models. We select three representative models: neural spline flows, conditional flow matching, and denoising diffusion probabilistic models, and examine their accuracy, computational cost, and generation speed across datasets with tunable dimensionality, complexity, and modal asymmetry. Our findings are varied, with no one framework being the best for all purposes. In a nutshell, (i) neural spline flows do best at capturing mode asymmetry present in low-dimensional data, (ii) conditional flow matching outperforms other models for high-dimensional data with low complexity, and (iii) denoising diffusion probabilistic models appear the best for low-dimensional data with high complexity. Our datasets include a Gaussian mixture model and the dihedral torsion angle distribution of the Aib9 peptide, generated via a molecular dynamics simulation. We hope our taxonomy of probabilistic generative frameworks and numerical results may guide model selection for a wide range of molecular tasks.

Artificial intelligence↗

Molecular simulation of lignin-related aromatic compound permeation through gram-negative bacterial outer membranes

Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Effect of Water on the Work of Adhesion at Epoxy Interfaces by Molecular Dynamics Simulation

Molecular dynamics simulation can be used to explore the detailed effects of chemistry on properties of materials. In this paper, two different epoxies found in aerospace resins are modeled using molecular dynamics. The first material, an amine-cured tetrafunctional epoxy, represents a composite matrix resin, while the second represents a 177 C-cured adhesive. Surface energies are derived for both epoxies and the work of adhesion values calculated for the epoxy/epoxy interfaces agree with experiment. Adding water -- to simulate the effect of moisture exposure -- reduced the work of adhesion in one case, and increased it in the other. To explore the difference, the various energy terms that make up the net work of adhesion were compared and the location of the added water was examined.

Hinkley, J.A.↗

Molecular Simulations of the Chain Length Dependent Adsorption of C7‐C14 n ‐Alkanes in ZIF‐8

Recent experiments show that the diffusivities of C8 and C10 n-alkanes in ZIF-8 are higher than those of C7 and C9, respectively. We investigated this unusual ‘odd-even' effect by simulating the adsorption of C7-C14 n-alkanes in ZIF-8 using hybrid Monte Carlo molecular simulations. The resultant adsorption isotherms, guest-host energies, isosteric heats, and chain length distributions are analyzed for trends among the n-alkanes. ZIF-8 cages filled with n-alkanes are characterized using a combination of image processing and data science to quantify the differences in packing that occur with chain length. In conclusion, results indicate that packing changes drastically from C7 to C8 and from C9 to C10, which is consistent with the diffusion trends previously reported in experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Revealing a Hidden Intermediate of Rotatory Catalysis with X-ray Crystallography and Molecular Simulations

The mechanism of rotatory catalysis in ATPhydrolyzing molecular motors remains an unresolved puzzle in biological energy transfer. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, knowledge on how the coupling between the chemical and mechanical steps within motors enforces directional rotatory movements remains fragmentary. Even more contentious is to pinpoint the rate-limiting step of a multistep rotation process. Here, using vacuolar or V 1 -type hexameric ATPase as an exemplary rotational motor, we present a model of the complete 4-step conformational cycle involved in rotatory catalysis. First, using Xray crystallography, a new intermediate or “dwell” is identified, which enables the release of an inorganic phosphate (or P i ) after ATP hydrolysis. Using molecular dynamics simulations, this new dwell is placed in a sequence with three other crystal structures to derive a putative cyclic rotation path. Free-energy simulations are employed to estimate the rate of the hexameric protein transformations and delineate allosteric effects that allow new reactant ATP entry only after hydrolysis product exit. An analysis of transfer entropy brings to light how the side-chain-level interactions transcend into larger-scale reorganizations, highlighting the role of the ubiquitous arginine-finger residues in coupling chemical and mechanical information. An inspection of all known rates encompassing the 4-step rotation mechanism implicates the overcoming of the ADP interactions with V 1 -ATPase to be the rate-limiting step of motor action.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer ice-philicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface–ice contact at the expense of surface–water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Furthermore, our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

LongBondEliminator: A Molecular Simulation Tool to Remove Ring Penetrations in Biomolecular Simulation Systems

We develop a workflow, implemented as a plugin to the molecular visualization program VMD, that can fix ring penetrations with minimal user input. LongBondEliminator, detects ring piercing artifacts by the long, strained bonds that are the local minimum energy conformation during minimization for some assembled simulation system. The LongBondEliminator tool then automatically treats regions near these long bonds using multiple biases applied through NAMD. By combining biases implemented through the collective variables module, density-based forces, and alchemical techniques in NAMD, LongBondEliminator will iteratively alleviate long bonds found within molecular simulation systems. Through three concrete examples with increasing complexity, a lignin polymer, an viral capsid assembly, and a large, highly glycosylated protein aggrecan, we demonstrate the utility for this method in eliminating ring penetrations from classical MD simulation systems. The tool is available via gitlab as a VMD plugin, and has been developed to be generically useful across a variety of biomolecular simulations.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

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↗

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↗

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↗

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↗

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↗

Investigating Lignin Aggregation and Interactions with Solvents during γ-Valerolactone (GVL) Pretreatment: A Combined Small Angle Neutron Scattering and Molecular Simulations Study

The strong tendency of lignin to aggregate in solution, coupled with limited understanding of how its molecular structure governs this behavior, hinders its effective utilization in biorefineries. Here, in this study, we investigated the solution behavior of lignin extracted from poplar using γ-valerolactone/water (GVL/H2O, 9:1 wt/wt) through combined small-angle neutron scattering (SANS) and molecular dynamics (MD) simulations. Lignin samples obtained at 100 °C (L100) and 120 °C (L120) differed in β–O–4 content, hydroxyl distribution, and S/G ratio, enabling direct assessment of how molecular composition governs solvation and aggregation. SANS showed that L120 formed rigid and elongated cylindrical aggregates at 25 °C that transitioned to more flexible spheroidal structures by 50 °C and remained stable up to 80 °C, whereas L100 adopted globular aggregates that progressively collapsed with increasing temperature. MD simulations reinforced these observations by showing that S-rich (L120-like) oligomers had larger radii of gyration, stronger solvent coordination driven by methoxy groups, and fewer lignin–lignin contacts. In contrast, G-rich (L100-like) oligomers displayed persistent aggregation and lower solubility. Collectively, these results reveal that increased aromatic methoxylation enhances lignin–solvent interactions and suppresses self-association, whereas reduced methoxylation and higher β–O–4 content promote persistent aggregation into colloid-like structures with restricted solvent penetration into the aggregate interior.

biorefinery pretreatment↗

Sulfonated polybenzimidazole membrane with graphene oxide additive for 2,3-butanediol/water separation: A molecular simulation

Membrane separation for 2,3-butanediol (2,3-BDO) recovery from fermentation broth is highly valued for sustainable and renewable processes, but it requires efficient membrane materials. Here, this work evaluates the sulfonated polybenzimidazole (sPBI) and its graphene oxide (GO) doped composite membrane for separating 2,3-BDO and water via atomistic simulations. Density functional theory calculations are applied to identify various forms of sPBI structures and quantify their binding interactions with 2,3-BDO and water. Classical molecular dynamic simulations are used to evaluate the structural changes, diffusivity, and selectivity of 2,3-BDO and water in different sPBI models, GO surfaces, and GO-doped sPBI composite models. Our results suggest that sPBI slightly increases the crystallinity of the membrane structures, enhances the adsorption strength for both 2,3-BDO and water, and improves the water/2,3-BDO selectivity by 2–3 times. The GO surfaces display a maximum selectivity at a surface coverage of 0.1–0.15 for both hydroxyl and epoxy surface groups. The addition of GO flakes to sPBI creates new interaction sites for 2,3-BDO and water at the interface of sPBI and GO, and the water/2,3-BDO selectivity of GO-doped sPBI models is further increased up to 3 times. This work illustrates how the integrated addition of sPBI and GO flakes offers a promising approach to selective separation of 2,3-BDO and water, providing theoretical guidance for polybenzimidazole-based membranes in the potential application of 2,3-BDO recovery.

2,3-butanediol↗

Molecular simulation of the structural and thermodynamic properties of n-alkane/brine interfacial systems with nonionic surfactants

All-atom molecular dynamics (MD) simulations are used to study the structural and thermodynamic properties of water + n-heptane interfacial systems in the presence of salinity (NaCl, NaI, KCl, and KI) and two different nonionic surfactants. Excellent qualitative and quantitative agreement with interfacial tension experimental data is obtained, but this requires applying a charge scaling factor to the ionic species, suggesting the likely role of polarizability effects. Tensoactive behavior is observed for the NaI and KI systems, coinciding with layering of the iodide species near the interface. Further, there are significant interaction differences between the surfactants and the different ionic species, but the structural and thermodynamic behavior of the surfactant molecules is relatively unaffected by the specific salt species. Our analysis suggest that the surfactant-surfactant interactions play a critical role in determining the interfacial behavior, with a much smaller impact associated with the composition of the brine phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗