[Pw_forum] PW+W90 spin polarized calculations
Gianluca Giovannetti
gianluca.giovannetti at gmail.com
Wed Feb 24 09:50:42 CET 2010
Dear All,
i`m performing PW+W90 spin polarized calculations.
I first do a scf run:
---------------------------------------------------------------------------------------
FeSe
&control
calculation = 'scf'
restart_mode = 'from_scratch'
prefix = 'fese'
tprnfor = .true.
pseudo_dir = './'
outdir = './'
iprint = 2
/
&system
ibrav = 8
celldm(1) = 10.0733642480265
celldm(2) = 1.00000000000000
celldm(3) = 1.0291721306383
nat = 8
ntyp = 5
ecutwfc = 25.0
ecutrho = 300.0
nbnd = 80
occupations = 'smearing'
smearing = 'methfessel-paxton'
degauss = 0.01
nspin = 2
starting_magnetization(1) = 1.0
starting_magnetization(2) = 1.0
starting_magnetization(3) = -1.0
starting_magnetization(4) = -1.0
starting_magnetization(5) = 0.0
lda_plus_u =.true.
Hubbard_U(1) = 1.d-8
Hubbard_U(2) = 1.d-8
Hubbard_U(3) = 1.d-8
Hubbard_U(4) = 1.d-8
Hubbard_alpha(1) = 1.d-8
Hubbard_alpha(2) = 1.d-8
Hubbard_alpha(3) = 1.d-8
Hubbard_alpha(4) = 1.d-8
/
&electrons
conv_thr = 1.0d-8
diagonalization = 'cg'
/
ATOMIC_SPECIES
Fe1 55.845 Fe.pbe-sp-van.UPF
Fe2 55.845 Fe.pbe-sp-van.UPF
Fe3 55.845 Fe.pbe-sp-van.UPF
Fe4 55.845 Fe.pbe-sp-van.UPF
Se 78.96 Se.pbe-van.UPF
ATOMIC_POSITIONS { angstrom }
Fe1 0.00000000 2.66529759 2.74305000
Fe2 0.00000000 0.00000000 2.74305000
Fe3 2.66529759 2.66529759 2.74305000
Fe4 2.66529759 0.00000000 2.74305000
Se 1.33264880 3.99794639 4.02515157
Se 3.99794639 3.99794639 1.46094843
Se 3.99794639 1.33264880 4.02515157
Se 1.33264880 1.33264880 1.46094843
K_POINTS {automatic}
12 12 8 0 0 0
---------------------------------------------------------------------------------------
then a nscf run:
---------------------------------------------------------------------------------------
FeSe
&control
calculation = 'nscf'
! restart_mode = 'from_scratch'
prefix = 'fese'
tprnfor = .true.
pseudo_dir = './'
outdir = './'
iprint = 2
wf_collect = .true.
/
&system
ibrav = 8
celldm(1) = 10.0733642480265
celldm(2) = 1.00000000000000
celldm(3) = 1.0291721306383
nat = 8
ntyp = 5
ecutwfc = 25.0
ecutrho = 300.0
nbnd = 80
occupations = 'smearing'
smearing = 'methfessel-paxton'
degauss = 0.01
nspin = 2
starting_magnetization(1) = 1.0
starting_magnetization(2) = 1.0
starting_magnetization(3) = -1.0
starting_magnetization(4) = -1.0
starting_magnetization(5) = 0.0
lda_plus_u =.true.
Hubbard_U(1) = 1.d-8
Hubbard_U(2) = 1.d-8
Hubbard_U(3) = 1.d-8
Hubbard_U(4) = 1.d-8
Hubbard_alpha(1) = 1.d-8
Hubbard_alpha(2) = 1.d-8
Hubbard_alpha(3) = 1.d-8
Hubbard_alpha(4) = 1.d-8
nosym=.true.
noinv=.true.
/
&electrons
conv_thr = 1.0d-8
diagonalization = 'cg'
/
ATOMIC_SPECIES
Fe1 55.845 Fe.pbe-sp-van.UPF
Fe2 55.845 Fe.pbe-sp-van.UPF
Fe3 55.845 Fe.pbe-sp-van.UPF
Fe4 55.845 Fe.pbe-sp-van.UPF
Se 78.96 Se.pbe-van.UPF
ATOMIC_POSITIONS { angstrom }
Fe1 0.00000000 2.66529759 2.74305000
Fe2 0.00000000 0.00000000 2.74305000
Fe3 2.66529759 2.66529759 2.74305000
Fe4 2.66529759 0.00000000 2.74305000
Se 1.33264880 3.99794639 4.02515157
Se 3.99794639 3.99794639 1.46094843
Se 3.99794639 1.33264880 4.02515157
Se 1.33264880 1.33264880 1.46094843
K_POINTS crystal
1152
[...here i removed the list....]
---------------------------------------------------------------------------------------
then i start with Wannier90.
Let look only at the spin channel "up".
I run at first:
/work/aha/giovanne/ESPRESSO/espresso-4.0.5/W90/wannier90.x -pp FeSe_up
the FeSe_up input file looks like:
---------------------------------------------------------------------------------------
num_wann = 20
num_iter = 5000
dis_num_iter = 500
iprint = 3
exclude_bands : 1-32,53-80
write_r2mn=.true.
hr_plot=.true.
kmesh_tol=0.00001
spin = up
begin atoms_cart
ang
Fe 0.00000000 2.66529759 2.74305000
Fe 0.00000000 0.00000000 2.74305000
Fe 2.66529759 2.66529759 2.74305000
Fe 2.66529759 0.00000000 2.74305000
Se 1.33264880 3.99794639 4.02515157
Se 3.99794639 3.99794639 1.46094843
Se 3.99794639 1.33264880 4.02515157
Se 1.33264880 1.33264880 1.46094843
end atoms_cart
begin kpoint_path
G 0.0000 0.0000 0.000 X 0.5000 0.0000 0.0000
X 0.5000 0.0000 0.000 M 0.5000 0.5000 0.000
M 0.5000 0.5000 0.000 Y 0.0000 0.5000 0.0000
Y 0.0000 0.5000 0.000 G 0.0000 0.0000 0.0000
G 0.0000 0.0000 0.000 Z 0.0000 0.0000 0.5000
end kpoint_path
bands_plot =T
!fermi_surface_plot = .true.
!To plot the WF
!wannier_plot = T
!wannier_plot_supercell = 2
!wannier_plot_list = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
!restart = plot
guiding_centres = T
begin projections
Fe:dxy;dxz;dyz;dx2-y2;dz2
end projections
begin unit_cell_cart
bohr
10.0733642480265 0.0000 0.0000
0.0000 10.0733642480265 0.0000
0.0000 0.0000 10.3672257458371
end unit_cell_cart
mp_grid : 12 12 8
begin kpoints
[....]
end kpoints
---------------------------------------------------------------------------------------
then:
mpirun -np 8 /work/aha/giovanne/ESPRESSO/espresso-4.0.5/bin/pw2wannier90.x <
FeSe_up.pw2wan
with FeSe_up.pw2wan:
---------------------------------------------------------------------------------------
&inputpp
outdir = './'
prefix = 'fese'
seedname = 'FeSe_up'
spin_component = 'up'
write_mmn = .true.
write_amn = .true.
write_unk = .true.
/
---------------------------------------------------------------------------------------
once i re-run W90 as:
/work/aha/giovanne/ESPRESSO/espresso-4.0.5/W90/wannier90.x FeSe_up
i get the following error:
---------------------------------------------------------------------------------------
Found a mismatch in FeSe_up.eig
Wanted band : 1 found band : 21
Wanted kpoint: 2 found kpoint: 1
A common cause of this error is using the wrong
number of bands. Check your input files.
If your pseudopotentials have shallow core states remember
to account for these electrons.
Exiting.......
param_read: mismatch in FeSe_up.eig
---------------------------------------------------------------------------------------
i think i should have 20 eigenvalues in FeSe_up.eig but they are 40.
Could you please help me to understand it?
thank you.
Gianluca
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20100224/cc91e261/attachment.html>
More information about the users
mailing list