<html dir="ltr"> <head> <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> <style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style> </head> <body ocsi="0" fpstyle="1"> <div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear Jacopo,<br> <br> the problem comes from the fact that you are trying to read molecule names from AMBER topology files. This kind of files only store information about residue names, and therefore CP2K assigns new (different) MOLNAMEs. <br> <br> Have a look at the previous discussion about this issue in the mailing list (Title:<span id="thread_subject_site"> "Constraint warning with AMBER PDB file</span>") where you can find a way to circumvent the problem based on the dumping of new PSF and PDB files which you can modify to give the MOLNAME information expected by cp2k.<br> <br> Cheers,<br> Pablo.<br> <br> <br> <div style="font-family: Times New Roman; color: rgb(0, 0, 0); font-size: 16px;"> <hr tabindex="-1"> <div style="direction: ltr;" id="divRpF649622"><font color="#000000" face="Tahoma" size="2"><b>From:</b> cp...@googlegroups.com [cp...@googlegroups.com] on behalf of Jacopo Sgrignani [sgr...@gmail.com]<br> <b>Sent:</b> Monday, October 17, 2011 3:46 PM<br> <b>To:</b> cp...@googlegroups.com<br> <b>Subject:</b> Re: [CP2K:3533] restraint dihedral<br> </font><br> </div> <div></div> <div>Dear Pablo<br> The input files are in the tar archive. <br> Actually, I found a post in this list about a similar problem with AMBER prmtop files.<br> Then I tried to put the number (456) of the molecule analysing the output of the MOLECULE print_key (it is the .mol file in the tar archive).<br> However the message is that the constrainted atoms are not in the 456 molecule.<br> Maybe I'm reading the molecule number in a wrong place.<br> <br> <br> Thanks again<br> <br> Jacopo<br> <div class="gmail_quote">2011/10/17 Campomanes Ramos Pablo <span dir="ltr"><<a href="mailto:pablo.ca...@epfl.ch" target="_blank">pablo.ca...@epfl.ch</a>></span><br> <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,<br> <br> I think that the message is clear. The code cannot find the MOLNAME defined for constraints, and therefore I assume that constraints are not applied. However, it is difficult to guess where your problem is (pdb, topology file...) without having a look at your input(s) file(s)... Could you please post them here?<br> <br> Regards,<br> Pablo.<br> <br> <div style="font-family: Times New Roman; color: rgb(0, 0, 0); font-size: 16px;"> <hr> <div style="direction: ltr;"><font color="#000000" face="Tahoma" size="2"> <div><b>From:</b> <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>]<br> </div> <b>Sent:</b> Monday, October 17, 2011 1:44 PM<br> <b>To:</b> <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a><br> <b>Subject:</b> Re: [CP2K:3530] restraint dihedral<br> </font><br> </div> <div> <div></div> <div> <div></div> <div>Thanks Pablo for your suggestion, I'm trying this solution.<br> 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" target="_blank"> MOLECULE</a>, actually<br> I tried both the keyword MOLECULE with the resid number from the PDB input file, and<br> the keyword MOLNAME with the resname fromthe pdb, but they don't seem to be the right solutions.<br> Actually using MOLECULE with the resname the simulation start but in the output I find this warning:<br> <br> <br> *** 13:22:30 WARNING in topology_constraint_util:print_warning_molname ***<br> *** MOLNAME (ASD) was defined for constraints, but this molecule ***<br> *** name is not defined. Please check carefully your PDB, PSF (has ***<br> *** priority over PDB) or input driven CP2K coordinates. In case you may ***<br> *** not find the reason for this warning it may be a good idea to print ***<br> *** all molecule information (including kind name) activating the ***<br> *** print_key MOLECULES specific of the SUBSYS%PRINT section.<br> Then I used also the MOLECULES print_key but the number of the molecule where the constrainted atoms are, it is <br> the same of the resid number in the PDB.<br> <br> Can somebody help me? Is it this waring important or the costraint is applied?<br> <br> <div class="gmail_quote">Thanks again<br> <br> Jacopo<br> <br> 2011/10/16 Campomanes Ramos Pablo <span dir="ltr"><<a href="mailto:pablo.ca...@epfl.ch" target="_blank">pablo.ca...@epfl.ch</a>></span><br> <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,<br> <br> I think you can easily solve your problem if you have a look at:<br> <br> <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><br> <br> Best regards,<br> Pablo.<br> <br> -----------------------------------------------<br> <b style="font-weight: normal;">Pablo Campomanes </b> <p><br> Tel:<a href="tel:%2B%2B41-%280%2921-693-0326" value="+41216930326" target="_blank">++41-(0)21-693-0326</a> <br> Fax: <a href="tel:%2B%2B41-%280%2921-693-0320" value="+41216930320" target="_blank"> ++41-(0)21-693-0320</a> <br> Email: <a href="mailto:pablo.ca...@epfl.ch" target="_blank">pablo.ca...@epfl.ch</a></p> <p><br> </p> <p>Laboratory of Computational<br> Chemistry and Biochemistry <br> BCH 4118 EPF Lausanne<br> CH-1015 Lausanne </p> ------------------------------------------------<br> <br> <div style="font-family: Times New Roman; color: rgb(0, 0, 0); font-size: 16px;"> <hr> <div style="direction: ltr;"><font color="#000000" face="Tahoma" size="2"><b>From:</b> <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>]<br> <b>Sent:</b> Sunday, October 16, 2011 4:09 PM<br> <b>To:</b> <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a><br> <b>Subject:</b> [CP2K:3528] restraint dihedral<br> </font><br> </div> <div></div> <div> <div> <div></div> <div> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB">Dear All</span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB">I would like to know the best way to restrict a dihedral angle during a simulation.</span><span style="font-size: 10pt; font-family: Times;" lang="EN-GB"></span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB">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.</span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB"> </span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB">However, I'm not sure that this is the correct solution to my problem.</span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB"> </span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB">Can somebody suggest me other solutions?</span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB"> </span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB">Thanks</span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB"> </span></p> <p class="MsoNormal"><span style="font-family: Arial; color: black;" lang="EN-GB">Jacopo</span></p> <div><br> </div> <div><br> </div> <div><br> </div> <p></p> </div> </div> -- <br> You received this message because you are subscribed to the Google Groups "cp2k" group.<br> To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank"> cp...@googlegroups.com</a>.<br> To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank"> cp2k+unsu...@googlegroups.com</a>.<br> 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>.<br> </div> </div> </div> </div> <font color="#888888"> <p></p> -- <br> You received this message because you are subscribed to the Google Groups "cp2k" group.<br> To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank"> cp...@googlegroups.com</a>.<br> To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank"> cp2k+unsu...@googlegroups.com</a>.<br> 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>.<br> </font></blockquote> </div> <br> <p></p> -- <br> You received this message because you are subscribed to the Google Groups "cp2k" group.<br> To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank"> cp...@googlegroups.com</a>.<br> To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank"> cp2k+unsu...@googlegroups.com</a>.<br> 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>.<br> </div> </div> </div> </div> </div> </div> <div> <div></div> <div> <p></p> -- <br> You received this message because you are subscribed to the Google Groups "cp2k" group.<br> To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank"> cp...@googlegroups.com</a>.<br> To unsubscribe from this group, send email to <a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank"> cp2k+unsu...@googlegroups.com</a>.<br> 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>.<br> </div> </div> </blockquote> </div> <br> <p></p> -- <br> You received this message because you are subscribed to the Google Groups "cp2k" group.<br> To post to this group, send email to cp...@googlegroups.com.<br> To unsubscribe from this group, send email to cp2k+unsu...@googlegroups.com.<br> For more options, visit this group at http://groups.google.com/group/cp2k?hl=en.<br> </div> </div> </div> </body> </html>