scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Polymorphic transitions in single crystals: A new molecular dynamics method

01 Dec 1981-Journal of Applied Physics (American Institute of PhysicsAIP)-Vol. 52, Iss: 12, pp 7182-7190
TL;DR: In this paper, a new Lagrangian formulation is introduced to make molecular dynamics (MD) calculations on systems under the most general externally applied, conditions of stress, which is well suited to the study of structural transformations in solids under external stress and at finite temperature.
Abstract: A new Lagrangian formulation is introduced. It can be used to make molecular dynamics (MD) calculations on systems under the most general, externally applied, conditions of stress. In this formulation the MD cell shape and size can change according to dynamical equations given by this Lagrangian. This new MD technique is well suited to the study of structural transformations in solids under external stress and at finite temperature. As an example of the use of this technique we show how a single crystal of Ni behaves under uniform uniaxial compressive and tensile loads. This work confirms some of the results of static (i.e., zero temperature) calculations reported in the literature. We also show that some results regarding the stress‐strain relation obtained by static calculations are invalid at finite temperature. We find that, under compressive loading, our model of Ni shows a bifurcation in its stress‐strain relation; this bifurcation provides a link in configuration space between cubic and hexagonal close packing. It is suggested that such a transformation could perhaps be observed experimentally under extreme conditions of shock.
Citations
More filters
Journal ArticleDOI
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

Journal ArticleDOI
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

Journal ArticleDOI
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

Journal ArticleDOI
TL;DR: A range of new simulation algorithms and features developed during the past 4 years are presented, leading up to the GROMACS 4.5 software package, which provides extremely high performance and cost efficiency for high-throughput as well as massively parallel simulations.
Abstract: Motivation: Molecular simulation has historically been a low-throughput technique, but faster computers and increasing amounts of genomic and structural data are changing this by enabling large-scale automated simulation of, for instance, many conformers or mutants of biomolecules with or without a range of ligands. At the same time, advances in performance and scaling now make it possible to model complex biomolecular interaction and function in a manner directly testable by experiment. These applications share a need for fast and efficient software that can be deployed on massive scale in clusters, web servers, distributed computing or cloud resources. Results: Here, we present a range of new simulation algorithms and features developed during the past 4 years, leading up to the GROMACS 4.5 software package. The software now automatically handles wide classes of biomolecules, such as proteins, nucleic acids and lipids, and comes with all commonly used force fields for these molecules built-in. GROMACS supports several implicit solvent models, as well as new free-energy algorithms, and the software now uses multithreading for efficient parallelization even on low-end systems, including windows-based workstations. Together with hand-tuned assembly kernels and state-of-the-art parallelization, this provides extremely high performance and cost efficiency for high-throughput as well as massively parallel simulations. Availability: GROMACS is an open source and free software available from http://www.gromacs.org. Contact: erik.lindahl@scilifelab.se Supplementary information:Supplementary data are available at Bioinformatics online.

6,029 citations


Cites methods from "Polymorphic transitions in single c..."

  • ...For pressure control, the existing algorithm by Parrinello and Rahman (Nosé and Klein, 1983; Parrinello and Rahman, 1981) is not correct using leapfrog, although it has been verified to give a correct distribution of volumes within statistical noise in many situations....

    [...]

Journal ArticleDOI
TL;DR: In this paper, a general all-atom force field for atomistic simulation of common organic molecules, inorganic small molecules, and polymers was developed using state-of-the-art ab initio and empirical parametrization techniques.
Abstract: A general all-atom force field for atomistic simulation of common organic molecules, inorganic small molecules, and polymers was developed using state-of-the-art ab initio and empirical parametrization techniques. The valence parameters and atomic partial charges were derived by fitting to ab initio data, and the van der Waals (vdW) parameters were derived by conducting MD simulations of molecular liquids and fitting the simulated cohesive energies and equilibrium densities to experimental data. The combined parametrization procedure significantly improves the quality of a general force field. Validation studies based on large number of isolated molecules, molecular liquids and molecular crystals, representing 28 molecular classes, show that the present force field enables accurate and simultaneous prediction of structural, conformational, vibrational, and thermophysical properties for a broad range of molecules in isolation and in condensed phases. Detailed results of the parametrization and validation f...

4,722 citations

References
More filters
Journal ArticleDOI

23,110 citations

Book
01 Jan 1953
TL;DR: In this paper, the Hartree-Fock Approximation of many-body techniques and the Electron Gas Polarons and Electron-phonon Interaction are discussed.
Abstract: Mathematical Introduction Acoustic Phonons Plasmons, Optical Phonons, and Polarization Waves Magnons Fermion Fields and the Hartree-Fock Approximation Many-body Techniques and the Electron Gas Polarons and the Electron-phonon Interaction Superconductivity Bloch Functions - General Properties Brillouin Zones and Crystal Symmetry Dynamics of Electrons in a Magnetic Field: de Haas-van Alphen Effect and Cyclotron Resonance Magnetoresistance Calculation of Energy Bands and Fermi Surfaces Semiconductor Crystals I: Energy Bands, Cyclotron Resonance, and Impurity States Semiconductor Crystals II: Optical Absorption and Excitons Electrodynamics of Metals Acoustic Attenuation in Metals Theory of Alloys Correlation Functions and Neutron Diffraction by Crystals Recoilless Emission Green's Functions - Application to Solid State Physics Appendix: Perturbation Theory and the Electron Gas Index.

21,954 citations

Book
01 Jan 1934
TL;DR: The theory of the slipline field is used in this article to solve the problem of stable and non-stressed problems in plane strains in a plane-strain scenario.
Abstract: Chapter 1: Stresses and Strains Chapter 2: Foundations of Plasticity Chapter 3: Elasto-Plastic Bending and Torsion Chapter 4: Plastic Analysis of Beams and Frames Chapter 5: Further Solutions of Elasto-Plastic Problems Chapter 6: Theory of the Slipline Field Chapter 7: Steady Problems in Plane Strain Chapter 8: Non-Steady Problems in Plane Strain

20,724 citations

Book
01 Jan 1976
TL;DR: In this article, the authors present a mathematical model for time-dependent correlation functions and response functions in liquid solvers, based on statistical mechanics and molecular distribution functions, and show that these functions are related to time correlation functions in Ionic and Ionic liquids.
Abstract: Introduction. Statistical Mechanics and Molecular Distribution Functions. Computer "Experiments" on Liquids. Diagrammatic Expansions. Distribution Function Theories. Perturbation Theories. Time-dependent Correlation Functions and Response Functions. Hydrodynamics And Transport Coefficients. Microscopic Theories of Time-Correlation Functions. Ionic Liquids. Simple Liquid Metals. Molecular Liquids. Appendices. References. Index.

9,144 citations

Journal ArticleDOI
TL;DR: In this paper, it is shown that time averages of properties of the simulated fluid are equal to averages over the isoenthalpic-isobaric, canonical, and isothermal-isboric ensembles.
Abstract: In the molecular dynamics simulation method for fluids, the equations of motion for a collection of particles in a fixed volume are solved numerically. The energy, volume, and number of particles are constant for a particular simulation, and it is assumed that time averages of properties of the simulated fluid are equal to microcanonical ensemble averages of the same properties. In some situations, it is desirable to perform simulations of a fluid for particular values of temperature and/or pressure or under conditions in which the energy and volume of the fluid can fluctuate. This paper proposes and discusses three methods for performing molecular dynamics simulations under conditions of constant temperature and/or pressure, rather than constant energy and volume. For these three methods, it is shown that time averages of properties of the simulated fluid are equal to averages over the isoenthalpic–isobaric, canonical, and isothermal–isobaric ensembles. Each method is a way of describing the dynamics of a certain number of particles in a volume element of a fluid while taking into account the influence of surrounding particles in changing the energy and/or density of the simulated volume element. The influence of the surroundings is taken into account without introducing unwanted surface effects. Examples of situations where these methods may be useful are discussed.

4,998 citations