-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathepw.txt
63 lines (57 loc) · 2.22 KB
/
epw.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
script python
=>save
2/ dossier epw a coté de phonon
3/ dans epw SNO2.in scf.in scf nscf
4/ SCF:
smearing = 'np' (petit => augmenter point k (666 ok) + degauss, tot_charge = 0.01
5/ nscf: shifter k (uniquement positif) 1/6 ... conv_tresh = 10d-12
nband = 28
6/ EPW: iverbosiry: 0 Tnk 3 Tnvnue (par phonon)/fishstick: borne enrgie outour fermi.
degausk = 0.1(en peu plus haut q'avant)
filkf: un point k gamma
nqf1: converngence la dessus
point q au bout (x cartesian puis points)
nbnd = 13
numiter 3000
wanier (frozen: point ou il doit passer)
linewidth
-----
Script crée dossier save et copie colle dossiers en les renommant pour être utilisés par EPW
0/ Phonon EPW
1/ Lancer script Python: Python pp.py
2/ Créer dossier PhononEPW (scf.in, ph.in) et EPW (+epw.in, scf.in, nscf.in)
- scf = même que phononEPW + occupation ‘smearing’/smearing = ‘mp’ = type de smearing, degauss (largeur de fermi-dirac) =1meV =0.001 petit plus précis /tot_charge =0.01
nband = 28 au moins pour inclure bande de conduction
- nscf = même que scf
Pour que k soit positifs (dans EPW) - shift zone de brilloin
Spécifier 6*6*6 = 216, specifier tous les points k 0 1/6 2/6
conv_thr = 1.0d-12
nband = 28 au moins pour inclure bande de conduction
3/ input EPW
dvscf_dir ../ph/save
iverbosity = 0 -> tau_nk
= 3 ->tau_nk\mu = donne en fonction des modes de phonons
elecselfen = .true. -> calcule sigma pour tau
phononselfen = .false.
a2f =
parallel k = true (q false)
fsthick = 1eV 😊 intervalle raisonnable de scattering considéré)
eptemp = 300K
degaussw = 0.02 (un peu plus) que pour scf (calcul moins précis mais plus smooth (?)
fermi_energy = (voir output nscf)
nk = 6
nq = 3
filukk = ‘’ (fichier séparé1 seul point k)
nqf1 = 3 puis convergence là-dessus 12*12*12 24 36 (points sur lesquels on prend l’interpolation)
à la fin
wannierize true
nombre de bandes wannierisees (13)
nombre skipper
niter pour minimiser extension spatiale wannier
dis-win_min = 0.8 (display bandes
dis_win_
dis_froz = window sur laquelle wannier passe par points calculés
proj = fonctions initiales données pour wannier, typiquement orbitales moléculaires hybridées
8 cartesian
Liste de points q donnée par phonon
line width