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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Reduced methane recovery at high pressure due to methane trapping in shale nanopores

By 2050, shale gas production is expected to exceed three-quarters of total US natural gas production. However, current unconventional hydrocarbon gas recovery rates are only around 20%. Maximizing production of this natural resource thus necessitates improved understanding of the fundamental mechanisms underlying hydrocarbon retention within the nanoporous shale matrix. In this study, we integrated molecular simulation with high-pressure small-angle neutron scattering (SANS), an experimental technique uniquely capable of characterizing methane behavior in situ within shale nanopores at elevated pressures. Samples were created using Marcellus shale, a gas-generative formation comprising the largest natural gas field in the United States. Our results demonstrate that, contrary to the conventional wisdom that elevated drawdown pressure increases methane recovery, a higher peak pressure led to the trapping of dense, liquid-like methane in sub-2 nm radius nanopores, which comprise more than 90% of the measured nanopore volume, due to irreversible deformation of the kerogen matrix. These findings have critical implications for pressure management strategies to maximize hydrocarbon recovery, as well as broad implications for fluid behavior under confinement.

58 GEOSCIENCES↗
DOE OSTI2023

Crystallization of binary nanocrystal superlattices and the relevance of short-range attraction

The synthesis of binary nanocrystal superlattices (BNSLs) enables the targeted integration of orthogonal physical properties, such as photoluminescence and magnetism, into a single superstructure, unlocking a vast design space for multifunctional materials. However, the formation mechanism of BNSLs remains poorly understood, restricting the prediction of the structure and properties of superlattices. Here we use a combination of in situ scattering and molecular simulation to elucidate the self-assembly of two common BNSLs (AlB 2 and NaZn 13 ) through emulsion templating. Our self-assembly experiments reveal that no intermediate structures precede the formation of the final binary phases, indicating that their formation proceeds through classical nucleation. Using simulations, we find that, despite the formation of AlB 2 and NaZn 13 typically being attributed to entropy, their self-assembly is most consistent with the nanocrystals possessing short-range interparticle attraction, which we find can accelerate nucleation kinetics in BNSLs. Here, we also find homogeneous, classical nucleation in simulations, corroborating our experiments. These results establish a robust correspondence between experiment and theory, paving the way towards prediction of BNSLs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Alkyl decorated metal–organic frameworks for selective trapping of ethane from ethylene above ambient pressures

Here the trapping of paraffins is beneficial compared to selective olefin adsorption for adsorptive olefin purification from a process engineering point of view. Here we demonstrate the use of a series of Zn 2 (X-bdc) 2 (dabco) (where X-bdc 2- is bdc 2- = 1,4-benzenedicarboxylate with substituting groups X, DM-bdc 2- = 2,5-dimethyl-1,4-benzenedicarboxylate or TM-bdc 2- = 2,3,5,6-tetramethyl-1,4-benzenedicarboxylate and dabco = diazabicyclo[2.2.2.]octane) metal-organic frameworks (MOFs) for the adsorptive removal of ethane from ethylene streams. The best performing material from this series is Zn2(TM-bdc)2(dabco) (DMOF-TM), which shows a high ethane uptake of 5.31 mmol g -1 at 110 kPa, with a good IAST selectivity of 1.88 towards ethane over ethylene. Through breakthrough measurements a high productivity of 13.1 L/kg per breakthrough is revealed with good reproducibility over five consecutive cycles. Molecular simulations show that the methyl groups of DMOF-TM are forming a van der Waals trap with the methylene groups from dabco, snuggly fitting the ethane. Further, rarely used high pressure coadsorption measurements, in pressure regimes that most scientific studies on hydrocarbon separation on MOFs ignore, reveal an increase in ethane capacity and selectivity for binary mixtures with increased pressures. The coadsorption measurements reveal good selectivity of 1.96 at 1000 kPa, which is verified also through IAST calculations up to 3000 kPa. This study overall showcases the opportunities that pore engineering by alkyl group incorporation and pressure increase offer to improve hydrocarbon separation in reticular materials.

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

Efficient separation of xylene isomers by using a robust calcium-based metal–organic framework through a synergetic thermodynamically and kinetically controlled mechanism

Adsorptive separation of physically and chemically similar molecules and understanding the underlying host–guest interactions at the molecular level are of significant scientific and practical importance. Here we report the development of a novel calcium-based metal–organic framework, formulated as Ca 3 (Htcpp) 2 (H 4 tcpp = 2,3,5,6-tetrakis(4-carboxyphenyl)-pyrazine) featuring microporosity and high stability. This compound shows distinct adsorption behavior toward xylene isomers and is capable of separating them efficiently at an industrially relevant temperature. The selective adsorption is attributed to a synergetic thermodynamic and kinetic effect. Here, the host–guest interactions were probed directly by single-crystal X-ray diffraction analysis and the adsorption affinity was evaluated through computational molecular simulations.

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

Deep learning potential of mean force between polymer grafted nanoparticles

Grafting polymer chains on the surfaces of nanoparticles is a well-known route to control their self-assembly and distribution in a polymer matrix. A wide variety of self-assembled structures are achieved by changing the grafting patterns on the surface of an individual nanoparticle. However, an accurate estimation of the effective potential of mean force between a pair of grafted nanoparticles that determines their assembly and distribution in a polymer matrix is an outstanding challenge in nanoscience. We address this problem via deep learning. As a proof of concept, here we report a deep learning framework that learns the interaction between a pair of single-chain grafted spherical nanoparticles from their molecular dynamics trajectory. Subsequently, we carry out the deep learning potential of mean force-based molecular simulation that predicts the self-assembly of a large number of single-chain grafted nanoparticles into various anisotropic superstructures, including percolating networks and bilayers depending on the nanoparticle concentration in three-dimensions. The deep learning potential of mean force-predicted self-assembled superstructures are consistent with the actual superstructures of single-chain polymer grafted spherical nanoparticles. This deep learning framework is very generic and extensible to more complex systems including multiple-chain grafted nanoparticles. Finally, we expect that this deep learning approach will accelerate the characterization and prediction of the self-assembly and phase behaviour of polymer-grafted and unfunctionalized nanoparticles in free space or a polymer matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Curvature-controlled geometrical lensing behavior in self-propelled colloidal particle systems

In many biological systems, the curvature of the surfaces cells live on influences their collective properties. Curvature should likewise influence the behavior of active colloidal particles. We show using molecular simulation of self-propelled active particles on surfaces of Gaussian curvature (both positive and negative) how curvature sign and magnitude can alter the system's collective behavior. Curvature acts as a geometrical lens and shifts the critical density of motility-induced phase separation (MIPS) to lower values for positive curvature and higher values for negative curvature, which we explain theoretically by the nature of parallel lines in spherical and hyperbolic space. Curvature also fluidizes dense MIPS clusters due to the emergence of defect patterns disrupting the crystalline order inside the clusters. Using our findings, we engineer three confining surfaces that strategically combine regions of different curvature to produce a host of novel dynamical behaviors, including cyclic MIPS on spherocylinders, directionally biased cyclic MIPS on spherocones, and position dependent cluster fluctuations on metaballs.

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

Confinement-induced clustering of H 2 and CO 2 gas molecules in hydrated nanopores

Gas molecule clustering within nanopores holds significance in the fields of nanofluidics, biology, gas adsorption/desorption, and geological gas storage. However, the intricate roles of nanoconfinement and surface chemistry that govern the formation of gas clusters remain inadequately explored. In this study, through free energy calculation in molecular simulations, we systematically compared the tendencies of H 2 and CO 2 molecules to aggregate within hydrated hydrophobic pyrophyllite and hydrophilic gibbsite nanopores. The results indicate that nanoconfinement enhances gas dimer formation in the nanopores, irrespective of surface chemistry. However, surface hydrophilicity prohibits the formation of gas clusters larger than dimers, while large gas clusters form easily in hydrophobic nanopores. Despite H 2 and CO 2 both being non-polar, the larger quadrupole moment of CO 2 leads to a stronger preference for dimer/cluster formation compared to H 2 . Here, our results also indicate that gases prefer to enter the nanopores as individual molecules, but exit the nanopores as dimers/clusters. This investigation provides a mechanistic understanding of gas cluster formation within nanopores, which is relevant to various applications, including geological gas storage.

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

Iodide oxidation by ozone at the surface of aqueous microdroplets

The kinetics of iodide oxidation by ozone are studied in single microdroplets to explore chemical reactivity at the air–water interface. Molecular simulations and kinetic modeling are employed to understand the underlying multiphase mechanism.

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

Low hydrogen solubility in clay interlayers limits gas loss in hydrogen geological storage

Gas intercalation into clay interlayers may result in hydrogen loss in the geological storage of hydrogen; a phenomenon that has not been fully understood and quantified. Here we use metadynamics molecular simulations to calculate the free energy landscape of H 2 intercalation into montmorillonite interlayers and the H 2 solubility in the confined water; in comparison with results obtained for CO 2 . Furthermore, the results indicate that H 2 intercalation into hydrated interlayers is thermodynamically unfavorable while CO 2 intercalation can be favorable. H 2 solubility in hydrated clay interlayers is in the same order of magnitude as that in bulk water and therefore no over-solubility effect due to nanoconfinement is observed – in striking contrast with CO 2 . These results indicate that H 2 loss and leakage through hydrated interlayers due to intercalation in a subsurface storage system, if any, is limited.

58 GEOSCIENCES↗
DOE OSTI2024

Impact of surface hydrophilicity on the ordering and transport properties of bicontinuous microemulsions

Microemulsions (MEs) have many industrial applications, where recent developments have shown that MEs can be utilized for electrochemical applications, including potentially in redox flow batteries. However, understanding the structure and dynamics of these systems, including at a surface, is needed to direct and rationally control their electrochemical behavior. While bulk solution measurements have provided insight into their structure, their assembly at an interface also impacts the electron (to the electrode) and ion (across the surfactant) charge transfer processes in the system. To address this shortcoming, neutron reflectivity experiments and molecular simulations have been completed that document the near surface structure of a series of deuterated water (D 2 O)/polysorbate-20/toluene MEs on hydrophilic and amphiphilic surfaces. These results show that the microemulsions form complex layered structures near a hard electrode surface, where most layers are not purely one component. Decreasing the D 2 O concentration in the ME increases the number of and purity of the layers established on the solid surface. These lamellar-type layers transition from the surface to the bulk microemulsion as a series of mixed layers (i.e., containing oil, water, and surfactant) that are consistent with perforated lamellae. Additionally, these mixed lamellae appear to become more perforated with oil and water pathways on an amphiphilic surface. Furthermore, the purity and thickness of these layers will influence the accessibility of an electrode by redox active species, as well as ion transport required to satisfy the electroneutrality condition.

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

Near-cryogenic direct air capture using adsorbents

Direct air capture (DAC) of CO 2 is a key component in the portfolio of negative emissions technologies for mitigating global warming. However, even with the most potent amine sorbents, large-scale DAC deployment remains limited by high energy and capital costs. Recently, adsorbents relying on weak interactions with CO 2 have emerged as a potential alternative, thanks to their rapid adsorption kinetics and superior long-term stability, particularly under sub-ambient conditions (∼253 K). Despite these advantages, their use is hindered by the need for a water-removal process, location-specific constraints, and insufficient working capacity even in cold climates. In this study, we hypothesized that further reducing the adsorption temperature to a near-cryogenic range (160–220 K) could enable cost-effective DAC by utilizing the full potential of physisorbents. We primarily consider integrating DAC with a relatively untapped source of cold energy—liquified natural gas (LNG) regasification—to perform near-cryogenic DAC. From large-scale molecular simulations, Zeolite 13X and CALF-20 were identified as promising candidates. These materials were subsequently examined through experiments, including breakthrough analyses at 195 K. Their high CO 2 sorption capacity (4.5–5.5 mmol g −1 ), combined with a low desorption enthalpy and robust long-term stability, led to a threefold reduction in the levelized cost of capture (down to 68.2 USD per tonne CO 2 ). Estimates of the global LNG regasification resource suggest that LNG–DAC coupling could potentially enable the capture of 103–142 megatonnes of CO 2 annually as of 2050.

Kim, Seo-Yul [Georgia Institute of Technology, Atl↗
DOE OSTI2021

NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields

The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar–600 bar and 450 K–750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.

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

A Gaussian field approach to the planar electric double layer structures in electrolyte solutions

Here, in this work, the planar, electric, double-layer structures of non-polarizable electrodes in electrolyte solutions are studied with Gaussian field theory. A response function with two Yukawa functions is used to capture the electrostatic response of the electrolyte solution, from which the modified response function in the planar symmetry is derived analytically. The modified response function is further used to evaluate the induced charge density and the electrostatic potential near an electrode. The Gaussian field theory, combined with a two-Yukawa response function, can reproduce the oscillatory decay behavior of the electric potentials in concentrated electrolyte solutions. When the exact sum rules for the bulk electrolyte solutions and the electric double layers are used as constraints to determine the parameters of the response function, the Gaussian field theory could at least partly capture the nonlinear response effect of the surface charge density. Comparison with results for a planar electrode with fixed surface charge densities from molecular simulations demonstrates the validity of Gaussian field theory.

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

Extended molecular eigenmodes treatment of dipole–dipole NMR relaxation in real fluids

Traditional models of NMR relaxation fail to account for the complex, multi-exponential behavior of the autocorrelation function in realistic systems characterized by soft-interactions and molecules that are chemically and physically complex. Here, in this study, we describe the relative diffusion of the spin dipoles by means of a Fokker–Planck equation that includes an interaction potential of mean force to account for the response of the physical/chemical environment around the dipoles. By numerically solving the Fokker–Planck equation for the diffusion propagator, we estimate dipole–dipole NMR relaxation for like- and unlike-spin systems via its eigenmode solution. We test the model against molecular simulations of diffusing dipoles with harmonic potentials and also validate using experimental longitudinal relaxation data from real systems, including Gd(III)–aqua and Gd(III)–DO3A–butrol complexes, the latter being an important MRI contrast agent. Using this novel approach, we predict both the inner- and outer-shell contributions to the relaxivity rates with excellent accuracy at frequencies relevant to MRI. We also show that, under the appropriate assumptions, our framework naturally recovers the Bloembergen–Purcell–Pound, the Solomon–Bloembergen–Morgan, and the Hwang–Freed models. Our implementation is general and publicly available for application to a broad range of systems.

Pinheiro dos Santos, Thiago J. [Rice Univ., Housto↗
DOE OSTI2020

TDP-43 α-helical structure tunes liquid–liquid phase separation and function

Liquid–liquid phase separation (LLPS) is involved in the formation of membraneless organelles (MLOs) associated with RNA processing. The RNA-binding protein TDP-43 is present in several MLOs, undergoes LLPS, and has been linked to the pathogenesis of amyotrophic lateral sclerosis (ALS). While some ALS-associated mutations in TDP-43 disrupt self-interaction and function, here we show that designed single mutations can enhance TDP-43 assembly and function via modulating helical structure. Using molecular simulation and NMR spectroscopy, we observe large structural changes upon dimerization of TDP-43. Two conserved glycine residues (G335 and G338) are potent inhibitors of helical extension and helix–helix interaction, which are removed in part by variants at these positions, including the ALS-associated G335D. Substitution to helix-enhancing alanine at either of these positions dramatically enhances phase separation in vitro and decreases fluidity of phase-separated TDP-43 reporter compartments in cells. Furthermore, G335A increases TDP-43 splicing function in a minigene assay. Therefore, the TDP-43 helical region serves as a short but uniquely tunable module where application of biophysical principles can precisely control assembly and function in cellular and synthetic biology applications of LLPS.

60 APPLIED LIFE SCIENCES↗