There is, I think, something ethereal about i —the square root of minus one. I remember first hearing about it at school. It seemed an odd beast at that time—an intruder hovering on the edge of reality.

Usually familiarity dulls this sense of the bizarre, but in the case of i it was the reverse: over the years the sense of its surreal nature intensified. It seemed that it was impossible to write mathematics that described the real world in …

I and i

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

http://kth.diva-portal.org/smash/get/diva2:1303357/FULLTEXT01

GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers

[신간의 별자리x] 우리/미술, 그리고 ‘슬픔의 박물관’

Using molecular dynamics simulations, we show that the thermal conductivity (κ) of Si0.5Ge0.5 can be reduced by more than one order of magnitude by etching nanometer-sized holes in the material, and it becomes almost constant as a function of temperature between 300 and 1100 K for samples with 1 nm wide pores. In nanoporous SiGe, thermal conduction is largely determined by mass disorder and boundary scattering, and thus the dependence of κ on pore distance and on structural, atomistic disorder is much weaker than in the case of nanoporous Si. This indicates that one may minimize κ of the alloy with less stringent morphological constraints than for pure Si.

Morphology and Temperature Dependence of the Thermal Conductivity of Nanoporous SiGe

The Adaptive Resolution Scheme (AdResS) is a hybrid scheme that allows to treat a molecular system with different levels of resolution depending on the location of the molecules. The construction of a Hamiltonian based on the this idea (H-AdResS) allows one to formulate the usual tools of ensembles and statistical mechanics. We present a number of exact and approximate results that provide a statistical mechanics foundation for this simulation method. We also present simulation results that illustrate the theory.

/pdf/statistical-mechanics-of-hamiltonian-adaptive-resolution-4ivowhwe2w.pdf

Statistical mechanics of Hamiltonian adaptive resolution simulations

Water–metal interfaces are ubiquitous and play a key role in many chemical processes, from catalysis to corrosion. Whereas water adlayers on atomically flat transition metal surfaces have been investigated in depth, little is known about the chemistry of water on stepped surfaces, commonly occurring in realistic situations. Using first-principles simulations, we study the adsorption of water on a stepped platinum surface. We find that water adsorbs preferentially at the step edge, forming linear clusters or chains, stabilized by the cooperative effect of chemical bonds with the substrate and hydrogen bonds. In contrast with their behavior on flat Pt, at steps water molecules dissociate, forming mixed hydroxyl/water structures, through an autocatalytic mechanism promoted by H-bonding. Nuclear quantum effects contribute to stabilize partially dissociated cluster and chains. Together with the recently demonstrated behavior of water chains adsorbed on stepped Pt surfaces to transfer protons via thermally acti...

/pdf/autocatalytic-and-cooperatively-stabilized-dissociation-of-2jepdvhezb.pdf

Autocatalytic and Cooperatively Stabilized Dissociation of Water on a Stepped Platinum Surface

We outline a method to investigate the role of nuclear quantum effects in liquid water making use of a force field derived from ab initio simulations. Starting from a first-principles molecular dynamics simulation, we obtain an effective force field for bulk liquid water using the force-matching technique. After validating that our effective model reproduces the key structural and dynamic properties of the reference system, we use it to perform path integral simulations to investigate the role played by nuclear quantum effects on bulk water, probing radial distribution functions, vibrational spectra, and hydrogen bond fluctuations. Our approach offers a practical route to derive ab initio quality molecular models to study quantum effects at a low computational cost.

/pdf/nuclear-quantum-effects-in-water-a-multiscale-study-x9169xa8wm.pdf

Nuclear Quantum Effects in Water: A Multiscale Study.

We use classical molecular dynamics combined with the recently developed metadynamics method [A. Laio and M. Parrinello, Procs. Natl. Acad. Sci. USA 99, 20 (2002)] to study the process of bulk melting in hexagonal ice. Our simulations show that bulk melting is mediated by the formation of topological defects which preserve the coordination of the tetrahedral network. Such defects cluster to form a defective region involving about 50 molecules with a surprisingly long life-time. The subsequent formation of coordination defects triggers the transition to the liquid state.

Davide Donadio

Papers

Morphology and Temperature Dependence of the Thermal Conductivity of Nanoporous SiGe

Statistical mechanics of Hamiltonian adaptive resolution simulations

Autocatalytic and Cooperatively Stabilized Dissociation of Water on a Stepped Platinum Surface

Nuclear Quantum Effects in Water: A Multiscale Study.

Topological defects and bulk melting of hexagonal ice