<div><br></div>Dear Lukas, <div><br></div><div>Did you  first  equilibrate the system at room temperature?  </div><div>Why did you choose the Langevin thermostat? How does the cell fluctuate? Which is the energy of the barostat? The stress tensor?</div><div>Is the total energy conserved? </div><div>Without more information it is hard to guess what the problem might be. </div><div>If you have an initially equilibrated sample (stress tensor around zero, conserved total energy, etc. ), try annealing it progressively to higher temperatures.</div><div><br></div><div>The model looks quite small and it will be anyway difficult to get a disordered structure due to the strong spatial correlations. </div><div><br></div><div>Kind regards</div><div>Marcella </div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, November 16, 2020 at 12:03:49 PM UTC+1 Lukas C wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi everyone!<br><br>I am trying to see signs of melting for a given SiO2 bulk. In order to achieve this, I ran MD calculations for a NPT_F (constant pressure and temperature, flexible cell) ensemble which is connected to a thermal bath (adaptive langevin thermostat) and to a barostat which are utilized to keep T and p at constant values for each simulation respectively. The temperature was set to T=600K,900K,1200K,1500K,1800K,2100K and 3000K. <br><br>However, i am experiencing multiple problems.<br>1. The calculations abort after 20 steps (every time) with the error:<br><br> *******************************************************************************<br> *   ___                                                                       *<br> *  /   \                                                                      *<br> * [ABORT]                                                                     *<br> *  \___/                  invalid value for enumeration:   104                *<br> *    |                                                                        *<br> *  O/|                                                                        *<br> * /| |                                                                        *<br> * / \                                     input/input_enumeration_types.F:188 *<br><div> *******************************************************************************</div><div><br></div><div>2. by comparing the output files, i found out that the simulation seems to trigger the same reactions (altough the temperature is very different). cp2k did recognice the temperature (thus, the kinetic energy is changed) but the curves that represent the change of kinetic energy resemble each other perfectly. to my understanding, this means that exactly the same reactions were triggered, which should not be the case for different temperatures. Each simulation was started with a .wfn file as an initial guess, which was taken from an single-point ENERGY calculation that was conducted beforehand. Thus the potential energy is the same for each simulation at t=0. The plot is attached together with the input file. I already tried to use various optimizers, preconditioners and so on. am i doing something seriously wrong? Thank you, your help is appreciated!<br><br>Regards,<br>Lukas<br></div></blockquote></div>