[gmx-users] problem with thermodynamic integration simulation

[Parag Mukhopadhyay]

I am trying to run thermodynamic integration
simulation to calculate the relative free energy of
slovation between 2-methyl-3-phenylpropylamine( state
A) and 3-phenylpropylamine(state B) by changing carbon
of CH3 (in A) to hydrogen and methyl hydrogen atoms to
dummy atoms (in B)

 For the simulation in water at lamda=1.00 the
molecule if breaking apart. The log file shows
relative constarint deviation after LINCS. However,
the simulation in vacuo (at lamda= 1.00) works fine
and I don't have any problem with any of the
intermediate lamda (and lamda = 0.00) values in either
water or vacuo simulations.

I have tried energy minimizing and then used the
output structure for the md run but it does not help.
I am using OPLS parameters for the molecule (.top file
attached). I have also attached the mdp options I am
using. I will really appreciate any suggestion on how
to fix the problem.


		
__________________________________________________________
How much free photo storage do you get? Store your friends 'n family snaps for FREE with Yahoo! Photos http://in.photos.yahoo.com
-------------- next part --------------
A non-text attachment was scrubbed...
Name: psa-ppa-new.top
Type: application/octet-stream
Size: 9155 bytes
Desc: 1210261089-psa-ppa-new.top
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20050721/21ac0fe1/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: equ-l0.00.mdp
Type: application/octet-stream
Size: 1183 bytes
Desc: 3643838528-equ-l0.00.mdp
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20050721/21ac0fe1/attachment-0001.obj>