<div dir="ltr"><p class="MsoNormal">Hello,</p><p class="MsoNormal">I have encountered a bug (feature?) within CP2K MC (Monte
Carlo) code version 5. The bug showed up when I tried to run GEMC_NVT (two
boxes) in order to get a vapor-liquid coexistence curve.<br></p><p class="MsoNormal">To demonstrate the bug, I have prepared two identical boxes
containing 1 water molecule each. Any MC move performed on water molecules should
50% end up in box 1, and 50% in box2.<br></p><p class="MsoNormal">The bug is related to ‘PMTRAION’ and ‘PMTRAION_MOL’ flags (probability
of attempting a conformational change)<br></p><p class="MsoNormal">MOTION/MC/MOVE_PROBABILITIES/<wbr>MOL_PROBABILITIES/PMTRAION_MOL<br></p><p class="MsoNormal">I observed that when I performed 5000 PMTRAION_MOL moves, about 75% of the
moves ended up in box1, and 25% ended up in box 2. This did not happen with
other moves like PMSWAP (swapping between 2 boxes) or PMTRANS (translational
move) which gave ~2500 moves for each box. I don’t think this is a problem with
my compilation since the precompiled executable also shows the same problem. (see attached)<br></p><p class="MsoNormal">Could anyone clarify if this behavior is intended?<br></p><p class="MsoNormal">Regards,<br></p><p class="MsoNormal">Cho</p></div>