<br><div>Dear CP2K experts,</div><div><br></div><div>I am currently trying to reproduce <a href="https://pubs.acs.org/doi/abs/10.1021/acs.jpcc.9b06348">this article</a> [1], where the authors use Monte Carlo simulations, and AVBMC steps in particular, to build a water cluster from an (ideal) water reservoir. <br></div><div><br></div><div>Since CP2K is the only freely available package that I could find with AVBMC implemented, I decided to give it a try. In other words, I am new to the package. <br></div><div><br></div><div>In the gas reservoir there is no force between the molecules. In the box where the water cluster is being built, the molecules are interacting with a force-field (e.g. TIP3P). Only particles that belong to the cluster should be accepted. We define a maximum allowed distance between the particles in a cluster: If the distances between a particle A and all particles in a cluster is larger than this maximum distance, the particle A does not belong to the cluster (Stillinger cluster criterion). When moving particles from the reservoir to the cluster, we want to reject all moves where a particle is not being attached to the cluster.</div><div><br></div><div>So here comes my main question: Is it possible to define such an acceptance criterion in CP2K?</div><div><br></div><div>I started from <a href="https://www.cp2k.org/howto:gemc">the CP2K GEMC tutorial </a>[2], but cannot find a way to continue. I also have a question about the output files. The trajectory files for the two boxes are both written to the file mc_coordinates. When visualizing it in Ovito, it shows the two boxes every other frame. I want a trajectory file for each box. I tried this:</div><div><br></div><div>&MOTION</div><div>  &PRINT<br>    &TRAJECTORY ON<br>      FILENAME box1.trj<br>      FORMAT XYZ<br>    &END TRAJECTORY<br>  &END PRINT<br>&END MOTION</div><div><br></div><div>but no box1.trj file was written.</div><div><br></div><div>I appreciate any help I can get, I feel quite lost at the moment.<br></div><div><br></div><div>Links:<br></div><div><div>[1] https://pubs.acs.org/doi/abs/10.1021/acs.jpcc.9b06348</div><div>[2] https://www.cp2k.org/howto:gemc</div></div><div><br></div><div>Best regards,</div><div>Even M. Nordhagen</div><div><div>PhD candidate in Computational Physics,</div> <div>Department of Physics, University of Oslo</div></div>