33
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 ...
- 10.7k
24
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}}(...
- 10.7k
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 ...
- 10.7k
16
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 ...
- 5,065
14
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-...
- 526
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 ...
- 10.7k
12
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., ...
- 216
10
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 ...
- 10.7k
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 ...
- 2,232
7
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. ...
- 10.7k
7
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 ...
- 1,416
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 ...
- 10.7k
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, ...
- 675
6
votes
Accepted
How to solve this Phonopy-related problem?
This is basic bash-script which will be quite useful, just start with scratch
...
- 3,834
6
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 ...
- 3,834
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 ...
- 8,612
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 ...
- 3,834
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 ...
- 6,469
5
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 ...
- 10.7k
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}\...
- 10.7k
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 ...
- 846
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 ...
- 1,015
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 (...
- 20.7k
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 ...
- 336
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 ...
- 10.7k
2
votes
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 (...
- 1,708
1
vote
How to calculate the lattice parameter of a supercell using Phonopy? (For Quantum ESPRESSO)
The process you described for obtaining a supercell can be found here.
If the 1*1*1 unit cell has a lattice parameter of a, then a 2*2*2 supercell will have a lattice parameter of 2a, a 3*3*3 ...
- 457
1
vote
k-path for phonon band structure
The question was:
"but why do some articles use different path"
However, no article was provided that uses a different path. ProfM said:
"I agree there is no reason why the two paths ...
1
vote
Phonon branch-resolved harmonic sampling
OP:
"Hello! Perhaps, the most general formalism, which would allow to perform such analysis is based on the Berry phase theory. There are a bunch of papers doing this, for example, journals.aps....
Only top scored, non community-wiki answers of a minimum length are eligible