<div dir="ltr">Hi,<br><br>I get the error:<br>Fatal error:<br>Step 0: The X-size
(0.000000) times the triclinic skew factor (1.000000) is smaller than
the smallest allowed cell size (1.200000) for domain decomposition grid
cell 0 0 0<br>



<br>I don&#39;t quite understand the initialization of cell_x1, cell_x0, thus have problems to find out the cause. <br><br>This is with RF and mdrun options -nosum -dlb<br><br>I
tested both revision 1.145 and 1.148 of domdec.c. (all other files
today version) Thus it is not the recent fatal-fix or 2d PME
preparation.<br>

<br>The simulation runs find without dlb.<br>
<br>Where should I start looking why the size is 0?<br><br>The below md.log file is for 1M atoms and 256 cores but it crashes the same with a subsystem of 4200 atoms and 16 cores.<br><br>Thanks<br>Roland<br><br>the md.log:<br>

<br><br>parameters of the run:<br>&nbsp;&nbsp; integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = md<br>&nbsp;&nbsp; nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>&nbsp;&nbsp; init_step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Grid<br>&nbsp;&nbsp; nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br>&nbsp;&nbsp; ndelta&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2<br>&nbsp;&nbsp; nstcomm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; comm_mode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Linear<br>&nbsp;&nbsp; nstcheckpoint&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>&nbsp;&nbsp; nstlog&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>



&nbsp;&nbsp; nstxout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; nstvout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; nstfout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; nstenergy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; nstxtcout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; init_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; delta_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.002<br>



&nbsp;&nbsp; xtcprec&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>&nbsp;&nbsp; nkx&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 72<br>&nbsp;&nbsp; nky&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 72<br>&nbsp;&nbsp; nkz&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 72<br>&nbsp;&nbsp; pme_order&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 4<br>&nbsp;&nbsp; ewald_rtol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1e-05<br>&nbsp;&nbsp; ewald_geometry&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>



&nbsp;&nbsp; epsilon_surface&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; optimize_fft&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = TRUE<br>&nbsp;&nbsp; ePBC&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = xyz<br>&nbsp;&nbsp; bPeriodicMols&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>&nbsp;&nbsp; bContinuation&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>&nbsp;&nbsp; bShakeSOR&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>&nbsp;&nbsp; etc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = No<br>



&nbsp;&nbsp; epc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = No<br>&nbsp;&nbsp; epctype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Isotropic<br>&nbsp;&nbsp; tau_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>&nbsp;&nbsp; ref_p (3x3):<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ref_p[&nbsp;&nbsp;&nbsp; 0]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ref_p[&nbsp;&nbsp;&nbsp; 1]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>



&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ref_p[&nbsp;&nbsp;&nbsp; 2]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>&nbsp;&nbsp; compress (3x3):<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; compress[&nbsp;&nbsp;&nbsp; 0]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; compress[&nbsp;&nbsp;&nbsp; 1]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>



&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; compress[&nbsp;&nbsp;&nbsp; 2]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>&nbsp;&nbsp; refcoord_scaling&nbsp;&nbsp;&nbsp;&nbsp; = No<br>&nbsp;&nbsp; posres_com (3):<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; posres_com[0]= 0.00000e+00<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; posres_com[1]= 0.00000e+00<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; posres_com[2]= 0.00000e+00<br>



&nbsp;&nbsp; posres_comB (3):<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; posres_comB[0]= 0.00000e+00<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; posres_comB[1]= 0.00000e+00<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; posres_comB[2]= 0.00000e+00<br>&nbsp;&nbsp; andersen_seed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 815131<br>&nbsp;&nbsp; rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.2<br>&nbsp;&nbsp; rtpi&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.05<br>



&nbsp;&nbsp; coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Reaction-Field<br>&nbsp;&nbsp; rcoulomb_switch&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.2<br>&nbsp;&nbsp; vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Cut-off<br>&nbsp;&nbsp; rvdw_switch&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.2<br>&nbsp;&nbsp; epsilon_r&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>



&nbsp;&nbsp; epsilon_rf&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 78<br>&nbsp;&nbsp; tabext&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>&nbsp;&nbsp; implicit_solvent&nbsp;&nbsp;&nbsp;&nbsp; = No<br>&nbsp;&nbsp; gb_algorithm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Still<br>&nbsp;&nbsp; gb_epsilon_solvent&nbsp;&nbsp; = 80<br>&nbsp;&nbsp; nstgbradii&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>&nbsp;&nbsp; rgbradii&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2<br>



&nbsp;&nbsp; gb_saltconc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; gb_obc_alpha&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>&nbsp;&nbsp; gb_obc_beta&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.8<br>&nbsp;&nbsp; gb_obc_gamma&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 4.85<br>&nbsp;&nbsp; sa_surface_tension&nbsp;&nbsp; = 2.092<br>&nbsp;&nbsp; DispCorr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = No<br>&nbsp;&nbsp; free_energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = no<br>



&nbsp;&nbsp; init_lambda&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; sc_alpha&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; sc_power&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; sc_sigma&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.3<br>&nbsp;&nbsp; delta_lambda&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; nwall&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; wall_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 9-3<br>



&nbsp;&nbsp; wall_atomtype[0]&nbsp;&nbsp;&nbsp;&nbsp; = -1<br>&nbsp;&nbsp; wall_atomtype[1]&nbsp;&nbsp;&nbsp;&nbsp; = -1<br>&nbsp;&nbsp; wall_density[0]&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; wall_density[1]&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; wall_ewald_zfac&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 3<br>&nbsp;&nbsp; pull&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = no<br>&nbsp;&nbsp; disre_weighting&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Conservative<br>



&nbsp;&nbsp; disre_mixed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>&nbsp;&nbsp; dr_fc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>&nbsp;&nbsp; dr_tau&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; nstdisreout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; orires_fc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; orires_tau&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; nstorireout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 100<br>



&nbsp;&nbsp; dihre-fc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>&nbsp;&nbsp; em_stepsize&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.01<br>&nbsp;&nbsp; em_tol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 100<br>&nbsp;&nbsp; niter&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 20<br>&nbsp;&nbsp; fc_stepsize&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; nstcgsteep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>&nbsp;&nbsp; nbfgscorr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br>



&nbsp;&nbsp; ConstAlg&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Lincs<br>&nbsp;&nbsp; shake_tol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1e-04<br>&nbsp;&nbsp; lincs_order&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 8<br>&nbsp;&nbsp; lincs_warnangle&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 30<br>&nbsp;&nbsp; lincs_iter&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>&nbsp;&nbsp; bd_fric&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; ld_seed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1993<br>



&nbsp;&nbsp; cos_accel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; deform (3x3):<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; deform[&nbsp;&nbsp;&nbsp; 0]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; deform[&nbsp;&nbsp;&nbsp; 1]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; deform[&nbsp;&nbsp;&nbsp; 2]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>



&nbsp;&nbsp; userint1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; userint2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; userint3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; userint4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; userreal1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; userreal2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; userreal3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; userreal4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>



grpopts:<br>&nbsp;&nbsp; nrdf:&nbsp;&nbsp;&nbsp; 2.16e+06<br>&nbsp;&nbsp; ref_t:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>&nbsp;&nbsp; tau_t:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>anneal:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; No<br>ann_npoints:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>&nbsp;&nbsp; acc:&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>&nbsp;&nbsp; nfreeze:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N<br>



&nbsp;&nbsp; energygrp_flags[&nbsp; 0]: 0 0<br>&nbsp;&nbsp; energygrp_flags[&nbsp; 1]: 0 0<br>&nbsp;&nbsp; efield-x:<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>&nbsp;&nbsp; efield-xt:<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>&nbsp;&nbsp; efield-y:<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>&nbsp;&nbsp; efield-yt:<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>&nbsp;&nbsp; efield-z:<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>



&nbsp;&nbsp; efield-zt:<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>&nbsp;&nbsp; bQMMM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>&nbsp;&nbsp; QMconstraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; QMMMscheme&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>&nbsp;&nbsp; scalefactor&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>qm_opts:<br>&nbsp;&nbsp; ngQM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br><br>Initializing Domain Decomposition on 256 nodes<br>



Will use dynamic load balancing<br>Will sort the charge groups at every domain (re)decomposition<br>Scaling the initial minimum size with 1/0.8 (option -dds) = 1.25<br>Optimizing the DD grid for 256 cells with a minimum initial size of 0.750 nm<br>



The maximum allowed number of cells is: X 24 Y 24 Z 37<br>Domain decomposition grid 8 x 4 x 8, separate PME nodes 0<br>Domain decomposition nodeid 0, coordinates 0 0 0<br><br>Using two step summing over 64 groups of on average 4.0 processes<br>



<br>There are 0 dihedral restraints<br>Table routines are used for coulomb: FALSE<br>Table routines are used for vdw:&nbsp;&nbsp;&nbsp;&nbsp; FALSE<br>Cut-off&#39;s:&nbsp;&nbsp; NS: 1.2&nbsp;&nbsp; Coulomb: 1.2&nbsp;&nbsp; LJ: 1.2<br>Reaction-Field:<br>epsRF = 78, rc = 1.2, krf = 0.283823, crf = 1.24204, epsfac = 138.935<br>



The electrostatics potential has its minimum at r = 1.20774<br>System total charge: -0.000<br>Generated table with 1100 data points for 1-4 COUL.<br>Tabscale = 500 points/nm<br>Generated table with 1100 data points for 1-4 LJ6.<br>



Tabscale = 500 points/nm<br>Generated table with 1100 data points for 1-4 LJ12.<br>Tabscale = 500 points/nm<br><br>Enabling SPC water optimization for 240000 molecules.<br><br>Configuring nonbonded kernels...<br>Testing x86_64 SSE support... present.<br>



<br><br>Removing pbc first time<br>Done rmpbc<br><br>Initializing LINear Constraint Solver<br><br>++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>B. Hess and H. Bekker and H. J. C. Berendsen and J. G. E. M. Fraaije<br>



LINCS: A Linear Constraint Solver for molecular simulations<br>J. Comp. Chem. 18 (1997) pp. 1463-1472<br>-------- -------- --- Thank You --- -------- --------<br><br>The number of constraints is 144000<br><br>++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>



S. Miyamoto and P. A. Kollman<br>SETTLE: An Analytical Version of the SHAKE and RATTLE Algorithms for Rigid<br>Water Models<br>J. Comp. Chem. 13 (1992) pp. 952-962<br>-------- -------- --- Thank You --- -------- --------<br>



<br><br>Linking all bonded interactions to atoms<br>There are 1092000 inter charge-group exclusions,<br>will use an extra communication step for exclusion forces for Reaction-Field<br><br>The maximum number of communication pulses is: X 1 Y 1 Z 1<br>



The minimum size for domain decomposition cells is 1.200 nm<br>The requested allowed shrink of DD cells (option -dds) is: 0.80<br>The allowed shrink of domain decomposition cells is: X 0.52 Y 0.26 Z 0.34<br><br>The maximum allowed distance for atoms involved in interactions is:<br>



&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; two-body bonded interactions&nbsp; (-rdd)&nbsp;&nbsp; 1.200 nm<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; multi-body bonded interactions&nbsp; (-rdd)&nbsp;&nbsp; 1.200 nm<br><br><br>Making 3D domain decomposition grid 8 x 4 x 8, home cell index 0 0 0<br><br><br>Will not sum the energies at every step,<br>



therefore the energy file does not contain exact averages and fluctuations.<br><br>WARNING:<br>Because of the no energy summing option setting nstcomm (was 0) to nstlist (10)<br><br>There are: 1008000 Atoms<br>Charge
group distribution at step 0: 1262 1274 1251 1282 1277 1242 1289 1260
1264 1271 1277 1274 1249 1285 1263 1263 1261 1274 1250 1281 1277 1240
1289 1259 1265 1271 1278 1275 1249 1287 1263 1264 1292 1267 1299 1275
1277 1301 1265 1284 1251 1283 1269 1272 1280 1270 1265 1289 1292 1267
1299 1275 1277 1301 1265 1284 1251 1283 1269 1272 1280 1270 1265 1289
1239 1263 1266 1234 1250 1265 1240 1257 1282 1252 1277 1264 1266 1281
1260 1255 1239 1263 1266 1234 1250 1265 1240 1257 1282 1252 1277 1264
1266 1281 1260 1255 1249 1279 1263 1254 1275 1255 1272 1272 1254 1252
1246 1245 1263 1236 1252 1260 1249 1279 1262 1254 1275 1254 1272 1271
1254 1252 1247 1245 1263 1237 1252 1261 1262 1274 1251 1282 1277 1242
1289 1260 1264 1271 1277 1274 1249 1285 1263 1263 1261 1274 1250 1281
1277 1240 1289 1259 1265 1271 1278 1275 1249 1287 1263 1264 1291 1267
1299 1274 1277 1301 1264 1284 1251 1283 1269 1272 1280 1270 1265 1289
1291 1267 1299 1274 1277 1301 1264 1284 1251 1283 1269 1272 1280 1270
1265 1289 1240 1263 1266 1235 1250 1265 1241 1257 1282 1252 1277 1264
1266 1281 1260 1255 1240 1263 1266 1235 1250 1265 1241 1257 1282 1252
1277 1264 1266 1281 1260 1255 1249 1279 1263 1254 1275 1255 1272 1272
1254 1252 1246 1245 1263 1236 1252 1260 1249 1279 1262 1254 1275 1254
1272 1271 1254 1252 1247 1245 1263 1237 1252 1261<br>


<br>------------------------------<div dir="ltr">-------------------------<br>Program mdrun, VERSION 3.3.99_development_20080718<br>Source code file: domdec.c, line: 2220<br><br>Fatal error:<br>Step
0: The X-size (0.000000) times the triclinic skew factor (1.000000) is
smaller than the smallest allowed cell size (1.200000) for domain
decomposition grid cell 0 0 0<br>


-------------------------------------------------------<br><br>Thanx for Using GROMACS - Have a Nice Day<br></div></div>