<div dir="ltr">Dear CP2K developers,<br><br><br>This is crossposted from the PLUMED mailing list since they think the issue may be in the interface rather than PLUMED itself. <br>I'm trying to do a steered MD simulation with PLUMED(2.1.2) + CP2K (latest version, attached are some minimal inputs).<br><br>However the simulation crashes for  virtual atoms (COM) constructed from a group of atoms of a certain size. <br>The number of atoms should not have any significant effect on the resulting CV (at least not initially). <br>It is difficult to identify which combations work and which don't.<br><br>(in a system of 63 atoms in total..)<br><br><span style="font-family:courier new,monospace">surface: COM ATOMS=1                       <wbr>                  [works]<br>surface: COM ATOMS=1,2,3,4,5,6,7,8,9,10,11,<wbr>12,13,14          [crashes]<br>surface: COM ATOMS=1,2,3,4,5,6,7,8,9,10,11,<wbr>12,13,14,15       [works]<br>surface: COM ATOMS=1,2,3,4,5,6,7,8,9,10,11,<wbr>12,13,14,15,16    [crashes]<br>surface: COM ATOMS=1-60                               <wbr>       [crashes]</span><br><br><span style="font-family:courier new,monospace">d_gr-pr: DISTANCE ATOMS=surface,61 COMPONENTS<br>restraint: ... <br>  MOVINGRESTRAINT <br>   ARG=d_gr-pr.z KAPPA0=0.1<br>    AT0=+3  STEP0=0000<br>    AT1=-3  STEP1=2000  <br>...<br></span><br><br>This only happens for parallel runs (pure MPI, pure openmp, or a mix of the two).<br>So I guess it is related to some MPI communication problem.<br><br>Did you ever see something similar or have any idea what could go wrong?<br><br>Thanks for any suggestions.<br><br><br>Best,<br>Jaap</div>