<div dir="ltr">Dear CP2K users,<div><br></div><div>I am trying to use G3X3 constraints on my water molecules during a CELL_OPT motion. Needless to say, that it does not work. The "manual" is not very clear about the indexes that need to be provided. Should I use the index of the atoms in the coordinate file, for example 15761 15762 15763, or the indexes within the water molecule itself: 1 2 3. Does such constraint apply to all my water molecules provided they have the same molname "SWM4" or do I need an entry for each water molecule?</div><div><br></div><div>What about VIRTUAL_SITES?</div><div><br></div><div>Finally, the "CONTRAINT_INFO ON" section does not print a single information in the output to help me out.</div><div><br></div><div>I attached my input file for clarity.</div><div><br></div><div>Regards,</div><div>Pierre</div><div><br></div><div><br></div></div>