He, Yang wrote:
Hi all users,
When I use the command "mdrun -tablep table.xvg " to run my model in gromacs,
it always shows that:
starting mdrun 'DNA in water'
10000 steps, 20.0 ps.
This is orders of magnitude too short to expect a large conformational
change to happen. Think nanoseconds at least. Depending on the number of
base pairs, think much more than that. You should be aware of how long
this might take, since you've read some literature about MD on DNA, right?
Warning: 1-4 interaction between 11 and 19 at distance 1.437 which is larger
than the 1-4 table size 1.000 nm
These are ignored for the rest of the simulation
This usually means your system is exploding,
if not, you should increase table-extension in your mdp file
Don't ignore warnings unless you know what they mean and it's OK. See
http://wiki.gromacs.org/index.php/Errors#1-4_interaction_not_within_cut-off
And no more information . Furthermore, I use the command ngmx to look at the trajectory ,
I found that only when I set the value of nsteps by 10, will there be little animation
about my DNA two strands' disassociation . While I try to reduce the value of dt ,
nstdisreout , nstorireout and nstdihreout , it also did not work and I can not see any
disassociation of DNA strands using the command "ngmx".
Probably, you only wrote one frame before it crashed. Does the log file
look like it terminated normally after 10000 steps? Use gmxcheck to see
how many frames you have in the trajectory file.
Also, what I got confused is that when I set the value of nsteps larger than
10, there is just static figure about DNA model. This is my part of mdp file
1) You need the simulation not to crash
2) You need to simulate for long enough (nsteps) to see a change
3) You need to write output (nstxout or nstxtcout) sufficiently often to
keep yourself happy.
Mark
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
