[gmx-users] so difficult problem
Amit Choubey
kgp.amit at gmail.com
Tue Jan 4 10:08:44 CET 2011
There is something wrong with your initial configuration. May be you forgot
to take care of periodicity, how did you get your initial configuration?
Also notice that these kind of problems have been discussed previously.
Amit
On Tue, Jan 4, 2011 at 12:33 AM, mohsen ramezanpour <
ramezanpour.mohsen at gmail.com> wrote:
> Dear All
> I am using this .mdp file and I recived the below warnings,I can't solve
> that.
>
>
> title = NPT
> define =
> integrator = md
> tinit = 0
> dt = 0.002
> nsteps = 500000
> nstcomm = 1
> comm-grps = protein non-protein
> niter = 20
> nstxout = 5000
> nstvout = 5000
> nstfout = 0
> nstlog = 5000
> nstenergy = 250
> nstxtcout = 250
> xtc-precision = 1000
> xtc_grps = protein non-protein
> energygrps = Protein non-protein
> nstlist = 5
> ns_type = grid
> pbc = xyz
> rlist = 1.4
> domain-decomposition = no
> coulombtype = PME
> rcoulomb-switch = 0
> rcoulomb = 1.4
> epsilon-r = 1
> vdw-type = Cut-off
> rvdw-switch = 0
> rvdw = 1.4
> DispCorr = EnerPres
> fourierspacing = 0.12
> fourier_nx = 0
> fourier_ny = 0
> fourier_nz = 0
> pme_order = 4
> ewald_rtol = 1e-05
> epsilon_surface = 0
> optimize_fft = no
> tcoupl = Nose-hoover
> tc-grps = Protein non-protein
> tau_t = 0.1 0.1
> ref_t = 300 300
> Pcoupl = Parrinello-Rahman
> Pcoupltype = Isotropic
> tau_p = 1.0
> compressibility = 4.5e-5
> ref_p = 1.0
> annealing = no
> gen_vel = yes
> gen_temp = 310
> gen_seed = 173529
> constraints = all-bonds
> constraint-algorithm = Lincs
> unconstrained-start = no
> lincs-order = 4
> lincs-warnangle = 30
> morse = no
>
> my sytem is protein-ligand,I want to generate a NPT.
> the result was:
>
> Step 0, time 0 (ps) LINCS WARNING
> relative constraint deviation after LINCS:
> rms 0.002537, max 0.119994 (between atoms 5293 and 5294)
> bonds that rotated more than 30 degrees:
> atom 1 atom 2 angle previous, current, constraint length
> starting mdrun 'Protein in water'
> 500000 steps, 1000.0 ps.
>
> Step 0, time 0 (ps) LINCS WARNING
> relative constraint deviation after LINCS:
> rms 1503.633433, max 55362.878906 (between atoms 5289 and 5290)
> bonds that rotated more than 30 degrees:
> atom 1 atom 2 angle previous, current, constraint length
> 5264 5263 41.7 0.1143 0.1261 0.1140
> 5293 5294 82.7 0.1221 2612.9744 0.1090
> 5291 5293 77.2 0.1517 6082.6147 0.1390
> 5291 5292 101.0 0.1117 1518.6106 0.1090
> 5289 5291 87.7 0.1353 6891.8911 0.1390
> 5289 5290 88.4 0.1437 7529.4878 0.1360
> 5289 5287 88.4 0.1456 7540.4873 0.1390
> 5287 5288 89.7 0.1092 381.5114 0.1090
> 5285 5287 92.3 0.1395 354.9220 0.1390
> 5285 5286 123.8 0.1101 44.8448 0.1090
> 5284 5293 89.0 0.1478 4605.1763 0.1390
> 5284 5285 79.3 0.1402 173.5627 0.1390
> 5278 5279 78.9 0.1529 1.1270 0.1530
> 5277 5278 122.9 0.1545 13.1423 0.1530
> 5276 5284 84.9 0.1422 159.0098 0.1390
> 5276 5277 68.4 0.1543 90.8248 0.1530
> 5276 5275 39.5 0.1439 72.5216 0.1430
> 5274 5275 103.7 0.1436 21.6028 0.1430
> 5273 5274 108.5 0.1394 3.3421 0.1390
> 5276 5272 45.7 0.1402 80.6795 0.1390
> 5272 5273 127.3 0.1340 14.2401 0.1330
> 5272 5270 69.8 0.1336 12.6250 0.1330
> 5270 5271 113.3 0.1092 0.2530 0.1090
> 5268 5270 104.6 0.1391 0.2718 0.1390
> 5268 5269 44.4 0.1091 0.1636 0.1090
> 5273 5266 58.5 0.1342 14.0125 0.1330
> 5265 5268 36.5 0.1394 0.1769 0.1390
>
> Back Off! I just backed up step0b.pdb to ./#step0b.pdb.2#
>
> Back Off! I just backed up step0c.pdb to ./#step0c.pdb.2#
> Wrote pdb files with previous and current coordinates
> step 0Warning: 1-4 interaction between 5264 and 5273 at distance 13.965
> 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
>
> Please let me know the solution.
> Thanks in advance
>
>
>
>
>
>
> --
> gmx-users mailing list gmx-users at 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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110104/edb9d95a/attachment.html>
More information about the gromacs.org_gmx-users
mailing list