Wall clock time and date at job start Wed Apr 1 2020 03:36:36 AMSOL-version 7.1 by G. D. Hawkins, D. J. Giesen, G. C. Lynch, C. C. Chambers, I. Rossi, J. W. Storer, J. Li, J. D. Thompson, P. Winget, B. J. Lynch, D. Rinaldi, D. A. Liotard, C. J. Cramer, and D. G. Truhlar Copyright 2004 by Regents of the University of Minnesota. All rights reserved. Notice: recipients of this code are asked to comply with the user agreement in Section 1 of the documentation file. ******************************************************************************* * 1SCF - SCF CALCULATION WITHOUT GEOMETRY OPTIMIZATION * - USE EF ROUTINE FOR MINIMUM SEARCH (DEFAULT) * GEO-OK - OVERRIDE INTERATOMIC DISTANCE CHECK * TLIMIT= - A TIME OF 15. SECONDS REQUESTED * CHARGE - CHARGE ON SYSTEM= -1 * AM1 - THE AM1 HAMILTONIAN TO BE USED * DEV - DEVELOPER OPTIONS ARE ALLOWED * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS USER-SPECIFIED * DIELEC - THE SOLVENT DIELECTRIC CONSTANT IS 2.06 * IOFR - THE SOLVENT INDEX OF REFRACTION IS 1.4345 * ALPHA - THE SOLVENT ALPHA IS 0.00 * BETA - THE SOLVENT BETA IS 0.00 * GAMMA - THE SOLVENT MACROSCOPIC SURFACE TENSION IS * 38.93 CAL MOL^-1 ANGSTROM^-2 * FACARB - THE FRACTION OF AROMATIC CARBONS IS 0.00 * FEHALO - THE FRACTION OF ELECTRONEGATIVE HALOGENS IS 0.00 ******************************************************************************* Atom NO. Chemical Bond length Bond angle Dihedral angle number (k) symbol (angstroms) (degrees) (degrees) (I) NA:I NB:NA:I NC:NB:NA:I NA NB NC 1 1 C 2 2 C 1.52998 * 1 3 3 H 1.08999 * 109.47357 * 2 1 4 4 O 1.42904 * 109.47264 * 119.99865 * 2 1 3 5 5 C 1.53001 * 109.47290 * 239.99930 * 2 1 3 6 6 C 1.53003 * 109.47040 * 174.99882 * 5 2 1 7 7 C 1.53001 * 109.47040 * 180.02562 * 6 5 2 8 8 N 1.46897 * 109.47274 * 180.02562 * 7 6 5 9 9 C 1.46894 * 111.00430 * 180.02562 * 8 7 6 10 10 H 1.08998 * 109.46828 * 325.00585 * 9 8 7 11 11 C 1.53006 * 109.47224 * 84.99711 * 9 8 7 12 12 C 1.52997 * 109.47164 * 180.02562 * 11 9 8 13 13 C 1.53002 * 109.46978 * 299.99874 * 12 11 9 14 14 H 1.08992 * 109.47056 * 300.00064 * 13 12 11 15 15 N 1.46505 * 109.46958 * 180.02562 * 13 12 11 16 16 C 1.34777 * 119.99797 * 84.99979 * 15 13 12 17 17 O 1.21278 * 120.00140 * 0.02562 * 16 15 13 18 18 C 1.50697 * 119.99981 * 179.97438 * 16 15 13 19 19 S 1.80998 * 109.47393 * 180.02562 * 18 16 15 20 20 C 1.76191 * 100.00112 * 180.02562 * 19 18 16 21 21 H 1.08003 * 120.00195 * 0.02562 * 20 19 18 22 22 C 1.38805 * 120.00294 * 179.97438 * 20 19 18 23 23 N 1.31627 * 120.00522 * 180.02562 * 22 20 19 24 24 Si 1.86803 * 119.99596 * 0.02562 * 22 20 19 25 25 H 1.48503 * 109.46964 * 89.99853 * 24 22 20 26 26 H 1.48496 * 109.47163 * 209.99738 * 24 22 20 27 27 H 1.48506 * 109.47253 * 330.00125 * 24 22 20 28 28 C 1.52996 * 109.46891 * 60.00410 * 13 12 11 29 29 C 1.52997 * 109.47515 * 205.00020 * 9 8 7 30 30 H 1.09000 * 109.46879 * 179.97438 * 1 2 3 31 31 H 1.08998 * 109.47530 * 300.00011 * 1 2 3 32 32 H 1.09001 * 109.46948 * 60.00434 * 1 2 3 33 33 H 0.96706 * 113.99641 * 59.99885 * 4 2 1 34 34 H 1.08998 * 109.47159 * 55.00019 * 5 2 1 35 35 H 1.09001 * 109.47307 * 294.99994 * 5 2 1 36 36 H 1.08998 * 109.47033 * 60.00060 * 6 5 2 37 37 H 1.09001 * 109.47162 * 300.00201 * 6 5 2 38 38 H 1.09000 * 109.46678 * 59.99939 * 7 6 5 39 39 H 1.09001 * 109.47307 * 300.00112 * 7 6 5 40 40 H 1.00902 * 110.99811 * 303.95484 * 8 7 6 41 41 H 1.09001 * 109.47002 * 60.00120 * 11 9 8 42 42 H 1.09000 * 109.47050 * 300.00387 * 11 9 8 43 43 H 1.09002 * 109.47045 * 180.02562 * 12 11 9 44 44 H 1.09001 * 109.47375 * 59.99966 * 12 11 9 45 45 H 0.97001 * 119.99745 * 265.00017 * 15 13 12 46 46 H 1.08998 * 109.47252 * 300.00615 * 18 16 15 47 47 H 1.09005 * 109.46927 * 60.00159 * 18 16 15 48 48 H 0.97001 * 119.99838 * 180.02562 * 23 22 20 49 49 H 1.09004 * 109.47368 * 179.97438 * 28 13 12 50 50 H 1.09001 * 109.47372 * 59.99964 * 28 13 12 51 51 H 1.09002 * 109.47153 * 59.99762 * 29 9 8 52 52 H 1.09001 * 109.46832 * 299.99953 * 29 9 8 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 6 0.0000 0.0000 0.0000 2 6 1.5300 0.0000 0.0000 3 1 1.8934 1.0276 0.0000 4 8 2.0064 -0.6736 1.1668 5 6 2.0400 -0.7213 -1.2492 6 6 3.5649 -0.6132 -1.3136 7 6 4.0749 -1.3340 -2.5632 8 7 5.5389 -1.2308 -2.6248 9 6 6.0652 -1.9125 -3.8148 10 1 5.4606 -2.7947 -4.0252 11 6 6.0134 -0.9614 -5.0121 12 6 6.5610 -1.6716 -6.2517 13 6 8.0094 -2.0949 -5.9988 14 1 8.6142 -1.2128 -5.7886 15 7 8.5339 -2.7744 -7.1860 16 6 9.0569 -2.0548 -8.1985 17 8 9.0937 -0.8447 -8.1266 18 6 9.5970 -2.7538 -9.4194 19 16 10.2254 -1.5246 -10.5899 20 6 10.7775 -2.5871 -11.8824 21 1 10.6683 -3.6572 -11.7856 22 6 11.3519 -2.0489 -13.0257 23 7 11.7639 -2.8426 -13.9916 24 14 11.5402 -0.1980 -13.1936 25 1 12.8531 0.2184 -12.6384 26 1 11.4675 0.1789 -14.6281 27 1 10.4498 0.4802 -12.4476 28 6 8.0611 -3.0460 -4.8016 29 6 7.5136 -2.3354 -3.5616 30 1 -0.3633 -1.0277 -0.0005 31 1 -0.3634 0.5138 0.8899 32 1 -0.3633 0.5138 -0.8900 33 1 1.7210 -1.5953 1.2329 34 1 1.6054 -0.2618 -2.1370 35 1 1.7517 -1.7715 -1.2056 36 1 3.9995 -1.0727 -0.4259 37 1 3.8532 0.4370 -1.3573 38 1 3.6402 -0.8740 -3.4507 39 1 3.7865 -2.3842 -2.5200 40 1 5.8348 -0.2664 -2.6004 41 1 6.6180 -0.0791 -4.8017 42 1 4.9815 -0.6601 -5.1925 43 1 6.5245 -0.9939 -7.1046 44 1 5.9562 -2.5537 -6.4620 45 1 8.5048 -3.7423 -7.2433 46 1 10.4059 -3.4232 -9.1268 47 1 8.8005 -3.3303 -9.8900 48 1 12.1656 -2.4664 -14.7903 49 1 9.0929 -3.3479 -4.6216 50 1 7.4563 -3.9281 -5.0118 51 1 7.5505 -3.0130 -2.7086 52 1 8.1182 -1.4531 -3.3512 RHF calculation, no. of doubly occupied orbitals= 64 REFERENCES FOR PARAMETERS IN GAS-PHASE HAMILTONIAN: H: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985) C: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985) N: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985) O: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985) Si: (AM1): M.J.S.DEWAR, C. JIE, ORGANOMETALLICS, 6, 1486-1490 (1987). S: (AM1): M.J.S. DEWAR, Y.-C. YUAN, INORGANIC CHEM., 29, 589 (1990) REFERENCE FOR CHARGE MODEL 2: J. Li, J. Xing, C. J. Cramer, and D. G. Truhlar, J. Chem. Phys. 111 (1999) 885. REFERENCE FOR THE PARAMETERIZATION OF THE SM5.42R SOLVATION MODEL: J. Li, T. Zhu, G. D. Hawkins, P. Winget, D. A. Liotard, C. J. Cramer, and D. G. Truhlar, Theor. Chem. Acc. 103 (1999) 9-63 P. Winget, J. D. Thompson C. J. Cramer,and D. G. Truhlar, J. Phys. Chem. B. submitted. CHARGE=-1 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=GENORG IOFR=1.4345 ALPHA=0.00 BETA=0.00 GAMMA=38.93 & DIELEC=2.06 FACARB=0.00 FEHALO=0.00 DEV ZINC001495442992.mol2 52 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Wed Apr 1 2020 03:36:36 Heat of formation + Delta-G solvation = -127.322602 kcal Electronic energy + Delta-G solvation = -27005.105686 eV Core-core repulsion = 23128.919125 eV Total energy + Delta-G solvation = -3876.186562 eV No. of doubly occupied orbitals = 64 Molecular weight (most abundant/longest-lived isotopes) = 344.192 amu Computer time = 0.79 seconds Orbital eigenvalues (eV) -40.01985 -38.78964 -38.32793 -36.13117 -34.10678 -32.66041 -32.43999 -31.04447 -30.89963 -29.53988 -27.59051 -25.63276 -25.24534 -23.03173 -22.29979 -21.91526 -21.67537 -20.44003 -20.06000 -17.86231 -16.91389 -16.70383 -16.57489 -16.39991 -15.92058 -15.36825 -15.13222 -14.99196 -14.92739 -14.57417 -14.35487 -14.00402 -13.75003 -13.72387 -13.41660 -13.26179 -13.02192 -12.80084 -12.67683 -12.49339 -12.38374 -12.19010 -12.02458 -11.92297 -11.50281 -11.44590 -11.33899 -11.26342 -11.25114 -11.14357 -11.01271 -10.74367 -10.64893 -10.53127 -10.36067 -9.85782 -9.48711 -9.01726 -8.93827 -8.71567 -8.50225 -6.91978 -6.29327 -4.00939 2.76887 2.98266 3.29933 3.33717 3.50054 3.59592 3.70196 3.89907 4.22310 4.32815 4.38842 4.44532 4.48231 4.55816 4.66138 4.71554 4.78962 4.90698 4.92644 4.95427 4.98080 5.06119 5.09822 5.11543 5.20274 5.22317 5.23037 5.29292 5.33154 5.37880 5.47302 5.50668 5.53262 5.60630 5.73727 5.77387 5.79197 5.83106 5.84892 5.90159 5.92719 6.01632 6.11712 6.16091 6.60531 6.70295 7.01793 7.13941 7.64420 7.93751 8.41246 8.51354 9.16971 10.79937 Molecular weight = 344.19amu Principal moments of inertia in cm(-1) A = 0.029870 B = 0.001436 C = 0.001402 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 937.172883 B =19496.420736 C =19962.569851 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.187 4.187 2 C 0.101 3.899 3 H 0.097 0.903 4 O -0.563 6.563 5 C -0.151 4.151 6 C -0.115 4.115 7 C 0.058 3.942 8 N -0.686 5.686 9 C 0.087 3.913 10 H 0.076 0.924 11 C -0.144 4.144 12 C -0.126 4.126 13 C 0.147 3.853 14 H 0.087 0.913 15 N -0.734 5.734 16 C 0.531 3.469 17 O -0.520 6.520 18 C -0.125 4.125 19 S -0.142 6.142 20 C -0.522 4.522 21 H 0.088 0.912 22 C 0.096 3.904 23 N -0.881 5.881 24 Si 0.857 3.143 25 H -0.267 1.267 26 H -0.274 1.274 27 H -0.236 1.236 28 C -0.125 4.125 29 C -0.101 4.101 30 H 0.058 0.942 31 H 0.063 0.937 32 H 0.066 0.934 33 H 0.374 0.626 34 H 0.070 0.930 35 H 0.066 0.934 36 H 0.077 0.923 37 H 0.060 0.940 38 H 0.028 0.972 39 H 0.069 0.931 40 H 0.341 0.659 41 H 0.067 0.933 42 H 0.060 0.940 43 H 0.070 0.930 44 H 0.062 0.938 45 H 0.398 0.602 46 H 0.096 0.904 47 H 0.096 0.904 48 H 0.271 0.729 49 H 0.072 0.928 50 H 0.067 0.933 51 H 0.073 0.927 52 H 0.065 0.935 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -20.533 0.229 28.683 35.276 Note: The Mulliken population analysis charges presented below were not used in the solvation calculation but are provided for completeness. The chosen solvation model uses CM2 partial charges for calculating solvation energies. Net atomic charges, atomic populations, and dipole contributions using Mulliken population analysis Atom NO. Type Charge No. of electrons 1 C -0.244 4.244 2 C 0.042 3.958 3 H 0.115 0.885 4 O -0.368 6.368 5 C -0.189 4.189 6 C -0.152 4.152 7 C -0.071 4.071 8 N -0.315 5.315 9 C -0.023 4.023 10 H 0.095 0.905 11 C -0.184 4.184 12 C -0.165 4.165 13 C 0.042 3.958 14 H 0.105 0.895 15 N -0.387 5.387 16 C 0.317 3.683 17 O -0.394 6.394 18 C -0.275 4.275 19 S 0.084 5.916 20 C -0.669 4.669 21 H 0.106 0.894 22 C -0.112 4.112 23 N -0.517 5.517 24 Si 0.681 3.319 25 H -0.192 1.192 26 H -0.199 1.199 27 H -0.160 1.160 28 C -0.164 4.164 29 C -0.140 4.140 30 H 0.077 0.923 31 H 0.082 0.918 32 H 0.085 0.915 33 H 0.221 0.779 34 H 0.089 0.911 35 H 0.085 0.915 36 H 0.096 0.904 37 H 0.078 0.922 38 H 0.047 0.953 39 H 0.088 0.912 40 H 0.159 0.841 41 H 0.086 0.914 42 H 0.079 0.921 43 H 0.089 0.911 44 H 0.081 0.919 45 H 0.232 0.768 46 H 0.115 0.885 47 H 0.115 0.885 48 H 0.081 0.919 49 H 0.091 0.909 50 H 0.086 0.914 51 H 0.092 0.908 52 H 0.084 0.916 Dipole moment (debyes) X Y Z Total from point charges -20.010 2.274 29.151 35.431 hybrid contribution -0.473 -1.133 -1.502 1.941 sum -20.483 1.141 27.649 34.428 Atomic orbital electron populations 1.22519 0.96546 1.02562 1.02816 1.21929 0.92623 0.96803 0.84449 0.88486 1.86606 1.77174 1.29051 1.43980 1.22285 0.96584 1.01299 0.98722 1.21506 0.94362 1.00687 0.98685 1.22006 0.87758 1.00649 0.96720 1.60599 1.09730 1.20658 1.40560 1.21642 0.96378 0.96702 0.87588 0.90531 1.22212 1.01048 0.99299 0.95817 1.21860 0.95890 1.00716 0.98035 1.20918 0.94378 0.96281 0.84229 0.89493 1.46056 1.63781 1.10219 1.18633 1.20281 0.76684 0.87762 0.83564 1.90690 1.55250 1.13885 1.79617 1.23290 1.03933 1.03522 0.96796 1.84777 1.78906 1.11594 1.16303 1.22325 1.46176 0.95417 1.03029 0.89408 1.24707 0.93536 0.98651 0.94349 1.69938 1.39113 1.35394 1.07273 0.86738 0.81078 0.85411 0.78641 1.19193 1.19937 1.15955 1.21814 1.01987 0.98402 0.94198 1.21424 0.94731 1.00788 0.97032 0.92317 0.91806 0.91547 0.77923 0.91111 0.91531 0.90384 0.92158 0.95333 0.91246 0.84075 0.91435 0.92089 0.91081 0.91877 0.76778 0.88532 0.88533 0.91945 0.90903 0.91423 0.90840 0.91575 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. Note: The number of geometries may not correspond to the number of cycles due to rejected geometry changes. In the following table subtotal= G_P + SS G_CDS. Atom Chem. CM2 G_P Area Sigma k SS G_CDS Subtotal M number symbol chg. (kcal) (Ang**2) cal/(Ang**2) (kcal) (kcal) value 1 C -0.19 -0.70 9.53 37.16 0.35 -0.35 16 2 C 0.10 0.48 3.59 -26.73 -0.10 0.39 16 3 H 0.10 0.51 8.14 -51.93 -0.42 0.09 16 4 O -0.56 -3.44 12.88 -35.23 -0.45 -3.89 16 5 C -0.15 -0.65 4.89 -26.73 -0.13 -0.78 16 6 C -0.11 -0.64 5.37 -26.73 -0.14 -0.79 16 7 C 0.06 0.36 5.48 -3.83 -0.02 0.34 16 8 N -0.69 -4.85 6.25 -106.73 -0.67 -5.51 16 9 C 0.09 0.64 2.47 -67.71 -0.17 0.48 16 10 H 0.08 0.58 7.80 -51.93 -0.40 0.18 16 11 C -0.14 -1.22 5.17 -26.73 -0.14 -1.36 16 12 C -0.13 -1.31 5.49 -26.73 -0.15 -1.46 16 13 C 0.15 1.61 2.45 -67.93 -0.17 1.44 16 14 H 0.09 1.15 7.57 -51.93 -0.39 0.76 16 15 N -0.73 -9.16 5.08 -53.81 -0.27 -9.44 16 16 C 0.53 9.33 7.81 -10.99 -0.09 9.25 16 17 O -0.52 -11.10 15.52 5.56 0.09 -11.02 16 18 C -0.12 -2.40 6.16 36.00 0.22 -2.18 16 19 S -0.14 -3.61 18.55 -107.50 -1.99 -5.60 16 20 C -0.52 -15.07 10.04 30.94 0.31 -14.76 16 21 H 0.09 2.64 8.03 -52.48 -0.42 2.22 16 22 C 0.10 2.82 5.50 -17.43 -0.10 2.72 16 23 N -0.88 -26.60 13.97 55.49 0.78 -25.83 16 24 Si 0.86 21.83 27.74 -169.99 -4.72 17.12 16 25 H -0.27 -6.72 7.11 56.52 0.40 -6.32 16 26 H -0.27 -6.72 7.11 56.52 0.40 -6.31 16 27 H -0.24 -6.00 6.74 56.52 0.38 -5.62 16 28 C -0.13 -1.00 5.39 -26.71 -0.14 -1.15 16 29 C -0.10 -0.76 5.53 -26.71 -0.15 -0.91 16 30 H 0.06 0.18 8.14 -51.93 -0.42 -0.24 16 31 H 0.06 0.27 8.14 -51.93 -0.42 -0.16 16 32 H 0.07 0.22 8.14 -51.93 -0.42 -0.20 16 33 H 0.37 1.04 9.05 45.56 0.41 1.45 16 34 H 0.07 0.29 8.14 -51.93 -0.42 -0.13 16 35 H 0.07 0.26 8.14 -51.93 -0.42 -0.16 16 36 H 0.08 0.54 7.82 -51.93 -0.41 0.14 16 37 H 0.06 0.35 8.14 -51.93 -0.42 -0.07 16 38 H 0.03 0.18 7.44 -51.93 -0.39 -0.20 16 39 H 0.07 0.44 8.06 -51.93 -0.42 0.03 16 40 H 0.34 2.32 8.47 -39.29 -0.33 1.99 16 41 H 0.07 0.61 8.13 -51.93 -0.42 0.19 16 42 H 0.06 0.50 7.45 -51.93 -0.39 0.12 16 43 H 0.07 0.94 8.14 -51.93 -0.42 0.52 16 44 H 0.06 0.61 8.14 -51.93 -0.42 0.18 16 45 H 0.40 3.86 8.60 -40.82 -0.35 3.51 16 46 H 0.10 1.73 8.14 -51.92 -0.42 1.31 16 47 H 0.10 1.73 8.14 -51.91 -0.42 1.31 16 48 H 0.27 7.60 8.31 -40.82 -0.34 7.26 16 49 H 0.07 0.58 8.14 -51.93 -0.42 0.16 16 50 H 0.07 0.48 8.14 -51.93 -0.42 0.05 16 51 H 0.07 0.55 8.14 -51.93 -0.42 0.12 16 52 H 0.07 0.53 8.14 -51.93 -0.42 0.11 16 LS Contribution 424.69 15.07 6.40 6.40 Total: -1.00 -34.20 424.69 -10.45 -44.65 By element: Atomic # 1 Polarization: 11.27 SS G_CDS: -9.01 Total: 2.26 kcal Atomic # 6 Polarization: -8.53 SS G_CDS: -0.60 Total: -9.13 kcal Atomic # 7 Polarization: -40.62 SS G_CDS: -0.17 Total: -40.78 kcal Atomic # 8 Polarization: -14.54 SS G_CDS: -0.37 Total: -14.91 kcal Atomic # 14 Polarization: 21.83 SS G_CDS: -4.72 Total: 17.12 kcal Atomic # 16 Polarization: -3.61 SS G_CDS: -1.99 Total: -5.60 kcal Total LS contribution 6.40 Total: 6.40 kcal Total: -34.20 -10.45 -44.65 kcal The number of atoms in the molecule is 52 The average number of expansion shells was 16.00 The maximum number of expansion shells was 16 The minimum number of expansion shells was 16 **** NOTA BENE **** This is the net solvation energy for this exact molecular structure (nuclear and electronic)! The standard-state solvation energy should be obtained as the difference between the heat of formation plus delta-G solvation for the relaxed solvated system and that for the relaxed gas-phase system. ZINC001495442992.mol2 52 1SCF run This is a breakdown of the solvation energy calculated without geometric relaxation in solution: (1) E-EN(sol) electronic-nuclear energy of solute -82.675 kcal (2) G-P(sol) polarization free energy of solvation -34.201 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -116.876 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.447 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.647 kcal (6) G-S(sol) free energy of system = (1) + (5) -127.323 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.79 seconds