Wall clock time and date at job start Sat Apr 11 2020 06:13:22 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.52996 * 1 3 3 H 1.09005 * 109.50222 * 2 1 4 4 C 1.53048 * 109.50091 * 239.90007 * 2 1 3 5 5 C 1.53032 * 109.54088 * 178.68199 * 4 2 1 6 6 C 1.53191 * 109.31120 * 61.36717 * 5 4 2 7 7 N 1.46930 * 108.77838 * 305.36195 * 6 5 4 8 8 C 1.34771 * 120.63623 * 233.63138 * 7 6 5 9 9 O 1.21277 * 120.00543 * 186.40752 * 8 7 6 10 10 C 1.50701 * 119.99838 * 6.40481 * 8 7 6 11 11 C 1.52998 * 109.47359 * 180.02562 * 10 8 7 12 12 O 1.42973 * 109.46671 * 294.99544 * 11 10 8 13 13 C 1.43087 * 113.76924 * 181.20460 * 12 11 10 14 14 C 1.53260 * 108.94368 * 57.97094 * 13 12 11 15 15 C 1.53254 * 108.84675 * 305.47754 * 14 13 12 16 16 O 1.42969 * 109.47184 * 175.00365 * 11 10 8 17 17 C 1.46929 * 118.73158 * 53.61029 * 7 6 5 18 18 H 1.08998 * 109.70754 * 66.24911 * 17 7 6 19 19 C 1.53003 * 109.58389 * 186.60680 * 17 7 6 20 20 C 1.52997 * 109.47433 * 65.83821 * 19 17 7 21 21 N 1.46498 * 109.47412 * 175.74376 * 20 19 17 22 22 C 1.36931 * 120.00067 * 180.02562 * 21 20 19 23 23 H 1.08005 * 120.00408 * 359.97438 * 22 21 20 24 24 C 1.38816 * 119.99847 * 180.02562 * 22 21 20 25 25 N 1.31617 * 120.00056 * 359.97438 * 24 22 21 26 26 Si 1.86807 * 119.99883 * 180.02562 * 24 22 21 27 27 H 1.48501 * 109.46888 * 359.97438 * 26 24 22 28 28 H 1.48495 * 109.47154 * 119.99757 * 26 24 22 29 29 H 1.48501 * 109.46568 * 239.99511 * 26 24 22 30 30 H 1.08993 * 109.47127 * 299.93667 * 1 2 3 31 31 H 1.09000 * 109.46756 * 59.93616 * 1 2 3 32 32 H 1.08997 * 109.47005 * 179.97438 * 1 2 3 33 33 H 1.08999 * 109.45594 * 298.70833 * 4 2 1 34 34 H 1.09001 * 109.45294 * 58.66534 * 4 2 1 35 35 H 1.08998 * 109.50331 * 181.31971 * 5 4 2 36 36 H 1.09004 * 109.50179 * 301.38521 * 5 4 2 37 37 H 1.08996 * 109.58899 * 65.15378 * 6 5 4 38 38 H 1.09007 * 109.58922 * 185.43604 * 6 5 4 39 39 H 1.09010 * 109.47302 * 299.99934 * 10 8 7 40 40 H 1.09001 * 109.47584 * 59.99505 * 10 8 7 41 41 H 1.08994 * 109.47259 * 54.99818 * 11 10 8 42 42 H 1.09005 * 109.56126 * 298.13708 * 13 12 11 43 43 H 1.08995 * 109.64803 * 177.86729 * 13 12 11 44 44 H 1.08999 * 109.57362 * 185.66922 * 14 13 12 45 45 H 1.09003 * 109.67515 * 65.35006 * 14 13 12 46 46 H 1.09001 * 109.56022 * 294.68307 * 15 14 13 47 47 H 1.08995 * 109.64866 * 174.41157 * 15 14 13 48 48 H 1.09005 * 109.47006 * 305.84110 * 19 17 7 49 49 H 1.08999 * 109.47465 * 185.84599 * 19 17 7 50 50 H 1.08998 * 109.47112 * 295.74495 * 20 19 17 51 51 H 1.08997 * 109.46800 * 55.74087 * 20 19 17 52 52 H 0.96998 * 119.99737 * 359.97438 * 21 20 19 53 53 H 0.97002 * 120.00020 * 180.02562 * 25 24 22 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.8939 1.0275 0.0000 4 6 2.0409 -0.7235 -1.2481 5 6 3.5708 -0.6954 -1.2658 6 6 4.1052 -1.4267 -0.0304 7 7 3.5046 -0.8288 1.1699 8 6 4.2811 -0.3849 2.1781 9 8 3.7737 -0.0009 3.2105 10 6 5.7799 -0.3659 2.0226 11 6 6.4186 0.1779 3.3021 12 8 6.0395 1.5449 3.4800 13 6 6.6110 2.1555 4.6411 14 6 8.1382 2.0647 4.5502 15 6 8.5356 0.5993 4.3421 16 8 7.8418 0.0876 3.2000 17 6 2.0413 -0.7235 1.2498 18 1 1.6032 -1.7203 1.2998 19 6 1.6517 0.0713 2.4978 20 6 2.1346 1.5159 2.3540 21 7 1.6739 2.3014 3.5015 22 6 1.9871 3.6314 3.5909 23 1 2.5741 4.1013 2.8155 24 6 1.5510 4.3756 4.6786 25 7 0.8362 3.8028 5.6237 26 14 1.9776 6.1902 4.8001 27 1 2.7836 6.5894 3.6185 28 1 0.7277 6.9909 4.8407 29 1 2.7611 6.4325 6.0381 30 1 -0.3633 0.5128 0.8905 31 1 -0.3633 0.5148 -0.8894 32 1 -0.3633 -1.0276 -0.0005 33 1 1.6967 -1.7577 -1.2347 34 1 1.6582 -0.2251 -2.1388 35 1 3.9335 -1.1903 -2.1667 36 1 3.9159 0.3385 -1.2548 37 1 3.8376 -2.4818 -0.0860 38 1 5.1899 -1.3267 0.0128 39 1 6.1392 -1.3787 1.8397 40 1 6.0500 0.2726 1.1814 41 1 6.0776 -0.4078 4.1557 42 1 6.2670 1.6340 5.5344 43 1 6.3099 3.2018 4.6907 44 1 8.5817 2.4356 5.4743 45 1 8.4920 2.6620 3.7099 46 1 8.2645 0.0189 5.2240 47 1 9.6107 0.5310 4.1767 48 1 2.1140 -0.3807 3.3754 49 1 0.5677 0.0601 2.6120 50 1 3.2238 1.5327 2.3158 51 1 1.7322 1.9437 1.4357 52 1 1.1465 1.8795 4.1977 53 1 0.5311 4.3229 6.3836 RHF calculation, no. of doubly occupied orbitals= 66 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). 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 ZINC001488993178.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Apr 11 2020 06:13:22 Heat of formation + Delta-G solvation = -90.813824 kcal Electronic energy + Delta-G solvation = -32993.726908 eV Core-core repulsion = 28866.826667 eV Total energy + Delta-G solvation = -4126.900242 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 6.38 seconds Orbital eigenvalues (eV) -40.81131 -39.67916 -37.51871 -35.47932 -35.16587 -33.64001 -33.09897 -32.55077 -30.74314 -30.19530 -28.14204 -26.37575 -25.86608 -25.22130 -23.62528 -22.96012 -21.95947 -21.17042 -21.07611 -19.77577 -18.62232 -17.44799 -16.85930 -16.59184 -16.40349 -15.86794 -15.64031 -15.28104 -14.89135 -14.74053 -14.59800 -14.52501 -14.10210 -13.79447 -13.42488 -13.17848 -13.13437 -13.02977 -12.97470 -12.65525 -12.42973 -12.42675 -11.88295 -11.79254 -11.63863 -11.47031 -11.25375 -11.03603 -11.01410 -10.98842 -10.84958 -10.73772 -10.60320 -10.40122 -10.12265 -10.06152 -9.88983 -9.70940 -9.49623 -9.44791 -8.63494 -8.48257 -8.19286 -6.54205 -5.52333 -2.75661 2.88345 3.01188 3.54503 3.59740 3.62504 3.67447 3.85203 4.60980 4.68887 4.74297 4.77422 4.79291 4.82168 4.87372 4.94201 5.01662 5.08851 5.15517 5.25906 5.31644 5.33290 5.43640 5.49760 5.55449 5.61779 5.66148 5.68388 5.71793 5.72281 5.77473 5.80337 5.85454 5.93390 6.02098 6.08547 6.11656 6.12894 6.24624 6.40215 6.43770 6.59984 6.67246 6.78274 6.88809 6.88842 7.32553 7.36464 7.45642 7.66288 7.92115 8.57511 8.68076 9.77957 10.44412 10.78668 11.71800 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.006806 B = 0.005955 C = 0.003458 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 4113.051934 B = 4700.973077 C = 8094.413614 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.144 4.144 2 C -0.110 4.110 3 H 0.087 0.913 4 C -0.113 4.113 5 C -0.136 4.136 6 C 0.108 3.892 7 N -0.612 5.612 8 C 0.508 3.492 9 O -0.475 6.475 10 C -0.126 4.126 11 C 0.263 3.737 12 O -0.364 6.364 13 C 0.061 3.939 14 C -0.193 4.193 15 C 0.062 3.938 16 O -0.401 6.401 17 C 0.163 3.837 18 H 0.064 0.936 19 C -0.139 4.139 20 C 0.162 3.838 21 N -0.724 5.724 22 C -0.212 4.212 23 H 0.065 0.935 24 C -0.023 4.023 25 N -0.938 5.938 26 Si 0.899 3.101 27 H -0.283 1.283 28 H -0.295 1.295 29 H -0.295 1.295 30 H 0.067 0.933 31 H 0.053 0.947 32 H 0.050 0.950 33 H 0.060 0.940 34 H 0.062 0.938 35 H 0.071 0.929 36 H 0.069 0.931 37 H 0.066 0.934 38 H 0.082 0.918 39 H 0.107 0.893 40 H 0.104 0.896 41 H 0.062 0.938 42 H 0.054 0.946 43 H 0.102 0.898 44 H 0.082 0.918 45 H 0.079 0.921 46 H 0.052 0.948 47 H 0.096 0.904 48 H 0.108 0.892 49 H 0.052 0.948 50 H 0.038 0.962 51 H -0.001 1.001 52 H 0.380 0.620 53 H 0.250 0.750 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 12.163 -14.540 -11.678 22.265 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.202 4.202 2 C -0.129 4.129 3 H 0.105 0.895 4 C -0.151 4.151 5 C -0.175 4.175 6 C -0.015 4.015 7 N -0.348 5.348 8 C 0.295 3.705 9 O -0.348 6.348 10 C -0.167 4.167 11 C 0.167 3.833 12 O -0.284 6.284 13 C -0.015 4.015 14 C -0.231 4.231 15 C -0.015 4.015 16 O -0.322 6.322 17 C 0.060 3.940 18 H 0.082 0.918 19 C -0.177 4.177 20 C 0.035 3.965 21 N -0.367 5.367 22 C -0.344 4.344 23 H 0.083 0.917 24 C -0.217 4.217 25 N -0.589 5.589 26 Si 0.731 3.269 27 H -0.208 1.208 28 H -0.222 1.222 29 H -0.222 1.222 30 H 0.086 0.914 31 H 0.072 0.928 32 H 0.069 0.931 33 H 0.079 0.921 34 H 0.081 0.919 35 H 0.089 0.911 36 H 0.087 0.913 37 H 0.084 0.916 38 H 0.101 0.899 39 H 0.125 0.875 40 H 0.123 0.877 41 H 0.080 0.920 42 H 0.073 0.927 43 H 0.120 0.880 44 H 0.101 0.899 45 H 0.098 0.902 46 H 0.071 0.929 47 H 0.114 0.886 48 H 0.126 0.874 49 H 0.071 0.929 50 H 0.055 0.945 51 H 0.017 0.983 52 H 0.213 0.787 53 H 0.060 0.940 Dipole moment (debyes) X Y Z Total from point charges 11.934 -14.146 -11.433 21.754 hybrid contribution 0.183 0.996 0.589 1.172 sum 12.117 -13.150 -10.844 20.912 Atomic orbital electron populations 1.21781 0.94483 1.01376 1.02555 1.21679 0.96188 1.01143 0.93935 0.89504 1.21587 0.94871 1.00543 0.98090 1.21927 0.96594 1.01597 0.97356 1.21689 0.98201 0.96859 0.84739 1.47974 1.06664 1.62109 1.18076 1.20034 0.88117 0.78172 0.84187 1.90872 1.72915 1.43162 1.27890 1.21514 0.93429 1.03214 0.98494 1.22844 0.78878 0.84722 0.96860 1.88271 1.69913 1.24958 1.45233 1.22588 0.90610 1.00052 0.88286 1.22827 0.98617 0.98008 1.03688 1.22592 0.97202 0.93205 0.88485 1.88002 1.19804 1.82289 1.42142 1.21256 0.78248 0.98387 0.96088 0.91776 1.22453 0.99787 0.97771 0.97656 1.20974 0.97462 0.88216 0.89835 1.42249 1.63746 1.07179 1.23496 1.19256 1.25505 0.85934 1.03750 0.91747 1.24641 1.01846 0.97640 0.97620 1.69677 1.42493 1.34782 1.11986 0.86204 0.77358 0.86246 0.77082 1.20807 1.22218 1.22224 0.91423 0.92793 0.93113 0.92148 0.91900 0.91059 0.91266 0.91597 0.89942 0.87511 0.87733 0.92039 0.92726 0.87957 0.89929 0.90199 0.92930 0.88627 0.87424 0.92947 0.94487 0.98309 0.78715 0.94036 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.14 -1.81 8.97 37.16 0.33 -1.47 16 2 C -0.11 -1.43 1.99 -90.50 -0.18 -1.61 16 3 H 0.09 1.37 5.93 -51.93 -0.31 1.06 16 4 C -0.11 -1.11 5.17 -26.69 -0.14 -1.25 16 5 C -0.14 -1.20 6.08 -26.62 -0.16 -1.36 16 6 C 0.11 0.98 6.34 -3.71 -0.02 0.96 16 7 N -0.61 -8.01 2.94 -165.22 -0.49 -8.50 16 8 C 0.51 7.90 6.17 -10.99 -0.07 7.83 16 9 O -0.48 -9.52 9.28 -8.85 -0.08 -9.60 16 10 C -0.13 -1.62 4.14 -27.88 -0.12 -1.73 16 11 C 0.26 3.94 2.80 37.16 0.10 4.04 16 12 O -0.36 -6.69 9.77 -64.76 -0.63 -7.32 16 13 C 0.06 0.96 6.96 37.29 0.26 1.22 16 14 C -0.19 -2.27 6.35 -26.47 -0.17 -2.44 16 15 C 0.06 0.68 6.96 37.29 0.26 0.94 16 16 O -0.40 -5.42 10.54 -59.03 -0.62 -6.04 16 17 C 0.16 2.32 2.66 -67.60 -0.18 2.14 16 18 H 0.06 0.83 8.14 -51.93 -0.42 0.41 16 19 C -0.14 -2.62 4.20 -26.73 -0.11 -2.73 16 20 C 0.16 3.71 4.94 -4.04 -0.02 3.69 16 21 N -0.72 -20.34 5.60 -59.91 -0.34 -20.68 16 22 C -0.21 -6.46 9.99 -15.06 -0.15 -6.61 16 23 H 0.06 1.94 7.83 -52.48 -0.41 1.53 16 24 C -0.02 -0.71 5.50 -17.43 -0.10 -0.80 16 25 N -0.94 -30.27 13.06 55.49 0.72 -29.54 16 26 Si 0.90 23.84 29.55 -169.99 -5.02 18.82 16 27 H -0.28 -7.28 7.11 56.52 0.40 -6.88 16 28 H -0.29 -7.80 7.11 56.52 0.40 -7.40 16 29 H -0.29 -7.79 7.11 56.52 0.40 -7.39 16 30 H 0.07 1.01 6.75 -51.93 -0.35 0.65 16 31 H 0.05 0.62 8.14 -51.93 -0.42 0.20 16 32 H 0.05 0.55 8.14 -51.93 -0.42 0.13 16 33 H 0.06 0.54 8.14 -51.93 -0.42 0.11 16 34 H 0.06 0.58 8.14 -51.93 -0.42 0.16 16 35 H 0.07 0.48 8.14 -51.93 -0.42 0.05 16 36 H 0.07 0.70 8.14 -51.93 -0.42 0.28 16 37 H 0.07 0.50 8.14 -51.93 -0.42 0.08 16 38 H 0.08 0.61 5.94 -51.93 -0.31 0.31 16 39 H 0.11 1.02 7.68 -51.92 -0.40 0.62 16 40 H 0.10 1.28 7.95 -51.93 -0.41 0.87 16 41 H 0.06 0.96 7.95 -51.93 -0.41 0.55 16 42 H 0.05 0.85 8.14 -51.93 -0.42 0.43 16 43 H 0.10 1.68 8.14 -51.93 -0.42 1.25 16 44 H 0.08 0.78 8.14 -51.93 -0.42 0.36 16 45 H 0.08 0.95 8.14 -51.93 -0.42 0.53 16 46 H 0.05 0.53 8.14 -51.93 -0.42 0.11 16 47 H 0.10 0.82 8.14 -51.93 -0.42 0.40 16 48 H 0.11 2.25 5.58 -53.72 -0.30 1.95 16 49 H 0.05 0.99 6.93 -51.93 -0.36 0.63 16 50 H 0.04 0.91 5.04 -51.93 -0.26 0.65 16 51 H 0.00 -0.01 5.97 -51.93 -0.31 -0.32 16 52 H 0.38 11.32 7.67 -40.82 -0.31 11.01 16 53 H 0.25 7.76 8.31 -40.82 -0.34 7.42 16 LS Contribution 394.81 15.07 5.95 5.95 Total: -1.00 -36.18 394.81 -10.16 -46.34 By element: Atomic # 1 Polarization: 18.96 SS G_CDS: -9.20 Total: 9.76 kcal Atomic # 6 Polarization: 1.27 SS G_CDS: -0.46 Total: 0.82 kcal Atomic # 7 Polarization: -58.62 SS G_CDS: -0.10 Total: -58.72 kcal Atomic # 8 Polarization: -21.63 SS G_CDS: -1.34 Total: -22.97 kcal Atomic # 14 Polarization: 23.84 SS G_CDS: -5.02 Total: 18.82 kcal Total LS contribution 5.95 Total: 5.95 kcal Total: -36.18 -10.16 -46.34 kcal The number of atoms in the molecule is 53 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. ZINC001488993178.mol2 53 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 -44.473 kcal (2) G-P(sol) polarization free energy of solvation -36.175 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -80.649 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.165 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -46.340 kcal (6) G-S(sol) free energy of system = (1) + (5) -90.814 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 6.38 seconds