[Pw_forum] scf convergence calculation

Francesca Costanzo costanzo at ms.fci.unibo.it
Wed Jan 23 12:00:25 CET 2008


Dear PWSCF users,

I am doing a scf calculation on si(100) slab surface.
I have tested two pseudo potentials fo silicon to see
which performs better.

I am wondering why  using
the   Si.pbe-n-van.UPF and Si.pw91-n-van.UPF pseudoes
in the first case I can reach quite easily ( after 400 
iterations) the convergence while in the latter case after 1200 iterations 
I still don't reach the convergence.
the mixing beta is 0.5.

Is there something wrong?

Anyway the input file is below:

   # self-consistent calculation
     cat > si-ex01.in << EOF
&control
     calculation='scf'
     prefix='si',
     tstress = .true.
     tprnfor = .true.
     restart_mode='restart',
     pseudo_dir ='/sfs/sanfs/home/userincm/incbo407/PWscf/pseudo/',
     outdir='/scratch/incbo407/au-PWscf/tmp-1/',
     etot_conv_thr = 1.d-4
   /
  &system
     ibrav = 8, celldm(1) = 29.02243, celldm(2) =1.0, celldm(3)= 0.89771, 
nat=56, ntyp=2,  ecutwfc=18.0,
  /
  &electrons
     mixing_mode = 'plain'
     mixing_beta = 0.5
     conv_thr =  1.0d-8
  /
  &ions
  /
ATOMIC_SPECIES
Si 28.086 Si.pw91-n-van.UPF
H   1     H.pw91-van_ak.UPF
ATOMIC_POSITIONS {angstrom}
  Si -.486986953827 -.486915661371 .021837699340
  Si 3.254734080912 3.254554580350 -.715466406181
  Si .424995921464 5.989287334956 -.032510900988
  Si 1.908574447668 7.596205581111 -.896329755432
  Si 5.990175519344 .425574367605 -.031655667742
  Si 7.597289525085 1.909448358802 -.895741285772
  Si 4.863694686691 4.863621862511 -.928089502092
  Si 9.019734606435 9.020000340300 -.836996926760
  Si 1.332006439043 4.282396930595 -1.483511378893
  Si 4.282762677873 1.331640602610 -1.482720710139
  Si 1.215421638705 9.726643706866 -1.436019240017
  Si 4.069065518611 6.940868463489 -1.558148325068
  Si 6.941730748367 4.070143770542 -1.556740812199
  Si 9.727400143869 1.215885487447 -1.436139812232
  Si 6.926671693616 9.825186228268 -1.480679150046
  Si 9.825040961721 6.927019069098 -1.481122837636
  Si 2.752199364097 -.023268345504 -2.733402882614
  Si -.022185141426 2.752036538539 -2.734578963613
  Si 2.755879272416 5.496899501655 -2.970326566691
  Si -.004779806850 8.317568614429 -2.969974442472
  Si 8.317316174588 -.004017804657 -2.969192076227
  Si 5.498360672151 2.755878723458 -2.968474960855
  Si 8.346177303503 5.558752009571 -2.782863187750
  Si 5.557877655641 8.346312021303 -2.782400281184
  Si 1.389966664098 1.388183764950 -4.139518840077
  Si 1.366337053897 6.940229376004 -4.258916306178
  Si 4.151574127028 4.149476803004 -4.278767316434
  Si 4.138929512730 9.694066210006 -4.160002340210
  Si 6.940071554385 1.366577579401 -4.258586275398
  Si 9.693861222560 4.139131304691 -4.160296524367
  Si 6.925553051516 6.926722417990 -4.154730726999
  Si 9.700291030289 9.701542440060 -4.273212145658
  Si .000000000000 .000000000000 -5.538591342479
  Si 2.769296729594 2.769296729594 -5.538591342479
  Si .000000000000 5.538591342479 -5.538591342479
  Si 2.769296729594 8.307892305490 -5.538591342479
  Si 5.538591342479 .000000000000 -5.538591342479
  Si 8.307892305490 2.769296729594 -5.538591342479
  Si 5.538591342479 5.538591342479 -5.538591342479
  Si 8.307892305490 8.307892305490 -5.538591342479
  H .847478422628 -.846884685755 -6.504440365581
  H -.847316494390 .847694326946 -6.503794769337
  H 3.610082118376 1.928241460414 -6.515886469477
  H 1.917014435899 3.621363118971 -6.495428477030
  H .852174341536 4.685717428620 -6.494724671289
  H -.840461532306 6.379649786722 -6.516156349874
  H 3.616235391428 7.460033933598 -6.504011732009
  H 1.921872283045 9.155369669764 -6.503953522512
  H 6.383803828126 -.844617690420 -6.508975414605
  H 4.685555500382 .852282293694 -6.494454790892
  H 9.156068183733 1.921926259125 -6.503307926268
  H 7.460679529841 3.616289367507 -6.504710245978
  H 6.385587155455 4.690845156163 -6.504170485184
  H 4.690683227924 6.386073998525 -6.503524888940
  H 9.160391561857 7.455445966849 -6.494835798512
  H 7.462679819843 9.153210626588 -6.508329818361
K_POINTS {automatic}
4   4   1  0  0  0
EOF


thank you.
 		Francesca Costanzo
-- 
Francesca Costanzo, Ph.D.
Dipartimento di Chimica Fisica ed Inorganica
Viale Risorgimento 4
40136 Bologna
tel.0039-051-2093710
fax 0039-051-2093690


More information about the Pw_forum mailing list