TY - JOUR

T1 - Thermodynamics of ion solvation in dipolar solvent using Monte Carlo mean reaction field simulation

AU - Bandura, A. V.

AU - Lvov, S. N.

AU - Macdonald, D. D.

PY - 1998/4/21

Y1 - 1998/4/21

N2 - The mean reaction field approximation has been developed for Monte Carlo simulation of ion solvation thermodynamics and structure in ion-dipolar systems. Calculations have been carried out for systems having one hard-sphere ion and N - 1 hardsphere dipoles of equal diameter, 256 ≤ N ≤ 1372, where N is the number of particles. The reduced charge on the ion, q*, was varied from 0 to 14.15, and two values of the reduced dipole moment, μ*, (1.14 and 1.50) were used for the simulations. The mean reaction field corrections were introduced to take into account the changes in long-range dipole-dipole interactions due to the presence of the ion. The corrections were performed for molecules both inside and outside the truncation sphere, which has an ion at the center, resulting in a small dependence of solvation energy on simulation cell size. Expressions have been developed within the framework of thermodynamic integration and perturbation techniques for calculating the free-energy change on ion charging. The applicability of this approach for simulating the solvation internal energy has been investigated and compared with other methods. The influence of system size and cut-off radius on the thermodynamic solvation parameters and the radial distribution functions were studied. It has been shown that for N ≥ 500, variation of the simulated system size and cut-off radius does not drastically affect the calculated thermodynamic quantities and radial distribution functions when the mean reaction field corrections were applied. The simulated thermodynamic solvation properties were compared with results found in the literature for various molecular statistical approaches [mean spherical approximation (MSA), linearized hypemetted-chain (LHNC), etc.].

AB - The mean reaction field approximation has been developed for Monte Carlo simulation of ion solvation thermodynamics and structure in ion-dipolar systems. Calculations have been carried out for systems having one hard-sphere ion and N - 1 hardsphere dipoles of equal diameter, 256 ≤ N ≤ 1372, where N is the number of particles. The reduced charge on the ion, q*, was varied from 0 to 14.15, and two values of the reduced dipole moment, μ*, (1.14 and 1.50) were used for the simulations. The mean reaction field corrections were introduced to take into account the changes in long-range dipole-dipole interactions due to the presence of the ion. The corrections were performed for molecules both inside and outside the truncation sphere, which has an ion at the center, resulting in a small dependence of solvation energy on simulation cell size. Expressions have been developed within the framework of thermodynamic integration and perturbation techniques for calculating the free-energy change on ion charging. The applicability of this approach for simulating the solvation internal energy has been investigated and compared with other methods. The influence of system size and cut-off radius on the thermodynamic solvation parameters and the radial distribution functions were studied. It has been shown that for N ≥ 500, variation of the simulated system size and cut-off radius does not drastically affect the calculated thermodynamic quantities and radial distribution functions when the mean reaction field corrections were applied. The simulated thermodynamic solvation properties were compared with results found in the literature for various molecular statistical approaches [mean spherical approximation (MSA), linearized hypemetted-chain (LHNC), etc.].

UR - http://www.scopus.com/inward/record.url?scp=33749639521&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33749639521&partnerID=8YFLogxK

U2 - 10.1039/a706091b

DO - 10.1039/a706091b

M3 - Article

AN - SCOPUS:33749639521

SN - 0956-5000

VL - 94

SP - 1063

EP - 1072

JO - Journal of the Chemical Society - Faraday Transactions

JF - Journal of the Chemical Society - Faraday Transactions

IS - 8

ER -