normal Bond/angle distribution not matching for standard C1 chain

  • yogi@martini
  • yogi@martini's Avatar Topic Author
  • Offline
  • Junior Boarder
2 years 1 month ago #7470 by yogi@martini
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 -

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

2 years 1 month ago #7472 by Pim
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.093 seconds