How to obtain factors for the virtual site constructions? #25
-
|
In the parameter file for caffeine, [constraints] [dihedrals] [virtual_sites2] [virtual_sites3] In the [virtual_site2], how do we get the factor 'a'? Is factor 'a' equal to 0.38 (i.e., the distance of bead 3 and bead 5) divided by 0.599 (i.e., the distance of bead 3 and bead 6)? And for the factors 'a' and 'b' in the [virtual_site3], I would very much like to know how to calculate them too. |
Beta Was this translation helpful? Give feedback.
Replies: 2 comments
-
|
Hi @Clare-xiong and thanks for the question. The factors for the virtual sites constructions are obtained in the same way as the bond length and other bonded parameters, i.e., by comparing the mapped distributions coming from all-atom simulations to the CG distributions, as described in tutorial to parametrize new small-molecules (around section 7)). Namely, in deriving the factors for building the virtual sites, you should consider them as actual particles and hence have bond lengths between the particles and the virtual sites and match such bond lengths (even though you won't have such bond lengths in the final model). Actually, a chapter we recently published describes in quite some detail exactly the making of the caffeine model - see in particular Fig 1.2a and Fig1.c and associated discussions. In principles, given that the virtual sites are built as linear combinations of the particle positions, you can probably devise a cleverer way to obtain those factors. However, you'll still need to extract bond lengths between the real particles as described in the tutorial to parametrize new small-molecules. |
Beta Was this translation helpful? Give feedback.
-
|
I know how to get the factors a and b. Thanks! |
Beta Was this translation helpful? Give feedback.
Hi @Clare-xiong and thanks for the question. The factors for the virtual sites constructions are obtained in the same way as the bond length and other bonded parameters, i.e., by comparing the mapped distributions coming from all-atom simulations to the CG distributions, as described in tutorial to parametrize new small-molecules (around section 7)). Namely, in deriving the factors for building the virtual sites, you should consider them as actual particles and hence have bond lengths between the particles and the virtual sites and match such bond lengths (even though you won't have such bond lengths in the final model). Actually, a chapter we recently published describes in quite some de…