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.

The optoelectronic properties of amorphous conjugated polymers are sensitive to conformational disorder and spectroscopy provides the means for structural characterization of the fragments of the chain which interact with light - "chromophores". A faithful interpretation of spectroscopic conformational signatures, however, presents a key challenge. We investigate the relationship between the ground state optical gaps, the properties of the excited states, and the structural features of chromophores of a single molecule poly(3-hexyl)-thiophene (P3HT), using quantum-classical atomistic simulations. Our results demonstrate that chromophoric disorder reflects an interplay between excited state de-localization and electron-hole polarization, and is controlled by torsional disorder that is specifically associated with the presence of side chains. Within this conceptual framework, we predict and explain a counter-intuitive spectral signature of P3HT: a red-shifted absorption, despite shortening of chromophores, with increasing temperature.

Relating Chromophoric and Structural Disorder in Conjugated Polymers

The role of specific solvent modes in the non-radiative relaxation of an excess electron in methanol

The rapidly advancing area of computer simulation for quantum-mechanical systems now permits the direct evaluation of spectroscopic observables for condensed phase systems, which increasingly overlaps frontier ultra-fast measurements. In this paper recent theoretical techniques are outlined and representative results from application to the hydrated electron are described that exemplify the close connection between theory and experiment that can now be made. Results focus, first, on the decomposition of the equilibrium optical spectrum into the underlying electronic excitations and on simulation of a proposed photophysical hole-burning experiment designed to test the veracity of this decomposition. Second, the author examines aspects of the transient relaxation of initially nonequilibrium photogenerated electrons by simulation of the spectroscopic signature of proposed early time states and of the relaxation of these states to equilibrium.

Quantum simulation of ultrafast spectroscopy

Coarse-grained modeling of conjugated polymers has become an increasingly popular route to investigate the physics of organic optoelectronic materials. While ultraviolet (UV)-vis spectroscopy remains one of the key experimental methods for the interrogation of these materials, a rigorous bridge between simulated coarse-grained structures and spectroscopy has not been established. Here, we address this challenge by developing a method that can predict spectra of conjugated polymers directly from coarse-grained representations while avoiding repetitive procedures such as ad hoc back-mapping from coarse-grained to atomistic representations followed by spectral computation using quantum chemistry. Our approach is based on a generative deep-learning model: the long-short-term memory recurrent neural network (LSTM-RNN). The latter is suggested by the apparent similarity between natural languages and the mathematical structure of perturbative expansions of, in our case, excited-state energies perturbed by conformational fluctuations. We also use this model to explore the level of sensitivity of spectra to the coarse-grained representation back-mapping protocol. Our approach presents a tool uniquely suited for improving postsimulation analysis protocols, as well as, potentially, for including spectral data as input in the refinement of coarse-grained potentials.

https://www.pnas.org/content/pnas/117/25/13945.full.pdf

Predicting optical spectra for optoelectronic polymers using coarse-grained models and recurrent neural networks.

We have applied the extended reference interaction site model (RISM) theory with the hypernetted chain (HNC) closure to evaluate the solvent structure and potential of mean force (PMF) for the sodium dimethyl phosphate (NaDMP) ion pair in «simple point charge» (SPC) water. Comparison is made between the free energies obtained from the RISM-HNC theory, from the so-called Gaussian bath theory using RISM solvent structure, and from computer simulation. Both the extended RISM theory and Gaussian bath theory predict the locations of PMF extrema in reasonable accord with those calculated by molecular dynamics simulation, although the PMFs obtained from the extended RISM theory show smaller amplitudes. Overall, the Gaussian bath theory yields the PMF curves in better accord with simulation

Peter J. Rossky

Papers

Relating Chromophoric and Structural Disorder in Conjugated Polymers

The role of specific solvent modes in the non-radiative relaxation of an excess electron in methanol

Quantum simulation of ultrafast spectroscopy

Predicting optical spectra for optoelectronic polymers using coarse-grained models and recurrent neural networks.

Potential of mean force for a sodium dimethyl phosphate ion pair in aqueous solution: a further test of the extended RISM theory