Michael L. Klein

Bio: Michael L. Klein is an academic researcher from Temple University. The author has contributed to research in topics: Molecular dynamics & Solvation. The author has an hindex of 104, co-authored 745 publications receiving 78805 citations. Previous affiliations of Michael L. Klein include National Research Council & University of Tübingen.

Comparison of simple potential functions for simulating liquid water

Abstract: Classical Monte Carlo simulations have been carried out for liquid water in the NPT ensemble at 25 °C and 1 atm using six of the simpler intermolecular potential functions for the water dimer: Bernal–Fowler (BF), SPC, ST2, TIPS2, TIP3P, and TIP4P. Comparisons are made with experimental thermodynamic and structural data including the recent neutron diffraction results of Thiessen and Narten. The computed densities and potential energies are in reasonable accord with experiment except for the original BF model, which yields an 18% overestimate of the density and poor structural results. The TIPS2 and TIP4P potentials yield oxygen–oxygen partial structure functions in good agreement with the neutron diffraction results. The accord with the experimental OH and HH partial structure functions is poorer; however, the computed results for these functions are similar for all the potential functions. Consequently, the discrepancy may be due to the correction terms needed in processing the neutron data or to an effect uniformly neglected in the computations. Comparisons are also made for self‐diffusion coefficients obtained from molecular dynamics simulations. Overall, the SPC, ST2, TIPS2, and TIP4P models give reasonable structural and thermodynamic descriptions of liquid water and they should be useful in simulations of aqueous solutions. The simplicity of the SPC, TIPS2, and TIP4P functions is also attractive from a computational standpoint.

Constant pressure molecular dynamics algorithms

Abstract: Modularly invariant equations of motion are derived that generate the isothermal–isobaric ensemble as their phase space averages. Isotropic volume fluctuations and fully flexible simulation cells as well as a hybrid scheme that naturally combines the two motions are considered. The resulting methods are tested on two problems, a particle in a one‐dimensional periodic potential and a spherical model of C60 in the solid/fluid phase.

Nosé-Hoover chains : the canonical ensemble via continuous dynamics

Abstract: Nose has derived a set of dynamical equations that can be shown to give canonically distributed positions and momenta provided the phase space average can be taken into the trajectory average, i.e., the system is ergodic [S. Nose, J. Chem. Phys. 81, 511 (1984), W. G. Hoover, Phys. Rev. A 31, 1695 (1985)]. Unfortunately, the Nose–Hoover dynamics is not ergodic for small or stiff systems. Here a modification of the dynamics is proposed which includes not a single thermostat variable but a chain of variables, Nose–Hoover chains. The ‘‘new’’ dynamics gives the canonical distribution where the simple formalism fails. In addition, the new method is easier to use than an extension [D. Kusnezov, A. Bulgac, and W. Bauer, Ann. Phys. 204, 155 (1990)] which also gives the canonical distribution for stiff cases.

Constant pressure molecular dynamics for molecular systems

Abstract: Technical aspects of the constant pressure molecular dynamics (MD) method proposed by Andersen and extended by Parrinello and Rahman to allow changes in the shape of the MD cell are discussed. The new MD method is extended to treat molecular systems and to include long range charge-charge interactions. Results on the conservation laws, the frequency of oscillation of the MD cell, and the equations which constrain the shape of the MD cell are also given. An additional constraint is introduced to stop the superfluous MD cell rotation which would otherwise complicate the analysis of crystal structures. The method is illustrated by examining the behaviour of solid nitrogen at high pressure.

Explicit reversible integrators for extended systems dynamics

Abstract: Explicit reversible integrators, suitable for use in large-scale computer simulations, are derived for extended systems generating the canonical and isothermal-isobaric ensembles. The new methods are compared with the standard implicit (iterative) integrators on some illustrative example problems. In addition, modification of the proposed algorithms for multiple time step integration is outlined.

Comparison of simple potential functions for simulating liquid water

Abstract: Classical Monte Carlo simulations have been carried out for liquid water in the NPT ensemble at 25 °C and 1 atm using six of the simpler intermolecular potential functions for the water dimer: Bernal–Fowler (BF), SPC, ST2, TIPS2, TIP3P, and TIP4P. Comparisons are made with experimental thermodynamic and structural data including the recent neutron diffraction results of Thiessen and Narten. The computed densities and potential energies are in reasonable accord with experiment except for the original BF model, which yields an 18% overestimate of the density and poor structural results. The TIPS2 and TIP4P potentials yield oxygen–oxygen partial structure functions in good agreement with the neutron diffraction results. The accord with the experimental OH and HH partial structure functions is poorer; however, the computed results for these functions are similar for all the potential functions. Consequently, the discrepancy may be due to the correction terms needed in processing the neutron data or to an effect uniformly neglected in the computations. Comparisons are also made for self‐diffusion coefficients obtained from molecular dynamics simulations. Overall, the SPC, ST2, TIPS2, and TIP4P models give reasonable structural and thermodynamic descriptions of liquid water and they should be useful in simulations of aqueous solutions. The simplicity of the SPC, TIPS2, and TIP4P functions is also attractive from a computational standpoint.

Fast parallel algorithms for short-range molecular dynamics

Abstract: 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.

Molecular dynamics with coupling to an external bath.

Abstract: In molecular dynamics (MD) simulations the need often arises to maintain such parameters as temperature or pressure rather than energy and volume, or to impose gradients for studying transport properties in nonequilibrium MD A method is described to realize coupling to an external bath with constant temperature or pressure with adjustable time constants for the coupling The method is easily extendable to other variables and to gradients, and can be applied also to polyatomic molecules involving internal constraints The influence of coupling time constants on dynamical variables is evaluated A leap‐frog algorithm is presented for the general case involving constraints with coupling to both a constant temperature and a constant pressure bath

QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials

Abstract: QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

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

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