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