[Pw_forum] relax vs multiple SCF calculations inconsistency
sarashs
sarashs at ece.ubc.ca
Tue Apr 18 03:55:21 CEST 2017
Hi,
I've posted this question before but wasn't successful in getting an
answer. I have a molecule ((OH)3ZrOSi(OH)3) and I want to find the
equilibrium angle of ZrOSi. The equilibrium angle that I get with a
"relax" calculation differs a lot from the result I get by increasing
the angle and performing an "SCF" calculation. I tried various things
including increasing the ecutwfc and ecutrho and fixing the number of
bands for both calculations but had no success. Is there any chance
someone can tell me what I'm missing?
Here's a sample code my code:
&CONTROL
calculation = "scf",
outdir = "outputs",
prefix ="ZrSi3",
restart_mode="from_scratch"
pseudo_dir = "pseudos",
/
&SYSTEM
nosym = .TRUE.,
ibrav = 0,
nat = 15,
ntyp = 4,
ecutwfc = 600,
ecutrho = 6000,
vdw_corr = "grimme-d2",
/
&ELECTRONS
conv_thr = 1.D-11,
mixing_beta = 0.5D0,
/
ATOMIC_SPECIES
O 15.9994 O.pz-mt.UPF
H 1.00794 H.pz-vbc.UPF
Si 28.0855 Si.pz-vbc.UPF
Zr 91.224 Zr.pz-hgh.UPF
ATOMIC_POSITIONS { angstrom }
H -0.217652040 2.144165004 -0.967842103
O -0.644926428 1.299751587 -1.173517450
H -0.763682000 -2.102234442 -0.637701677
O -0.631971073 -1.266428593 -1.106959840
H -2.244617658 0.059203828 1.041286276
O -1.287708239 0.028876445 1.183017241
Si -0.341076411 0.043743634 -0.163179170
O 1.183131599 0.097692502 0.401792709
Zr 2.959296806 0.743741575 0.268340451
H 3.161927401 3.483669731 0.077078720
O 2.777805281 2.601338899 0.023924414
H 4.228044154 0.110551898 2.670358608
O 4.003457303 0.315517992 1.757071324
H 3.597999032 -0.225640653 -2.168546302
O 3.766434111 -0.029092443 -1.238766538
K_POINTS { gamma }
CELL_PARAMETERS { angstrom }
13 0 0
0 13 0
0 0 13
Regards,
S. Arash Sheikholeslam
PhD candidate at electrical engineering department, University of
British columbia
More information about the users
mailing list