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

NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This article, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

/pdf/scalable-molecular-dynamics-with-namd-51o1z7rwp8.pdf

Scalable molecular dynamics with NAMD

/pdf/scalable-molecular-dynamics-with-namd-3j7xyc70g4.pdf

Scalable Molecular Dynamics with NAMD

CASTEP Computer program / Density functional theory / Pseudopotentials / ab initio study / Plane-wave method / Computational crystallography Abstract. The CASTEP code for first principles electro- nic structure calculations will be described. A brief, non- technical overview will be given and some of the features and capabilities highlighted. Some features which are un- ique to CASTEP will be described and near-future devel- opment plans outlined.

/pdf/first-principles-methods-using-castep-2ke0q7ufxk.pdf

First principles methods using CASTEP

Calcium carbonate is an abundant substance that can be created in several mineral forms by the reaction of dissolved carbon dioxide in water with calcium ions. Through biomineralization, organisms can harness and control this process to form various functional materials that can act as anything from shells through to lenses. The early stages of calcium carbonate formation have recently attracted attention as stable prenucleation clusters have been observed, contrary to classical models. Here we show, using computer simulations combined with the analysis of experimental data, that these mineral clusters are made of an ionic polymer, composed of alternating calcium and carbonate ions, with a dynamic topology consisting of chains, branches and rings. The existence of a disordered, flexible and strongly hydrated precursor provides a basis for explaining the formation of other liquid-like amorphous states of calcium carbonate, in addition to the non-classical behaviour during growth of amorphous calcium carbonate.

/pdf/stable-prenucleation-mineral-clusters-are-liquid-like-ionic-1929lhc1t2.pdf

Stable prenucleation mineral clusters are liquid-like ionic polymers

The performance of existing force-field models for the calcium carbonate - water system has been critically assessed with particular reference to the thermodynamic consequences. It is demonstrated that all currently available parametrizations fail to describe the calcite-aragonite phase transition, and the free energies of solvation for the calcium cation are also considerably in error leading to a poor description of the dissolution enthalpy for calcite. A new force-field, based on rigid carbonate ions, has been developed that corrects these deficiencies and accurately describes the thermodynamics of the aqueous calcium carbonate system within molecular dynamics simulations. Not only does this new model lead to quantitative changes in the properties of the calcite (1014) surface in contact with water, but also significant qualitative differences. With this more accurate model it is found that calcium ions do not adsorb at the pristine basal plane of calcite, while carbonate ions only weakly bind. Carbon...

Derivation of an accurate force-field for simulating the growth of calcium carbonate from aqueous solution : a new model for the calcite-water interface

Nucleation of gas hydrates remains a poorly understood phenomenon, despite its importance as a critical step in understanding the performance and mode of action of low dosage hydrate inhibitors. We present here a detailed analysis of the structural and mechanistic processes by which gas hydrates nucleate in a molecular dynamics simulation of dissolved methane at a methane/water interface. It was found that hydrate initially nucleates into a phase consistent with none of the common bulk crystal structures, but containing structural units of all of them. The process of water cage formation has been found to correlate strongly with the collective arrangement of methane molecules.

Gas hydrate nucleation and cage formation at a water/methane interface

Nucleation and crystal growth are important in material synthesis, climate modeling, biomineralization, and pharmaceutical formulation. Despite tremendous efforts, the mechanisms and kinetics of nucleation remain elusive to both theory and experiment. Here we investigate sodium chloride (NaCl) nucleation from supersaturated brines using seeded atomistic simulations, polymorph-specific order parameters, and elements of classical nucleation theory. We find that NaCl nucleates via the common rock salt structure. Ion desolvation-not diffusion-is identified as the limiting resistance to attachment. Two different analyses give approximately consistent attachment kinetics: diffusion along the nucleus size coordinate and reaction-diffusion analysis of approach-to-coexistence simulation data from Aragones et al. ( J. Chem. Phys. 2012, 136, 244508 ). Our simulations were performed at realistic supersaturations to enable the first direct comparison to experimental nucleation rates for this system. The computed and measured rates converge to a common upper limit at extremely high supersaturation. However, our rate predictions are between 15 and 30 orders of magnitude too fast. We comment on possible origins of the large discrepancy.

Nucleation of NaCl from Aqueous Solution: Critical Sizes, Ion-Attachment Kinetics, and Rates.

The metadynamics method for accelerating rate events in molecular simulations is applied to the problem of ice freezing. We demonstrate homogeneous nucleation and growth of ice at 180 K in the isothermal-isobaric ensemble without the presence of external fields or surfaces. This result represents the first report of continuous and dynamic ice nucleation in a system of freely evolving density. Simulations are conducted using a variety of periodic simulation domains. In all cases the cubic polymorph ice Ic is grown. The influence of boundary effects on estimates of the nucleation free energy barrier are discussed in relation to differences between this and earlier work.

David Quigley

Papers

Stable prenucleation mineral clusters are liquid-like ionic polymers

Derivation of an accurate force-field for simulating the growth of calcium carbonate from aqueous solution : a new model for the calcite-water interface

Gas hydrate nucleation and cage formation at a water/methane interface

Nucleation of NaCl from Aqueous Solution: Critical Sizes, Ion-Attachment Kinetics, and Rates.

Metadynamics simulations of ice nucleation and growth