Biomedical subjects

Christopher J Fennell

Publications and source records attributed to Christopher J Fennell.

3 recordsLinked to original sources

PubMed2006

Is the Ewald summation still necessary? Pairwise alternatives to the accepted standard for long-range electrostatics.

We investigate pairwise electrostatic interaction methods and show that there are viable computationally efficient (O(N)) alternatives to the Ewald summation for typical modern molecular simulations. These methods are extended from the damped and cutoff-neutralized Coulombic sum originally proposed by Wolf et al. [J. Chem. Phys. 110, 8255 (1999)]. One of these, the damped shifted force method, shows a remarkable ability to reproduce the energetic and dynamic characteristics exhibited by simulations employing lattice summation techniques. Comparisons were performed with this and other pairwise methods against the smooth particle-mesh Ewald summation to see how well they reproduce the energetics and dynamics of a variety of molecular simulations.

Journal Article↗

PubMed2005

OOPSE: an object-oriented parallel simulation engine for molecular dynamics.

OOPSE is a new molecular dynamics simulation program that is capable of efficiently integrating equations of motion for atom types with orientational degrees of freedom (e.g. "sticky" atoms and point dipoles). Transition metals can also be simulated using the embedded atom method (EAM) potential included in the code. Parallel simulations are carried out using the force-based decomposition method. Simulations are specified using a very simple C-based meta-data language. A number of advanced integrators are included, and the basic integrator for orientational dynamics provides substantial improvements over older quaternion-based schemes.

Journal Article↗

PubMed2004

On the structural and transport properties of the soft sticky dipole and related single-point water models.

The density maximum and temperature dependence of the self-diffusion constant were investigated for the soft sticky dipole (SSD) water model and two related reparametrizations of this single-point model. A combination of microcanonical and isobaric-isothermal molecular dynamics simulations was used to calculate these properties, both with and without the use of reaction field to handle long-range electrostatics. The isobaric-isothermal simulations of the melting of both ice-Ih and ice-Ic showed a density maximum near 260 K. In most cases, the use of the reaction field resulted in calculated densities which were significantly lower than experimental densities. Analysis of self-diffusion constants shows that the original SSD model captures the transport properties of experimental water very well in both the normal and supercooled liquid regimes. We also present our reparametrized versions of SSD for use both with the reaction field or without any long-range electrostatic corrections. These are called the SSD/RF and SSD/E models, respectively. These modified models were shown to maintain or improve upon the experimental agreement with the structural and transport properties that can be obtained with either the original SSD or the density-corrected version of the original model (SSD1). Additionally, a novel low-density ice structure is presented which appears to be the most stable ice structure for the entire SSD family.

Journal Article↗