Hello,
I´m using grimaces 2016 to pull a dihedral angle.
With the angle I can force the flipping in and out of a single bulge in a RNA
system.
I tried with different k.
I do a simulation of 180 ns with a rate of 1 degree each ns. I always put the
reference at 0 degree
I have obtained different results:
Sometime when I see the file pullx.xvg it says that my COM is changing
but in reality it dose not.
When I see the trajectory with VMD I can see that the bulge is not
moving.
Some other time instead (changing randomly k), the bulge moves, but it dose a
rotation of 180 degrees with few steps and for the rest of the simulations stay
steady.
Some other time it continues 'to be pulled’ and break the structure.
Any suggestion?
This is the part for the pulling:
pull-coord1-geometry = dihedral
pull-coord1-groups = 1 2 3 4 5 6
pull-group1-pbcatom = 0
pull_group1_name = phospate_6
pull_group2_name = phospate_8
pull_group3_name = phospate_8
pull_group4_name = U7
pull_group5_name = U7
pull_group6_name = bulge
pull_coord1_rate = 1
pull_coord1_k = 100. (I modify this)
pull-coord1-start = yes
pull-coord1-init = -59.903
Thanks for the help!
Luca
--
Gromacs Users mailing list
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a
mail to gmx-users-requ...@gromacs.org.
