Dear Users,<br><br>Could anyone please let me know how to write a new constraint on the centre of mass of a molecule.<br>I
want to allow an approaching molecule onto a substrate to move along XY
plane while its Z-component of centre of mass is arrested.<br>
One can arrest the Z-component of an atom by setting &#39; if_pos(nat)(3)=0 &#39; for the atom with index number &#39;nat&#39;. <br>Similarly, can we arrest the Z-component of the centre of mass of a molecule ?<br><br>Thanks and Regards,<br>


<br>Sivaramakrishna J<br>Member of Technical Staff<br>Computational Research Labs, Pune.<br>India