<div dir="ltr">Hi guys,<div><br></div><div>As Amandine Chouanard mentioned in a <a href="https://groups.google.com/forum/#!topic/cp2k/mwieBkpEjns">recent post </a>the conserved quantity shows discontinuities when a flexible cell MD simulation is restarted. I used her input file and found that restart works fine with NVT and NPT_I simulations, but NOT with NPT_F (see attached figure). The discontinuities in NVT_F the conserved quantity are larger if restarts are less frequent (i.e. trajectory pieces are longer).</div><div><br></div><div>It seems that the conserved quantities in NPT_I and NPT_F simulations are computed by the same subroutine in motion/md_conserved_quantities.F:390. So it seems that the restart file contains ALL data on the current state of the barostat and thermostat to ensure a smooth restart.</div><div><br></div><div>Is it possible that this has something to do with the CELL_REF setting? Something along the lines that the grids before and after the restart are different? I tried printing the GRID info but it seems to keep the number of points constant...</div><div><br></div><div>I will appreciate any suggestions/comments,</div><div>Rustam</div><div><br></div><div><br></div>







</div>