 MCCCS Towhee - Version 7.0.1   (August 18 2001)    
 Copyright (C) <2011> Marcus G. Martin
 Code Contributors for 2011: Fred James; Marcus G. Martin; Suleiman Oloriegbe;
 Copyright (C) <2010> Marcus G. Martin
 Code Contributors for 2010: Rene Haber; Andrej Lajovic; Marcus G. Martin; Loukas Peristeras;
 Copyright (C) <2009> Marcus G. Martin
 Code Contributors for 2009: Rene Haber; Marcus G. Martin; Loukas Peristeras;
 Copyright (C) <2008> Marcus G. Martin
 Code Contributors for 2008: Rene Haber; Iyad A. Hijazi; Marcus G. Martin; Loukas Peristeras; Craig Tenney;
 Copyright (C) <2007> Marcus G. Martin
 Code Contributors for 2007: Bernhard Eckl; Frank Heilmann; Jianhui Li; Marcus G. Martin; Craig M. Tenney; Ozgur Yazaydin
 Copyright (C) <2006> Marcus G. Martin
 Code Contributors for 2006: Alan A. Chen; Arben Jusufi; Christian D. Lorenz; Marcus G. Martin; Jeffrey J. Potoff; Vincent K. Shen; Matthew A. Wyczalkowski; Ozgur Yazaydin
 Copyright (C) <2005> Marcus G. Martin
 Code Contributors for 2005: Marcus G. Martin; Alan A. Chen; Christian D. Lorenz; Matthew A. Wyczalkowski
 Copyright (C) <2004> Marcus G. Martin
 Code Contributors for 2004: Marcus G. Martin; Christian D. Lorenz; Peter A. Schultz; Aidan P. Thompson; Jeffrey J. Potoff; Steve R. Lustig; Paul S. Crozier; Nicholas du Preez
 Copyright (C) <2003> Marcus G. Martin
 Code Contributors for 2003: Marcus G. Martin; Christian D. Lorenz; David B. Ritch; Michael L. Greenfield; Yogesh Srivastava; Evangelos A Coutsias
 Copyright (C) <2002> Marcus G. Martin
 Code Contributors for 2002: Marcus G. Martin; Christian D. Lorenz; Aidan P. Thompson; Allen G. Sault
 Copyright (C) <2000-2001> Marcus G. Martin
 Copyright (C) <1999> J. Ilja Siepmann and Marcus G. Martin
 Code Contributors for 1999: J. Ilja Siepmann; Marcus G. Martin; Bin Chen; Collin D. Wick; John Stubbs

 Direct comments about this code to Marcus G. Martin (marcus_martin@users.sourceforge.net)
 See the MCCCS Towhee web site for more information and the users manual.
 http://towhee.sourceforge.net

 We would like to thank the following agencies for providing funding used to develop MCCCS Towhee
 -Department of Energy Industrial Technologies Program (MGM)
 -National Science Foundation (JIS)
 -Department of Energy Computational Science Graduate Fellowship (MGM,CDW)

 This program is free software; you can redistribute it and/or modify
 it under the terms of the GNU General Public License as published by
 the Free Software Foundation; either version 2 of the License, or
 (at your option) any later version.
 This program is distributed in the hope that it will be useful,
 but WITHOUT ANY WARRANTY; without even the implied warranty of
 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 GNU General Public License for more details.

 You should have received a copy of the GNU General Public License
 along with this program; if not, write to the Free Software
 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

 Reading from towhee_input file: towhee_input                                                                                        
 in directory: current directory
 inputformat: Towhee         
 ensemble: nvt                           
 temperature:  400.000    
 nmolty:            1
 nmolectyp:          160
 numboxes:            2
 stepstyle: cycles  
 nstep:            1
 printfreq:         1000
 blocksize:         2000
 moviefreq:       100000
 backupfreq:         1000
 runoutput: full                
Full output of updates and block averages
 pdb_output_freq:        10000
 pressurefreq:           20
 trmaxdispfreq:         1000
 volmaxdispfreq:         1000
 potentialstyle: internal                      
 ffnumber:            1
 ff_filename:
/towheebase/ForceFields/towhee_ff_SMMKnaip                                                          
 classical_potential: Lennard-Jones                 
 READCLASSICAL: pot_num:     1 potential name: Lennard-Jones                 
 classical_mixrule: Lorentz-Berthelot             
 lshift:  F
 ltailc:  T
 rmin:  2.00000    
 rcut:  14.0000    
 rcutin:  5.00000    
 electrostatic_form: none                                              
 Setting up force field parameters from files
 opening forcefield file:    1
 Lorentz-Berthelot Mixing rules
 Arithmetic mean of sigma terms
 Geometric mean of epsilon term
default max_bond_length:  3.00
 No solvation model used
 linit:  T
 initboxtype: dimensions          
initstyle Box:    1
 full cbmc           
initstyle Box:    2
 full cbmc           
Box:  1 initlattice: simple cubic        
Box:  2 initlattice: simple cubic        
Box:  1 initmol:       120
Box:  2 initmol:        40
Box:  1 inix,iniy,iniz:         5          5          5
Box:  2 inix,iniy,iniz:         4          3          4
Box idim hmatrix:    1 1   31.00000    0.00000    0.00000
Box idim hmatrix:    1 2    0.00000   31.00000    0.00000
Box idim hmatrix:    1 3    0.00000    0.00000   31.00000
Box idim hmatrix:    2 1   90.00000    0.00000    0.00000
Box idim hmatrix:    2 2    0.00000   90.00000    0.00000
Box idim hmatrix:    2 3    0.00000    0.00000   90.00000
itest:  1 pairbox:  1   2
 pmvol: 0.200000E-02
 pmvlpr:  1.00000    
 rmvol: 0.100000    
 tavol: 0.500000    
 pmcell:  0.00000    
 pmcellpr:  1.00000    
pairbox: 1  2 pmcellpt: 0.500000
 rmcell:  1.00000    
 tacell: 0.500000    
 pm2boxrbswap:  0.00000    
 pm2rbswmt:  1.00000    
 pm2rbswpr:  1.00000    
 pm2boxcbswap: 0.100000E-01
 pm2cbswmt:  1.00000    
 pm2cbswpr:  1.00000    
 pm1boxcbswap:  0.00000    
 pm1cbswmt:  1.00000    
 pmavb1: 0.200000E-01
 pmavb1in: 0.500000    
 pmavb1mt:  1.00000    
moltyp:  1 pmavb1ct: 1.000000
 avb1rad:  8.00000    
 pmavb2:  0.00000    
 pmavb2in: 0.500000    
 pmavb2mt:  1.00000    
moltyp:  1 pmavb2ct: 1.000000
 avb2rad:  8.00000    
 pmavb3:  0.00000    
 pmavb3mt:  1.00000    
moltyp:  1 pmavb3ct: 1.000000
 avb3rad:  8.00000    
 pmcb: 0.330000    
 pmcbmt:  1.00000    
 pmall:  0.00000    
 pmback:  0.00000    
 pmbkmt:  1.00000    
 pmpivot:  0.00000    
 pmpivmt:  1.00000    
 pmconrot:  0.00000    
 pmcrmt:  1.00000    
 pmcrback:  0.00000    
 pmcrbmt:  1.00000    
 pmplane:  0.00000    
 pmplanebox: 0.500000     1.00000    
 planewidth:  3.00000    
 pmrow:  0.00000    
 pmrowbox: 0.500000     1.00000    
 rowwidth:  3.00000    
 pmtraat:  0.00000    
 pmtamt:  1.00000    
 rmtraa: 0.500000    
 tatraa: 0.500000    
 pmtracm: 0.670000    
 pmtcmt:  1.00000    
 rmtrac: 0.500000    
 tatrac: 0.500000    
 pmrotate:  1.00000    
 pmromt:  1.00000    
 rmrot: 0.500000E-01
 tarot: 0.500000    
cbmc_formulation: Martin and Siepmann 1999 + Martin and Thompson 2004         
 Original Coupled-Decoupled formulation from
 M.G. Martin; J.I. Siepmann; J. Phys. Chem. B103 2977-2980 (1999)
 with a decoupled bond length selection from
 M.G. Martin; A.P. Thompson; Fluid Phase Equilib. 217 105-110 (2004)
cbmc_setting_style: default ideal                 
 input_style: basic connectivity map                            
 nunit:            7
 nmaxcbmc:            7
 lpdbnames:  F
   using the SMMKnaip   force field
 charge_assignment: none                          
   Building the input file for molecule type:     1
unit:    1 name:CH3       
unit:    2 name:CH2       
unit:    3 name:CH        
unit:    4 name:CH2       
unit:    5 name:CH3       
unit:    6 name:CH2       
unit:    7 name:CH3       
 Verifying input structures are consistent
Determining CBMC bond distributions
     for molecule type     1
 Determining cyclic subunits for molecule type          1
Default total charge on molecule   1 is   0.00000
Total charge in the simulation system:   0.00000
Bond Types
Type:      1 Style: Fixed  Length: 1.5400
Angle Types
Type:   1 Style: Standard Harmonic  Angle:    114.000 Constant:    31250.0
Type:   2 Style: Standard Harmonic  Angle:    112.000 Constant:    31250.0
Torsion Types
Type:   2 Style: Old UA OPLS Cosine Series
          k0:    -251.1 k1:     428.7 k2:    -111.8 k3:     441.3
Improper Torsion Types
   No Improper Types
 Canonical Gibbs ensemble
 3-dimensional periodic box
 Additional Center-of-Mass cutoff
 Dual Cutoff Configurational-bias Monte Carlo
 Coupled-decoupled Configurational-bias MC

Molecular mass for molecule type     1 is   100.1334 g/mol
 Calling initconf
Selected random number generator DX-1597-2-7 and single integer seed    1302002
Testing random number generator using the single integer seed
    0.812479  0.687663  0.080298  0.200888  0.516914
10 million RNG sum (approximately 5 million):   4999978.2264809748
Box idim hinverse:    1 1    0.03226    0.00000    0.00000
Box idim hinverse:    1 2    0.00000    0.03226    0.00000
Box idim hinverse:    1 3    0.00000    0.00000    0.03226
Box idim hinverse:    2 1    0.01111    0.00000    0.00000
Box idim hinverse:    2 2    0.00000    0.01111    0.00000
Box idim hinverse:    2 3    0.00000    0.00000    0.01111
INITCONF: building template for molecule type:    1
Initial Structure Growth Trial:    1
 Finished initconf
Energies exclusively from internal potentials
Nonbonded Force Field
Lennard-Jones 12-6 potential
    with tail corrections
u(r) = 4*epsilon[(sigma/r)^12 - (sigma/r)^6] - shift
Num. Atom(i)    Num. Atom(j)         sigma    epsilon      shift     1-4sig     1-4eps
   1 CH3           1 CH3            3.7700    98.1000     0.0000     0.0000     0.0000
   1 CH3           2 CH2            3.8500    67.9021     0.0000     0.0000     0.0000
   1 CH3           3 CH             3.9350    34.3103     0.0000     0.0000     0.0000
   2 CH2           2 CH2            3.9300    47.0000     0.0000     0.0000     0.0000
   2 CH2           3 CH             4.0150    23.7487     0.0000     0.0000     0.0000
   3 CH            3 CH             4.1000    12.0000     0.0000     0.0000     0.0000

Number of MC cycles:                    1
Number of molecules:                  160
Temperature [K]:       400.00000

Initial Energies for Box     1
Total molecules in this box        120
Molecules of type   1 :        120
 total vibration            0.000 [K]         0.00000 [kcal/mol]
    regular                 0.000 [K]         0.00000 [kcal/mol]
    bond-bond(1-2)          0.000 [K]         0.00000 [kcal/mol]
 total angle           305435.783 [K]       606.96683 [kcal/mol]
    regular            305435.783 [K]       606.96683 [kcal/mol]
    angle-angle             0.000 [K]         0.00000 [kcal/mol]
 total torsion         112318.185 [K]       223.20048 [kcal/mol]
    regular            112318.185 [K]       223.20048 [kcal/mol]
    improper                0.000 [K]         0.00000 [kcal/mol]
 total nonbond        1166616.252 [K]      2318.31830 [kcal/mol]
    intramolecular     -31361.609 [K]       -62.32229 [kcal/mol]
    2-body nonbond    1211983.389 [K]      2408.47259 [kcal/mol]
    3-body nonbond          0.000 [K]         0.00000 [kcal/mol]
    tail correct.      -14005.528 [K]       -27.83201 [kcal/mol]
 external field             0.000 [K]         0.00000 [kcal/mol]
 solvation                  0.000 [K]         0.00000 [kcal/mol]
 total classical           1584370.2194 [K]          3148.485601 [kcal/mol]
Initial Energies for Box     2
Total molecules in this box         40
Molecules of type   1 :         40
 total vibration            0.000 [K]         0.00000 [kcal/mol]
    regular                 0.000 [K]         0.00000 [kcal/mol]
    bond-bond(1-2)          0.000 [K]         0.00000 [kcal/mol]
 total angle           101811.928 [K]       202.32228 [kcal/mol]
    regular            101811.928 [K]       202.32228 [kcal/mol]
    angle-angle             0.000 [K]         0.00000 [kcal/mol]
 total torsion          37439.395 [K]        74.40016 [kcal/mol]
    regular             37439.395 [K]        74.40016 [kcal/mol]
    improper                0.000 [K]         0.00000 [kcal/mol]
 total nonbond         -10517.464 [K]       -20.90047 [kcal/mol]
    intramolecular     -10453.870 [K]       -20.77410 [kcal/mol]
    2-body nonbond          0.000 [K]         0.00000 [kcal/mol]
    3-body nonbond          0.000 [K]         0.00000 [kcal/mol]
    tail correct.         -63.594 [K]        -0.12637 [kcal/mol]
 external field             0.000 [K]         0.00000 [kcal/mol]
 solvation                  0.000 [K]         0.00000 [kcal/mol]
 total classical            128733.8591 [K]           255.821964 [kcal/mol]
 initial virial pressure in box  1 =     4908938.04
 initial virial pressure in box  2 =         300.62

 +++++ start of markov chain +++++

Cycle      Box   Energy [K]  Volume [A^3] Press. [kPa] Molecules

 +++++ end of markov chain +++++

 Final hmatrix (general box dimensions) 
Box:     1
  hmatrix(1,x)      31.00000       0.00000       0.00000
  hmatrix(2,x)       0.00000      31.00000       0.00000
  hmatrix(3,x)       0.00000       0.00000      31.00000
Box:     2
  hmatrix(1,x)      90.00000       0.00000       0.00000
  hmatrix(2,x)       0.00000      90.00000       0.00000
  hmatrix(3,x)       0.00000       0.00000      90.00000

* Aggregation-Volume-Bias-1 Moves *
 Box Moved Around From To  Attempted  Accepted  Percent
   1     1      1  out out         1         0     0.00

* Configurational-Bias REGROWTH Moves *
 Molecule type:    1 Box:    1
     Length  Attempts  Regrown  Accepted  %Regrown   %Accep.
        1         9         9         4    100.00     44.44
        4        10        10         5    100.00     50.00
        5         4         4         1    100.00     25.00
        6        14        14         0    100.00      0.00
 Molecule type:    1 Box:    2
     Length  Attempts  Regrown  Accepted  %Regrown   %Accep.
        1         2         2         1    100.00     50.00
        4         3         3         2    100.00     66.67
        6         2         2         1    100.00     50.00

* COM Translation Moves *
Molecule:  1 Box: 1 Attempts:             49. Accepted:             22. Accepted:  44.898 %
Molecule:  1 Box: 2 Attempts:             16. Accepted:             16. Accepted: 100.000 %

* Rotation Moves *
Molecule:  1 Box: 1 Attempts:             36. Accepted:             18. Accepted:  50.000 %
Molecule:  1 Box: 2 Attempts:             14. Accepted:             14. Accepted: 100.000 %

Final Energies for Box     1
Total molecules in this box        120
Molecules of type   1 :        120
 total vibration            0.000 [K]         0.00000 [kcal/mol]
    regular                 0.000 [K]         0.00000 [kcal/mol]
    bond-bond(1-2)          0.000 [K]         0.00000 [kcal/mol]
 total angle           293148.516 [K]       582.54938 [kcal/mol]
    regular            293148.516 [K]       582.54938 [kcal/mol]
    angle-angle             0.000 [K]         0.00000 [kcal/mol]
 total torsion         115631.631 [K]       229.78501 [kcal/mol]
    regular            115631.631 [K]       229.78501 [kcal/mol]
    improper                0.000 [K]         0.00000 [kcal/mol]
 total nonbond         920687.072 [K]      1829.60394 [kcal/mol]
    intramolecular     -27671.071 [K]       -54.98839 [kcal/mol]
    2-body nonbond     962363.671 [K]      1912.42433 [kcal/mol]
    3-body nonbond          0.000 [K]         0.00000 [kcal/mol]
    tail correct.      -14005.528 [K]       -27.83201 [kcal/mol]
 external field             0.000 [K]         0.00000 [kcal/mol]
 solvation                  0.000 [K]         0.00000 [kcal/mol]
 total classical           1329467.2200 [K]          2641.938322 [kcal/mol]
Final Energies for Box     2
Total molecules in this box         40
Molecules of type   1 :         40
 total vibration            0.000 [K]         0.00000 [kcal/mol]
    regular                 0.000 [K]         0.00000 [kcal/mol]
    bond-bond(1-2)          0.000 [K]         0.00000 [kcal/mol]
 total angle            98850.297 [K]       196.43688 [kcal/mol]
    regular             98850.297 [K]       196.43688 [kcal/mol]
    angle-angle             0.000 [K]         0.00000 [kcal/mol]
 total torsion          37578.764 [K]        74.67712 [kcal/mol]
    regular             37578.764 [K]        74.67712 [kcal/mol]
    improper                0.000 [K]         0.00000 [kcal/mol]
 total nonbond          -9926.948 [K]       -19.72699 [kcal/mol]
    intramolecular      -9863.354 [K]       -19.60061 [kcal/mol]
    2-body nonbond          0.000 [K]         0.00000 [kcal/mol]
    3-body nonbond          0.000 [K]         0.00000 [kcal/mol]
    tail correct.         -63.594 [K]        -0.12637 [kcal/mol]
 external field             0.000 [K]         0.00000 [kcal/mol]
 solvation                  0.000 [K]         0.00000 [kcal/mol]
 total classical            126502.1136 [K]           251.387004 [kcal/mol]
Averages               Units Type       Box  1      Box  2
 Volume                 nm^3       0.29791E+02 0.72900E+03
 Volume^2               nm^6       0.88750E+03 0.53144E+06
 Molecule Number                1      120.000      40.000
 Molar Volume         ml/mol       0.14950E+03 0.10975E+05
 Specific Density       g/ml        0.66977839  0.00912363
 Number Density         nm-3    1      4.02806     0.05487
 Mole Fraction                  1    1.0000000   1.0000000
 Radius of Gyration        A    1    1.9715061   1.9705753
 Ideal Pressure          kPa       0.22252E+05 0.30311E+03
 Ideal p_i <N/V>kT       kPa    1  0.22252E+05 0.30311E+03
 Total Classical           K        0.1477E+07  0.1271E+06
 Inter vdw                 K        0.1095E+07 -0.6359E+02
 Angle                     K        0.2984E+06  0.9992E+05
 Torsion                   K        0.1137E+06  0.3715E+05
 Intra vdw                 K       -0.3062E+05 -0.9889E+04
 External Field            K        0.0000E+00  0.0000E+00
 Vibration                 K        0.0000E+00  0.0000E+00
 Coulomb                   K        0.0000E+00  0.0000E+00
 Tail vdw                  K       -0.1401E+05 -0.6359E+02
 Solvation                 K        0.0000E+00  0.0000E+00
 G: Sum{<u_i><N_i>}   kJ/mol        0.0000E+00  0.0000E+00
 U                    kJ/mol        0.1228E+05  0.1057E+04
 H_vap (pV/n=RT)      kJ/mol       -0.7256E+02
Please see towhee_citations for a list of suggested citations for this simulation
