[Pw_forum] VC_relax
Ari P Seitsonen
Ari.P.Seitsonen at iki.fi
Thu Mar 31 17:47:51 CEST 2016
Dear Vishal Gupta,
Your cut-off energy looks _very_ low, no matter which elements/pseudo
potentials you do employ. Probably your system becomes metallic upon the
relaxation, so I would include some kind of broadening of the occupation
numbers, in addition possibly using the mixing mode 'local-TF', if you
have vacuum in your system.
Greetings,
apsi
-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-
Ari Paavo Seitsonen / Ari.P.Seitsonen at iki.fi / http://www.iki.fi/~apsi/
Ecole Normale Supérieure (ENS), Département de Chimie, Paris
Mobile (F) : +33 789 37 24 25 (CH) : +41 79 71 90 935
On Thu, 31 Mar 2016, Vishal Gupta wrote:
> Hi all,I've been running a vc_relax calculation on a 40 atom system. The calculation works fine for some
> cycles. The total force had come down to 0.108. After some time. the force started increasing and now the
> iterations don't converge even in 400 steps. The input file is
> /
> &SYSTEM
> ibrav = 0,
> nat = 40,
> ntyp = 1,
> ecutwfc = 12 ,
> /
> &ELECTRONS
> mixing_beta = 0.3 ,
> diagonalization = 'david' ,
> electron_maxstep = 200 ,
> /
>
> [&Ions/]
> /
> &CELL
> cell_dynamics= 'sd' ,
> cell_dofree='xy' ,
> /
>
> CELL_PARAMETERS angstrom
> 22.265358047804796 0.000000000000000 0.000000000000000
> 0.000000000000000 11.131216979724510 0.000000000000000
> 0.000000000000000 0.000000000000000 10.000000000000000
>
> P.S. Ive tried changing cell dynamics to 'damp-pr' , 'sd' and 'damp-w', but the same problem persists.
> Can anybody Please help me identifying the mistake ?
> Thank You.
> Best Regards,
> Vishal Gupta
>
> B.Tech. 3rd year Mechanical
>
> Indian Institute of Technology Ropar
> Rupnagar (140001), Punjab, India.
> Email :- vishal.gupta at iitrpr.ac.in
> RMML, IIT Ropar
>
>
>
> [487a748cedfb942f92401d52262f336d11784524388949.png]
>
>
More information about the users
mailing list