Thanks Pablo for your suggestion, I'm trying this solution. However I have a problem with the keyword MOLNAME or <a href="http://cp2k.berlios.de/manual/CP2K_INPUT/MOTION/CONSTRAINT/COLLECTIVE.html#list_MOLECULE" id="desc_MOLECULE">MOLECULE</a>, actually I tried both the keyword MOLECULE with the resid number from the PDB input file, and the keyword MOLNAME with the resname fromthe pdb, but they don't seem to be the right solutions. Actually using MOLECULE with the resname the simulation start but in the output I find this warning: *** 13:22:30 WARNING in topology_constraint_util:print_warning_molname *** *** MOLNAME (ASD) was defined for constraints, but this molecule *** *** name is not defined. Please check carefully your PDB, PSF (has *** *** priority over PDB) or input driven CP2K coordinates. In case you may *** *** not find the reason for this warning it may be a good idea to print *** *** all molecule information (including kind name) activating the *** *** print_key MOLECULES specific of the SUBSYS%PRINT section. Then I used also the MOLECULES print_key but the number of the molecule where the constrainted atoms are, it is the same of the resid number in the PDB. Can somebody help me? Is it this waring important or the costraint is applied? <div class="gmail_quote">Thanks again Jacopo 2011/10/16 Campomanes Ramos Pablo <<a href="mailto:pablo.ca...@epfl.ch">pablo.ca...@epfl.ch</a>> <blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"> <div> <div style="direction: ltr; font-family: Tahoma; color: rgb(0, 0, 0); font-size: 10pt;">Dear Jacopo, I think you can easily solve your problem if you have a look at: <a href="http://cp2k.berlios.de/manual/CP2K_INPUT/MOTION/CONSTRAINT/COLLECTIVE/RESTRAINT.html" target="_blank">http://cp2k.berlios.de/manual/CP2K_INPUT/MOTION/CONSTRAINT/COLLECTIVE/RESTRAINT.html</a> Best regards, Pablo. ----------------------------------------------- Pablo Campomanes Tel:<a href="tel:%2B%2B41-%280%2921-693-0326" value="+41216930326" target="_blank">++41-(0)21-693-0326</a> Fax: <a href="tel:%2B%2B41-%280%2921-693-0320" value="+41216930320" target="_blank">++41-(0)21-693-0320</a> Email: <a href="mailto:pablo.ca...@epfl.ch" target="_blank">pablo.ca...@epfl.ch</a> Laboratory of Computational Chemistry and Biochemistry BCH 4118 EPF Lausanne CH-1015 Lausanne ------------------------------------------------ <div style="font-family: Times New Roman; color: rgb(0, 0, 0); font-size: 16px;"> <hr> <div style="direction: ltr;">From: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> [<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>] on behalf of Jacopo Sgrignani [<a href="mailto:sgr...@gmail.com" target="_blank">sgr...@gmail.com</a>] Sent: Sunday, October 16, 2011 4:09 PM To: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> Subject: [CP2K:3528] restraint dihedral </div> <div></div> <div><div><div></div><div class="h5"> Dear All I would like to know the best way to restrict a dihedral angle during a simulation. Without any experience about this, I tried to use a fictitious Colvar putting two walls to define the space that the dihedral can freely explore. Then, the program asks for a value of the SCALE variable during the metadynamics calculation, and I set this to 0. However, I'm not sure that this is the correct solution to my problem. Can somebody suggest me other solutions? Thanks Jacopo <div> </div> <div> </div> <div> </div> </div></div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>. To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>. For more options, visit this group at <a href="http://groups.google.com/group/cp2k?hl=en" target="_blank">http://groups.google.com/group/cp2k?hl=en</a>. </div> </div> </div> </div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>. To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>. For more options, visit this group at <a href="http://groups.google.com/group/cp2k?hl=en" target="_blank">http://groups.google.com/group/cp2k?hl=en</a>. </blockquote></div>