<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On Jun 22, 2011, at 8:44 AM, marco.stenta wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><br>The
problem is that I have 47883 water molecules: do I have to repeat this
CONSTRAINT section (or the &G3X3 subsection ) 50 thousand time? the
input would be pretty big.<br>Is there a clever way to SHAKE water molecules in my case?<br></blockquote>yes (see below)</div><div><br><blockquote type="cite">you write in your previous answer <br>MOLECULE 1..15<br><br>but
from the manual I understand MOLECULE takes only one integer: what does
your line mean? can I specify for a range instead then for a single
molecule?<br></blockquote><div>that's right.. I'll see how tough it is to convert this keyword to a list..</div><div><br></div><blockquote type="cite">Is it possible to have the code recognizing a MOLNAME also when one uses an amber topology?<br></blockquote>not easily, because AMBER topology has to be heavily digested unfortunately and molnames have to be changed.. Clearly one can always hack this procedure, but it is not my intention do so at the moment..</div><div><br><blockquote type="cite">I
am forced to use amber topology as I did not find the ambertocharmm
converter (leap2fist ?) anymore and I have to shake water, of course.<br></blockquote>I decided to embed the leap2fist in cp2k: less maintenance.. you can still dump the PSF and PDB after digestion and use them to re-run cp2k with the proper styles..</div><div>At this stage you can decide to give the name you like more to the residues in the PSF, and then use just 1 constraint for WATER (assuming all water have the molname WATER).</div><div><br><blockquote type="cite">More in general: what do you mean by "molecule"?<br>each residue of a protein is assigned to a "molecule": what is the meaning of that?<br>TOPOLOGY > <a href="http://cp2k.berlios.de/manual/CP2K_INPUT/FORCE_EVAL/SUBSYS/TOPOLOGY.html#list_PARA_RES" id="desc_PARA_RES">PARA_RES</a> is set to TRUE<br>what if I set it to FALSE?<br></blockquote>PARA_RES distributes the residues instead of molecules.</div><div><br><blockquote type="cite">I
get the term molecule has nothing to do with its chemical meaning, but
my protein part is composed by more than one sub-units: does that
conflicts with <a href="http://cp2k.berlios.de/manual/CP2K_INPUT/FORCE_EVAL/SUBSYS/TOPOLOGY.html#list_PARA_RES" id="desc_PARA_RES">PARA_RES</a> FALSE?<br>Then I get some strange molecules as well like <br> 1. Molecule kind: _QM_MOL0 Number of atoms: 10657<br>this is not the QM part, nor the QM + link and I am lost <br></blockquote><div><br></div>no input.. no idea..</div><div>cheers,</div><div>Teo</div><div><br><blockquote type="cite"><br>thanks for the help<br>m<br><br><br><div><br class="webkit-block-placeholder"></div>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msg/cp2k/-/YYO4fTrQUIwJ">https://groups.google.com/d/msg/cp2k/-/YYO4fTrQUIwJ</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
To unsubscribe from this group, send email to <a href="mailto:cp2k+uns...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
For more options, visit this group at <a href="http://groups.google.com/group/cp2k?hl=en">http://groups.google.com/group/cp2k?hl=en</a>.<br>
</blockquote></div><br></body></html>