<div dir="ltr">Dear All:<div><br></div><div>I am running an AIMD simulation with CP2K for a system of water molecules inside carbon nanotube (CNT). All the CNT atoms are fixed with CONSTRAINT/FIXED_ATOMS module. Then I added GLE Thermostat to the movable atoms with THERMOSTAT/DEFINE_REGION module. The corresponding parts of the input file are as the following.</div><div><br></div><div><div> &CONSTRAINT</div><div>  &FIXED_ATOMS</div><div>  Â COMPONENTS_TO_FIX XYZ</div><div>  Â LIST 841..1008</div><div>  &END FIXED_ATOMS</div><div> &END CONSTRAINT</div></div><div><br></div><div><div>  Â &THERMOSTAT</div><div>  Â  Â REGION DEFINED</div><div>  Â  Â &DEFINE_REGION</div><div>  Â  Â  Â LIST 1..840</div><div>  Â  Â &END DEFINE_REGION</div><div>  Â  Â TYPE GLE</div></div><div>  Â  Â ...</div><div>  Â  Â ...</div><div>  Â Â &END THERMOSTAT</div><div><br></div><div>The whole set of input files are attached.</div><div><br></div><div>I got the following error: </div><div><br></div><div><div>*******************************************************************************</div><div> * Â  ___ Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  *</div><div> * Â / Â  \ Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â *</div><div> * [ABORT] Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  *</div><div> * Â \___/ Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  CPASSERT failed Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â *</div><div> * Â  Â | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â *</div><div> * Â O/| Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â *</div><div> * /| | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â *</div><div> * / \ Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â motion/thermostat/thermostat_types.F:282 *</div><div> *******************************************************************************</div><div><br></div><div><br></div><div> ===== Routine Calling Stack =====</div><div><br></div><div>  Â  Â  Â  Â  Â  2 qs_mol_dyn_low</div><div>  Â  Â  Â  Â  Â  1 CP2K</div></div><div><br></div><div><br></div><div>When I remove the FIXED_ATOMS constraint, seems the simulation could run. Could you please help me with the setup of the Thermostat part with fixed atom system? Thanks in advance ~</div><div><br></div><div>Tao</div></div>