About small displacement method for phonons

Click For Summary
SUMMARY

The small displacement method for calculating phonons requires that the wave vector k be orthogonal to the supercell size, meaning not all k vectors are permissible. Each choice of k vector necessitates a full electronic calculation, and the supercell must accommodate the periodicity of the k vector. Additionally, a sufficiently large supercell is essential to minimize the Hellmann-Feynman force outside the supercell, as the dynamical matrix involves an infinite summation of force constants.

PREREQUISITES
  • Understanding of phonon calculations in solid-state physics
  • Familiarity with the small displacement method
  • Knowledge of supercell construction in computational materials science
  • Experience with electronic structure calculations
NEXT STEPS
  • Research the implementation of the small displacement method in Quantum ESPRESSO
  • Learn about the construction and optimization of supercells in VASP
  • Study the Hellmann-Feynman theorem and its implications in phonon calculations
  • Explore the use of k-point sampling in density functional theory (DFT) calculations
USEFUL FOR

Researchers and practitioners in computational materials science, particularly those focused on phonon calculations and electronic structure methods.

Hyla Brook
Messages
22
Reaction score
0
Dear All,

I have trouble to understand the calculation of phonon using small displacement method. I found people said the limitation of this method was that it requires the wave vector k orthogonal to the supercell size(?). Does it mean the wave vector is not any you want? why? I got really confused. Any idea is highly appreciated. Thank you!

H.B
 
Physics news on Phys.org
Yes. You have to do a full electronic calculation for each choice of the k vector of the displacement. The super-cell has to fit the period of the k-vector.
 
DrDu said:
Yes. You have to do a full electronic calculation for each choice of the k vector of the displacement. The super-cell has to fit the period of the k-vector.

Thank you for reply Dr. Du. But in the program (e.g. phon), the number of k point is variable. That means we can choose any k wave vector we want in the BZ.

The tutorial of this program says the limitation of this method is we need a large supercell to make Hellmann-Feynman force negligible outside the supercell. I think this is easy to understand because the element of the dynamical matrix is a summation involving the force constants, which in principle has infinite terms. I asked that question because that explanation was mentioned as equivalent as this one. I cannot even understand that well.

Could you explain it more?
 
Last edited:

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 11 ·
Replies
11
Views
4K
  • · Replies 1 ·
Replies
1
Views
5K
  • · Replies 5 ·
Replies
5
Views
508
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 6 ·
Replies
6
Views
4K