[Pw_forum] conversence problem in vc-relax (Atom Mass in

Cesar R.S. da Silva cesards at msi.umn.edu
Mon Sep 10 19:33:29 CEST 2007


Hi Yukihiro,

if you are using

>     etot_conv_thr = 1.0D-5
>     forc_conv_thr = 1.0D-4

you should set a tighter convergence for SCF. My suggestion is

conv_thr = 1.0d-8


Yours,

Cesar.

>
> Message: 2
> Date: Sun, 09 Sep 2007 01:02:49 +0900
> From: Yukihiro Okuno <okunoyukihiro2000 at yahoo.co.jp>
> Subject: [Pw_forum] conversence problem in vc-relax (Atom Mass in
> 	relaxation)
> To: pw_forum at pwscf.org
> Message-ID: <200709081602.AA00341 at Einstein.yahoo.co.jp>
> Content-Type: text/plain; charset=us-ascii
>
> Dear PWScf Users.
>
> I'm Yukihiro Okuno.
>
> I had posted about the Atom Mass in vc-relax, and thank you for the
> responces.
> We set same atom mass 20.0 ( average of atom mass), and do vc-relax
> calculation as
> the following input file,
>
>
> &control
>     calculation  = 'vc-relax',
>     restart_mode = 'from_scratch',
>     pseudo_dir   = '/pseudo/',
>     outdir       = './',
>     prefix='vc-test'
>     nstep = 100
>     etot_conv_thr = 1.0D-5
>     forc_conv_thr = 1.0D-4
>     dt = 80.0
>     tstress = .true.
>     tprnfor = .true.
>   /
>  &system
>     ibrav=6
>     celldm(1)=6.83994507291
>     celldm(3)=1.09801321247
>     nat=5
>     ntyp=3
>     nbnd=26
>     ecutwfc=60.0
>     ecutrho = 600.0
>     occupations = 'fixed'
>     degauss=0.00
>  /
>  &electrons
>     conv_thr = 1e-10,
>     mixing_beta=0.3,
>  /
> &IONS
>  ion_dynamics='damp'
>  bfgs_ndim=3
>  pot_extrapolation = "second_order"
>  wfc_extrapolation = "second_order"
> /
> &CELL
>   cell_dynamics  = 'damp-w',
>   press = 0.0
>   wmass = 0.02
> /
>
>
> But the calculation is not conversed and oscillate. The system has only
> five atoms
> and simple P4mm perovskite structure.
> I had change wmass as default value, but the convergence is not achieved.
> Is it too hard to converse to set  forc_conv_thr = 1.0D-4 in VC-relax
> calculation ?
> ( I want to get severe internal atomic coordinate, and
> forc_conv_thr=1.0D-3 is not
> sufficient.) Simple atomic internal relaxation, I can get converged
> result.
>
>
> The stress and total force of output file is as follow.
>
> stress values ,,
>
>     entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -0.13
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -2.09
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -1.17
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -0.60
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 0.07
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 0.37
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 0.84
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 0.95
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 1.07
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 1.22
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 1.36
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 1.47
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 1.55
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 0.96
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 0.65
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 0.39
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 0.20
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> 0.04
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -0.08
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -0.13
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -0.13
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -0.63
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -0.89
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -1.05
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -1.04
>      entering subroutine stress ...
>           total   stress  (Ry/bohr**3)                   (kbar)     P=
> -1.13
>
>
> and  The total force,
>
>      Total force =     0.000512     Total SCF correction =     0.000057
>      Total force =     0.000499     Total SCF correction =     0.000063
>      Total force =     0.000472     Total SCF correction =     0.000024
>      Total force =     0.000411     Total SCF correction =     0.000022
>      Total force =     0.000385     Total SCF correction =     0.000053
>      Total force =     0.000527     Total SCF correction =     0.000012
>      Total force =     0.000487     Total SCF correction =     0.000188
>      Total force =     0.000572     Total SCF correction =     0.000010
>      Total force =     0.000537     Total SCF correction =     0.000016
>      Total force =     0.000484     Total SCF correction =     0.000016
>      Total force =     0.000440     Total SCF correction =     0.000007
>      Total force =     0.000407     Total SCF correction =     0.000017
>      Total force =     0.000369     Total SCF correction =     0.000016
>      Total force =     0.000391     Total SCF correction =     0.000045
>      Total force =     0.000552     Total SCF correction =     0.000015
>      Total force =     0.000569     Total SCF correction =     0.000130
>      Total force =     0.000683     Total SCF correction =     0.000026
>      Total force =     0.000623     Total SCF correction =     0.000025
>      Total force =     0.000544     Total SCF correction =     0.000020
>      Total force =     0.000404     Total SCF correction =     0.000015
>      Total force =     0.000280     Total SCF correction =     0.000028
>      Total force =     0.000219     Total SCF correction =     0.000011
>      Total force =     0.000348     Total SCF correction =     0.000039
>      Total force =     0.000552     Total SCF correction =     0.000015
>      Total force =     0.000688     Total SCF correction =     0.000019
>      Total force =     0.000741     Total SCF correction =     0.000013
>      Total force =     0.000723     Total SCF correction =     0.000013
>      Total force =     0.000638     Total SCF correction =     0.000018
>      Total force =     0.000504     Total SCF correction =     0.000017
>      Total force =     0.000368     Total SCF correction =     0.000029
>      Total force =     0.000256     Total SCF correction =     0.000020
>      Total force =     0.001173     Total SCF correction =     0.000016
>      Total force =     0.001841     Total SCF correction =     0.000271
>      Total force =     0.001924     Total SCF correction =     0.000220
>      Total force =     0.002328     Total SCF correction =     0.000057
>      Total force =     0.002448     Total SCF correction =     0.000019
>      Total force =     0.002318     Total SCF correction =     0.000016
>      Total force =     0.002004     Total SCF correction =     0.000018
>      Total force =     0.001547     Total SCF correction =     0.000022
>      Total force =     0.001029     Total SCF correction =     0.000026
>      Total force =     0.000669     Total SCF correction =     0.000015
>      Total force =     0.002834     Total SCF correction =     0.000206
>      Total force =     0.003057     Total SCF correction =     0.000170
>      Total force =     0.001864     Total SCF correction =     0.000273
>      Total force =     0.001107     Total SCF correction =     0.000101
>      Total force =     0.000870     Total SCF correction =     0.000021
>      Total force =     0.000355     Total SCF correction =     0.000068
>      Total force =     0.000613     Total SCF correction =     0.000114
>      Total force =     0.000484     Total SCF correction =     0.000038
>      Total force =     0.000338     Total SCF correction =     0.000180
>      Total force =     0.000372     Total SCF correction =     0.000125
>      Total force =     0.000460     Total SCF correction =     0.000028
>      Total force =     0.000482     Total SCF correction =     0.000014
>      Total force =     0.000528     Total SCF correction =     0.000018
>      Total force =     0.000652     Total SCF correction =     0.000008
>      Total force =     0.000729     Total SCF correction =     0.000017
>      Total force =     0.000740     Total SCF correction =     0.000013
>      Total force =     0.000684     Total SCF correction =     0.000005
>      Total force =     0.000573     Total SCF correction =     0.000005
>      Total force =     0.000446     Total SCF correction =     0.000012
>      Total force =     0.000369     Total SCF correction =     0.000010
>      Total force =     0.000408     Total SCF correction =     0.000005
>      Total force =     0.000602     Total SCF correction =     0.000044
>      Total force =     0.000897     Total SCF correction =     0.000015
>      Total force =     0.001145     Total SCF correction =     0.000020
>      Total force =     0.001247     Total SCF correction =     0.000004
>      Total force =     0.001348     Total SCF correction =     0.000050
>
> Sincerely.
>
>
> ----
> Yukihiro Okuno  okunoyukihiro2000 at yahoo.co.jp
>
>




More information about the Pw_forum mailing list