Log file opened on Fri Dec 19 12:01:52 2014 Host: etna.ii.uib.no pid: 19140 nodeid: 0 nnodes: 1 Gromacs version: VERSION 4.6.6 Precision: single Memory model: 64 bit MPI library: thread_mpi OpenMP support: enabled GPU support: disabled invsqrt routine: gmx_software_invsqrt(x) CPU acceleration: AVX_256 FFT library: fftpack (built-in) Large file support: enabled RDTSCP usage: enabled Built on: Wed Aug 13 11:30:34 CEST 2014 Built by: venken@etna.ii.uib.no [CMAKE] Build OS/arch: Linux 2.6.32-431.11.2.el6.x86_64 x86_64 Build CPU vendor: GenuineIntel Build CPU brand: Intel(R) Xeon(R) CPU E5-2609 v2 @ 2.50GHz Build CPU family: 6 Model: 62 Stepping: 4 Build CPU features: aes apic avx clfsh cmov cx8 cx16 f16c htt lahf_lm mmx msr nonstop_tsc pcid pclmuldq pdcm pdpe1gb popcnt pse rdrnd rdtscp sse2 sse3 sse4.1 sse4.2 ssse3 tdt x2apic C compiler: /usr/bin/cc GNU cc (GCC) 4.4.7 20120313 (Red Hat 4.4.7-4) C compiler flags: -mavx -Wextra -Wno-missing-field-initializers -Wno-sign-compare -Wall -Wno-unused -Wunused-value -fomit-frame-pointer -funroll-all-loops -O3 -DNDEBUG :-) G R O M A C S (-: GROtesk MACabre and Sinister :-) VERSION 4.6.6 (-: Contributions from Mark Abraham, Emile Apol, Rossen Apostolov, Herman J.C. Berendsen, Aldert van Buuren, Pär Bjelkmar, Rudi van Drunen, Anton Feenstra, Gerrit Groenhof, Christoph Junghans, Peter Kasson, Carsten Kutzner, Per Larsson, Pieter Meulenhoff, Teemu Murtola, Szilard Pall, Sander Pronk, Roland Schulz, Michael Shirts, Alfons Sijbers, Peter Tieleman, Berk Hess, David van der Spoel, and Erik Lindahl. Copyright (c) 1991-2000, University of Groningen, The Netherlands. Copyright (c) 2001-2012,2013, The GROMACS development team at Uppsala University & The Royal Institute of Technology, Sweden. check out http://www.gromacs.org for more information. This program is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation; either version 2.1 of the License, or (at your option) any later version. :-) mdrun (-: ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++ B. Hess and C. Kutzner and D. van der Spoel and E. Lindahl GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation J. Chem. Theory Comput. 4 (2008) pp. 435-447 -------- -------- --- Thank You --- -------- -------- ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++ D. van der Spoel, E. Lindahl, B. Hess, G. Groenhof, A. E. Mark and H. J. C. Berendsen GROMACS: Fast, Flexible and Free J. Comp. Chem. 26 (2005) pp. 1701-1719 -------- -------- --- Thank You --- -------- -------- ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++ E. Lindahl and B. Hess and D. van der Spoel GROMACS 3.0: A package for molecular simulation and trajectory analysis J. Mol. Mod. 7 (2001) pp. 306-317 -------- -------- --- Thank You --- -------- -------- ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++ H. J. C. Berendsen, D. van der Spoel and R. van Drunen GROMACS: A message-passing parallel molecular dynamics implementation Comp. Phys. Comm. 91 (1995) pp. 43-56 -------- -------- --- Thank You --- -------- -------- Input Parameters: integrator = steep nsteps = 10 init-step = 0 cutoff-scheme = Group ns_type = Grid nstlist = 10 ndelta = 2 nstcomm = 100 comm-mode = Linear nstlog = 100 nstxout = 0 nstvout = 0 nstfout = 0 nstcalcenergy = 100 nstenergy = 1000 nstxtcout = 100 init-t = 0 delta-t = 0.02 xtcprec = 10 fourierspacing = 0.12 nkx = 0 nky = 0 nkz = 0 pme-order = 4 ewald-rtol = 1e-05 ewald-geometry = 0 epsilon-surface = 0 optimize-fft = FALSE ePBC = xyz bPeriodicMols = FALSE bContinuation = FALSE bShakeSOR = FALSE etc = No bPrintNHChains = FALSE nsttcouple = -1 epc = No epctype = Isotropic nstpcouple = -1 tau-p = 1 ref-p (3x3): ref-p[ 0]={ 0.00000e+00, 0.00000e+00, 0.00000e+00} ref-p[ 1]={ 0.00000e+00, 0.00000e+00, 0.00000e+00} ref-p[ 2]={ 0.00000e+00, 0.00000e+00, 0.00000e+00} compress (3x3): compress[ 0]={ 0.00000e+00, 0.00000e+00, 0.00000e+00} compress[ 1]={ 0.00000e+00, 0.00000e+00, 0.00000e+00} compress[ 2]={ 0.00000e+00, 0.00000e+00, 0.00000e+00} refcoord-scaling = No posres-com (3): posres-com[0]= 0.00000e+00 posres-com[1]= 0.00000e+00 posres-com[2]= 0.00000e+00 posres-comB (3): posres-comB[0]= 0.00000e+00 posres-comB[1]= 0.00000e+00 posres-comB[2]= 0.00000e+00 verlet-buffer-drift = 0.005 rlist = 1.4 rlistlong = 1.4 nstcalclr = 0 rtpi = 0.05 coulombtype = Shift coulomb-modifier = None rcoulomb-switch = 0 rcoulomb = 1.2 vdwtype = Shift vdw-modifier = None rvdw-switch = 0.9 rvdw = 1.2 epsilon-r = 15 epsilon-rf = inf tabext = 1 implicit-solvent = No gb-algorithm = Still gb-epsilon-solvent = 80 nstgbradii = 1 rgbradii = 1 gb-saltconc = 0 gb-obc-alpha = 1 gb-obc-beta = 0.8 gb-obc-gamma = 4.85 gb-dielectric-offset = 0.009 sa-algorithm = Ace-approximation sa-surface-tension = 2.05016 DispCorr = No bSimTemp = FALSE free-energy = no nwall = 0 wall-type = 9-3 wall-atomtype[0] = -1 wall-atomtype[1] = -1 wall-density[0] = 0 wall-density[1] = 0 wall-ewald-zfac = 3 pull = no rotation = FALSE disre = No disre-weighting = Conservative disre-mixed = FALSE dr-fc = 1000 dr-tau = 0 nstdisreout = 100 orires-fc = 0 orires-tau = 0 nstorireout = 100 dihre-fc = 0 em-stepsize = 0.01 em-tol = 10 niter = 20 fc-stepsize = 0 nstcgsteep = 1000 nbfgscorr = 10 ConstAlg = Lincs shake-tol = 0.0001 lincs-order = 4 lincs-warnangle = 30 lincs-iter = 1 bd-fric = 0 ld-seed = 1993 cos-accel = 0 deform (3x3): deform[ 0]={ 0.00000e+00, 0.00000e+00, 0.00000e+00} deform[ 1]={ 0.00000e+00, 0.00000e+00, 0.00000e+00} deform[ 2]={ 0.00000e+00, 0.00000e+00, 0.00000e+00} adress = FALSE userint1 = 0 userint2 = 0 userint3 = 0 userint4 = 0 userreal1 = 0 userreal2 = 0 userreal3 = 0 userreal4 = 0 grpopts: nrdf: 105 ref-t: 0 tau-t: 0 anneal: No ann-npoints: 0 acc: 0 0 0 nfreeze: N N N energygrp-flags[ 0]: 0 efield-x: n = 0 efield-xt: n = 0 efield-y: n = 0 efield-yt: n = 0 efield-z: n = 0 efield-zt: n = 0 bQMMM = FALSE QMconstraints = 0 QMMMscheme = 0 scalefactor = 1 qm-opts: ngQM = 0 Using 1 MPI thread Detecting CPU-specific acceleration. Present hardware specification: Vendor: GenuineIntel Brand: Intel(R) Xeon(R) CPU E5-2609 v2 @ 2.50GHz Family: 6 Model: 62 Stepping: 4 Features: aes apic avx clfsh cmov cx8 cx16 f16c htt lahf_lm mmx msr nonstop_tsc pcid pclmuldq pdcm pdpe1gb popcnt pse rdrnd rdtscp sse2 sse3 sse4.1 sse4.2 ssse3 tdt x2apic Acceleration most likely to fit this hardware: AVX_256 Acceleration selected at GROMACS compile time: AVX_256 Table routines are used for coulomb: TRUE Table routines are used for vdw: TRUE Using shifted Lennard-Jones, switch between 0.9 and 1.2 nm Cut-off's: NS: 1.4 Coulomb: 1.2 LJ: 1.2 System total charge: 1.000 Generated table with 1200 data points for Shift. Tabscale = 500 points/nm Generated table with 1200 data points for LJ6Shift. Tabscale = 500 points/nm Generated table with 1200 data points for LJ12Shift. Tabscale = 500 points/nm Potential shift: LJ r^-12: 0.000 r^-6 0.000 Removing pbc first time Initializing LINear Constraint Solver ++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++ B. Hess P-LINCS: A Parallel Linear Constraint Solver for molecular simulation J. Chem. Theory Comput. 4 (2008) pp. 116-122 -------- -------- --- Thank You --- -------- -------- The number of constraints is 36 12 constraints are involved in constraint triangles, will apply an additional matrix expansion of order 4 for couplings between constraints inside triangles Initiating Steepest Descents Max number of connections per atom is 19 Total number of connections is 382 Max number of graph edges per atom is 3 Total number of graph edges is 102 Rel. Constraint Deviation: RMS MAX between atoms Before LINCS 0.270931 0.900270 3 4 After LINCS 0.000303 0.000727 44 45 Started Steepest Descents on node 0 Fri Dec 19 12:01:52 2014 Steepest Descents: Tolerance (Fmax) = 1.00000e+01 Number of steps = 10 Grid: 3 x 3 x 3 cells Step Time Lambda 0 0.00000 0.00000 Energies (kJ/mol) Bond G96Angle Proper Dih. Improper Dih. LJ (SR) 5.41972e+01 1.12410e+02 2.00628e+01 1.53785e+00 -1.56968e+02 Coulomb (SR) Potential Pressure (bar) Constr. rmsd 1.23559e-01 3.13632e+01 -2.92177e+01 3.03048e-04 Step Time Lambda 1 1.00000 0.00000 Energies (kJ/mol) Bond G96Angle Proper Dih. Improper Dih. LJ (SR) 5.06771e+01 1.01206e+02 1.39458e+01 1.56662e+00 -2.34590e+02 Coulomb (SR) Potential Pressure (bar) Constr. rmsd 1.22060e-01 -6.70727e+01 -4.18976e+01 1.31100e-06 Step Time Lambda 2 2.00000 0.00000 Energies (kJ/mol) Bond G96Angle Proper Dih. Improper Dih. LJ (SR) 4.50046e+01 8.75763e+01 8.00557e+00 1.75323e+00 -3.02143e+02 Coulomb (SR) Potential Pressure (bar) Constr. rmsd 1.19195e-01 -1.59685e+02 -5.24144e+01 2.77924e-06 Step Time Lambda 3 3.00000 0.00000 Energies (kJ/mol) Bond G96Angle Proper Dih. Improper Dih. LJ (SR) 3.58355e+01 7.51218e+01 4.77488e+00 2.12411e+00 -3.56110e+02 Coulomb (SR) Potential Pressure (bar) Constr. rmsd 1.13520e-01 -2.38140e+02 -6.26012e+01 5.85488e-06 Step Time Lambda 4 4.00000 0.00000 Energies (kJ/mol) Bond G96Angle Proper Dih. Improper Dih. LJ (SR) 2.13765e+01 6.08987e+01 2.82120e+00 2.48456e+00 -3.96307e+02 Coulomb (SR) Potential Pressure (bar) Constr. rmsd 1.02088e-01 -3.08623e+02 -6.98308e+01 1.30111e-05 Step Time Lambda 5 5.00000 0.00000 Energies (kJ/mol) Bond G96Angle Proper Dih. Improper Dih. LJ (SR) 9.66966e+00 4.80460e+01 1.02360e+01 2.87272e+00 -4.20156e+02 Coulomb (SR) Potential Pressure (bar) Constr. rmsd 8.74061e-02 -3.49244e+02 -6.74774e+01 2.08672e-05 Step Time Lambda 6 6.00000 0.00000 Step Time Lambda 7 7.00000 0.00000 Energies (kJ/mol) Bond G96Angle Proper Dih. Improper Dih. LJ (SR) 6.11962e+00 4.52793e+01 3.35540e+00 2.67977e+00 -4.25966e+02 Coulomb (SR) Potential Pressure (bar) Constr. rmsd 8.58974e-02 -3.68446e+02 -6.71071e+01 2.67643e-05 Step Time Lambda 8 8.00000 0.00000 Step Time Lambda 9 9.00000 0.00000 Energies (kJ/mol) Bond G96Angle Proper Dih. Improper Dih. LJ (SR) 6.61139e+00 4.13129e+01 3.87155e+00 2.80572e+00 -4.32188e+02 Coulomb (SR) Potential Pressure (bar) Constr. rmsd 8.38512e-02 -3.77503e+02 -6.73752e+01 5.77355e-06 Step Time Lambda 10 10.00000 0.00000 Energies (kJ/mol) Bond G96Angle Proper Dih. Improper Dih. LJ (SR) 4.32672e+00 4.20438e+01 4.14474e+00 2.71118e+00 -4.34519e+02 Coulomb (SR) Potential Pressure (bar) Constr. rmsd 8.09832e-02 -3.81211e+02 -6.49873e+01 8.50705e-06 Energy minimization reached the maximum number of steps before the forces reached the requested precision Fmax < 10. Steepest Descents did not converge to Fmax < 10 in 11 steps. Potential Energy = -3.8121140e+02 Maximum force = 4.7244736e+02 on atom 20 Norm of force = 1.9088266e+02 M E G A - F L O P S A C C O U N T I N G NB=Group-cutoff nonbonded kernels NxN=N-by-N cluster Verlet kernels RF=Reaction-Field VdW=Van der Waals QSTab=quadratic-spline table W3=SPC/TIP3p W4=TIP4p (single or pairs) V&F=Potential and force V=Potential only F=Force only Computing: M-Number M-Flops % Flops ----------------------------------------------------------------------------- NB VdW [V&F] 0.497390 0.497 49.2 NB Elec. [V&F] 0.004048 0.004 0.4 NS-Pairs 0.011854 0.249 24.6 Reset In Box 0.000528 0.002 0.2 Shift-X 0.000528 0.003 0.3 CG-CoM 0.000528 0.002 0.2 Bonds 0.000165 0.010 1.0 Angles 0.000473 0.079 7.9 Propers 0.000154 0.035 3.5 Impropers 0.000044 0.009 0.9 Virial 0.001023 0.018 1.8 Lincs 0.000792 0.048 4.7 Lincs-Mat 0.011352 0.045 4.5 Constraint-Vir 0.000396 0.010 0.9 ----------------------------------------------------------------------------- Total 1.011 100.0 ----------------------------------------------------------------------------- R E A L C Y C L E A N D T I M E A C C O U N T I N G Computing: Nodes Th. Count Wall t (s) G-Cycles % ----------------------------------------------------------------------------- Neighbor search 1 1 11 0.000 0.002 2.8 Force 1 1 11 0.000 0.001 1.5 Constraints 1 1 21 0.000 0.000 0.7 Rest 1 0.012 0.060 95.1 ----------------------------------------------------------------------------- Total 1 0.013 0.063 100.0 ----------------------------------------------------------------------------- Core t (s) Wall t (s) (%) Time: 0.000 0.013 0.0 (steps/hour) Performance: 3109800.4 Finished mdrun on node 0 Fri Dec 19 12:01:52 2014