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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Interfacial Ice Density Fluctuations Inform Surface Ice-Philicity

The propensity of a surface to nucleate ice or bind to ice is governed by its ice-philicity─its relative preference for ice over liquid water. However, the relationship between the features of a surface and its ice-philicity is not well understood, and for surfaces with chemical or topographical heterogeneity, such as proteins, their ice-philicity is not even well-defined. In the analogous problem of surface hydrophobicity, it has been shown that hydrophobic surfaces display enhanced low water-density (vapor-like) fluctuations in their vicinity. To interrogate whether enhanced ice-like fluctuations are similarly observed near ice-philic surfaces, here we use molecular simulations and enhanced sampling techniques. Using a family of model surfaces for which the wetting coefficient, k , has previously been characterized, we show that the free energy of observing rare interfacial ice-density fluctuations decreases monotonically with increasing k . By utilizing this connection, we investigate a set of fcc systems and find that the (110) surface is more ice-philic than the (111) or (100) surfaces. By additionally analyzing the structure of interfacial ice, we find that all surfaces prefer to bind to the basal plane of ice, and the topographical complementarity of the (110) surface to the basal plane explains its higher ice-philicity. Using enhanced interfacial ice-like fluctuations as a measure of surface ice-philicity, we then characterize the ice-philicity of chemically heterogeneous and topologically complex systems. In particular, we study the spruce budworm antifreeze protein (sbwAFP), which binds to ice using a known ice-binding site (IBS) and resists engulfment using nonbinding sites of the protein (NBSs). We find that the IBS displays enhanced interfacial ice-density fluctuations and is therefore more ice-philic than the two NBSs studied. We also find the two NBSs are similarly ice-phobic. By establishing a connection between interfacial ice-like fluctuations and surface ice-philicity, our findings thus provide a way to characterize the ice-philicity of heterogeneous surfaces.

crystals↗
DOE OSTI2020

High Pressure Infiltration/Expulsion of Aqueous NaCl in Planar Hydrophobic Nanopores

Pressure-driven permeation of water in a poorly wettable material results in a conversion of mechanical work into surface free energy representing a new form of energy storage, or absorption. When water is replaced by a concentrated electrolyte solution, the storage capacity of a nanoporous medium becomes comparable to high-end supercapacitors. The addition of salt can also reduce the hysteresis of the infiltration/expulsion cycle. Our molecular simulations provide a theoretical perspective into the mechanisms involved in the process, and underlying structures and interactions in compressed nanoconfined solutions. Specifically, we consider aqueous NaCl in planar confinements of widths of 1.0 nm and 1.64 nm and pressures of up to 3 kbar. Open ensemble Monte Carlo simulations utilizing fractional exchanges of molecules for efficient additions/removal of ions have been utilized in conjunction with pressure-dependent chemical potentials to model bulk phases under pressure. Confinements open to these pressurized bulk, aqueous electrolyte phases show reversibility at narrow pore sizes and irreversibility in wider ones, consistent with experiment, as well as strong hysteresis at both pore size. The addition of salt results in significant increases in the solid/liquid interfacial tension in narrower pores and associated infiltration and expulsion pressures. Furthermore, these changes are consistent with strong desalination effects at the lower pore size observed irrespective of external pressure and initial concentration.

25 ENERGY STORAGE↗
DOE OSTI2022

Depletion-Driven Assembly of Polymer-Coated Nanocrystals

Depletion-driven assembly has been widely studied for micron-sized colloids, but questions remain at the nanoscale where the governing physics are impacted by stabilizing surface ligands or wrapping polymers, whose dimensions are similar to those of the colloidal core and the depletant. Here, we probe how wrapping colloidal tin-doped indium oxide nanocrystals with polymers affects their depletion-induced interactions and assembly in solutions of polyethylene glycol. Copolymers of polyacrylic acid grafted with polyethylene oxide provide nanocrystal wrappings with different effective polymer graft densities and molecular weights. (Ultra) small angle X-ray scattering, coarse-grained molecular dynamics simulation, and molecular thermodynamic theory were combined to analyze how depletant size and polymer wrapping characteristics affect depletion interactions, structure, and phase behavior. The results show how depletant molecular weight, as well as surface density and molecular weight of polymer grafts, set thresholds for assembly. These signatures are unique to depletion-driven assembly of nanoscale colloids and mirror phase behaviors of grafted nanoparticle–polymer composites. Optical and rheological responses of depletion-driven assemblies of nanocrystals with different polymer shell architectures were probed to learn how their structural differences impact properties. Finally, we discuss how these wrinkles for depletion-driven assembly at the nanoscale may provide fresh opportunities for designing responsive depletion interactions and dynamically reconfigurable materials.

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

Post-Hydration Crosslinking of Ion Exchange Membranes to Control Water Content

The utilization of crosslinking as a technique to limit excess water uptake in ion exchange membranes is a well-established approach. However, traditional methods involve performing the crosslinking reaction before the membrane self-segregates into hydrophilic and hydrophobic domains upon exposure to water. In this study, we employ numerical and molecular simulations to explore the potential of post-hydration crosslinking of the polymer at subsaturated conditions. This process allows for the development of a nanosegregated morphology, which establishes well-connected but narrow water channels, while limiting additional sorption of water. The significant reduction in membrane swelling observed in our study indicates that post-hydration crosslinking could be a promising approach to regulate water uptake in ion exchange membranes. Here based on our analyses, we propose chemical crosslinking strategies that may yield high membrane gel fractions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Degradation of Anion Exchange Membranes by Cation Elimination: Impact on Water Uptake, Nanostructure, and Ionic Mobility

Anion exchange membranes (AEMs) are an attractive platform for fuel cell and electrolysis technologies since they enable the use of cheaper, nonplatinum group metal electrodes and components. However, the widespread adoption of AEM-based devices is limited by chemical degradation of the AEM in the highly alkaline medium of operation. Experimental studies report three major pathways of degradation, including substitution (S N 2), elimination (E2), and backbone cleavage. The decline in membrane performance is likely due to a cumulative effect of several pathways, which has proven to be difficult to disentangle through experiments. Here we use coarse-grained molecular simulations to isolate the impact of E2 degradation, where cationic sites are replaced by an alkene group, on AEM water uptake, nanoscale morphology, structure of the ion-conducting channels, water mobility, and ionic conductivity. Our studies focus on the well-studied model AEM polyphenylene oxide with tetraalkylammonium cationic sites (PPO-TMA). We find that about half of the hydrophobic groups resulting from the degradation retract into the hydrophobic domains of the membrane. The reduction in ion exchange capacity (IEC) and the hydrophobicity of the alkene groups decrease both the equilibrium water uptake (%WU) and the number of water molecules per cation (λ) of the membrane. Interestingly, the evolution of λ with the IEC seems to follow that of undegraded PPO-TMA membranes. Additionally, we ascribe this to the similarity between the structure of the degraded monomer and the neutral one in the polymer. We conclude that most of the performance losses originate in the lower hydration of the degraded membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting Partial Atomic Charges in Metal–Organic Frameworks: An Extension to Ionic MOFs

Molecular simulation is an invaluable tool to predict and understand the usage of metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges, commonly obtained from density functional theory (DFT) calculations, are often required to model the electrostatic interactions between the MOF and adsorbates, especially when the adsorbates have dipole or quadrupole moments, such as water and CO 2 . Machine learning (ML) models have been previously employed to predict partial charges and avoid the computational cost associated with DFT calculations. However, previous ML models suffer from small training data sets, which limit their scope of application. In this work, we introduce two novel machine learning models, PACMOF2-neutral and PACMOF2-ionic, aimed at predicting the density-derived electrostatic and chemical (DDEC6) partial atomic charges for both neutral and ionic MOFs. These models not only yield DFT-level accuracy at a fraction of the computational cost but also demonstrate a remarkable improvement in prediction of adsorption, as validated with grand canonical Monte Carlo simulations. Furthermore, the robustness and fast computational time of the PACMOF2 models, along with their transferability to other porous materials such as covalent organic frameworks and zeolites, underscores their potential in high-throughput screening of MOFs for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hydrophobic Nanoconfinement Enhances CO 2 Conversion to H 2 CO 3

Understanding the formation of H 2 CO 3 in water from CO 2 is important in environmental and industrial processes. Although numerous investigations have studied this reaction, the conversion of CO 2 to H 2 CO 3 in nanopores, and how it differs from that in bulk water, has not been understood. We use ReaxFF metadynamics molecular simulations to demonstrate striking differences in the free energy of CO 2 conversion to H 2 CO 3 in bulk and nanoconfined aqueous environments. We find that nanoconfinement not only reduces the energy barrier but also reverses the reaction from endothermic in bulk water to exothermic in nanoconfined water. Also, charged intermediates are observed more often under nanoconfinement than in bulk water. Stronger solvation and more favorable proton transfer with increasing nanoconfinement enhance the thermodynamics and kinetics of the reaction. Here our results provide a detailed mechanistic understanding of an important step in the carbonation process, which depends intricately on confinement, surface chemistry, and CO 2 concentration.

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

Nuclear Quantum Effects on the Electronic Structure of Water and Ice

The electronic properties and optical response of ice and water are intricately shaped by their molecular structure, including the quantum mechanical nature of the hydrogen atoms. Despite numerous previous studies, a comprehensive understanding of the nuclear quantum effects (NQEs) on the electronic structure of water and ice at finite temperatures remains elusive. Here, we utilize molecular simulations that harness efficient machine-learning potentials and many-body perturbation theory to assess how NQEs impact the electronic bands of water and hexagonal ice. By comparing path-integral and classical simulations, we find that NQEs lead to a larger renormalization of the fundamental gap of ice, compared to that of water, ultimately yielding similar bandgaps in the two systems, consistent with experimental estimates. Our calculations suggest that the increased quantum mechanical delocalization of protons in ice, relative to water, is a key factor leading to the enhancement of NQEs on the electronic structure of ice.

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

Nanoscale hydration in layered manganese oxides

Birnessite is a layered MnO2 mineral with a nanoscale interlayer region that accommodates layers of intercalated water. Variable distributions of Mn sites with oxidation states of II, III and IV are responsible for the high catalytic reactivity of birnessite in nature and as a valuable candidate for energy storage solutions. We here report water loading capabilities and the vibrational spectral signatures of two forms of birnessite of strongly contrasting particle size. Using X-ray diffraction we find that potassium-birnessite accommodates no more than one monolayer (1W) of water in its interlayer region. Molecular simulations show that this is an energetically favorable hydration state where interlayer potassium-water and direct water-birnessite interactions are greater than at other hydration levels. Simulations also suggest a stable 2W state but that is not achieved experimentally by contact with water vapor. Finally, this work provides a means to predict the distribution of adsorbed and interlayer water molecules using a recently developed composite model.

Cheng, Wei↗
DOE OSTI2022

Exponential Water Uptake in Ionomer Membranes Results from Polymer Plasticization

Water content is the most influential factor in the performance of ion-exchange membranes (IEMs), controlling their mechanical rigidity, ionic conductivity, and degradation rates. Membranes absorb water exponentially at high water activity (a w ), making that region of the sorption isotherm the most influential on membrane properties. Plasticization of the polymer by water has been proposed to cause this exponential uptake at high a w . However, an integrated microscopic picture of the structure, thermodynamic, and mechanical properties of IEMs as a function of a w has remained elusive. In this work, we use large-scale molecular simulations validated with experimental measurements to compute the sorption isotherms, Young’s modulus, polymer dynamics, and structure of IEMs. The simulations unveil that the exponential increase in water uptake (WU) coincides in all cases with the glass to rubber transition of the membrane, as measured through its Young’s modulus and segmental polymer dynamics. Functionalization of the polymer with alkyl groups further contributes to its plasticization, increasing the WU at a given a w and ion-exchange capacity (IEC). The alkyl chains act synergistically with water to plasticize the polymer matrix and allow water penetration in the membrane. The simulations reveal that the width of the water channels depends on the ratio λ of water to ions in the polymer but is independent of its IEC. We conclude that differences in the polymer matrix–not the water channels–are responsible for the distinct uptake response of IEMs to the thermodynamic driving force of water activity.

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

Control of Structural Hydrophobicity and Cation Solvation on Interlayer Water Transport during Clay Dehydration

Swelling clay hydration/dehydration is important to many environmental and industrial processes. Experimental studies usually probe equilibrium hydration states in an averaged manner and thus cannot capture the fast water transport and structural change in interlayers during hydration/dehydration. Using molecular simulations and thermogravimetric analyses, we observe a two-stage dehydration process. The first stage is controlled by evaporation at the edges: water molecules near hydrophobic sites and the first few water molecules of the hydration shell of cations move fast to particle edges for evaporation. The second stage is controlled by slow desorption of the last 1–2 water molecules from the cations and slow transport through the interlayers. The two-stage dehydration is strongly coupled with interlayer collapse and the coordination number changes of cations, all of which depend on layer charge distribution. Furthermore, this mechanistic interpretation of clay dehydration can be key to the coupled chemomechanical behavior in natural/engineered barriers.

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

Hydrogen Adsorption in Ultramicroporous Metal–Organic Frameworks Featuring Silent Open Metal Sites

In this work, we utilized an ultramicroporous metal–organic framework (MOF) named [Ni 3 (pzdc) 2 (ade) 2 (H 2 O) 4 ]·2.18H 2 O (where H 3 pzdc represents pyrazole-3,5-dicarboxylic acid and ade represents adenine) for hydrogen (H 2 ) adsorption. Upon activation, [Ni 3 (pzdc) 2 (ade) 2 ] was obtained, and in situ carbon monoxide loading by transmission infrared spectroscopy revealed the generation of open Ni(II) sites. The MOF displayed a Brunauer–Emmett–Teller (BET) surface area of 160 m 2 /g and a pore size of 0.67 nm. Hydrogen adsorption measurements conducted on this MOF at 77 K showed a steep increase in uptake (up to 1.93 mmol/g at 0.04 bar) at low pressure, reaching a H 2 uptake saturation at 2.11 mmol/g at ~0.15 bar. The affinity of this MOF for H 2 was determined to be 9.7 ± 1.0 kJ/mol. In situ H 2 loading experiments supported by molecular simulations confirmed that H 2 does not bind to the open Ni(II) sites of [Ni 3 (pzdc) 2 (ade) 2 ], and the high affinity of the MOF for H 2 is attributed to the interplay of pore size, shape, and functionality.

08 HYDROGEN↗
DOE OSTI2020

Crossover from Rouse to Reptation Dynamics in Salt-Free Polyelectrolyte Complex Coacervates

Considerable interest in the dynamics and rheology of polyelectrolyte complex coacervates has been motivated by their industrial application as viscosity modifiers. A central question is the extent to which classical Rouse and reptation models can be applied to systems where electrostatic interactions play a critical role on the thermodynamics. By relying on molecular simulations, we present a direct analysis of the crossover from Rouse to reptation dynamics in salt-free complex coacervates as a function of chain length. This crossover shifts to shorter chain lengths as electrostatic interactions become stronger, which corresponds to the formation of denser coacervates. To distinguish the roles of Coulomb interactions and density, here we compare the dynamics of coacervates to those of neutral, semidilute solutions at the same density. Both systems exhibit a universal dynamical behavior in the connectivity-dominated (subdiffusion and normal diffusion) regimes, but the monomer relaxation time in coacervates is much longer and increases with increasing Bjerrum length. This is similar to the cage effect observed in glass-forming polymers, but the local dynamical slowdown is caused here by strong Coulomb attractions (ion pairing) between oppositely charged monomers. Our findings provide a microscopic framework for the quantitative understanding of coacervate dynamics and rheology.

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

Polymer Architecture-Induced Trade-off between Conductivities and Transference Numbers in Salt-Doped Polymeric Ionic Liquids

Recent experiments have demonstrated that polymeric ionic liquids that share the same cation and anion but possess different architectures can exhibit markedly different conductivity and transference number characteristics when doped with lithium salt. In this study, we used atomistic molecular simulations on polymer chemistries inspired by the experiments to probe the mechanistic origins underlying the competition between conductivity and transference numbers. Our results indicate that the architecture of the polycationic ionic liquid plays a subtle but crucial role in modulating the anion–cation interactions, especially their dynamical coordination characteristics. Chemistries leading to longer-lived anion–cation coordinations relative to lithium–anion coordinations lead to lower conductivities and higher transference numbers. Furthermore, our results suggest that higher conductivities are accompanied by lower transference numbers and vice versa, revealing that alternative approaches may need to be considered to break this trade-off in salt-doped polyILs.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Elastomer Mechanics of Cross-Linked Linear-Ring Polymer Blends

Cross-linking a blend of linear and ring polymers creates a new topology-based dual-network elastomer in which the two components differ significantly in their topology. We use molecular simulations and topological analysis to examine key mechanical properties as functions of ring polymer volume fraction Φ R . For Φ R < Φ R *, where the rings begin to overlap, the network shear modulus G and the maximum stretch ratio λ p are weakly dependent on Φ R . For Φ R > Φ R *, entanglements trapped in the network are diluted as the rings overlap, leading to a significant decrease in G and an increase in λ p with increasing Φ R . Here, the peak tensile stress, σ p , exhibits a maximum around Φ R *, indicating an enhancement of network strength due to the stronger cohesion from the entanglements between linear and ring polymers.

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

The Effects of Morphology and Hydration on Anion Transport in Self-Assembled Nanoporous Membranes

Ordered nanoporous polymer membranes offer opportunities for systematically probing the mechanisms of ion transport under confinement and for realizing useful materials for electrochemical devices. Here, we examine the impact of morphology and ion hydration on the transport of hydroxide and bromide anions in nanostructured polymer membranes with 1 nm scale pores. We use aqueous lyotropic self-assembly of an amphiphilic monomer, with a polymerizable surfactant to create direct hexagonal (H I ) and gyroid mesophases. UV-induced cross-linking leads to the formation of nanoporous polymers with water continuous channels. The membranes are mechanically robust and chemically durable, resisting degradation during extended exposure to 1 M NaOH solutions. We use a combination of electrochemical impedance spectroscopy, pulsed-field gradient NMR spectroscopy, and molecular simulations to elucidate anion and water transport. The as-prepared hexagonal systems display higher conductivity and lower activation energies for both anions relative to the gyroid system. When compared at equivalent hydration, however, gyroid and hexagonal membranes show similar activation energies, with nearly identical conductivities at ambient temperatures. Both ionic conductivity and water diffusivity increase with increasing hydration. Here, the water uptake as a function of relative humidity for the hexagonal and gyroid mesophases ultimately dictates the water diffusion and magnitude of the ionic conductivity, with the hexagonal system showing overall higher capacity for hydration and thus faster ion transport. The durability of these materials under aggressive alkaline conditions and their relatively high hydroxide ion conductivity suggest that these nanostructured polymers could be of interest as membranes for alkaline fuel cells.

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

Can Multiple Ions in an Ionic Liquid Improve the Biomass Pretreatment Efficacy?

we report over the last few decades, efforts to transition the global production of fuels and chemicals toward renewable carbon feedstocks have accelerated. A large portion of these efforts have focused on valorization of one of the most abundant renewable carbon sources, lignocellulose. Pretreatment of lignocellulose is the first critical step in this process. In this study, novel ionic liquid (IL) systems consisting of multiple ions known to be effective at biomass pretreatment were tested on woody and grassy biomass. Molecular simulations and experimental results established the synergistic advantages of combining specific individual components in these systems. For pine (woody) biomass, pretreatment with the combination of imidazolium, cholinium, acetate, and lysinate ions achieved 80% glucose and 70% xylose yields at high biomass loading. For sorghum biomass, an IL system comprising cholinium, lysinate, and palmitate ions not only enabled a 98% glucose yield but was also found to be biocompatible in a one-pot configuration, producing the biofuel precursor bisabolene using an engineered strain of the yeast Rhodosporidium toruloides.

09 BIOMASS FUELS↗
DOE OSTI2025

Molecular Insights into Geochemical Reactions of Iron-Bearing Minerals: Implications for Hydrogen Geo-Storage

Here, this study investigates the reaction of hydrogen (H 2 ) with pyrite (FeS 2 ), focusing on how temperature and the presence of water influence the reaction pathways and kinetics. Utilizing computational molecular simulations and kinetic analyses, we explore the impact of these factors on the formation of hydrogen sulfide (H 2 S) and related species. First, grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations reveal that physical H 2 adsorption occurs in distinct layers on the pyrite surface. In addition, increased temperatures reduce the absolute adsorption capacities. Reactive MD simulations demonstrate that H 2 interacts differently with pyrite under varying conditions. At 298 K, H 2 reacts with pyrite to form HS – , leading to the formation of HS – through covalent bonding with sulfur of pyrite. However, no H 2 S is produced at this temperature, suggesting that a kinetic barrier (i.e., activation energy) may prevent this reaction. At higher temperatures, H 2 S production significantly increases. The presence of water introduces additional complexity to the reaction mechanism. Unlike dry conditions, water enhances H 2 S generation, even at low temperatures. Water also facilitates the formation of additional products, such as SOH, indicating a more intricate chemical environment on the pyrite surface. Our findings identify the association of HS – ions to form H 2 S as the rate-limiting step, with temperature influencing this process. This finding suggests that while the presence of water can create a more dynamic reaction environment, the overall mechanisms leading to H 2 S formation remain consistent. These outcomes suggest the need for developing targeted strategies to manage and control H 2 S emissions within the context of underground hydrogen storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗