Wall clock time and date at job start Tue Mar 31 2020 05:46:18 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.53007 * 1 3 3 C 1.52995 * 109.47104 * 2 1 4 4 H 1.09005 * 109.47155 * 298.14978 * 3 2 1 5 5 C 1.50705 * 109.46970 * 178.15570 * 3 2 1 6 6 H 1.08001 * 119.99429 * 246.80980 * 5 3 2 7 7 C 1.37875 * 120.00384 * 66.81008 * 5 3 2 8 8 N 1.31519 * 119.99613 * 9.83518 * 7 5 3 9 9 Si 1.86794 * 120.00194 * 189.84227 * 7 5 3 10 10 H 1.48493 * 109.47052 * 269.99701 * 9 7 5 11 11 H 1.48504 * 109.47244 * 29.99820 * 9 7 5 12 12 H 1.48504 * 109.47297 * 149.99629 * 9 7 5 13 13 N 1.46499 * 109.47080 * 58.15356 * 3 2 1 14 14 C 1.34770 * 120.00132 * 261.46907 * 13 3 2 15 15 O 1.21280 * 120.00048 * 3.53902 * 14 13 3 16 16 C 1.50701 * 120.00411 * 183.53169 * 14 13 3 17 17 C 1.52994 * 109.47526 * 185.55802 * 16 14 13 18 18 C 1.53002 * 115.54966 * 280.72313 * 17 16 14 19 19 N 1.46496 * 109.46831 * 300.06804 * 18 17 16 20 20 C 1.34781 * 119.99607 * 182.93802 * 19 18 17 21 21 O 1.21511 * 120.00071 * 5.08452 * 20 19 18 22 22 O 1.34642 * 119.99959 * 185.08623 * 20 19 18 23 23 C 1.45193 * 117.00030 * 180.02562 * 22 20 19 24 24 C 1.52995 * 109.47344 * 59.99733 * 23 22 20 25 25 C 1.53002 * 109.47352 * 179.97438 * 23 22 20 26 26 C 1.53006 * 109.47129 * 299.99926 * 23 22 20 27 27 C 1.53000 * 117.50039 * 135.01477 * 17 16 14 28 28 C 1.53003 * 117.49861 * 66.39126 * 17 16 14 29 29 H 1.08996 * 109.46887 * 54.00726 * 1 2 3 30 30 H 1.08997 * 109.47044 * 174.01217 * 1 2 3 31 31 H 1.09001 * 109.46903 * 294.00792 * 1 2 3 32 32 H 1.09006 * 109.46506 * 239.99934 * 2 1 3 33 33 H 1.09000 * 109.47398 * 120.00390 * 2 1 3 34 34 H 0.96992 * 119.99969 * 356.77268 * 8 7 5 35 35 H 0.97001 * 119.99538 * 81.47492 * 13 3 2 36 36 H 1.09002 * 109.46808 * 305.55751 * 16 14 13 37 37 H 1.08996 * 109.46711 * 65.55524 * 16 14 13 38 38 H 1.09000 * 109.47133 * 60.06686 * 18 17 16 39 39 H 1.09002 * 109.46855 * 180.02562 * 18 17 16 40 40 H 0.96995 * 120.00307 * 2.94153 * 19 18 17 41 41 H 1.09002 * 109.46989 * 60.00180 * 24 23 22 42 42 H 1.09004 * 109.47532 * 179.97438 * 24 23 22 43 43 H 1.08997 * 109.47602 * 299.99997 * 24 23 22 44 44 H 1.08998 * 109.46593 * 59.99599 * 25 23 22 45 45 H 1.08992 * 109.47476 * 179.97438 * 25 23 22 46 46 H 1.08999 * 109.46849 * 300.00605 * 25 23 22 47 47 H 1.09001 * 109.47085 * 59.99725 * 26 23 22 48 48 H 1.08992 * 109.47158 * 180.02562 * 26 23 22 49 49 H 1.08996 * 109.46814 * 300.00114 * 26 23 22 50 50 H 1.09002 * 117.49735 * 359.97438 * 27 17 16 51 51 H 1.08998 * 117.49976 * 145.02314 * 27 17 16 52 52 H 1.08999 * 117.50237 * 214.98059 * 28 17 16 53 53 H 1.08994 * 117.49786 * 359.97438 * 28 17 16 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.5301 0.0000 0.0000 3 6 2.0400 1.4425 0.0000 4 1 1.7040 1.9466 -0.9061 5 6 3.5464 1.4427 0.0457 6 1 4.0600 1.8203 0.9175 7 6 4.2688 0.9608 -1.0253 8 7 3.6598 0.6824 -2.1572 9 14 6.1122 0.7033 -0.8685 10 1 6.8212 1.9455 -1.2678 11 1 6.4479 0.3727 0.5398 12 1 6.5331 -0.4130 -1.7528 13 7 1.5157 2.1458 1.1733 14 6 0.4333 2.9399 1.0549 15 8 -0.1452 3.0221 -0.0078 16 6 -0.0586 3.7255 2.2432 17 6 -1.3683 4.4292 1.8827 18 6 -1.2236 5.7143 1.0649 19 7 -0.4260 6.6851 1.8182 20 6 -0.2164 7.9186 1.3169 21 8 -0.7640 8.2567 0.2862 22 8 0.5998 8.7761 1.9582 23 6 0.7785 10.0875 1.3612 24 6 -0.5746 10.7962 1.2725 25 6 1.7317 10.9157 2.2252 26 6 1.3661 9.9328 -0.0430 27 6 -2.5089 4.3637 2.9004 28 6 -2.5900 3.5502 1.6070 29 1 -0.3633 0.6039 -0.8315 30 1 -0.3633 -1.0220 -0.1072 31 1 -0.3633 0.4181 0.9388 32 1 1.8933 -0.5139 -0.8901 33 1 1.8935 -0.5139 0.8899 34 1 2.6963 0.7714 -2.2242 35 1 1.9477 2.0399 2.0353 36 1 -0.2282 3.0490 3.0810 37 1 0.6885 4.4686 2.5218 38 1 -0.7268 5.4902 0.1210 39 1 -2.2109 6.1310 0.8653 40 1 -0.0487 6.4402 2.6775 41 1 -0.9932 10.9064 2.2729 42 1 -0.4404 11.7809 0.8247 43 1 -1.2540 10.2064 0.6572 44 1 2.6957 10.4109 2.2887 45 1 1.8662 11.8999 1.7768 46 1 1.3127 11.0257 3.2253 47 1 0.6867 9.3430 -0.6583 48 1 1.5007 10.9173 -0.4910 49 1 2.3301 9.4280 0.0202 50 1 -2.3263 3.8188 3.8266 51 1 -3.1564 5.2362 2.9876 52 1 -3.2908 3.8876 0.8434 53 1 -2.4604 2.4706 1.6821 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 ZINC000408521737.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 05:46:18 Heat of formation + Delta-G solvation = -115.921565 kcal Electronic energy + Delta-G solvation = -31077.377487 eV Core-core repulsion = 26949.388492 eV Total energy + Delta-G solvation = -4127.988995 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.85 seconds Orbital eigenvalues (eV) -40.83269 -39.65481 -38.02886 -36.99496 -35.27890 -34.55759 -33.71472 -32.46953 -31.22247 -29.27083 -27.60407 -27.52028 -27.20397 -26.19563 -23.14909 -22.66264 -22.19998 -20.94331 -20.17839 -19.69971 -19.10335 -18.12992 -17.37499 -16.68983 -16.33623 -16.18206 -15.83082 -15.16020 -14.91792 -14.69143 -14.46833 -14.42015 -13.91004 -13.80081 -13.53950 -13.35159 -13.03296 -12.90111 -12.72196 -12.67969 -12.30672 -12.11704 -12.07829 -12.04275 -11.80715 -11.59191 -11.57277 -11.46616 -11.40399 -11.20982 -10.77487 -10.61474 -10.49886 -10.45880 -10.38738 -10.24288 -10.18863 -10.00437 -9.93658 -9.19501 -8.97869 -8.90667 -8.31906 -7.90142 -6.45662 -4.18725 2.09692 2.83601 3.34047 3.38524 3.41647 3.53834 3.57066 4.09563 4.20985 4.46831 4.51093 4.54248 4.58478 4.70844 4.72974 4.78037 4.91824 5.26597 5.27436 5.34774 5.37868 5.45758 5.48909 5.53883 5.55463 5.58198 5.61008 5.70003 5.74252 5.79635 5.91535 5.96634 6.02022 6.05131 6.09704 6.09962 6.11281 6.16779 6.30354 6.32847 6.47427 6.71150 6.85490 7.02588 7.13834 7.39142 7.40646 7.53884 7.86034 8.32404 8.35848 8.52964 8.75877 9.06209 9.53265 10.92165 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.009829 B = 0.003036 C = 0.002469 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2847.928525 B = 9220.390069 C =11335.967204 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.141 4.141 2 C -0.129 4.129 3 C 0.256 3.744 4 H 0.048 0.952 5 C -0.617 4.617 6 H 0.054 0.946 7 C 0.027 3.973 8 N -0.902 5.902 9 Si 1.072 2.928 10 H -0.288 1.288 11 H -0.288 1.288 12 H -0.283 1.283 13 N -0.706 5.706 14 C 0.503 3.497 15 O -0.544 6.544 16 C -0.109 4.109 17 C -0.121 4.121 18 C 0.178 3.822 19 N -0.706 5.706 20 C 0.650 3.350 21 O -0.568 6.568 22 O -0.367 6.367 23 C 0.136 3.864 24 C -0.183 4.183 25 C -0.140 4.140 26 C -0.185 4.185 27 C -0.175 4.175 28 C -0.138 4.138 29 H 0.054 0.946 30 H 0.032 0.968 31 H 0.036 0.964 32 H 0.067 0.933 33 H 0.045 0.955 34 H 0.274 0.726 35 H 0.388 0.612 36 H 0.089 0.911 37 H 0.090 0.910 38 H 0.086 0.914 39 H 0.071 0.929 40 H 0.404 0.596 41 H 0.058 0.942 42 H 0.069 0.931 43 H 0.088 0.912 44 H 0.067 0.933 45 H 0.079 0.921 46 H 0.065 0.935 47 H 0.089 0.911 48 H 0.069 0.931 49 H 0.061 0.939 50 H 0.096 0.904 51 H 0.095 0.905 52 H 0.097 0.903 53 H 0.098 0.902 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -9.040 14.555 13.294 21.686 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.198 4.198 2 C -0.168 4.168 3 C 0.154 3.846 4 H 0.066 0.934 5 C -0.655 4.655 6 H 0.072 0.928 7 C -0.179 4.179 8 N -0.537 5.537 9 Si 0.900 3.100 10 H -0.215 1.215 11 H -0.214 1.214 12 H -0.209 1.209 13 N -0.354 5.354 14 C 0.287 3.713 15 O -0.422 6.422 16 C -0.149 4.149 17 C -0.122 4.122 18 C 0.056 3.944 19 N -0.365 5.365 20 C 0.403 3.597 21 O -0.457 6.457 22 O -0.281 6.281 23 C 0.100 3.900 24 C -0.240 4.240 25 C -0.198 4.198 26 C -0.243 4.243 27 C -0.212 4.212 28 C -0.175 4.175 29 H 0.073 0.927 30 H 0.051 0.949 31 H 0.055 0.945 32 H 0.086 0.914 33 H 0.063 0.937 34 H 0.086 0.914 35 H 0.220 0.780 36 H 0.108 0.892 37 H 0.109 0.891 38 H 0.104 0.896 39 H 0.090 0.910 40 H 0.239 0.761 41 H 0.077 0.923 42 H 0.088 0.912 43 H 0.106 0.894 44 H 0.086 0.914 45 H 0.098 0.902 46 H 0.084 0.916 47 H 0.107 0.893 48 H 0.088 0.912 49 H 0.080 0.920 50 H 0.114 0.886 51 H 0.114 0.886 52 H 0.116 0.884 53 H 0.116 0.884 Dipole moment (debyes) X Y Z Total from point charges -8.313 14.379 11.113 19.984 hybrid contribution -0.842 0.288 1.706 1.924 sum -9.155 14.668 12.819 21.524 Atomic orbital electron populations 1.21763 0.97977 0.99184 1.00916 1.21986 0.93111 1.00281 1.01397 1.18852 0.93445 0.87687 0.84653 0.93373 1.22072 0.90374 1.47178 1.05887 0.92807 1.25758 1.06797 0.93581 0.91801 1.70748 1.12475 1.45592 1.24910 0.83511 0.74134 0.76999 0.75386 1.21472 1.21416 1.20856 1.46234 1.28886 1.46709 1.13525 1.20293 0.81111 0.80231 0.89631 1.90652 1.56795 1.65562 1.29171 1.21258 0.96153 0.99459 0.98063 1.21130 0.90589 0.97502 1.02968 1.20819 0.95298 0.84110 0.94189 1.44732 1.56495 1.10581 1.24737 1.18059 0.78372 0.81434 0.81862 1.90900 1.53069 1.76970 1.24798 1.86293 1.62399 1.25781 1.53671 1.22045 0.96182 0.78279 0.93527 1.22511 0.97143 1.01318 1.03007 1.21858 0.99021 0.98978 0.99908 1.22530 1.02269 1.03425 0.96033 1.23205 1.00898 1.03043 0.94097 1.22833 0.98291 0.98147 0.98273 0.92737 0.94920 0.94467 0.91409 0.93655 0.91414 0.78050 0.89227 0.89132 0.89605 0.91031 0.76128 0.92318 0.91227 0.89355 0.91369 0.90206 0.91576 0.89256 0.91221 0.91975 0.88585 0.88631 0.88435 0.88380 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 -2.81 8.76 37.16 0.33 -2.48 16 2 C -0.13 -2.91 5.37 -26.73 -0.14 -3.06 16 3 C 0.26 6.22 1.90 -69.08 -0.13 6.09 16 4 H 0.05 1.37 6.50 -51.93 -0.34 1.03 16 5 C -0.62 -16.40 8.68 -34.44 -0.30 -16.70 16 6 H 0.05 1.40 7.83 -52.49 -0.41 0.99 16 7 C 0.03 0.74 5.22 -17.82 -0.09 0.64 16 8 N -0.90 -26.55 13.13 55.43 0.73 -25.82 16 9 Si 1.07 24.36 29.91 -169.99 -5.09 19.28 16 10 H -0.29 -6.42 7.11 56.52 0.40 -6.01 16 11 H -0.29 -6.17 7.11 56.52 0.40 -5.77 16 12 H -0.28 -6.44 7.11 56.52 0.40 -6.04 16 13 N -0.71 -13.83 4.55 -54.93 -0.25 -14.08 16 14 C 0.50 9.10 6.54 -10.99 -0.07 9.02 16 15 O -0.54 -11.71 13.30 5.56 0.07 -11.63 16 16 C -0.11 -1.37 5.55 -27.89 -0.15 -1.52 16 17 C -0.12 -1.43 1.10 -154.51 -0.17 -1.60 16 18 C 0.18 2.37 5.46 -4.04 -0.02 2.35 16 19 N -0.71 -8.71 5.33 -65.28 -0.35 -9.06 16 20 C 0.65 9.24 8.05 54.06 0.44 9.68 16 21 O -0.57 -9.37 11.88 -19.28 -0.23 -9.60 16 22 O -0.37 -4.68 9.94 -40.33 -0.40 -5.08 16 23 C 0.14 1.47 1.13 -90.62 -0.10 1.37 16 24 C -0.18 -1.83 8.37 37.16 0.31 -1.52 16 25 C -0.14 -1.11 8.85 37.16 0.33 -0.78 16 26 C -0.19 -2.12 8.37 37.16 0.31 -1.81 16 27 C -0.18 -1.61 9.89 -26.73 -0.26 -1.88 16 28 C -0.14 -1.60 9.26 -26.73 -0.25 -1.85 16 29 H 0.05 1.17 7.65 -51.93 -0.40 0.77 16 30 H 0.03 0.58 8.14 -51.93 -0.42 0.15 16 31 H 0.04 0.69 7.12 -51.93 -0.37 0.32 16 32 H 0.07 1.69 6.60 -51.93 -0.34 1.35 16 33 H 0.04 0.99 8.14 -51.93 -0.42 0.56 16 34 H 0.27 8.14 6.05 -40.82 -0.25 7.90 16 35 H 0.39 6.62 8.74 -40.82 -0.36 6.27 16 36 H 0.09 0.92 7.96 -51.93 -0.41 0.51 16 37 H 0.09 1.03 7.06 -51.93 -0.37 0.66 16 38 H 0.09 1.46 7.48 -51.93 -0.39 1.07 16 39 H 0.07 0.93 8.10 -51.93 -0.42 0.51 16 40 H 0.40 4.09 7.69 -40.82 -0.31 3.78 16 41 H 0.06 0.54 8.14 -51.93 -0.42 0.12 16 42 H 0.07 0.56 8.14 -51.93 -0.42 0.14 16 43 H 0.09 1.10 5.88 -51.93 -0.31 0.80 16 44 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16 45 H 0.08 0.48 8.14 -51.93 -0.42 0.06 16 46 H 0.07 0.48 8.14 -51.93 -0.42 0.06 16 47 H 0.09 1.26 5.88 -51.93 -0.31 0.96 16 48 H 0.07 0.64 8.14 -51.93 -0.42 0.21 16 49 H 0.06 0.72 8.14 -51.93 -0.42 0.30 16 50 H 0.10 0.76 8.09 -51.93 -0.42 0.34 16 51 H 0.10 0.81 8.14 -51.93 -0.42 0.39 16 52 H 0.10 1.11 8.14 -51.93 -0.42 0.69 16 53 H 0.10 1.21 8.14 -51.93 -0.42 0.79 16 LS Contribution 418.22 15.07 6.30 6.30 Total: -1.00 -32.25 418.22 -8.46 -40.71 By element: Atomic # 1 Polarization: 22.28 SS G_CDS: -9.26 Total: 13.02 kcal Atomic # 6 Polarization: -4.05 SS G_CDS: 0.01 Total: -4.04 kcal Atomic # 7 Polarization: -49.09 SS G_CDS: 0.13 Total: -48.95 kcal Atomic # 8 Polarization: -25.76 SS G_CDS: -0.56 Total: -26.32 kcal Atomic # 14 Polarization: 24.36 SS G_CDS: -5.09 Total: 19.28 kcal Total LS contribution 6.30 Total: 6.30 kcal Total: -32.25 -8.46 -40.71 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. ZINC000408521737.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 -75.212 kcal (2) G-P(sol) polarization free energy of solvation -32.248 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -107.461 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.461 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.709 kcal (6) G-S(sol) free energy of system = (1) + (5) -115.922 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.85 seconds