AMBER Archive (2007)Subject: Re: AMBER: Minimization amber9 segmentation fault
From: Adrian Roitberg (roitberg_at_qtp.ufl.edu)
Date: Sat Dec 01 2007 - 10:17:15 CST
Francesco,
Seems like a straight vdw clash.
Strat by loading the system onto vmd and look around atom "O
17223" and
see if it is too close to someone else.
The gradient there is HUGE, so I suspect that. There might be other
clashes after you look at that one
NSTEP ENERGY RMS GMAX NAME NUMBER
1 2.8852E+18 1.4602E+18 3.9863E+20 O 17223
BOND = 149932.7945 ANGLE = 46211.6820 DIHED =
9161.8416
VDWAALS = ************* EEL = -246440.0364 HBOND =
0.0000
1-4 VDW = 29126.1421 1-4 EEL = 23905.1720 RESTRAINT =
0.0000
Then, after ONLY 12 steps minimization , you get
NSTEP ENERGY RMS GMAX NAME NUMBER
> 12 2.8815E+11 3.6643E+10 1.0002E+13 O 17223
>
> BOND = 149938.5018 ANGLE = 46211.8441 DIHED =
9161.8506
> VDWAALS = ************* EEL = -246163.7303 HBOND =
0.0000
> 1-4 VDW = 29126.0930 1-4 EEL = 23905.1432 RESTRAINT =
0.0000
> EAMBER = *************
So, your total energy decreased by 7 orders of magnitude and gradient by
eight orders.
I suspect the segmentation fault is due to numerical errors of handling
these large numbers, but cannot say for sure.
Adrian
Francesco Pietra wrote:
> Does this huge energy in trial minimization mean steric clashes or simply a bad
> min.in? I am posting as I can't understand better the segmentation fault.
>
> The system is a protein with non-covalently-bound large non-polymeric ligand,
> immersed in a POPC hydrated membrane, the whole hydrated (all water is TIP3P).
>
>>From ambmask:
> 1-76 POPC
> 77-520 protein
> 521 ligand
> 522-9999 WAT
> **** WAT
>
> (the latter corresponds to 0-4759 WAT from pdb obtained with ambpdb, otherwise
> residue number and name are the same).
>
> If relevant, should the last line of min.in end with a comma, as I did? There
> are different versions to this regard on the web.
>
> Thanks for help
> francesco pietra
>
> -------------------------------------------------------
> Amber 9 SANDER 2006
> -------------------------------------------------------
>
> | Run on 12/01/2007 at 16:29:59
> [-O]verwriting output
>
> File Assignments:
> | MDIN: box_min1.in
> | MDOUT: box_min1.out
> |INPCRD: protein_ligand_pop_box.inpcrd
>
> | PARM: protein_ligand_pop_box.prmtop
>
> |RESTRT: box_min1.rst
> | REFC: protein_gam_pop_box.inpcrd
>
> | MDVEL: mdvel
> | MDEN: mden
> | MDCRD: mdcrd
> |MDINFO: mdinfo
> |INPDIP: inpdip
> |RSTDIP: rstdip
>
>
> Here is the input file:
>
> 1st minimization protein_ligand_pop_box80x80
>
> &cntrl
> imin=1, maxcyc=20, ncyc=10,
> cut=10, ntb=1, ntpr=1,
> ntr=1, restraintmask=":77-520, 521",
> restraint_wt=30,
> /
>
>
>
> --------------------------------------------------------------------------------
> 1. RESOURCE USE:
> --------------------------------------------------------------------------------
>
> | Flags: MPI
>
> getting new box info from bottom of inpcrd
> | INFO: Old style inpcrd file read
>
> | peek_ewald_inpcrd: Box info found
> |Largest sphere to fit in unit cell has radius = 46.090
> | New format PARM file being parsed.
> | Version = 1.000 Date = 11/28/07 Time = 12:12:33
> NATOM = 89900 NTYPES = 18 NBONH = 82478 MBONA = 7441
> NTHETH = 22945 MTHETA = 9352 NPHIH = 38729 MPHIA = 20393
> NHPARM = 0 NPARM = 0 NNB = 192544 NRES = 24759
> NBONA = 7441 NTHETA = 9352 NPHIA = 20393 NUMBND = 62
> NUMANG = 130 NPTRA = 58 NATYP = 44 NPHB = 1
> IFBOX = 1 NMXRS = 134 IFCAP = 0 NEXTRA = 0
> NCOPY = 0
>
>
> | Memory Use Allocated
> | Real 6108229
> | Hollerith 564161
> | Integer 3173435
> | Max Pairs 12945600
> | nblistReal 1078800
> | nblist Int 3319286
> | Total 134283 kbytes
> | Duplicated 0 dihedrals
> | Duplicated 0 dihedrals
>
> BOX TYPE: RECTILINEAR
>
> --------------------------------------------------------------------------------
> 2. CONTROL DATA FOR THE RUN
> --------------------------------------------------------------------------------
>
>
>
>
> General flags:
> imin = 1, nmropt = 0
>
> Nature and format of input:
> ntx = 1, irest = 0, ntrx = 1
>
> Nature and format of output:
> ntxo = 1, ntpr = 1, ntrx = 1, ntwr = 500
> iwrap = 0, ntwx = 0, ntwv = 0, ntwe = 0
> ioutfm = 0, ntwprt = 0, idecomp = 0, rbornstat= 0
>
> Potential function:
> ntf = 1, ntb = 1, igb = 0, nsnb = 25
> ipol = 0, gbsa = 0, iesp = 0
> dielc = 1.00000, cut = 10.00000, intdiel = 1.00000
> scnb = 2.00000, scee = 1.20000
>
> Frozen or restrained atoms:
> ibelly = 0, ntr = 1
>
> Energy minimization:
> maxcyc = 20, ncyc = 10, ntmin = 1
> dx0 = 0.01000, drms = 0.00010
>
> Ewald parameters:
> verbose = 0, ew_type = 0, nbflag = 1, use_pme = 1
> vdwmeth = 1, eedmeth = 1, netfrc = 0
> Box X = 112.533 Box Y = 109.793 Box Z = 92.179
> Alpha = 90.000 Beta = 90.000 Gamma = 90.000
> NFFT1 = 120 NFFT2 = 120 NFFT3 = 96
> Cutoff= 10.000 Tol =0.100E-04
> Ewald Coefficient = 0.27511
> Interpolation order = 4
>
> LOADING THE CONSTRAINED ATOMS AS GROUPS
>
>
> 5. REFERENCE ATOM COORDINATES
>
>
> Mask :77-520, 521; matches 7002 atoms
>
> --------------------------------------------------------------------------------
> 3. ATOMIC COORDINATES AND VELOCITIES
> --------------------------------------------------------------------------------
>
>
>
> begin time read from input coords = 0.000 ps
>
> Number of triangulated 3-point waters found: 24238
> | Atom division among processors:
> | 0 22475 44951 67427 89900
>
> Sum of charges from parm topology file = -0.00074086
> Forcing neutrality...
> | Running AMBER/MPI version on 4 nodes
>
>
> --------------------------------------------------------------------------------
> 4. RESULTS
> --------------------------------------------------------------------------------
>
> ---------------------------------------------------
> APPROXIMATING switch and d/dx switch using CUBIC SPLINE INTERPOLATION
> using 5000.0 points per unit in tabled values
> TESTING RELATIVE ERROR over r ranging from 0.0 to cutoff
> | CHECK switch(x): max rel err = 0.2738E-14 at 2.422500
> | CHECK d/dx switch(x): max rel err = 0.8314E-11 at 2.736960
> ---------------------------------------------------
> | Local SIZE OF NONBOND LIST = 6454031
> | TOTAL SIZE OF NONBOND LIST = 25821257
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 1 2.8852E+18 1.4602E+18 3.9863E+20 O 17223
>
> BOND = 149932.7945 ANGLE = 46211.6820 DIHED = 9161.8416
> VDWAALS = ************* EEL = -246440.0364 HBOND = 0.0000
> 1-4 VDW = 29126.1421 1-4 EEL = 23905.1720 RESTRAINT = 0.0000
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 2 8.2905E+17 3.7817E+17 1.0324E+20 O 17223
>
> BOND = 149932.8374 ANGLE = 46211.6729 DIHED = 9161.8416
> VDWAALS = ************* EEL = -246403.3807 HBOND = 0.0000
> 1-4 VDW = 29126.1396 1-4 EEL = 23905.1707 RESTRAINT = 0.0000
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 3 2.1637E+17 8.8247E+16 2.4091E+19 O 17223
>
> BOND = 149932.9075 ANGLE = 46211.6637 DIHED = 9161.8416
> VDWAALS = ************* EEL = -246367.9129 HBOND = 0.0000
> 1-4 VDW = 29126.1367 1-4 EEL = 23905.1692 RESTRAINT = 0.0000
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 4 5.1492E+16 1.8633E+16 5.0866E+18 O 17223
>
> BOND = 149933.0186 ANGLE = 46211.6549 DIHED = 9161.8417
> VDWAALS = ************* EEL = -246334.1279 HBOND = 0.0000
> 1-4 VDW = 29126.1334 1-4 EEL = 23905.1675 RESTRAINT = 0.0000
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 5 1.1230E+16 3.5795E+15 9.7717E+17 O 17223
>
> BOND = 149933.1907 ANGLE = 46211.6479 DIHED = 9161.8419
> VDWAALS = ************* EEL = -246302.4146 HBOND = 0.0000
> 1-4 VDW = 29126.1294 1-4 EEL = 23905.1654 RESTRAINT = 0.0000
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 6 2.2583E+15 6.2972E+14 1.7191E+17 O 17223
>
> BOND = 149933.4532 ANGLE = 46211.6443 DIHED = 9161.8422
> VDWAALS = ************* EEL = -246273.0505 HBOND = 0.0000
> 1-4 VDW = 29126.1249 1-4 EEL = 23905.1629 RESTRAINT = 0.0000
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 7 4.2146E+14 1.0218E+14 2.7895E+16 O 17223
>
> BOND = 149933.8485 ANGLE = 46211.6471 DIHED = 9161.8428
> VDWAALS = ************* EEL = -246246.1983 HBOND = 0.0000
> 1-4 VDW = 29126.1198 1-4 EEL = 23905.1600 RESTRAINT = 0.0000
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 8 7.3514E+13 1.5407E+13 4.2061E+15 O 17223
>
> BOND = 149934.4390 ANGLE = 46211.6606 DIHED = 9161.8438
> VDWAALS = ************* EEL = -246221.9135 HBOND = 0.0000
> 1-4 VDW = 29126.1140 1-4 EEL = 23905.1567 RESTRAINT = 0.0000
> EAMBER = *************
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 9 1.2074E+13 2.1752E+12 5.9383E+14 O 17223
>
> BOND = 149935.3146 ANGLE = 46211.6910 DIHED = 9161.8452
> VDWAALS = ************* EEL = -246200.1611 HBOND = 0.0000
> 1-4 VDW = 29126.1076 1-4 EEL = 23905.1528 RESTRAINT = 0.0000
> EAMBER = *************
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 10 1.8860E+12 2.8967E+11 7.9077E+13 O 17223
>
> BOND = 149936.6064 ANGLE = 46211.7479 DIHED = 9161.8474
> VDWAALS = ************* EEL = -246180.8298 HBOND = 0.0000
> 1-4 VDW = 29126.1006 1-4 EEL = 23905.1483 RESTRAINT = 0.0000
> EAMBER = *************
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 11 2.8815E+11 3.6643E+10 1.0002E+13 O 17223
>
> BOND = 149938.5018 ANGLE = 46211.8441 DIHED = 9161.8506
> VDWAALS = ************* EEL = -246163.7303 HBOND = 0.0000
> 1-4 VDW = 29126.0930 1-4 EEL = 23905.1432 RESTRAINT = 0.0000
> EAMBER = *************
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 12 2.8815E+11 3.6643E+10 1.0002E+13 O 17223
>
> BOND = 149938.5018 ANGLE = 46211.8441 DIHED = 9161.8506
> VDWAALS = ************* EEL = -246163.7303 HBOND = 0.0000
> 1-4 VDW = 29126.0930 1-4 EEL = 23905.1432 RESTRAINT = 0.0000
> EAMBER = *************
>
>
>
> ____________________________________________________________________________________
> Never miss a thing. Make Yahoo your home page.
> http://www.yahoo.com/r/hs
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
--
Dr. Adrian E. Roitberg
Associate Professor
Quantum Theory Project and Department of Chemistry
University of Florida PHONE 352 392-6972
P.O. Box 118435 FAX 352 392-8722
Gainesville, FL 32611-8435 Email adrian_at_qtp.ufl.edu
============================================================================
To announce that there must be no criticism of the president,
or that we are to stand by the president right or wrong,
is not only unpatriotic and servile, but is morally treasonable
to the American public."
-- Theodore Roosevelt
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|