Wall clock time and date at job start Sun Mar 7 2021 05:16:52 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= 0 * AM1 - THE AM1 HAMILTONIAN TO BE USED * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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 O 1.45200 * 1 3 3 C 1.34705 * 116.99807 * 2 1 4 4 O 1.21469 * 120.00383 * 359.97438 * 3 2 1 5 5 C 1.47732 * 120.00155 * 180.02562 * 3 2 1 6 6 C 1.39607 * 120.15295 * 180.24960 * 5 3 2 7 7 C 1.37962 * 119.84885 * 180.26874 * 6 5 3 8 8 C 1.50690 * 119.92870 * 179.72992 * 7 6 5 9 9 H 1.09000 * 112.90879 * 348.48016 * 8 7 6 10 10 C 1.53786 * 113.61518 * 119.72668 * 8 7 6 11 11 C 1.53781 * 87.17350 * 140.02616 * 10 8 7 12 12 H 1.08998 * 113.61196 * 219.97759 * 11 10 8 13 13 N 1.46506 * 113.61613 * 89.15811 * 11 10 8 14 14 C 1.34774 * 119.99981 * 154.99882 * 13 11 10 15 15 O 1.21275 * 120.00090 * 0.02562 * 14 13 11 16 16 C 1.50696 * 119.99541 * 180.02562 * 14 13 11 17 17 N 1.46505 * 109.47028 * 180.02562 * 16 14 13 18 18 C 1.36828 * 125.06122 * 265.02489 * 17 16 14 19 19 C 1.34246 * 109.93680 * 179.97438 * 18 17 16 20 20 C 1.50697 * 126.49906 * 179.97438 * 19 18 17 21 21 O 1.42906 * 109.46729 * 255.02161 * 20 19 18 22 22 Si 1.86296 * 109.47328 * 15.01959 * 20 19 18 23 23 C 1.46491 * 106.99691 * 359.97438 * 19 18 17 24 24 C 1.39596 * 134.03437 * 180.02562 * 23 19 18 25 25 C 1.36662 * 119.75958 * 179.97438 * 24 23 19 26 26 C 1.38790 * 120.53304 * 0.02562 * 25 24 23 27 27 C 1.37785 * 120.66312 * 359.97438 * 26 25 24 28 28 C 1.37697 * 125.06059 * 84.71379 * 17 16 14 29 29 C 1.53781 * 87.08070 * 334.52830 * 11 10 8 30 30 C 1.38349 * 120.14739 * 359.46105 * 7 6 5 31 31 C 1.38384 * 120.30413 * 0.29385 * 30 7 6 32 32 C 1.37926 * 120.15042 * 359.97438 * 31 30 7 33 33 H 1.08994 * 109.47586 * 59.99976 * 1 2 3 34 34 H 1.09002 * 109.46989 * 179.97438 * 1 2 3 35 35 H 1.09000 * 109.47236 * 299.99703 * 1 2 3 36 36 H 1.08003 * 120.07620 * 0.02562 * 6 5 3 37 37 H 1.09002 * 113.61357 * 25.43533 * 10 8 7 38 38 H 1.09002 * 113.61608 * 254.61695 * 10 8 7 39 39 H 0.97003 * 119.99276 * 335.00814 * 13 11 10 40 40 H 1.09007 * 109.47140 * 299.99355 * 16 14 13 41 41 H 1.08993 * 109.47822 * 59.99828 * 16 14 13 42 42 H 1.07995 * 125.03441 * 359.95924 * 18 17 16 43 43 H 1.09005 * 109.47205 * 135.02660 * 20 19 18 44 44 H 0.96693 * 113.99951 * 60.00518 * 21 20 19 45 45 H 1.48505 * 109.99761 * 179.97438 * 22 20 19 46 46 H 1.48496 * 110.00265 * 301.32025 * 22 20 19 47 47 H 1.07992 * 120.12202 * 359.95062 * 24 23 19 48 48 H 1.08005 * 119.73099 * 179.97438 * 25 24 23 49 49 H 1.08002 * 119.66580 * 179.97438 * 26 25 24 50 50 H 1.08001 * 120.11077 * 180.02562 * 27 26 25 51 51 H 1.09005 * 113.61471 * 139.98321 * 29 11 10 52 52 H 1.09003 * 113.61469 * 270.88099 * 29 11 10 53 53 H 1.07999 * 119.85070 * 180.27628 * 30 7 6 54 54 H 1.08006 * 119.92620 * 180.02562 * 31 30 7 55 55 H 1.08000 * 120.07562 * 180.02562 * 32 31 30 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 8 1.4520 0.0000 0.0000 3 6 2.0635 1.2003 0.0000 4 8 1.4020 2.2190 -0.0005 5 6 3.5388 1.2777 -0.0006 6 6 4.1758 2.5199 0.0047 7 6 5.5539 2.5849 -0.0015 8 6 6.2434 3.9248 -0.0020 9 1 5.5618 4.7526 -0.1977 10 6 7.1521 4.1505 1.2180 11 6 8.1507 4.8151 0.2556 12 1 9.1907 4.5640 0.4637 13 7 7.9354 6.2545 0.0884 14 6 8.9585 7.0505 -0.2806 15 8 10.0572 6.5752 -0.4746 16 6 8.7370 8.5312 -0.4520 17 7 9.9930 9.1695 -0.8540 18 6 10.8827 9.7613 0.0006 19 6 11.9264 10.2474 -0.6897 20 6 13.1193 10.9707 -0.1199 21 8 13.0079 12.3676 -0.4001 22 14 13.1812 10.7102 1.7237 23 6 11.6857 9.9413 -2.1019 24 6 12.3928 10.1823 -3.2812 25 6 11.8845 9.7562 -4.4761 26 6 10.6702 9.0859 -4.5269 27 6 9.9559 8.8376 -3.3752 28 6 10.4567 9.2578 -2.1476 29 6 7.5225 3.9627 -0.8596 30 6 6.3044 1.4226 -0.0016 31 6 5.6794 0.1879 -0.0017 32 6 4.3024 0.1085 -0.0011 33 1 -0.3634 0.5138 -0.8899 34 1 -0.3633 -1.0277 -0.0005 35 1 -0.3634 0.5138 0.8900 36 1 3.5912 3.4280 0.0096 37 1 7.5230 3.2271 1.6629 38 1 6.7326 4.8363 1.9540 39 1 7.0565 6.6346 0.2433 40 1 7.9824 8.6992 -1.2205 41 1 8.3970 8.9594 0.4909 42 1 10.7629 9.8278 1.0718 43 1 14.0310 10.5809 -0.5728 44 1 12.2179 12.7818 -0.0267 45 1 14.3537 11.4181 2.2976 46 1 11.9301 11.2086 2.3494 47 1 13.3379 10.7039 -3.2495 48 1 12.4318 9.9432 -5.3882 49 1 10.2814 8.7557 -5.4789 50 1 9.0121 8.3149 -3.4250 51 1 7.4009 4.4973 -1.8018 52 1 7.9946 2.9880 -0.9828 53 1 7.3829 1.4791 -0.0020 54 1 6.2712 -0.7155 -0.0014 55 1 3.8159 -0.8558 -0.0007 There are 77 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 77 No. of singly occupied orbitals= 1 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=0 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=WATER s_22_15508394_9818708.mol2 55 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sun Mar 7 2021 05:16:52 Heat of formation + Delta-G solvation = -18.547554 kcal Electronic energy + Delta-G solvation = -38705.706263 eV Core-core repulsion = 33641.261033 eV Total energy + Delta-G solvation = -5064.445229 eV No. of doubly occupied orbitals = 77 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 421.177 amu Computer time = 1.84 seconds Orbital eigenvalues (eV) -42.54241 -41.81807 -41.30918 -40.18240 -38.77385 -38.14025 -37.65964 -37.56666 -34.73024 -34.02143 -32.40125 -32.14310 -31.33457 -30.23650 -29.02081 -28.67892 -26.68354 -26.39877 -25.61232 -25.05993 -23.77470 -22.95385 -22.78277 -22.02985 -21.72544 -20.68075 -20.20562 -18.91282 -18.61776 -18.41477 -18.22264 -17.72963 -17.47630 -17.34849 -17.06634 -16.80199 -16.41274 -16.30493 -16.20529 -16.05787 -15.89465 -15.47788 -15.38913 -14.97314 -14.84845 -14.71446 -14.66567 -14.63221 -14.48629 -14.28133 -14.06669 -14.03291 -13.64441 -13.48343 -13.27959 -13.15223 -13.00068 -12.86509 -12.72191 -12.68090 -12.56592 -12.39437 -11.99879 -11.90674 -11.84280 -11.79547 -11.73190 -11.69980 -11.22644 -11.12534 -10.54789 -10.44600 -10.13534 -10.05824 -9.74002 -9.15635 -8.63041 -4.85968 -0.38643 0.14569 0.30695 0.64487 1.35697 1.40613 1.41436 1.47564 1.58888 1.61443 1.68692 1.97925 2.20661 2.68509 2.95545 3.45394 3.54072 3.58453 3.61635 3.62820 3.66217 3.71263 3.77299 3.84665 3.88802 3.97472 4.04083 4.13069 4.14049 4.15424 4.22276 4.23235 4.26697 4.28310 4.30326 4.41816 4.50508 4.59434 4.64335 4.66384 4.72288 4.76261 4.79288 4.79912 4.86756 4.90446 4.97205 4.98563 5.04104 5.14241 5.14361 5.18603 5.26283 5.40114 5.47112 5.51353 5.63617 5.66921 5.74000 5.75683 5.83996 6.19416 6.23917 6.49133 6.55993 6.60700 7.16257 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.029 3.971 2 O -0.358 6.358 3 C 0.505 3.495 4 O -0.540 6.540 5 C -0.147 4.147 6 C -0.061 4.061 7 C -0.083 4.083 8 C -0.068 4.068 9 H 0.110 0.890 10 C -0.123 4.123 11 C 0.118 3.882 12 H 0.087 0.913 13 N -0.696 5.696 14 C 0.496 3.504 15 O -0.565 6.565 16 C 0.086 3.914 17 N -0.466 5.466 18 C 0.075 3.925 19 C -0.215 4.215 20 C -0.051 4.051 21 O -0.553 6.553 22 Si 0.766 3.234 23 C -0.116 4.116 24 C -0.088 4.088 25 C -0.154 4.154 26 C -0.107 4.107 27 C -0.134 4.134 28 C 0.092 3.908 29 C -0.142 4.142 30 C -0.077 4.077 31 C -0.129 4.129 32 C -0.069 4.069 33 H 0.062 0.938 34 H 0.133 0.867 35 H 0.063 0.937 36 H 0.131 0.869 37 H 0.104 0.896 38 H 0.089 0.911 39 H 0.419 0.581 40 H 0.141 0.859 41 H 0.161 0.839 42 H 0.197 0.803 43 H 0.104 0.896 44 H 0.383 0.617 45 H -0.302 1.302 46 H -0.260 1.260 47 H 0.119 0.881 48 H 0.136 0.864 49 H 0.144 0.856 50 H 0.132 0.868 51 H 0.067 0.933 52 H 0.100 0.900 53 H 0.152 0.848 54 H 0.156 0.844 55 H 0.146 0.854 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -4.946 -4.938 3.436 7.788 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.064 4.064 2 O -0.274 6.274 3 C 0.352 3.648 4 O -0.431 6.431 5 C -0.148 4.148 6 C -0.080 4.080 7 C -0.083 4.083 8 C -0.086 4.086 9 H 0.128 0.872 10 C -0.161 4.161 11 C 0.012 3.988 12 H 0.105 0.895 13 N -0.349 5.349 14 C 0.281 3.719 15 O -0.442 6.442 16 C -0.038 4.038 17 N -0.148 5.148 18 C -0.051 4.051 19 C -0.222 4.222 20 C -0.159 4.159 21 O -0.361 6.361 22 Si 0.668 3.332 23 C -0.120 4.120 24 C -0.106 4.106 25 C -0.173 4.173 26 C -0.126 4.126 27 C -0.154 4.154 28 C -0.013 4.013 29 C -0.180 4.180 30 C -0.095 4.095 31 C -0.146 4.146 32 C -0.087 4.087 33 H 0.081 0.919 34 H 0.152 0.848 35 H 0.081 0.919 36 H 0.149 0.851 37 H 0.122 0.878 38 H 0.108 0.892 39 H 0.257 0.743 40 H 0.159 0.841 41 H 0.179 0.821 42 H 0.214 0.786 43 H 0.122 0.878 44 H 0.232 0.768 45 H -0.230 1.230 46 H -0.185 1.185 47 H 0.137 0.863 48 H 0.154 0.846 49 H 0.162 0.838 50 H 0.150 0.850 51 H 0.085 0.915 52 H 0.118 0.882 53 H 0.169 0.831 54 H 0.174 0.826 55 H 0.163 0.837 Dipole moment (debyes) X Y Z Total from point charges -4.173 -4.879 3.887 7.505 hybrid contribution -1.275 1.128 -0.657 1.824 sum -5.447 -3.751 3.231 7.361 Atomic orbital electron populations 1.23581 0.74307 1.05921 1.02603 1.86532 1.25198 1.31918 1.83780 1.21647 0.89946 0.81394 0.71764 1.90818 1.67661 1.37024 1.47631 1.19187 0.88187 0.96150 1.11305 1.21219 0.94755 0.96369 0.95619 1.19761 0.91996 0.96124 1.00406 1.21663 0.95349 0.93849 0.97721 0.87219 1.23387 0.94871 1.01909 0.95953 1.21800 1.00053 0.79785 0.97117 0.89489 1.45515 1.13079 1.07127 1.69216 1.21334 0.84450 0.89860 0.76228 1.90804 1.26055 1.76684 1.50663 1.21975 0.85401 0.89039 1.07399 1.43432 1.19085 1.48145 1.04143 1.22012 0.89435 0.97050 0.96605 1.20703 0.99595 1.10087 0.91861 1.24585 0.97205 0.81065 1.13058 1.86562 1.49026 1.21758 1.78771 0.94050 0.77781 0.90239 0.71092 1.19338 0.95717 1.02952 0.93974 1.20451 0.99466 0.99336 0.91382 1.20914 0.97339 1.01059 0.97945 1.21123 0.94808 0.98474 0.98227 1.20235 1.01694 1.03924 0.89565 1.18385 0.95221 1.01619 0.86084 1.23618 0.95508 1.01153 0.97752 1.21651 1.01580 0.91528 0.94759 1.21497 0.93925 0.99290 0.99910 1.21662 0.94680 0.97772 0.94556 0.91943 0.84848 0.91865 0.85097 0.87809 0.89245 0.74320 0.84125 0.82105 0.78612 0.87842 0.76764 1.22950 1.18539 0.86346 0.84629 0.83775 0.85010 0.91452 0.88189 0.83067 0.82606 0.83651 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.03 -0.03 10.56 113.37 1.20 1.16 16 2 O -0.36 -1.35 9.82 -56.15 -0.55 -1.90 16 3 C 0.51 3.44 7.90 70.33 0.56 4.00 16 4 O -0.54 -6.55 16.14 19.94 0.32 -6.23 16 5 C -0.15 -0.41 5.87 -20.10 -0.12 -0.53 16 6 C -0.06 -0.03 9.25 22.53 0.21 0.18 16 7 C -0.08 0.33 5.29 -19.91 -0.11 0.22 16 8 C -0.07 0.54 3.87 -10.92 -0.04 0.50 16 9 H 0.11 -1.06 8.09 -2.39 -0.02 -1.08 16 10 C -0.12 1.16 7.28 31.12 0.23 1.39 16 11 C 0.12 -0.74 4.32 45.52 0.20 -0.55 16 12 H 0.09 -0.12 7.46 -2.39 -0.02 -0.14 16 13 N -0.70 5.89 5.27 -440.04 -2.32 3.57 16 14 C 0.50 -1.27 7.90 87.66 0.69 -0.58 16 15 O -0.56 -4.85 16.02 -3.01 -0.05 -4.90 16 16 C 0.09 -0.66 6.48 85.63 0.55 -0.11 16 17 N -0.47 -0.57 2.27 -519.99 -1.18 -1.75 16 18 C 0.07 0.07 9.89 83.00 0.82 0.89 16 19 C -0.22 -1.18 5.58 -18.62 -0.10 -1.28 16 20 C -0.05 -0.20 3.06 71.23 0.22 0.01 16 21 O -0.55 -3.66 12.88 -148.98 -1.92 -5.58 16 22 Si 0.77 -0.37 30.82 68.60 2.11 1.74 16 23 C -0.12 -1.11 6.35 -20.20 -0.13 -1.24 16 24 C -0.09 -0.83 10.05 22.23 0.22 -0.60 16 25 C -0.15 -1.01 9.97 22.03 0.22 -0.79 16 26 C -0.11 -0.52 10.01 22.29 0.22 -0.29 16 27 C -0.13 -0.67 10.00 22.37 0.22 -0.45 16 28 C 0.09 0.64 6.80 39.03 0.27 0.90 16 29 C -0.14 0.92 7.26 31.21 0.23 1.14 16 30 C -0.08 0.43 8.61 22.33 0.19 0.63 16 31 C -0.13 0.62 10.03 22.23 0.22 0.85 16 32 C -0.07 0.09 9.62 22.54 0.22 0.31 16 33 H 0.06 -0.03 8.14 -2.39 -0.02 -0.05 16 34 H 0.13 -0.98 8.14 -2.39 -0.02 -1.00 16 35 H 0.06 -0.04 8.14 -2.39 -0.02 -0.06 16 36 H 0.13 0.12 7.65 -2.91 -0.02 0.10 16 37 H 0.10 -1.00 7.72 -2.39 -0.02 -1.02 16 38 H 0.09 -0.98 8.14 -2.39 -0.02 -0.99 16 39 H 0.42 -6.27 8.70 -92.71 -0.81 -7.08 16 40 H 0.14 -1.71 8.14 -2.38 -0.02 -1.73 16 41 H 0.16 -2.49 8.12 -2.39 -0.02 -2.51 16 42 H 0.20 -1.15 6.27 -2.91 -0.02 -1.17 16 43 H 0.10 0.52 7.96 -2.38 -0.02 0.50 16 44 H 0.38 -0.01 8.70 -74.06 -0.64 -0.65 16 45 H -0.30 -2.17 7.11 99.48 0.71 -1.46 16 46 H -0.26 0.13 7.00 99.48 0.70 0.83 16 47 H 0.12 1.00 8.06 -2.91 -0.02 0.98 16 48 H 0.14 0.34 8.06 -2.91 -0.02 0.32 16 49 H 0.14 -0.03 8.06 -2.91 -0.02 -0.05 16 50 H 0.13 0.04 8.06 -2.91 -0.02 0.01 16 51 H 0.07 -0.29 8.14 -2.38 -0.02 -0.30 16 52 H 0.10 -0.71 6.31 -2.39 -0.02 -0.72 16 53 H 0.15 -1.48 6.65 -2.91 -0.02 -1.50 16 54 H 0.16 -1.48 8.06 -2.91 -0.02 -1.50 16 55 H 0.15 -0.50 7.61 -2.91 -0.02 -0.52 16 Total: 0.00 -32.24 463.64 2.13 -30.10 By element: Atomic # 1 Polarization: -20.34 SS G_CDS: -0.47 Total: -20.82 kcal Atomic # 6 Polarization: -0.42 SS G_CDS: 6.19 Total: 5.77 kcal Atomic # 7 Polarization: 5.31 SS G_CDS: -3.50 Total: 1.82 kcal Atomic # 8 Polarization: -16.41 SS G_CDS: -2.20 Total: -18.61 kcal Atomic # 14 Polarization: -0.37 SS G_CDS: 2.11 Total: 1.74 kcal Total: -32.24 2.13 -30.10 kcal The number of atoms in the molecule is 55 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. s_22_15508394_9818708.mol2 55 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 11.556 kcal (2) G-P(sol) polarization free energy of solvation -32.237 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -20.681 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 2.133 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -30.104 kcal (6) G-S(sol) free energy of system = (1) + (5) -18.548 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.84 seconds