SEARCH · Search NASA

Results for “hydration shell”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2024

DRIFT: Diffusivity Regulation of Isotopes by nanoconFinementT

Spontaneous isotope fractionation has been reported under nanoconfinement conditions in naturally occurring systems, but the origin of this phenomena is currently unknown. Two existing hypotheses have been proposed, one based on changes in the solvation environment of the isotopes that reduces the non-mass dependent hydrodynamics contribution to diffusion. The other is that isotopes have mass-dependent surface adsorption, varying their total diffusion through nanoconfined channels. To investigate these hypotheses, benchtop experiments, nuclear magnetic resonance (NMR) spectroscopy, and molecule scale modeling were applied. Classical molecular dynamics simulations identified that the Na + and Cl - hydration shells across the three different salt solutions ( 22 Na 35 Cl, 23 Na 35 Cl, 24 Na 35 Cl) did not vary as a function of the Na + isotope, but that there was a significant pore size effect, with larger hydration shells at larger pore sizes. Additionally, while total adsorption times did not vary as a function of the Na + isotope or pore size, the free ion concentration, or those adsorbed on the surface for <5% of the simulation time did exhibit isotope dependence. Experimentally, challenges occurred developing a repeatable experiment, but NMR characterization of water diffusion rates through ordered alumina membranes was able to identify the existence of two distinct water environments associated with water inside and outside the pore. Further NMR studies could be used to confirm variation in hydration shells and diffusion rates of dissolved ions in water. Ultimately, mass-dependence adsorption is a primary driver of variations in isotope diffusion rates, rather than variation in hydration shells that occur under nanoconfinement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗

DOE OSTI2023

DEHP− extractant binding to trivalent lanthanide Er3+: Fast binding accompanied by concerted angular motions of hydration water

Solvent extraction of trivalent rare earth metal ions by organophosphorus extractants proceeds via binding of phosphoric acid headgroups to the metal ion. Water molecules in the tightly bound first hydration shell of the metal ions must be displaced by oxygen atoms from phosphoric acid headgroups. Here, we use classical molecular dynamics simulations to explore the event in which a fully hydrated Er3+ binds to its first phosphoric acid headgroup. Approach of the headgroup into the region between the first and second hydration shells leads to a fast ejection of a water molecule that is accompanied by reordering of the hydration water molecules, including discretization of their angular positions and collective rotation about the metal ion. The water molecule ejected from the first shell is located diametrically opposite from the binding oxygen. Headgroup binding places a headgroup oxygen closer to Er3+ than its first hydration shell and creates a loosely bound water that subsequently exchanges between the first shell and its environment. This second exchange of water also occurs at discrete angular positions. This geometrical aspect of binding may be of relevance to understanding the binding and transport of ion–extractant complexes that are expected to occur at the organic–aqueous liquid–liquid interface used in solvent extraction processes.

Chemistry↗

DOE OSTI2024

The effect of low temperature on poly(3-methyl- N -vinylcaprolactam)- b -poly( N -vinylpyrrolidone) diblock copolymer nanovesicles assembled from all-aqueous media

Nanosized polymeric vesicles (polymersomes) self-assembled from double hydrophilic copolymers of poly(3-methyl-N-vinylcaprolactam) n -b-poly(N-vinylpyrrolidone) m (PMVC n -b-PVPON m ) using all aqueous media are a promising platform for biomedical applications, because of their superior stability over liposomes in vivo and high loading capacity. Herein, we explored the temperature-sensitive behavior of PMVC 58 -b-PVPON 65 vesicles using transmission electron microscopy (TEM), dynamic light scattering (DLS), atomic force microscopy (AFM), and small-angle neutron scattering (SANS) in response to lowering the solution temperature from 37 to 25, 20, 14 and 4 °C. The copolymer vesicles with an average size of 350 nm at 37 °C were assembled from the diblock copolymer dissolved in aqueous solution at 4 °C. We show that while the polymersome's size gradually decreases upon the temperature decrease from 37 to 4 °C, the average shell thickness increases from 17 nm to 25 nm, respectively. SANS study revealed that the PMVC 58 -b-PVPON 65 vesicle undergoes a gradual structure evolution from a dense-shell vesicle at 37–25 °C to a highly-hydrated shell vesicle at 20–14 °C to molecular chain aggregates at 4 °C. From SANS contrast matching study, this vesicle behavior is found to be driven by the gradual rehydration of PMVC block at 37–14 °C. The shell hydration at 20–14 °C also correlated with the 4.4-fold decrease in the relative fluorescence intensity from vesicle-encapsulated fluorescent dye, indicating ~80% of the dye release within 12 hours after the vesicle exposure to 14 °C. No significant (<5%) dye release was observed for the vesicle solutions at 37–20 °C, indicating excellent cargo retention inside the vesicles. Our study provides new fundamental insights on temperature-sensitive polymer vesicles and demonstrates that the copolymer assembly into polymersomes can be achieved by decreasing a copolymer aqueous solution temperature below 14 °C followed by solution exposure to ≥20 °C. This type of all-aqueous assembly, instead of nanoprecipitation from organic solvents or solvent exchange, can be highly desirable for encapsulating a wide range of biological molecules, including proteins, peptides, and nucleic acids, into stable polymer vesicles without a need for organic solvents for dissolution of the copolymers that are amphiphilic at physiologically relevant temperatures of 20–37 °C.

36 MATERIALS SCIENCE↗

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 OSTI2016

Tuning the tetrahedrality of the hydrogen-bonded network of water: Comparison of the effects of pressure and added salts

Experiments and simulations demonstrate some intriguing equivalences in the effect of pressure and electrolytes on the hydrogen-bonded network of water. Here, we examine the extent and nature of equivalence effects between pressure and salt concentration using relationships between structure, entropy, and transport properties based on two key ideas: first, the approximation of the excess entropy of the fluid by the contribution due to the atom-atom pair correlation functions and second, Rosenfeld-type excess entropy scaling relations for transport properties. We perform molecular dynamics simulations of LiCl–H{sub 2}O and bulk SPC/E water spanning the concentration range 0.025–0.300 molefraction of LiCl at 1 atm and pressure range from 0 to 7 GPa, respectively. The temperature range considered was from 225 to 350 K for both the systems. To establish that the time-temperature-transformation behaviour of electrolyte solutions and water is equivalent, we use the additional observation based on our simulations that the pair entropy behaves as a near-linear function of pressure in bulk water and of composition in LiCl–H{sub 2}O. This allows for the alignment of pair entropy isotherms and allows for a simple mapping of pressure onto composition. Rosenfeld-scaling implies that pair entropy is semiquantitatively related to the transport properties. At a given temperature, equivalent state points in bulk H{sub 2}O and LiCl–H{sub 2}O (at 1 atm) are defined as those for which the pair entropy, diffusivity, and viscosity are nearly identical. The microscopic basis for this equivalence lies in the ability of both pressure and ions to convert the liquid phase into a pair-dominated fluid, as demonstrated by the O–O–O angular distribution within the first coordination shell of a water molecule. There are, however, sharp differences in local order and mechanisms for the breakdown of tetrahedral order by pressure and electrolytes. Increasing pressure increases orientational disorder within the first neighbour shell while addition of ions shifts local orientational order from tetrahedral to close-packed as water molecules get incorporated in ionic hydration shells. The variations in local order within the first hydration shell may underlie ion-specific effects, such as the Hofmeister series.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Hydrophilic and Apolar Hydration in Densely Grafted Cationic Brushes and Counterions with Large Mobilities

We employ all-atom molecular dynamics (MD) simulation framework to unravel water microstructure and ion properties for cationic [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] (PMETAC) brushes with chloride ions as counterions. First, we identify locally separate water domains (or first hydration shells) each around the {N(CH 3 ) 3 } + and the C=O functional groups of the PMETAC chain and one around the Cl - ion. These first hydration shells around the respective moieties overlap and the extent of the overlap depends on the nature of the species triggering it. Second, despite the overlap, the water molecules in these domains demonstrate disparate properties dictated by the properties of the atoms/groups around which they are located. For example, the presence of the methyl groups make the {N(CH 3 ) 3 } + group trigger apolar hydration as evidenced by the corresponding orientation of the dipole of the water molecules around the {N(CH 3 ) 3 } + moiety. These water molecules around the N(CH 3 ) 3 } + group also have enhanced tetrahedrality as compared to the water molecules constituting the hydration layer around the C=O group and the Cl - counterion. Our simulations also identify that there is an intervening water layer between the Cl - ion and {N(CH 3 ) 3 } + group: this layer prevents the Cl - ion from coming very close to the {N(CH 3 ) 3 } + group. As a consequence, there is a significantly large mobility of the Cl - ions inside the PMETAC brush layer. Furthermore, the C=O group of the PE chain, due to the partial negative charge on the oxygen atom and the specific structure of the PMETAC brush system, demonstrates strongly hydrophilic behavior and enforces a specific dipole response of water molecules analogous to that experienced by water around anionic species of high charge density. As a result, our findings confirm that PMETAC brushes undergo hydrophilic hydration at one site, apolar hydration at another site, and ensures large mobility of the supported Cl- counterions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1992

Damage to cellular DNA from particulate radiations, the efficacy of its processing and the radiosensitivity of mammalian cells. Emphasis on DNA double strand breaks and chromatin breaks

For several years, it has been evident that cellular radiation biology is in a necessary period of consolidation and transition (Lett 1987, 1990; Lett et al. 1986, 1987). Both changes are moving apace, and have been stimulated by studies with heavy charged particles. From the standpoint of radiation chemistry, there is now a consensus of opinion that the DNA hydration shell must be distinguished from bulk water in the cell nucleus and treated as an integral part of DNA (chromatin) (Lett 1987). Concomitantly, sentiment is strengthening for the abandonment of the classical notions of "direct" and "indirect" action (Fielden and O'Neill 1991; O'Neill 1991; O'Neill et al. 1991; Schulte-Frohlinde and Bothe 1991 and references therein). A layer of water molecules outside, or in the outer edge of, the DNA (chromatin) hydration shell influences cellular radiosensitivity in ways not fully understood. Charge and energy transfer processes facilitated by, or involving, DNA hydration must be considered in rigorous theories of radiation action on cells. The induction and processing of double stand breaks (DSBs) in DNA (chromatin) seem to be the predominant determinants of the radiotoxicity of normally radioresistant mammalian cells, the survival curves of which reflect the patterns of damage induced and the damage present after processing ceases, and can be modelled in formal terms by the use of reaction (enzyme) kinetics. Incongruities such as sublethal damage are neither scientifically sound nor relevant to cellular radiation biology (Calkins 1991; Lett 1990; Lett et al. 1987a). Increases in linear energy transfer (LET infinity) up to 100-200 keV micron-1 cause increases in the extents of neighboring chemical and physical damage in DNA denoted by the general term DSB. Those changes are accompanied by decreasing abilities of cells normally radioresistant to sparsely ionizing radiations to process DSBs in DNA and chromatin and to recover from radiation exposure, so they make significant contributions to the relative biological effectiveness (RBE) of a given radiation.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

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 OSTI2022

Vibrational signatures of HNO 3 acidity when complexed with microhydrated alkali metal ions, M + ∙(HNO 3 )(H 2 O) n =5 (M = Li, K, Na, Rb, Cs), at 20 K

The speciation of strong acids like HNO 3 under conditions of restricted hydration is an important factor in the rates of chemical reactions at the air-water interface. Here we explore the trade-offs in play when HNO 3 is attached to alkali ions (Li + -Cs + ) with four water molecules in their primary hydration shells. This is achieved by analyzing the vibrational spectra of the M + ∙(HNO 3 )(H 2 O) 5 clusters cooled to about 20 K in a cryogenic photofragmentation mass spectrometer. The local acidity of the acidic OH group is estimated by the extent of the red shift in its stretching frequency when attached to a single water molecule. Here, the persistence of this local structural motif (HNO 3 -H 2 O) in all of these alkali metal clusters enables us to determine the competition between the effect of the direct complexation of the acid with the cation, which acts to enhance acidity, and the role of the water network in the first hydration shell around the ions, which acts to counter (screen) the intrinsic effect of the ion. Analysis of the vibrational features associated with the acid molecule as well as those of the water network reveal how cooperative interactions in the micro hydration regime conspire to effectively offset the intrinsic enhancement of HNO 3 acidity afforded by attachment to the smaller cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Exploring the impact of ions on oxygen K-edge X-ray absorption spectroscopy in NaCl solution using the GW-Bethe-Salpeter-equation approach

X-ray absorption spectroscopy (XAS) is a powerful experimental tool to probe the local structure in materials with the core hole excitations. Here, the oxygen K-edge XAS spectra of the NaCl solution and pure water are computed by using a recently developed GW-Bethe-Salpeter equation approach, based on configurations modeled by path-integral molecular dynamics with the deep-learning technique. The neural network is trained on ab initio data obtained with strongly constrained and appropriately normed density functional theory. The observed changes in the XAS features of the NaCl solution, compared to those of pure water, are in good agreement between experimental and theoretical results. We provided detailed explanations for these spectral changes that occur when NaCl is solvated in pure water. Specifically, the presence of solvating ion pairs leads to localization of electron-hole excitons. As a result, our theoretical XAS results support the theory that the effects of the solvating ions on the H-bond network are mainly confined within the first hydration shell of ions, however beyond the shell the arrangement of water molecules remains to be comparable to that observed in pure water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Description of Hydration Water in Protein (Green Fluorescent Protein) Solution

The structurally and dynamically perturbed hydration shells that surround proteins and biomolecules have a substantial influence upon their function and stability. This makes the extent and degree of water perturbation of practical interest for general biological study and industrial formulation. Here, we present an experimental description of the dynamical perturbation of hydration water around green fluorescent protein in solution. Less than two shells (~5.5 Å) were perturbed, with dynamics a factor of 2–10 times slower than bulk water, depending on their distance from the protein surface and the probe length of the measurement. Furthermore, this dependence on probe length demonstrates that hydration water undergoes subdiffusive motions (τ ∝ q –2.5 for the first hydration shell, τ ∝ q –2.3 for perturbed water in the second shell), an important difference with neat water, which demonstrates diffusive behavior (τ ∝ q –2 ). Our results help clarify the seemingly conflicting range of values reported for hydration water retardation as a logical consequence of the different length scales probed by the analytical techniques used.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Mechanisms of interfacial reactivity in near surface and extreme environments

The local water structure surrounding ions in aqueous solutions greatly affects their chemical properties such as reaction rates, ion association, and proton and electron transport. These properties result in the behavior of ions in natural aqueous environments. For example ore transport is facilitated by chloride ion pair formation and the reaction of ions in an interface is strongly dependent on the dehydration of the ion hydration shell. We are developing the use of high-resolution XAFS observations and 1st principles based MD-XAFS analysis (spectra simulated using 1st principle methods with no adjustable parameters, AIMD) to interpret the solution properties of strongly interacting aqueous solutes under arbitrary pressure and temperature conditions. In the 1st principle MD-XAFS method density functional theory (DFT) based MD simulations(Car and Parrinello, 1985) are used to generate a large ensemble of structural snap shots of the hydration region. These are then used to generate scattering intensities. I emphasize three points about this novel approach to analyzing XAFS data. 1st: As illustrated in Figure 1, the level of agreement between the calculated and observed intensities is considerably higher than has been obtained by any XAFS analysis to date (note 2nd shell region, R> 2 Å). 2nd: This result was obtained from a parameter free simulation with no fitting of the interaction potentials to any data. This supports the use of these methods for more difficult environments and more complex solutes (polyions). 3rd: New information about the shell structure (Figure 1) is now available because of this more detailed agreement. We note also that both multiple scattering and second shell features are well represented in the analysis. As far as we know this is the 1st analysis of second shell structure and multiple scattering. Excellent agreement has been obtained for most of the third row metal ions: Ca 2+ , Zn 2+ , Cu 2+ , Ni 2+ , Co 2+ , Mn 2+ , Fe 3+ , Cr 3+ . Calculations on these systems are demanding because of their open electronic shells, and high ionic charge. Principal Investigator: Professor John Weare (University of California, San Diego) The prediction of the interactions of geochemical fluids with minerals, nanoparticles, and colloids under extreme near surface conditions of temperature (T) and pressure (P) is a grand challenge research need in geosciences (U.S. DOE 2007, Basic Research Needs for Geosciences: Facilitating the 21st Energy Systems.). To evaluate the impact of these processes on energy production and management strategies it is necessary to have a high level of understanding of the interaction between complex natural fluids and mineral formations. This program emphasizes 1st principle parameter free simulations of complex chemical processes in solutions, in the mineral phase, and in the interfaces between these phases The development of new computational tools (with emphasis on oxide materials and reaction dynamics) tailored to treat wide range of conditions and time scales experienced in such geochemical applications is has been developed. Because of the sensitivity of the interaction in these systems to electronic structure and local bonding environments, and of the need to describe bond breaking/formation, our simulations are based on interactions calculated at the electronic structure level (ab-initio molecular dynamics, AIMD). The progress in the computational aspects of program may be summarized in terms of the following themes (objectives); Development of efficient parameter free dynamical simulation technology based on 1st principles force and energy calculations especially adapted for geochemical applications (e.g., mineral, interfaces and aqueous solutions) (continuing program); Calculation of the dynamics of water structure of in the surface-water interface of transition metal oxides and oxihydroxides; and Development of improved (beyond DFT+GGA) electronic structure calculations for minerals and the interface region that more accurately calculate electron correlation, spin density, and localization. The focus of the program is also on the iron oxide and oxihydroxide minerals and Fe 2+ (aq)/Fe 3+ (cr) oxidation in the mineral solution interface region. These methods included the development of model Hamilitonian methods that can be solved to near convergence for single site models (DMFT) and many-body perturbation methods (MP2, GW); Development of time decomposition methods to extend time scales of molecular dynamics (MD) simulations and support the use of high complexity electronic structure calculations (MP2, CCSD(T)) of forces for use in dynamical simulations where very high chemical accuracy is required (microsolvated reactions in absorbed surface layers); and The development of a new linear scaling finite element solver for eigenvalue problem that supports solution of quantum problems with unusual potential and boundary values. Application progress of the above new simulation technology to problems of geochemical interests includes; The prediction of metal oxide surface structure and the reduction/oxidation of Fe 3+ (cr)/Fe 2+ (aq) in metal oxide (hematite, goethite)/solution interfaces. Result: water interacts strongly with the 001 Hematite surface; interaction of water with the 100 goethite is weak; The study of ion solvation and the composition of ion hydration shells under extreme conditions (focus on Fe 3+/2+ , Al 3+ and Mg 2+ and their hydroxide speciation). Result: Ion association in water solutions can be calculated from 1st principle methods. Efficient sampling of the free energy requires more development; The continued development of new high resolution analysis of XAFS scattering of disordered systems (particularly Al, Mg) and of XANES calculations for aqueous ions. Result: EXAFS spectra can be calculated to high accuracy with DFT level dynamic simulations; The exploration of electron localization and electron transport in metal oxides (highly correlated materials). Result: Proper description of electron localization requires levels of calculation beyond DFT; and Localization of electrons in DFT type Hamiltonians was studied. Result: For very Dirac high exchange new solutions (New unphysical bifrucations) to the eigenvalue problem are found. The program was highly collaborative involving faculty and students in mathematics, physics and computer science departments as well as coworkers at the Pacific Northwest National Laboratories (PNNL). The students in this program had the opportunity to develop skills in the development of methods, the implementation of method on high performance parallel computers and the application of these methods to problem in geochemical science. Much of the software that was developed was incorporated in the NWchem software package maintained by PNNL.

58 GEOSCIENCES↗

DOE OSTI2020

Microsecond Molecular Dynamics Simulations of Proteins Using a Quasi-Equilibrium Solvation Shell Model

In this work, we describe the development and implementation of a quasi-equilibrium hydration shell model of biomolecular solvation with adaptive boundaries. Applying the model to microsecond-long molecular dynamics simulations of several protein systems of varying complexity, we find that the model simulation results are of comparable quality to those obtained from simulations of fully solvated systems, but at a reduced computational cost. We discuss the dominant sources of error in the model and outline directions for future improvements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Reversible (De)Intercalation of Hydrated Zn 2+ in Mg 2+ -Stabilized V 2 O 5 Nanobelts with High Areal Capacity

The rechargeable aqueous zinc ion battery (ZIB) is regarded as one of the most promising candidates for large-scale energy storage applications due to its low-cost and eco-friendly properties. However, the development of a suitable cathode operating with high areal capacity and uncovering the relevant reaction mechanisms remain challenging. Herein, the application of Mg 0.26 V 2 O 5 ∙0.73H 2 O (MVO) nanobelts as a ZIB cathode is demonstrated. In situ FT-IR reveals the shift of OH stretching from 3350 cm -1 to 3200 cm -1 , corresponding to the hydration shell of Zn 2+ , while in situ Raman suggests the interlayer charges creening effect, which would boost the intercalation of hydrated Zn 2+ . Density function theory reveals that the hydrated Zn 2+ can lower the Coulombic repulsion at the electrode-electrolyte interface and circumvents the desolvation penalty of hydrated Zn 2+ during the (de)intercalation process. Additionally, excellent structure stability and large interlayer spacing guarantee the highly reversible (de)intercalation of hydrated Zn 2+ . Therefore, the MVO nanobelts exhibit a high areal capacity of 2.12 mAh cm -2 at 0.05 A g -1 , outstanding cycling stability of 2500 cycles at 10 A g -1 with a mass loading of 5 mg cm -2 . Finally, it is believed that the use of hydrated intercalation charge carriers will boost further studies in other multivalent rechargeable batteries.

25 ENERGY STORAGE↗

DOE OSTI2015

Structure of the ordered hydration of amino acids in proteins: analysis of crystal structures

The hydration of protein crystal structures was studied at the level of individual amino acids. The dependence of the number of water molecules and their preferred spatial localization on various parameters, such as solvent accessibility, secondary structure and side-chain conformation, was determined. Crystallography provides unique information about the arrangement of water molecules near protein surfaces. Using a nonredundant set of 2818 protein crystal structures with a resolution of better than 1.8 Å, the extent and structure of the hydration shell of all 20 standard amino-acid residues were analyzed as function of the residue conformation, secondary structure and solvent accessibility. The results show how hydration depends on the amino-acid conformation and the environment in which it occurs. After conformational clustering of individual residues, the density distribution of water molecules was compiled and the preferred hydration sites were determined as maxima in the pseudo-electron-density representation of water distributions. Many hydration sites interact with both main-chain and side-chain amino-acid atoms, and several occurrences of hydration sites with less canonical contacts, such as carbon–donor hydrogen bonds, OH–π interactions and off-plane interactions with aromatic heteroatoms, are also reported. Information about the location and relative importance of the empirically determined preferred hydration sites in proteins has applications in improving the current methods of hydration-site prediction in molecular replacement, ab initio protein structure prediction and the set-up of molecular-dynamics simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2018

Molecular Dynamics Simulation of Hydration and Swelling of Mixed-Layer Clays

Swelling of clay minerals is important to a broad class of problems in science and engineering. While the problem has been extensively studied experimentally, past molecular modeling of the phenomenon was focused on pure clays of one type or another. In practice, however, there is a diverse class of mixed-layer clays (MLCs) in sedimentary rock with intermixed stacking sequence of two or more types of distinct layers within a single crystal. In fact, more than 60% of sedimentary rocks in the U.S. contain various types of MLCs. In this work, we present the results, to our knowledge, of the first molecular dynamics simulation of hydration energetics and swelling of illite–montmorillonite (I–MMT) MLCs, the most common type of mixed clays. The swelling is studied as a function of the water concentration with three combinations of interlayer cations, namely, Na + and K + . The hydration energies, the radial distribution functions, and the density profiles in the interlayer region are computed. For regular Na–MMT with layer charge concentrated in the octahedral sheet, weak cation–surface interaction results in fully hydrated ions and significant swelling. In the asymmetric interlayer of the MLC, however, the illite sheet with stronger interaction of surface and ions causes adsorption of the cations deep in the ditrigonal cavities of the siloxane surface. Given that the hydration enthalpy of K + is smaller than that of Na + , its hydration shell is quite unstable compared with that of Na+. Therefore, swelling is inhibited as the ratio K + /Na + increases. Finally, the results demonstrate the significant differences between the hydration and swelling properties of pure clays and the mixed ones, which have important implications in practice, particularly for sequestration of CO 2 in sedimentary rock.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Unconventional Complex Coacervation between Neutral Polymer and Inorganic Polyoxometalate in Aqueous Solution via Direct Water Mediation

Water, specifically in a hydration shell, is critical for many biological supramolecular aggregations in nature, where water can directly mediate intermolecular association via hydrogen bonding and is regarded as “structured water”. Conversely, little has been reported on the biomimetic water-mediated supramolecular assembly with adequately high water content to date, because of the competing thermodynamic processes of water hydration and water as a building block to participate in self-assembly. In this work, we explore water-mediated complexation based on entropy-driven biphasic coacervate formation using highly hydrophilic neutral polymer and inorganic mineral-analogous nanoclusters. For the first time (to the best of our knowledge), nonelectrostatic liquid–liquid separating coacervate formation is demonstrated between polyethylene glycol (PEG) and polyoxometalate (POM) nanoclusters in aqueous solutions of varied PEG and POM concentrations, POM types, and aqueous medium conditions. Comprehensive characterization using fluorescence microscopy, small-angle X-ray scattering, calorimetry, and other techniques has confirmed that the compositions, microstructure, and thermodynamics of PEG–POM complex coacervation are highly similar to entropy-driven complex coacervation between oppositely charged polyelectrolytes in aqueous solution. However, the effect of heavy water on critical POM concentration for the onset of coacervate formation suggests that water, instead of the counter ions as commonly debated for polyelectrolyte complex coacervation, is responsible for PEG–POM coacervate formation. Specifically, structured water works as a hydrogen bond donor for both highly hydrated PEG and POM to directly mediate the PEG–water–POM association, resulting in the release of excess hydrated water for entropy-driven PEG–POM complex coacervation. Furthermore, water-mediated complex coacervation could be developed as a general and simple strategy to build biomimetic hybrid nanomaterials with high water content for various applications from energy-related functional nanomaterials to biomedical ramification.

36 MATERIALS SCIENCE↗