Dear Pwscf users,<br>

<br>

I want to perform a vc-relax calculation for a low-symmetry structure
(P-1), but the calculation is not always convergenced. Maybe this
structure is far away from its equilibrium. From the help files I know
the time step flag "dt"  is a important parameter for variable
cell relaxing calculation. In the manual, it is said the good values of
"dt" is 50/70 or 100/150, but the default value in INPUT_PW file is
20.0D0. How to understand them?  Here, a part of my input file is
attached. The initial atomic position is determined by performing
"relax" calculation.<br>

<br>

&CONTROL<br>

   calculation =   "vc-relax"   ,<br>

   restart_mode = 'from_scratch' ,<br>

   outdir='~/tmp/' ,<br>

   pseudo_dir = '~/espresso32/pseudo' ,<br>

   disk_io = 'default' ,<br>

   verbosity = 'default' ,<br>

   tstress = .true. ,<br>

   tprnfor = .true. ,<br>

   nstep =  55  ,<br>

               etot_conv_thr = 1.0E-5  ,<br>

               forc_conv_thr = 1.0D-4 ,<br>
                         
dt = 50 ,<br>

/<br>

&SYSTEM<br>

  ibrav     = 0,<br>

  celldm(1) = 12.2409,<br>

  nosym     = .true.,<br>

  nat       = 30,<br>

  ntyp      = 3,<br>

  ecutwfc   = 25.D0,<br>

  ecutrho   = 150.D0,<br>

/<br>

&ELECTRONS<br>

  electron_maxstep = 500,<br>

  conv_thr    = 1.D-7,<br>

  mixing_beta = 0.7D0,<br>

  diagonalization   =   'david'   ,<br>

/<br>

&IONS<br>

/<br>

&CELL<br>

  cell_dynamics = 'damp-w',<br>

          press = 0.00,<br>

      wmass = 0.007,<br>

/<br>

...<br>

...<br>

K_POINTS automatic<br>

5 5 5 1 1 1<br><br>
Best regards,<br>
C. H. Hu<br clear="all"><br>-- <br>============================<br>C.H. Hu<br>Postdoctoral fellow<br>Chimie et Physico-Chimie appliquees<br>Institut Francais du Petrole (IFP)<br>Rueil-Malmaison, France<br>E_mail: <a href="mailto:chaohao.mse@gmail.com">
chaohao.mse@gmail.com</a><br>============================