<div dir="auto">Dear all,<div dir="auto"><br></div><div dir="auto">I am carrying out high pressure NPT MD simulations using cp.x for a system I am working on. As these simulations are being done for long times they require multiple restarts. While restarting a simulation, I noticed that there is a discontinuity in the fictitious K.E in the new simulation. I started looking into this when, in a particular simulation, there was a very large discontinuity of fictitious K.E (more than 10 meV/atom) on restart and this caused large fluctuation in other parameters such as the enthalpy. The following is a plot of the same (The restart step is marked with a vertical line).</div><div dir="auto"><br></div><div dir="auto"><img src="cid:ii_195921a95aff79b6ec62" style="max-width:100%;height:auto"><br></div><div dir="auto"><br></div><div dir="auto">When I further investigated several other simulations, I observed that a discontinuity in fictitious K.E occurs whenever a simulation is restarted, although in most cases the magnitude of such discontinuity is negligible. Therefore, it went unnoticed previously. An example is shown below:</div><div dir="auto"><br><img src="cid:ii_195921d14acf78d57457" style="max-width:100%;height:auto"><br><br></div><div dir="auto">However, in some cases, the magnitude of this discontinuity is higher, which causes sudden fluctuations in the same and dies down very fast as below:</div><div dir="auto"><br><img src="cid:ii_195921cabaff7a986475" style="max-width:100%;height:auto"><br><br></div><div dir="auto">It can be seen that the magnitude of the discontinuity here is very low compared to the first case. I haven't yet identified the cause of the same and am unsure whether such a discontinuity is to be expected whenever there is a restart. It would be really helpful if someone has an insight on this issue. Also, I would appreciate some help to get an idea on how I can proceed to investigate whether the above behaviour would have affected the dynamics and if so by how much.</div><div dir="auto"><br></div><div dir="auto">The relevant part of the cp.x restart input file is shown below for reference:</div><div dir="auto">&control</div><div dir="auto">    title = 'mol',</div><div dir="auto">    calculation = 'vc-cp',</div><div dir="auto">    restart_mode = 'restart',</div><div dir="auto">    ndr = 77,</div><div dir="auto">    ndw = 78,</div><div dir="auto">    dt    = 5.0d0,</div><div dir="auto">    nstep  = 100000,</div><div dir="auto">    iprint = 5,</div><div dir="auto">    isave  = 100,</div><div dir="auto">    etot_conv_thr = 1.0D-8,</div><div dir="auto">    forc_conv_thr = 1.0D-5,</div><div dir="auto">    ekin_conv_thr = 1.0D-8,</div><div dir="auto">    tstress = .true.,</div><div dir="auto">    tprnfor = .true.,</div><div dir="auto">    prefix = 'mol',</div><div dir="auto">    outdir = './',</div><div dir="auto">    pseudo_dir = '/home/rashidr/pseud/',</div><div dir="auto">    verbosity = 'medium',</div><div dir="auto">    max_seconds = 258700,</div><div dir="auto">/</div><div dir="auto">&system</div><div dir="auto">    ibrav = 0,</div><div dir="auto">    celldm(1) = 15.54772162,</div><div dir="auto">    ntyp = 4,</div><div dir="auto">    nat  = 80,</div><div dir="auto">    ecutwfc = 62,</div><div dir="auto">    ecutrho = 480,</div><div dir="auto">    occupations="fixed",</div><div dir="auto">    vdw_corr= 'grimme-d2'</div><div dir="auto">    nr1b = 30,</div><div dir="auto">    nr2b = 30,</div><div dir="auto">    nr3b = 30,</div><div dir="auto">    ecfixed = 80,</div><div dir="auto">    qcutz = 150,</div><div dir="auto">    q2sigma = 2.0,</div><div dir="auto">/</div><div dir="auto">&electrons</div><div dir="auto">    electron_dynamics = 'verlet',</div><div dir="auto">    emass = 200.d0,</div><div dir="auto">    emass_cutoff = 1.d0,</div><div dir="auto">    electron_velocities = 'default',</div><div dir="auto">    electron_maxstep = 2000,</div><div dir="auto">    conv_thr = 1.0D-10,</div><div dir="auto">    electron_temperature = 'not_controlled',</div><div dir="auto">/</div><div dir="auto">&ions</div><div dir="auto">    ion_dynamics = 'verlet',</div><div dir="auto">    ion_temperature = 'nose',</div><div dir="auto">    tempw= 300.0,</div><div dir="auto">    fnosep = 55.0,</div><div dir="auto">    ion_velocities = 'default',</div><div dir="auto">    ion_positions = 'default',</div><div dir="auto">/</div><div dir="auto">&cell</div><div dir="auto">    cell_dynamics = 'pr',</div><div dir="auto">    cell_parameters = 'default'</div><div dir="auto">    cell_velocities = 'default'</div><div dir="auto">    press = 400.0,</div><div dir="auto">/</div><div dir="auto"><br></div><div dir="auto">Any help regarding the above would be appreciated.</div><div dir="auto"><br></div><div dir="auto">Thanks and Regards,</div><div dir="auto">Rashid Rafeek</div><div dir="auto">PhD student</div><div dir="auto">Department of Chemistry</div><div dir="auto">IISER Bhopal</div></div>