Phonon Dispersions and Density of States Calculation on Grid¶
This tutorial page explains how to calculate the Phonon Dispersion Curves and Phonon Density of States of materials based on Density Functional Theory. We will be studying crystalline Silicon in the standard cubicdiamond crystal structure, and we will use Quantum ESPRESSO as our simulation engine.
Quantum ESPRESSO version considered in this tutorial
The present tutorial is written for Quantum ESPRESSO at versions 5.2.1, 5.4.0, 6.0.0 or 6.3.
What sets the present tutorial apart from the other tutorial on phonon calculations is the employment of the "Grid Method" for computing the vibrational properties of materials, which is reviewed in the subsequent paragraph. This method is based on a map type workflow, where multiple branches are executed in parallel as separate independent Jobs, with the consequent gain in computational efficiency and overall speed of the phonon calculation. More information on this method, together with a demonstration of its application and results on a sample set of materials, can be found in Ref. ^{1}.
The Grid Method for Phonon Calculations¶
The Grid Method allows for an efficient parallelization of the tasks for calculating the individual vibrational modes. This method optimizes the corresponding workflow in order to obtain the frequencies for each individual symmetryirreducible representation ^{2} of the phonon lattice perturbations in parallel.
Steps Involved in Grid Method¶
We thus implement a gridparallel workflow for the calculation of the phonon dynamical matrices, initially explained in the first reference cited in this page. During the actual phonon calculation part of the workflow, the following steps happen:

First, the irreducible representations for the vibrational modes (irreps) are generated, based on the sampling grid in the reciprocal space (qpoint grid).

Second, a separate calculation is prepared and submitted for execution in parallel for each pair or irreducible representation and qpoint, via the implementation of a Map for performing a distributed calculation ("map" stage).

Finally, after the calculations for all irrepsq points pairs are finished, the dynamical matrices are collected and aggregated together, and phonon dispersions and density of states are calculated ("reduce" stage).
Thus, we employ a "mapreduce" general type of logic and scenario within the overall workflow, where the individual calculation tasks are performed independently and in parallel with one another. This allows for an improved efficiency and speedup of phonon calculations compared to the more traditional serial approach, such that the limiting factor within the overall calculation is the longest run per individual irreducible representationq point pair.
Schematic Visualization of Grid Method¶
A schematic summary of the above workflow procedure is offered in the figure below. This flowchart depicts all different approaches to the calculation of the phonon dispersions and lattice vibrations. The approach used in the present tutorial is the rightmost one. Here, ”SCF” stands for the self consistent field preliminary calculation, ”ph.x” denotes the phonon calculations by means of Density Functional Perturbation Theory, and "irrep" is an irreducible representation of a vibrational phonon mode.
Workflow Structure¶
We review now the different steps involved in implementing the above general Grid Phonon theoretical framework in an actual Workflow deployable on our platform. We shall consider the example case of the Quantum ESPRESSO modeling engine.
1. Preliminary SCF Calculation¶
The first subworkflow step in the overall "Phonons Grid" Workflow is a standard selfconsistent field (scf) total energy calculation, providing the ensuing steps of the workflow with the wavefunctions of the material structure under investigation.
For the sake of this example, we can set the grid of special kpoints to 6 x 6 x 6, under Important Settings.
2. Qpoints and Irrep Generation¶
The second subworkflow step ("phinitqpoints") is composed of a single unit, which consists in generating the grid of qpoints over which the vibrational phonon modes calculations will be performed, for each irreducible representation of such modes.
The size of this qgrid should be a divisor of the size of the abovementioned kgrid. Hence, a qgrid of 2 x 2 x 2 should suffice in this case, which can be set under Important Settings.
3. Extraction of qpoints/irrep pairs¶
The "espressoxmlgetqptirr" subworkflow comprises a main python script, whose role is to parse and extract qpoints and irreducible representations from the previouslygenerated Quantum ESPRESSO XML data. Each distinct combination of qpoint and irrep is added as a separate entry to what follows.
4. Map Distributed Phonon Calculation¶
The ensuing subworkflow consists in the Map for performing the distributed parallel phonon calculations across each qpoint/irrep pair extracted previously. The list of qpoints/irrep pairs can be inspected under the "Data" tab of the Map editor interface.
Care should be taken to set the qgrid under the "Important Settings" of the "phsingleirrqpt" map subworkflow again to 2 x 2 x 2, as in the previous steps.
5. Reduce and Aggregate Results¶
The final "Reduce" subworkflow collects together the results of the previous calculations over each independent qpoint/irrep pair, via the "ph_grid_restart" unit. Here, the size of the qgrid under Important Settings should once again be set to the 2 x 2 x 2 value being considered in the present example.
These combined results are then used to complete the phonon dispersion and density of states calculation, through the help of the Quantum ESPRESSO "q2r" and "matdyn" executables.
Creating and Executing Job¶
"Phonon Map" workflows can readily be imported into the accountowned collection from the Workflows Bank.
Apart from this, the same procedural instructions as in the other phonons calculation tutorial should be followed for creating and launching the corresponding gridbased phonon Job through our Web Interface, and for inspecting the associated results.
Animation¶
In the video animation below, we outline the procedure for creating and executing a phonon calculation job via the Grid Method. We conclude by inspecting the corresponding results for the Phonon Dispersion Curves and Density of States, considering crystalline silicon as our demonstrative sample operated in conjunction with the Quantum ESPRESSO simulation engine.
Computational cost of phonon calculation
Phonon calculations are in general quite computationally demanding. We therefore recommend the employment of at least 8 computing cores. For larger calculations, OF queues will have faster turnaround than the OR queues considered in the video.