[gmx-users] Implementing Radial Distance Constraints

[David van der Spoel]

Bob Johnson wrote:
> Hello everyone,
> My system consists of a carbon nanotube and various chemicals that adsorb to the
> tube. I would like to use radial distance constraints to enforce adsorbates to
> maintain a fixed distance away from the nanotube. These restraints must obey
> the relation R = sqrt (x^2 + y^2) = user defined constant. With this type of
> restraint the adsorbates are free to move along the z and circumferential
> directions (something that can not be done with normal position or distance
> restraints).
> 
> I would like to add a restraint potential given by V = 1/2*k*(R-R0)^2 = 1/2*k*(
> sqrt(x^2+y^2) - R0)^2. Here R0 is the fixed radial distance. Is this difficult
> to implement in Gromacs? What file in the source code would I edit?
> Thanks,
> Bob Johnson
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
please be more clear.

if we assume your cnt is centered along the z-axis, given roughly by 
x^2+y^2 = d0. Am I correct in assuming that you want all other particles 
to have x^2+y^2 <= d1, where d1 < d0 ?

-- 
David.
________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  	75124 Uppsala, Sweden
phone:	46 18 471 4205		fax: 46 18 511 755
spoel at xray.bmc.uu.se	spoel at gromacs.org   http://folding.bmc.uu.se
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++