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.
In this paper identities are derived which allow the computation of the Coulomb energy associated with N charges in a central cell and all their periodic images.These identities are all consequences of one basic identity which is obtained in a simple and straightforward way. It is possible to extend the results to other types of potent i a l s a s w ell.
A summation formula is proven for a lattice sum occurring frequently in molecular dynamics calculations. It has a much faster convergence rate than the original sum. An important application is in the case of Coulomb forces.
We study a general class of scalar reaction/interacting population diffusion equations in two space dimensions: convective terms, due to wind, are included. We consider boundary conditions which include a measure of the hostility to the species in the exterior of the domain. The main point of the paper is to obtain estimates for the minimum domain size which can sustain spatially heterogeneous structures and indicate the type of patterns which could appear.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.