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 unpolarized absorption and circular dichroism spectra of the fundamental vibrational transitions of the chiral molecule, 4-methyl-2-oxetanone, are calculated ab initio. Harmonic force fields are obtained using Density Functional Theory (DFT), MP2, and SCF methodologies and a 5S4P2D/3S2P (TZ2P) basis set. DFT calculations use the Local Spin Density Approximation (LSDA), BLYP, and Becke3LYP (B3LYP) density functionals. Mid-IR spectra predicted using LSDA, BLYP, and B3LYP force fields are of significantly different quality, the B3LYP force field yielding spectra in clearly superior, and overall excellent, agreement with experiment. The MP2 force field yields spectra in slightly worse agreement with experiment than the B3LYP force field. The SCF force field yields spectra in poor agreement with experiment.The basis set dependence of B3LYP force fields is also explored: the 6-31G* and TZ2P basis sets give very similar results while the 3-21G basis set yields spectra in substantially worse agreements with experiment. jg

Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields

Advanced polyimide materials: Syntheses, physical properties and applications

VASPKIT: A user-friendly interface facilitating high-throughput computing and analysis using VASP code

An activity coefficient model using molecular solvation based on the COSMO-RS method is proposed. In this model, quantum mechanical COSMO calculations are performed to obtain the screening charges for molecules in a perfect conductor. A statistical mechanical model that considers molecules to be a collection of surface segments is developed for the calculation of segment activity coefficients using these screening charges. Activity coefficients for molecules are then obtained by summing the contributions of the segments. This model requires only a single radius for each atom in the COSMO solvation calculations, one universal parameter to discern hydrogen-bond acceptors and donors, and two universal parameters to determine segment interactions. This is a significantly fewer number of parameters for phase equilibrium calculations than group contribution methods such as the UNIFAC (168 parameters) and modified UNIFAC (612 parameters) models. The resulting completely a priori prediction method results in abso...

A Priori Phase Equilibrium Prediction from a Segment Contribution Solvation Model

To validate this two phase model for the thermodynamics of liquids, we applied it to pure Lennard-Jones systems for a range of reduced temperatures from 0.9 to 1.8 and reduced densities from 0.05 to 1.10. These conditions cover the gas, liquid, crystal, metastable, and unstable states in the phase diagram. Our results compare quite well with accurate Monte Carlo calculations of the phase diagram for classical Lennard-Jones particles throughout the entire phase diagram. Thus the two-phase thermodynamics approach provides an efficient means for extracting thermodynamic properties of liquids ~and gases and solids! .© 2003 American Institute of Physics. @DOI: 10.1063/1.1624057#

/pdf/the-two-phase-model-for-calculating-thermodynamic-properties-1vku2doglh.pdf

The two-phase model for calculating thermodynamic properties of liquids from molecular dynamics: Validation for the phase diagram of Lennard-Jones fluids

We report various structural and conformational properties of generations 4, 5, and 6 PAMAM (polyamidoamine) dendrimer [EDA (ethylenediamine) core)] at various protonation levels through extensive molecular dynamics (MD) simulations in explicit solvent. The presence of solvent leads to swelling of the dendrimer (by 33% for G5 compared to the case of no solvent). We find that decreasing the solution from high pH (∼10, no protonation) to neutral (∼7, only primary amines protonated) to low pH (∼4, tertiary amines also protonated) changes the radius of gyration of G5 from 21 to 22 to 25 A, respectively. We also report such other structural quantities as radial density, distribution of terminal groups, solvent accessible surface area and volume, shape, and structure factors (to compare with SAXS and SANS experiments) at various pH conditions. We find significant back-folding of the outer subgenerations in the interior of the molecules at all levels of pH, contrary to original expectations and some SANS experiments but in agreement with other SANS experiments. We find significant water penetration inside the dendrimer, with ∼3 water/tertiary amine for high pH and ∼6 water/tertiary amine for low pH (all for G5). This indicates that the interior of the dendrimer is quite open with internal cavities available for accommodating guest molecules, suggesting using PAMAM dendrimer for guest−host applications. This estimate of internal waters suggests that sufficient water is available to facilitate metal ion binding.

Effect of Solvent and pH on the Structure of PAMAM Dendrimers

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Shiang-Tai Lin

Papers

A Priori Phase Equilibrium Prediction from a Segment Contribution Solvation Model

The two-phase model for calculating thermodynamic properties of liquids from molecular dynamics: Validation for the phase diagram of Lennard-Jones fluids

Effect of Solvent and pH on the Structure of PAMAM Dendrimers

Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Improvements of COSMO-SAC for vapor-liquid and liquid-liquid equilibrium predictions