Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

The previously developed particle mesh Ewald method is reformulated in terms of efficient B‐spline interpolation of the structure factors This reformulation allows a natural extension of the method to potentials of the form 1/rp with p≥1 Furthermore, efficient calculation of the virial tensor follows Use of B‐splines in place of Lagrange interpolation leads to analytic gradients as well as a significant improvement in the accuracy We demonstrate that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N) For biomolecular systems with many thousands of atoms this method permits the use of Ewald summation at a computational cost comparable to that of a simple truncation method of 10 A or less

A smooth particle mesh Ewald method

CHARMM (Chemistry at HARvard Macromolecular Mechanics) is a highly flexible computer program which uses empirical energy functions to model macromolecular systems. The program can read or model build structures, energy minimize them by first- or second-derivative techniques, perform a normal mode or molecular dynamics simulation, and analyze the structural, equilibrium, and dynamic properties determined in these calculations. The operations that CHARMM can perform are described, and some implementation details are given. A set of parameters for the empirical energy function and a sample run are included.

CHARMM: A program for macromolecular energy, minimization, and dynamics calculations

6.2.2. Definition of Effective Properties 3064 6.3. Response Properties to Magnetic Fields 3066 6.3.1. Nuclear Shielding 3066 6.3.2. Indirect Spin−Spin Coupling 3067 6.3.3. EPR Parameters 3068 6.4. Properties of Chiral Systems 3069 6.4.1. Electronic Circular Dichroism (ECD) 3069 6.4.2. Optical Rotation (OR) 3069 6.4.3. VCD and VROA 3070 7. Continuum and Discrete Models 3071 7.1. Continuum Methods within MD and MC Simulations 3072

/pdf/quantum-mechanical-continuum-solvation-models-3edv963g14.pdf

Quantum mechanical continuum solvation models.

Chemistry in Supercritical Water: Insights from Theory and Simulation

http://link.springer.com/content/pdf/10.1007/978-3-642-85176-6_13.pdf

Electronic-Relaxation Dynamics in Solution

Discretized path integral simulation methods, supplemented by a wavefunction-based representation, has been applied to the determination of aspects of the structure and spectroscopy in two quantum mechanical aqueous systems. The first of these applications is the determination of the consequences of quantizing the molecular degrees of freedom of the water molecules in the pure room temperature liquid. The results provide a quantitative estimate of the significance of approximating such a system as classical and also of the size of isotope effects on the liquid structure. The latter exhibit quantitative differences from recent experiments, indicating limitations in the model used for the solvent. Second, aspects of the structure and spectroscopy of the hydrated electron are considered. Here, we treat the water classically but treat the electron quantum mechanically. The excess electron density and solvent distribution are shown to exhibit structural similarities to ionic solvation. However, it is found that the electronic state fluctuates in response to the fluid, and that these fluctuations have essential consequences for the spectroscopy of the species. In both examples, the high frequency librational motions, and the associated fluctuating fields, that are manifestations of the light hydrogen mass, polarity, and strong restoring forces characteristic of molecular displacements in hydrogen bonded liquids are found to be key ingredients in the description of the phenomena.

Quantum Effects in Hydrogen Bonded Liquids

We investigate the thermodynamics of hydrophobic oligomer collapse using a water-explicit, three-dimensional lattice model. The model captures several aspects of protein thermodynamics, including the emergence of cold- and thermal-unfolding, as well as unfolding at high solvent density (a phenomenon akin to pressure-induced denaturation). We show that over a range of conditions spanning a ≈14% increase in solvent density, the oligomer transforms into a compact, strongly water-penetrated conformation at low temperature. This contrasts with thermal unfolding at high temperature, where the system “denatures” into an extended random coil conformation. We report a phase diagram for hydrophobic collapse that correctly captures qualitative aspects of cold and thermal unfolding at low to intermediate solvent densities.

Phase Behavior of a Lattice Hydrophobic Oligomer in Explicit Water

We examine the range of validity of the Gaussian model for various water-like liquids whose intermolecular potentials differ from SPC/E water, to provide insight into the temperature dependence of the hydrophobic effect for small hard sphere solutes. We find that low compressibility liquids that have more close-packed network structures show much larger deviations from Gaussian fluctuations for low or zero occupancies relative to more compressible fluids with more open networks. Water appears to be a unique molecular fluid in possessing equilibrium density fluctuations that are faithfully described by the Gaussian theory. We ascribe this success to the fact, shown here, that the orientational correlations near a small hard sphere solute involve remarkably little reorganization from the bulk, which is a consequence of water's low solvent reorganization enthalpy and entropy.

/pdf/predicting-cavity-formation-free-energy-how-far-is-the-2xbkxsa3rf.pdf

Peter J. Rossky

Papers

Chemistry in Supercritical Water: Insights from Theory and Simulation

Electronic-Relaxation Dynamics in Solution

Quantum Effects in Hydrogen Bonded Liquids

Phase Behavior of a Lattice Hydrophobic Oligomer in Explicit Water

Predicting cavity formation free energy: how far is the Gaussian approximation valid?