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 ' if_pos(nat)(3)=0 ' for the atom with index number 'nat'. <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