normal Bond/angle distribution not matching for standard C1 chain

  • yogi@martini
  • yogi@martini's Avatar Topic Author
  • Offline
  • Junior Boarder
More
9 months 3 weeks ago #7470 by yogi@martini
yogi@martini created the topic: Bond/angle distribution not matching for standard C1 chain
Hello experts,

I am parameterizing new lipid with hair pin structure, two tails attached to one head.
I have atomistic trajectories for lipid bilayer (simulation using gromacs with gromos FF). Now I am developing martini model. After 4-1 mapping, I kept standard bond length 0.47 and 1280 force.c, and 180 angle with 25 force.c. for tails. And assigned bond/angle for headgroup beads from atomistic simulation. I did self assembly simulation. I am getting good agreement between AA and CG for headgroup beads but not for tail atoms which are simple C1-C1-C1... tail. The bond angle distribution which I am getting can be found here - https://imgur.com/a/Yx7GY

Please Log in or Create an account to join the conversation.

More
9 months 2 weeks ago #7472 by Pim
Pim replied the topic: Bond/angle distribution not matching for standard C1 chain
Hey, quick reply:
I'd say that your CG results, although the bond length is a bit lower than expected maybe, are what is to be expected from such a simulation. I'm a bit surprised by the very narrow atomistic distributions to be honest. Atomistic simulations of hexadecane, for example, give average angles of 137 deg. according to the 2004 Martini paper (probably also not very hard to check and I'm sure there's more experts here that can comment). So it's a bit stereotypical to say, but are you sure your gromos simulations are accurate? If you want to match those you'll have to deviate from the 0.47/1250 and 180/25 distributions for this lipid.

Please Log in or Create an account to join the conversation.

Time to create page: 0.096 seconds