Abstract:A method to measure the phonon dispersion of a crystal based on molecular dynamics simulation is proposed and implemented as an extension to an open source classical molecular dynamics simulation code LAMMPS. In the proposed method, the dynamical matrix is constructed by observing the displacements of atoms during molecular dynamics simulation, making use of the fluctuation-dissipation theory. The dynamical matrix can then be employed to compute the phonon spectra by evaluating its eigenvalues. It is found that the proposed method is capable of yielding the phonon dispersion accurately, while taking into account the anharmonic effect on phonons simultaneously. The implementation is done in the style of fix of LAMMPS, which is designed to run in parallel and to exploit the functions provided by LAMMPS; the measured dynamical matrices could be passed to an auxiliary postprocessing code to evaluate the phonons.Program summaryProgram title: FixPhonon, version 1.0Catalogue identifier: AEJB_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEJB_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU General Public licenseNo. of lines in distributed program, including test data, etc.: 105 393No. of bytes in distributed program, including test data, etc.: 3 231 800Distribution format: tar.gzProgramming language: C++Computer: AllOperating system: LinuxHas the code been vectorized or parallelized?: Yes. 1 to N processors may be usedRAM: Depends on problem, approximate to 1 kB to several MBClassification: 7.8External routines: MPI, FIT, LAMMPS version 15, January 2010 (http://lammps.sandia.gov/)Nature of problem: Atoms in solids make ceaseless vibrations about their equilibrium positions, and a collective vibration forms a wave of allowed wavelength and amplitude. The quantum of such lattice vibration is called the phonon, and the so-called "lattice dynamics" is the field of study to find the normal modes of these vibrations. In other words, lattice dynamics examines the relationship between the frequencies of phonons and the wave vectors, i.e., the phonon dispersion. The evaluation of the phonon dispersion requires the construction of the dynamical matrix. In atomic scale modeling, the dynamical matrices are usually constructed by deriving the derivatives of the force field employed, which cannot account for the effect of temperature on phonons, with an exception of the tedious "quasi-harmonic" procedure.Solution method: We propose here a method to construct the dynamical matrix directly from molecular dynamics simulations, simply by observing the displacements of atoms in the system thus making the constructing of the dynamical matrix a straightforward task. Moreover, the anharmonic effect was taken into account in molecular dynamics simulations naturally, the resultant phonons therefore reflect the finite temperature effect simultaneously.Restrictions: A well defined lattice is necessary to employ the proposed method as well as the implemented code to evaluate the phonon dispersion. In other words, the system under study should be in solid state where atoms vibrate about their equilibrium positions. Besides, no drifting of the lattice is expected. The method is best suited for periodic systems, although non-periodic system with a supercell approach is also possible, it will however become inefficient when the unit cell contains too many atoms. Additional comments: The readers are encouraged to visit http://code.google.com/p/fix-phonon for subsequent update of the code as well as the associated postprocessing code, so as to keep up with the latest version of LAMMPS.Running time: Running time depends on the system size, the numbers of processors used, and the complexity of the force field, like a typical molecular dynamics simulation. For the third example shown in this paper, it took about 2.5 hours on an Intel Xeon X3220 architecture (2.4G, quadcore).

InterPhon: Ab initio interface phonon calculations within a 3D electronic structure framework

PHON: A program to calculate phonons using the small displacement method

API Phonons: Python Interfaces for Phonon Transport Modeling

Phq: A Fortran Code to Compute Phonon Quasiparticle Properties and Dispersions.

Phonon Dispersion Measured Directly from Molecular Dynamics Simulations

First-principles Phonon Calculations with Phonopy and Phono3py

Phoebe: a High-Performance Framework for Solving Phonon and Electron Boltzmann Transport Equations

A framework for solving atomistic phonon-structure scattering problems in the frequency domain using Perfectly Matched Layer boundaries

Ab initio calculation of electron-phonon linewidths and molecular dynamics with electronic friction at metal surfaces with numeric atom-centered orbitals

Exploring Phononic Properties of Two-Dimensional Materials using Machine Learning Interatomic Potentials

EPIq : an open-source software for the calculation of electron-phonon interaction related properties

First principles phonon calculations in materials science

Phonons from Density-Functional Perturbation Theory using the All-Electron Full-Potential Linearized Augmented Plane-Wave Method FLEUR

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

Accelerating High-Throughput Phonon Calculations via Machine Learning Universal Potentials

Perturbo: A software package for ab initio electron–phonon interactions, charge transport and ultrafast dynamics

Probing precise interatomic potentials by nonadiabatic nonlinear phonons

Panoramic Mapping of Phonon Transport from Ultrafast Electron Diffraction and Scientific Machine Learning

Primitive to conventional geometry projection for efficient phonon transport calculations

Interfacial Phonon Transport with Frequency-Dependent Transmissivity by Monte Carlo Simulation

Adaptively compressed polarizability operator for accelerating large scale \textit{ab initio} phonon calculations