Troubleshooting electronic convergence: Difference between revisions
No edit summary |
No edit summary |
||
Line 3: | Line 3: | ||
Try to create a minimal [[INCAR]] file with as few tags as possible. | Try to create a minimal [[INCAR]] file with as few tags as possible. | ||
If the calculation converges then gradually add them back until you find which one was causing the | If the calculation converges then gradually add them back until you find which one was causing the problem. | ||
Check if you have enough bands. You can do this by looking at the [[OUTCAR]] file and checking that there are enough 'empty' states (i.e. states with zero occupation). When using an iterative solver the last states might not be accurately described, if these are occupied them convergence is | Check if you have enough bands. You can do this by looking at the [[OUTCAR]] file and checking that there are enough 'empty' states (i.e. states with zero occupation). When using an iterative solver the last states might not be accurately described, if these are occupied them convergence is likely to fail. | ||
Systems with f-orbitals are in particular difficult to converge. | Systems with f-orbitals are in particular difficult to converge. |
Revision as of 16:28, 23 November 2021
There can be many reasons why the SCF convergence is not reached. Naturally, there are also many different possible solutions.
Try to create a minimal INCAR file with as few tags as possible. If the calculation converges then gradually add them back until you find which one was causing the problem.
Check if you have enough bands. You can do this by looking at the OUTCAR file and checking that there are enough 'empty' states (i.e. states with zero occupation). When using an iterative solver the last states might not be accurately described, if these are occupied them convergence is likely to fail.
Systems with f-orbitals are in particular difficult to converge. Here are a couple of tips for a magnetic calculation with LDA+U:
- give initial magnetization only to the magnetic atoms
- use spin-polarized calculation
- perform the calculation in 3 steps (always starting from the previous WAVECAR):