This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
creating_inputs [2016/12/08 08:29] krejcio |
creating_inputs [2017/02/10 15:23] krejcio |
||
---|---|---|---|
Line 4: | Line 4: | ||
[[http://www.fireball-dft.org]] | [[http://www.fireball-dft.org]] | ||
- | Unfortunately official version cannot write out input files. | + | Official version of //FireballTG// can write out input files now. |
- | A working executable for creation of necessary input files is at: | + | A working executable for creation of necessary input files is also at: |
/storage/praha1/home/krejcio/bin_fireball_stable/fireball.x | /storage/praha1/home/krejcio/bin_fireball_stable/fireball.x | ||
Input files can computed with McWEDA functional as well as with computing XC on a grid. Cluster systems as well as systems with Periodic Boundary Conditions (PBC) can be computed. | Input files can computed with McWEDA functional as well as with computing XC on a grid. Cluster systems as well as systems with Periodic Boundary Conditions (PBC) can be computed. | ||
- | A fireball.in file for calculations with McWEDA: | + | A __fireball.in__ file for calculations with McWEDA: |
&OPTION | &OPTION | ||
Line 28: | Line 28: | ||
&END | &END | ||
- | A fireball.in file for calculations with XC on a grid computations: | + | A __fireball.in__ file for calculations with XC on a grid computations: |
&OPTION | &OPTION | ||
Line 51: | Line 51: | ||
&END | &END | ||
- | In case of PBC calculations phik_0001_s.dat, phik_0001_py.dat, ... files are produced by the Fireball. In case of cluster calculations phik_s.dat, phik_py.dat, ... are outputs of the Fireball calculations. They serve as inputs for the PP-STM calculations. Inside they look like: | + | In case of PBC calculations __phik_0001_s.dat__, __phik_0001_py.dat__, ... files are produced by the Fireball. In case of cluster calculations __phik_s.dat__, __phik_py.dat__, ... are outputs of the //Fireball// calculations. They serve as inputs for the PP-STM calculations. Inside they look like: |
38 280 -5.37896401 Number of atoms Number of states (Molecular orbitals) The Fermi Level | 38 280 -5.37896401 Number of atoms Number of states (Molecular orbitals) The Fermi Level | ||
Line 61: | Line 61: | ||
[[http://wiki.fysik.dtu.dk/gpaw/]] | [[http://wiki.fysik.dtu.dk/gpaw/]] | ||
- | Even though the GPAW is mainly used for representing the wave-function on a grid it can work in LCAO mode as well. For the purpose of making inputs for the PP-STM calculations the LCAO mode is necessary. Both - default or double-zeta (basis='dzp'; for more information look at the GPAW web page) - basis sets can be used. The PP-STM code reads the stored *.gpw binary produced by the GPAW calculations. Here is an example of some GPAW script for the calculations of the input: | + | Even though the //GPAW// is mainly used for representing the wave-function on a grid it can work in LCAO mode as well. For the purpose of making inputs for the PP-STM calculations the LCAO mode is necessary. Both - default or double-zeta (basis='dzp'; for more information look at the //GPAW// web page http://wiki.fysik.dtu.dk/gpaw/) - basis sets can be used. The PP-STM code reads the stored *.gpw binary produced by the //GPAW// calculations. Here is an example of some //GPAW// script for the calculations of the input: |
from ase import * | from ase import * | ||
Line 81: | Line 81: | ||
calc.write('out_LCAO_'+xc+'.gpw',mode='all') # saves the calculation into binary 'out_LCAO_LDA.gpw' file | calc.write('out_LCAO_'+xc+'.gpw',mode='all') # saves the calculation into binary 'out_LCAO_LDA.gpw' file | ||
- | The results of the GPAW calculations is stored in binary file 'out_LCAO_LDA.gpw' | + | The results of the //GPAW// calculations is stored in binary file '__out_LCAO_LDA.gpw__' |
Line 89: | Line 89: | ||
Works for PBC calculations, just add: | Works for PBC calculations, just add: | ||
output eigenvectors | output eigenvectors | ||
- | output band 0 0 0 0.5 0.5 0.0 3 G K | + | output band 0 0 0 0.5 0.5 0.0 2 G K |
- | into control.in. | + | into __control.in__. |
+ | In the case of hybrid functionals (B3LYP, PBE0, HSE) add also: | ||
+ | exx_band_structure_version 1 ## for hybrid functionals only ## | ||
+ | into your __control.in__ file. | ||
- | The calculations produces: KS_eigenvectors.band_1.kpt_1.out for spin-restricted calculations or KS_eigenvectors_dn.band_1.kpt_1.out & KS_eigenvectors_up.band_1.kpt_1.out in the case of spin-polarized calculations. | + | The calculations produces: __KS_eigenvectors.band_1.kpt_1.out__ for spin-restricted calculations or __KS_eigenvectors_dn.band_1.kpt_1.out__ & __KS_eigenvectors_up.band_1.kpt_1.out__ in the case of spin-polarized calculations. |
Note: For cluster calculations Mathematica scripts have to be used for creating PP-STM inputs. | Note: For cluster calculations Mathematica scripts have to be used for creating PP-STM inputs. | ||