Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydration free energy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)-an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Finally, analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration-physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

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

Consequences of the failure of equipartition for the p – V behavior of liquid water and the hydration free energy components of a small protein

Earlier we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step δ t ≤ 0.5 fs to ensure equipartition between translational and rotational modes. Here we extend that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examine integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small δ t is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, and the dielectric constant, a measure of dipole moment fluctuations, are also seen to be sensitive to δ t . Using the mean volume estimated from the NpT simulation, we examine the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to δ t . In going from δ t = 2 fs to δ t = 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein.

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

Dataset for manuscript "Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein"

Previously, we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step dt that is less than or equal to 0.5 fs to ensure equipartition between translational and rotational modes. We extended that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examined integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small time-step, dt, is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, is seen to be sensitive to dt. Using the mean volume estimated from the NpT simulation, we examined the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to dt. In going from a time-step of 2 fs to a time-step of 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein. The data-set contains the simulation metadata and log files that support the claims noted above.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

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

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

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

Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F – over Cl – in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory–experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

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

Dialing in Direct Air Capture of CO 2 by Crystal Engineering of Bis-iminoguanidines

Direct air capture (DAC) technologies that extract carbon dioxide from the atmosphere via chemical processes have the potential to restore the atmospheric CO 2 concentration to an optimal level. This study elucidates structure-property relationships in DAC by crystallization of bis-iminoguanidine (BIG) carbonate salts. Here, their crystal structures are analyzed by X-ray and neutron diffraction to accurately measure key structural parameters including molecular conformations, hydrogen bonding, and π-stacking. Experimental measurements of key properties, such as aqueous solubilities and regeneration energies and temperatures, are complemented by first-principles calculations of lattice and hydration free energies, as well as free energies of reactions with CO 2 , and BIG regenerations. Minor structural modifications in the molecular structure of the BIGs are found to result in major changes in the crystal structures and the aqueous solubilities within the series, leading to enhanced DAC.

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

Interplay of physically different properties leading to challenges in separating lanthanide cations – an ab initio molecular dynamics and experimental study

Lanthanide elements have well-documented similarities in their chemical behavior, which make the valuable trivalent lanthanide cations (Ln 3+ ) particularly difficult to separate from each other in water. Here, we apply ab initio molecular dynamics simulations to compare the free energies (Δ G ads ) associated with the adsorption of lanthanide cations to silica surfaces at a pH condition where SiO – groups are present. The predicted Δ G ads for lutetium (Lu 3+ ) and europium (Eu 3+ ) are similar within statistical uncertainties; this is in qualitative agreement with our batch adsorption measurements on silica. This finding is remarkable because the two cations exhibit hydration free energies (Δ G hyd ) that differ by >2 eV, different hydration numbers, and different hydrolysis behavior far from silica surfaces. We observe that the similarity in Lu 3+ and Eu 3+ Δ G ads is the result of a delicate cancellation between the difference in Eu 3+ and Lu 3+ hydration (Δ G hyd ), and their difference in binding energies to silica. We propose that disrupting this cancellation at the two end points, either for adsorbed or completely desorbed lanthanides ( e.g. , via nanoconfinment or mixed solvents), will lead to effective Ln 3+ separation.

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

The H•/H– Redox Couple and Absolute Hydration Energy of H–

A supermolecule-continuum approach with water clusters up to n = 16 H2O molecules has been used to predict the absolute hydration free energies at 298 K ($\Delta$G hyd ) of both hydrogen (H•) and hydride (H - ) to be 4.6 and -78 ± 2 kcal/mol, respectively. These values are combined with a high accuracy prediction of the gas phase electron affinity ($\Delta$G gas,298K = -16.9 kcal/mol) to determine the aqueous electron affinity of H• of 99.5 ± 2 kcal/mol, which yields a reduction potential for H• vs SHE of -0.03 ± 0.15 V. This value is in agreement within 0.2 V with most estimates obtained using a wide variety of approaches. The results of this study can be used to improve the absolute hydricity scale in water which provides additional insights into how a putative hydride interacts with solvent, but do not change the ability to predict the relative reactivity of two species using relative hydricity scales.

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

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗
DOE OSTI2023

Ion solvation as a predictor of lanthanide adsorption structures and energetics in alumina nanopores

Adsorption reactions at solid-water interfaces define elemental fate and transport and enable contaminant clean-up, water purification, and chemical separations. For nanoparticles and nanopores, nanoconfinement may lead to unexpected and hard-to-predict products and energetics of adsorption, compared to analogous unconfined surfaces. Here we use X-ray absorption fine structure spectroscopy and operando flow microcalorimetry to determine nanoconfinement effects on the energetics and local coordination environment of trivalent lanthanides adsorbed on Al 2 O 3 surfaces. We show that the nanoconfinement effects on adsorption become more pronounced as the hydration free energy, ΔG hydr , of a lanthanide decreases. Neodymium (Nd 3+ ) has the least exothermic ΔG hydr (-3336 kJ·mol -1 ) and forms mostly outer-sphere complexes on unconfined Al 2 O 3 surfaces but shifts to inner-sphere complexes within the 4 nm Al 2 O 3 pores. Lutetium (Lu 3+ ) has the most exothermic ΔG hydr (-3589 kJ·mol -1 ) and forms inner-sphere adsorption complexes regardless of whether Al 2 O 3 surfaces are nanoconfined. Importantly, the energetics of adsorption is exothermic in nanopores only, and becomes endothermic with increasing surface coverage. Changes to the energetics and products of adsorption in nanopores are ion-specific, even within chemically similar trivalent lanthanide series, and can be predicted by considering the hydration energies of adsorbing ions.

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

Impact of Ligands on the Ion–Ion Selectivity of Ligand-Appended-Pillar[ n ]arene Channels

Ligand-appended pillar[5]arene (LAP) channels are emerging as a promising platform for ion-selective membranes. In this study, we investigated the ion selectivities of three LAP channels functionalized with distinct ligands: diglycolamine, carbamoylmethyl phosphine oxide, and propionamide phosphonic acid. These ligands were chosen based on their demonstrated affinities for lanthanides in solvent extraction studies. We examined the selectivity of each channel toward a series of monovalent, divalent, and trivalent ions: Li + , Na + , K + , Mg 2+ , Ca 2+ , La 3+ , Eu 3+ , and Yb 3+ . To quantify ion-channel interactions and the energetics of translocation, we computed the potential of mean force (PMF) profiles for each ion. These PMFs were subsequently recast into permeability coefficients, enabling a direct evaluation of ion permselectivity. Our results reveal distinct selectivity patterns governed by the chemical nature of the appended ligands. The diglycolamine-functionalized channel exhibits strong interactions with monovalent ions, leading to replication of the Li + bulk hydration shell within the channel. This promotes preferential transport Li + ions while effectively excluding divalent and trivalent species. In contrast, the carbamoylmethyl phosphine oxide-functionalized channel displays a reduced selectivity between monovalent ions while similarly exhibiting an effective rejection of divalent and trivalent ions. This arises from steric hindrance around channel-lining oxygens imposed by the ligand's bulky phenyl groups which limits heavy ion coordination and promotes retention of hydration shells effectively increasing effective ion size. The propionamide phosphonic acid-functionalized channel exhibits a different trend, favoring ions with lower hydration free energies, with K + ions showing enhanced permeation. These findings highlight the critical role of ligand-ion interactions in modulating ion selectivity within LAP channels. In particular, the strong coordination exhibited by diglycolamine and propionamide phosphonic acid ligands significantly influences the ion transport energetics and selectivity profiles.

hydration↗
DOE OSTI2023

Cation–polymer interactions and local heterogeneity determine the relative order of alkali cation diffusion coefficients in PEGDA hydrogels

Current research efforts are focused on endowing polymer membranes with ion–ion selectivity by incorporating ion–polymer interactions into materials to bias the selective partitioning and or diffusivity of one species over another. However, little is known about the impact of such interactions on the mechanisms of ion transport. In this study, we probe the influence of cation–polymer interactions on cation, anion, and salt diffusivity in a model membrane material, poly(ethylene glycol) diacrylate (PEGDA) by modeling concentrated polyethylene oxide solutions via molecular dynamics simulations. These results are compared to published experimental data for LiCl, NaCl, and KCl diffusion in PEGDA. Experimentally, the order of salt and cation diffusion coefficients for LiCl, NaCl, and KCl deviate from the order in aqueous solutions. Here, simulations identify these deviations to arise from cation–polymer coordination in the membrane. Both the fraction of bound cations and the average binding lifetime increases with decreasing cation hydration free energy (moving down the alkali series), leading to different diffusivity trends in the membrane compared to solution. However, to recover the experimentally observed order of diffusivities cations and salt in our simulations, we needed to incorporate membrane heterogeneity explicitly via a polymer charge scaling procedure. Together, our results indicate that cation–polymer interactions, as well as spatial heterogeneity within the membrane, play a critical role in dictating the observed order of alkali cation and salt diffusion coefficients in membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Lanthanide transport in angstrom-scale MoS 2 -based two-dimensional channels

Rare earth elements (REEs), critical to modern industry, are difficult to separate and purify, given their similar physicochemical properties originating from the lanthanide contraction. Here, we systematically study the transport of lanthanide ions (Ln 3+ ) in artificially confined angstrom-scale two-dimensional channels using MoS 2 -based building blocks in an aqueous environment. The results show that the uptake and permeability of Ln 3+ assume a well-defined volcano shape peaked at Sm 3+ . This transport behavior is rooted from the tradeoff between the barrier for dehydration and the strength of interactions of lanthanide ions in the confinement channels, reminiscent of the Sabatier principle. Molecular dynamics simulations reveal that Sm 3+ , with moderate hydration free energy and intermediate affinity for channel interaction, exhibit the smallest dehydration degree, consequently resulting in the highest permeability. Our work not only highlights the distinct mass transport properties under extreme confinement but also demonstrates the potential of dialing confinement dimension and chemistry for greener REEs separation.

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

Hydration Mimicry by Membrane Ion Channels

Ions transiting biomembranes might pass readily from water through ion-specific membrane proteins if these protein channels provide environments similar to the aqueous solution hydration environment. In fact, bulk aqueous solution is an important reference condition for the ion permeation process. Assessment of this hydration mimicry concept depends on understanding the hydration structure and free energies of metal ions in water in order to provide a comparison for the membrane channel environment. To refine these considerations, we review local hydration structures of ions in bulk water and the molecular quasi-chemical theory that provides hydration free energies. In doing so, we note some current views of ion binding to membrane channels and suggest new physical chemical calculations and experiments that might further clarify the hydration mimicry concept.

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

Beyond Thermodynamics: Assessing the Dynamical Softness of Hydrated Ions from First Principles

Ion (de)hydration is a key rate-determining step in interfacial processes from corrosion to electrochemical energy storage. However, predicting the kinetics of ion (de)hydration remains challenging, prompting the use of static proxies such as hydration energy and valence. While useful for assessing thermodynamic preferences, such descriptors cannot fully capture the dynamical softness of the hydration shell that dictates kinetics. Accordingly, we use first-principles molecular dynamics to analyze hydration shell softness for a diverse set of metal cations. Additionally, three dynamic metrics are introduced to intuitively describe the bond rigidity, shape deformability, and exchange fluidity of the solvation shell. Together, these metrics capture the relevant physics in the static descriptors, while offering a far more complete and efficient representation for the overall propensity for (de)hydration. Application to the hydrated ion set demonstrates a weak connection between dynamical softness and hydration energy, confirming that dynamical descriptors of hydration are key for correctly describing ion transfer processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structure–Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca)

Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities-and their impacts on environmental processes-requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca–Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ~0.6 nm between alginate strands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗