https://www.sandia.gov/~sjplimp/papers/jcompphys95.pdf

Fast parallel algorithms for short-range molecular dynamics

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.

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

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

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

Scalable Molecular Dynamics with NAMD

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

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

The topological constraints of nonconcatenated ring polymers force them to form compact loopy globular conformations with much lower entropy than unconstrained ideal rings. The closed-loop structure of ring polymers also enables them to be threaded by linear polymers in ring/linear blends, resulting in less compact ring conformations with higher entropy. This conformational entropy increase promotes mixing rings with linear polymers. Here, using molecular dynamics simulations for bead-spring chains, ring/linear blends are shown to be significantly more miscible than linear/linear blends and that there is an entropic mixing, negative χ, for ring/linear blends compared to linear/linear and ring/ring blends. In analogy with small angle neutron scattering, the static structure function S(q) is measured, and the resulting data are fit to the random phase approximation model to determine χ. In the limit that the two components are the same, χ = 0 for the linear/linear and ring/ring blends as expected, while χ < 0 for the ring/linear blends. With increasing chain stiffness, χ for the ring/linear blends becomes more negative, varying reciprocally with the number of monomers between entanglements. Ring/linear blends are also shown to be more miscible than either ring/ring or linear/linear blends and stay in single phase for a wider range of increasing repulsion between the two components.

https://doi.org/10.1021/acspolymersau.2c00050

Entropic Mixing of Ring/Linear Polymer Blends

A parallel algorithm for the concurrent atomistic-continuum methodology

The final report for a Laboratory Directed Research and Development project entitled, ''From Atom-Picoseconds to Centimeter-Years in Simulation and Experiment'' is presented. In this project, separate modeling methods at the atomic scale were used to bridge gaps in time and space with higher scales. For understanding of continuum mechanics quantities at various scales atomistic simulations that ranged from nanometers to microns were performed and experiments from centimeters to millimeters were performed. Certain continuum mechanical quantities were clearly defined as a function of size scale, for example, the yield stress. Several techniques were used to extend the time scale of simulations, including calculating prefactors and activation energies for diffusion events and mapping complex atomic motions onto more tractable lattice models. In the case of transport of small molecules in polymeric and nanoporous materials, new Monte Carlo methods for sampling transition rates were developed.

/pdf/from-atom-picoseconds-to-centimeter-years-in-simulation-and-2cneoruv0g.pdf

From Atom-Picoseconds to Centimeter-Years in Simulation and Experiment

Many of the most important and hardest-to-solve problems related to the synthesis, performance, and aging of materials involve diffusion through the material or along surfaces and interfaces. These diffusion processes are driven by motions at the atomic scale, but traditional atomistic simulation methods such as molecular dynamics are limited to very short timescales on the order of the atomic vibration period (less than a picosecond), while macroscale diffusion takes place over timescales many orders of magnitude larger. We have completed an LDRD project with the goal of developing and implementing new simulation tools to overcome this timescale problem. In particular, we have focused on two main classes of methods: accelerated molecular dynamics methods that seek to extend the timescale attainable in atomistic simulations, and so-called 'equation-free' methods that combine a fine scale atomistic description of a system with a slower, coarse scale description in order to project the system forward over long times.

/pdf/accelerated-molecular-dynamics-and-equation-free-methods-for-1vwhx24pvr.pdf

Accelerated molecular dynamics and equation-free methods for simulating diffusion in solids.

This paper describes a method for incorporating a diffusion field modeling oxygen usage and dispersion in a multi-scale model of Mycobacterium tuberculosis (Mtb) infection mediated granuloma formation. We implemented this method over a floating-point field to model oxygen dynamics in host tissue during chronic phase response and Mtb persistence. The method avoids the requirement of satisfying the Courant-Friedrichs-Lewy (CFL) condition, which is necessary in implementing the explicit version of the finite-difference method, but imposes an impractical bound on the time step. Instead, diffusion is modeled by a matrix-based, steady state approximate solution to the diffusion equation. Presented in figure 1 is the evolution of the diffusion profiles of a containment granuloma over time.

Steven J. Plimpton

Papers

Entropic Mixing of Ring/Linear Polymer Blends

A parallel algorithm for the concurrent atomistic-continuum methodology

From Atom-Picoseconds to Centimeter-Years in Simulation and Experiment

Accelerated molecular dynamics and equation-free methods for simulating diffusion in solids.

A method for modeling oxygen diffusion in an agent-based model with application to host-pathogen infection