[gmx-users] mdrun Segmentation fault (after install make_hole_mdrun)

Minghu Song songm at rpi.edu
Sat May 31 03:45:01 CEST 2003


Dear gmx users,

Recently I am using the modified make_hole version mdrun (thanks Dr. 
Smith) to make the
hole inthe lipid bilayer for membrane proteins.  The process of making 
hole itself
(mdrun -hole -holep *.mdp .....) worked  smoothly.  But today I found 
that I cannot run the
normal mdrun (except em process) anymore . For instance, the following is
the  results of simple MD simulation for 216 SPC water mentioned in 
tutorial:
 
mdrun  -s water.tpr -o water.trr  -c water_out.gro -v water.log
..................................................................................
Back Off! I just backed up md.log to ./#md.log.6#
Getting Loaded...
Reading file water.tpr, VERSION 3.1.4 (single precision)
Reading file water.tpr, VERSION 3.1.4 (single precision)
Loaded with Money
Back Off! I just backed up ener.edr to ./#ener.edr.13#
starting mdrun 'Pure Water'
5000 steps,     10.0 ps.
Back Off! I just backed up water.trr to ./#water.trr.7#
Back Off! I just backed up insidesurf.pdb to ./#insidesurf.pdb.11#
Segmentation fault
-------------------------

It seems that the line ".... insidesurf.pdb to ./#insidesurf.pdb.11#" 
should not apprear because I didn't the option -hole.
I re-configuration and compile them again, but the error still exists. 
 This is how I re-compiled today:
------------------------
make distclean
./configure
make
make install
make links
------------------------
 
Also, I have another question about the results from mdrun (make_hole 
version).  In the mdp file, I set up NPT ensemble (pcoupletype = 
anisotropic).
There is no any error reported after the calculation, but when I check 
the log file I found that the presssure is quite different as the setup 
one (1 bar):
mdrun -v -hole -holep hp.mdp  -deffnm insert
input run.mdp file:
----------------------------------------------------------------------
.....
Pcoupl                      =  berendsen
Pcoupltype               =  anisotropic
tau_p                        =  5
compressibility       =  4.5e-5  4.5e-5  4.5e-5  0   0   0
ref_p                        =  1.0       1.0       1.0       0   0   0
-----------------------------------------------------------------------

output in log file
-----------------------------------------------------------------------------------------------------------
Initial step 0:
Energies (kJ/mol)
          Angle    Proper Dih. Ryckaert-Bell.  Improper Dih.          LJ-14
    3.16511e+03    4.01188e+03    4.21941e+03    1.88532e+02    2.33083e+03
     Coulomb-14        LJ (SR)   Coulomb (SR)   Coulomb (LR) Position Rest.
    1.47721e+04   -4.50831e+03   -3.00000e+05   -2.99713e+03    2.42487e-02
      Potential    Kinetic En.        Total Energy       Temperature 
 Pressure (bar)
   -2.78818e+05    4.09423e+04   -2.37875e+05    2.99016e+02   
-3.39316e+03  (1bar ?)
                                                                        
                                     
last step 10000:
   Energies (kJ/mol)
          Angle    Proper Dih. Ryckaert-Bell.  Improper Dih.          LJ-14
    1.20003e+04    3.73106e+03    5.73505e+03    3.63259e+02    3.46746e+03
     Coulomb-14        LJ (SR)   Coulomb (SR)   Coulomb (LR) Position Rest.
    1.00511e+04    5.84592e+03   -2.99810e+05   -1.21922e+04    5.76862e+02
      Potential    Kinetic En.   Total Energy    Temperature       
Pressure (bar)
   -2.70231e+05    4.16222e+04   -2.28609e+05    3.03982e+02    4.94660e+02
--------------------------------------------------------------------------------------------------------
 
I just began to learn the MD as well as gromacs.  Maybe I made some 
"silly" mistakes duringing the installation step
or ask a very simple question.  Anyone help is greatly appreciated. 
Thanks inadvance.




More information about the gromacs.org_gmx-users mailing list