<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Dear QE users,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">        </div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">            I want to run cp.x for long time (nearly 20-30ps) for a heavy system like LiFePO4. Starting form the minimized structure (after electron minimization & randomization) I checked the dynamics run in microcanonical ensemble which runs very fine (total constant of motion nearly constant asthe criteria suggests). But when I run the simulation using Nose-hoover thermostat for 10-20 ps the electron fictitious kinetic energy increases rapidly which suggests that the system is not proceeding in the right direction. As suggested I minimized the electron ke  again and restated the simulation using NH thermostat. But results after re minimization is not proceeding in the correct way. So, help me how to restart the dynamics for long time step. </div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><div class="gmail_default"><br></div><div class="gmail_default">&control</div><div class="gmail_default">    calculation = 'cp'</div><div class="gmail_default">    nstep  = 50000,</div><div class="gmail_default">    iprint = 10,   isave  = 100,</div><div class="gmail_default">    dt    = 4.0d0,</div><div class="gmail_default">    etot_conv_thr = 1.d-9,   ekin_conv_thr = 1.d-5,</div><div class="gmail_default"> /</div><div class="gmail_default"> &system</div><div class="gmail_default">    nat=  55,  ntyp= 8,</div><div class="gmail_default">    ecutwfc =30.0, ecutrho = 250.0,</div><div class="gmail_default">    nr1b=12,nr2b=12,nr3b=12,</div><div class="gmail_default"> /</div><div class="gmail_default"> &electrons</div><div class="gmail_default">    emass = 400.d0, emass_cutoff = 2.5d0,</div><div class="gmail_default">    orthogonalization = 'ortho',</div><div class="gmail_default">    ortho_eps = 5.d-8, ortho_max = 20,</div><div class="gmail_default">    electron_dynamics = 'verlet', </div><div class="gmail_default">    !electron_damping = 0.1,</div><div class="gmail_default">    electron_velocities ='zero'</div><div class="gmail_default"> /</div><div class="gmail_default"> &ions</div><div class="gmail_default">    ion_dynamics ='verlet',</div><div class="gmail_default">    ion_velocities='zero', ion_temperature = 'nose', fnosep = 24., tempw= 2000.</div><div class="gmail_default"> /</div><div class="gmail_default"> &cell</div><div class="gmail_default">    cell_dynamics = 'none', press = 0.0d0,</div></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Any help would be appreciated.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div><div class="gmail_signature"><div dir="ltr"><div><div><div><div>Shashi Bhusan Mishra<br></div>Research Scholar,<br></div>Dept. Of Physics,<br></div>Room No- HSB209A<br></div>IIT Madras-36<br>9087218085<br></div></div></div>
</div>