Wall clock time and date at job start Sat Apr 11 2020 06:12:09 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 O 1.45206 * 1 3 3 C 1.34224 * 117.00077 * 2 1 4 4 O 1.20833 * 120.00103 * 359.97438 * 3 2 1 5 5 C 1.50696 * 120.00302 * 180.02562 * 3 2 1 6 6 H 1.09009 * 109.46972 * 299.99362 * 5 3 2 7 7 N 1.46904 * 109.47493 * 180.02562 * 5 3 2 8 8 C 1.46895 * 111.00065 * 84.99579 * 7 5 3 9 9 C 1.53004 * 109.47348 * 180.02562 * 8 7 5 10 10 C 1.50698 * 109.47386 * 180.02562 * 9 8 7 11 11 O 1.21280 * 120.00098 * 359.97438 * 10 9 8 12 12 N 1.34779 * 119.99971 * 179.97438 * 10 9 8 13 13 C 1.46497 * 119.99661 * 0.02562 * 12 10 9 14 14 C 1.46498 * 119.99913 * 179.97438 * 12 10 9 15 15 C 1.50706 * 109.46767 * 90.00404 * 14 12 10 16 16 H 1.08000 * 120.00080 * 0.02562 * 15 14 12 17 17 C 1.37876 * 119.99833 * 179.97438 * 15 14 12 18 18 N 1.31513 * 119.99888 * 0.02562 * 17 15 14 19 19 Si 1.86798 * 120.00007 * 179.97438 * 17 15 14 20 20 H 1.48504 * 109.47003 * 90.00004 * 19 17 15 21 21 H 1.48497 * 109.47370 * 210.00240 * 19 17 15 22 22 H 1.48496 * 109.47019 * 330.00459 * 19 17 15 23 23 C 1.53000 * 109.47495 * 59.99416 * 5 3 2 24 24 C 1.53006 * 109.46588 * 60.00415 * 23 5 3 25 25 C 1.53000 * 109.46588 * 179.97438 * 24 23 5 26 26 C 1.52999 * 109.47002 * 60.00305 * 25 24 23 27 27 C 1.52998 * 109.47310 * 299.99119 * 26 25 24 28 28 C 1.52999 * 109.47002 * 300.00720 * 23 5 3 29 29 H 1.09000 * 109.47070 * 300.00162 * 1 2 3 30 30 H 1.09001 * 109.46626 * 60.00180 * 1 2 3 31 31 H 1.09001 * 109.47153 * 179.97438 * 1 2 3 32 32 H 1.00893 * 110.99542 * 321.04266 * 7 5 3 33 33 H 1.09000 * 109.47443 * 300.00296 * 8 7 5 34 34 H 1.09000 * 109.47459 * 60.00245 * 8 7 5 35 35 H 1.09000 * 109.46728 * 300.00290 * 9 8 7 36 36 H 1.09000 * 109.47004 * 59.99522 * 9 8 7 37 37 H 1.08999 * 109.47839 * 89.99873 * 13 12 10 38 38 H 1.09002 * 109.46572 * 209.99739 * 13 12 10 39 39 H 1.09003 * 109.47438 * 329.99533 * 13 12 10 40 40 H 1.08995 * 109.47056 * 210.00073 * 14 12 10 41 41 H 1.08995 * 109.47010 * 330.00543 * 14 12 10 42 42 H 0.96998 * 119.99879 * 179.97438 * 18 17 15 43 43 H 1.08994 * 109.47081 * 180.02562 * 23 5 3 44 44 H 1.09004 * 109.46892 * 300.00052 * 24 23 5 45 45 H 1.08998 * 109.47144 * 60.00447 * 24 23 5 46 46 H 1.08994 * 109.47539 * 299.99631 * 25 24 23 47 47 H 1.09001 * 109.46747 * 179.97438 * 25 24 23 48 48 H 1.09008 * 109.47133 * 179.97438 * 26 25 24 49 49 H 1.08994 * 109.47529 * 59.99875 * 26 25 24 50 50 H 1.08994 * 109.47441 * 300.00810 * 27 26 25 51 51 H 1.08998 * 109.47254 * 179.97438 * 27 26 25 52 52 H 1.09001 * 109.46780 * 299.99257 * 28 23 5 53 53 H 1.08994 * 109.47529 * 59.99875 * 28 23 5 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.4521 0.0000 0.0000 3 6 2.0614 1.1959 0.0000 4 8 1.4034 2.2093 -0.0005 5 6 3.5663 1.2749 -0.0006 6 1 3.9564 0.7806 0.8893 7 7 3.9829 2.6836 0.0002 8 6 4.0180 3.2215 1.3666 9 6 4.4545 4.6874 1.3267 10 6 4.4912 5.2390 2.7286 11 8 4.1967 4.5297 3.6673 12 7 4.8512 6.5206 2.9393 13 6 5.2063 7.3775 1.8054 14 6 4.8863 7.0570 4.3021 15 6 6.2557 6.8394 4.8924 16 1 7.0232 6.3509 4.3105 17 6 6.5291 7.2650 6.1751 18 7 5.5946 7.8604 6.8835 19 14 8.2263 6.9947 6.9070 20 1 9.0892 8.1640 6.6012 21 1 8.1111 6.8361 8.3790 22 1 8.8279 5.7693 6.3227 23 6 4.1130 0.5815 -1.2501 24 6 3.5655 1.2753 -2.4990 25 6 4.1126 0.5823 -3.7485 26 6 3.6793 -0.8850 -3.7489 27 6 4.2269 -1.5789 -2.5001 28 6 3.6791 -0.8857 -1.2501 29 1 -0.3633 0.5139 0.8900 30 1 -0.3632 0.5138 -0.8900 31 1 -0.3633 -1.0277 -0.0005 32 1 3.3827 3.2383 -0.5915 33 1 3.0252 3.1504 1.8108 34 1 4.7256 2.6475 1.9650 35 1 5.4473 4.7588 0.8824 36 1 3.7466 5.2611 0.7285 37 1 4.3126 7.8804 1.4360 38 1 5.9353 8.1215 2.1267 39 1 5.6355 6.7682 1.0099 40 1 4.6663 8.1243 4.2783 41 1 4.1420 6.5461 4.9129 42 1 5.7868 8.1595 7.7860 43 1 5.2015 0.6382 -1.2502 44 1 2.4769 1.2182 -2.4986 45 1 3.8746 2.3205 -2.4990 46 1 5.2011 0.6398 -3.7487 47 1 3.7222 1.0763 -4.6383 48 1 4.0688 -1.3787 -4.6394 49 1 2.5909 -0.9425 -3.7488 50 1 5.3153 -1.5214 -2.5002 51 1 3.9178 -2.6241 -2.5001 52 1 4.0692 -1.3800 -0.3604 53 1 2.5907 -0.9428 -1.2498 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 ZINC000542619548.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:12:09 Heat of formation + Delta-G solvation = -142.333442 kcal Electronic energy + Delta-G solvation = -30750.517724 eV Core-core repulsion = 26621.383424 eV Total energy + Delta-G solvation = -4129.134300 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.63 seconds Orbital eigenvalues (eV) -41.20094 -39.02146 -38.70630 -37.15765 -35.74413 -34.74746 -32.82758 -32.42436 -31.02515 -29.80085 -28.33832 -27.64271 -26.00542 -25.44023 -23.53867 -22.77872 -22.61361 -21.63469 -20.82714 -19.90343 -18.35511 -17.83112 -17.18632 -16.80072 -16.53303 -16.27532 -15.82740 -15.34030 -15.16272 -14.76736 -14.60701 -14.38793 -14.21168 -14.16078 -13.95204 -13.50303 -13.43516 -13.26082 -12.66215 -12.46192 -12.42100 -12.24637 -12.16030 -12.12523 -11.84108 -11.75121 -11.58260 -11.50496 -11.40255 -11.17487 -11.09477 -11.01406 -10.85650 -10.66530 -10.57130 -10.53549 -10.36028 -10.26843 -10.13529 -9.60035 -8.86719 -8.79938 -8.41559 -7.60383 -6.02023 -4.06711 1.95475 2.46969 3.31498 3.37803 3.41516 3.49806 3.64423 4.17686 4.34018 4.38796 4.47124 4.51121 4.57568 4.83409 4.83872 4.91509 4.95269 5.03055 5.08586 5.15531 5.24173 5.24573 5.25935 5.30126 5.32053 5.33716 5.43492 5.45027 5.61712 5.65405 5.78208 5.82981 5.84361 5.87598 5.94421 5.96365 5.99353 6.04731 6.07856 6.08889 6.13243 6.26041 6.30755 6.68652 7.06739 7.11461 7.47355 7.48877 7.57750 7.94551 8.26150 8.45870 8.92912 9.13391 9.58265 11.06082 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.016806 B = 0.002189 C = 0.002128 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1665.656158 B =12788.169919 C =13151.986620 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.351 6.351 3 C 0.436 3.564 4 O -0.508 6.508 5 C 0.099 3.901 6 H 0.118 0.882 7 N -0.658 5.658 8 C 0.068 3.932 9 C -0.135 4.135 10 C 0.506 3.494 11 O -0.545 6.545 12 N -0.603 5.603 13 C 0.084 3.916 14 C 0.244 3.756 15 C -0.524 4.524 16 H 0.057 0.943 17 C 0.059 3.941 18 N -0.896 5.896 19 Si 0.899 3.101 20 H -0.280 1.280 21 H -0.285 1.285 22 H -0.274 1.274 23 C -0.068 4.068 24 C -0.117 4.117 25 C -0.116 4.116 26 C -0.121 4.121 27 C -0.116 4.116 28 C -0.129 4.129 29 H 0.062 0.938 30 H 0.061 0.939 31 H 0.099 0.901 32 H 0.351 0.649 33 H 0.044 0.956 34 H 0.084 0.916 35 H 0.096 0.904 36 H 0.084 0.916 37 H 0.043 0.957 38 H 0.073 0.927 39 H 0.059 0.941 40 H 0.031 0.969 41 H 0.043 0.957 42 H 0.263 0.737 43 H 0.096 0.904 44 H 0.052 0.948 45 H 0.073 0.927 46 H 0.065 0.935 47 H 0.060 0.940 48 H 0.059 0.941 49 H 0.058 0.942 50 H 0.062 0.938 51 H 0.060 0.940 52 H 0.065 0.935 53 H 0.084 0.916 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -6.184 -17.256 -23.418 29.739 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.065 4.065 2 O -0.266 6.266 3 C 0.277 3.723 4 O -0.395 6.395 5 C -0.013 4.013 6 H 0.136 0.864 7 N -0.289 5.289 8 C -0.060 4.060 9 C -0.176 4.176 10 C 0.292 3.708 11 O -0.424 6.424 12 N -0.335 5.335 13 C -0.061 4.061 14 C 0.125 3.875 15 C -0.563 4.563 16 H 0.075 0.925 17 C -0.142 4.142 18 N -0.535 5.535 19 Si 0.726 3.274 20 H -0.206 1.206 21 H -0.211 1.211 22 H -0.199 1.199 23 C -0.087 4.087 24 C -0.155 4.155 25 C -0.154 4.154 26 C -0.158 4.158 27 C -0.154 4.154 28 C -0.167 4.167 29 H 0.081 0.919 30 H 0.079 0.921 31 H 0.118 0.882 32 H 0.172 0.828 33 H 0.063 0.937 34 H 0.103 0.897 35 H 0.115 0.885 36 H 0.103 0.897 37 H 0.062 0.938 38 H 0.092 0.908 39 H 0.077 0.923 40 H 0.049 0.951 41 H 0.061 0.939 42 H 0.073 0.927 43 H 0.114 0.886 44 H 0.071 0.929 45 H 0.091 0.909 46 H 0.083 0.917 47 H 0.078 0.922 48 H 0.077 0.923 49 H 0.077 0.923 50 H 0.080 0.920 51 H 0.079 0.921 52 H 0.084 0.916 53 H 0.103 0.897 Dipole moment (debyes) X Y Z Total from point charges -5.941 -17.348 -22.730 29.204 hybrid contribution -0.294 0.290 0.696 0.809 sum -6.235 -17.058 -22.033 28.554 Atomic orbital electron populations 1.23279 0.77169 1.03551 1.02537 1.86640 1.22924 1.33436 1.83632 1.25320 0.93645 0.81615 0.71702 1.90684 1.67028 1.36391 1.45385 1.21708 0.84713 0.92694 1.02234 0.86418 1.60927 1.55002 1.00687 1.12320 1.21815 1.01373 0.94800 0.88050 1.21567 1.03750 0.95966 0.96268 1.20236 0.77996 0.82334 0.90258 1.90663 1.49419 1.57652 1.44632 1.48150 1.65373 1.10436 1.09581 1.21641 0.99884 0.95393 0.89152 1.19283 0.97929 0.94049 0.76232 1.21437 0.97187 1.38433 0.99283 0.92460 1.24951 0.97401 0.95121 0.96765 1.69937 1.32936 1.39906 1.10685 0.86486 0.84046 0.78510 0.78321 1.20650 1.21107 1.19871 1.21134 1.01716 0.94173 0.91663 1.21616 0.99700 1.00090 0.94078 1.21494 1.01058 0.95385 0.97461 1.21552 1.00711 0.95747 0.97811 1.21504 1.00538 0.98902 0.94414 1.21892 1.02781 0.94010 0.97982 0.91895 0.92071 0.88210 0.82814 0.93716 0.89748 0.88544 0.89746 0.93834 0.90837 0.92262 0.95061 0.93902 0.92748 0.88626 0.92921 0.90856 0.91663 0.92168 0.92267 0.92277 0.91972 0.92084 0.91603 0.89741 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.22 10.56 101.05 1.07 1.29 16 2 O -0.35 -3.45 7.48 -39.24 -0.29 -3.74 16 3 C 0.44 5.02 5.76 36.00 0.21 5.23 16 4 O -0.51 -7.27 15.06 -18.76 -0.28 -7.56 16 5 C 0.10 1.00 2.37 -68.86 -0.16 0.84 16 6 H 0.12 1.17 7.80 -51.92 -0.40 0.76 16 7 N -0.66 -7.44 5.62 -107.95 -0.61 -8.05 16 8 C 0.07 0.99 4.67 -3.82 -0.02 0.98 16 9 C -0.13 -2.03 4.51 -27.88 -0.13 -2.16 16 10 C 0.51 10.74 7.67 -10.99 -0.08 10.66 16 11 O -0.55 -14.01 15.56 5.56 0.09 -13.92 16 12 N -0.60 -13.15 2.86 -182.96 -0.52 -13.67 16 13 C 0.08 1.35 10.24 59.85 0.61 1.96 16 14 C 0.24 6.57 5.16 -5.19 -0.03 6.54 16 15 C -0.52 -15.15 9.39 -34.44 -0.32 -15.47 16 16 H 0.06 1.65 7.81 -52.49 -0.41 1.24 16 17 C 0.06 1.73 5.46 -17.82 -0.10 1.63 16 18 N -0.90 -27.05 13.29 55.43 0.74 -26.31 16 19 Si 0.90 22.13 29.54 -169.99 -5.02 17.11 16 20 H -0.28 -6.68 7.11 56.52 0.40 -6.28 16 21 H -0.29 -6.96 7.11 56.52 0.40 -6.56 16 22 H -0.27 -6.72 7.11 56.52 0.40 -6.31 16 23 C -0.07 -0.55 2.21 -90.62 -0.20 -0.75 16 24 C -0.12 -0.93 4.69 -26.73 -0.13 -1.06 16 25 C -0.12 -0.78 5.92 -26.73 -0.16 -0.94 16 26 C -0.12 -0.76 5.92 -26.73 -0.16 -0.92 16 27 C -0.12 -0.73 5.93 -26.70 -0.16 -0.89 16 28 C -0.13 -0.94 4.69 -26.70 -0.13 -1.06 16 29 H 0.06 0.54 8.13 -51.93 -0.42 0.12 16 30 H 0.06 0.47 8.13 -51.93 -0.42 0.04 16 31 H 0.10 0.49 8.14 -51.93 -0.42 0.07 16 32 H 0.35 3.98 7.00 -39.30 -0.28 3.71 16 33 H 0.04 0.75 7.58 -51.93 -0.39 0.36 16 34 H 0.08 1.30 8.01 -51.93 -0.42 0.88 16 35 H 0.10 1.30 7.40 -51.93 -0.38 0.92 16 36 H 0.08 1.12 8.11 -51.93 -0.42 0.70 16 37 H 0.04 0.61 8.14 -51.93 -0.42 0.19 16 38 H 0.07 1.22 8.07 -51.93 -0.42 0.80 16 39 H 0.06 0.78 6.44 -51.93 -0.33 0.44 16 40 H 0.03 0.85 8.07 -51.93 -0.42 0.43 16 41 H 0.04 1.30 7.42 -51.93 -0.39 0.91 16 42 H 0.26 7.48 8.31 -40.82 -0.34 7.14 16 43 H 0.10 0.77 8.14 -51.93 -0.42 0.35 16 44 H 0.05 0.47 7.22 -51.93 -0.37 0.10 16 45 H 0.07 0.61 7.21 -51.93 -0.37 0.24 16 46 H 0.06 0.43 8.14 -51.93 -0.42 0.01 16 47 H 0.06 0.38 8.14 -51.93 -0.42 -0.04 16 48 H 0.06 0.33 8.14 -51.92 -0.42 -0.09 16 49 H 0.06 0.39 8.14 -51.93 -0.42 -0.03 16 50 H 0.06 0.39 8.14 -51.93 -0.42 -0.04 16 51 H 0.06 0.35 8.14 -51.93 -0.42 -0.07 16 52 H 0.07 0.47 8.14 -51.93 -0.42 0.04 16 53 H 0.08 0.68 5.32 -51.93 -0.28 0.40 16 LS Contribution 415.32 15.07 6.26 6.26 Total: -1.00 -34.57 415.32 -9.02 -43.59 By element: Atomic # 1 Polarization: 9.92 SS G_CDS: -9.50 Total: 0.42 kcal Atomic # 6 Polarization: 5.75 SS G_CDS: 0.12 Total: 5.87 kcal Atomic # 7 Polarization: -47.63 SS G_CDS: -0.39 Total: -48.03 kcal Atomic # 8 Polarization: -24.73 SS G_CDS: -0.49 Total: -25.22 kcal Atomic # 14 Polarization: 22.13 SS G_CDS: -5.02 Total: 17.11 kcal Total LS contribution 6.26 Total: 6.26 kcal Total: -34.57 -9.02 -43.59 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. ZINC000542619548.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 -98.744 kcal (2) G-P(sol) polarization free energy of solvation -34.571 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -133.315 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.018 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.589 kcal (6) G-S(sol) free energy of system = (1) + (5) -142.333 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.63 seconds