Pushing convergence in parallel
2D hBN
In the previous session we focused on the scaling properties of the code. We now focus more on the physical aspect.
The system we are working on is the 2D-hBN sheet we already used the first day.
If you wan a published reference of some properties of hBN you can look to the following arxiv, https://arxiv.org/abs/1511.00129 and references therein were the QP correction are computed. In that case the convergence with the number of k-points is studied in detail (for hBN and other materials). The focus is on a proper description of the q->0 limit of the (screened) electron-electron interaction entering the GW self-energy. Here you can see an image extracted from the arxiv. As you can see the QP corrections converge slowly with the number of k-points.
In the present tutorial you will start using a 12x12 grid. Only in the last part we will re-run pw.x to generate bigger k-grids
Notice that to run convergence you will need to manually change and re-run the input file. For the k-points case in particular you will need to re-run also pw.x Tomorrow you will learn how to do that in a more automatic way using the yambopy scripting tools, however running yambo in serial. While yambopy offers an easier interface for the use, in the present parallel tutorial you will be able to push much more the convergence. In a real life calculation you will likely take advantage of both.
Convergence on the number bands in X and G
Convergence with the energy cut-off on the response function
Convergence with the number of k-points (optional)
| Prev: Tutorials Home | Now: Tutorials Home --> GW Convergence | Next: Parallel tutorials are over |