Wall clock time and date at job start Sun Mar 7 2021 07:06:34 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.42905 * 1 3 3 C 1.35964 * 116.99979 * 2 1 4 4 C 1.39449 * 120.12454 * 180.02562 * 3 2 1 5 5 O 1.34776 * 133.23678 * 359.95041 * 4 3 2 6 6 C 1.34930 * 110.83703 * 180.02562 * 5 4 3 7 7 C 1.46997 * 125.13852 * 180.02562 * 6 5 4 8 8 O 1.21645 * 119.99894 * 179.97438 * 7 6 5 9 9 N 1.34777 * 119.99895 * 359.96908 * 7 6 5 10 10 C 1.46091 * 122.76464 * 0.02562 * 9 7 6 11 11 C 1.46013 * 114.87864 * 78.52434 * 10 9 7 12 12 C 1.56521 * 108.96926 * 60.18122 * 11 10 9 13 13 Si 1.83376 * 108.24137 * 291.77998 * 12 11 10 14 14 C 1.88359 * 111.35353 * 322.39847 * 13 12 11 15 15 C 1.56249 * 99.08913 * 221.45851 * 14 13 12 16 16 C 1.53713 * 109.68512 * 39.76719 * 15 14 13 17 17 C 1.55531 * 111.35697 * 311.79605 * 16 15 14 18 18 C 1.87777 * 113.96572 * 86.62842 * 13 12 11 19 19 C 1.44791 * 122.76574 * 179.70817 * 9 7 6 20 20 C 1.35809 * 109.72080 * 0.02562 * 6 5 4 21 21 C 1.40608 * 119.13650 * 179.97438 * 4 3 2 22 22 C 1.39773 * 120.12168 * 359.97438 * 21 4 3 23 23 C 1.50697 * 120.09087 * 179.97438 * 22 21 4 24 24 O 1.42900 * 109.47080 * 75.00442 * 23 22 21 25 25 Si 1.86310 * 109.47013 * 195.00183 * 23 22 21 26 26 H 1.48494 * 109.46823 * 300.00126 * 25 23 22 27 27 H 1.48492 * 109.46984 * 60.00082 * 25 23 22 28 28 H 1.48502 * 109.46945 * 180.02562 * 25 23 22 29 29 C 1.36625 * 119.81151 * 0.02562 * 22 21 4 30 30 C 1.38264 * 120.11839 * 359.72613 * 3 2 1 31 31 H 1.09001 * 109.47103 * 180.02562 * 1 2 3 32 32 H 1.09004 * 109.47054 * 299.99742 * 1 2 3 33 33 H 1.09000 * 109.47187 * 59.99797 * 1 2 3 34 34 H 1.09008 * 108.40298 * 199.92203 * 10 9 7 35 35 H 1.08997 * 108.40613 * 317.12260 * 10 9 7 36 36 H 1.09003 * 109.55353 * 180.02562 * 11 10 9 37 37 H 1.09004 * 109.55330 * 300.22933 * 11 10 9 38 38 H 1.08999 * 109.71284 * 51.46689 * 12 11 10 39 39 H 1.09003 * 109.70750 * 172.05138 * 12 11 10 40 40 H 1.08996 * 111.47567 * 338.92487 * 14 13 12 41 41 H 1.08996 * 111.47295 * 103.99027 * 14 13 12 42 42 H 1.09001 * 109.43816 * 279.69914 * 15 14 13 43 43 H 1.08997 * 109.39443 * 159.81528 * 15 14 13 44 44 H 1.08994 * 109.11579 * 191.29589 * 16 15 14 45 45 H 1.09001 * 109.02528 * 72.23796 * 16 15 14 46 46 H 1.08992 * 111.11135 * 270.25654 * 17 16 15 47 47 H 1.09000 * 111.15220 * 146.14078 * 17 16 15 48 48 H 1.09004 * 110.13287 * 80.07215 * 18 13 12 49 49 H 1.08999 * 110.13837 * 201.81460 * 18 13 12 50 50 H 1.09004 * 110.03754 * 50.61654 * 19 9 7 51 51 H 1.09000 * 110.03903 * 172.11224 * 19 9 7 52 52 H 1.07999 * 126.95167 * 180.02562 * 20 6 5 53 53 H 1.09003 * 109.47399 * 315.00246 * 23 22 21 54 54 H 0.96707 * 114.00115 * 179.97438 * 24 23 22 55 55 H 1.08006 * 119.74903 * 180.02562 * 29 22 21 56 56 H 1.08002 * 119.66509 * 0.02562 * 30 3 2 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.4291 0.0000 0.0000 3 6 2.0463 1.2114 0.0000 4 6 3.4387 1.2875 -0.0005 5 8 4.4141 0.3574 -0.0003 6 6 5.6317 0.9388 -0.0005 7 6 6.9132 0.2186 -0.0009 8 8 7.9595 0.8390 -0.0015 9 7 6.9287 -1.1291 -0.0003 10 6 5.7095 -1.9339 0.0001 11 6 5.0516 -2.0528 -1.2980 12 6 6.0389 -2.7256 -2.3091 13 14 7.4126 -1.5553 -2.6350 14 6 6.7956 0.2228 -2.7096 15 6 7.6365 0.7249 -3.9270 16 6 7.7023 -0.3742 -4.9995 17 6 8.0858 -1.7482 -4.3797 18 6 8.8231 -1.6722 -1.4010 19 6 8.1552 -1.8986 -0.0068 20 6 5.4944 2.2899 -0.0004 21 6 4.0554 2.5511 -0.0011 22 6 3.2765 3.7117 -0.0016 23 6 3.9381 5.0657 -0.0016 24 8 4.4609 5.3403 -1.3028 25 14 2.6861 6.3686 0.4523 26 1 1.5797 6.3583 -0.5381 27 1 2.1428 6.0832 1.8045 28 1 3.3380 7.7029 0.4517 29 6 1.9136 3.6151 -0.0011 30 6 1.2957 2.3726 -0.0057 31 1 -0.3633 -1.0277 0.0005 32 1 -0.3633 0.5138 0.8900 33 1 -0.3633 0.5139 -0.8900 34 1 5.9580 -2.9349 0.3529 35 1 5.0048 -1.4907 0.7036 36 1 4.1566 -2.6671 -1.1985 37 1 4.7758 -1.0621 -1.6594 38 1 6.4344 -3.6461 -1.8796 39 1 5.5168 -2.9489 -3.2395 40 1 5.7266 0.2658 -2.9183 41 1 7.0441 0.7698 -1.8003 42 1 8.6457 0.9686 -3.5950 43 1 7.1684 1.6148 -4.3478 44 1 8.4468 -0.0997 -5.7468 45 1 6.7283 -0.4595 -5.4813 46 1 9.1668 -1.8866 -4.3702 47 1 7.5920 -2.5672 -4.9027 48 1 9.4697 -2.5140 -1.6489 49 1 9.3967 -0.7453 -1.3978 50 1 8.8213 -1.5463 0.7808 51 1 7.9333 -2.9563 0.1354 52 1 6.2874 3.0231 -0.0009 53 1 4.7503 5.0733 0.7254 54 1 4.9016 6.1980 -1.3767 55 1 1.3127 4.5126 -0.0019 56 1 0.2174 2.3120 -0.0098 RHF calculation, no. of doubly occupied orbitals= 73 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_16608542_6672420.mol2 56 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 07:06:34 Heat of formation + Delta-G solvation = -73.875551 kcal Electronic energy + Delta-G solvation = -39425.929523 eV Core-core repulsion = 34812.746462 eV Total energy + Delta-G solvation = -4613.183061 eV No. of doubly occupied orbitals = 73 Molecular weight (most abundant/longest-lived isotopes) = 403.182 amu Computer time = 1.64 seconds Orbital eigenvalues (eV) -43.84674 -41.55389 -40.50351 -38.97395 -38.20967 -37.10107 -36.26152 -34.15968 -33.65134 -32.18960 -31.64826 -30.48807 -29.77068 -28.54207 -27.34634 -26.32181 -25.76610 -24.75503 -24.01426 -23.41235 -22.73105 -21.72880 -21.13799 -20.85298 -19.19300 -18.91894 -18.50687 -18.32087 -17.81556 -17.52659 -17.04426 -16.79455 -16.36095 -16.28515 -16.07448 -15.81132 -15.69490 -15.58419 -15.34695 -14.95818 -14.83758 -14.66198 -14.62358 -14.40201 -14.09493 -14.01115 -13.93571 -13.61000 -13.50480 -13.39306 -13.22497 -12.90408 -12.63348 -12.53760 -12.48084 -12.37858 -12.35278 -12.12418 -12.08461 -12.00583 -11.84891 -11.69510 -11.64063 -11.47734 -11.39195 -10.98677 -10.55604 -10.45329 -10.32878 -10.21753 -9.95487 -9.72834 -8.99218 -0.66291 0.07852 0.60489 1.11215 1.18535 1.46021 1.51747 1.63773 1.70280 1.82994 1.88275 1.94788 2.02703 2.12235 2.35699 2.65287 3.02479 3.36455 3.56259 3.71894 3.76405 3.83906 3.90216 4.00230 4.02605 4.05578 4.11880 4.13076 4.21032 4.28989 4.31540 4.36440 4.37058 4.40688 4.43281 4.51342 4.53424 4.60153 4.60716 4.74899 4.76635 4.82312 4.83549 4.91148 4.93192 4.96435 4.99846 5.01976 5.07250 5.09423 5.20507 5.23835 5.27854 5.29632 5.31836 5.37934 5.42156 5.45550 5.55599 5.64643 6.11706 6.24207 6.75125 6.77782 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.014 3.986 2 O -0.326 6.326 3 C 0.141 3.859 4 C 0.019 3.981 5 O -0.161 6.161 6 C -0.025 4.025 7 C 0.618 3.382 8 O -0.537 6.537 9 N -0.598 5.598 10 C 0.090 3.910 11 C -0.128 4.128 12 C -0.402 4.402 13 Si 1.021 2.979 14 C -0.443 4.443 15 C -0.118 4.118 16 C -0.118 4.118 17 C -0.414 4.414 18 C -0.425 4.425 19 C 0.112 3.888 20 C -0.121 4.121 21 C -0.095 4.095 22 C -0.071 4.071 23 C -0.091 4.091 24 O -0.562 6.562 25 Si 0.986 3.014 26 H -0.254 1.254 27 H -0.252 1.252 28 H -0.255 1.255 29 C -0.102 4.102 30 C -0.172 4.172 31 H 0.119 0.881 32 H 0.073 0.927 33 H 0.070 0.930 34 H 0.121 0.879 35 H 0.074 0.926 36 H 0.098 0.902 37 H 0.053 0.947 38 H 0.104 0.896 39 H 0.085 0.915 40 H 0.059 0.941 41 H 0.074 0.926 42 H 0.054 0.946 43 H 0.049 0.951 44 H 0.059 0.941 45 H 0.059 0.941 46 H 0.077 0.923 47 H 0.085 0.915 48 H 0.104 0.896 49 H 0.063 0.937 50 H 0.070 0.930 51 H 0.133 0.867 52 H 0.170 0.830 53 H 0.083 0.917 54 H 0.397 0.603 55 H 0.153 0.847 56 H 0.178 0.822 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -5.616 -3.217 2.020 6.780 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.080 4.080 2 O -0.240 6.240 3 C 0.095 3.905 4 C -0.029 4.029 5 O -0.062 6.062 6 C -0.075 4.075 7 C 0.407 3.593 8 O -0.416 6.416 9 N -0.329 5.329 10 C -0.031 4.031 11 C -0.166 4.166 12 C -0.491 4.491 13 Si 1.212 2.788 14 C -0.525 4.525 15 C -0.155 4.155 16 C -0.155 4.155 17 C -0.499 4.499 18 C -0.512 4.512 19 C -0.009 4.009 20 C -0.141 4.141 21 C -0.096 4.096 22 C -0.072 4.072 23 C -0.199 4.199 24 O -0.371 6.371 25 Si 0.808 3.192 26 H -0.178 1.178 27 H -0.176 1.176 28 H -0.179 1.179 29 C -0.119 4.119 30 C -0.190 4.190 31 H 0.138 0.862 32 H 0.092 0.908 33 H 0.089 0.911 34 H 0.139 0.861 35 H 0.092 0.908 36 H 0.116 0.884 37 H 0.072 0.928 38 H 0.123 0.877 39 H 0.104 0.896 40 H 0.077 0.923 41 H 0.092 0.908 42 H 0.073 0.927 43 H 0.068 0.932 44 H 0.078 0.922 45 H 0.078 0.922 46 H 0.096 0.904 47 H 0.104 0.896 48 H 0.122 0.878 49 H 0.082 0.918 50 H 0.089 0.911 51 H 0.151 0.849 52 H 0.187 0.813 53 H 0.101 0.899 54 H 0.247 0.753 55 H 0.170 0.830 56 H 0.195 0.805 Dipole moment (debyes) X Y Z Total from point charges -5.887 -3.359 2.180 7.120 hybrid contribution 0.611 1.177 -0.415 1.390 sum -5.276 -2.182 1.765 5.976 Atomic orbital electron populations 1.23480 0.75521 1.06100 1.02905 1.86139 1.24065 1.25362 1.88395 1.17861 0.86182 0.86675 0.99781 1.19196 0.88417 0.87365 1.07917 1.84288 1.12276 1.33674 1.76009 1.21965 0.84435 0.92556 1.08572 1.16282 0.85174 0.82144 0.75664 1.90902 1.31322 1.67576 1.51844 1.47648 1.10285 1.04198 1.70818 1.20599 0.85981 0.97389 0.99156 1.20955 0.99378 1.01482 0.94749 1.26310 1.09910 1.09648 1.03206 0.79473 0.67459 0.65617 0.66285 1.27743 1.01992 1.22465 1.00269 1.21918 0.99459 0.98559 0.95605 1.21793 1.01422 0.94306 0.98022 1.26650 1.05872 0.97982 1.19371 1.27359 1.12468 1.05164 1.06178 1.21940 0.84279 0.97158 0.97552 1.22762 0.95119 0.94557 1.01706 1.19544 0.92533 0.92793 1.04714 1.18709 0.92999 0.92800 1.02667 1.24948 1.05813 1.00055 0.89115 1.86496 1.78260 1.43274 1.29042 0.85201 0.77722 0.76276 0.80031 1.17831 1.17568 1.17883 1.20801 0.93219 0.98849 0.99077 1.21305 1.01751 0.89692 1.06279 0.86235 0.90839 0.91099 0.86050 0.90778 0.88400 0.92806 0.87731 0.89633 0.92256 0.90788 0.92713 0.93174 0.92217 0.92189 0.90397 0.89582 0.87756 0.91763 0.91142 0.84937 0.81255 0.89874 0.75335 0.83012 0.80521 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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.01 -0.08 9.81 113.37 1.11 1.03 16 2 O -0.33 -0.93 10.50 -93.14 -0.98 -1.91 16 3 C 0.14 0.70 6.70 22.57 0.15 0.85 16 4 C 0.02 0.17 7.35 23.14 0.17 0.34 16 5 O -0.16 -1.33 5.81 -13.43 -0.08 -1.40 16 6 C -0.02 -0.24 5.32 24.00 0.13 -0.11 16 7 C 0.62 5.56 4.11 86.54 0.36 5.92 16 8 O -0.54 -7.79 15.20 -4.17 -0.06 -7.85 16 9 N -0.60 -1.39 1.63 -837.47 -1.36 -2.75 16 10 C 0.09 -0.17 4.41 84.31 0.37 0.20 16 11 C -0.13 0.15 5.25 29.68 0.16 0.30 16 12 C -0.40 1.34 4.83 73.21 0.35 1.69 16 13 Si 1.02 0.54 8.77 68.60 0.60 1.14 16 14 C -0.44 -3.62 4.43 73.11 0.32 -3.29 16 15 C -0.12 -1.01 5.96 31.96 0.19 -0.82 16 16 C -0.12 -0.65 6.09 31.71 0.19 -0.45 16 17 C -0.41 -0.83 5.84 72.85 0.43 -0.41 16 18 C -0.43 -0.30 5.03 73.11 0.37 0.07 16 19 C 0.11 -0.13 5.78 87.83 0.51 0.38 16 20 C -0.12 -1.21 10.68 23.79 0.25 -0.96 16 21 C -0.09 -0.87 6.16 -20.25 -0.12 -0.99 16 22 C -0.07 -0.41 5.30 -19.87 -0.11 -0.52 16 23 C -0.09 -0.15 2.70 71.24 0.19 0.04 16 24 O -0.56 -1.79 12.86 -148.98 -1.92 -3.70 16 25 Si 0.99 -1.57 25.50 68.60 1.75 0.18 16 26 H -0.25 0.12 6.74 99.48 0.67 0.79 16 27 H -0.25 0.22 7.11 99.48 0.71 0.93 16 28 H -0.25 0.39 7.11 99.48 0.71 1.10 16 29 C -0.10 -0.24 8.62 22.02 0.19 -0.05 16 30 C -0.17 -0.13 9.03 22.40 0.20 0.07 16 31 H 0.12 -1.09 8.14 -2.39 -0.02 -1.11 16 32 H 0.07 -0.64 7.67 -2.38 -0.02 -0.65 16 33 H 0.07 -0.57 7.65 -2.39 -0.02 -0.59 16 34 H 0.12 -1.02 7.19 -2.38 -0.02 -1.04 16 35 H 0.07 0.03 6.19 -2.39 -0.01 0.02 16 36 H 0.10 -0.32 8.14 -2.39 -0.02 -0.34 16 37 H 0.05 0.21 4.69 -2.38 -0.01 0.20 16 38 H 0.10 -0.81 7.93 -2.39 -0.02 -0.83 16 39 H 0.08 -0.28 7.95 -2.39 -0.02 -0.30 16 40 H 0.06 0.58 7.12 -2.39 -0.02 0.57 16 41 H 0.07 0.84 3.11 -2.39 -0.01 0.83 16 42 H 0.05 0.52 8.14 -2.39 -0.02 0.50 16 43 H 0.05 0.50 8.14 -2.39 -0.02 0.48 16 44 H 0.06 0.32 8.14 -2.39 -0.02 0.30 16 45 H 0.06 0.30 8.14 -2.39 -0.02 0.28 16 46 H 0.08 0.13 8.03 -2.39 -0.02 0.11 16 47 H 0.09 0.00 8.03 -2.39 -0.02 -0.02 16 48 H 0.10 -0.37 8.00 -2.38 -0.02 -0.39 16 49 H 0.06 0.42 7.50 -2.39 -0.02 0.40 16 50 H 0.07 0.07 8.14 -2.38 -0.02 0.06 16 51 H 0.13 -1.18 7.19 -2.39 -0.02 -1.19 16 52 H 0.17 1.51 8.06 -2.91 -0.02 1.48 16 53 H 0.08 -0.05 7.96 -2.39 -0.02 -0.07 16 54 H 0.40 -2.36 8.74 -74.05 -0.65 -3.00 16 55 H 0.15 -0.30 4.84 -2.91 -0.01 -0.32 16 56 H 0.18 -1.13 6.26 -2.91 -0.02 -1.15 16 Total: 0.00 -20.31 415.70 4.36 -15.96 By element: Atomic # 1 Polarization: -3.96 SS G_CDS: 0.99 Total: -2.97 kcal Atomic # 6 Polarization: -2.11 SS G_CDS: 5.41 Total: 3.31 kcal Atomic # 7 Polarization: -1.39 SS G_CDS: -1.36 Total: -2.75 kcal Atomic # 8 Polarization: -11.83 SS G_CDS: -3.04 Total: -14.86 kcal Atomic # 14 Polarization: -1.03 SS G_CDS: 2.35 Total: 1.32 kcal Total: -20.31 4.36 -15.96 kcal The number of atoms in the molecule is 56 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_16608542_6672420.mol2 56 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 -57.920 kcal (2) G-P(sol) polarization free energy of solvation -20.314 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -78.234 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 4.358 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -15.956 kcal (6) G-S(sol) free energy of system = (1) + (5) -73.876 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.65 seconds