Pushing convergence in parallel: Difference between revisions
(First version) |
No edit summary |
||
| Line 5: | Line 5: | ||
The system we are working on is the 2D-hBN sheet we already used the first day. | The system we are working on is the 2D-hBN sheet we already used the first day. | ||
[[File:HBN2.png| | [[File:HBN2.png|500px|center|Atomic structure of 2D hBN]] | ||
If you wan a published reference of some properties of hBN you can look to the following arxiv, | If you wan a published reference of some properties of hBN you can look to the following arxiv, | ||
| Line 13: | Line 13: | ||
interaction entering the GW self-energy. Here you can see an image extracted from the arxiv. | 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. | As you can see the QP corrections converge slowly with the number of k-points. | ||
[[File:QP convergence hBN from literature.png| | [[File:QP convergence hBN from literature.png|400px|center|Yambo tutorial image]] | ||
In the present tutorial you will start using a 12x12 grid. | 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 | Only in the last part we will re-run pw.x to generate bigger k-grids | ||
| Line 24: | Line 24: | ||
will be able to push much more the convergence. In a real life calculation you will likely take advantage | will be able to push much more the convergence. In a real life calculation you will likely take advantage | ||
of both. | 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) == | |||
<br> | <br> | ||
Revision as of 22:35, 25 April 2017
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 |