<div dir="ltr">Hi Alin,<div><br></div><div>Thanks for the example. I checked some of the regtest but not the one you mentioned. This quite helpful. Yet, if you check my input file, this exactly what I am doing, so I am puzzled as for why it seems the constraints (at least g3x3) are not enforced. This leads to distorted water molecules (I checked the restart file) and then to the bad contacts causing cp2k to crash.</div><div><br></div><div>Regards,</div><div>Pierre <br><br>On Wednesday, June 17, 2020 at 1:02:19 PM UTC+1, Alin Marin Elena wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi Pierre,
<br>
<br>there are some example in tests/Fist/regtest-1-3/h2o-32-<wbr>g3x3-mult.inp
<br>
<br>
<br>Regards,
<br>Alin
<br>
<br>Without Questions there are no Answers!
<br>______________________________<wbr>______________________________<wbr>__________
<br>Dr. Alin Marin ELENA
<br><a href="http://alin.elena.space/" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Falin.elena.space%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNH4RbkuB6jbzC5uP2isWSj179KFxg';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Falin.elena.space%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNH4RbkuB6jbzC5uP2isWSj179KFxg';return true;">http://alin.elena.space/</a>
<br>______________________________<wbr>______________________________<wbr>__________
<br>
<br>On Wed, 17 Jun 2020 at 12:54, Pierre-André Cazade
<br><<a href="javascript:" target="_blank" gdf-obfuscated-mailto="f1Y3z0w8BwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">pie...@gmail.com</a><wbr>> wrote:
<br>>
<br>> Hi Azade,
<br>>
<br>> Thanks for your answer. The molname is defined in the PSF file, in the topology section. The molname seems to be identified correctly as it is present in the output file (attached) in the description of the topology.
<br>>
<br>> Regards,
<br>> Pierre
<br>>
<br>> On Wednesday, June 17, 2020 at 12:20:32 PM UTC+1, Azade Yazdan Yar wrote:
<br>>>
<br>>> Hi Pierre,
<br>>> I know some of the answers to your first question.
<br>>> What you use in the ATOMS should be the index of atoms within that molecule. Constraints will be applied to all molecules with the same MOLNAME or MOLECULE number.
<br>>> From your input file, I cannot see how the SWM4 molecules are defined in the coordinates section. You have to define the MOLNAME in that part as below (what I have as H2O, or MOLECULE 1). Note that if you use MOLECULE, it is the 'molecule kind number' and not molecule number.
<br>>>
<br>>> Hope this helps,
<br>>> Azade
<br>>>
<br>>> &G3X3
<br>>>        MOLECULE  1
<br>>>        ATOMS  1..3
<br>>>        DISTANCES     a b c
<br>>>  &END G3X3
<br>>>
<br>>> O    5.1692824047308950E+00    5.3185376653354242E+00    2.1483827179189252E+01 H2O 1
<br>>>
<br>>> On Wed, Jun 17, 2020 at 12:41 PM Pierre-André Cazade <<a>pie...@gmail.com</a>> wrote:
<br>>>>
<br>>>> Dear CP2K users,
<br>>>>
<br>>>> 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?
<br>>>>
<br>>>> What about VIRTUAL_SITES?
<br>>>>
<br>>>> Finally, the "CONTRAINT_INFO ON" section does not print a single information in the output to help me out.
<br>>>>
<br>>>> I attached my input file for clarity.
<br>>>>
<br>>>> Regards,
<br>>>> Pierre
<br>>>>
<br>>>>
<br>>>> --
<br>>>> You received this message because you are subscribed to the Google Groups "cp2k" group.
<br>>>> To unsubscribe from this group and stop receiving emails from it, send an email to <a>c...@googlegroups.com</a>.
<br>>>> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0a1b04bc-ad4c-4cd0-babc-7f694e836288o%40googlegroups.com" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/0a1b04bc-ad4c-4cd0-babc-7f694e836288o%40googlegroups.com';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/0a1b04bc-ad4c-4cd0-babc-7f694e836288o%40googlegroups.com';return true;">https://groups.google.com/d/<wbr>msgid/cp2k/0a1b04bc-ad4c-4cd0-<wbr>babc-7f694e836288o%<wbr>40googlegroups.com</a>.
<br>>
<br>> --
<br>> You received this message because you are subscribed to the Google Groups "cp2k" group.
<br>> To unsubscribe from this group and stop receiving emails from it, send an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="f1Y3z0w8BwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.<wbr>com</a>.
<br>> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/eb929574-947c-45dd-8f2c-dd767482b7f1o%40googlegroups.com" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/eb929574-947c-45dd-8f2c-dd767482b7f1o%40googlegroups.com';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/eb929574-947c-45dd-8f2c-dd767482b7f1o%40googlegroups.com';return true;">https://groups.google.com/d/<wbr>msgid/cp2k/eb929574-947c-45dd-<wbr>8f2c-dd767482b7f1o%<wbr>40googlegroups.com</a>.
<br></blockquote></div></div>