Dear CP2K users,<br><br>I am trying to do a run using a flexible version
 of a 4-site water model.  I am trying to use the &CONSTRAINT module
 to describe a virtual atom (MW) as a function of the O-H distances but I
 do not want the whole water molecule treated as a rigid body. I have 
specified the virtual atom constraint using the &VIRTUAL_SITE 
definition but I have commented other constraints regarding the water 
molecule. However, I am unable to make it work. Could you please have a 
look at it and let me know what I am doing wrong or whether it is 
possible or not?  <br><br>Thanks very much for any help,<br><br>Carlos<br><br> I have attached the input files and below is the error I get<br><br>*********************************************************<br> *** ERROR in get_molecule (MODULE molecule_types_new) ***<br> *********************************************************<br><br> *** The pointer lci is not associated ***<br><br> *** Program stopped at line number 279 of MODULE molecule_types_new ***<br><br> ===== Routine Calling Stack ===== <br><br>            3 velocity_verlet<br>            2 qs_mol_dyn_low<br>            1 CP2K<br>