This issue was moved to a discussion.
You can continue the conversation there. Go to discussion →
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Methyl-benzoate bonding #9
Comments
Hi and thanks for the question. Yep, the structure looks like your top one. I've uploaded some more sketches to clarify further: https://github.com/ricalessandri/Martini3-mappings/blob/main/methyl-benzoate.pdf This is because the -COOMe group is fairly rigid and this construction is very numerically stable (see, e.g., minute 11:16 of this talk). Feel free to follow up, I'd be happy to clarify the choice further. |
Hi, Thank you very much for the quick response. It makes sense now :). Thank you, Vaidyanathan M. S |
I'm a little confused about the three TC5 beads. The carbon atom immediately connected to the -COOMe group (call it "alpha carbon"). Why do the TC5 map 2 non-hydrogen atoms like this :(1C_beta + 1/2C_alpha + 1/2 * C_gamma)? Why is it not like a benzene(TC5 bead map the -CH=CH- group)? When I map the molecule to the CG model, how should I deal with the benzene? Should I make the TC5 bead map like methyl-benzoate, or the benzene? |
The mapping is chosen to preserve the symmetry of the molecule. You can see for yourself by trying to draw mappings such as |
If I understand correctly, you mean that there can be more than one way to map while guaranteeing that TC5 beads will represent two carbon atoms. The mapping method will be chosen based on considerations of symmetry. Am I right in my understanding? |
Yes, exactly. |
This issue was moved to a discussion.
You can continue the conversation there. Go to discussion →
In the parameter file for methyl-benzoate,
; id type resnr residu atom cgnr charge
1 N4a 1 MBZOA O1 1 0
2 TC5 1 MBZOA R2 2 0
3 TC5 1 MBZOA R3 3 0
4 TC5 1 MBZOA R4 4 0
[bonds]
; i j funct length force.c.
1 2 1 0.389 5000 ; cog
1 3 1 0.389 5000 ; cog
Why is there a bond between N4a and two different TC5 atoms?
Wouldn't the structure then look like
Instead of
The text was updated successfully, but these errors were encountered: