<!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">
    <blockquote
      cite="mid:BANLkTi=jvg-U-0vAkuFQPttR7JaZeQXbAw@mail.gmail.com"
      type="cite">Hello gmx users,<br>
      <br>
      My system for NVT equilbration runs into segmentation fault as
      soon as I try to run it.<br>
      It does not give any warning or hint of what might be going wrong.<br>
      Since I am a new user I am having difficulty exploring the
      plausible reasons.<br>
      <br>
      System: Protein( polyhistidine), 20&nbsp; 2,5-dihydroxybenzoic acid
      anions, 1:1 water: methanol (~3000 molecules of each) in 8 nm cube<br>
      <br>
      <br>
      I had had EM of the system using steepest decent. Outcome:<br>
      <br>
      Steepest Descents converged to machine precision in 15 steps, but
      did not reach the requested Fmax &lt; 1000.<br>
      Potential Energy&nbsp; =&nbsp; 1.5354981e+19<br>
      Maximum force&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inf on atom 651<br>
      Norm of force&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inf<br>
    </blockquote>
    <br>
    So that's broken already - you have enormous positive energy and
    infinite forces. Stop there and fix it. Either your starting
    configuration has severe atomic overlaps (go and visualize it with
    the periodic box), or some of your topology is broken (try a box of
    methanol on its own, try the protein in vacuum, try a single acid in
    vacuum)<br>
    <br>
    Mark<br>
    <br>
    <blockquote
      cite="mid:BANLkTi=jvg-U-0vAkuFQPttR7JaZeQXbAw@mail.gmail.com"
      type="cite"><b>The minim.mdp </b>is:<br>
      ; minim.mdp - used as input into grompp to generate em.tpr<br>
      ; Parameters describing what to do, when to stop and what to save<br>
      integrator&nbsp; = steep&nbsp;&nbsp;&nbsp;&nbsp; ; Algorithm (steep = steepest descent
      minimization)<br>
      emtol&nbsp;&nbsp;&nbsp; = 1000.0&nbsp;&nbsp;&nbsp; ; Stop minimization when the maximum force
      &lt; 1000.0 kJ/mol/nm<br>
      emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.02&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; Energy step size<br>
      nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 50000&nbsp;&nbsp;&nbsp;&nbsp; ; Maximum number of (minimization) steps
      to perform<br>
      <br>
      ; Parameters describing how to find the neighbors of each atom and
      how to calculate the interactions<br>
      nstlist&nbsp;&nbsp;&nbsp;&nbsp; = 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; Frequency to update the neighbor list and
      long range forces<br>
      ns_type&nbsp;&nbsp;&nbsp;&nbsp; = grid&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; Method to determine neighbor list
      (simple, grid)<br>
      rlist&nbsp;&nbsp;&nbsp; = 1.0&nbsp;&nbsp;&nbsp; ; Cut-off for making neighbor list (short range
      forces)<br>
      coulombtype = PME&nbsp;&nbsp;&nbsp; ; Treatment of long range electrostatic
      interactions<br>
      rcoulomb = 1.0&nbsp;&nbsp;&nbsp; ; Short-range electrostatic cut-off<br>
      rvdw&nbsp;&nbsp;&nbsp;&nbsp; = 1.0&nbsp;&nbsp;&nbsp; ; Short-range Van der Waals cut-off<br>
      pbc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = xyz&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; Periodic Boundary Conditions (yes/no)<br>
      constraints = none<br>
      <br>
      <br>
      <br>
      <br>
      <b>The nvt.mdp</b>:<br>
      <br>
      title&nbsp;&nbsp;&nbsp; = hist NVT equilibration<br>
      define&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -DPOSRES&nbsp; ; position restrain the protein<br>
      ; Run parameters<br>
      integrator&nbsp; = md&nbsp;&nbsp;&nbsp;&nbsp; ; leap-frog integrator<br>
      nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 50000&nbsp;&nbsp;&nbsp;&nbsp; ; 2 * 50000 = 100 ps<br>
      dt&nbsp;&nbsp;&nbsp; = 0.002 &nbsp;&nbsp;&nbsp; ; 2 fs<br>
      ; Output control<br>
      nstxout&nbsp;&nbsp;&nbsp;&nbsp; = 100&nbsp;&nbsp;&nbsp; ; save coordinates every 0.2 ps<br>
      nstvout&nbsp;&nbsp;&nbsp;&nbsp; = 100&nbsp;&nbsp;&nbsp; ; save velocities every 0.2 ps<br>
      nstenergy&nbsp;&nbsp; = 100&nbsp;&nbsp;&nbsp; ; save energies every 0.2 ps<br>
      nstlog&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 100&nbsp;&nbsp;&nbsp; ; update log file every 0.2 ps<br>
      ; Bond parameters<br>
      continuation&nbsp;&nbsp; = no&nbsp;&nbsp;&nbsp;&nbsp; ; first dynamics run<br>
      constraint_algorithm = lincs&nbsp; ; holonomic constraints<br>
      constraints = none ;<br>
      lincs_iter&nbsp; = 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; accuracy of LINCS<br>
      lincs_order = 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; also related to accuracy<br>
      ; Neighborsearching<br>
      ns_type&nbsp;&nbsp;&nbsp;&nbsp; = grid&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; search neighboring grid cells<br>
      nstlist&nbsp;&nbsp;&nbsp;&nbsp; = 5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; 10 fs<br>
      rlist&nbsp;&nbsp;&nbsp; = 1.0&nbsp;&nbsp;&nbsp; ; short-range neighborlist cutoff (in nm)<br>
      rcoulomb = 1.0&nbsp;&nbsp;&nbsp; ; short-range electrostatic cutoff (in nm)<br>
      rvdw&nbsp;&nbsp;&nbsp;&nbsp; = 1.0&nbsp;&nbsp;&nbsp; ; short-range van der Waals cutoff (in nm)<br>
      ; Electrostatics<br>
      coulombtype = PME&nbsp;&nbsp;&nbsp; ; Particle Mesh Ewald for long-range
      electrostatics<br>
      pme_order&nbsp;&nbsp; = 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; cubic interpolation<br>
      fourierspacing = 0.16&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; grid spacing for FFT<br>
      ; Temperature coupling is on<br>
      tcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = V-rescale ; modified Berendsen thermostat<br>
      tc-grps&nbsp;&nbsp;&nbsp;&nbsp; = Protein Non-Protein&nbsp;&nbsp; ; two coupling groups - more
      accurate<br>
      tau_t&nbsp;&nbsp;&nbsp; = 0.1 0.1&nbsp;&nbsp; ; time constant, in ps<br>
      ref_t&nbsp;&nbsp;&nbsp; = 300&nbsp;&nbsp;&nbsp; 300&nbsp;&nbsp; ; reference temperature, one for each
      group, in K<br>
      ; Pressure coupling is off<br>
      pcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = no&nbsp;&nbsp;&nbsp;&nbsp; ; no pressure coupling in NVT<br>
      ; Periodic boundary conditions<br>
      pbc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = xyz&nbsp;&nbsp;&nbsp; ; 3-D PBC<br>
      ; Dispersion correction<br>
      DispCorr = EnerPres&nbsp; ; account for cut-off vdW scheme<br>
      ; Velocity generation<br>
      gen_vel&nbsp;&nbsp;&nbsp;&nbsp; = yes&nbsp;&nbsp;&nbsp; ; assign velocities from Maxwell distribution<br>
      gen_temp = 300&nbsp;&nbsp;&nbsp; ; temperature for Maxwell distribution<br>
      gen_seed = -1&nbsp;&nbsp;&nbsp;&nbsp; ; generate a random seed<br>
      <br>
      <br>
      I tried to decrease the step size, that also runs into seg fault
      error.<br>
      <br>
      <br>
      Kindly guide.<br>
      <br>
      Thanks,<br>
      SN<br>
      <br>
      <br>
      <pre wrap=""></pre>
    </blockquote>
    <br>
  </body>
</html>