1
|
Zhang J, Searles DJ, Duignan T. A Method for Efficiently Predicting the Radial Distribution Function and Osmotic Coefficients of Aqueous Electrolyte Solutions. J Chem Theory Comput 2024. [PMID: 39099091 DOI: 10.1021/acs.jctc.4c00363] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/06/2024]
Abstract
The prediction of the structural and thermodynamic properties of electrolyte solutions is critical for a huge range of practical situations where these solutions play a vital role. Theoretical models, such as the continuum solvent model, attempt to explain the behavior of solutions using a coarse-grained description of the interactions of species in the solution, whereas molecular simulations aim to directly compute the behavior of the solution, including the interactions between all ions and molecules in the system. Both methods have limitations: theoretical models are generally less accurate because they rely on assumptions, while molecular simulations require significant computational resources, particularly if higher accuracy is desired. To address these issues, we propose an affordable and effective method that combines the advantages of the modified Poisson-Boltzmann equation (MPBE) with classical molecular dynamics (MD) simulations to predict the radial distribution functions and thermodynamic properties of electrolyte solutions. We demonstrate a method of using the MPBE to compute the short-range potential of mean force (PMF) from the radial distribution functions (RDFs) and vice versa. Furthermore, we provide insights into the relationship between the RDFs and the short-range PMF based on the MPBE. Our analysis reveals that the effective short-range PMFs can be approximately calculated using low concentration simulations but the short-range PMFs are slightly concentration-dependent in simulations at higher concentrations. Additionally, we demonstrate that for concentrated solutions, osmotic coefficients can be calculated in agreement with experiment using a virial approach. This is based on the effective short-range PMFs and RDFs obtained from the MPBE method. Our proposed MPBE can therefore accelerate the calculation of the structural and thermodynamic properties of electrolyte solutions.
Collapse
Affiliation(s)
- Junji Zhang
- School of Chemical Engineering, The University of Queensland, Brisbane, QLD 4072, Australia
| | - Debra J Searles
- Australian Institute for Bioengineering and Nanotechnology, The University of Queensland, Brisbane, QLD 4072, Australia
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, QLD 4072, Australia
| | - Timothy Duignan
- School of Chemical Engineering, The University of Queensland, Brisbane, QLD 4072, Australia
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan, QLD 4111, Australia
| |
Collapse
|
2
|
Li X, Bourg IC. Hygroscopic Growth of Adsorbed Water Films on Smectite Clay Particles. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2024; 58:1109-1118. [PMID: 38164899 PMCID: PMC10795194 DOI: 10.1021/acs.est.3c08253] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/05/2023] [Revised: 12/07/2023] [Accepted: 12/12/2023] [Indexed: 01/03/2024]
Abstract
Hygroscopic growth of adsorbed water films on clay particles underlies a number of environmental science questions, from the air quality and climate impacts of mineral dust aerosols to the hydrology and mechanics of unsaturated soils and sedimentary rocks. Here, we use molecular dynamics (MD) simulations to establish the relation between adsorbed water film thickness (h) and relative humidity (RH) or disjoining pressure (Π), which has long been uncertain due to factors including sensitivity to particle shape, surface roughness, and aqueous chemistry. We present a new MD simulation approach that enables precise quantification of Π in films up to six water monolayers thick. We find that the hygroscopicity of phyllosilicate mineral surfaces increases in the order mica < K-smectite < Na-smectite. The relationship between Π and h on clay surfaces follows a double exponential decay with e-folding lengths of 2.3 and 7.5 Å. The two decay length scales are attributed to hydration repulsion and osmotic phenomena in the electrical double layer (EDL) at the clay-water interface.
Collapse
Affiliation(s)
- Xiaohan Li
- Department
of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ian C. Bourg
- Department
of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, United States
- High
Meadows Environmental Institute, Princeton
University, Princeton, New Jersey 08544, United States
| |
Collapse
|
3
|
Reinhardt A, Chew PY, Cheng B. A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals. J Chem Phys 2023; 159:184110. [PMID: 37962445 DOI: 10.1063/5.0173341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Accepted: 10/20/2023] [Indexed: 11/15/2023] Open
Abstract
Computing the solubility of crystals in a solvent using atomistic simulations is notoriously challenging due to the complexities and convergence issues associated with free-energy methods, as well as the slow equilibration in direct-coexistence simulations. This paper introduces a molecular-dynamics workflow that simplifies and robustly computes the solubility of molecular or ionic crystals. This method is considerably more straightforward than the state-of-the-art, as we have streamlined and optimised each step of the process. Specifically, we calculate the chemical potential of the crystal using the gas-phase molecule as a reference state, and employ the S0 method to determine the concentration dependence of the chemical potential of the solute. We use this workflow to predict the solubilities of sodium chloride in water, urea polymorphs in water, and paracetamol polymorphs in both water and ethanol. Our findings indicate that the predicted solubility is sensitive to the chosen potential energy surface. Furthermore, we note that the harmonic approximation often fails for both molecular crystals and gas molecules at or above room temperature, and that the assumption of an ideal solution becomes less valid for highly soluble substances.
Collapse
Affiliation(s)
- Aleks Reinhardt
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Pin Yu Chew
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Bingqing Cheng
- Institute of Science and Technology Austria, Am Campus 1, 3400 Klosterneuburg, Austria
| |
Collapse
|
4
|
Zhang J, Clennell MB, Dewhurst DN. Transport Properties of NaCl in Aqueous Solution and Hydrogen Solubility in Brine. J Phys Chem B 2023; 127:8900-8915. [PMID: 37794729 PMCID: PMC10591480 DOI: 10.1021/acs.jpcb.3c03863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Revised: 08/17/2023] [Indexed: 10/06/2023]
Abstract
Ion transport properties and hydrogen solubility in brine play pivotal roles in various engineering and scientific scopes including chemical, physical, geochemical, and geothermal domains. Molecular dynamics simulations were performed to obtain transport properties of NaCl in the binary H2O + NaCl system using different force fields. Brine density, ion diffusivity, molar conductivity, conductivity, and hydrogen solubilities were obtained as functions of temperature and salt concentration. We compared the performance of different force fields against the experimental correlation model and developed three mathematical models. The first was the modified brine density model based on the simulated brine density over a wide range of salinity levels, and the second and third analytical mathematical models were derived for the ion diffusivity and molar conductivity as a function of salinity and temperature. The results of this study illustrated that the modified brine density model not only produced the same results of the previous model for lower salinity levels but also applied well to predict the brine density for a higher salinity level. The derived mathematical models indicated that the ion diffusivity and molar conductivity decreased linearly with salinity, and the slope and y-intercept of the lines of diffusivity and molar conductivity versus temperature were third-order polynomials of temperature. The developed models provided the mechanism for the behavior of decreasing molar conductivity with increasing salinity and increasing conductivity with increasing salinity. The directions of the effect of salinity on the molar conductivity and conductivity were opposite. The molar conductivity increased with a decreasing salinity level. However, the conductivity increased with increasing salinity, as the contribution of the ion concentration or salinity level to conductivity dominated over that of the ion movement.
Collapse
Affiliation(s)
- Junfang Zhang
- CSIRO Energy, 26 Dick Perry Avenue, Kensington, Western Australia 6151, Australia
| | - Michael B. Clennell
- CSIRO Energy, 26 Dick Perry Avenue, Kensington, Western Australia 6151, Australia
| | - David N. Dewhurst
- CSIRO Energy, 26 Dick Perry Avenue, Kensington, Western Australia 6151, Australia
| |
Collapse
|
5
|
Li M, Lv L, Fang T, Hao L, Li S, Dong S, Wu Y, Dong X, Liu H. Self-Consistent Implementation of a Solvation Free Energy Framework to Predict the Salt Solubilities of Six Alkali Halides. J Chem Theory Comput 2023; 19:5586-5601. [PMID: 37471389 DOI: 10.1021/acs.jctc.3c00083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2023]
Abstract
To assess the salt solubilities of six alkali halides in aqueous systems, we proposed a thermodynamic cycle and an efficient molecular modeling methodology. The Gibbs free energy changes for vaporization, dissociation, and dissolution were calculated using the experimental data of ionic thermodynamic properties obtained from the NBS tables. Additionally, the Marcus' and Tissandier's solvation free energy data for Li+, Na+, K+, Cl-, and Br- ions were compared with the conventional solvation free energies by substituting into our self-consistent thermodynamic cycle. Furthermore, Tissandier's absolute solvation free energy data were used as the training set to refit the Lennard-Jones parameters of OPLS-AA force field for ions. To predict salt solubilities, an assumption of a pseudo-solvent was proposed to characterize the coupling work of a solute with its environment from infinitely diluted to saturated solutions, indicating that the Gibbs energy change of solvation process is a function of ionic strength. Following the self-consistency of the cycle, the newly derived formulas were used to determine the salt solubilities by interpolating the intersection of Gibbs free energy of dissolution and the zero free energy line. The refined ion parameters can also predict the structure and thermodynamic properties of aqueous electrolyte solutions, such as densities, pair correlation functions, hydration numbers, mean activity coefficients, vapor pressures, and the radial dependences of the net charge at 298.15 K and 1 bar. Our method can be used to characterize the solid-liquid equilibria of ions or charged particles in aqueous systems. Furthermore, for highly concentrated strong electrolyte systems, it is essential to introduce accurate water models and polarizable force fields.
Collapse
Affiliation(s)
- Miyi Li
- School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 10081, China
| | - Liqiang Lv
- College of Chemical Engineering, Shijiazhuang University, Hebei, Shijiazhuang 050035, China
| | - Tao Fang
- Beijing Institute of Aerospace Testing Technology, Beijing 100074, China
| | - Long Hao
- Beijing Institute of Aerospace Testing Technology, Beijing 100074, China
| | - Shenhui Li
- School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 10081, China
| | - Shoulong Dong
- School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 10081, China
| | - Yufeng Wu
- School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 10081, China
| | - Xiao Dong
- School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 10081, China
| | - Helei Liu
- School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 10081, China
- Key Laboratory of Low-Carbon Conversion Science & Engineering, Shanghai Advanced Research Institute, Chinese Academy of Sciences (Shanghai Advanced Research Institute, Chinese Academy of Sciences), Shanghai 201210, China
| |
Collapse
|
6
|
Sanchez-Burgos I, Espinosa JR. Direct Calculation of the Interfacial Free Energy between NaCl Crystal and Its Aqueous Solution at the Solubility Limit. PHYSICAL REVIEW LETTERS 2023; 130:118001. [PMID: 37001068 DOI: 10.1103/physrevlett.130.118001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 11/07/2022] [Accepted: 01/30/2023] [Indexed: 06/19/2023]
Abstract
Salty water is the most abundant electrolyte aqueous mixture on Earth, however, very little is known about the NaCl-saturated solution interfacial free energy (γ_{s}). Here, we provide the first direct estimation of γ_{s} for several NaCl crystallographic planes by means of the mold integration technique, a highly efficient computational method to evaluate interfacial free energies with anisotropic crystal resolution. Making use of the JC-SPC/E model, one of the most benchmarked force fields for NaCl water solutions, we measure γ_{s} of four different crystal planes, (100), (110), (111), and (112[over ¯]) with the saturated solution at normal conditions. We find high anisotropy between the different crystal orientations with values ranging from 100 to 150 mJ m^{-2}, and the average value of the distinct planes being γ[over ¯]_{s}=137(20) mJ m^{-2}. This value for the coexistence interfacial free energy is in reasonable agreement with previous extrapolations from nucleation studies. Our Letter represents a milestone in the computational calculation of interfacial free energies between ionic crystals and aqueous solutions.
Collapse
Affiliation(s)
- Ignacio Sanchez-Burgos
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
7
|
Gupta S, Elliott JR, Anderko A, Crosthwaite J, Chapman WG, Lira CT. Current Practices and Continuing Needs in Thermophysical Properties for the Chemical Industry. Ind Eng Chem Res 2023. [DOI: 10.1021/acs.iecr.2c03153] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2023]
Affiliation(s)
- Sumnesh Gupta
- The Dow Chemical Company, 1254 Enclave Parkway, Houston, Texas 77077, United States
| | - J. Richard Elliott
- Chemical, Biomolecular, and Corrosion Engineering Department, University of Akron, Akron, Ohio 44325-3906, United States
| | - Andrzej Anderko
- OLI Systems, Inc., 2 Gatehall Drive, Suite 1D, Parsippany, New Jersey 07054, United States
| | - Jacob Crosthwaite
- The Dow Chemical Company, 1897 Building, Midland, Michigan 48667, United States
| | - Walter G. Chapman
- Chemical and Biomolecular Engineering Department, Rice University, 6100 Main Street, Houston, Texas 77005, United States
| | - Carl T. Lira
- Chemical Engineering & Materials Science, Michigan State University, East Lansing, Michigan 48824-2288, United States
| |
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
|
Tong J, Peng B, Kontogeorgis GM, Liang X. Behavior of the aqueous sodium chloride solutions from molecular simulations and theories. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2022.121086] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
10
|
Izarra AD, Coudert FX, Fuchs AH, Boutin A. Alchemical Osmostat for Monte Carlo Simulation: Sampling Aqueous Electrolyte Solution in Open Systems. J Phys Chem B 2023; 127:766-776. [PMID: 36634303 DOI: 10.1021/acs.jpcb.2c07902] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Molecular simulations involving electrolytes are usually performed at a fixed amount of salt ions in the simulation box, reproducing macroscopic concentration. Although this statement is valid in the bulk, the concentration of an electrolyte confined in nanoporous materials such as MOFs or zeolites is greatly affected and remains a priori unknown. The nanoporous material in equilibrium with the bulk electrolyte exchange water and ions at a given chemical potential Δμ in the semi-grand-canonical ensemble, that must be calibrated in order to determine the concentration in the nanoporous material. In this work, we propose an algorithm based on nonequilibrium candidate Monte Carlo (NCMC) moves to ultimately perform MC simulations in contact with a saline reservoir. First, we adapt the Widom insertion technique to calibrate the chemical potential by alchemically transmuting water molecules into ions by using NCMC moves. The chemical potential defines a Monte Carlo osmostat in the semi-grand-constant volume and temperature ensemble (Δμ, N, V, T) to be added in a Monte Carlo simulation where the number of ions fluctuates. In order to validate the method, we adapted the NCMC move to determine the free energy of water solvation and subsequently explore thermodynamics of electrolyte solvation at infinite dilution in water. Finally, we implemented the osmostat in MC simulations initialized with bulk water that are driven toward electrolyte solutions of similar concentration as the saline reservoir. Our results demonstrate that alchemical osmostat for MC simulation is a promising tool for use to sample electrolyte insertion in nanoporous materials.
Collapse
Affiliation(s)
- Ambroise de Izarra
- PASTEUR, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France.,Chimie ParisTech, PSL University, CNRS, Institut de Recherche de Chimie Paris, Paris75005, France
| | - François-Xavier Coudert
- Chimie ParisTech, PSL University, CNRS, Institut de Recherche de Chimie Paris, Paris75005, France
| | - Alain H Fuchs
- Chimie ParisTech, PSL University, CNRS, Institut de Recherche de Chimie Paris, Paris75005, France
| | - Anne Boutin
- PASTEUR, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France
| |
Collapse
|
11
|
Polarizable force fields for accurate molecular simulations of aqueous solutions of electrolytes, crystalline salts, and solubility: Li+, Na+, K+, Rb+, F−, Cl−, Br−, I−. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119659] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
12
|
Dawass N, Langeveld J, Ramdin M, Pérez-Gallent E, Villanueva AA, Giling EJM, Langerak J, van den Broeke LJP, Vlugt TJH, Moultos OA. Solubilities and Transport Properties of CO 2, Oxalic Acid, and Formic Acid in Mixed Solvents Composed of Deep Eutectic Solvents, Methanol, and Propylene Carbonate. J Phys Chem B 2022; 126:3572-3584. [PMID: 35507866 PMCID: PMC9125562 DOI: 10.1021/acs.jpcb.2c01425] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
![]()
Recently, deep eutectic
solvents (DES) have been considered as
possible electrolytes for the electrochemical reduction of CO2 to value-added products such as formic and oxalic acids.
The applicability of pure DES as electrolytes is hindered by high
viscosities. Mixtures of DES with organic solvents can be a promising
way of designing superior electrolytes by exploiting the advantages
of each solvent type. In this study, densities, viscosities, diffusivities,
and ionic conductivities of mixed solvents comprising DES (i.e., reline
and ethaline), methanol, and propylene carbonate were computed using
molecular simulations. To provide a quantitative assessment of the
affinity and mass transport of CO2 and oxalic and formic
acids in the mixed solvents, the solubilities and self-diffusivities
of these solutes were also computed. Our results show that the addition
of DES to the organic solvents enhances the solubilities of oxalic
and formic acids, while the solubility of CO2 in the ethaline-containing
mixtures are in the same order of magnitude with the respective pure
organic components. A monotonic increase in the densities and viscosities
of the mixed solvents is observed as the mole fraction of DES in the
mixture increases, with the exception of the density of ethaline-propylene
carbonate which shows the opposite behavior due to the high viscosity
of the pure organic component. The self-diffusivities of all species
in the mixtures significantly decrease as the mole fraction of DES
approaches unity. Similarly, the self-diffusivities of the dissolved
CO2 and the oxalic and formic acids also decrease by at
least 1 order of magnitude as the composition of the mixture shifts
from the pure organic component to pure DES. The computed ionic conductivities
of all mixed solvents show a maximum value for mole fractions of DES
in the range from 0.2 to 0.6 and decrease as more DES is added to
the mixtures. Since for most mixtures studied here no prior experimental
measurements exist, our findings can serve as a first data set based
on which further investigation of DES-containing electrolyte solutions
can be performed for the electrochemical reduction of CO2 to useful chemicals.
Collapse
Affiliation(s)
- Noura Dawass
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Jilles Langeveld
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Mahinder Ramdin
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Elena Pérez-Gallent
- Department of Sustainable Process and Energy Systems, TNO, Delft, Zuid-Holland 2628CA, The Netherlands
| | - Angel A Villanueva
- Department of Sustainable Process and Energy Systems, TNO, Delft, Zuid-Holland 2628CA, The Netherlands
| | - Erwin J M Giling
- Department of Sustainable Process and Energy Systems, TNO, Delft, Zuid-Holland 2628CA, The Netherlands
| | - Jort Langerak
- Research and Development Department, DMT Environmental Technology, Yndustrywei 3, 8501SN Joure, The Netherlands
| | - Leo J P van den Broeke
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & 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 & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| |
Collapse
|
13
|
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
|
14
|
Saravi SH, Panagiotopoulos AZ. Activity Coefficients and Solubilities of NaCl in Water-Methanol Solutions from Molecular Dynamics Simulations. J Phys Chem B 2022; 126:2891-2898. [PMID: 35411772 DOI: 10.1021/acs.jpcb.2c00813] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We obtain activity coefficients and solubilities of NaCl in water-methanol solutions at 298.15 K and 1 bar from molecular dynamics (MD) simulations with the Joung-Cheatham, SPC/E, and OPLS-AA force fields for NaCl, water, and methanol, respectively. The Lorentz-Berthelot combining rules were adopted for the unlike-pair interactions of Na+, Cl-, and the oxygen site in SPC/E water, and geometric combining rules were utilized for the remainder of the cross interactions. We found that the selection of appropriate combining rules is important in obtaining physically realistic solubilities. The solvent compositions studied range from pure water to pure methanol. Several salt concentrations were investigated at each solvent composition, from the lowest concentrations permitted by the system size used up to the experimental solubilities. We first calculated individual ion activity coefficients (IIACs) for Na+ and Cl- from the free energy change due to the gradual insertion of a single cation or anion into the solution, accompanied by a neutralizing background. We obtained the salt solubilities by comparing the chemical potentials in solution with solid NaCl chemical potentials calculated previously using the Einstein crystal method. Mean ionic activity coefficients obtained from the IIACs are in reasonable agreement with experimental data, with deviations increasing for solutions of higher methanol content. Predictions for the salt solubility are in surprisingly good agreement with experimental data, despite well-known challenges in the simultaneous calculation of activity coefficients and solubilities with classical MD simulations. The present study demonstrates that good predictions for these two important phase equilibrium properties can be obtained for mixed-solvent electrolyte solutions using existing nonpolarizable models and further suggests that the previously proposed single ion insertion technique can be extended to complex mixed-solvent solutions as well.
Collapse
Affiliation(s)
- Sina Hassanjani Saravi
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | | |
Collapse
|
15
|
Lamas CP, Vega C, Noya EG. Freezing point depression of salt aqueous solutions using the Madrid-2019 model. J Chem Phys 2022; 156:134503. [PMID: 35395902 DOI: 10.1063/5.0085051] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Salt aqueous solutions are relevant in many fields, ranging from biological systems to seawater. Thus, the availability of a force-field that is able to reproduce the thermodynamic and dynamic behavior of salt aqueous solutions would be of great interest. Unfortunately, this has been proven challenging, and most of the existing force-fields fail to reproduce much of their behavior. In particular, the diffusion of water or the salt solubility are often not well reproduced by most of the existing force-fields. Recently, the Madrid-2019 model was proposed, and it was shown that this force-field, which uses the TIP4P/2005 model for water and non-integer charges for the ions, provides a good description of a large number of properties, including the solution densities, viscosities, and the diffusion of water. In this work, we assess the performance of this force-field on the evaluation of the freezing point depression. Although the freezing point depression is a colligative property that at low salt concentrations depends solely on properties of pure water, a good model for the electrolytes is needed to accurately predict the freezing point depression at moderate and high salt concentrations. The coexistence line between ice and several salt aqueous solutions (NaCl, KCl, LiCl, MgCl2, and Li2SO4) up to the eutectic point is estimated from direct coexistence molecular dynamics simulations. Our results show that this force-field reproduces fairly well the experimentally measured freezing point depression with respect to pure water freezing for all the salts and at all the compositions considered.
Collapse
Affiliation(s)
- Cintia P Lamas
- Departamento de Química-Física I (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química-Física I (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eva G Noya
- Instituto de Química Física Rocasolano, CSIC, C/ Serrano 119, 28006 Madrid, Spain
| |
Collapse
|
16
|
Bianco V, Conde MM, Lamas CP, Noya EG, Sanz E. Phase diagram of the NaCl–water system from computer simulations. J Chem Phys 2022; 156:064505. [DOI: 10.1063/5.0083371] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- V. Bianco
- Departamento de Química Física (Unidad de I+D+i asociada al CSIC), Facultad de 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. P. Lamas
- Departamento de Química Física (Unidad de I+D+i asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
- Instituto de Química Física Rocasolano, Consejo Superior de Investigaciones Científicas, CSIC, Calle Serrano 119, 28006 Madrid, Spain
| | - E. G. Noya
- Instituto de Química Física Rocasolano, Consejo Superior de Investigaciones Científicas, CSIC, Calle Serrano 119, 28006 Madrid, Spain
| | - E. Sanz
- Departamento de Química Física (Unidad de I+D+i asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
17
|
Berkowitz ML. Molecular Simulations of Aqueous Electrolytes: Role of Explicit Inclusion of Charge Transfer into Force Fields. J Phys Chem B 2021; 125:13069-13076. [PMID: 34807628 DOI: 10.1021/acs.jpcb.1c08383] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We describe here simulations of aqueous salt solutions that are performed using an explicit charge transfer force field. The emphasis of the discussion is on the calculation of a dynamical property of the solutions: self-diffusion of water. While force fields that are based on pairwise additive potentials or on potentials with explicit inclusion of polarization or with scaled charges can provide at best a qualitative agreement with experiments, force fields with explicit inclusion of charge transfer can produce quantitative agreement with experiment for NaCl and KCl solutions. We argue that a force field with explicit charge transfer contains new physics absent in the previously used force fields described in recent reviews of molecular simulations of aqueous electrolytes.
Collapse
Affiliation(s)
- Max L Berkowitz
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
18
|
Saravi SH, Panagiotopoulos AZ. Activity coefficients of aqueous electrolytes from implicit-water molecular dynamics simulations. J Chem Phys 2021; 155:184501. [PMID: 34773944 DOI: 10.1063/5.0064963] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We obtain activity coefficients in NaCl and KCl solutions from implicit-water molecular dynamics simulations, at 298.15 K and 1 bar, using two distinct approaches. In the first approach, we consider ions in a continuum with constant relative permittivity (ɛr) equal to that of pure water; in the other approach, we take into account the concentration-dependence of ɛr, as obtained from explicit-water simulations. Individual ion activity coefficients (IIACs) are calculated using gradual insertion of single ions with uniform neutralizing backgrounds to ensure electroneutrality. Mean ionic activity coefficients (MIACs) obtained from the corresponding IIACs in simulations with constant ɛr show reasonable agreement with experimental data for both salts. Surprisingly, large systematic negative deviations are observed for both IIACs and MIACs in simulations with concentration-dependent ɛr. Our results suggest that the absence of hydration structure in implicit-water simulations cannot be compensated by correcting for the concentration-dependence of the relative permittivity ɛr. Moreover, even in simulations with constant ɛr for which the calculated MIACs are reasonable, the relative positioning of IIACs of anions and cations is incorrect for NaCl. We conclude that there are severe inherent limitations associated with implicit-water simulations in providing accurate activities of aqueous electrolytes, a finding with direct relevance to the development of electrolyte theories and to the use and interpretation of implicit-solvent simulations.
Collapse
Affiliation(s)
- Sina Hassanjani Saravi
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | | |
Collapse
|
19
|
Saravi SH, Panagiotopoulos AZ. Individual Ion Activity Coefficients in Aqueous Electrolytes from Explicit-Water Molecular Dynamics Simulations. J Phys Chem B 2021; 125:8511-8521. [PMID: 34319101 DOI: 10.1021/acs.jpcb.1c04019] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
We compute individual ion activity coefficients (IIACs) in aqueous NaCl, KCl, NaF, and KF solutions from explicit-water molecular dynamics simulations. Free energy changes are obtained from insertion of single ions-accompanied by uniform neutralizing backgrounds-into solution by gradually turning on first Lennard-Jones interactions, followed by Coulombic interactions using Ewald electrostatics. Simulations are performed at multiple system sizes, and all results are extrapolated to the thermodynamic limit, thus eliminating any possible artifacts from the neutralizing backgrounds. Because of controversies associated with measurements of IIACs from electrochemical cells with ion-selective electrodes, the reported experimental data are not widely accepted; thus there remains a knowledge gap with respect to the contributions of individual ions to solution nonidealities. Our results are in good qualitative agreement with these reported measurements, though significantly larger in magnitude. In particular, the relative positioning for the activity coefficients of anions and cations matches the experimental ordering for all four systems. This work establishes a robust thermodynamic framework, without a need to invoke extra hypotheses, that sheds light on the behavior of individual ions and their contributions to nonidealities of aqueous electrolyte solutions.
Collapse
Affiliation(s)
- Sina Hassanjani Saravi
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | | |
Collapse
|
20
|
Baryiames CP, Garrett P, Baiz CR. Bursting the bubble: A molecular understanding of surfactant-water interfaces. J Chem Phys 2021; 154:170901. [PMID: 34241044 DOI: 10.1063/5.0047377] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Surfactant science has historically emphasized bulk, thermodynamic measurements to understand the microemulsion properties of greatest industrial significance, such as interfacial tensions, phase behavior, and thermal stability. Recently, interest in the molecular properties of surfactants has grown among the physical chemistry community. This has led to the application of cutting-edge spectroscopic methods and advanced simulations to understand the specific interactions that give rise to the previously studied bulk characteristics. In this Perspective, we catalog key findings that describe the surfactant-oil and surfactant-water interfaces in molecular detail. We emphasize the role of ultrafast spectroscopic methods, including two-dimensional infrared spectroscopy and sum-frequency-generation spectroscopy, in conjunction with molecular dynamics simulations, and the role these techniques have played in advancing our understanding of interfacial properties in surfactant microemulsions.
Collapse
Affiliation(s)
- Christopher P Baryiames
- Department of Chemistry, University of Texas at Austin, 105 E 24th St. Stop A5300, Austin, Texas 78712-1224, USA
| | - Paul Garrett
- Department of Chemistry, University of Texas at Austin, 105 E 24th St. Stop A5300, Austin, Texas 78712-1224, USA
| | - Carlos R Baiz
- Department of Chemistry, University of Texas at Austin, 105 E 24th St. Stop A5300, Austin, Texas 78712-1224, USA
| |
Collapse
|
21
|
Abstract
Many important processes affecting the earth's climate are determined by the physical properties of seawater. In addition, desalination of seawater is a significant source of drinking water for the human population living in coastal areas. Since the physical properties of seawater governing these processes depend on the molecular interactions among its components, a deeper knowledge of seawater at the molecular level would contribute to a better understanding of these phenomena. However, in strong contrast with the situation in other areas such as biomolecules or materials science, molecular simulation studies reporting the physical properties of seawater are currently lacking. This is probably due to the usual perception of the seawater composition being too complex to approach. This point of view ignores the fact that physical properties of seawater are dependent on a single parameter representing the composition, namely the salinity. This is because the relative proportions of any two major constituents of seasalt are always the same. Another obstacle to performing molecular simulations of seawater could have been the unavailability of a satisfactory force field representing the interactions between water molecules and dissolved substances. However, this drawback has recently been overcome with the proposal of the Madrid-2019 force field. In this work we show for the first time that molecular simulation of seawater is feasible. We have performed molecular dynamics simulations of a system, the composition of which is close to the average composition of standard seawater and with the molecular interactions given by the Madrid-2019 force field. In this way we are able to provide quantitative or semiquantitative predictions for a number of relevant physical properties of seawater for temperatures and salinities from the oceanographic range to those relevant to desalination processes. The computed magnitudes include static (density), dynamical (viscosity and diffusion coefficients), structural (ionic hydration, ion-ion distribution functions), and interfacial (surface tension) properties.
Collapse
Affiliation(s)
- Ivan M Zeron
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Miguel A Gonzalez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Edoardo Errani
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jose L F Abascal
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
22
|
Rozsa VF, Galli G. Molecular Polarizabilities in Aqueous Systems from First-Principles. J Phys Chem B 2021; 125:2183-2192. [DOI: 10.1021/acs.jpcb.0c10732] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Viktor F. Rozsa
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Giulia Galli
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division, Argonne National Laboratory, Chicago, Illinois 60439, United States
| |
Collapse
|
23
|
Celebi AT, Vlugt TJH, Moultos OA. Thermal conductivity of aqueous solutions of reline, ethaline, and glyceline deep eutectic solvents; a molecular dynamics simulation study. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1876263] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Alper T. Celebi
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| |
Collapse
|
24
|
Malloum A, Fifen JJ, Conradie J. Determination of the absolute solvation free energy and enthalpy of the proton in solutions. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.114919] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
|
25
|
Baryiames CP, Ma E, Baiz CR. Ions Slow Water Dynamics at Nonionic Surfactant Interfaces. J Phys Chem B 2020; 124:11895-11900. [DOI: 10.1021/acs.jpcb.0c09086] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Christopher P. Baryiames
- Department of Chemistry, University of Texas at Austin, 105 East 24th Street, Stop A5300, Austin, Texas 78712-1224, United States
| | - Emily Ma
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Carlos R. Baiz
- Department of Chemistry, University of Texas at Austin, 105 East 24th Street, Stop A5300, Austin, Texas 78712-1224, United States
| |
Collapse
|
26
|
Noroozi J, Smith WR. Accurately Predicting CO2 Reactive Absorption Properties in Aqueous Alkanolamine Solutions by Molecular Simulation Requiring No Solvent Experimental Data. Ind Eng Chem Res 2020. [DOI: 10.1021/acs.iecr.0c03738] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Javad Noroozi
- Department of Chemical Engineering, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - William R. Smith
- Department of Chemical Engineering, University of Waterloo, Waterloo, ON N2L 3G1, Canada
- Department of Mathematics and Statistics, University of Guelph, Guelph, ON N1G 2W1, Canada
- Faculty of Science, University of Ontario Institute of Technology, Oshawa, ON L1H 7K4, Canada
| |
Collapse
|
27
|
Panagiotopoulos AZ. Simulations of activities, solubilities, transport properties, and nucleation rates for aqueous electrolyte solutions. J Chem Phys 2020; 153:010903. [DOI: 10.1063/5.0012102] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
|
28
|
Le Breton G, Joly L. Molecular modeling of aqueous electrolytes at interfaces: Effects of long-range dispersion forces and of ionic charge rescaling. J Chem Phys 2020; 152:241102. [PMID: 32610967 DOI: 10.1063/5.0011058] [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/13/2023] Open
Abstract
Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions-described by a truncated Lennard-Jones (LJ) potential-and electrostatic interactions-described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges [electronic continuum correction (ECC)], possibly complemented with rescaling of LJ parameters [ECC rescaled (ECCR)], have shown promising results in bulk, but their performance at interfaces has been less explored. Here, we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk accurately describes image charge repulsion at the liquid-vapor interface, and the rescaling of LJ parameters in ECCR models-aimed at capturing local ion-ion and ion-water interactions in bulk- describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.
Collapse
Affiliation(s)
- Guillaume Le Breton
- Département de Physique, École Normale Supérieure de Lyon, 46 Allée d'Italie, Lyon Cedex 07, France
| | - Laurent Joly
- Univ. Lyon, Univ. Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| |
Collapse
|
29
|
Shao Y, Shigenobu K, Watanabe M, Zhang C. Role of Viscosity in Deviations from the Nernst-Einstein Relation. J Phys Chem B 2020; 124:4774-4780. [PMID: 32412758 PMCID: PMC7497660 DOI: 10.1021/acs.jpcb.0c02544] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Revised: 05/05/2020] [Indexed: 11/30/2022]
Abstract
Deviations from the Nernst-Einstein relation are commonly attributed to ion-ion correlation and ion pairing. Despite the fact that these deviations can be quantified by either experimental measurements or molecular dynamics simulations, there is no rule of thumb to tell the extent of deviations. Here, we show that deviations from the Nernst-Einstein relation are proportional to the inverse viscosity by exploring the finite-size effect on transport properties under periodic boundary conditions. This conclusion is in accord with the established experimental results of ionic liquids.
Collapse
Affiliation(s)
- Yunqi Shao
- Department
of Chemistry-Ångström Laboratory, Uppsala University, Lägerhyddsvägen 1, P.O. Box 538, 75121 Uppsala, Sweden
| | - Keisuke Shigenobu
- Department
of Chemistry and Biotechnology, Yokohama
National University, Yokohama 240-8501, Japan
| | - Masayoshi Watanabe
- Department
of Chemistry and Biotechnology, Yokohama
National University, Yokohama 240-8501, Japan
| | - Chao Zhang
- Department
of Chemistry-Ångström Laboratory, Uppsala University, Lägerhyddsvägen 1, P.O. Box 538, 75121 Uppsala, Sweden
| |
Collapse
|
30
|
Dočkal J, Lísal M, Moučka F. Molecular Force Field Development for Aqueous Electrolytes: 2. Polarizable Models Incorporating Crystalline Chemical Potential and Their Accurate Simulations of Halite, Hydrohalite, Aqueous Solutions of NaCl, and Solubility. J Chem Theory Comput 2020; 16:3677-3688. [DOI: 10.1021/acs.jctc.0c00161] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jan Dočkal
- Department of Physics, Faculty of Science, J. E. Purkyně University, 400 96 Ústí n. Lab., Czech Republic
| | - Martin Lísal
- Department of Physics, Faculty of Science, J. E. Purkyně University, 400 96 Ústí n. Lab., Czech Republic
- Department of Molecular and Mesoscopic Modelling, Institute of Chemical Process Fundamentals of the CAS, v. v. i., 165 02 Prague 6-Suchdol, Czech Republic
| | - Filip Moučka
- Department of Physics, Faculty of Science, J. E. Purkyně University, 400 96 Ústí n. Lab., Czech Republic
- Department of Molecular and Mesoscopic Modelling, Institute of Chemical Process Fundamentals of the CAS, v. v. i., 165 02 Prague 6-Suchdol, Czech Republic
| |
Collapse
|
31
|
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
|
32
|
Xin X, Niu X, Liu W, Wang D. Hybrid Solvation Model with First Solvation Shell for Calculation of Solvation Free Energy. Chemphyschem 2020; 21:762-769. [PMID: 32154979 DOI: 10.1002/cphc.202000039] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2020] [Revised: 02/15/2020] [Indexed: 02/03/2023]
Abstract
We present a hybrid solvation model with first solvation shell to calculate solvation free energies. This hybrid model combines the quantum mechanics and molecular mechanics methods with the analytical expression based on the Born solvation model to calculate solvation free energies. Based on calculated free energies of solvation and reaction profiles in gas phase, we set up a unified scheme to predict reaction profiles in solution. The predicted solvation free energies and reaction barriers are compared with experimental results for twenty bimolecular nucleophilic substitution reactions. These comparisons show that our hybrid solvation model can predict reliable solvation free energies and reaction barriers for chemical reactions of small molecules in aqueous solution.
Collapse
Affiliation(s)
- Xin Xin
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| | - Xiao Niu
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| | - Wanqi Liu
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong, 250014, China
| |
Collapse
|
33
|
Lbadaoui-Darvas M, Takahama S. Water Activity from Equilibrium Molecular Dynamics Simulations and Kirkwood-Buff Theory. J Phys Chem B 2019; 123:10757-10768. [DOI: 10.1021/acs.jpcb.9b06735] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Mária Lbadaoui-Darvas
- ENAC/IIE Swiss Federal Institute of Technology Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Satoshi Takahama
- ENAC/IIE Swiss Federal Institute of Technology Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
34
|
Seal S, Doblhoff-Dier K, Meyer J. Dielectric Decrement for Aqueous NaCl Solutions: Effect of Ionic Charge Scaling in Nonpolarizable Water Force Fields. J Phys Chem B 2019; 123:9912-9921. [PMID: 31647235 PMCID: PMC6875873 DOI: 10.1021/acs.jpcb.9b07916] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We investigate the dielectric constant and the dielectric decrement of aqueous NaCl
solutions by means of molecular dynamic simulations. We thereby compare the performance
of four different force fields and focus on disentangling the origin of the dielectric
decrement and the influence of scaled ionic charges, as often used in nonpolarizable
force fields to account for the missing dynamic polarizability in the shielding of
electrostatic ion interactions. Three of the force fields showed excessive contact ion
pair formation, which correlates with a reduced dielectric decrement. In spite of the
fact that the scaling of charges only weakly influenced the average polarization of
water molecules around an ion, the rescaling of ionic charges did influence the
dielectric decrement, and a close-to-linear relation of the slope of the dielectric
constant as a function of concentration with the ionic charge was found.
Collapse
Affiliation(s)
- Sayan Seal
- Gorlaeus Laboratories, Leiden Institute of Chemistry , Leiden University , PO Box 9502, 2300 RA Leiden , The Netherlands
| | - Katharina Doblhoff-Dier
- Gorlaeus Laboratories, Leiden Institute of Chemistry , Leiden University , PO Box 9502, 2300 RA Leiden , The Netherlands
| | - Jörg Meyer
- Gorlaeus Laboratories, Leiden Institute of Chemistry , Leiden University , PO Box 9502, 2300 RA Leiden , The Netherlands
| |
Collapse
|
35
|
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: 131] [Impact Index Per Article: 26.2] [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
|
36
|
Dočkal J, Svoboda M, Lísal M, Moučka F. A general hydrogen bonding definition based on three-dimensional spatial distribution functions and its extension to quantitative structural analysis of solutions and general intermolecular bonds. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.02.036] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
37
|
Noroozi J, Smith WR. An Efficient Molecular Simulation Methodology for Chemical Reaction Equilibria in Electrolyte Solutions: Application to CO2 Reactive Absorption. J Phys Chem A 2019; 123:4074-4086. [DOI: 10.1021/acs.jpca.9b00302] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Javad Noroozi
- Department of Chemical Engineering, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - William R. Smith
- Department of Chemical Engineering, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
- Department of Mathematics and Statistics, University of Guelph, Guelph, Ontario N1G 2W1, Canada
- Faculty of Science, University of Ontario Institute of Technology, Oshawa, Ontario L1H 7K4, Canada
| |
Collapse
|
38
|
Jiang H, Debenedetti PG, Panagiotopoulos AZ. Nucleation in aqueous NaCl solutions shifts from 1-step to 2-step mechanism on crossing the spinodal. J Chem Phys 2019; 150:124502. [DOI: 10.1063/1.5084248] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Hao Jiang
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Pablo G. Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | | |
Collapse
|
39
|
Bellucci MA, Gobbo G, Wijethunga TK, Ciccotti G, Trout BL. Solubility of paracetamol in ethanol by molecular dynamics using the extended Einstein crystal method and experiments. J Chem Phys 2019; 150:094107. [PMID: 30849885 DOI: 10.1063/1.5086706] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Li and co-workers [Li et al., J. Chem. Phys. 146, 214110 (2017)] have recently proposed a methodology to compute the solubility of molecular compounds from first principles, using molecular dynamics simulations. We revise and further explore their methodology that was originally applied to naphthalene in water at low concentration. In particular, we compute the solubility of paracetamol in an ethanol solution at ambient conditions. For the simulations, we used a force field that we previously reparameterized to reproduce certain thermodynamic properties of paracetamol but not explicitly its solubility in ethanol. In addition, we have determined the experimental solubility by performing turbidity measurements using a Crystal16 over a range of temperatures. Our work serves a dual purpose: (i) methodologically, we clarify how to compute, with a relatively straightforward procedure, the solubility of molecular compounds and (ii) applying this procedure, we show that the solubility predicted by our force field (0.085 ± 0.014 in mole ratio) is in good agreement with the experimental value obtained from our experiments and those reported in the literature (average 0.0585 ± 0.004), considering typical deviations for predictions from first principle methods. The good agreement between the experimental and the calculated solubility also suggests that the method used to reparameterize the force field can be used as a general strategy to optimize force fields for simulations in solution.
Collapse
Affiliation(s)
- Michael A Bellucci
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Gianpaolo Gobbo
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Tharanga K Wijethunga
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | | | - Bernhardt L Trout
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
40
|
Sayer T, Sprik M, Zhang C. Finite electric displacement simulations of polar ionic solid-electrolyte interfaces: Application to NaCl(111)/aqueous NaCl solution. J Chem Phys 2019; 150:041716. [DOI: 10.1063/1.5054843] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Thomas Sayer
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Michiel Sprik
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Chao Zhang
- Department of Chemistry-Ångström Laboratory, Uppsala University, Lägerhyddsvägen 1, P.O. Box 538, 75121 Uppsala, Sweden
| |
Collapse
|
41
|
Nezbeda I. Vapour–liquid equilibria from molecular simulations: some issues affecting reliability and reproducibility. Mol Phys 2019. [DOI: 10.1080/00268976.2018.1562125] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- I. Nezbeda
- Faculty of Science, J.E. Purkinje University, Ústí nad Labem, Czech Republic
| |
Collapse
|
42
|
Espinosa JR, Wand CR, Vega C, Sanz E, Frenkel D. Calculation of the water-octanol partition coefficient of cholesterol for SPC, TIP3P, and TIP4P water. J Chem Phys 2018; 149:224501. [PMID: 30553262 DOI: 10.1063/1.5054056] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a numerical study of the relative solubility of cholesterol in octanol and water. Our calculations allow us to compare the accuracy of the computed values of the excess chemical potential of cholesterol for several widely used water models (SPC, TIP3P, and TIP4P). We compute the excess solvation free energies by means of a cavity-based method [L. Li et al., J. Chem. Phys. 146(21), 214110 (2017)] which allows for the calculation of the excess chemical potential of a large molecule in a dense solvent phase. For the calculation of the relative solubility ("partition coefficient," log10 P o / w ) of cholesterol between octanol and water, we use the OPLS/AA force field in combination with the SPC, TIP3P, and TIP4P water models. For all water models studied, our results reproduce the experimental observation that cholesterol is less soluble in water than in octanol. While the experimental value for the partition coefficient is log10 P o / w = 3.7, SPC, TIP3P, and TIP4P give us a value of log10 P o / w = 4.5, 4.6, and 2.9, respectively. Therefore, although the results for the studied water models in combination with the OPLS/AA force field are acceptable, further work to improve the accuracy of current force fields is needed.
Collapse
Affiliation(s)
- Jorge R Espinosa
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Charlie R Wand
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Carlos Vega
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Daan Frenkel
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
43
|
Svoboda M, Moučka F, Lísal M. Saturated aqueous NaCl solution and pure water in Na-montmorillonite clay at thermodynamic conditions of hydraulic fracturing: Thermodynamics, structure and diffusion from molecular simulations. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.08.144] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
44
|
Prasetyo N, Hünenberger PH, Hofer TS. Single-Ion Thermodynamics from First Principles: Calculation of the Absolute Hydration Free Energy and Single-Electrode Potential of Aqueous Li + Using ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations. J Chem Theory Comput 2018; 14:6443-6459. [PMID: 30284829 DOI: 10.1021/acs.jctc.8b00729] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
A recently proposed thermodynamic integration (TI) approach formulated in the framework of quantum mechanical/molecular mechanical molecular dynamics (QM/MM MD) simulations is applied to study the structure, dynamics, and absolute intrinsic hydration free energy Δs GM+,wat◦ of the Li+ ion at a correlated ab initio level of theory. Based on the results, standard values (298.15 K, ideal gas at 1 bar, ideal solute at 1 molal) for the absolute intrinsic hydration free energy [Formula: see text] of the proton, the surface electric potential jump χwat◦ upon entering bulk water, and the absolute single-electrode potential [Formula: see text] of the reference hydrogen electrode are calculated to be -1099.9 ± 4.2 kJ·mol-1, 0.13 ± 0.08 V, and 4.28 ± 0.04 V, respectively, in excellent agreement with the standard values recommended by Hünenberger and Reif on the basis of an extensive evaluation of the available experimental data (-1100 ± 5 kJ·mol-1, 0.13 ± 0.10 V, and 4.28 ± 0.13 V). The simulation results for Li+ are also compared to those for Na+ and K+ from a previous study in terms of relative hydration free energies ΔΔs GM+,wat◦ and relative electrode potentials [Formula: see text]. The calculated values are found to agree extremely well with the experimental differences in standard conventional hydration free energies ΔΔs GM+,wat• and redox potentials [Formula: see text]. The level of agreement between simulation and experiment, which is quantitative within error bars, underlines the substantial accuracy improvement achieved by applying a highly demanding QM/MM approach at the resolution-of-identity second-order Møller-Plesset perturbation (RIMP2) level over calculations relying on purely molecular mechanical or density functional theory (DFT) descriptions. A detailed analysis of the structural and dynamical properties of the Li+ hydrate indicates that a correct description of the solvation structure and dynamics is achieved as well at this level of theory. Consideration of the QM/MM potential-energy components also shows that the partitioning into QM and MM zones does not induce any significant energetic artifact for the system considered.
Collapse
Affiliation(s)
- Niko Prasetyo
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry , University of Innsbruck , Innrain 80-82 , A-6020 Innsbruck , Austria.,Austria-Indonesia Centre (AIC) for Computational Chemistry , Universitas Gadjah Mada , Sekip Utara , Yogyakarta 55281 , Indonesia.,Department of Chemistry, Faculty of Mathematics and Natural Sciences , Universitas Gadjah Mada , Sekip Utara , Yogyakarta 55281 , Indonesia
| | - Philippe H Hünenberger
- Laboratorium für Physikalische Chemie , ETH Zürich, ETH-Hönggerberg , HCI Building , CH-8093 Zürich , Switzerland
| | - Thomas S Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry , University of Innsbruck , Innrain 80-82 , A-6020 Innsbruck , Austria
| |
Collapse
|
45
|
Škvára J, Nezbeda I. Surface of aqueous solutions of alkali halides: layer by layer analysis. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1540871] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Jiří Škvára
- Faculty of Science, J. E. Purkinje University, Ústí nad Labem, Czech Republic
| | - Ivo Nezbeda
- Faculty of Science, J. E. Purkinje University, Ústí nad Labem, Czech Republic
| |
Collapse
|
46
|
Patel LA, Kindt JT. Simulations of NaCl Aggregation from Solution: Solvent Determines Topography of Free Energy Landscape. J Comput Chem 2018; 40:135-147. [DOI: 10.1002/jcc.25554] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2018] [Revised: 07/10/2018] [Accepted: 07/11/2018] [Indexed: 11/10/2022]
Affiliation(s)
- Lara A. Patel
- Department of Chemistry; Emory University; 1515 Dickey Drive, Atlanta Georgia 30322
| | - James T. Kindt
- Department of Chemistry; Emory University; 1515 Dickey Drive, Atlanta Georgia 30322
| |
Collapse
|
47
|
Li L, Totton T, Frenkel D. Computational methodology for solubility prediction: Application to sparingly soluble organic/inorganic materials. J Chem Phys 2018; 149:054102. [PMID: 30089373 DOI: 10.1063/1.5040366] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The solubility of a crystalline material can be estimated from the absolute free energy of the solid and the excess solvation free energy. In the earlier work, we presented a general-purpose molecular-dynamics-based methodology enabling solubility predictions of crystalline compounds, yielding accurate estimates of the aqueous solubilities of naphthalene at various pressures and temperatures. In the present work, we investigate a number of prototypical complex materials, including phenanthrene, calcite, and aragonite over a range of temperatures and pressures. Our results provide stronger evidence for the power of the methodology for universal solubility predictions.
Collapse
Affiliation(s)
- Lunna Li
- Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| | - Tim Totton
- BP Exploration Operating Co. Ltd., Sunbury-on-Thames TW16 7LN, United Kingdom
| | - Daan Frenkel
- Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| |
Collapse
|
48
|
Wand CR, Totton TS, Frenkel D. Addressing hysteresis and slow equilibration issues in cavity-based calculation of chemical potentials. J Chem Phys 2018; 149:014105. [PMID: 29981554 DOI: 10.1063/1.5036963] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we explore the strengths and weaknesses of a cavity-based method to calculate the excess chemical potential of a large molecular solute in a dense liquid solvent. Use of the cavity alleviates some technical problems associated with the appearance of (integrable) divergences in the integrand during alchemical particle growth. The excess chemical potential calculated using the cavity-based method should be independent of the cavity attributes. However, the performance of the method (equilibration time and the robustness) does depend on the cavity attributes. To illustrate the importance of a suitable choice of the cavity attributes, we calculate the partition coefficient of pyrene in toluene and heptane using a coarse-grained model. We find that a poor choice for the functional form of the cavity may lead to hysteresis between growth and shrinkage of the cavity. Somewhat unexpectedly, we find that, by allowing the cavity to move as a pseudo-particle within the simulation box, the decay time of fluctuations in the integrand of the thermodynamic integration can be reduced by an order of magnitude, thereby increasing the statistical accuracy of the calculation.
Collapse
Affiliation(s)
- C R Wand
- Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| | - T S Totton
- BP Exploration Operating Co. Ltd., Sunbury-on-Thames TW16 7LN, United Kingdom
| | - D Frenkel
- Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| |
Collapse
|
49
|
Soria GD, Espinosa JR, Ramirez J, Valeriani C, Vega C, Sanz E. A simulation study of homogeneous ice nucleation in supercooled salty water. J Chem Phys 2018; 148:222811. [DOI: 10.1063/1.5008889] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Guiomar D. Soria
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jorge R. Espinosa
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jorge Ramirez
- Departamento de Ingenieria Quimica Industrial y Medio Ambiente, Escuela Tecnica Superior de Ingenieros Industriales, Universidad Politecnica de Madrid, 28006 Madrid, Spain
| | - Chantal Valeriani
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
- Departamento de Fisica Aplicada I, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
50
|
Svoboda M, Lísal M. Concentrated aqueous sodium chloride solution in clays at thermodynamic conditions of hydraulic fracturing: Insight from molecular dynamics simulations. J Chem Phys 2018; 148:222806. [DOI: 10.1063/1.5017166] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Martin Svoboda
- Department of Molecular and Mesoscopic Modelling, Institute of Chemical Process Fundamentals of the CAS, v. v. i., Prague, Czech Republic
- Department of Physics, Faculty of Science, J. E. Purkinje University, Ústí nad Labem, Czech Republic
| | - Martin Lísal
- Department of Molecular and Mesoscopic Modelling, Institute of Chemical Process Fundamentals of the CAS, v. v. i., Prague, Czech Republic
- Department of Physics, Faculty of Science, J. E. Purkinje University, Ústí nad Labem, Czech Republic
| |
Collapse
|