1
|
Blazquez S, de Lucas M, Vega C, Troncoso J, Gámez F. The temperature of maximum in density of aqueous solutions of nitrate and ammonium salts: Testing the Madrid-2019 force field. J Chem Phys 2024; 161:046103. [PMID: 39056391 DOI: 10.1063/5.0217827] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Accepted: 07/08/2024] [Indexed: 07/28/2024] Open
Abstract
The shift in the temperature of maximum in density (TMD) at room pressure of aqueous solutions of a set of five salts containing NO3- and/or NH4+ groups is studied both through experiments and through molecular dynamics simulations using the Madrid-2019 force field for ions and the TIP4P/2005 model for water. The experiments demonstrate the potential transferability and limitations of the Madrid-2019 force field for nitrate and ammonium ions recently developed by our group at different temperatures and add updated information to the reported datasets of TMDs for strong electrolytes. By using the Despretz law, individual ion contributions are extracted for predictive purposes from the experimental values of the shift in the TMD. Interesting findings for the behavior of the shift in the TMD in nitrate salts expose that this property might be particularly challenging for modelization approaches when dealing with polyatomic species.
Collapse
Affiliation(s)
- S Blazquez
- Departamento de Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, E 28040 Madrid, Spain
| | - M de Lucas
- Departamento de Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, E 28040 Madrid, Spain
| | - C Vega
- Departamento de Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, E 28040 Madrid, Spain
| | - J Troncoso
- Departamento de Física Aplicada, Universidade de Vigo, Escola de Enxeñaría Aeronaútica e do Espazo, E 32004 Ourense, Spain
| | - F Gámez
- Departamento de Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, E 28040 Madrid, Spain
| |
Collapse
|
2
|
Sedano LF, Blazquez S, Vega C. Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant? J Chem Phys 2024; 161:044505. [PMID: 39046346 DOI: 10.1063/5.0211871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2024] [Accepted: 06/30/2024] [Indexed: 07/25/2024] Open
Abstract
The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.
Collapse
Affiliation(s)
- L F Sedano
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - S Blazquez
- 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
|
3
|
Cruces Chamorro V, Jungwirth P, Martinez-Seara H. Building Water Models Compatible with Charge Scaling Molecular Dynamics. J Phys Chem Lett 2024; 15:2922-2928. [PMID: 38451169 PMCID: PMC10945568 DOI: 10.1021/acs.jpclett.4c00344] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Revised: 02/29/2024] [Accepted: 03/04/2024] [Indexed: 03/08/2024]
Abstract
Charge scaling has proven to be an efficient way to account in a mean-field manner for electronic polarization by aqueous ions in force field molecular dynamics simulations. However, commonly used water models with dielectric constants over 50 are not consistent with this approach leading to "overscaling", i.e., generally too weak ion-ion interactions. Here, we build water models fully compatible with charge scaling, i.e., having the correct low-frequency dielectric constant of about 45. To this end, we employ advanced optimization and machine learning schemes in order to explore the vast parameter space of four-site water models efficiently. As an a priori unwarranted positive result, we find a sizable range of force field parameters that satisfy the above dielectric constant constraint providing at the same time accuracy with respect to experimental data comparable with the best existing four-site water models such as TIP4P/2005, TIP4P-FB, or OPC. The present results thus open the way to the development of a consistent charge scaling force field for modeling ions in aqueous solutions.
Collapse
Affiliation(s)
- Victor Cruces Chamorro
- Institute of Organic Chemistry and
Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Pavel Jungwirth
- Institute of Organic Chemistry and
Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| | - Hector Martinez-Seara
- Institute of Organic Chemistry and
Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
| |
Collapse
|
4
|
Troncoso J, González-Salgado D. The temperature of maximum density for aqueous solutions. J Chem Phys 2024; 160:100902. [PMID: 38465676 DOI: 10.1063/5.0180094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2023] [Accepted: 12/26/2023] [Indexed: 03/12/2024] Open
Abstract
Experimental and theoretical advances for understanding the temperature of maximum density (TMD) of aqueous solutions are outlined. The main equations that relate the TMD behavior to key thermodynamic properties are stated. The experimental TMD data are classified as a function of the nature of the solute (inorganic electrolytes, non-electrolytes, organic salts and ionic liquids, and amino acids and proteins). In addition, the experimental results that explore the effect of pressure are detailed. These experimental data are rationalized by making use of qualitative and semi-quantitative arguments based on the thermodynamics of aqueous systems. The main theoretical and simulation advances in TMD for aqueous solutions are also shown-including new calculations in the context of the scaled particle theory-and their ability to reproduce the experimental data is evaluated. Finally, new experiments and theoretical and simulation developments, which could give important insights into the problem of TMD for aqueous solutions, are proposed.
Collapse
Affiliation(s)
- Jacobo Troncoso
- Instituto de Física e Ciencias Aeroespaciais da Universidade de Vigo and Unidad MSMN Asociada al CSIC por el IQF Blas Cabrera, Ourense 32004, Spain
| | - Diego González-Salgado
- Instituto de Física e Ciencias Aeroespaciais da Universidade de Vigo and Unidad MSMN Asociada al CSIC por el IQF Blas Cabrera, Ourense 32004, Spain
| |
Collapse
|
5
|
Blazquez S, Bourg IC, Vega C. Madrid-2019 force field: An extension to divalent cations Sr2+ and Ba2. J Chem Phys 2024; 160:046101. [PMID: 38270237 DOI: 10.1063/5.0186233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2023] [Accepted: 12/27/2023] [Indexed: 01/26/2024] Open
Abstract
In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Ian C Bourg
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, USA
- High Meadows Environmental Institute, Princeton University, Princeton, New Jersey 08544, USA
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
6
|
Blazquez S, Abascal JLF, Lagerweij J, Habibi P, Dey P, Vlugt TJH, Moultos OA, Vega C. Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property. J Chem Theory Comput 2023; 19:5380-5393. [PMID: 37506381 PMCID: PMC10448725 DOI: 10.1021/acs.jctc.3c00562] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Indexed: 07/30/2023]
Abstract
In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.
Collapse
Affiliation(s)
- Samuel Blazquez
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jose L. F. Abascal
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jelle Lagerweij
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Parsa Habibi
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD Delft, The Netherlands
| | - Poulumi Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD Delft, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Carlos Vega
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
7
|
Perin L, Gallo P. Phase Diagram of Aqueous Solutions of LiCl: a Study of Concentration Effects on the Anomalies of Water. J Phys Chem B 2023; 127:4613-4622. [PMID: 37167579 DOI: 10.1021/acs.jpcb.3c00703] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
We perform molecular dynamics simulations in order to study thermodynamics and the structure of supercooled aqueous solutions of lithium chloride (LiCl) at concentrations c = 0.678 and 2.034 mol/kg. We model the solvent using the TIP4P/2005 potential and the ions using the Madrid-2019 force field, a force field particularly suited for studying this solution. We find that, for c = 0.678 mol/kg, the behavior of the equation of state, studied in the P-T plane, indicates the presence of a liquid-liquid phase transition, similar to what was previously found for bulk water. We estimate the position of the liquid-liquid critical point to be at Tc ≈ 174 K, Pc ≈ 1775 bar, and ρc ≈ 1.065 g/cm3. When the concentration is tripled to c = 2.034 mol/kg, no critical point is observed, indicating its possible disappearance at this concentration. We also study the water-water and water-ions structure in the two solutions, and we find that at the concentrations examined the effect of ions on the water-water structure is not strong, and all the features found in bulk water are preserved. We also calculate the hydration number of the Li and Cl ions, and in line with experiments, we find the value of 4 for Li+ and between 5.5 and 6 for Cl-, confirming the good performances of the Madrid-2019 force field.
Collapse
Affiliation(s)
- Leonardo Perin
- Dipartimento di Fisica, Università Roma Tre, Via della Vasca Navale 84, I-00146 Roma, Italy
| | - Paola Gallo
- Dipartimento di Fisica, Università Roma Tre, Via della Vasca Navale 84, I-00146 Roma, Italy
| |
Collapse
|
8
|
Blazquez S, Conde MM, Vega C. Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water. J Chem Phys 2023; 158:054505. [PMID: 36754806 DOI: 10.1063/5.0136498] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na+ and Cl- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
9
|
Building a Hofmeister-like series for the maximum in density temperature of aqueous electrolyte solutions. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
|
10
|
Cerdeiriña CA. Water's Unusual Thermodynamics in the Realm of Physical Chemistry. J Phys Chem B 2022; 126:6608-6613. [PMID: 36001372 PMCID: PMC9797112 DOI: 10.1021/acs.jpcb.2c05274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 08/10/2022] [Indexed: 12/31/2022]
Abstract
While it is known since the early work by Edsall, Frank and Evans, Kauzmann, and others that the thermodynamics of solvation of nonpolar solutes in water is unusual and has implications for the thermodynamics of protein folding, only recently have its connections with the unusual temperature dependence of the density of solvent water been illuminated. Such density behavior is, in turn, one of the manifestations of a nonstandard thermodynamic pattern contemplating a second, liquid-liquid critical point at conditions of temperature and pressure at which water exists as a deeply supercooled liquid. Recent experimental and computational work unambiguously points toward the existence of such a critical point, thereby providing concrete answers to the questions posed by the 1976 pioneering experiments by Speedy and Angell and the associated "liquid-liquid transition hypothesis" posited in 1992 by Stanley and co-workers. Challenges of this phenomenology to the branch of Statistical Mechanics remain.
Collapse
Affiliation(s)
- Claudio A. Cerdeiriña
- Departamento de Física Aplicada, Universidad de Vigo—Campus del Agua, Ourense 32004, Spain
| |
Collapse
|