A new approach for efficient simulation of Coulomb interactions in ionic fluids
- PMID: 18376910
- DOI: 10.1063/1.2894478
A new approach for efficient simulation of Coulomb interactions in ionic fluids
Abstract
We propose a simplified version of local molecular field (LMF) theory to treat Coulomb interactions in simulations of ionic fluids. LMF theory relies on splitting the Coulomb potential into a short-ranged part that combines with other short-ranged core interactions and is simulated explicitly. The averaged effects of the remaining long-ranged part are taken into account through a self-consistently determined effective external field. The theory contains an adjustable length parameter sigma that specifies the cutoff distance for the short-ranged interaction. This can be chosen to minimize the errors resulting from the mean-field treatment of the complementary long-ranged part. Here we suggest that in many cases an accurate approximation to the effective field can be obtained directly from the equilibrium charge density given by the Debye theory of screening, thus eliminating the need for a self-consistent treatment. In the limit sigma-->0, this assumption reduces to the classical Debye approximation. We examine the numerical performance of this approximation for a simple model of a symmetric ionic mixture. Our results for thermodynamic and structural properties of uniform ionic mixtures agree well with similar results of Ewald simulations of the full ionic system. In addition, we have used the simplified theory in a grand-canonical simulation of a nonuniform ionic mixture where an ion has been fixed at the origin. Simulations using short-ranged truncations of the Coulomb interactions alone do not satisfy the exact condition of complete screening of the fixed ion, but this condition is recovered when the effective field is taken into account. We argue that this simplified approach can also be used in the simulations of more complex nonuniform systems.
Similar articles
-
Accurate thermodynamics for short-ranged truncations of Coulomb interactions in site-site molecular models.J Chem Phys. 2009 Dec 28;131(24):244108. doi: 10.1063/1.3276729. J Chem Phys. 2009. PMID: 20059055
-
Structure of nonuniform hard sphere fluids from shifted linear truncations of functional expansions.J Phys Chem B. 2005 Apr 14;109(14):6892-901. doi: 10.1021/jp0446884. J Phys Chem B. 2005. PMID: 16851776
-
Electrostatic correlations in inhomogeneous charged fluids beyond loop expansion.J Chem Phys. 2012 Sep 14;137(10):104902. doi: 10.1063/1.4750044. J Chem Phys. 2012. PMID: 22979885
-
Double-layer in ionic liquids: paradigm change?J Phys Chem B. 2007 May 24;111(20):5545-57. doi: 10.1021/jp067857o. Epub 2007 May 1. J Phys Chem B. 2007. PMID: 17469864
-
Local Molecular Field Theory for Coulomb Interactions in Aqueous Solutions.J Phys Chem B. 2023 Feb 2;127(4):809-821. doi: 10.1021/acs.jpcb.2c06988. Epub 2023 Jan 20. J Phys Chem B. 2023. PMID: 36669139 Review.
Cited by
-
Isotropic periodic sum for multipole interactions and a vector relation for calculation of the Cartesian multipole tensor.J Chem Phys. 2016 Oct 28;145(16):164110. doi: 10.1063/1.4966019. J Chem Phys. 2016. PMID: 27802614 Free PMC article.
-
Isotropic periodic sum of electrostatic interactions for polar systems.J Chem Phys. 2009 Jul 14;131(2):024107. doi: 10.1063/1.3160730. J Chem Phys. 2009. PMID: 19603970 Free PMC article.
-
Interplay of local hydrogen-bonding and long-ranged dipolar forces in simulations of confined water.Proc Natl Acad Sci U S A. 2008 Dec 9;105(49):19136-41. doi: 10.1073/pnas.0807623105. Epub 2008 Dec 8. Proc Natl Acad Sci U S A. 2008. PMID: 19064931 Free PMC article.
-
Long-ranged contributions to solvation free energies from theory and short-ranged models.Proc Natl Acad Sci U S A. 2016 Mar 15;113(11):2819-26. doi: 10.1073/pnas.1521570113. Epub 2016 Feb 29. Proc Natl Acad Sci U S A. 2016. PMID: 26929375 Free PMC article.
LinkOut - more resources
Full Text Sources