[gmx-users] LINCS WARNING relative constraint deviation
Ali Alizadeh
ali.alizadehmojarad at gmail.com
Mon Oct 1 13:10:26 CEST 2012
Dear all users
1- I am trying to simulate a system that contains 1182 methane
molecules(5910 atoms) and 1830 water molecules(5490), I
keep getting these errors.
2- Gromacs does not give me a .xtc file, Why?
My command lines for grompp and mdrun in terminal window:
grompp -f npt.mdp -c wwfull.top.gro -p topol.top -o topol.tpr -n ww.ndx
mdrun -s topol.tpr -o wwfull.trr -x wwfull.xtc
These are errors and my .mdp file:
...........................................
Steepest Descents:
Tolerance (Fmax) = 1.00000e-04
Number of steps = 1000
Step 50, time 0.01 (ps) LINCS WARNING
relative constraint deviation after LINCS:
rms 0.000589, max 0.014000 (between atoms 2636 and 2638)
bonds that rotated more than 30 degrees:
atom 1 atom 2 angle previous, current, constraint length
2636 2637 39.4 0.1090 0.1092 0.1090
Step 51, time 0.0102 (ps) LINCS WARNING
relative constraint deviation after LINCS:
rms 0.001802, max 0.029832 (between atoms 2026 and 2029)
bonds that rotated more than 30 degrees:
atom 1 atom 2 angle previous, current, constraint length
2331 2334 43.1 0.1090 0.1093 0.1090
2026 2028 49.6 0.1090 0.1095 0.1090
5071 5074 41.7 0.1090 0.1093 0.1090
4806 4808 47.0 0.1090 0.1095 0.1090
Step 52, time 0.0104 (ps) LINCS WARNING
relative constraint deviation after LINCS:
rms 0.003002, max 0.049492 (between atoms 156 and 159)
bonds that rotated more than 30 degrees:
atom 1 atom 2 angle previous, current, constraint length
471 474 51.7 0.1090 0.1096 0.1090
2331 2334 32.4 0.1093 0.1078 0.1090
5071 5074 31.3 0.1093 0.1081 0.1090
156 160 55.3 0.1090 0.1095 0.1090
916 919 41.0 0.1090 0.1093 0.1090
5886 5890 61.1 0.1090 0.1092 0.1090
Step 53, time 0.0106 (ps) LINCS WARNING
relative constraint deviation after LINCS:
rms 0.034606, max 0.843079 (between atoms 5886 and 5890)
bonds that rotated more than 30 degrees:
atom 1 atom 2 angle previous, current, constraint length
471 475 31.5 0.1049 0.1325 0.1090
471 474 92.9 0.1096 0.1519 0.1090
471 473 38.6 0.1044 0.1402 0.1090
471 472 92.7 0.1048 0.1568 0.1090
2331 2335 48.8 0.1082 0.1083 0.1090
2331 2334 48.4 0.1078 0.1070 0.1090
2331 2333 88.8 0.1110 0.1564 0.1090
2331 2332 88.8 0.1110 0.1586 0.1090
2026 2030 89.4 0.1103 0.1177 0.1090
2026 2028 40.0 0.1091 0.1191 0.1090
2026 2027 38.2 0.1089 0.1207 0.1090
5071 5075 41.9 0.1084 0.1048 0.1090
5071 5074 41.6 0.1081 0.1043 0.1090
5071 5073 88.6 0.1109 0.1376 0.1090
5071 5072 88.7 0.1110 0.1392 0.1090
4806 4810 89.3 0.1104 0.1176 0.1090
4806 4809 32.7 0.1104 0.0927 0.1090
4806 4808 35.2 0.1089 0.1140 0.1090
4806 4807 35.0 0.1088 0.1153 0.1090
156 160 92.7 0.1095 0.1324 0.1090
156 159 80.8 0.1036 0.0523 0.1090
156 158 90.0 0.1036 0.1338 0.1090
156 157 90.4 0.1037 0.1489 0.1090
2276 2280 38.3 0.1091 0.1114 0.1090
2276 2277 33.0 0.1089 0.1134 0.1090
2286 2290 37.8 0.1090 0.1102 0.1090
916 920 90.1 0.1068 0.1704 0.1090
916 919 45.9 0.1093 0.1237 0.1090
916 917 65.0 0.1067 0.1285 0.1090
1296 1300 31.2 0.1090 0.1094 0.1090
5886 5890 90.9 0.1092 0.2009 0.1090
5886 5889 47.9 0.1043 0.1326 0.1090
5886 5888 85.9 0.1045 0.1282 0.1090
5886 5887 91.4 0.1046 0.1208 0.1090
Back Off! I just backed up step53b_n0.pdb to ./#step53b_n0.pdb.1#
Back Off! I just backed up step53c_n0.pdb to ./#step53c_n0.pdb.1#
Wrote pdb files with previous and current coordinates
..
..
..
This is my .mdp file:
title = Yo
cpp = /usr/bin/cpp
include =
define =
; RUN CONTROL PARAMETERS
integrator = steep
; Start time and timestep in ps
tinit = 0
dt = 0.0002
nsteps = 1000
; For exact run continuation or redoing part of a run
init_step = 0
; mode for center of mass motion removal
comm_mode = Linear
; number of steps for center of mass motion removal
nstcomm = 10
; ENERGY MINIMIZATION OPTIONS
; Force tolerance and initial step-size
emtol = 0.0001
emstep = 0.00001
; Max number of iterations in relax_shells
niter = 100
; Step size (ps^2) for minimization of flexible constraints
fcstep = 0
; Frequency of steepest descents steps when doing CG
nstcgsteep = 1000
nbfgscorr = 10
; OUTPUT CONTROL OPTIONS
; Output frequency for coords (x), velocities (v) and forces (f)
nstxout = 1
nstvout = 100
nstfout = 100
nstcalcenergy = -1
; Checkpointing helps you continue after crashes
nstcheckpoint = 1000
; Output frequency for energies to log file and energy file
nstlog = 250
nstenergy = 60
; Output frequency and precision for xtc file
nstxtcout = 1
xtc_precision = 1000
; This selects the subset of atoms for the xtc file. You can
; select multiple groups. By default all atoms will be written.
xtc_grps =
; NEIGHBORSEARCHING PARAMETERS
; nblist update frequency
nstlist = 20
; ns algorithm (simple or grid)
ns_type = grid
; Periodic boundary conditions: xyz (default), no (vacuum)
; or full (infinite systems only)
pbc = xyz
; nblist cut-off
rlist = 0.35
rlistlong = 0.35 ; cut off distance for
short range neighbor list
; OPTIONS FOR ELECTROSTATICS AND VDW
; Method for doing electrostatics
coulombtype = PME
rcoulomb_switch = 0
rcoulomb = 0.35
; Method for doing Van der Waals
vdw_type = Cut-off
; cut-off lengths
rvdw_switch = 0
rvdw = 0.35
; Apply long range dispersion corrections for Energy and Pressure
DispCorr = EnerPres
; Extension of the potential lookup tables beyond the cut-off
table_extension = 1
; Spacing for the PME/PPPM FFT grid
fourierspacing = 0.1
; FFT grid size, when a value is 0 fourierspacing will be used
fourier_nx =
fourier_ny =
fourier_nz =
; EWALD/PME/PPPM parameters
pme_order = 4
ewald_rtol = 1e-05
ewald_geometry = 3d
epsilon_surface = 0
optimize_fft = no
; OPTIONS FOR WEAK COUPLING ALGORITHMS
; Temperature coupling
Tcoupl = berendsen
; Groups to couple separately
nsttcouple = -1
tc_grps = System
; Time constant (ps) and reference temperature (K)
tau_t = 1
ref_t = 240
; Pressure coupling
Pcoupl = berendsen
Pcoupltype = isotropic
nstpcouple = -1
; Time constant (ps), compressibility (1/bar) and reference P (bar)
tau_p = 5
compressibility = 5e-5
ref_p = 300
; Random seed for Andersen thermostat
andersen_seed = 815131
; GENERATE VELOCITIES FOR STARTUP RUN
gen_vel = no
gen_temp =
gen_seed =
; OPTIONS FOR BONDS
constraints = all-bonds
; Type of constraint algorithm
constraint_algorithm = Lincs
; Do not constrain the start configuration
continuation = no
; Use successive overrelaxation to reduce the number of shake iterations
Shake_SOR = no
; Relative tolerance of shake
shake_tol = 1e-04
; Highest order in the expansion of the constraint coupling matrix
lincs_order = 4
; Number of iterations in the final step of LINCS. 1 is fine for
; normal simulations, but use 2 to conserve energy in NVE runs.
; For energy minimization with constraints it should be 4 to 8.
lincs_iter = 1
; Lincs will write a warning to the stderr if in one step a bond
; rotates over more degrees than
lincs_warnangle = 30
; Convert harmonic bonds to morse potentials
morse = no
--
Sincerely
More information about the gromacs.org_gmx-users
mailing list