1
|
Wei Z, Xu H, Xu X, Feng G, Zheng W, Li T. Solvation of magnesium chloride dimer in water: The case of anionic and neutral clusters. J Chem Phys 2023; 158:2888211. [PMID: 37140000 DOI: 10.1063/5.0146319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2023] [Accepted: 04/24/2023] [Indexed: 05/05/2023] Open
Abstract
The structures of magnesium chloride dimer-water clusters, (MgCl2)2(H2O)n-/0, were investigated with size-selected anion photoelectron spectroscopy and theoretical calculations to understand the dissolution of magnesium chloride in water. The most stable structures were confirmed by comparing vertical detachment energies (VDEs) with the experimental measurements. A dramatic drop of VDE at n = 3 has been observed in the experiment, which is in accordance with the structural change of (MgCl2)2(H2O)n-. Compared to the neutral clusters, the excess electron induces two significant phenomena in (MgCl2)2(H2O)n-. First, the planar D2h geometry can be converted into a C3v structure at n = 0, making the Mg-Cl bonds easier to be broken by water molecules. More importantly, a negative charge-transfer-to-solvent process occurs after adding three water molecules (i.e., at n = 3), which leads to an obvious deviation in the evolution of the clusters. Such electron transfer behavior was noticed at n = 1 in monomer MgCl2(H2O)n-, indicating that the dimerization between two MgCl2 molecules can make the cluster more capable of binding electron. In neutral (MgCl2)2(H2O)n, this dimerization provides more sites for the added water molecules, which can stabilize the entire cluster and maintain its initial structure. Specifically, filling the coordination number to be 6 for Mg atoms can be seen as a link between structural preferences in the dissolution of the monomers, dimers, and extended bulk-state of MgCl2. This work represents an important step forward into fully understanding the solvation of MgCl2 crystals and other multivalent salt oligomers.
Collapse
Affiliation(s)
- Zhiyou Wei
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325001, China
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hongguang Xu
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xiling Xu
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Gang Feng
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 401331, China
| | - Weijun Zheng
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Tao Li
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325001, China
- Oujiang Laboratory (Zhejiang Lab for Regenerative Medicine, Vision and Brain Health), Wenzhou, Zhejiang 325001, China
| |
Collapse
|
2
|
Fan K, Zhang Y, Qiu Y, Zhang H. Impacts of targeting different hydration free energy references on the development of ion potentials. Phys Chem Chem Phys 2022; 24:16244-16262. [PMID: 35758314 DOI: 10.1039/d2cp01237e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H+). The choice of ΔG°(H+) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K+, Mg2+, Al3+, and Cl- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H+) for the cation and an increase in ΔG°(H+) for the anion. The increase in ΔG°(H+) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H+) value of -259.8 kcal mol-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H+) for Mg2+ model is better for a reasonable description of the metal binding network.
Collapse
Affiliation(s)
- Kun Fan
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| |
Collapse
|
3
|
Venetsanos F, Anogiannakis SD, Theodorou DN. Mixing Thermodynamics and Flory–Huggins Interaction Parameter of Polyethylene Oxide/Polyethylene Oligomeric Blends from Kirkwood–Buff Theory and Molecular Simulations. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c00642] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Fotis Venetsanos
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| | - Stefanos D. Anogiannakis
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| | - Doros N. Theodorou
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| |
Collapse
|
4
|
Wang X, Clegg SL, Di Tommaso D. Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions. J Chem Phys 2022; 156:024502. [PMID: 35032987 DOI: 10.1063/5.0074970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Chemical thermodynamic models of solvent and solute activities predict the equilibrium behavior of aqueous solutions. However, these models are semi-empirical. They represent micro-scale ion and solvent behaviors controlling the macroscopic properties using small numbers of parameters whose values are obtained by fitting to activities and other partial derivatives of the Gibbs energy measured for the bulk solutions. We have conducted atomistic simulations of aqueous electrolyte solutions (MgCl2 and CaCl2) to determine the parameters of thermodynamic hydration models. We have implemented a cooperative hydration model to categorize the water molecules in electrolyte solutions into different subpopulations. The value of the electrolyte-specific parameter, k, was determined from the ion-affected subpopulation with the lowest absolute value of the free energy of removing the water molecule. The other equilibrium constant parameter, K1, associated with the first degree of hydration, was computed from the free energy of hydration of hydrated clusters. The hydration number, h, was determined from a reorientation dynamic analysis of the water subpopulations compared to bulk-like behavior. The reparameterized models [R. H. Stokes and R. H. Robinson, J. Solution Chem. 2, 173 (1973) and Balomenos et al., Fluid Phase Equilib. 243, 29 (2006)] using the computed values of the parameters lead to the osmotic coefficients of MgCl2 solutions that are consistent with measurements. Such an approach removes the dependence on the availability of experimental data and could lead to aqueous thermodynamic models capable of estimating the values of solute and solvent activities as well as thermal and volumetric properties for a wide range of compositions and concentrations.
Collapse
Affiliation(s)
- Xiangwen Wang
- Department of Chemistry, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | - Simon L Clegg
- School of Environmental Sciences, University of East Anglia, Norwich NR4 7TJ, United Kingdom
| | - Devis Di Tommaso
- Department of Chemistry, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| |
Collapse
|
5
|
Rubinson KA. Measurement of Enthalpies and Entropies of Activation as a Function of Pairwise Distance for the Pairwise Relative Diffusion of SrI 2 in Water over Lengthscales from 6 Å to 40 Å. NEW J CHEM 2022. [DOI: 10.1039/d2nj00712f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Quasi-elastic neutron scattering (QENS) has been used to provide insight into the motions within solutions for many decades, and coherent QENS (CQENS) determines the motions of one ion relative to...
Collapse
|
6
|
Wei ZY, Yang LJ, Gong SY, Xu HG, Xu XL, Gao YQ, Zheng WJ. Comparison of the Microsolvation of CaX 2 (X = F, Cl, Br, I) in Water: Size-Selected Anion Photoelectron Spectroscopy and Theoretical Calculations. J Phys Chem A 2021; 125:3288-3306. [PMID: 33872010 DOI: 10.1021/acs.jpca.1c00573] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To understand the microsolvation of alkaline-earth dihalides in water and provide information about the dependence of solvation processes on different halides, we investigated CaBr2(H2O)n-, CaI2(H2O)n-, and CaF2(H2O)n- (n = 0-6) clusters using size-selected anion photoelectron spectroscopy and conducted theoretical calculations on these clusters and their neutrals. The results are compared with those of CaCl2(H2O)n-/0 clusters reported previously. It is found that the vertical detachment energies (VDEs) of CaCl2(H2O)n-, CaBr2(H2O)n-, and CaI2(H2O)n- show a similar trend with increasing cluster size, while the VDEs of CaF2(H2O)n- show a different trend. The VDEs of CaF2(H2O)n- are much lower than those of CaCl2(H2O)n-, CaBr2(H2O)n-, and CaI2(H2O)n-. A detailed probing of the structures shows that a significant increase of the Ca-X distance (separation of Ca2+-X- ion pair) in CaCl2(H2O)n-/0, CaBr2(H2O)n-/0, and CaI2(H2O)n-/0 clusters occurred at about n = 5. However, for CaF2(H2O)n-/0, no abrupt change of the Ca-F distance with the increasing cluster size has been observed. In CaCl2(H2O)6-/0, CaBr2(H2O)6-/0, and CaI2(H2O)6-/0, the Ca atom coordinates directly with 5 H2O molecules. However, in CaF2(H2O)n-/0, the Ca atom coordinates directly with only 2 or 3 H2O molecules. The similarity or differences in the structures and coordination numbers are consistent with the fact that CaCl2, CaBr2, and CaI2 have similar solubility, while CaF2 has much lower solubility.
Collapse
Affiliation(s)
- Zhi-You Wei
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Li-Jiang Yang
- Beijing National Laboratory for Molecular Science, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Shi-Yan Gong
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hong-Guang Xu
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xi-Ling Xu
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yi Qin Gao
- Beijing National Laboratory for Molecular Science, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.,Biomedical Pioneering Innovation Center, Peking University, Beijing 100871, China
| | - Wei-Jun Zheng
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.,University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
7
|
Ploetz EA, Karunaweera S, Bentenitis N, Chen F, Dai S, Gee MB, Jiao Y, Kang M, Kariyawasam NL, Naleem N, Weerasinghe S, Smith PE. Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20. J Chem Theory Comput 2021; 17:2964-2990. [PMID: 33878263 DOI: 10.1021/acs.jctc.1c00075] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.
Collapse
Affiliation(s)
- Elizabeth A Ploetz
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Sadish Karunaweera
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nikolaos Bentenitis
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Feng Chen
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Shu Dai
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Moon B Gee
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Yuanfang Jiao
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Myungshim Kang
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nilusha L Kariyawasam
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nawavi Naleem
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | | | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| |
Collapse
|
8
|
Characterization of a Bacillus megaterium strain with metal bioremediation potential and in silico discovery of novel cadmium binding motifs in the regulator, CadC. Appl Microbiol Biotechnol 2021; 105:2573-2586. [PMID: 33651131 DOI: 10.1007/s00253-021-11193-2] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2021] [Revised: 02/11/2021] [Accepted: 02/17/2021] [Indexed: 12/13/2022]
Abstract
Bioremediation of toxic metal ions using bacterial strains is a promising tool. Metal binding motifs in microbial proteins are involved in the regulation and transport of such toxic metals for metal detoxification. A bacterial strain designated TWSL_4 with metal (Cu, Cd, and Pb) resistance and removal ability was isolated and identified as a Bacillus megaterium strain using 16S rRNA gene analysis. An operon with 2 open reading frames (ORFs) was identified, cloned, and sequenced. ORF1 and ORF2 were identical to the cadmium efflux system accessory protein (CadC) and cadmium-translocating P-type ATPases (CadA) of B. megaterium strain YC4-R4 respectively. A protein homology search using Swiss model retrieved no crystal structures for CadC and CadA of Bacillus sp.. CadC of TWSL_4 had a sequence identity of 53% to the CadC (121aa) protein and 51.69% to the CadC crystal structure (1U2W.1.B; GMQE=0.75) of Staphylococcus sp. pI258. Molecular dynamic simulation studies revealed the presence of three metal binding regions in CadC of TWSL_4, [ASP7-TYR9], [ASP100-HIS102], and [LYS113-ASP116]. This is the first report showing evidence for the presence of Cd2+ and Zn2+ metal binding motifs in the CadC regulator of the Bacillus megaterium cad operon. The bacterial strain TWSL_4 was also found to contain two different P type ATPases encoding genes, cadA and zosA involved in metal resistance. Furthermore, the metal bioremediation potential of strain TWSL_4 was confirmed using an industrial effluent. KEY POINTS: • Isolation of a metal-resistant bacterial strain with potential for industrial bioremediation. • Discovery of novel Cd binding sites in CadC of the cad operon from B. megaterium. • Involvement of aspartic acid in the coordination of metal ions (Cd2+).
Collapse
|
9
|
Wei ZY, Yang LJ, Xu HG, Farooq U, Xu XL, Gao YQ, Zheng WJ. Hydration processes of barium chloride: Size-selected anion photoelectron spectroscopy and theoretical calculations of BaCl 2-water clusters. J Chem Phys 2020; 153:134301. [PMID: 33032412 DOI: 10.1063/5.0021991] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In order to understand the hydration processes of BaCl2, we investigated BaCl2(H2O)n - (n = 0-5) clusters using size-selected anion photoelectron spectroscopy and theoretical calculations. The structures of neutral BaCl2(H2O)n clusters up to n = 8 were also investigated by theoretical calculations. It is found that in BaCl2(H2O)n -/0, the Ba-Cl distances increase very slowly with the cluster size. The hydration process is not able to induce the breaking of a Ba-Cl bond in the cluster size range (n = 0-8) studied in this work. In small BaCl2(H2O)n clusters with n ≤ 5, the Ba atom has a coordination number of n + 2; however, in BaCl2(H2O)6-8 clusters, the Ba atom coordinates with two Cl atoms and (n - 1) water molecules, and it has a coordination number of n + 1. Unlike the previously studied MgCl2(H2O)n - and CaCl2(H2O)n -, negative charge-transfer-to-solvent behavior has not been observed for BaCl2(H2O)n -, and the excess electron of BaCl2(H2O)n - is mainly localized on the Ba atom rather on the water molecules. No observation of Ba2+-Cl- separation in current work is consistent with the lower solubility of BaCl2 compared to MgCl2 and CaCl2. Considering the BaCl2/H2O mole ratio in the saturated solution, one would expect that about 20-30 H2O molecules are needed to break the first Ba-Cl bond in BaCl2.
Collapse
Affiliation(s)
- Zhi-You Wei
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Li-Jiang Yang
- Beijing National Laboratory for Molecular Science, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Hong-Guang Xu
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Umar Farooq
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Xi-Ling Xu
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Yi-Qin Gao
- Beijing National Laboratory for Molecular Science, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Wei-Jun Zheng
- Beijing National Laboratory for Molecular Sciences (BNLMS), State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
10
|
Okamoto R, Koga K, Onuki A. Theory of electrolytes including steric, attractive, and hydration interactions. J Chem Phys 2020; 153:074503. [PMID: 32828079 DOI: 10.1063/5.0015446] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
We present a continuum theory of electrolytes composed of a waterlike solvent and univalent ions. First, we start with a density functional F for the coarse-grained solvent, cation, and anion densities, including the Debye-Hückel free energy, the Coulombic interaction, and the direct interactions among these three components. These densities fluctuate obeying the distribution ∝exp(-F/kBT). Eliminating the solvent density deviation in F, we obtain the effective non-Coulombic interactions among the ions, which consist of the direct ones and the solvent-mediated ones. We then derive general expressions for the ion correlation, the apparent partial volume, and the activity and osmotic coefficients up to linear order in the average salt density ns. Second, we perform numerical analysis using the Mansoori-Carnahan-Starling-Leland model [J. Chem. Phys. 54, 1523 (1971)] for three-component hardspheres. The effective interactions sensitively depend on the cation and anion sizes due to competition between the steric and hydration effects, which are repulsive between small-large ion pairs and attractive between symmetric pairs. These agree with previous experiments and Collins' rule [Biophys. J. 72, 65 (1997)]. We also give simple approximate expressions for the ionic interaction coefficients valid for any ion sizes.
Collapse
Affiliation(s)
- Ryuichi Okamoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Kenichiro Koga
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Akira Onuki
- Department of Physics, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
11
|
Rubinson KA, Mountain RD. Ion and water transport reasonably involves rotation and pseudorotation: measurement and modeling the temperature dependence of small-angle neutron scattering from aqueous SrI 2. Phys Chem Chem Phys 2020; 22:13479-13488. [PMID: 32525150 DOI: 10.1039/d0cp02088e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
X-ray and neutron scattering have provided insight into the short range (<8 Å) structures of ionic solutions for over a century. For longer distances, single scattering bands have, however, been seen. For the non-hydrolyzing salt SrI2 in aqueous (D2O) solution, a structure sufficient to scatter slow neutrons has been seen to persist down to a concentration of 0.1 mol L-1 where the measured average spacing between scatterers is over 20 Å. Theoretical studies of such long distance solution structures are difficult, and these difficulties are discussed. The width of the distribution in distances between the scatterers (ions, ion pairs, etc.) remains less than 10 Å, which approximates the average size of the ions and their first hydration shell. Here, we measure the temperature dependence from 10 °C to 90 °C of the small angle neutron scattering (SANS) by a 0.5 molar SrI2 solution in D2O and find that this surprisingly narrow distribution of the distances remains constant within experimental uncertainty. This structure of the ions in the solution appears to endure because changes in interion distances along any single spatial dimension require displacements near the size of a water molecule. Together, the experimental measurements support a rotatory mechanism for simultaneous ion transport and water countertransport. Since rotation minimizes displacement of the solution framework, it is suggested that water transport alone also involves rotation of multimolecular structures, and that the interpretation of single-molecule water rotation is confounded by pseudorotation that results from paired picosecond proton exchanges. It is pointed out that NMR-determined millisecond to microsecond proton exchange times of chelated-metal-ion bound waters and the much faster chelate rotational correlation times around 10 picoseconds, both of which require making and breaking of hydrogen bonds, are difficult to impossible to reconcile.
Collapse
Affiliation(s)
- Kenneth A Rubinson
- Department of Biochemistry and Molecular Biology, Wright State University, Dayton, OH 45435, USA. and NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA.
| | - Raymond D Mountain
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA.
| |
Collapse
|
12
|
Fetisov EO, Mundy CJ, Schenter GK, Benmore CJ, Fulton JL, Kathmann SM. Nanometer-Scale Correlations in Aqueous Salt Solutions. J Phys Chem Lett 2020; 11:2598-2604. [PMID: 32163289 DOI: 10.1021/acs.jpclett.0c00322] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The intermediate-range (1-4 nm) correlation of cations, anions, and water in aqueous alkaline earth salt solutions is measured using synchrotron X-ray diffraction. We differentiate from the entire solution structure factor, SX(Q), a separate region at low Q (<1.5 Å-1) containing local diffraction maxima (prepeaks) that indicate nanometer-scale oscillatory behavior. These features are quantitatively reproduced with classical molecular dynamics simulations. At high concentrations, the prepeaks emerge from correlations arising from the existence of small quasi close-packed lattice-like structures comprised of cation hydration spheres. We also analyze the concentration dependence of the prepeak and discuss the overall results in light of the rich literature dealing with intermediate-range correlations underlying universal phenomena in concentrated electrolytes.
Collapse
Affiliation(s)
- Evgenii O Fetisov
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christopher J Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Gregory K Schenter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Chris J Benmore
- Advanced Photon Source, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - John L Fulton
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Shawn M Kathmann
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
13
|
Polák J, Ondo D, Heyda J. Thermodynamics of N-Isopropylacrylamide in Water: Insight from Experiments, Simulations, and Kirkwood-Buff Analysis Teamwork. J Phys Chem B 2020; 124:2495-2504. [PMID: 32118440 DOI: 10.1021/acs.jpcb.0c00413] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The behavior of thermoresponsive polymer poly(N-isopropylacrylamide) (PNiPAM), an essential building block in the design of smart soft materials, in aqueous solutions has attracted much interest, which contrasts with our knowledge of N-isopropylacrylamide (NiPAM) monomer. Strikingly, the physicochemical properties of aqueous NiPAM are similarly rich, and their understanding is far from being complete. This stems from the lack of accurate thermodynamic data and quantitative model for atomistic simulations. In this joint study, we have probed the thermodynamic behavior of aqueous NiPAM by experimental methods, molecular dynamics (MD) simulations, and Kirkwood-Buff (KB) analysis at ambient conditions. From the partial molar volumes and simultaneously correlated osmotic coefficients, with excess partial molar enthalpies of NiPAM in water, the concentration and temperature dependence of KB integrals was determined. For the purpose of this work, we have developed and employed a novel NiPAM force field, which not only reproduces KB integrals (Gij) and adequately captures macroscopic thermodynamic quantities but also provides more accurate structural insight than the original force fields. We revealed in the vicinity of NiPAM the competing effect of amide hydration with interaction between nonpolar regions. This microscopic picture is reflected in the experimentally observed NiPAM-NiPAM association, which is present from highly dilute conditions up to the solubility limit and is evidenced by G22. From intermediate concentrations, it is accompanied by the existence of apparent dense-water regions, as indicated by positive G11 values. The here-employed KB-based framework provided a mutually consistent thermodynamic and microscopic insight into the NiPAM solution and may be further extended for ion-specific effects. Moreover, our findings contribute to the understanding of thermodynamic grounds behind PNiPAM collapse transition.
Collapse
Affiliation(s)
- Jakub Polák
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| | - Daniel Ondo
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| | - Jan Heyda
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| |
Collapse
|
14
|
Tomaník L, Muchová E, Slavíček P. Solvation energies of ions with ensemble cluster-continuum approach. Phys Chem Chem Phys 2020; 22:22357-22368. [DOI: 10.1039/d0cp02768e] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
An alternative cluster-continuum approach for the calculation of solvation free energies of ions.
Collapse
Affiliation(s)
- Lukáš Tomaník
- Department of Physical Chemistry
- University of Chemistry and TechnologyTechnická 5
- 16628 Prague 6
- Czech Republic
| | - Eva Muchová
- Department of Physical Chemistry
- University of Chemistry and TechnologyTechnická 5
- 16628 Prague 6
- Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry
- University of Chemistry and TechnologyTechnická 5
- 16628 Prague 6
- Czech Republic
| |
Collapse
|
15
|
Chialvo AA. On the Solute-Induced Structure-Making/Breaking Effect: Rigorous Links among Microscopic Behavior, Solvation Properties, and Solution Non-Ideality. J Phys Chem B 2019; 123:2930-2947. [PMID: 30794414 DOI: 10.1021/acs.jpcb.9b00364] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
We studied the solute-induced perturbation of the solvent environment around a solute species from a microscopic viewpoint and propose a novel approach to the understanding of the structure-making/breaking process, regardless of the type and nature of the solute-solvent interactions. Based on the Kirkwood-Buff fluctuation formalism, we present a rigorous statistical mechanics description of the evolution of the solvent structure around the solute, analyze its response to small perturbations of the ( TP) state conditions and composition of the system, and make direct connections between a few equivalent micro- and macroscopic manifestations as probes for, and targets of, experimental measurements. We illustrate the analysis with theoretical results from integral equation calculations of model fluids and experimental evidence from available data for a variety of aqueous electrolyte and nonelectrolyte real fluid solutions. Finally, we provide a critical discussion about the inadequacy underlying a widely used de facto criterion for the classification of structure-making/breaking solutes.
Collapse
|
16
|
Petris PC, Anogiannakis SD, Tzounis PN, Theodorou DN. Thermodynamic Analysis of n-Hexane-Ethanol Binary Mixtures Using the Kirkwood-Buff Theory. J Phys Chem B 2019; 123:247-257. [PMID: 30516991 DOI: 10.1021/acs.jpcb.8b10425] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A complete thermodynamic analysis of mixtures consisting of molecules with complex chemical constitution can be rather demanding. The Kirkwood-Buff theory of solutions allows the estimation of thermodynamic properties, which cannot be directly extracted from atomistic simulations, such as the Gibbs energy of mixing (Δmix G). In this work, we perform molecular dynamics simulations of n-hexane-ethanol binary mixtures in the liquid state under two temperature-pressure conditions and at various mole fractions. On the basis of the recently published methodology of Galata [ Fluid Phase Equilib. 2018 , 470 , 25 - 37 ] , we first calculate the Kirkwood-Buff integrals in the isothermal-isobaric ( NpT) ensemble, identifying how system size affects their estimation. We then extract the activity coefficients, excess Gibbs energy, excess enthalpy, and excess entropy for the n-hexane-ethanol binary mixtures we simulate. We employ two approaches for quantifying composition fluctuations: one based on counting molecular centers of mass and a second one based on counting molecular segments. Results from the two approaches are practically indistinguishable. We compare our results against predictions of vapor-liquid equilibria obtained in a previous simulation work using the same force field, as well as with experimental data, and find very good agreement. In addition, we develop a simple methodology to identify the hydrogen bonds between ethanol molecules and analyze their effects on mixing properties.
Collapse
Affiliation(s)
- Panagiotis C Petris
- School of Chemical Engineering , National Technical University of Athens , GR 15780 Athens , Greece
| | - Stefanos D Anogiannakis
- School of Chemical Engineering , National Technical University of Athens , GR 15780 Athens , Greece
| | | | - Doros N Theodorou
- School of Chemical Engineering , National Technical University of Athens , GR 15780 Athens , Greece
| |
Collapse
|