Wall clock time and date at job start Sat Apr 11 2020 06:10:45 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.53001 * 1 3 3 H 1.09003 * 109.47008 * 2 1 4 4 C 1.50697 * 109.47048 * 240.00821 * 2 1 3 5 5 H 1.08007 * 120.00033 * 240.00064 * 4 2 1 6 6 C 1.37883 * 120.00382 * 59.99547 * 4 2 1 7 7 N 1.31512 * 119.99785 * 0.02562 * 6 4 2 8 8 Si 1.86799 * 120.00159 * 180.02562 * 6 4 2 9 9 H 1.48508 * 109.47012 * 90.00054 * 8 6 4 10 10 H 1.48494 * 109.47202 * 209.99919 * 8 6 4 11 11 H 1.48504 * 109.47158 * 330.00070 * 8 6 4 12 12 N 1.46491 * 109.47267 * 120.00606 * 2 1 3 13 13 C 1.34780 * 120.00254 * 274.99731 * 12 2 1 14 14 O 1.21284 * 119.99517 * 0.02562 * 13 12 2 15 15 C 1.50692 * 120.00439 * 180.02562 * 13 12 2 16 16 H 1.08997 * 110.88928 * 36.32440 * 15 13 12 17 17 C 1.54910 * 111.04109 * 160.22924 * 15 13 12 18 18 N 1.47417 * 104.83540 * 155.09637 * 17 15 13 19 19 C 1.34783 * 125.65056 * 156.13868 * 18 17 15 20 20 O 1.21508 * 120.00215 * 179.69017 * 19 18 17 21 21 O 1.34633 * 119.99491 * 359.69122 * 19 18 17 22 22 C 1.45200 * 116.99621 * 179.97438 * 21 19 18 23 23 C 1.53006 * 109.47087 * 60.00413 * 22 21 19 24 24 C 1.53002 * 109.47039 * 179.97438 * 22 21 19 25 25 C 1.52990 * 109.47491 * 300.00085 * 22 21 19 26 26 C 1.47024 * 108.70594 * 335.87420 * 18 17 15 27 27 C 1.54325 * 107.26914 * 1.05604 * 26 18 17 28 28 C 1.52995 * 110.72062 * 140.55429 * 27 26 18 29 29 C 1.53002 * 110.72284 * 263.83246 * 27 26 18 30 30 H 1.09003 * 109.47008 * 180.02562 * 1 2 3 31 31 H 1.09005 * 109.46912 * 299.99258 * 1 2 3 32 32 H 1.08996 * 109.47410 * 59.99054 * 1 2 3 33 33 H 0.96994 * 120.00152 * 179.97438 * 7 6 4 34 34 H 0.96999 * 120.00180 * 95.00165 * 12 2 1 35 35 H 1.09011 * 110.36206 * 36.25751 * 17 15 13 36 36 H 1.08998 * 110.37026 * 273.94082 * 17 15 13 37 37 H 1.08994 * 109.46996 * 59.99528 * 23 22 21 38 38 H 1.09002 * 109.47410 * 179.97438 * 23 22 21 39 39 H 1.09003 * 109.46748 * 300.00017 * 23 22 21 40 40 H 1.09002 * 109.46812 * 60.00605 * 24 22 21 41 41 H 1.09002 * 109.46920 * 180.02562 * 24 22 21 42 42 H 1.09002 * 109.47154 * 300.00333 * 24 22 21 43 43 H 1.09004 * 109.47129 * 60.00144 * 25 22 21 44 44 H 1.09002 * 109.47730 * 180.02562 * 25 22 21 45 45 H 1.09002 * 109.47382 * 300.00021 * 25 22 21 46 46 H 1.08994 * 109.88365 * 120.47022 * 26 18 17 47 47 H 1.08999 * 109.88557 * 241.48762 * 26 18 17 48 48 H 1.09003 * 109.47390 * 61.45963 * 28 27 26 49 49 H 1.09008 * 109.47296 * 181.45440 * 28 27 26 50 50 H 1.08999 * 109.47160 * 301.45696 * 28 27 26 51 51 H 1.08994 * 109.47103 * 64.54506 * 29 27 26 52 52 H 1.09003 * 109.46917 * 184.54666 * 29 27 26 53 53 H 1.08997 * 109.46922 * 304.54567 * 29 27 26 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.8933 1.0277 0.0000 4 6 2.0323 -0.7102 -1.2305 5 1 2.6533 -1.5884 -1.1316 6 6 1.6992 -0.2391 -2.4828 7 7 0.9434 0.8304 -2.6034 8 14 2.3213 -1.1199 -4.0081 9 1 1.3426 -2.1613 -4.4121 10 1 2.4864 -0.1409 -5.1124 11 1 3.6293 -1.7579 -3.7125 12 7 2.0183 -0.6907 1.1960 13 6 2.1472 -0.0185 2.3571 14 8 1.8583 1.1582 2.4117 15 6 2.6500 -0.7288 3.5873 16 1 3.4138 -1.4619 3.3279 17 6 3.1829 0.2813 4.6340 18 7 3.0568 -0.4140 5.9277 19 6 3.7771 -0.1486 7.0357 20 8 3.5838 -0.7846 8.0528 21 8 4.7107 0.8212 7.0152 22 6 5.4477 1.0525 8.2447 23 6 6.1905 -0.2248 8.6421 24 6 6.4567 2.1825 8.0295 25 6 4.4735 1.4432 9.3576 26 6 2.0109 -1.4407 5.8113 27 6 1.4817 -1.3934 4.3624 28 6 1.2172 -2.8055 3.8364 29 6 0.2175 -0.5360 4.2755 30 1 -0.3633 -1.0277 0.0005 31 1 -0.3633 0.5137 0.8901 32 1 -0.3634 0.5139 -0.8898 33 1 0.7094 1.1621 -3.4843 34 1 2.2491 -1.6318 1.1524 35 1 2.5767 1.1873 4.6296 36 1 4.2268 0.5213 4.4324 37 1 6.8846 -0.5031 7.8491 38 1 6.7442 -0.0510 9.5648 39 1 5.4720 -1.0301 8.7952 40 1 5.9279 3.0926 7.7462 41 1 7.0103 2.3558 8.9523 42 1 7.1508 1.9037 7.2366 43 1 3.7550 0.6379 9.5109 44 1 5.0268 1.6173 10.2805 45 1 3.9444 2.3532 9.0746 46 1 2.4319 -2.4241 6.0203 47 1 1.2026 -1.2274 6.5108 48 1 0.4556 -3.2855 4.4510 49 1 0.8693 -2.7497 2.8048 50 1 2.1379 -3.3873 3.8784 51 1 0.4613 0.4970 4.5233 52 1 -0.1837 -0.5819 3.2631 53 1 -0.5258 -0.9124 4.9783 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 ZINC000878533582.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:10:45 Heat of formation + Delta-G solvation = -129.011288 kcal Electronic energy + Delta-G solvation = -31634.358429 eV Core-core repulsion = 27505.801821 eV Total energy + Delta-G solvation = -4128.556608 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 2.46 seconds Orbital eigenvalues (eV) -41.07966 -40.05568 -38.04939 -36.68259 -34.82476 -34.47417 -33.46017 -32.36087 -30.33102 -28.70269 -27.60741 -27.57836 -27.33563 -26.78437 -25.29817 -24.16313 -22.09291 -21.00458 -19.83871 -19.49624 -19.03643 -17.79101 -17.18058 -16.80200 -16.36496 -15.95560 -15.39125 -15.17959 -14.77201 -14.70881 -14.37108 -14.14244 -14.05377 -13.77559 -13.41291 -13.21478 -13.07261 -12.69338 -12.64513 -12.47142 -12.41221 -12.26663 -12.09143 -12.06699 -11.97640 -11.70045 -11.63075 -11.55658 -11.46389 -11.15213 -11.06376 -10.97817 -10.87442 -10.73342 -10.57130 -10.50986 -10.34710 -10.26482 -10.14522 -9.10234 -9.06518 -8.76535 -8.58326 -7.95781 -6.15503 -4.24864 1.97738 2.75822 3.21427 3.23189 3.29357 3.31148 3.65569 4.06704 4.28647 4.36704 4.49114 4.51725 4.56101 4.69840 4.74683 5.05439 5.13899 5.17097 5.19622 5.25143 5.30322 5.32214 5.47896 5.48979 5.51538 5.56778 5.63025 5.64633 5.65357 5.69178 5.78111 5.87573 5.92208 5.99811 6.03144 6.05696 6.08013 6.15702 6.20141 6.35369 6.41371 6.78721 6.92578 7.11224 7.38029 7.41273 7.51806 7.65776 7.87969 8.11814 8.24682 8.61541 8.74210 9.06672 9.42889 10.90371 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.018282 B = 0.002512 C = 0.002383 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1531.225561 B =11142.324407 C =11745.561537 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.146 4.146 2 C 0.265 3.735 3 H 0.049 0.951 4 C -0.534 4.534 5 H 0.061 0.939 6 C 0.067 3.933 7 N -0.883 5.883 8 Si 0.901 3.099 9 H -0.279 1.279 10 H -0.285 1.285 11 H -0.274 1.274 12 N -0.704 5.704 13 C 0.519 3.481 14 O -0.550 6.550 15 C -0.121 4.121 16 H 0.102 0.898 17 C 0.122 3.878 18 N -0.602 5.602 19 C 0.669 3.331 20 O -0.572 6.572 21 O -0.359 6.359 22 C 0.135 3.865 23 C -0.183 4.183 24 C -0.141 4.141 25 C -0.183 4.183 26 C 0.138 3.862 27 C -0.059 4.059 28 C -0.141 4.141 29 C -0.142 4.142 30 H 0.026 0.974 31 H 0.015 0.985 32 H 0.089 0.911 33 H 0.264 0.736 34 H 0.388 0.612 35 H 0.091 0.909 36 H 0.082 0.918 37 H 0.059 0.941 38 H 0.069 0.931 39 H 0.087 0.913 40 H 0.067 0.933 41 H 0.079 0.921 42 H 0.067 0.933 43 H 0.087 0.913 44 H 0.069 0.931 45 H 0.060 0.940 46 H 0.074 0.926 47 H 0.074 0.926 48 H 0.063 0.937 49 H 0.071 0.929 50 H 0.061 0.939 51 H 0.066 0.934 52 H 0.068 0.932 53 H 0.058 0.942 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 6.003 -3.746 19.979 21.195 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.205 4.205 2 C 0.164 3.836 3 H 0.067 0.933 4 C -0.573 4.573 5 H 0.079 0.921 6 C -0.137 4.137 7 N -0.520 5.520 8 Si 0.728 3.272 9 H -0.206 1.206 10 H -0.211 1.211 11 H -0.199 1.199 12 N -0.352 5.352 13 C 0.302 3.698 14 O -0.429 6.429 15 C -0.143 4.143 16 H 0.120 0.880 17 C -0.001 4.001 18 N -0.342 5.342 19 C 0.425 3.575 20 O -0.463 6.463 21 O -0.274 6.274 22 C 0.098 3.902 23 C -0.240 4.240 24 C -0.198 4.198 25 C -0.240 4.240 26 C 0.016 3.984 27 C -0.060 4.060 28 C -0.198 4.198 29 C -0.199 4.199 30 H 0.045 0.955 31 H 0.034 0.966 32 H 0.108 0.892 33 H 0.074 0.926 34 H 0.221 0.779 35 H 0.109 0.891 36 H 0.101 0.899 37 H 0.078 0.922 38 H 0.088 0.912 39 H 0.106 0.894 40 H 0.087 0.913 41 H 0.098 0.902 42 H 0.086 0.914 43 H 0.106 0.894 44 H 0.088 0.912 45 H 0.079 0.921 46 H 0.092 0.908 47 H 0.092 0.908 48 H 0.082 0.918 49 H 0.090 0.910 50 H 0.080 0.920 51 H 0.085 0.915 52 H 0.087 0.913 53 H 0.077 0.923 Dipole moment (debyes) X Y Z Total from point charges 5.756 -2.898 20.267 21.267 hybrid contribution 0.481 -0.923 -0.909 1.382 sum 6.237 -3.821 19.358 20.693 Atomic orbital electron populations 1.22073 0.98447 0.99044 1.00898 1.18740 0.91000 0.93963 0.79934 0.93350 1.21212 1.32489 1.13225 0.90336 0.92114 1.24882 0.94811 0.94922 0.99055 1.69895 1.35863 1.22362 1.23913 0.86462 0.79173 0.79289 0.82257 1.20553 1.21078 1.19863 1.46535 1.64782 1.13475 1.10358 1.19767 0.78724 0.87220 0.84068 1.90665 1.51405 1.15926 1.84883 1.22645 0.98625 0.98598 0.94463 0.88015 1.22207 1.02502 0.95499 0.79905 1.47093 1.38345 1.38587 1.10221 1.17418 0.78448 0.79240 0.82397 1.90880 1.70933 1.53619 1.30824 1.86344 1.48305 1.52417 1.40297 1.22069 0.90697 0.96095 0.81311 1.22515 1.01314 0.98329 1.01830 1.21861 0.97855 0.97129 1.02977 1.22511 1.00524 1.03238 0.97733 1.22082 0.89135 0.90658 0.96485 1.21537 0.95034 0.96841 0.92625 1.21857 1.02131 0.93626 1.02203 1.21906 0.95219 0.99888 1.02928 0.95521 0.96633 0.89210 0.92559 0.77947 0.89057 0.89919 0.92170 0.91238 0.89419 0.91350 0.90213 0.91399 0.89425 0.91237 0.92123 0.90804 0.90801 0.91815 0.91041 0.92014 0.91518 0.91296 0.92311 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.15 -3.46 9.26 37.16 0.34 -3.12 16 2 C 0.27 6.42 2.56 -69.09 -0.18 6.24 16 3 H 0.05 1.37 7.58 -51.93 -0.39 0.98 16 4 C -0.53 -13.66 8.73 -34.44 -0.30 -13.96 16 5 H 0.06 1.48 7.74 -52.48 -0.41 1.08 16 6 C 0.07 1.75 5.30 -17.82 -0.09 1.65 16 7 N -0.88 -24.34 11.32 55.43 0.63 -23.71 16 8 Si 0.90 20.21 29.54 -169.99 -5.02 15.19 16 9 H -0.28 -6.20 7.11 56.52 0.40 -5.80 16 10 H -0.29 -6.43 7.11 56.52 0.40 -6.03 16 11 H -0.27 -6.05 7.11 56.52 0.40 -5.65 16 12 N -0.70 -13.53 5.07 -54.94 -0.28 -13.81 16 13 C 0.52 9.36 5.55 -10.99 -0.06 9.29 16 14 O -0.55 -12.00 14.92 5.55 0.08 -11.92 16 15 C -0.12 -1.63 2.75 -89.62 -0.25 -1.88 16 16 H 0.10 1.21 8.04 -51.93 -0.42 0.79 16 17 C 0.12 1.67 5.90 -2.53 -0.01 1.66 16 18 N -0.60 -7.97 3.33 -176.30 -0.59 -8.56 16 19 C 0.67 9.50 7.93 54.06 0.43 9.93 16 20 O -0.57 -8.78 12.21 -19.28 -0.24 -9.01 16 21 O -0.36 -4.76 9.42 -40.31 -0.38 -5.13 16 22 C 0.14 1.41 1.13 -90.62 -0.10 1.31 16 23 C -0.18 -1.79 8.37 37.16 0.31 -1.48 16 24 C -0.14 -1.13 8.85 37.16 0.33 -0.80 16 25 C -0.18 -1.82 8.37 37.15 0.31 -1.51 16 26 C 0.14 1.59 6.04 -3.06 -0.02 1.57 16 27 C -0.06 -0.69 0.76 -152.67 -0.12 -0.80 16 28 C -0.14 -1.38 8.49 37.16 0.32 -1.06 16 29 C -0.14 -1.84 7.41 37.16 0.28 -1.57 16 30 H 0.03 0.58 8.14 -51.93 -0.42 0.16 16 31 H 0.01 0.34 8.14 -51.93 -0.42 -0.08 16 32 H 0.09 2.35 6.07 -51.93 -0.32 2.04 16 33 H 0.26 6.94 8.31 -40.82 -0.34 6.60 16 34 H 0.39 6.47 8.54 -40.82 -0.35 6.12 16 35 H 0.09 1.41 6.71 -51.92 -0.35 1.06 16 36 H 0.08 1.11 8.04 -51.93 -0.42 0.69 16 37 H 0.06 0.57 8.14 -51.93 -0.42 0.15 16 38 H 0.07 0.53 8.14 -51.93 -0.42 0.11 16 39 H 0.09 1.05 5.88 -51.93 -0.31 0.75 16 40 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16 41 H 0.08 0.47 8.14 -51.93 -0.42 0.05 16 42 H 0.07 0.54 8.14 -51.93 -0.42 0.12 16 43 H 0.09 1.07 5.88 -51.93 -0.31 0.76 16 44 H 0.07 0.54 8.14 -51.93 -0.42 0.12 16 45 H 0.06 0.60 8.14 -51.93 -0.42 0.17 16 46 H 0.07 0.79 7.82 -51.93 -0.41 0.39 16 47 H 0.07 0.84 7.82 -51.93 -0.41 0.43 16 48 H 0.06 0.54 8.14 -51.93 -0.42 0.12 16 49 H 0.07 0.78 8.14 -51.92 -0.42 0.36 16 50 H 0.06 0.53 8.13 -51.93 -0.42 0.11 16 51 H 0.07 0.97 7.05 -51.93 -0.37 0.61 16 52 H 0.07 1.03 7.58 -51.93 -0.39 0.64 16 53 H 0.06 0.63 8.10 -51.93 -0.42 0.21 16 LS Contribution 413.36 15.07 6.23 6.23 Total: -1.00 -30.25 413.36 -7.83 -38.09 By element: Atomic # 1 Polarization: 16.62 SS G_CDS: -9.46 Total: 7.16 kcal Atomic # 6 Polarization: 4.29 SS G_CDS: 1.18 Total: 5.47 kcal Atomic # 7 Polarization: -45.84 SS G_CDS: -0.24 Total: -46.08 kcal Atomic # 8 Polarization: -25.53 SS G_CDS: -0.53 Total: -26.06 kcal Atomic # 14 Polarization: 20.21 SS G_CDS: -5.02 Total: 15.19 kcal Total LS contribution 6.23 Total: 6.23 kcal Total: -30.25 -7.83 -38.09 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. ZINC000878533582.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 -90.924 kcal (2) G-P(sol) polarization free energy of solvation -30.253 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -121.177 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.834 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -38.087 kcal (6) G-S(sol) free energy of system = (1) + (5) -129.011 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 2.46 seconds