[Pw_forum] Wrong effective charges in AFM insulator using "zue"

xirainbow nkxirainbow at gmail.com
Wed Jan 19 09:27:29 CET 2011


In 4.2.1 manual, it is said that:*If .zue=true.is an alternative
algorithm, different from the default one (if trans .and. epsil ) The
results should be the same within numerical noise.*

*However, I get wrong "effective charges" in zue=true in AFM insulator
CaMnO3 with **4.2.1*
***These two Mn atoms are the same except antiparallel spin.*
*In "epsil+trans", effective charges of Mn1=7.79, Mn2=7.78;*
*In "zue"             , effective charges of Mn1=3.1,   Mn2=51.2;*

I will describe my calculation procedure using *4.2.1* in short.
1:  I do scf calculation in scf.in
*restart_mode  = 'from_scratch'*
*occupations   = 'smearing'*
*smearing      = 'gauss'*
*degauss       = 0.002*
*starting_magnetization(1) = 1*
*starting_magnetization(2) = -1*
*starting_magnetization(3) = 0*
*starting_magnetization(4) = 0*

2: I do a scf calculation  in scf.in
*restart_mode  = 'restart'*
*occupations               = 'fixed'*
*nspin                     = 2*
*tot_magnetization         = 0*

3: I do ph calculation with
*
zue         = .TRUE. !the same with epsil and trans = TRUE
(In this condition, effective charges of Mn1=3.1,   Mn2=51.2;)
or with:
epsil       = .TRUE. , trans       = .TRUE.
*
(In this condition *, *effective charges of Mn1=7.79, Mn2=7.78;)

(The reason is why I restart scf is below:
"starting_magnetization" must be used with "smearing". However, "smearing"
will give a metal result. Effective charges can not be calculated for metal.
But, CaMnO3 is really a insulator AFM.
"occupations=fixed" can get a insulator result. However it must used with
"tot_magnetization". If I use "tot_magnetization=0", I will get a PM no AFM
state.
In the end, I use "smearing"+"starting_magnetization" to get a AFM metal.
Then I use "fixed"+"tot_magnetization" to get a AFM insulator.
*Is there any other good method to get a AFM insulator?*)


The following my script for "*zue= .TRUE *"
And I have attached output file of "*zue= .TRUE *" in accessary.  Thanks:)

#*********************************************
# scf-antimag-smearing
cat>scf.in<<EOF
for vc-relax
&CONTROL
calculation   = 'scf'
title         = 'G-10-scf'
verbosity     = 'high'
restart_mode  = 'from_scratch'
wf_collect    = .FALSE.
tstress       = .TRUE.
tprnfor       = .TRUE.
prefix        = 'G-10-scf'
etot_conv_thr = 1.0e-4
forc_conv_thr = 1.0e-3
disk_io       = 'low'
pseudo_dir    = './'
/

&SYSTEM
ibrav         = 2
celldm        = 14.20419765
nat           = 10
ntyp          = 4
nbnd          = 50
ecutwfc       = 70
ecutrho       = 560
nosym         = .FALSE.
!nosym_evc     = .TRUE.
!noinv         = .TRUE.
occupations   = 'smearing'
smearing      = 'gauss'
degauss       = 0.002
starting_magnetization(1) = 1
starting_magnetization(2) = -1
starting_magnetization(3) = 0
starting_magnetization(4) = 0
!occupations               = 'fixed'
nspin                     = 2
!tot_magnetization         = 0
/

&ELECTRONS
electron_maxstep  = 100
conv_thr          = 1.0e-7
mixing_mode       = 'plain'
mixing_beta       = 0.7
mixing_ndim       = 8
diagonalization   = 'david'
diago_david_ndim  = 4
/

ATOMIC_SPECIES
 Mn1 55  Mn.pbe-sp-van.UPF
 Mn2 55  Mn.pbe-sp-van.UPF
 Ca  40  Ca.pbe-nsp-van.UPF
 O  16  O.pbe-rrkjus.UPF

ATOMIC_POSITIONS {crystal}
 Mn1 0.00    0.00    0.00
 Mn2 0.50    0.50    0.50
 Ca  0.25    0.25    0.25
 Ca  0.75    0.75    0.75
 O   0.75    0.25    0.75
 O   0.25    0.75    0.25
 O   0.75    0.75    0.25
 O   0.25    0.25    0.75
 O   0.25    0.75    0.75
 O   0.75    0.25    0.25


K_POINTS {automatic}
6 6 6 0 0 0
EOF
./dopw <scf.in >G-scf-10-smearing.out

# scf-antimag-insulator
cat>scf.in<<EOF
for vc-relax
&CONTROL
calculation   = 'scf'
title         = 'G-10-scf'
verbosity     = 'high'
restart_mode  = 'restart'
wf_collect    = .FALSE.
tstress       = .TRUE.
tprnfor       = .TRUE.
prefix        = 'G-10-scf'
etot_conv_thr = 1.0e-4
forc_conv_thr = 1.0e-3
disk_io       = 'low'
pseudo_dir    = './'
/

&SYSTEM
ibrav         = 2
celldm        = ${lattic_parameter} !bohr
nat           = 10
ntyp          = 4
nbnd          = 50
ecutwfc       = 70
ecutrho       = 560
nosym         = .FALSE.
!nosym_evc     = .TRUE.
!noinv         = .TRUE.
!occupations   = 'smearing'
!smearing      = 'gauss'
!degauss       = 0.002
!starting_magnetization(1) = 1
!starting_magnetization(2) = -1
!starting_magnetization(3) = 0
!starting_magnetization(4) = 0
occupations               = 'fixed'
nspin                     = 2
tot_magnetization         = 0
/

&ELECTRONS
electron_maxstep  = 100
conv_thr          = 1.0e-7
mixing_mode       = 'plain'
mixing_beta       = 0.7
mixing_ndim       = 8
diagonalization   = 'david'
diago_david_ndim  = 4
/

ATOMIC_SPECIES
 Mn1 55  Mn.pbe-sp-van.UPF
 Mn2 55  Mn.pbe-sp-van.UPF
 Ca  40  Ca.pbe-nsp-van.UPF
 O  16  O.pbe-rrkjus.UPF

ATOMIC_POSITIONS {crystal}
 Mn1 0.00    0.00    0.00
 Mn2 0.50    0.50    0.50
 Ca  0.25    0.25    0.25
 Ca  0.75    0.75    0.75
 O   0.75    0.25    0.75
 O   0.25    0.75    0.25
 O   0.75    0.75    0.25
 O   0.25    0.25    0.75
 O   0.25    0.75    0.75
 O   0.75    0.25    0.25


K_POINTS {automatic}
6 6 6 0 0 0
EOF
./dopw <scf.in >G-scf-10-insulator.out

#for phonon-antimag
cat>ph.in <<EOF
for phonon
&INPUTPH
amass(1)    = 55
amass(2)    = 55
amass(3)    = 40
amass(4)    = 16
outdir      = "./"
prefix      = 'G-10-scf' !must be the same with scf
ldisp       = .FALSE.
niter_ph    = 100
tr2_ph      = 1.0e-12
alpha_mix(1)= 0.7
nmix_ph     = 4
iverbosity  = 1
fildyn      = 'matdyn'
!epsil       = .TRUE. !Do not set epsil to .true. for metallic system or
q/=0
!trans       = .TRUE. !if trans .and. epsil effective charges are calculated
zue         = .TRUE. !the same with epsil and trans = TRUE
/
0.0 0.0 0.0
EOF

./doph <ph.in >G-ph-10.out
done
#*********************************************

-- 
____________________________________
Hui Wang
School of physics, Fudan University, Shanghai, China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.democritos.it/pipermail/pw_forum/attachments/20110119/0bff5c39/attachment-0001.htm 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: output.rar
Type: application/rar
Size: 28147 bytes
Desc: not available
Url : http://www.democritos.it/pipermail/pw_forum/attachments/20110119/0bff5c39/attachment-0001.rar 


More information about the Pw_forum mailing list