Author
Shuichi Nosé
Other affiliations: National Research Council
Bio: Shuichi Nosé is an academic researcher from Keio University. The author has contributed to research in topics: Molecular dynamics & Glass transition. The author has an hindex of 16, co-authored 37 publications receiving 23339 citations. Previous affiliations of Shuichi Nosé include National Research Council.
Papers
More filters
TL;DR: In this article, the authors compared the canonical distribution in both momentum and coordinate space with three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 77, 63 (1983); and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983).
Abstract: Three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 48, 1818 (1982)], and Evans and Morriss [Chem. Phys. 77, 63 (1983)]; and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983)] are examined analytically via calculating the equilibrium distribution functions and comparing them with that of the canonical ensemble. Except for effects due to momentum and angular momentum conservation, method (1) yields the rigorous canonical distribution in both momentum and coordinate space. Method (2) can be made rigorous in coordinate space, and can be derived from method (1) by imposing a specific constraint. Method (3) is not rigorous and gives a deviation of order N−1/2 from the canonical distribution (N the number of particles). The results for the constant temperature–constant pressure ensemble are similar to the canonical ensemble case.
13,921 citations
TL;DR: In this paper, a molecular dynamics simulation method which can generate configurations belonging to the canonical (T, V, N) ensemble or the constant temperature constant pressure ensemble was proposed, which is tested for an atomic fluid (Ar) and works well.
Abstract: A molecular dynamics simulation method which can generate configurations belonging to the canonical (T, V, N) ensemble or the constant temperature constant pressure (T, P, N) ensemble, is proposed. The physical system of interest consists of N particles (f degrees of freedom), to which an external, macroscopic variable and its conjugate momentum are added. This device allows the total energy of the physical system to fluctuate. The equilibrium distribution of the energy coincides with the canonical distribution both in momentum and in coordinate space. The method is tested for an atomic fluid (Ar) and works well.
8,110 citations
TL;DR: 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 is discussed in this paper.
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.
2,562 citations
TL;DR: In this article, the canonical distribution is derived analytically as a stationary solution of a generalized Liouville's equation, which expresses the conservation of probability in a phase space.
Abstract: How the canonical distribution is realized in simulations based on deterministic dynamical equations is explained in this review. Basic formulations and their recent extensions of two constant temperature molecular dynamics methods; the constraint and the extended system methods, are discussed. In both methods, the canonical distribution is derived analytically as a stationary solution of a generalized Liouville's equation which expresses the conservation of probability in a phase space. In the constraint method, the total kinetic energy of a system is kept to a constant by imposing a constraint
1,066 citations
TL;DR: In this paper, a molecular dynamics simulation method which can generate configurations belonging to the canonical (T, V, N) ensemble or the constant temperature constant pressure ensemble was proposed, which is tested for an atomic fluid (Ar) and works well.
Abstract: A molecular dynamics simulation method which can generate configurations belonging to the canonical (T, V, N) ensemble or the constant temperature constant pressure (T, P, N) ensemble, is proposed The physical system of interest consists of N particles (f degrees of freedom), to which an external, macroscopic variable and its conjugate momentum are added This device allows the total energy of the physical system to fluctuate The equilibrium distribution of the energy coincides with the canonical distribution both in momentum and in coordinate space The method is tested for an atomic fluid (Ar) and works well
300 citations
Cited by
More filters
TL;DR: In this paper, a method is described to realize coupling to an external bath with constant temperature or pressure with adjustable time constants for the coupling, which can be easily extendable to other variables and to gradients, and can be applied also to polyatomic molecules involving internal constraints.
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
25,256 citations
TL;DR: It is demonstrated that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N), which is comparable to that of a simple truncation method of 10 A or less.
Abstract: 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
17,897 citations
TL;DR: The simulation allows us to study in detail the changes in the structure-property relationship through the metal-semiconductor transition, and a detailed analysis of the local structural properties and their changes induced by an annealing process is reported.
Abstract: We present ab initio quantum-mechanical molecular-dynamics simulations of the liquid-metal--amorphous-semiconductor transition in Ge. Our simulations are based on (a) finite-temperature density-functional theory of the one-electron states, (b) exact energy minimization and hence calculation of the exact Hellmann-Feynman forces after each molecular-dynamics step using preconditioned conjugate-gradient techniques, (c) accurate nonlocal pseudopotentials, and (d) Nos\'e dynamics for generating a canonical ensemble. This method gives perfect control of the adiabaticity of the electron-ion ensemble and allows us to perform simulations over more than 30 ps. The computer-generated ensemble describes the structural, dynamic, and electronic properties of liquid and amorphous Ge in very good agreement with experiment. The simulation allows us to study in detail the changes in the structure-property relationship through the metal-semiconductor transition. We report a detailed analysis of the local structural properties and their changes induced by an annealing process. The geometrical, bonding, and spectral properties of defects in the disordered tetrahedral network are investigated and compared with experiment.
16,744 citations
TL;DR: In this article, the authors compared the canonical distribution in both momentum and coordinate space with three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 77, 63 (1983); and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983).
Abstract: Three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 48, 1818 (1982)], and Evans and Morriss [Chem. Phys. 77, 63 (1983)]; and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983)] are examined analytically via calculating the equilibrium distribution functions and comparing them with that of the canonical ensemble. Except for effects due to momentum and angular momentum conservation, method (1) yields the rigorous canonical distribution in both momentum and coordinate space. Method (2) can be made rigorous in coordinate space, and can be derived from method (1) by imposing a specific constraint. Method (3) is not rigorous and gives a deviation of order N−1/2 from the canonical distribution (N the number of particles). The results for the constant temperature–constant pressure ensemble are similar to the canonical ensemble case.
13,921 citations
TL;DR: The software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s is described, which is a very fast program for molecular dynamics simulation.
Abstract: This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.
13,116 citations