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

Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software

/pdf/advanced-capabilities-for-materials-modelling-with-quantum-40poaosrve.pdf

Advanced capabilities for materials modelling with Quantum ESPRESSO.

Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudo-potential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement theirs ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

/pdf/advanced-capabilities-for-materials-modelling-with-quantum-1mm65lmkt8.pdf

Advanced capabilities for materials modelling with Quantum ESPRESSO

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.

/pdf/advances-in-methods-and-algorithms-in-a-modern-quantum-29c5vso11m.pdf

Advances in methods and algorithms in a modern quantum chemistry program package

The Feynman-Kleinert linearized path integral molecular dynamics (FK-LPI), ring polymer molecular dynamics (RPMD), and centroid molecular dynamics (CMD) methods are applied to the simulation of normal liquid helium. Comparisons of the simulation results at the T = 4 K and rho = 0.01873 A-3 state point are presented. The calculated quantum correlation functions for the three methods show significant differences, both in the short time and in the intermediate regions of the spectrum. Our simulation results are also compared to the recent results of other approximate quantum simulation methods. We find that FK-LPI qualitatively agrees with other approximate quantum simulation results while CMD and RPMD predict a qualitatively different impulsive rebound in the velocity autocorrelation function. Frequency space analysis reveals that RPMD exhibits a broad high-frequency tail similar to that from quantum mode coupling theory and numerical analytic continuation approaches, while FK-LPI provides a somewhat more rapid decay at high frequency than any of these three methods. CMD manifests a high-frequency component that is greatly reduced compared with the other methods.

Comparison of approximate quantum simulation methods applied to normal liquid helium at 4 K.

Nuclear Quantum Effects in H$^+$ and OH$^-$ Diffusion Along Confined Water Wires from Ab Initio Path Integral Molecular Dyanmics

We shall use this introduction to the Faraday Discussion on quantum effects in complex systems to review the recent progress that has been made in using imaginary time path integral methods to calculate chemical reaction rates. As a result of this progress, it is now routinely possible to calculate accurate rate constants including quantum mechanical zero point energy and tunnelling effects for arbitrarily complex (anharmonic and multi-dimensional) systems. This is true in the adiabatic (Born–Oppenheimer) limit, in the non-adiabatic (Fermi Golden Rule) limit, and everywhere between these two limits in the normal Marcus regime. Quantum mechanical effects on reaction rates can be enormous, even at room temperature, and the problem of including these effects in simulations of a wide variety of chemical reactions in complex systems has now effectively been solved.

/pdf/path-integral-methods-for-reaction-rates-in-complex-systems-avexhdai9k.pdf

Path integral methods for reaction rates in complex systems.

Magnetic field effects on radical pair reactions arise due to the interplay of coherent electron spin dynamics and spin relaxation effects, a rigorous treatment of which requires the solution of the Liouville-von Neumann equation. However, it is often found that simple incoherent kinetic models of the radical pair singlet-triplet intersystem crossing provide an acceptable description of experimental measurements. In this paper, we outline the theoretical basis for this incoherent kinetic description, elucidating its connection to exact quantum mechanics. We show, in particular, how the finite lifetime of the radical pair spin states, as well as any additional spin-state dephasing, leads to incoherent intersystem crossing. We arrive at simple expressions for the radical pair spin state interconversion rates to which the functional form proposed recently by Steiner et al. [J. Phys. Chem. C 122, 11701 (2018)] can be regarded as an approximation. We also test the kinetic master equation against exact quantum dynamical simulations for a model radical pair and for a series of PTZ•+–Phn–PDI•− molecular wires.

Radical pair intersystem crossing: Quantum dynamics or incoherent kinetics?

The log derivative version of the Kohn variational principle is reviewed in the context of a general inelastic molecular collision. Particular emphasis is placed on the possibility of solving the resulting linear equations iteratively for a single initial state column of the scattering matrix, S, and several important practical observations are made in this regard. In particular, it is argued that the discrete variable representation proposed by Light and coworkers leads to an extremely sparse coefficient matrix, and so has distinct advantages over the more obvious variational basis approach. The resulting methodology is applied to the diffractive scattering of a beam of helium atoms from the (001) face of LiF. Here test calculations with up to 2601 coupled channels and 216 translational grid points clearly demonstrate the practical potential of the iterative technique. The implications of these tests for more general scattering problems are also briefly discussed.

David E. Manolopoulos

Papers

Comparison of approximate quantum simulation methods applied to normal liquid helium at 4 K.

Nuclear Quantum Effects in H$^+$ and OH$^-$ Diffusion Along Confined Water Wires from Ab Initio Path Integral Molecular Dyanmics

Path integral methods for reaction rates in complex systems.

Radical pair intersystem crossing: Quantum dynamics or incoherent kinetics?

Iterative solution in quantum scattering theory. The log derivative Kohn approach