[Pw_forum] All Lowdin charges are zero
Semen Yesylevskyy
yesint at mail.ru
Fri Feb 24 16:30:47 CET 2017
Dear All,
I'm new to QE so please forgive if I'm missing something very obvious. I'm trying to study the charge transfer between organic molecules and MoS2 monolayer. I optimized the system, did an scf and then tried to get Lowdin charges with projwcf. However, all charges are almost zero which is not what I can expect:
Lowdin Charges:
Atom # 1: total charge = 0.0000, s = 0.0000,
Atom # 1: total charge = 0.0000, p = 0.0000, pz= 0.0000, px= 0.0000, py= 0.0000,
Atom # 2: total charge = 0.0001, s = 0.0000,
Atom # 2: total charge = 0.0001, p = 0.0001, pz= 0.0000, px= 0.0000, py= 0.0000,
Atom # 3: total charge = 0.0001, s = 0.0000,
Atom # 3: total charge = 0.0001, p = 0.0001, pz= 0.0000, px= 0.0000, py= 0.0000,
Atom # 4: total charge = 0.0000, s = 0.0000,
Atom # 4: total charge = 0.0000, p = 0.0000, pz= 0.0000, px= 0.0000, py= 0.0000,
Atom # 5: total charge = 0.0001, s = 0.0000,
Atom # 5: total charge = 0.0001, p = 0.0001, pz= 0.0000, px= 0.0000, py= 0.0000,
Atom # 6: total charge = 0.0001, s = 0.0000,
....
Could you please give me some advice of what I'm doing wrong?
--------------------------------
The projwcf input:
&PROJWFC
prefix = './3cell-cb_scf' ,
outdir = '.' ,
ngauss = 0.1 ,
DeltaE = 0.01 ,
/
--------------------------------
The scf input:
&CONTROL
calculation = 'scf' ,
restart_mode = 'from_scratch' ,
outdir = '.' ,
wfcdir = './tmp' ,
pseudo_dir = '.' ,
prefix = '3cell-cb_scf' ,
nstep = 1000 ,
tstress = .false. ,
/
&SYSTEM
ibrav = 4 ,
A = 9.495 ,
B = 9.495 ,
C = 100.0 ,
cosAB = 0.5 ,
cosAC = 0 ,
cosBC = 0 ,
nat = 39,
ntyp = 5,
ecutwfc = 35 ,
ecutrho = 140 ,
nosym = .true. ,
tot_charge = 0.000000 ,
/
&ELECTRONS
electron_maxstep = 200 ,
scf_must_converge = .false. ,
diagonalization='cg' ,
/
&IONS
ion_dynamics = 'bfgs' ,
/
ATOMIC_SPECIES
Mo 95.96 Mo_ONCV_PBE-1.0.upf
S 32.06 S_pbe_v1.2.uspp.F.UPF
C 12.011 C_pbe_v1.2.uspp.F.UPF
H 1.008 H.pbe-rrkjus_psl.0.1.UPF
Cl 35.45 cl_pbe_v1.4.uspp.F.UPF
ATOMIC_POSITIONS angstrom
Mo 1.392123361 1.675146760 5.927158394
S 1.409356288 3.492127081 4.353797199
S 1.409009361 3.492914431 7.500696022
Mo 2.974783806 4.416004278 5.927352435
S 2.991861295 6.232688291 4.353114990
S 2.991733615 6.233889787 7.500216780
Mo 4.557309411 7.156655993 5.926482117
S 4.574550104 8.974187415 4.353268832
S 4.573865790 8.974308816 7.499997172
Mo 4.556900653 1.675497042 5.926861427
S 4.573955171 3.492114510 4.352949498
S 4.574377417 3.492858021 7.500504130
Mo 6.139333906 4.416035304 5.926363713
S 6.156932316 6.232998289 4.352535948
S 6.156651090 6.233324282 7.500124138
Mo 7.721896396 7.156797182 5.926687354
S 7.739160210 8.974228888 4.353457824
S 7.739048493 8.974459743 7.500654977
Mo 7.722031044 1.675217921 5.926454615
S 7.739331052 3.492171149 4.352920557
S 7.738377057 3.492639347 7.499947690
Mo 9.304280720 4.416123626 5.927097215
S 9.321652087 6.232837357 4.353059451
S 9.321594176 6.234116854 7.500174405
Mo 10.887047440 7.156845127 5.926347019
S 10.903965611 8.974283787 4.353226240
S 10.903829602 8.974680010 7.500128022
C 5.700416780 4.960942076 11.705590151
C 5.786661978 6.354516349 11.705273773
C 7.046391983 6.956535202 11.705217371
C 8.203237318 6.173116561 11.705190411
C 8.099168585 4.780492709 11.705023203
C 6.846678887 4.164104064 11.704869934
Cl 4.129012811 4.204270654 11.705957792
H 4.877613854 6.955359406 11.704377918
H 7.118320846 8.046628444 11.704165169
H 9.184749419 6.647034640 11.704328004
H 8.998646725 4.161998637 11.704137157
H 6.757143341 3.077849963 11.704290956
K_POINTS automatic
1 1 1 1 1 1
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20170224/c75f69f1/attachment.html>
More information about the users
mailing list