[Pw_forum] ecut

Matteo Cococcioni matteo at umn.edu
Mon May 7 16:20:28 CEST 2007


Hi,

one thing missing in your input file is ecutrho (the energy cutoff for 
the charge density) (read INPUT_PW).
If you use ultrasoft pseudopotential you need to specify this parameter 
besides the ecutwfc.
It is needed to reconstruct the "hard" (augmentation) part of the charge 
density that you won't obtain
from the square of (soft) wavefunction. Typically, for ultrasoft 
pseudopotential, this parameter is
8 to 12 times the cut-off for the wavefunctions but it can even be 
bigger. If you don't specify
its value the code probably assumes you are using norm conserving 
pseudos and puts
ecutrho = 4.0 * ecutwfc which is not big enough.
I hope this will help.

Matteo


ali kazempoor wrote:
> hi
> i run scf input file for MnAs that has NiAs -type
> (hexagonal)and i use ultrasoft pseudopotentials for Mn
> and As from espersso  and optimize k point but my run
> don't converge at 50 ,55 ,.. rydberg for ecut.please
> correct my input file.
>  &CONTROL
>                        title = MnAs ,
>                  calculation = 'scf' ,
>                 restart_mode = 'from_scratch' ,
>                       outdir =
> '/scratch/kazempour/tmp/MnAs/' ,
>                   pseudo_dir =
> '/home/akazem/espresso-3.2/pseudo/' ,
>                       prefix = MnAs ,
>  /
>  &SYSTEM
>                        ibrav = 4,
>                    celldm(1) = 7.0373,
>                    celldm(3) = 1.5318,
>                          nat = 4,
>                         ntyp = 2,
>                      ecutwfc = 20,
>                      occupations = 'smearing' ,
>                      degauss = 0.05 ,
>                     smearing = 'methfessel-paxton' ,
>                        nspin = 2 ,
>    starting_magnetization(1) = 1.0,
>    starting_magnetization(2) = 0.1,
>  /
>  &ELECTRONS
>             electron_maxstep = 100,
>                     conv_thr = .00001 ,
>  /
> ATOMIC_SPECIES
>    Mn  54.93   Mn.pbe-sp-van.UPF 
>    As  74.92   As.pbe-n-van.UPF 
> ATOMIC_POSITIONS crystal 
>    Mn      0.000000000    0.000000000    0.000000000  
>  
>    Mn      0.000000000    0.000000000    0.500000000  
>  
>    As      0.333333333    0.666666666    0.250000000  
>  
>    As      0.666666666    0.333333333    0.750000000  
>  
> K_POINTS automatic 
>   11 11 9   1 1 1 
>
>
>
>
>  
> ____________________________________________________________________________________
> No need to miss a message. Get email on-the-go 
> with Yahoo! Mail for Mobile. Get started.
> http://mobile.yahoo.com/mail 
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>   




More information about the Pw_forum mailing list