Hot answers tagged

29 votes

How big should a supercell be in phonon calculations?

Quick Summary: There's no way around performing a convergence test. However, it is possible to obtain convergence much faster than the Phonopy approach by using nondiagonal supercells [1]. The basic ...
user avatar
  • 10k
22 votes
Accepted

What do negative phonon frequencies signify?

Phonons are a measure of the curvature of the potential energy surface about a stationary point. In particular, the matrix of force constants is calculated as: $$ D_{i\alpha,i^{\prime}\alpha^{\prime}}(...
user avatar
  • 10k
18 votes
Accepted

Are there differences in accuracy and reliability between the frozen phonon method and Density Functional Perturbation Theory?

Short answer: Modern implementations of these two methods lead to similar accuracies. Longer answer: The calculation of phonons requires the calculation of the Hessian of the potential energy surface ...
user avatar
  • 10k
15 votes

How big should a supercell be in phonon calculations?

Ideally, a convergence test would be the best way to decide the required size of the supercell, but it can get expensive. When phonopy (or any similar calculation technique) finds displacements in the ...
user avatar
13 votes

How big should a supercell be in phonon calculations?

Considering this in terms of a 2x2x2 supercell matrix is the wrong way to think about this, as choice depends on the cell length and bonding type. Given that rigorous convergence testing is near-...
user avatar
12 votes
Accepted

What do phonon dispersion (or lattice dynamics) studies include?

TLDR: When you calculate phonons, you can describe electrons at different levels of theory, typically semilocal DFT, but also hybrids or dynamical mean-field theory. Phonons do include zero-point ...
user avatar
  • 10k
10 votes
Accepted

Reliability about the predicted two-dimensional crystal structures based on DFT high-throughput calculation?

[Disclaimer - I am one of the co-authors of the 2D database on Materials Cloud (What you call "2D structures and layered materials", publishing the data of this work: N. Mounet et al., ...
user avatar
9 votes
Accepted

Why do we need to build a supercell to use the DFPT in Phonopy interfaced with VASP?

Disclaimer: I have never used Phonopy. The advantage of using DFPT is that in principle it can be used to calculate a perturbation of finite wave vector $\mathbf{q}$ using the primitive cell. This ...
user avatar
  • 10k
7 votes
Accepted

How to determine LA/TA, ZA/ZO and LO/TO by PHONOPY

There are two major techniques used to model phonon interactions - The frozen phonon method and density functional perturbation theory (DFPT). Phonopy is used to carry out calculations in the frozen ...
user avatar
  • 2,172
6 votes

"Input forces are not enough to calculate force constants" Phonopy error?

I have never used Phonopy, so we will have to wait for an actual user for a full answer. However, I can make some general statements about phonon calculations that could plausibly explain your error. ...
user avatar
  • 10k
6 votes

Can we do anything to reduce the number of displacements in a frozen phonon lattice dynamics study?

ProfM's answer gets the core idea perfectly right: Symmetry really is your best friend here. However, symmetry analysis is often quite involved, especially for larger unit cells. I recently ...
user avatar
6 votes

Can we do anything to reduce the number of displacements in a frozen phonon lattice dynamics study?

You can (should) use symmetry to reduce the number of displacements needed to construct the matrix of force constants. A nice pratical description of how to do this can be found in the description of ...
user avatar
  • 10k
6 votes

Calculation of linear coefficient of thermal expansion using VASP-Phonopy for Hexagonal material

You can minimize the free energy F(a,c) as a function of (a,c), which is simply high school math. An example is provided in the section III (A) of the Supplementary Material of Phys. Rev. B 100, ...
user avatar
  • 675
5 votes

Temperature effect on elastic constant using VASP

It is always better to write your own code to get elastic constants. For example , In case of cubic system, we need three types of distorsion to unit cell. Now elastic constant is simply slope of ...
user avatar
  • 3,054
5 votes

Tools to Visualize Phonons

1. Interactive Phonon visualizer by Materials Cloud https://interactivephonon.materialscloud.io/ This tool allows to visualize the phonon dispersion calculated with an external software and display ...
user avatar
  • 8,362
5 votes
Accepted

How to solve this Phonopy-related problem?

This is basic bash-script which will be quite useful, just start with scratch ...
user avatar
  • 3,054
5 votes
Accepted

Which IBRION tag should be used for finite displacement method in calculating Phonons spectrum?

Phonopy creates different POSCAR file in which atoms are perturbed, as you are running each ...
user avatar
  • 3,054
5 votes

Phonon density of states in statistical thermodynamics

There are two possible things tripping you up here: Phonons are collective oscillations: they involve the motion of all the atoms together. Therefore it only makes sense to talk about the phonons of ...
user avatar
4 votes

Best practice(s) for q-points convergence in phonon density of states calculations

In a calculation of the phonon density of states, $\mathbf{q}$-points feature in two ways: Explicitly calculated $\mathbf{q}$-points. These are the $\mathbf{q}$-points for which you explicitly ...
user avatar
  • 10k
4 votes

Should all q-points in the first Brillouin zone be sampled in phonon density of states calculations?

Background. The phonon density of states $g$ is given by: $$ \tag{1} g(\omega)=\sum_{\nu}\int\frac{d\mathbf{q}}{(2\pi)^3}\delta(\omega-\omega_{\mathbf{q}\nu})\approx\frac{1}{N_{\mathbf{q}}}\sum_{\nu}\...
user avatar
  • 10k
4 votes

"Input forces are not enough to calculate force constants" Phonopy error?

This error (input forces are not enough to calculate force constants) may arise when you are trying to plot the band structure or the density of state plots for some systems. Sometimes setting a ...
user avatar
3 votes

Lattice Thermal conductivity calculation error

Partial Answers To remove this error Error in routine cdiaghg (211): S matrix not positive definite Use PAW pseudopotential files instead of USPP. Make sure the exchange correlational functional are ...
user avatar
3 votes

Calculating lattice thermal conductivity using Quantum ESPRESSO and phonopy

Some years ago I found the following steps in ResearchGate: You need: QuantumEspresso (QE) Phono3py ShengBTE & thirdorder.py First, you need to create a QE input file (...
user avatar
  • 18.5k
3 votes

Can we translate vibrational frequencies at q-points to vibrational frequencies of atoms?

As known [1], the eigenvalues of the dynamical Hessian matrix represent the phonon frequencies, whereas the eigenvectors represent the particular atomic displacement patterns contributing to the ...
user avatar
2 votes

Phonon frequencies and free energy with respect to polymorphs

Your interpretation of the results. I agree with you that if you find no imaginary frequencies in the cubic phase it means it is at a local minimum of the potential energy landscape, and that if you ...
user avatar
  • 10k
1 vote

Temperature effect on elastic constant using VASP

The following paper is quite nice, in my opinion: Shang et al., Temperature-dependent elastic stiffness constants of α- and θ-Al2O3 from first-principles calculations, J. Phys. Condens. Matter 22 (...
user avatar

Only top scored, non community-wiki answers of a minimum length are eligible