[Pw_forum] (no subject)
degironc
degironc at sissa.it
Wed Feb 15 09:23:03 CET 2006
stress converges slowly with cutoff and if you are not at full convergence
setting the variables qcutz=150., q2sigma=2.0, ecfixed=24.0, modifies
the kinetic energy
of the electrons beyond 24 ry adding a large penalty function and the
"effective" cutoff of your
calculation is something around 24 ry.
This affects the total energy and the stress.
Also having smearing instead of fixed occupation makes a difference but
if your system is gapped it should not be large.
In order to compare output of the two codes you should use the same
setting in both.
stefano
CP>
> &SYSTEM
> ibrav = 14,
> celldm(1) = 24.4946,
> celldm(2) = 1.0,celldm(3) = 1.0,celldm(4) = 0.0,celldm(5) = 0.0,celldm(6) = 0.0,
> nat = 64,
> ntyp = 2,
> nbnd=288,nelec=576,nspin=1,
> ecutwfc = 30.0 ,
> ecutrho = 240 ,
> occupations = 'fixed' ,
> xc_type = 'PBE',
>
> qcutz=150., q2sigma=2.0, ecfixed=24.0,
> /
>
PW>
> &SYSTEM
> ibrav = 1,
> celldm(1) = 24.4946,
> nat = 64,
> ntyp = 2,
> ecutwfc = 30.0 ,
> ecutrho = 240 ,
> occupations = 'smearing' ,
> degauss = 0.02 ,
> smearing = 'fermi-dirac' ,
> /
> &ELECTRONS
>
>
More information about the users
mailing list