[Pw_forum] ecut dependence of symmetry operations(?)

Madhura Marathe madhura at jncasr.ac.in
Tue Jan 5 05:24:29 CET 2010


 Dear all,

 While performing convergence calculations w.r.t energy cut-off for a bulk
system, I found that the number of symmetry operations detected by a code
is different depending on the value of ecut. This is very surprising
because for the same system, the symmetries should not depend on the ecut
value.

 For Ecut = 20, 25, 40, 50 and 55 Ry, the output says -
"24 Sym.Ops. (with inversion)";
whereas for ecut = 30, 35, 45 and 60 Ry, the output is -
"12 Sym.Ops. (no inversion)"

 I have performed the same calculations using versions 4.0.4 and 4.1 with
the same results.
The input file is as follows -
*****
 &control
    calculation = 'scf'
    verbosity = 'high'
    restart_mode = 'from_scratch',
    prefix = 'RuGGA',
    tstress=.true.
    tprnfor=.true.
    pseudo_dir = '/home/madhura-data/pseudo/',
    outdir = '/home/madhura-data/tmp/'
/
 &system
    ibrav = 4, celldm(1) = 5.1831, celldm(3) = 1.584,
    nat = 2, ntyp = 1,
    ecutwfc = 20, ecutrho = 160,                     # ecutrho = 8*ecutwfc
    occupations='smearing', smearing='mp',degauss=0.05
    nbnd = 20
/
 &electrons
    diagonalization = 'david', mixing_mode = 'plain',
    mixing_beta = 0.7, conv_thr = 1.0d-12
/
ATOMIC_SPECIES
 Ru 101.07 Ru.pbe-n-van.UPF
ATOMIC_POSITIONS {crystal}
 Ru 0.00       0.00       0.00
 Ru 0.66666667 0.33333333 0.50
K_POINTS automatic
 8 8 8 0 0 0
*****

 Can anybody explain what is happening? Or do I need to file a bug report?
Because except for this problem, energies and stress values seem to be
converging.

 Thanks and regards,
 Madhura.

-- 
Madhura Marathe,
PhD student, TSU,
JNCASR, Bangalore.
India.
Phone No: +91-80-22082835



More information about the users mailing list