<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body text="#000000" bgcolor="#ffffff">
    On 13/04/2011 2:10 AM, Anna Duncan wrote:
    <blockquote
      cite="mid:4ABEAAE5-3AD3-4AA5-888A-8EF4A7FDC005@mrc-mbu.cam.ac.uk"
      type="cite">
      <div>Hi,</div>
      <div><br>
      </div>
      <div>I've started using version 4.5.4 of gromacs, having
        previously been using version 4.0 on an older server.</div>
      <div><br>
      </div>
      <div>When I run a simulation on v4.0 everything seems to run fine.
        &nbsp;However, when I run the same simulation on v4.5.4, mdrun gives
        me the error message:</div>
      <div><br>
      </div>
      <div>########################################################################</div>
      <div><br>
      </div>
      <div>
        <div>Making 3D domain decomposition 4 x 3 x 2</div>
        <div>starting mdrun 'PROTEIN IN BILAYER'</div>
        <div>5000000 steps, 150000.0 ps.</div>
        <div>step 0</div>
        <div>Step 11 &nbsp;Warning: pressure scaling more than 1%, mu:
          1.02292 1.02292 1.0122</div>
      </div>
      <div><br>
      </div>
      <div>
        <div>Step 11 &nbsp;Warning: pressure scaling more than 1%, mu:
          1.02292 1.02292 1.0122</div>
      </div>
      <div><br>
      </div>
      <div>
        <div>.</div>
        <div>.</div>
        <div>.</div>
        <div>.</div>
        <div>Step 11 &nbsp;Warning: pressure scaling more than 1%, mu:
          1.02292 1.02292 1.0122</div>
        <div><br>
        </div>
        <div>Step 21 &nbsp;Warning: pressure scaling more than 1%, mu:
          1.09061 1.09061 1.02015</div>
        <div><br>
        </div>
        <div>Step 21 &nbsp;Warning: pressure scaling more than 1%, mu:
          1.09061 1.09061 1.02015</div>
        <div>.</div>
      </div>
      <div>.</div>
      <div>.</div>
      <div>
        <div>Step 21 &nbsp;Warning: pressure scaling more than 1%, mu:
          1.09061 1.09061 1.02015</div>
        <div><br>
        </div>
        <div>Step 25, time 0.75 (ps) &nbsp;LINCS WARNING</div>
        <div>relative constraint deviation after LINCS:</div>
        <div>rms 0.000004, max 0.000015 (between atoms 147 and 148)</div>
        <div>bonds that rotated more than 30 degrees:</div>
        <div>&nbsp;atom 1 atom 2 &nbsp;angle &nbsp;previous, current, constraint length</div>
        <div>&nbsp;&nbsp; &nbsp;178 &nbsp; &nbsp;179 &nbsp; 34.7 &nbsp; &nbsp;0.2600 &nbsp; 0.2600 &nbsp; &nbsp; &nbsp;0.2600</div>
        <div><br>
        </div>
        <div>Step 28, time 0.84 (ps) &nbsp;LINCS WARNING</div>
        <div>relative constraint deviation after LINCS:</div>
        <div>rms 19.241184, max 126.849724 (between atoms 35 and 36)</div>
        <div>bonds that rotated more than 30 degrees:</div>
        <div>&nbsp;atom 1 atom 2 &nbsp;angle &nbsp;previous, current, constraint length</div>
        <div><br>
        </div>
        <div>Step 28, time 0.84 (ps) &nbsp;LINCS WARNING</div>
        <div>relative constraint deviation after LINCS:</div>
        <div>rms 1.739422, max 8.610986 (between atoms 496 and 497)</div>
        <div>bonds that rotated more than 30 degrees:</div>
        <div>&nbsp;atom 1 atom 2 &nbsp;angle &nbsp;previous, current, constraint length</div>
        <div>&nbsp;&nbsp; &nbsp; 45 &nbsp; &nbsp; 46 &nbsp; 90.0 &nbsp; &nbsp;0.2650 &nbsp; 0.5083 &nbsp; &nbsp; &nbsp;0.2650</div>
        <div>&nbsp;&nbsp; &nbsp; 35 &nbsp; &nbsp; 36 &nbsp; 90.0 &nbsp; &nbsp;0.3100 &nbsp;39.6334 &nbsp; &nbsp; &nbsp;0.3100</div>
        <div>&nbsp;&nbsp; &nbsp; 42 &nbsp; &nbsp; 43 &nbsp; 90.0 &nbsp; &nbsp;0.2600 &nbsp; 1.2129 &nbsp; &nbsp; &nbsp;0.2600</div>
        <div>&nbsp;&nbsp; &nbsp; 45 &nbsp; &nbsp; 46 &nbsp; 90.0 &nbsp; &nbsp;0.2650 &nbsp; 0.5083 &nbsp; &nbsp; &nbsp;0.2650</div>
        <div>&nbsp;&nbsp; &nbsp;482 &nbsp; &nbsp;483 &nbsp; 90.0 &nbsp; &nbsp;0.2600 &nbsp; 0.8649 &nbsp; &nbsp; &nbsp;0.2600</div>
        <div>&nbsp;&nbsp; &nbsp;496 &nbsp; &nbsp;497 &nbsp; 90.0 &nbsp; &nbsp;0.2700 &nbsp; 2.5950 &nbsp; &nbsp; &nbsp;0.2700</div>
        <div>&nbsp;&nbsp; &nbsp;496 &nbsp; &nbsp;498 &nbsp; 90.0 &nbsp; &nbsp;0.2700 &nbsp; 1.4511 &nbsp; &nbsp; &nbsp;0.2700</div>
        <div>&nbsp;&nbsp; &nbsp;497 &nbsp; &nbsp;498 &nbsp; 90.0 &nbsp; &nbsp;0.2700 &nbsp; 2.5822 &nbsp; &nbsp; &nbsp;0.2700</div>
        <div>Wrote pdb files with previous and current coordinates</div>
        <div>Wrote pdb files with previous and current coordinates</div>
        <div>Segmentation fault</div>
      </div>
    </blockquote>
    <br>
    This looks like a standard case of <a
      href="http://www.gromacs.org/Documentation/Terminology/Blowing_Up">http://www.gromacs.org/Documentation/Terminology/Blowing_Up</a>,
    which can be quite situation dependent. In the absence of further
    information, I'd guess that 4.0 got lucky with the numerical
    integration, and 4.5.4 didn't.<br>
    <br>
    <blockquote
      cite="mid:4ABEAAE5-3AD3-4AA5-888A-8EF4A7FDC005@mrc-mbu.cam.ac.uk"
      type="cite">
      <div>
        <div>########################################################################</div>
      </div>
      <div><br>
      </div>
      <div><br>
      </div>
      <div>I've run g_gmxdump/gmxdump to get the input parameters for
        the simulations using the 2 different versions, </div>
    </blockquote>
    <br>
    gmxcheck is a good tool for comparing various file types.<br>
    <br>
    <blockquote
      cite="mid:4ABEAAE5-3AD3-4AA5-888A-8EF4A7FDC005@mrc-mbu.cam.ac.uk"
      type="cite">
      <div>which I &nbsp;can't attach because they're too big but
        which&nbsp;differ in the lines shown below (v4.0 on the left and
        v4.5.4 on the right):</div>
      <div><br>
      </div>
      <div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nstcomm &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</span></font></font><span
            class="Apple-tab-span" style="white-space: pre;"><font
              class="Apple-style-span" face="'Courier New'"><font
                class="Apple-style-span" size="3"><span
                  class="Apple-style-span" style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;|</span></font></font><span class="Apple-tab-span"
            style="white-space: pre;"><font class="Apple-style-span"
              face="'Courier New'"><font class="Apple-style-span"
                size="3"><span class="Apple-style-span"
                  style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nstcomm &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 10</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nstcalcenergy &nbsp; &nbsp; &nbsp; &nbsp;= 1</span></font></font><span
            class="Apple-tab-span" style="white-space: pre;"><font
              class="Apple-style-span" face="'Courier New'"><font
                class="Apple-style-span" size="3"><span
                  class="Apple-style-span" style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;|</span></font></font><span class="Apple-tab-span"
            style="white-space: pre;"><font class="Apple-style-span"
              face="'Courier New'"><font class="Apple-style-span"
                size="3"><span class="Apple-style-span"
                  style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nstcalcenergy &nbsp; &nbsp; &nbsp; &nbsp;= 10</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nsttcouple &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1</span></font></font><span
            class="Apple-tab-span" style="white-space: pre;"><font
              class="Apple-style-span" face="'Courier New'"><font
                class="Apple-style-span" size="3"><span
                  class="Apple-style-span" style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;|</span></font></font><span class="Apple-tab-span"
            style="white-space: pre;"><font class="Apple-style-span"
              face="'Courier New'"><font class="Apple-style-span"
                size="3"><span class="Apple-style-span"
                  style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nsttcouple &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nstpcouple &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1</span></font></font><span
            class="Apple-tab-span" style="white-space: pre;"><font
              class="Apple-style-span" face="'Courier New'"><font
                class="Apple-style-span" size="3"><span
                  class="Apple-style-span" style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;|</span></font></font><span class="Apple-tab-span"
            style="white-space: pre;"><font class="Apple-style-span"
              face="'Courier New'"><font class="Apple-style-span"
                size="3"><span class="Apple-style-span"
                  style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nstpcouple &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10</span></font></font></div>
      </div>
    </blockquote>
    <br>
    The heuristics by which these are set by default have changed. These
    values are not of themselves indicative of a problem.<br>
    <br>
    <blockquote
      cite="mid:4ABEAAE5-3AD3-4AA5-888A-8EF4A7FDC005@mrc-mbu.cam.ac.uk"
      type="cite">
      <div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                rgbradii &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 2</span></font></font><span
            class="Apple-tab-span" style="white-space: pre;"><font
              class="Apple-style-span" face="'Courier New'"><font
                class="Apple-style-span" size="3"><span
                  class="Apple-style-span" style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;|</span></font></font><span class="Apple-tab-span"
            style="white-space: pre;"><font class="Apple-style-span"
              face="'Courier New'"><font class="Apple-style-span"
                size="3"><span class="Apple-style-span"
                  style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                rgbradii &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                sa_surface_tension &nbsp; = 2.092</span></font></font><span
            class="Apple-tab-span" style="white-space: pre;"><font
              class="Apple-style-span" face="'Courier New'"><font
                class="Apple-style-span" size="3"><span
                  class="Apple-style-span" style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;|</span></font></font><span class="Apple-tab-span"
            style="white-space: pre;"><font class="Apple-style-span"
              face="'Courier New'"><font class="Apple-style-span"
                size="3"><span class="Apple-style-span"
                  style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                sa_surface_tension &nbsp; = 2.05016</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                sc_sigma_min &nbsp; &nbsp; &nbsp; &nbsp; = 0</span></font></font><span
            class="Apple-tab-span" style="white-space: pre;"><font
              class="Apple-style-span" face="'Courier New'"><font
                class="Apple-style-span" size="3"><span
                  class="Apple-style-span" style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;|</span></font></font><span class="Apple-tab-span"
            style="white-space: pre;"><font class="Apple-style-span"
              face="'Courier New'"><font class="Apple-style-span"
                size="3"><span class="Apple-style-span"
                  style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                sc_sigma_min &nbsp; &nbsp; &nbsp; &nbsp; = 0.3</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nstdhdl &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</span></font></font><span
            class="Apple-tab-span" style="white-space: pre;"><font
              class="Apple-style-span" face="'Courier New'"><font
                class="Apple-style-span" size="3"><span
                  class="Apple-style-span" style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;|</span></font></font><span class="Apple-tab-span"
            style="white-space: pre;"><font class="Apple-style-span"
              face="'Courier New'"><font class="Apple-style-span"
                size="3"><span class="Apple-style-span"
                  style="font-size: 11px;"> </span></font></font></span><font
            class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp;
                nstdhdl &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 10</span></font></font></div>
      </div>
    </blockquote>
    <br>
    These only matter if you're using the relevant algorithms, for which
    the defaults have apparently changed.<br>
    <br>
    <blockquote
      cite="mid:4ABEAAE5-3AD3-4AA5-888A-8EF4A7FDC005@mrc-mbu.cam.ac.uk"
      type="cite">
      <div>and also for cmap atomtypes values:</div>
      <div>in v4.0:</div>
      <div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;0]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;1]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;2]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;3]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;4]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;5]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;6]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;7]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;8]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;9]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 10]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 11]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 12]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 13]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 14]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 15]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 16]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 17]={radius= 0.00000e+00, volume=
                0.00000e+00, gb_radius= 0.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1, S_hct=
                0.00000e+00)}</span></font></font></div>
      </div>
      <div><br>
      </div>
      <div>and in v4.1</div>
      <div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;0]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;1]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;2]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;3]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;4]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;5]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;6]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;7]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;8]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ &nbsp;9]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 10]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 11]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 12]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 13]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 14]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 15]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 16]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
        <div><font class="Apple-style-span" face="'Courier New'"><font
              class="Apple-style-span" size="3"><span
                class="Apple-style-span" style="font-size: 11px;">&nbsp;&nbsp; &nbsp;
                &nbsp;atomtype[ 17]={radius=-1.00000e+00,
                volume=-1.00000e+00, gb_radius=-1.00000e+00,
                surftens=-1.00000e+00, atomnumber= &nbsp;-1,
                S_hct=-1.00000e+00)}</span></font></font></div>
      </div>
    </blockquote>
    <br>
    Since I don't think any native 4.0 forcefield used CMAP at all, I
    infer you aren't, so this doesn't matter either.<br>
    <br>
    <blockquote
      cite="mid:4ABEAAE5-3AD3-4AA5-888A-8EF4A7FDC005@mrc-mbu.cam.ac.uk"
      type="cite">
      <div>I haven't set any of these values and I don't understand why
        the default values for parameters such as nsttcouple differ by a
        factor of 10. &nbsp;I know that in v4.5.4 parameters such as
        nsttcouple are by default set to be equal to nstlist (which =10
        in my simulations) but I don't understand how the default is
        being set in v4.0. &nbsp;</div>
      <div><br>
      </div>
      <div>The larger nstpcouple default values seem to me like the
        cause of my warning messages and eventual crash of the system
        when run in v4.5.4 but I am also worried that I don't know how
        the cmap atomtypes values are being set and how these
        differences will affect the simulations. &nbsp;Do you have any
        pointers?</div>
    </blockquote>
    <br>
    We don't know whether you're even using T-coupling, P-coupling or
    CMAP. How about a quick description of your objective, simulation
    system, equilibration protocol and maybe an .mdp file?<br>
    <br>
    In any case, until you've exhausted the diagnosis steps on the above
    link, the smart money is betting that you have a somewhat unhappy
    starting configuration.<br>
    <br>
    Mark<br>
  </body>
</html>