### Abstract

The effects of Coulomb potential truncation schemes used in computer simulations of ionic and polar fluids are examined by use of integral equation techniques. A renormalized HNC type equation capable of describing both ionic and polar molecular fluids with truncated interactions is derived and applied to several model systems of interest. Good agreement is found between the integral equation results and Monte Carlo simulations of the same potential for dilute solutions of ions in a dielectric continuum. Very large effects on the distribution functions result from truncation of the electrostatic interaction in dilute systems. Even in comparatively dense systems, unrealistic pair correlations near the cutoff distance result from some of the proposed truncation schemes. The effect of Coulomb potential truncation for a molecular model of pure water is also studied. Significant errors appear in the second neighbor region for commonly used truncation schemes; a simple switching function that zeros the potential and its first derivative yields results closest to the Coulomb potential without truncation.

Original language | English (US) |
---|---|

Pages (from-to) | 5897-5908 |

Number of pages | 12 |

Journal | The Journal of Chemical Physics |

Volume | 83 |

Issue number | 11 |

State | Published - 1985 |

Externally published | Yes |

### Fingerprint

### ASJC Scopus subject areas

- Atomic and Molecular Physics, and Optics

### Cite this

*The Journal of Chemical Physics*,

*83*(11), 5897-5908.

**Structural and energetic effects of truncating long ranged interactions in ionic and polar fluids.** / Brooks, Charles L.; Pettitt, Bernard; Karplus, Martin.

Research output: Contribution to journal › Article

*The Journal of Chemical Physics*, vol. 83, no. 11, pp. 5897-5908.

}

TY - JOUR

T1 - Structural and energetic effects of truncating long ranged interactions in ionic and polar fluids

AU - Brooks, Charles L.

AU - Pettitt, Bernard

AU - Karplus, Martin

PY - 1985

Y1 - 1985

N2 - The effects of Coulomb potential truncation schemes used in computer simulations of ionic and polar fluids are examined by use of integral equation techniques. A renormalized HNC type equation capable of describing both ionic and polar molecular fluids with truncated interactions is derived and applied to several model systems of interest. Good agreement is found between the integral equation results and Monte Carlo simulations of the same potential for dilute solutions of ions in a dielectric continuum. Very large effects on the distribution functions result from truncation of the electrostatic interaction in dilute systems. Even in comparatively dense systems, unrealistic pair correlations near the cutoff distance result from some of the proposed truncation schemes. The effect of Coulomb potential truncation for a molecular model of pure water is also studied. Significant errors appear in the second neighbor region for commonly used truncation schemes; a simple switching function that zeros the potential and its first derivative yields results closest to the Coulomb potential without truncation.

AB - The effects of Coulomb potential truncation schemes used in computer simulations of ionic and polar fluids are examined by use of integral equation techniques. A renormalized HNC type equation capable of describing both ionic and polar molecular fluids with truncated interactions is derived and applied to several model systems of interest. Good agreement is found between the integral equation results and Monte Carlo simulations of the same potential for dilute solutions of ions in a dielectric continuum. Very large effects on the distribution functions result from truncation of the electrostatic interaction in dilute systems. Even in comparatively dense systems, unrealistic pair correlations near the cutoff distance result from some of the proposed truncation schemes. The effect of Coulomb potential truncation for a molecular model of pure water is also studied. Significant errors appear in the second neighbor region for commonly used truncation schemes; a simple switching function that zeros the potential and its first derivative yields results closest to the Coulomb potential without truncation.

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

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

M3 - Article

VL - 83

SP - 5897

EP - 5908

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 11

ER -