--------------------------- Missatge original ---------------------------- Assumpte: De: [EMAIL PROTECTED] Data: Tue, Novembre 25, 2008 9:52 am A: "Discussion list for GROMACS users" <gmx-users@gromacs.org> [EMAIL PROTECTED] --------------------------------------------------------------------------
OOOOps, I forgot the Subject in the mail! Justin A. Lemkul wrote: > > > [EMAIL PROTECTED] wrote: >> >> Hi everyone, >> >> I'm having problems running an md simulation (with gmx-4.0.2) using vsites >> and a time step of 4fs. I did generate the tpr with pdb2gmx. The problem >> seems to occur when the molecule crosses pbc, since the first sign of the >> simulation not working fine after ~3ns is >> >> " >> Warning: 1-4 interaction between 6 and 8 at distance 6.122 which is larger >> than the 1-4 table size 2.400 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 >> or with user tables increase the table size >> " >> > > This is a commonly reported problem. See here: > > http://wiki.gromacs.org/index.php/blowing_up > Thanks, I was aware of this entry in the wiki. > Also search the list archives; you will pull up several hundred posts describing the problems other users have faced and how they overcame them. > > Other things that would be helpful to know if you are still having problems: What does your system contain? Did it adequately minimize? What kind of equilibration procedure did you perform? It's dna, ions and water. After I placed the appropriate entries in the ddb and bonded itp, I was checking whether I could go for 4fs, and I prepared different runs with different time steps. The set-ups that made the simulation crash were the ones where the dna had the time/chance to cross pbc. If this did not occur, the simulation went fine (i.e. did not explode). Also, in the mean time I run a simulation in NVE with the same settings and I did no observe the crash after 5ns, despite having jumps over the pbc in the last few hundreds of ps. I energy minimized the system with steepest descend, repeated a couple of times. They went fine after some initial lincs warnings. I admit I did not equilibrate much, however the crash happens after more than 3ns... let's consider the first 2ns some sort of tempering. Although my nstxtcout was set low enough to resolve the frames just before the crash, I strongly suspect it has to do with crossing the pbc.... I will run a couple of simulations with output freq. 1 and double check. > > <snip> > >> Tcoupl = v-rescale >> tc_grps = system >> tau_t = 0.1 >> ref_t = 300.00 > > Maybe you want "Protein Non-Protein" instead of "system" here? But that's just general advice, without any knowledge of what's in your system. > Right, although I don't think this is the problem. Just wanted to see if the sim. would run fine for some ns before production runs. Just occured to me now, I will also try with the old berendsen, just to double check, since I did not have any problem in NVE. > -Justin Thanks for your time! Guillem > >> Pcoupl = Berendsen >> Pcoupltype = isotropic >> tau_p = 1.0 >> compressibility = 4.5e-5 >> ref_p = 1.0 >> constraints = all-bonds >> constraint-algorithm = Lincs >> lincs-order = 6 >> lincs-iter = 2 >> lincs-warnangle = 30 >> >> >> I've tried, perhaps naively, using both options for periodic_molecule with >> the same results. >> >> I guess it's hard to tell what did go wrong, but perhaps somebody has an >> idea.. I thought of submitting a bugzilla, but maybe there are some more >> tests I could do to pin point the problem before. Any ideas? >> >> >> All the best, >> >> >> Guillem >> >> Dr. Guillem Portella >> MMB - Institute for Research in Biomedicine >> Parc Cientific de Barcelona >> http://mmb.pcb.ub.es/~gportella >> >> >> >> >> >> >> >> >> >> _______________________________________________ >> 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 >> > -- Dr. Guillem Portella MMB - Institute for Research in Biomedicine Parc Cientific de Barcelona http://mmb.pcb.ub.es/~gportella _______________________________________________ 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