Electron-phonon interactions from Monte-Carlo sampling: Difference between revisions
No edit summary |
|||
| (21 intermediate revisions by 2 users not shown) | |||
| Line 1: | Line 1: | ||
'''Important''': This feature | '''Important''': This feature is only available from VASP6.0 or higher. | ||
First of all this method needs a sufficiently large super cell. | First of all this method needs a sufficiently large super cell. | ||
| Line 5: | Line 5: | ||
To enable electron-phonon interactions from MC methods {{TAG|PHON_LMC}}=''.TRUE.'' has to be set in the {{TAG|INCAR}} file. Also {{TAG|IBRION}}=6 has to be selected (the Monte-Carlo (MC) method is currently only implemented for {{TAG|IBRION}}=6). | To enable electron-phonon interactions from MC methods {{TAG|PHON_LMC}}=''.TRUE.'' has to be set in the {{TAG|INCAR}} file. Also {{TAG|IBRION}}=6 has to be selected (the Monte-Carlo (MC) method is currently only implemented for {{TAG|IBRION}}=6). | ||
The first implementation of electron-phonon interactions from MC sampling in VASP is found in Ref. {{cite|karsai:njp:2018}}. | |||
The original publication of the ZG configuration (one-shot method) is found in Ref. {{cite|zacharias:prb:2016}}. | |||
== Full MC sampling == | == Full MC sampling == | ||
| Line 29: | Line 33: | ||
The MC sampling code produces many {{TAG|POSCAR}} files with different distorted Wycoff positions but unchanged Brillouin matrix. The files are labeled as | The MC sampling code produces many {{TAG|POSCAR}} files with different distorted Wycoff positions but unchanged Brillouin matrix. The files are labeled as | ||
POSCAR.TEBEG.NUMBER | POSCAR.TEBEG.NUMBER | ||
where the NUMBER runs from 1 to {{TAGBL|PHON_NSTRUCT}}. | where the NUMBER runs from 1 to {{TAGBL|PHON_NSTRUCT}}. | ||
== ZG configuration (one-shot sampling) == | == ZG configuration (one-shot sampling) == | ||
Alternatively to the full MC sampling, a one shot method, introduced by M. Zacharias and F. Giustino{{cite|zacharias:prb:2016}} (named ZG configuration after the authors), is available. This method only uses a single distorted structure and hence it is orders of magnitude faster than the full MC sampling, while it retains an accuracy very close to the full MC sampling for converged super cell sizes. For example we showed that for the zero-point renormalization of the band gap the accuracy is within 5 meV between the ZG configurations and the full MC sampling. Hence we suggest to use this method preferably, where convergence of the super cell size is hard to achieve or the 5 meV accuracy is | Alternatively to the full MC sampling, a one shot method, introduced by M. Zacharias and F. Giustino{{cite|zacharias:prb:2016}} (named ZG configuration after the authors), is available. This method only uses a single distorted structure and hence it is orders of magnitude faster than the full MC sampling, while it retains an accuracy very close to the full MC sampling for converged super cell sizes. For example we showed that for the zero-point renormalization of the band gap the accuracy is within 5 meV between the ZG configurations and the full MC sampling{{cite|karsai:njp:2018}}. Hence we suggest to use this method preferably, where convergence of the super cell size is hard to achieve or the 5 meV accuracy is enough. | ||
=== Input === | === Input === | ||
| Line 40: | Line 43: | ||
To select the ZG configuration {{TAG|PHON_NSTRUCT}}=0 has to be set in the {{TAG|INCAR}} file. | To select the ZG configuration {{TAG|PHON_NSTRUCT}}=0 has to be set in the {{TAG|INCAR}} file. | ||
The temperatures (in K) have to be provided | The number of different temperatures and the list of temperatures (in K) have to be provided using the tags {{TAG|PHON_NTLIST}} and {{TAG|PHON_TLIST}}, respectively, in the {{TAG|INCAR}} file. An example would look like: | ||
PHON_TLIST = 0.0 100.0 200.0 . | {{TAGBL|PHON_NTLIST}} = 4 | ||
{{TAGBL|PHON_TLIST}} = 0.0 100.0 200.0 350.0 | |||
This makes the simultaneous calculation of the ZG configuration at several temperatures possible. | This makes the simultaneous calculation of the ZG configuration at several temperatures possible. | ||
An example {{TAG|INCAR}} file for a temperature range from 0-700 K (with step size of 100 K) is given as: | |||
{{TAGBL|System}} = DEFAULT | |||
{{TAGBL|PREC}} = Accurate | |||
{{TAGBL|ISMEAR}} = 0; {{TAGBL|SIGMA}} = 0.1; | |||
{{TAGBL|IBRION}} = 6 | |||
{{TAGBL|PHON_NTLIST}} = 8 | |||
{{TAGBL|PHON_TLIST}} = 0.0 100.0 200.0 300.0 400.0 500.0 600.0 700.0 | |||
{{TAGBL|PHON_NSTRUCT}} = 0 | |||
{{TAGBL|PHON_LMC}} = .TRUE. | |||
=== Output === | |||
Similar to the MC sampling the ZG configuration method produces several {{TAG|POSCAR}} files with different distorted Wycoff positions but unchanged Brillouin matrix. The files are labeled as | |||
POSCAR.TEMP | |||
where TEMP runs over all temperatures defined by {{TAG|PHON_TLIST}}. | |||
== Related tags and articles == | |||
{{TAG|PHON_LBOSE}}, | |||
{{TAG|PHON_LMC}}, | |||
{{TAG|PHON_NSTRUCT}}, | |||
{{TAG|PHON_NTLIST}}, | |||
{{TAG|PHON_TLIST}} | |||
== References == | == References == | ||
<references/> | <references/> | ||
---- | ---- | ||
[[Category: | [[Category:Phonons]][[Category:Electron-phonon interactions]][[Category:Howto]] | ||
Latest revision as of 14:43, 20 February 2024
Important: This feature is only available from VASP6.0 or higher.
First of all this method needs a sufficiently large super cell. It also involves phonon calculations for the point (see Phonons from finite differences). So many tags in the INCAR will be used from the phonon calculations.
To enable electron-phonon interactions from MC methods PHON_LMC=.TRUE. has to be set in the INCAR file. Also IBRION=6 has to be selected (the Monte-Carlo (MC) method is currently only implemented for IBRION=6).
The first implementation of electron-phonon interactions from MC sampling in VASP is found in Ref. [1].
The original publication of the ZG configuration (one-shot method) is found in Ref. [2].
Full MC sampling
Input
The tag PHON_NSTRUCT sets the number of structures generated due to the MC sampling. Convergence of the observable with respect to this number should be monitored.
The tag TEBEG=0 is also needed to choose the temperature at which the sampling is run.
Additionally the PHON_LBOSE can be set .TRUE. or .FALSE. (default PHON_LBOSE=.TRUE.), which selects Bose-Einstein or Maxwell-Boltzmann statistics, respectively.
A sample INCAR file for 0 K looks like the following:
System = DEFAULT PREC = Accurate ISMEAR = 0; SIGMA = 0.1; IBRION = 6 PHON_LMC = .TRUE. PHON_NSTRUCT = 100 TEBEG = 0.0
Output
The MC sampling code produces many POSCAR files with different distorted Wycoff positions but unchanged Brillouin matrix. The files are labeled as
POSCAR.TEBEG.NUMBER
where the NUMBER runs from 1 to PHON_NSTRUCT.
ZG configuration (one-shot sampling)
Alternatively to the full MC sampling, a one shot method, introduced by M. Zacharias and F. Giustino[2] (named ZG configuration after the authors), is available. This method only uses a single distorted structure and hence it is orders of magnitude faster than the full MC sampling, while it retains an accuracy very close to the full MC sampling for converged super cell sizes. For example we showed that for the zero-point renormalization of the band gap the accuracy is within 5 meV between the ZG configurations and the full MC sampling[1]. Hence we suggest to use this method preferably, where convergence of the super cell size is hard to achieve or the 5 meV accuracy is enough.
Input
To select the ZG configuration PHON_NSTRUCT=0 has to be set in the INCAR file.
The number of different temperatures and the list of temperatures (in K) have to be provided using the tags PHON_NTLIST and PHON_TLIST, respectively, in the INCAR file. An example would look like:
PHON_NTLIST = 4 PHON_TLIST = 0.0 100.0 200.0 350.0
This makes the simultaneous calculation of the ZG configuration at several temperatures possible.
An example INCAR file for a temperature range from 0-700 K (with step size of 100 K) is given as:
System = DEFAULT PREC = Accurate ISMEAR = 0; SIGMA = 0.1; IBRION = 6 PHON_NTLIST = 8 PHON_TLIST = 0.0 100.0 200.0 300.0 400.0 500.0 600.0 700.0 PHON_NSTRUCT = 0 PHON_LMC = .TRUE.
Output
Similar to the MC sampling the ZG configuration method produces several POSCAR files with different distorted Wycoff positions but unchanged Brillouin matrix. The files are labeled as
POSCAR.TEMP
where TEMP runs over all temperatures defined by PHON_TLIST.
Related tags and articles
PHON_LBOSE, PHON_LMC, PHON_NSTRUCT, PHON_NTLIST, PHON_TLIST