A generalized reaction field method for molecular dynamics simulations

Tironi, Ilario G.; Sperb, René; Smith, Paul E.; van Gunsteren, Wilfred F.
April 1995
Journal of Chemical Physics;4/1/1995, Vol. 102 Issue 13, p5451
Academic Journal
Molecular dynamics simulations of ionic systems require the inclusion of long-range electrostatic forces. We propose an expression for the long-range electrostatic forces based on an analytical solution of the Poisson–Boltzmann equation outside a spherical cutoff, which can easily be implemented in molecular simulation programs. An analytical solution of the linearized Poisson–Boltzmann (PB) equation valid in a spherical region is obtained. From this general solution special expressions are derived for evaluating the electrostatic potential and its derivative at the origin of the sphere. These expressions have been implemented for molecular dynamics (MD) simulations, such that the surface of the cutoff sphere around a charged particle is identified with the spherical boundary of the Poisson–Boltzmann problem. The analytical solution of the Poisson–Boltzmann equation is valid for the cutoff sphere and can be used for calculating the reaction field forces on the central charge, assuming a uniform continuum of given ionic strength beyond the cutoff. MD simulations are performed for a periodic system consisting of 2127 SPC water molecules with 40 NaCl ions (1 molar). We compare the structural and dynamical results obtained from MD simulations in which the long range electrostatic interactions are treated differently; using a cutoff radius, using a cutoff radius and a Poisson–Boltzmann generalized reaction field force, and using the Ewald summation. Application of the Poisson–Boltzmann generalized reaction field gives a dramatic improvement of the structure of the solution compared to a simple cutoff treatment, at no extra computational cost. © 1995 American Institute of Physics.


Related Articles

  • An ab initio molecular dynamics study of the S[sub N]2 reaction... Raugei, Simone; Cardini, Gianni // Journal of Chemical Physics;12/22/1999, Vol. 111 Issue 24, p10887 

    Describes an ab initio molecular dynamics study of the S[sub N]2 reaction Cl[sup -]CH[sub 3]Br...CH[sub 3]Cl+Br[sup -]. Potential energy and free energy profiles determined along the reaction coordinate; Analysis of the structural parameters along the reaction pathway; Results of impact studies.

  • Forward-backward initial value representation for the calculation of thermal rate constants for reactions in complex molecular systems. Wang, Haobin; Haobin Wang; Thoss, Michael; Miller, William H. // Journal of Chemical Physics;1/1/2000, Vol. 112 Issue 1 

    The semiclassical (SC) initial value representation (IVR) provides a potentially practical way for including quantum effects into classical molecular dynamics simulations. The forward-backward (FB) version of the IVR provides an especially attractive way for calculating time correlation...

  • Dynamic response of an irreversible catalytic reaction to periodic variation of the reactant's pressure. López, A. C.; Albano, E. V. // Journal of Chemical Physics;2/22/2000, Vol. 112 Issue 8 

    The dynamic response of the Ziff, Gulari, and Barshad model for the catalytic oxidation of CO, to a periodic variation of the external pressure is studied close to the CO-poisoning transition. Due to the applied pressure, a considerable enhancement in the catalytic activity is observed when the...

  • A new mechanism for the enhancement of activated bimolecular reactions by rotational excitation. Miklavc, Adolf; Perdih, Marko; Smith, Ian W. // Journal of Chemical Physics;5/22/2000, Vol. 112 Issue 20 

    The kinematic mass model (KMM), which has been developed to examine the dynamics of activated bimolecular reactions, has here been adapted to examine how orientational effects associated with reagent rotation influence the rotational state dependence of reaction cross-section. It is shown that,...

  • Probing the effect of the H[sub 2] rotational state in O([sup 1]D)+H[sub 2]→OH+H: Theoretical dynamics including nonadiabatic effects and a crossed molecular beam study. Gray, Stephen K.; Flikkema, Edwin; Balint-Kurti, Gabriel G.; Subbotin, Andrei; Schatz, George C.; ten Brinke, Gerrit; Lin, Jim J.; Liu, Xianghong; Harich, Steve; Yang, Xueming // Journal of Chemical Physics;11/1/2000, Vol. 113 Issue 17 

    Theoretical estimates of reactive cross sections for O([sup 1]D)+H[sub 2](X,v=0,j)→OH(X)+H([sup 2]S), with H[sub 2] rotational quantum numbers j=0 and 1, are obtained for a range of collision energies, E[sub col]. Crossed molecular beam measurements are also used to infer the ratio, r[sub...

  • Theoretical studies of intersystem crossing effects in the O+H[sub 2] reaction. Hoffmann, Mark R.; Schatz, George C. // Journal of Chemical Physics;12/1/2000, Vol. 113 Issue 21 

    We present a general procedure for studying intersystem crossing effects in bimolecular chemical reactions, along with an application of this to the O+H[sub 2] reaction. In this procedure, we use previously derived singlet and triplet potential energy surfaces that were based on high quality...

  • Reaction dynamics on a thermally fluctuating potential. Berezhkovskii, Alexander M.; Szabo, Attila; Weiss, George H.; Huan-Xiang Zhou // Journal of Chemical Physics;12/8/1999, Vol. 111 Issue 22, p9952 

    Analyzes the reaction dynamics on a thermally fluctuating potential. Jump rates for transitions between two states; Analysis of fast-fluctuation limit in terms of the potential of mean force; Turnover behavior.

  • A Study of Deflagration to Detonation Transition in a Model A-B System Using Molecular Dynamics. Fellows, Justin; Haskins, Peter J.; Cook, Malcolm D. // AIP Conference Proceedings;2002, Vol. 620 Issue 1, p391 

    Molecular Dynamics employing a model 2 dimensional molecular AB system has been used to study the transition from burning through to detonation. In these studies a REBO (Reactive Empirical Bond Order) potential was used to describe the interaction between the atoms. A small block of AB molecules...

  • Semiclassical initial value representation for the Boltzmann operator in thermal rate constants. Zhao, Yi; Miller, William H. // Journal of Chemical Physics;12/1/2002, Vol. 117 Issue 21, p9605 

    The thermal rate constant for a chemical reaction, k(T), can be expressed as the long time limit of the flux-side correlation C[sub fs](t)=tr[e[sup -βĤ/2]&Fcirc;e[sup -βĤ/2]e[sup iĤ/2]e[sup iĤt/h]ĥe[sup -iĤt/h]].. Previous work has focused on semiclassical (SC)...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics