Hello,
I am working with a DPPC lipid bilayer system in water, and I would like
to add a shear force to the bilayer. Preferably, this force would be a
function of the z coordinate and be maximal at the bilayer surface and decay
to zero at the box top and bottom. For this, I have been looking a the
cos_acceleration option, but it doesn't quite give me what I want. The
Manual says that the magnitude of the acceleration is calculated by
cos(2*pi*z/boxheight). For the size of the bilayer and box in my system,
things happen to work out nicely if the magnutude of the acceleration were
calculated as cos(2*pi*z/boxheight - pi/2). This would give me an
acceleration of zero at the box top and bottom as well as the center, and
maximal acceleration at roughly the bilayer top and bottom surface. I was
wondering if there was an easy to make this switch in Gromacs?
If there is no easy to do this, I am assuming that I have to modify the
source code. I found in update.c a line that reads:
cosz = cos(fac*x[n][ZZ])
Would this essentially be the part I am looking for, and I would just have to
change it to: cosz = cos(fac*x[n][ZZ] - PI/2)?
Are there other parts of the source code that I would have to change besides
update.c?
Thank you,
Mike Tomasini
_______________________________________________
gmx-users mailing list [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
