Dear Ruibin<br><br>The QM box in QMMM is not periodic. <br>However, for such a system I would rather suggest a multiple force environment strategy, where all the subsystems can be periodic.<div>I would set-up 3 force environments:</div><div>A full system MM</div><div>B water QM</div><div>C water MM</div><div>The total energy is then E_A +E_B -E_C, and the forces on water are F_A + F_B -F_C</div><div><br></div><div>regards</div><div>marcella</div><div><br>On Friday, November 16, 2012 12:22:40 AM UTC+1, Ruibin wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Dear CP2K group:
<br>     I am new to CP2K and I have a question about QM/MM PBC calculation. I want to simulate water chain inside a carbon nanotube. I  treat the water with QM and the nanotube carbon with MM, and set the length of both QM and MM box to be to be the same as the nanotube length (~30A). However, I noticed that the QM waters, which are restricted by the nanotube to move along the tube axis,  can't pass through the QM box boundary and goes to the neighboring box. It just gives me warning that I need to add wall to the QM box to hold the QM atoms within the boundary. This is not physically meaningful for my system. I need the QM water to freely move along the nanotube axis by crossing the box boundary back and forth. Is there any way to solve this in CP2K? Or should I turn to GROMACS-CPMD or other packages? Any help will be highly appreciated.
<br>
<br>Thank you
<br>Best wishes
<br>Ruibin</blockquote></div>