1
|
Sun Y, Zhang D, Bashir A, Li C, Fan Z. Scaling Solute-Solvent Distances to Improve Solubility and Ion Paring Predictions in Rigid Ion Models. J Phys Chem B 2023; 127:9575-9586. [PMID: 37906589 DOI: 10.1021/acs.jpcb.3c05993] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
Force fields based on the rigid ion model (RIM) have been developed to accurately predict the various physical and chemical properties of salts and water. However, the combined use of these models often fails to accurately predict the solubility of salts in water. To address this issue, several approaches, such as charge scaling or reparameterization, have been proposed. Nevertheless, these methods require laborious reparameterization of nonbonded force field parameters. In this article, we propose a scaling solute-solvent distance (SSSD) method to improve force fields in predicting salt solubility without changing the solute-solute and solvent-solvent interactions in the original force fields. This method can also tune the ion pairing of salt in water. One main advantage of the SSSD method is that reparameterization of the crystal and water models is not needed. We use two RIMs for the NaCl-water system (JC-SPC/E and SD-SPC/E) and the CHARMM force field for the KCl-water system to demonstrate the improved accuracy in predicting solubility by the SSSD method. Furthermore, we use the RDG-SPC/Fw force field to show that the SSSD method can also be used to tune the ion pairing of CaCO3 in water. Limitations of this method are also discussed.
Collapse
Affiliation(s)
- Ying Sun
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
| | - Dan Zhang
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
- Department of Chemistry, Shanghai University, Shanghai 200444, P. R. China
| | - Ayesha Bashir
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
- School of Nano-Tech and Nano-Bionics, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Chuncheng Li
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
| | - Zhaochuan Fan
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
- School of Nano-Tech and Nano-Bionics, University of Science and Technology of China, Hefei 230026, P. R. 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
|
Bernhardt MP, Nagata Y, van der Vegt NFA. Where Lennard-Jones Potentials Fail: Iterative Optimization of Ion-Water Pair Potentials Based on Ab Initio Molecular Dynamics Data. J Phys Chem Lett 2022; 13:3712-3717. [PMID: 35439420 DOI: 10.1021/acs.jpclett.2c00121] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The use of the Lennard-Jones (LJ) potential in computer simulations of aqueous electrolyte solutions is widespread. The standard approach is to parametrize LJ potential parameters against thermodynamic solution properties, but problems in representing the local structural and dynamic properties of ion hydration shells remain. The r-12-term in the LJ potential is responsible for this as it leads to overly repulsive ion-water interactions at short range. As a result, the LJ potential predicts blue-shifted vibrational peaks of the cations' rattling mode and too large negative ion hydration entropies. We demonstrate that cation-water effective pair potentials derived from ab initio MD data have softer short-range repulsions and represent hydration shell properties significantly better. Our findings indicate that replacing the LJ potential with these effective pair potentials offers a promising route to represent thermodynamic solution properties and local interactions of specific ions with nonpolarizable force field models.
Collapse
Affiliation(s)
- Marvin P Bernhardt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Straße 10, 64287 Darmstadt, Germany
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Nico F A van der Vegt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Straße 10, 64287 Darmstadt, Germany
| |
Collapse
|
4
|
Fu Y, Su S, Zhang N, Wang Y, Guo X, Xue J. Dehydration-Determined Ion Selectivity of Graphene Subnanopores. ACS APPLIED MATERIALS & INTERFACES 2020; 12:24281-24288. [PMID: 32349478 DOI: 10.1021/acsami.0c03932] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
Graphene membranes with subnanopores are considered to be the next-generation materials for water desalination and ion separation, while their performance is mainly determined by the relative ion selectivity of the pores. However, the origin of this phenomenon has been controversial in the past few years, which strongly limits the development of related applications. Here, using direct Au ion bombardment, we fabricated the desired subnanopores with average diameters of 0.8 ± 0.16 nm in monolayer graphene. The pores showed the ability to sieve K+, Na+, Li+, Cs+, Mg2+, and Ca2+ cations, and the observed K+/Mg2+ selectivity ratio was over 4. With further molecular dynamics simulations, we demonstrated that the ion selectivity is primarily attributed to the dehydration process of ions that can be quantitatively described by the ion-dependent free-energy barriers. Hopefully, this work is helpful in further enhancing the ion selectivity of graphene nanopores and also presenting a new paradigm for improving the performance of other nanoporous atomically thin membranes, such as MXenes and MoS2.
Collapse
Affiliation(s)
- Yanjun Fu
- State Key Laboratory of Nuclear Physics and Technology, School of Physics, Peking University, Beijing 100871, P. R. China
- CAPT, HEDPS, College of Engineering, Peking University, Beijing 100871, P. R. China
| | - Shihao Su
- State Key Laboratory of Nuclear Physics and Technology, School of Physics, Peking University, Beijing 100871, P. R. China
- CAPT, HEDPS, College of Engineering, Peking University, Beijing 100871, P. R. China
| | - Ning Zhang
- State Key Laboratory of Nuclear Physics and Technology, School of Physics, Peking University, Beijing 100871, P. R. China
- CAPT, HEDPS, College of Engineering, Peking University, Beijing 100871, P. R. China
| | - Yihan Wang
- State Key Laboratory of Nuclear Physics and Technology, School of Physics, Peking University, Beijing 100871, P. R. China
- CAPT, HEDPS, College of Engineering, Peking University, Beijing 100871, P. R. China
| | - Xun Guo
- State Key Laboratory of Nuclear Physics and Technology, School of Physics, Peking University, Beijing 100871, P. R. China
| | - Jianming Xue
- State Key Laboratory of Nuclear Physics and Technology, School of Physics, Peking University, Beijing 100871, P. R. China
- CAPT, HEDPS, College of Engineering, Peking University, Beijing 100871, P. R. China
| |
Collapse
|
5
|
Saric D, Kohns M, Vrabec J. Dielectric constant and density of aqueous alkali halide solutions by molecular dynamics: A force field assessment. J Chem Phys 2020; 152:164502. [PMID: 32357782 DOI: 10.1063/1.5144991] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The concentration dependence of the dielectric constant and the density of 11 aqueous alkali halide solutions (LiCl, NaCl, KCl, RbCl, CsCl, LiI, NaI, KI, CsI, KF, and CsF) is investigated by molecular simulation. Predictions using eight non-polarizable ion force fields combined with the TIP4P/ε water model are compared to experimental data. The influence of the water model and the temperature on the results for the NaCl brine are also addressed. The TIP4P/ε water model improves the accuracy of dielectric constant predictions compared to the SPC/E water model. The solution density is predicted well by most ion models. Almost all ion force fields qualitatively capture the decline of the dielectric constant with the increase of concentration for all solutions and with the increase of temperature for NaCl brine. However, the sampled dielectric constant is mostly in poor quantitative agreement with experimental data. These results are related to the microscopic solution structure, ion pairing, and ultimately the force field parameters. Ion force fields with excessive contact ion pairing and precipitation below the experimental solubility limit generally yield higher dielectric constant values. An adequate reproduction of the experimental solubility limit should therefore be a prerequisite for further investigations of the dielectric constant of aqueous electrolyte solutions by molecular simulation.
Collapse
Affiliation(s)
- Denis Saric
- Thermodynamics and Energy Technology, University of Paderborn, 33098 Paderborn, Germany
| | - Maximilian Kohns
- Laboratory of Engineering Thermodynamics, Technische Universität Kaiserslautern, 67633 Kaiserslautern, Germany
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technical University Berlin, 10587 Berlin, Germany
| |
Collapse
|
6
|
Döpke MF, Moultos OA, Hartkamp R. On the transferability of ion parameters to the TIP4P/2005 water model using molecular dynamics simulations. J Chem Phys 2020; 152:024501. [PMID: 31941316 DOI: 10.1063/1.5124448] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Countless molecular dynamics studies have relied on available ion and water force field parameters to model aqueous electrolyte solutions. The TIP4P/2005 model has proven itself to be among the best rigid water force fields, whereas many of the most successful ion parameters were optimized in combination with SPC/E, TIP3P, or TIP4P/Ew water. Many researchers have combined these ions with TIP4P/2005, hoping to leverage the strengths of both parameter sets. To assess if this widely used approach is justified and to provide a guide in selecting ion parameters, we investigated the transferability of various commonly used monovalent and multivalent ion parameters to the TIP4P/2005 water model. The transferability is evaluated in terms of ion hydration free energy, hydration radius, coordination number, and self-diffusion coefficient at infinite dilution. For selected ion parameters, we also investigated density, ion pairing, chemical potential, and mean ionic activity coefficients at finite concentrations. We found that not all ions are equally transferable to TIP4P/2005 without compromising their performance. In particular, ions optimized for TIP3P water were found to be poorly transferable to TIP4P/2005, whereas ions optimized for TIP4P/Ew water provided nearly perfect transferability. The latter ions also showed good overall agreement with experimental values. The one exception is that no combination of ion parameters and water model considered here was found to accurately reproduce experimental self-diffusion coefficients. Additionally, we found that cations optimized for SPC/E and TIP3P water displayed consistent underpredictions in the hydration free energy, whereas anions consistently overpredicted the hydration free energy.
Collapse
Affiliation(s)
- Max F Döpke
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Othonas A Moultos
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Remco Hartkamp
- Process & Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| |
Collapse
|
7
|
Tolmachev DA, Boyko OS, Lukasheva NV, Martinez-Seara H, Karttunen M. Overbinding and Qualitative and Quantitative Changes Caused by Simple Na+ and K+ Ions in Polyelectrolyte Simulations: Comparison of Force Fields with and without NBFIX and ECC Corrections. J Chem Theory Comput 2019; 16:677-687. [DOI: 10.1021/acs.jctc.9b00813] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- D. A. Tolmachev
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| | - O. S. Boyko
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| | - N. V. Lukasheva
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| | - H. Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 542/2, Prague 6 CZ166 10, Czech Republic
| | - Mikko Karttunen
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, St. Petersburg 199004, Russia
| |
Collapse
|
8
|
Zeron IM, Abascal JLF, Vega C. A force field of Li +, Na +, K +, Mg 2+, Ca 2+, Cl -, and SO 4 2- in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions. J Chem Phys 2019; 151:134504. [PMID: 31594349 DOI: 10.1063/1.5121392] [Citation(s) in RCA: 123] [Impact Index Per Article: 24.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
In this work, a force field for several ions in water is proposed. In particular, we consider the cations Li+, Na+, K+, Mg2+, and Ca2+ and the anions Cl- and SO4 2-. These ions were selected as they appear in the composition of seawater, and they are also found in biological systems. The force field proposed (denoted as Madrid-2019) is nonpolarizable, and both water molecules and sulfate anions are rigid. For water, we use the TIP4P/2005 model. The main idea behind this work is to further explore the possibility of using scaled charges for describing ionic solutions. Monovalent and divalent ions are modeled using charges of 0.85 and 1.7, respectively (in electron units). The model allows a very accurate description of the densities of the solutions up to high concentrations. It also gives good predictions of viscosities up to 3 m concentrations. Calculated structural properties are also in reasonable agreement with the experiment. We have checked that no crystallization occurred in the simulations at concentrations similar to the solubility limit. A test for ternary mixtures shows that the force field provides excellent performance at an affordable computer cost. In summary, the use of scaled charges, which could be regarded as an effective and simple way of accounting for polarization (at least to a certain extend), improves the overall description of ionic systems in water. However, for purely ionic systems, scaled charges will not adequately describe neither the solid nor the melt.
Collapse
Affiliation(s)
- I M Zeron
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J L F Abascal
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
9
|
Bachler J, Handle PH, Giovambattista N, Loerting T. Glass polymorphism and liquid-liquid phase transition in aqueous solutions: experiments and computer simulations. Phys Chem Chem Phys 2019; 21:23238-23268. [PMID: 31556899 DOI: 10.1039/c9cp02953b] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
One of the most intriguing anomalies of water is its ability to exist as distinct amorphous ice forms (glass polymorphism or polyamorphism). This resonates well with the possible first-order liquid-liquid phase transition (LLPT) in the supercooled state, where ice is the stable phase. In this Perspective, we review experiments and computer simulations that search for LLPT and polyamorphism in aqueous solutions containing salts and alcohols. Most studies on ionic solutes are devoted to NaCl and LiCl; studies on alcohols have mainly focused on glycerol. Less attention has been paid to protein solutions and hydrophobic solutes, even though they reveal promising avenues. While all solutions show polyamorphism and an LLPT only in dilute, sub-eutectic mixtures, there are differences regarding the nature of the transition. Isocompositional transitions for varying mole fractions are observed in alcohol but not in ionic solutions. This is because water can surround alcohol molecules either in a low- or high-density configuration whereas for ionic solutes, the water ion hydration shell is forced into high-density structures. Consequently, the polyamorphic transition and the LLPT are prevented near the ions, but take place in patches of water within the solutions. We highlight discrepancies and different interpretations within the experimental community as well as the key challenges that need consideration when comparing experiments and simulations. We point out where reinterpretation of past studies helps to draw a unified, consistent picture. In addition to the literature review, we provide original experimental results. A list of eleven open questions that need further consideration is identified.
Collapse
Affiliation(s)
- Johannes Bachler
- Institute of Physical Chemistry, University of Innsbruck, A-6020 Innsbruck, Austria.
| | | | | | | |
Collapse
|
10
|
Yuan R, Fayer MD. Dynamics of Water Molecules and Ions in Concentrated Lithium Chloride Solutions Probed with Ultrafast 2D IR Spectroscopy. J Phys Chem B 2019; 123:7628-7639. [DOI: 10.1021/acs.jpcb.9b06038] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Rongfeng Yuan
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Michael D. Fayer
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
11
|
Yllö A, Zhang C. Experimental and molecular dynamics study of the ionic conductivity in aqueous LiCl electrolytes. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.05.004] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
12
|
Pethes I. The structure of aqueous lithium chloride solutions at high concentrations as revealed by a comparison of classical interatomic potential models. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.05.044] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
|
13
|
Peng H, Nguyen AV. A link between viscosity and cation-anion contact pairs: Adventure on the concept of structure-making/breaking for concentrated salt solutions. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.04.145] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
14
|
Camisasca G, De Marzio M, Rovere M, Gallo P. High density liquid structure enhancement in glass forming aqueous solution of LiCl. J Chem Phys 2018; 148:222829. [DOI: 10.1063/1.5024375] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- G. Camisasca
- Dipartimento di Matematica e Fisica, Università “Roma Tre,” Via della Vasca Navale 84, 00146 Roma, Italy
| | - M. De Marzio
- Dipartimento di Matematica e Fisica, Università “Roma Tre,” Via della Vasca Navale 84, 00146 Roma, Italy
| | - M. Rovere
- Dipartimento di Matematica e Fisica, Università “Roma Tre,” Via della Vasca Navale 84, 00146 Roma, Italy
| | - P. Gallo
- Dipartimento di Matematica e Fisica, Università “Roma Tre,” Via della Vasca Navale 84, 00146 Roma, Italy
| |
Collapse
|
15
|
Kameda Y, Maeda S, Amo Y, Usuki T, Ikeda K, Otomo T. Neutron Diffraction Study on the Structure of Hydrated Li+ in Dilute Aqueous Solutions. J Phys Chem B 2018; 122:1695-1701. [DOI: 10.1021/acs.jpcb.7b12218] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yasuo Kameda
- Department
of Material and Biological Chemistry, Faculty of Science, Yamagata University, Yamagata, Yamagata 990-8560, Japan
| | - Shunya Maeda
- Department
of Material and Biological Chemistry, Faculty of Science, Yamagata University, Yamagata, Yamagata 990-8560, Japan
| | - Yuko Amo
- Department
of Material and Biological Chemistry, Faculty of Science, Yamagata University, Yamagata, Yamagata 990-8560, Japan
| | - Takeshi Usuki
- Department
of Material and Biological Chemistry, Faculty of Science, Yamagata University, Yamagata, Yamagata 990-8560, Japan
| | - Kazutaka Ikeda
- Institute
of Material Structure Science, KEK, Tsukuba, Ibaraki 305-080, Japan
| | - Toshiya Otomo
- Institute
of Material Structure Science, KEK, Tsukuba, Ibaraki 305-080, Japan
| |
Collapse
|
16
|
Conde MM, Rovere M, Gallo P. High precision determination of the melting points of water TIP4P/2005 and water TIP4P/Ice models by the direct coexistence technique. J Chem Phys 2017; 147:244506. [DOI: 10.1063/1.5008478] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- M. M. Conde
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| | - M. Rovere
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| | - P. Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy
| |
Collapse
|
17
|
Hirata M, Yagasaki T, Matsumoto M, Tanaka H. Phase Diagram of TIP4P/2005 Water at High Pressure. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:11561-11569. [PMID: 28796510 DOI: 10.1021/acs.langmuir.7b01764] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We report a new ice phase that forms spontaneously at the interface between ice VII and liquid water in molecular dynamics simulations of TIP4P/2005 water. The new phase is structurally quite similar to an ice phase originally found to be a precursor in the course of the homogeneous nucleation of ice VII in liquid water. A part of the water molecules in these ice phases can rotate easily because the number of hydrogen bonds in them is less than four, and thus they can be regarded as partial plastic phases. A rough estimate suggests that these phases are thermodynamically more stable than either ice VI or ice VII for 3 GPa < P < 18 GPa at T = 300 K. Although the partial plastic phases would be metastable states at any point in the phase diagram of real water, they might be realized experimentally with the aid of dopants and/or solid surfaces.
Collapse
Affiliation(s)
- Masanori Hirata
- Graduate School of Natural Science and Technology and ‡Research Institute for Interdisciplinary Science, Okayama University , Okayama 700-8530, Japan
| | - Takuma Yagasaki
- Graduate School of Natural Science and Technology and ‡Research Institute for Interdisciplinary Science, Okayama University , Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Graduate School of Natural Science and Technology and ‡Research Institute for Interdisciplinary Science, Okayama University , Okayama 700-8530, Japan
| | - Hideki Tanaka
- Graduate School of Natural Science and Technology and ‡Research Institute for Interdisciplinary Science, Okayama University , Okayama 700-8530, Japan
| |
Collapse
|
18
|
Benavides AL, Portillo MA, Chamorro VC, Espinosa JR, Abascal JLF, Vega C. A potential model for sodium chloride solutions based on the TIP4P/2005 water model. J Chem Phys 2017; 147:104501. [DOI: 10.1063/1.5001190] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- A. L. Benavides
- Departamento de Ingeniería Física, División de Ciencias e Ingenierías, Universidad de Guanajuato, Loma del Bosque 103, Col. Lomas del Campestre, CP 37150 León, Mexico
| | - M. A. Portillo
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - V. C. Chamorro
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. R. Espinosa
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. L. F. Abascal
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
19
|
Li RZ, Liu YY, Yang M. Microsolvation of lithium iodide dimer studied by ab initio calculations. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.06.005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
20
|
A comparison of classical interatomic potentials applied to highly concentrated aqueous lithium chloride solutions. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.07.076] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
21
|
Lanaro G, Patey GN. Crystal structures of model lithium halides in bulk phase and in clusters. J Chem Phys 2017; 146:154501. [PMID: 28433031 DOI: 10.1063/1.4979926] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
We employ lattice energy calculations and molecular dynamics simulations to compare the stability of wurtzite and rock salt crystal structures of four lithium halides (LiF, LiCl, LiBr, and LiI) modeled using the Tosi-Fumi and Joung-Cheatham potentials, which are models frequently used in simulation studies. Both infinite crystals and finite clusters are considered. For the Tosi-Fumi model, we find that all four salts prefer the wurtzite structure both at 0 K and at finite temperatures, in disagreement with experiments, where rock salt is the stable structure and wurtzite exists as a metastable state. For Joung-Cheatham potentials, rock salt is more stable for LiF and LiCl, but the wurtzite structure is preferred by LiBr and LiI. It is clear that the available lithium halide force fields need improvement to bring them into better accord with the experiment. Finite-size clusters that are more stable as rock salt in the bulk phase tend to solidify as small rock salt crystals. However, small clusters of salts that prefer the wurtzite structure as bulk crystals tend to form structures that have hexagonal motifs, but are not finite-size wurtzite crystals. We show that small wurtzite structures are unstable due to the presence of a dipole and rearrange into more stable, size-dependent structures. We also show that entropic contributions can act in favor of the wurtzite structure at higher temperatures. The possible relevance of our results for simulation studies of crystal nucleation from melts and/or aqueous solutions is discussed.
Collapse
Affiliation(s)
- G Lanaro
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| | - G N Patey
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| |
Collapse
|
22
|
Benavides A, Portillo M, Abascal J, Vega C. Estimating the solubility of 1:1 electrolyte aqueous solutions: the chemical potential difference rule. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1288939] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- A.L. Benavides
- División de Ciencias e Ingenierías, Universidad de Guanajuato, Guanajuato, Mexico
| | - M.A. Portillo
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J.L.F. Abascal
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
23
|
Unger I, Seidel R, Thürmer S, Pohl MN, Aziz EF, Cederbaum LS, Muchová E, Slavíček P, Winter B, Kryzhevoi NV. Observation of electron-transfer-mediated decay in aqueous solution. Nat Chem 2017. [DOI: 10.1038/nchem.2727] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
24
|
Benavides AL, Aragones JL, Vega C. Consensus on the solubility of NaCl in water from computer simulations using the chemical potential route. J Chem Phys 2016; 144:124504. [PMID: 27036458 DOI: 10.1063/1.4943780] [Citation(s) in RCA: 83] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
The solubility of NaCl in water is evaluated by using three force field models: Joung-Cheatham for NaCl dissolved in two different water models (SPC/E and TIP4P/2005) and Smith Dang NaCl model in SPC/E water. The methodology based on free-energy calculations [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] and [J. L. Aragones et al., J. Chem. Phys. 136, 244508 (2012)] has been used, except, that all calculations for the NaCl in solution were obtained by using molecular dynamics simulations with the GROMACS package instead of homemade MC programs. We have explored new lower molalities and made longer runs to improve the accuracy of the calculations. Exploring the low molality region allowed us to obtain an analytical expression for the chemical potential of the ions in solution as a function of molality valid for a wider range of molalities, including the infinite dilute case. These new results are in better agreement with recent estimations of the solubility obtained with other methodologies. Besides, two empirical simple rules have been obtained to have a rough estimate of the solubility of a certain model, by analyzing the ionic pairs formation as a function of molality and/or by calculating the difference between the NaCl solid chemical potential and the standard chemical potential of the salt in solution.
Collapse
Affiliation(s)
- A L Benavides
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J L Aragones
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
25
|
Chen H, Cox JR, Panagiotopoulos AZ. Force Fields for Carbohydrate-Divalent Cation Interactions. J Phys Chem B 2016; 120:5203-8. [PMID: 27210229 DOI: 10.1021/acs.jpcb.6b01438] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
We report molecular dynamics simulations to study intermolecular interactions for carbohydrate-divalent cation complexes. We observed that common force fields from literature with standard Lorentz-Berthelot combining rules are unable to reproduce the experimental stability constants for model carbohydrate monomer (α-d-Allopyranose) and alkali earth metal cation (Mg(2+), Ca(2+), Sr(2+), or Ba(2+)) complexes. A modified combining rule with rescaled effective cross-interaction radius between cations and the hydroxyl oxygens on the carbohydrates was introduced to reproduce the experimental stability constants, which the preferential carbohydrate-cation complexing structures through the ax-eq-ax sequence of O-1, O-2, and O-3 on α-d-Allopyranose were also observed. The effective radius scaling factor obtained from (α-d-Allopyranose)-Ca(2+) complexes was directly transferrable to the similar six-membered ring (α-d-Ribopyranose)-Ca(2+) complexes; however, reparameterization for the scaling factor may be necessary for the five-membered ring (α-d-Ribofuranose)-Ca(2+) complexes.
Collapse
Affiliation(s)
- Hsieh Chen
- Aramco Services Company: Aramco Research Center - Boston, Cambridge, Massachusetts 02139, United States
| | - Jason R Cox
- Aramco Services Company: Aramco Research Center - Boston, Cambridge, Massachusetts 02139, United States
| | - Athanassios Z Panagiotopoulos
- Department of Chemical and Biological Engineering, Princeton University , Princeton, New Jersey 08544, United States
| |
Collapse
|
26
|
Affiliation(s)
- Wen Jun Xie
- Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, and Biodynamic
Optical Imaging Center, Peking University, Beijing 100871, China
| | - Zhen Zhang
- Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, and Biodynamic
Optical Imaging Center, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, and Biodynamic
Optical Imaging Center, Peking University, Beijing 100871, China
| |
Collapse
|
27
|
Kohagen M, Pluhařová E, Mason PE, Jungwirth P. Exploring Ion-Ion Interactions in Aqueous Solutions by a Combination of Molecular Dynamics and Neutron Scattering. J Phys Chem Lett 2015; 6:1563-1567. [PMID: 26263314 DOI: 10.1021/acs.jpclett.5b00060] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Recent advances in computational and experimental techniques have allowed for accurate description of ion pairing in aqueous solutions. Free energy methods based on ab initio molecular dynamics, as well as on force fields accounting effectively for electronic polarization, can provide quantitative information about the structures and occurrences of individual types of ion pairs. When properly benchmarked against electronic structure calculations for model systems and against structural experiments, in particular neutron scattering, such force field simulations represent a powerful tool for elucidating interactions of salt ions in complex biological aqueous environments.
Collapse
Affiliation(s)
- Miriam Kohagen
- †Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 16610 Prague 6, Czech Republic
| | - Eva Pluhařová
- ‡Department of Chemistry, École Normale Supérieure, UMR ENS-CNRS-UPMC 8640, 24 rue Lhomond, 75005 Paris, France
| | - Philip E Mason
- †Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 16610 Prague 6, Czech Republic
| | - Pavel Jungwirth
- †Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 2, 16610 Prague 6, Czech Republic
| |
Collapse
|
28
|
|
29
|
Singh MB, Dalvi VH, Gaikar VG. Investigations of clustering of ions and diffusivity in concentrated aqueous solutions of lithium chloride by molecular dynamic simulations. RSC Adv 2015. [DOI: 10.1039/c4ra15124k] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The diffusion coefficient of Li+ ions decreases with increase in LiCl concentration which depends on the size of coordination structure of ions formed in solutions.
Collapse
Affiliation(s)
- Meena B. Singh
- Department of Chemical Engineering
- Institute of Chemical Technology
- Matunga
- India
| | - Vishwanath H. Dalvi
- Department of Chemical Engineering
- Institute of Chemical Technology
- Matunga
- India
| | - Vilas G. Gaikar
- Department of Chemical Engineering
- Institute of Chemical Technology
- Matunga
- India
| |
Collapse
|
30
|
Hartkamp R, Coasne B. Structure and transport of aqueous electrolytes: From simple halides to radionuclide ions. J Chem Phys 2014; 141:124508. [DOI: 10.1063/1.4896380] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
31
|
|