VCA: Difference between revisions
No edit summary |
No edit summary |
||
| Line 3: | Line 3: | ||
Description: {{TAG|VCA}} is short for the virtual crystal approximation; the tag allows to "weight" each species found in the POTCAR file. | Description: {{TAG|VCA}} is short for the virtual crystal approximation; the tag allows to "weight" each species found in the POTCAR file. | ||
---- | ---- | ||
The tag {{TAG|VCA}} | The tag {{TAG|VCA}} must be supplied for each atom type or species found in the {{TAG|POTCAR}} and {{TAG|POSCAR}} file, respectively . | ||
It weights the corresponding POTCAR files according to the values given in the {{TAG|INCAR}} file, with the default obviously | It weights the corresponding POTCAR files according to the values given in the {{TAG|INCAR}} file, with the default obviously | ||
being 1. For instance the formal valency found in the POTCAR files is multiplied by the supplied values, and likewise | being 1. For instance the formal valency found in the POTCAR files is multiplied by the supplied values, and likewise the local potential, | ||
the augmentation charges, and the non-local pseudopotential | the augmentation charges, and the non-local pseudopotential strength parameters are scaled by the supplied values. | ||
It is possible to use this flag to perform calculations in the framework of the virtual crystal approximation. | It is possible to use this flag to perform calculations in the framework of the virtual crystal approximation. | ||
| Line 25: | Line 25: | ||
0.25 0.25 0.25 | 0.25 0.25 0.25 | ||
If {{TAG| | If {{TAG|VCA}} is set to | ||
VCA = 0.99 0.01 | VCA = 0.99 0.01 | ||
the Ge atoms are weighted with a weight of 0.99, whereas the Sn atoms are weighted by 0.01 (see <ref name="eckhardt2014indirect"/>). | the Ge atoms are weighted with a weight of 0.99, whereas the Sn atoms are weighted by 0.01 (see <ref name="eckhardt2014indirect"/> for an example application). | ||
The implementation in VASP closely follows the | The implementation in VASP closely follows the methodology suggested by Bellaiche and Vanderbilt.<ref name="bellaiche2000VCA"/> | ||
Caveats: Unfortunately results of | Caveats: Unfortunately results of this kind of VCA calculations are often not very reliable. The problems are even apparent in the original publications by <ref name="bellaiche2000VCA"/>. The | ||
key point is that the used PAW potentials need to be constructed in such a | key point is that the used PAW potentials need to be constructed in such a way that the pseudo atomic waves are very similar for the potentials | ||
that are "mixed" (in the example above this would be the Ge and Sn potentials). This can be achieved by carefully optimizing the radial cutoffs. Furthermore, | that are "mixed" (in the example above this would be the Ge and Sn potentials). This can be achieved by carefully optimizing the radial cutoffs. Furthermore, | ||
the local | the local potentials of the two POTCAR files need to be very similar. This means that results for many standard potentials are not accurate. For instance, Vegard's law is | ||
often not even approximately observed (instead the volume is way too large at 50 %). The problem is particularly severe, if semi-core states | often not even approximately observed (instead the volume is way too large at 50 % mixing). The problem is particularly severe, if semi-core states | ||
are treated as valence states. For instance, for Ge and Sn the d electrons need to be treated as core electrons to obtain reasonable results. | are treated as valence states. For instance, for Ge and Sn the d electrons need to be treated as core electrons to obtain reasonable results. | ||
Revision as of 12:03, 16 April 2020
VCA = [real array]
Default: VCA = read from the POTCAR file
Description: VCA is short for the virtual crystal approximation; the tag allows to "weight" each species found in the POTCAR file.
The tag VCA must be supplied for each atom type or species found in the POTCAR and POSCAR file, respectively . It weights the corresponding POTCAR files according to the values given in the INCAR file, with the default obviously being 1. For instance the formal valency found in the POTCAR files is multiplied by the supplied values, and likewise the local potential, the augmentation charges, and the non-local pseudopotential strength parameters are scaled by the supplied values.
It is possible to use this flag to perform calculations in the framework of the virtual crystal approximation. Say you want to simulate Sn doping in a Ge lattice. Using a POTCAR file with a Ge and Sn data set and the following POSCAR file, this can be achieved:
cd:
1.00000000000000
2.82173 2.82173 0.00000
0.00000 2.82173 2.82173
2.82173 0.00000 2.82173
Ge Sn
2 2
Direct
0.00 0.00 0.00
0.25 0.25 0.25
0.00 0.00 0.00
0.25 0.25 0.25
If VCA is set to
VCA = 0.99 0.01
the Ge atoms are weighted with a weight of 0.99, whereas the Sn atoms are weighted by 0.01 (see [1] for an example application). The implementation in VASP closely follows the methodology suggested by Bellaiche and Vanderbilt.[2]
Caveats: Unfortunately results of this kind of VCA calculations are often not very reliable. The problems are even apparent in the original publications by [2]. The key point is that the used PAW potentials need to be constructed in such a way that the pseudo atomic waves are very similar for the potentials that are "mixed" (in the example above this would be the Ge and Sn potentials). This can be achieved by carefully optimizing the radial cutoffs. Furthermore, the local potentials of the two POTCAR files need to be very similar. This means that results for many standard potentials are not accurate. For instance, Vegard's law is often not even approximately observed (instead the volume is way too large at 50 % mixing). The problem is particularly severe, if semi-core states are treated as valence states. For instance, for Ge and Sn the d electrons need to be treated as core electrons to obtain reasonable results.