Dear Michael: I have uploaded them as http://redmine.gromacs.org/issues/1306
It does not crash immediately. The crash is stochastic, giving a segfault between 200 and 5000 integration steps. That made me think it was a simple exploding system problem, but there are other things (listed in my original post) that make me think otherwise. Most notably, the drug is fine in both water and vacuum. I have also built numerous systems and get the crash in each one. My actual production system also has a protein, but during debugging I found that the error persists in a simple and small water solution. Thank you for your assistance. Chris. -- original message -- Chris, can you post a redmine on this so I can look at the files? Also, does it crash immediately, or after a while? On Thu, Jul 18, 2013 at 2:45 PM, Christopher Neale <chris.neale at mail.utoronto.ca> wrote: > Dear Users: > > I have a system with water and a drug (54 total atoms; 27 heavy atoms). The > system is stable when I simulate it for 1 ns. However, Once I add the > following options to the .mdp file, the run dies after a few ps with a > segfault. > > free-energy = yes > init-lambda = 1 > couple-lambda0 = vdw-q > couple-lambda1 = none > couple-intramol = no > couple-moltype = drug > > I do not get any step files or any lincs warnings. If I look at the .xtc and > .edr files, there is no indication of something blowing up before the > segfault. I have also verified that the drug runs without any problem in > vacuum. I get the same behaviour if I remove constraints and use a timestep > of 0.5 fs. The segfault is reproducible with v4.6.1 and v4.6.3. I am using > the charmm FF, but I converted all UB angles in my drug to type-1 angles and > still got the segfault. I also get the segfault with particle decomposition > and/or while running a single thread. I am currently using the SD integrator, > but I get the same segfault with md and md-vv. Couple-intramol=yes doesn't > resolve it, neither does using separate T-coupling groups for the water and > drug. Neither does turning off pressure coupling. > > Here is the .mdp file that works fine, but gives me a segfault when I add the > free energy stuff (above): > > constraints = all-bonds > lincs-iter = 1 > lincs-order = 6 > constraint_algorithm = lincs > integrator = sd > dt = 0.002 > tinit = 0 > nsteps = 100000 > nstcomm = 1 > nstxout = 0 > nstvout = 0 > nstfout = 0 > nstxtcout = 500 > nstenergy = 500 > nstlist = 10 > nstlog=0 ; reduce log file size > ns_type = grid > vdwtype = cut-off > rlist = 0.8 > rvdw = 0.8 > rcoulomb = 0.8 > coulombtype = cut-off > tc_grps = System > tau_t = 1.0 > ld_seed = -1 > ref_t = 310 > gen_temp = 310 > gen_vel = yes > unconstrained_start = no > gen_seed = -1 > Pcoupl = berendsen > pcoupltype = isotropic > tau_p = 4 > compressibility = 4.5e-5 > ref_p = 1.0 > > I do realize that some of these settings are not ideal for a production run. > I started with the real Charmm cutoffs + PME, etc, (which also gives the > segfault) but this is what I am using right now for quick testing. > > The only thing keeping me from filing a redmine issue is that if I remove my > drug and do the FEP on one of the water molecules (using the FEP code listed > above), I have no segfault. Therefore it is clearly related to the drug, > whose parameters I built so I may have caused the problem somehow. > Nevertheless, the drug runs fine in water and in vacuum without the FEP code, > so I can't imagine what could be causing this segfault (also, the fact that > it's a segfault means that I don;t get any useful info from mdrun as to what > might be going wrong). > > Thank you, > Chris. -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! * Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists