<div>Hello all,</div><div> </div><div>I'm trying to perform vc-relax calculations using damped MD method. QE version 6.7. Here is the part of my input file:</div><div> </div><div><div>&CONTROL</div><div>  calculation = 'vc-relax'</div><div>  etot_conv_thr =   1.0000000000d-06</div><div>  forc_conv_thr =   1.0000000000d-04</div><div>  outdir = './tmp/'</div><div>  tprnfor = .true.</div><div>  tstress = .true.</div><div>  verbosity = 'low'</div><div>  nstep = 300</div><div>  dt = 50.0d0 ! 2fs timestep</div><div>/</div><div> </div><div>&SYSTEM</div><div>  ibrav = 0</div><div>  nat = 84</div><div>  nspin = 2</div><div>  ntyp = 6</div><div>  ecutwfc =     50.00</div><div>  ecutrho =    400.00</div><div>  occupations = 'smearing'</div><div>  smearing = 'gaussian'</div><div>  degauss = 0.015</div><div>  starting_magnetization(1) =   -1.0</div><div>  starting_magnetization(2) =    1.0</div><div>/</div><div> </div><div>&ELECTRONS</div><div>  conv_thr          =   1.0d-10</div><div>  scf_must_converge =  .false.</div><div>  electron_maxstep  =   50</div><div>  mixing_beta       =   0.2</div><div>  startingpot       = 'file'</div><div>/</div><div> </div><div>&ions</div><div> ion_dynamics    = 'damp'</div><div> ion_temperature = 'reduce-T'</div><div> tempw           = 1000.00d0 !initial temperature</div><div> delta_t          = 30.0d0</div><div> nraise          = 10</div><div> pot_extrapolation = 'second_order'</div><div> wfc_extrapolation = 'second_order'</div><div>/</div><div> </div><div>&cell</div><div> cell_dynamics  = 'damp-w'</div><div> press_conv_thr = 0.05</div><div> cell_dofree    = 'all'</div><div> press          = 10.0</div><div>/</div><div> </div><div>K_POINTS automatic</div><div>2  1  3      0 0 0</div><div>....</div><div> </div><div>So, the goal is to start with T=1000 K and go down by 10 K. When I look at my output file, I see that temperature is 0 and does not change at all:</div><div> </div><div><div>  Ekin =     0.00000000 Ry    T =    0.0 K  Etot =    -2876.41170851</div><div>     Ekin =     0.00000244 Ry    T =    0.0 K  Etot =    -2876.41170845</div><div>     Ekin =     0.00000654 Ry    T =    0.0 K  Etot =    -2876.41170862</div><div>     Ekin =     0.00001080 Ry    T =    0.0 K  Etot =    -2876.41170872</div><div>     Ekin =     0.00001342 Ry    T =    0.0 K  Etot =    -2876.41170901</div><div>     Ekin =     0.00000711 Ry    T =    0.0 K  Etot =    -2876.41171702</div><div>     Ekin =     0.00000840 Ry    T =    0.0 K  Etot =    -2876.41171760</div><div>     Ekin =     0.00000870 Ry    T =    0.0 K  Etot =    -2876.41171916</div><div>     Ekin =     0.00001002 Ry    T =    0.0 K  Etot =    -2876.41171991</div><div>     Ekin =     0.00001177 Ry    T =    0.0 K  Etot =    -2876.41172036</div><div>     Ekin =     0.00001377 Ry    T =    0.0 K  Etot =    -2876.41172073</div><div>     Ekin =     0.00001586 Ry    T =    0.0 K  Etot =    -2876.41172118</div><div>     Ekin =     0.00001820 Ry    T =    0.0 K  Etot =    -2876.41172153</div><div>     Ekin =     0.00002099 Ry    T =    0.0 K  Etot =    -2876.41172161</div><div>     Ekin =     0.00002386 Ry    T =    0.0 K  Etot =    -2876.41172176</div><div>     Ekin =     0.00002671 Ry    T =    0.0 K  Etot =    -2876.41172211</div><div>     Ekin =     0.00002980 Ry    T =    0.0 K  Etot =    -2876.41172240</div><div>     Ekin =     0.00003324 Ry    T =    0.0 K  Etot =    -2876.41172252</div><div>     Ekin =     0.00003683 Ry    T =    0.0 K  Etot =    -2876.41172268</div><div>     Ekin =     0.00004054 Ry    T =    0.0 K  Etot =    -2876.41172293</div><div>     Ekin =     0.00004397 Ry    T =    0.0 K  Etot =    -2876.41172366</div><div>     Ekin =     0.00004755 Ry    T =    0.0 K  Etot =    -2876.41172445</div><div>     Ekin =     0.00005169 Ry    T =    0.0 K  Etot =    -2876.41172487</div><div>     Ekin =     0.00005627 Ry    T =    0.0 K  Etot =    -2876.41172507</div><div>     Ekin =     0.00006112 Ry    T =    0.0 K  Etot =    -2876.41172519</div><div>     Ekin =     0.00006620 Ry    T =    0.0 K  Etot =    -2876.41172529</div><div>     Ekin =     0.00007149 Ry    T =    0.0 K  Etot =    -2876.41172538</div><div>     Ekin =     0.00007702 Ry    T =    0.0 K  Etot =    -2876.41172545</div><div>     Ekin =     0.00008276 Ry    T =    0.0 K  Etot =    -2876.41172550</div><div>     Ekin =     0.00008872 Ry    T =    0.0 K  Etot =    -2876.41172555</div><div>     Ekin =     0.00009486 Ry    T =    0.0 K  Etot =    -2876.41172562</div><div> </div><div>I'm confused. It seems that my temperature settings are ignored. I played with "dt" timestep (increased up to 100), and output shows again initial temperature 0 and after 4 MD steps it jumped to "*****" value and calculations have diverged.</div><div> </div><div>What am I doing wrong? Why the initial temperature was not inforced?</div><div> </div><div>Thanks,</div><div> Sergey</div><div>ISF</div></div></div>