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 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 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>============================