1
|
Marin TW, Janik I, Bartels DM, Chipman DM. Failure of molecular dynamics to provide appropriate structures for quantum mechanical description of the aqueous chloride ion charge-transfer-to-solvent ultraviolet spectrum. Phys Chem Chem Phys 2021; 23:9109-9120. [PMID: 33885094 DOI: 10.1039/d1cp00930c] [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
The lowest band in the charge-transfer-to-solvent ultraviolet absorption spectrum of aqueous chloride ion is studied by experiment and computation. Interestingly, the experiments indicate that at concentrations up to at least 0.25 M, where calculations indicate ion pairing to be significant, there is no notable effect of ionic strength on the spectrum. The experimental spectra are fitted to aid comparison with computations. Classical molecular dynamic simulations are carried out on dilute aqueous Cl-, Na+, and NaCl, producing radial distribution functions in reasonable agreement with experiment and, for NaCl, clear evidence of ion pairing. Clusters are extracted from the simulations for quantum mechanical excited state calculations. Accurate ab initio coupled-cluster benchmark calculations on a small number of representative clusters are carried out and used to identify and validate an efficient protocol based on time-dependent density functional theory. The latter is used to carry out quantum mechanical calculations on thousands of clusters. The resulting computed spectrum is in excellent agreement with experiment for the peak position, with little influence from ion pairing, but is in qualitative disagreement on the width, being only about half as wide. It is concluded that simulation by classical molecular dynamics fails to provide an adequate variety of structures to explain the experimental CTTS spectrum of aqueous Cl-.
Collapse
Affiliation(s)
- Timothy W Marin
- Department of Physical Sciences, Benedictine University, 5700 College Rd, Lisle, IL 60532, USA
| | | | | | | |
Collapse
|
2
|
Yamaguchi T, Koda S. Dielectrically-Consistent Stell Correction for Integral Equation Theory on Electrolyte Solutions. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2015. [DOI: 10.1246/bcsj.20150041] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Tsuyoshi Yamaguchi
- Department of Molecular Design and Engineering, Graduate School of Engineering, Nagoya University
| | - Shinobu Koda
- Department of Molecular Design and Engineering, Graduate School of Engineering, Nagoya University
| |
Collapse
|
3
|
Ratkova EL, Palmer DS, Fedorov MV. Solvation thermodynamics of organic molecules by the molecular integral equation theory: approaching chemical accuracy. Chem Rev 2015; 115:6312-56. [PMID: 26073187 DOI: 10.1021/cr5000283] [Citation(s) in RCA: 135] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Ekaterina L Ratkova
- †G. A. Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, Akademicheskaya Street 1, Ivanovo 153045, Russia.,‡The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany
| | - David S Palmer
- ‡The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany.,§Department of Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, United Kingdom
| | - Maxim V Fedorov
- ‡The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany.,∥Department of Physics, Scottish Universities Physics Alliance (SUPA), University of Strathclyde, John Anderson Building, 107 Rottenrow East, Glasgow G4 0NG, United Kingdom
| |
Collapse
|
4
|
Joung IS, Luchko T, Case DA. Simple electrolyte solutions: comparison of DRISM and molecular dynamics results for alkali halide solutions. J Chem Phys 2013; 138:044103. [PMID: 23387564 DOI: 10.1063/1.4775743] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Using the dielectrically consistent reference interaction site model (DRISM) of molecular solvation, we have calculated structural and thermodynamic information of alkali-halide salts in aqueous solution, as a function of salt concentration. The impact of varying the closure relation used with DRISM is investigated using the partial series expansion of order-n (PSE-n) family of closures, which includes the commonly used hypernetted-chain equation (HNC) and Kovalenko-Hirata closures. Results are compared to explicit molecular dynamics (MD) simulations, using the same force fields, and to experiment. The mean activity coefficients of ions predicted by DRISM agree well with experimental values at concentrations below 0.5 m, especially when using the HNC closure. As individual ion activities (and the corresponding solvation free energies) are not known from experiment, only DRISM and MD results are directly compared and found to have reasonably good agreement. The activity of water directly estimated from DRISM is nearly consistent with values derived from the DRISM ion activities and the Gibbs-Duhem equation, but the changes in the computed pressure as a function of salt concentration dominate these comparisons. Good agreement with experiment is obtained if these pressure changes are ignored. Radial distribution functions of NaCl solution at three concentrations were compared between DRISM and MD simulations. DRISM shows comparable water distribution around the cation, but water structures around the anion deviate from the MD results; this may also be related to the high pressure of the system. Despite some problems, DRISM-PSE-n is an effective tool for investigating thermodynamic properties of simple electrolytes.
Collapse
Affiliation(s)
- In Suk Joung
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey 08854, USA
| | | | | |
Collapse
|
5
|
Luchko T, Joung IS, Case DA. Integral Equation Theory of Biomolecules and Electrolytes. INNOVATIONS IN BIOMOLECULAR MODELING AND SIMULATIONS 2012. [DOI: 10.1039/9781849735049-00051] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
The so-called three-dimensional version (3D-RISM) can be used to describe the interactions of solvent components (here we treat water and ions) with a chemical or biomolecular solute of arbitrary size and shape. Here we give an overview of the current status of such models, describing some aspects of “pure” electrolytes (water plus simple ions) and of ionophores, proteins and nucleic acids in the presence of water and salts. Here we focus primarily on interactions with water and dissolved salts; as a practical matter, the discussion is mostly limited to monovalent ions, since studies of divalent ions present many difficult problems that have not yet been addressed. This is not a comprehensive review, but covers a few recent examples that illustrate current issues.
Collapse
Affiliation(s)
- Tyler Luchko
- Department of Chemistry and Chemical Biology and BioMaPS Institute Rutgers University Piscataway NJ 08854, USA
| | - In Suk Joung
- Department of Chemistry and Chemical Biology and BioMaPS Institute Rutgers University Piscataway NJ 08854, USA
| | - David A. Case
- Department of Chemistry and Chemical Biology and BioMaPS Institute Rutgers University Piscataway NJ 08854, USA
| |
Collapse
|
6
|
Molecular dynamics simulations of LiCl association and NaCl association in water by means of ABEEM/MM. ACTA ACUST UNITED AC 2008. [DOI: 10.1007/s11426-008-0129-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
7
|
Woelki S, Kohler HH, Krienke H, Schmeer G. Improvements of DRISM calculations: symmetry reduction and hybrid algorithms. Phys Chem Chem Phys 2008; 10:898-910. [DOI: 10.1039/b712306j] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
8
|
Smirnov PR, Trostin VN. Structures of the nearest surroundings of the K+, Rb+, and Cs+ ions in aqueous solutions of their salts. RUSS J GEN CHEM+ 2007. [DOI: 10.1134/s1070363207120043] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
9
|
Woelki S, Kohler HH, Krienke H. A Singlet-RISM Theory for Solid/Liquid Interfaces Part I: Uncharged Walls. J Phys Chem B 2007; 111:13386-97. [DOI: 10.1021/jp068998t] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Stefan Woelki
- Institute of Analytical Chemistry, Chemo- and Biosensors, University of Regensburg, D-93040 Regensburg, Germany, and Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Hans-Helmut Kohler
- Institute of Analytical Chemistry, Chemo- and Biosensors, University of Regensburg, D-93040 Regensburg, Germany, and Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Hartmut Krienke
- Institute of Analytical Chemistry, Chemo- and Biosensors, University of Regensburg, D-93040 Regensburg, Germany, and Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| |
Collapse
|
10
|
Smirnov PR, Trostin VN. Structure of the nearest surrounding of the Na+ ion in aqueous solutions of its salts. RUSS J GEN CHEM+ 2007. [DOI: 10.1134/s1070363207050052] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
11
|
Vatamanu J, Cann NM. Evaluation of site-site bridge diagrams for molecular fluids. J Chem Phys 2004; 121:6922-34. [PMID: 15473751 DOI: 10.1063/1.1789131] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The presence of bridge functions in formally exact integral equation theories is the primary obstacle preventing the extraction of exact fluid structure from these theories. The bridge functions are typically neglected but in many fluids their impact may be significant. Each bridge function can be subdivided into bridge diagrams, which are well defined but difficult to evaluate. The calculation of bridge diagrams for the Chandler-Silbey-Ladanyi (CSL) integral equation theory is the subject of this paper. In particular, we evaluate the diagrams required to yield an exact theory up to the first power in density [O(rho(1))] and provide algorithms that remain feasible for any molecule. Further, the bridge diagrams are evaluated and compared with the f-bond and h-bond formulations. Exact bridge diagrams are numerically evaluated for several chiral molecules, for two polar dimers, and for SPC/E water. The quality of the diagrams is assessed in two ways: First, the predicted interatomic distributions are compared with those obtained from Monte Carlo simulations. Second, the connectivity constraints are evaluated and the errors in satisfying these exact relationships are compared for the f-bond and h-bond formulations. For apolar fluids, a clear improvement in CSL theory is evident with the inclusion of O(rho(0)) and O(rho(1)) diagrams. In contrast, for polar fluids, the inclusion of bridge diagrams does not lead to improvement in the structural predictions.
Collapse
Affiliation(s)
- J Vatamanu
- Department of Chemistry, Queen's University Kingston, Ontario K7L3N6, Canada
| | | |
Collapse
|
12
|
Zhang Z, Duan Z. Lithium chloride ionic association in dilute aqueous solution: a constrained molecular dynamics study. Chem Phys 2004. [DOI: 10.1016/j.chemphys.2003.10.030] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
13
|
Bagatella-Flores N, González-Mozuelos P. Effective pair potentials in charged molecular systems: Water mediated interactions. J Chem Phys 2002. [DOI: 10.1063/1.1503333] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
|
14
|
Fedotova M, Oparin R, Trostin V. Structure formation of aqueous electrolyte solutions under extreme conditions by the extended RISM-approach. A possibility of predicting. J Mol Liq 2001. [DOI: 10.1016/s0167-7322(01)00154-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
15
|
Kovalenko A, Truong TN. Thermochemistry of solvation: A self-consistent three-dimensional reference interaction site model approach. J Chem Phys 2000. [DOI: 10.1063/1.1313388] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
|
16
|
Kovalenko A, Hirata F. Hydration free energy of hydrophobic solutes studied by a reference interaction site model with a repulsive bridge correction and a thermodynamic perturbation method. J Chem Phys 2000. [DOI: 10.1063/1.1305885] [Citation(s) in RCA: 141] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
|
17
|
Kovalenko A, Hirata F. Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach. J Chem Phys 2000. [DOI: 10.1063/1.481676] [Citation(s) in RCA: 233] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
|
18
|
Imai T, Kinoshita M, Hirata F. Salt Effect on Stability and Solvation Structure of Peptide: An Integral Equation Study. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2000. [DOI: 10.1246/bcsj.73.1113] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
19
|
|
20
|
Kovalenko A, Hirata F. Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model. J Chem Phys 1999. [DOI: 10.1063/1.478883] [Citation(s) in RCA: 533] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Georgalis Y, Umbach P, Saenger W, Ihmels B, Soumpasis DM. Ordering of Fractal Clusters in Crystallizing Lysozyme Solutions. J Am Chem Soc 1999. [DOI: 10.1021/ja982407y] [Citation(s) in RCA: 48] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yannis Georgalis
- Contribution from the Institut für Kristallographie, Freie Universität Berlin, Takustrasse 6, 14195 Berlin, Germany, and Max-Planck Institut fur Biophysikalische Chemie, Biocomputation Group, Postfach 2841, 37018 Göttingen, Germany
| | - Patrick Umbach
- Contribution from the Institut für Kristallographie, Freie Universität Berlin, Takustrasse 6, 14195 Berlin, Germany, and Max-Planck Institut fur Biophysikalische Chemie, Biocomputation Group, Postfach 2841, 37018 Göttingen, Germany
| | - Wolfram Saenger
- Contribution from the Institut für Kristallographie, Freie Universität Berlin, Takustrasse 6, 14195 Berlin, Germany, and Max-Planck Institut fur Biophysikalische Chemie, Biocomputation Group, Postfach 2841, 37018 Göttingen, Germany
| | - Bernd Ihmels
- Contribution from the Institut für Kristallographie, Freie Universität Berlin, Takustrasse 6, 14195 Berlin, Germany, and Max-Planck Institut fur Biophysikalische Chemie, Biocomputation Group, Postfach 2841, 37018 Göttingen, Germany
| | - Dikeos Mario Soumpasis
- Contribution from the Institut für Kristallographie, Freie Universität Berlin, Takustrasse 6, 14195 Berlin, Germany, and Max-Planck Institut fur Biophysikalische Chemie, Biocomputation Group, Postfach 2841, 37018 Göttingen, Germany
| |
Collapse
|
22
|
Kameda Y, Sugawara K, Usuki T, Uemura O. Hydration Structure of Na+in Concentrated Aqueous Solutions. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 1998. [DOI: 10.1246/bcsj.71.2769] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
23
|
Kinoshita M, Hirata F. Analysis of salt effects on solubility of noble gases in water using the reference interaction site model theory. J Chem Phys 1997. [DOI: 10.1063/1.473519] [Citation(s) in RCA: 72] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
|
25
|
Llano‐Restrepo M, Chapman WG. Monte Carlo simulation of the structural properties of concentrated aqueous alkali halide solutions at 25 °C using a simple civilized model. J Chem Phys 1994. [DOI: 10.1063/1.466777] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
26
|
Hummer G, Soumpasis DM, Neumann M. Computer simulations do not support Cl-Cl pairing in aqueous NaCl solution. Mol Phys 1994. [DOI: 10.1080/00268979400100771] [Citation(s) in RCA: 48] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
27
|
Hummer G, Soumpasis DM. Correlations and free energies in restricted primitive model descriptions of electrolytes. J Chem Phys 1993. [DOI: 10.1063/1.464600] [Citation(s) in RCA: 50] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Soumpasis DM. Formal Aspects of the Potential of Mean Force Approach. COMPUTATION OF BIOMOLECULAR STRUCTURES 1993:223-239. [DOI: 10.1007/978-3-642-77798-1_16] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
|
29
|
Perkyns J, Pettitt BM. A site–site theory for finite concentration saline solutions. J Chem Phys 1992. [DOI: 10.1063/1.463485] [Citation(s) in RCA: 249] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
|