Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvation”

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 757 records · Page 42

A hybrid calorimetry-simulation model of mixing enthalpy for molten salt

Calorimetric determination of enthalpies of mixing (ΔH mix ) in multicomponent molten salts is often interpreted using empirical models that lack physically meaningful parameters. However, for improving pyrochemical separation of spent nuclear fuel, where lanthanides are major fission products and critical elements, a deeper thermodynamic understanding of the link between excess thermodynamic properties and solvation structure is critically needed. In this work, we implement a hybrid and physics-informed framework, MIVM+Calorimetry+AIMD, which integrates experimentally measured ΔH mix (via high temperature drop calorimetry) with solvation structures from ab initio molecular dynamics (AIMD). This approach is demonstrated using LaCl 3 mixed with eutectic LiCl-KCl (58 mol% – 42 mol%) at 873 K and 1133 K. MIVM-derived parameters enable extrapolation of excess Gibbs energy and La 3+ activity across compositions. In contrast, direct ΔH mix predictions from AIMD and polarizable ion model simulations deviate significantly. By incorporating experimentally benchmarked solvation structures into an interpretable thermodynamic model, the MIVM+Calorimetry+AIMD formalism achieves higher accuracy and generalizable method for studying molten salts, offering a robust path for understanding and optimizing molten salt chemistry relevant to nuclear fuel cycles and separation science.

Goncharov, Vitaliy G. [Washington State Univ., Pul↗

A fundamental study on the [(μ-Cl) 3 Mg 2 (THF) 6 ] + dimer electrolytes for rechargeable Mg batteries

We present a fundamental study on [(μ-Cl) 3 Mg 2 (THF) 6 ] + dimer electrolytes using various physical methods including Subambient Pressure Ionization with Nanoelectrospray Mass spectrometry (SPIN-MS), Raman spectroscopy, 25 Mg{ 1 H} NMR, 27 Al{ 1 H} NMR and electrochemical analysis. For the first time, long time sought THF solvated [MgCl] + species was experimentally characterized by SPIN mass spectrometry in the solution of the Mgdimer containing electrolyte, confirming the mono-Cl- abstraction reaction between MgCl 2 and an Al Lewis acid. Solvated MgCl 2 in the electrolyte was confirmed by Raman spectroscopy. The experimental results establish the previously proposed dimerization equilibrium of solvated [MgCl] + and MgCl2 with [(μ-Cl) 3 Mg 2 (THF)6] + . 25 Mg{ 1 H} NMR, 27 Al{ 1 H} NMR and electrochemical analysis on chloration reaction of [(μ-Cl) 3 Mg 2 (THF) 6 ]AlPh 3 Cl with external Cl- led to further insights on the coordination chemistry of the dimer electrolyte. Finally, a comprehensive mechanism is proposed for the reversible electrochemical Mg deposition and stripping and Mg 2+ and Cl- ion transports of the Mg dimer electrolytes in rechargeable Mg batteries.

25 ENERGY STORAGE↗

A density functional theory based approach for predicting melting points of ionic liquids

Accurate prediction of melting points of ILs is important both from the fundamental point of view and from the practical perspective for screening ILs with low melting points and broadening their utilization in a wider temperature range. In this work, we present an ab initio approach to calculating melting points of ILs with known crystal structures and illustrate its application for a series of 11 ILs containing imidazolium/pyrrolidinium cations and halide/polyatomic fluoro-containing anions. The melting point is determined as a temperature at which the Gibbs free energy of fusion is zero. The Gibbs free energy of fusion can be expressed through the use of the Born-Fajans-Haber cycle via the lattice free energy of forming a solid IL from gaseous phase ions and the sum of the solvation free energies of ions comprising IL. Dispersion-corrected density functional theory (DFT) involving (semi)local (PBE-D3) and hybrid exchange-correlation (HSE06-D3) functionals is applied to estimate the lattice enthalpy, entropy, and free energy. The ions solvation free energies are calculated with the SMD-generic-IL solvation model at the M06-2X/6-31+G(d) level of theory under standard conditions. The melting points of ILs computed with the HSE06-D3 functional are in good agreement with the experimental data, with a mean absolute error of 30.5 K and a mean relative error of 8.5%. The model is capable of accurately reproducing the trends in melting points upon variation of alkyl substituents in organic cations and replacement one anion by another. The results verify that the lattice energies of ILs containing polyatomic fluoro-containing anions can be approximated reasonably well using the volume-based thermodynamic approach. However, there is no correlation of the computed lattice energies with molecular volume for ILs containing halide anions. Moreover, entropies of solid ILs follow two different linear relationships with molecular volume for halides and polyatomic fluoro-containing anions. As a result, continuous progress in predicting crystal structures of organic salts with halide anions will be a key factor for successful prediction of melting points with no prior knowledge of the crystal structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular mechanism and binding free energy of doxorubicin intercalation in DNA

The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of Δ G ele , Δ G vDW , Δ G polar , Δ G non-polar , T Δ S total , Δ G deform , Δ G con , and Δ G ion . To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (Δ G deform ), the energy cost from translational and rotational entropic contributions ( T Δ S tran+rot ), the total electrostatic interactions (Δ G polar-PB/GB + Δ G ele ) of incorporation, the intramolecular electrostatic interactions (Δ G ele ) and electrostatic polar solvation interactions (Δ G polar-PB/GB ) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (Δ G vDW ), the non-polar solvation interaction (Δ G non-polar ), the vibrational entropic contribution ( T Δ S vib ), and the standard concentration dependent free energies of DOX (Δ G con ) and the ionic solution (Δ G ion ). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX–DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. Here, we have calculated the total binding free energy (BFE) (Δ G t-sim ) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX–DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX–DNA1 and DOX–DNA2 in the 1 : 1 complex is –9.1 and –5.1 kcal mol –1 respectively. The same quantities from the modified MM-PBSA method are –12.74 and –8.35 kcal mol –1 respectively. The value of the total BFE Δ G t-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of –7.7 ± 0.3 kcal mol –1 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular dynamics simulations of alkaline earth metal ions binding to DNA reveal ion size and hydration effects

The identity of metal ions surrounding DNA is key to its biological function and materials applications. In this work, we compare atomistic molecular dynamics simulations of double strand DNA (dsDNA) with four alkaline earth metal ions (Mg 2+ , Ca 2+ , Sr 2+ , and Ba 2+ ) to elucidate the physical interactions that govern DNA–ion binding. Simulations accurately model the ion–phosphate distance of Mg 2+ and reproduce ion counting experiments for Ca 2+ , Sr 2+ , and Ba 2+ . Our analysis shows that alkaline earth metal ions prefer to bind at the phosphate backbone compared to the major groove and negligible binding occurs in the minor groove. Larger alkaline earth metal ions with variable first solvation shells (Ca 2+ , Sr 2+ , and Ba 2+ ) show both direct and indirect binding, where indirect binding increases with ion size. Mg 2+ does not fit this trend because the strength of its first solvation shell predicts indirect binding only. Ions bound to the phosphate backbone form fewer contacts per ion compared to the major groove. Within the major groove, metal ions preferentially bind to guanine–cystosine base pairs and form simultaneous contacts with the N7 and O6 atoms of guanine. Overall, we find that the interplay among ion size, DNA–ion interaction, and the size and flexibility of the first solvation shell are key to predicting how alkaline earth metal ions interact with DNA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis, modular composition, and electrochemical properties of lamellar iron sulfides

Transition metal chalcogenides with layered structures are promising materials for energy storage, catalysis, and electronics, among other areas. We have identified a new layered phase of iron sulfide containing interlayer solvated cations. We present an optimized synthesis for the Li + -containing material from an Fe(III) xanthate complex. Structure and composition data indicate the material consists of poorly-ordered iron sulfide layers separated by solvated cations. The lamellar spacing in these materials can be tuned by changing the identity of the cation. Furthermore, the lamellar spacing can also be reversibly tuned by the degree of solvation of the material. The material is electrically conductive and can serve as a pseudocapacitor with comparable performance to commercial materials such as MnO 2 . Finally, these materials also show promise as lithium or sodium ion battery cathodes with good capacity and reversibility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of diluent concentration in localized high concentration electrolytes: elucidation of hidden diluent-Li + interactions and Li + transport mechanism

Localized high concentration electrolytes (LHCE) offer a viable dilution strategy for high concentration electrolytes (HCE) as the dilution process barely impacts the enhanced reductive/oxidative behavior of the HCE formulation but significantly lowers the overall viscosity and, in most cases, increases the ionic conductivity. On the other hand, experimental studies indicate that fluorinated ether electrolytes such as 1,1,2,2-tetrafluoroethylene 2,2,3,3-tetrafluoropropyl ether (TTE) help grow enhanced passivation layers on Ni-rich NMC cathodes. In this work, we study LHCE formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethyl carbonate (DMC), and TTE as the diluent. Here we use molecular dynamics methodologies, and Raman spectra measurements, to address to what extent the diluent content impacts the coordination behavior of the aggregated structures in the electrolyte and to evaluate the Li + transport properties under the influence of an external electric field. In contrast to other fluorinated ethers, we find that TTE interacts with Li + via fluorine atoms, partially limiting the DMC–Li + interactions hence altering the Li + solvation coordination. This competitive interaction with Li+ between the organic solvent and the TTE diluent influences the electrolyte's reductive/oxidative behavior. Nevertheless, the bonding strength of the Li + –F TTE is much weaker than those of the Li + –O DMC and Li + –OFSI – . Therefore, the existence of Li + –F TTE is in a transient state rather than in a steady state. These results provide plausible guiding rules for future dilution strategies of HCE electrolytes. We also demonstrate that Li + ions drift under the electric field's influence via repeated ion dissociation/association processes following a hopping conduction mechanism. Li + ions jump between aggregated networks where Li–O interactions dominate via diluent-enriched phases, a process in which the solvation shells temporarily mutate to a Li–F dominated coordination structure. We expect our results to contribute an improved atomic-level understanding of the solvation structure and dynamics of LHCE electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) of zwitterionic organic HCl salts

In this work, we benchmark quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) for determining the crystal structures of two zwitterionic organic HCl salts, L-ornithine HCl ( Orn ) and L-histidine HCl·H 2 O ( Hist ). These salts present an interesting challenge for QNMRX-CSP, as gas-phase geometry optimizations used to generate starting structures for the organic zwitterionic fragments fail to capture their correct solid-state geometries. To overcome this limitation, geometry optimizations using the COSMO water-solvation model are employed to generate initial structural models. Using this approach, QNMRX-CSP yields structural models of the two zwitterionic organic HCl salts that closely match experimentally determined crystal structures. In addition, the application of QNMRX-CSP to Hist represents a further step toward the de novo structural determination of solvated organic HCl salts, as Hist is the first benchmark system of this type to include a water molecule as a component of its crystal structure. This work is significant for its potential application to the structural determination of active pharmaceutical ingredients, which often feature complex organic components and solvated solid forms.

Fleischer, Carl H. [Florida State Univ., Tallahass↗

Molecular density functional theory for water with liquid-gas coexistence and correct pressure

The solvation of hydrophobic solutes in water is special because liquid and gas are almost at coexistence. In the common hypernetted chain approximation to integral equations, or equivalently in the homogenous reference fluid of molecular density functional theory, coexistence is not taken into account. Hydration structures and energies of nanometer-scale hydrophobic solutes are thus incorrect. In this article, we propose a bridge functional that corrects this thermodynamic inconsistency by introducing a metastable gas phase for the homogeneous solvent. We show how this can be done by a third order expansion of the functional around the bulk liquid density that imposes the right pressure and the correct second order derivatives. Although this theory is not limited to water, we apply it to study hydrophobic solvation in water at room temperature and pressure and compare the results to all-atom simulations. The solvation free energy of small molecular solutes like n-alkanes and hard sphere solutes whose radii range from angstroms to nanometers is now in quantitative agreement with reference all atom simulations. The macroscopic liquid-gas surface tension predicted by the theory is comparable to experiments. This theory gives an alternative to the empirical hard sphere bridge correction used so far by several authors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

C–IBI: Targeting cumulative coordination within an iterative protocol to derive coarse-grained models of (multi-component) complex fluids

We present a coarse-graining strategy that we test for aqueous mixtures. The method uses pair-wise cumulative coordination as a target function within an iterative Boltzmann inversion (IBI) like protocol. We name this method coordination iterative Boltzmann inversion (C–IBI). While the underlying coarse-grained model is still structure based and, thus, preserves pair-wise solution structure, our method also reproduces solvation thermodynamics of binary and/or ternary mixtures. In addition, we observe much faster convergence within C–IBI compared to IBI. To validate the robustness, we apply C–IBI to study test cases of solvation thermodynamics of aqueous urea and a triglycine solvation in aqueous urea.

74 ATOMIC AND MOLECULAR PHYSICS↗

Stochastic level-set variational implicit-solvent approach to solute-solvent interfacial fluctuations

Recent years have seen the initial success of a variational implicit-solvent model (VISM), implemented with a robust level-set method, in capturing efficiently different hydration states and providing quantitatively good estimation of solvation free energies of biomolecules. The level-set minimization of the VISM solvation free-energy functional of all possible solute-solvent interfaces or dielectric boundaries predicts an equilibrium biomolecular conformation that is often close to an initial guess. In this work, we develop a theory in the form of Langevin geometrical flow to incorporate solute-solvent interfacial fluctuations into the VISM. Such fluctuations are crucial to biomolecular conformational changes and binding process. We also develop a stochastic level-set method to numerically implement such a theory. We describe the interfacial fluctuation through the “normal velocity” that is the solute-solvent interfacial force, derive the corresponding stochastic level-set equation in the sense of Stratonovich so that the surface representation is independent of the choice of implicit function, and develop numerical techniques for solving such an equation and processing the numerical data. We apply our computational method to study the dewetting transition in the system of two hydrophobic plates and a hydrophobic cavity of a synthetic host molecule cucurbit[7]uril. Numerical simulations demonstrate that our approach can describe an underlying system jumping out of a local minimum of the free-energy functional and can capture dewetting transitions of hydrophobic systems. In the case of two hydrophobic plates, we find that the wavelength of interfacial fluctuations has a strong influence to the dewetting transition. In addition, we find that the estimated energy barrier of the dewetting transition scales quadratically with the inter-plate distance, agreeing well with existing studies of molecular dynamics simulations. Our work is a first step toward the inclusion of fluctuations into the VISM and understanding the impact of interfacial fluctuations on biomolecular solvation with an implicit-solvent approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The influence of mesoscopic confinement on the dynamics of imidazolium-based room temperature ionic liquids in polyether sulfone membranes

The structural dynamics of a series of 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (C n mimNTf 2 , n = 2, 4, 6, 10: ethyl—Emim; butyl—Bmim; hexyl—Hmim; decyl—Dmim) room temperature ionic liquids confined in the pores of polyether sulfone (PES 200) membranes with an average pore size of ~350 nm and in the bulk liquids were studied. Time correlated single photon counting measurements of the fluorescence of the fluorophore coumarin 153 (C153) were used to observe the time-dependent Stokes shift (solvation dynamics). The solvation dynamics of C153 in the ionic liquids are multiexponential decays. The multiexponential functional form of the decays was confirmed as the slowest decay component of each bulk liquid matches the slowest component of the liquid dynamics measured by optical heterodyne-detected optical Kerr effect (OHD-OKE) experiments, which is single exponential. The fact that the slowest component of the Stokes shift matches the OHD-OKE data in all four liquids identifies this component of the solvation dynamics as arising from the complete structural randomization of the liquids. Although the pores in the PES membranes are large, confinement on the mesoscopic length scale results in substantial slowing of the dynamics, a factor of ~4, for EmimNTf2, with the effect decreasing as the chain length increases. By DmimNTf 2 , the dynamics are virtually indistinguishable from those in the bulk liquid. The rotation relaxation of C153 in the four bulk liquids was also measured and showed strong coupling between the C153 probe and its environment.

14 SOLAR ENERGY↗

Unraveling the Spectral Signatures of Solvent Ordering in K-edge XANES of Aqueous Na +

The aqueous solvation structure of the Na + ion is studied both by molecular dynamics (MD) simulations and K-edge X-ray absorption near edge spectroscopy (XANES). In this paper, we present a systematic study contrasting the differences in the predicted XANES spectra, using molecular configurations generated from classical and quantum MD. Using a set of suitable order parameters, we elucidate how the spectroscopic features are influenced by the local solvation structure around the Na + ion. Lastly, our findings suggest that XANES is sensitive to fluctuations between the first and second solvation shells and can distinguish between classical and quantum interaction representations used to generate molecular ensembles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water

Many physical phenomena must be accounted for to accurately model solution-phase optical spectra line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein (GFP) and photoactive yellow protein (PYP). By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we introduce how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We discover that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Radiolytic redox interplay defines nanomaterial synthesis in liquids

Irradiation of a liquid solution generates solvated electrons and radiolysis products, which can lead to material deposition or etching. The chemical environment dictates the dominant reactions. Radiolysis-induced reactions in salt solutions have substantially different results in pure water versus water-ammonia, which extends the lifetime of solvated electrons. We investigate the interplay between transport and solution chemistry via the example of solid silver formation from e-beam irradiation of silver nitrate solutions in water and water-ammonia. The addition of ammonia results in the formation of a secondary ring-shaped deposit tens of microns in diameter (formed over tens of seconds) around the primary point of deposition (formed over milliseconds). Simulations uncover the relative importance of oxidizing and reducing reactions and transport effects. Our explanation of this behavior involves mechanisms beyond ammonia’s role in extending solvated electron lifetimes.

36 MATERIALS SCIENCE↗

Impact of the LiPF6 Concentration on the Interfacial Charge Transfer and Fast-charging Capabilities of Lithium-Ion Batteries

Fast-charging lithium-ion batteries (LIB) demand optimized electrolyte formulations to balance ionic conductivity, viscosity, and interfacial charge transfer kinetics. This study examines how LiPF 6 concentration shapes solvation structure, desolvation energy, charge transfer activation energy, and solid electrolyte interphase (SEI) properties, which are critical for fast-charging performance. Using Raman spectroscopy, electrochemical cycling, X-ray photoelectron spectroscopy, and atomistic modeling, we analyze how varying LiPF 6 concentrations impact interfacial and bulk transport properties. Our findings show that increasing LiPF 6 concentration alters lithium solvation structures, reduces desolvation energy, and accelerates charge transfer at the electrode interface. Higher concentrations lower the activation energy for charge transfer and suppress excessive SEI growth, improving interfacial kinetics. However, concentrations above a certain threshold increase viscosity and reduce ionic conductivity, limiting transport efficiency. These results offer insights into electrolyte solvation and interfacial charge transfer mechanisms, providing guidelines for designing next-generation fast-charging LIB electrolytes with enhanced efficiency, stability, and longevity.

Son, Seoung-Bum [Argonne National Laboratory (ANL)↗

Transitioning NWChem to the Next Generation of Manycore Machines

The NorthWest chemistry (NWChem) modeling software is a popular molecular chemistry simulation software that was designed from the start to work on massively parallel processing supercomputers [1-3]. It contains an umbrella of modules that today includes self-consistent eld (SCF), second order Møller-Plesset perturbation theory (MP2), coupled cluster (CC), multiconguration self-consistent eld (MCSCF), selected conguration interaction (CI), tensor contraction engine (TCE) many body methods, density functional theory (DFT), time-dependent density functional theory (TDDFT), real-time time-dependent density functional theory, pseudopotential plane-wave density functional theory (PSPW), band structure (BAND), ab initio molecular dynamics (AIMD), Car-Parrinello molecular dynamics (MD), classical MD, hybrid quantum mechanics molecular mechanics (QM/MM), hybrid ab initio molecular dynamics molecular mechanics (AIMD/MM), gauge independent atomic orbital nuclear magnetic resonance (GIAO NMR), conductor like screening solvation model (COSMO), conductor-like screening solvation model based on density (COSMO-SMD), and reference interaction site model (RISM) solvation models, free energy simulations, reaction path optimization, parallel in time, among other capabilities [4]. Moreover, new capabilities continue to be added with each new release.

Bylaska, Eric J.↗

Molecular Basis for Microbial Adhesion and Geochemical Surface Reactions: A Study Across Scales

Computational chemistry was used to help provide a molecular level description of the interactions of Gram-negative microbial membranes with subsurface materials. The goal is to develop a better understanding of the molecular processes involved in microbial metal binding, microbial attachment to mineral surfaces, and, eventually, oxidation/reduction reactions (electron transfer) that can occur at these surfaces and are mediated by the bacterial exterior surface. The project focused on the interaction of the outer microbial membrane, which is dominated by an exterior lipopolysaccharide (LPS) portion, of Pseudomonas aeruginosa with the mineral goethite and with solvated ions in the environment. This was originally a collaborative project with T.P. Straatsma and B. Lowery of the Pacific Northwest National Laboratory. The University of Alabama effort used electronic structure calculations to predict the molecular behavior of ions in solution and the behavior of the sugars which form a critical part of the LPS. The interactions of the sugars with metal ions are expected to dominate much of the microscopic structure and transport phenomena in the LPS. This work, in combination with the molecular dynamics simulations of Straatsma and the experimental electrochemistry and microscopy measurements of Lowry, both at PNNL, is providing new insights into the detailed molecular behavior of these membranes in geochemical environments. The effort at The University of Alabama has three components: solvation energies and structures of ions in solution, prediction of the acidity of the critical groups in the sugars in the LPS, and binding of metal ions to the sugar anions. An important aspect of the structure of the LPS membrane as well as ion transport in the LPS is the ability of the sugar side groups such as the carboxylic acids and the phosphates to bind positively charged ions. We are studying the acidity of the acidic side groups in order to better understand the ability of these groups to bind metal ions. We need to understand the solvation properties of the metal ions in solution and their ability to bind not only to the sugars but to proteins and to other anions. Our goal is then to be able to predict the ability of the side groups to bind metal ions. One result from the earlier molecular dynamics simulations is the exclusion of water from the inner hydrophobic part of the membrane. We thus need to investigate the binding of the cations in media with different dielectric constants.

58 GEOSCIENCES↗