<div dir="ltr"><font face="georgia, serif">Hi everyone in cp2k community,</font><div><font face="georgia, serif"><br></font></div><div><font face="georgia, serif">I'm using cp2k4.1 to run a QMMM calculation, but the strange thing is that different result emerged under 1 message passing process and </font></div><div><font face="georgia, serif">72 </font><span style="font-family: georgia, serif;">message passing processes situation. I run same example using 1 core and 72 cores </span><font face="georgia, serif">simultaneously. But after 2000 steps, the program </font></div><div><font face="georgia, serif">crashed in </font><span style="font-family: georgia, serif;">72 cores' simulation. And the error report says</span></div><div><span style="font-family: georgia, serif;"><div class="prettyprint" style="background-color: rgb(250, 250, 250); border-color: rgb(187, 187, 187); border-style: solid; border-width: 1px; overflow-wrap: break-word;"><code class="prettyprint"><div class="subprettyprint"><div class="subprettyprint">#2  0x2AD0BBB5A24F</div><div class="subprettyprint">#3  0xAC14E1 in __input_cp2k_restarts_util_MOD_section_velocity_val_set</div><div class="subprettyprint">#0  0x2AF10E6856F7</div><div class="subprettyprint">#2  0x2B0AD401624F</div><div class="subprettyprint">#3  0xAC14E1 in __input_cp2k_restarts_util_MOD_section_velocity_val_set</div><div class="subprettyprint">#3  0xAC14E1 in __input_cp2k_restarts_util_MOD_section_velocity_val_set</div><div class="subprettyprint">#2  0x2AF10F55924F</div><div class="subprettyprint">#4  0x9EE04F in __input_restart_force_eval_MOD_update_subsys</div><div class="subprettyprint">#4  0x9EE04F in __input_restart_force_eval_MOD_update_subsys</div><div class="subprettyprint">#4  0x9EE04F in __input_restart_force_eval_MOD_update_subsys</div><div class="subprettyprint">#0  0x2B06FD3176F7</div><div class="subprettyprint">#2  0x2B06FE1EB24F</div><div class="subprettyprint">#2  0x2AE46DD2C24F</div><div class="subprettyprint">#5  0x9EF0B3 in __input_restart_force_eval_MOD_update_force_eval</div><div class="subprettyprint">#0  0x2AC0F9FB96F7</div><div class="subprettyprint">#1  0x2AC0F9FB9D3E</div><div class="subprettyprint">#5  0x9EF0B3 in __input_restart_force_eval_MOD_update_force_eval</div></div></code></div><br>Setting about restart part reads</span></div><div><span style="font-family: georgia, serif;"><br></span></div><div><div class="prettyprint" style="font-family: georgia, serif; background-color: rgb(250, 250, 250); border-color: rgb(187, 187, 187); border-style: solid; border-width: 1px; overflow-wrap: break-word;"><code class="prettyprint"><div class="subprettyprint"><div class="subprettyprint">&EXT_RESTART</div><div class="subprettyprint">  RESTART_FILE_NAME ../equi_qmmm/UREA-ZW-1_20000.restart</div><div class="subprettyprint">  RESTART_DEFAULT F</div><div class="subprettyprint">  RESTART_POS</div><div class="subprettyprint">  RESTART_VEL</div><div class="subprettyprint">&END</div></div></code></div><br><font face="georgia, serif">In this simulation, I choose GPW as<font size="2"> the </font></font><span style="color: rgb(0, 0, 0);"><font face="georgia, serif" style="" size="2">electronic structure method. Part of relative setting reads</font></span></div><div><span style="color: rgb(0, 0, 0);"><font face="georgia, serif" style="" size="2"><div class="prettyprint" style="background-color: rgb(250, 250, 250); border-color: rgb(187, 187, 187); border-style: solid; border-width: 1px; overflow-wrap: break-word;"><code class="prettyprint"><div class="subprettyprint"><div class="subprettyprint">  &DFT</div><div class="subprettyprint">    CHARGE 0</div><div class="subprettyprint">    BASIS_SET_FILE_NAME BASIS_SET</div><div class="subprettyprint">    POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div class="subprettyprint">    &QS</div><div class="subprettyprint">      METHOD GPW</div><div class="subprettyprint">      EPS_DEFAULT 1.0E-10</div><div class="subprettyprint">    &END QS</div><div class="subprettyprint">    &XC</div><div class="subprettyprint">        &XC_FUNCTIONAL PADE</div><div class="subprettyprint">        &END XC_FUNCTIONAL</div><div class="subprettyprint">    &END XC</div><div class="subprettyprint">    &SCF</div><div class="subprettyprint">      MAX_SCF 30</div><div class="subprettyprint">      EPS_SCF 1.0E-6</div><div class="subprettyprint">      SCF_GUESS ATOMIC</div><div class="subprettyprint">      &OT</div><div class="subprettyprint">        MINIMIZER DIIS</div><div class="subprettyprint">        PRECONDITIONER FULL_SINGLE_INVERSE</div><div class="subprettyprint">      &END</div><div class="subprettyprint">      &OUTER_SCF</div><div class="subprettyprint">        EPS_SCF 1.0E-6</div><div class="subprettyprint">        MAX_SCF 10</div><div class="subprettyprint">      &END</div></div></code></div><br><br></font></span></div><div><font face="georgia, serif">Prior to this simulation, I also run same example using PM6 as QM method. But there is no such problem even under MPI running. </font><br></div><div><font face="georgia, serif"><br></font></div><div><font face="georgia, serif">So where is the problem ? Could someone give any suggestion ?</font></div><div><span style="font-family: georgia, serif;"><br></span></div><div><span style="font-family: georgia, serif;"><br></span></div><div><font face="georgia, serif"><br></font><div><font face="georgia, serif"><br></font></div></div></div>