Wall clock time and date at job start Sat Apr 11 2020 06:13:24 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.52998 * 1 3 3 H 1.09003 * 109.54442 * 2 1 4 4 C 1.53097 * 109.50318 * 239.89221 * 2 1 3 5 5 O 1.43016 * 109.25478 * 183.21848 * 4 2 1 6 6 C 1.43015 * 113.73965 * 58.67152 * 5 4 2 7 7 C 1.53093 * 109.25810 * 301.32492 * 6 5 4 8 8 N 1.47011 * 109.50618 * 120.10963 * 2 1 3 9 9 C 1.46896 * 111.00160 * 302.32448 * 8 2 1 10 10 C 1.53002 * 109.47214 * 302.41102 * 9 8 2 11 11 N 1.46499 * 109.47041 * 179.97438 * 10 9 8 12 12 C 1.34777 * 120.00113 * 179.97438 * 11 10 9 13 13 O 1.21592 * 119.99835 * 4.97602 * 12 11 10 14 14 N 1.34774 * 120.00363 * 184.98105 * 12 11 10 15 15 C 1.47082 * 120.88098 * 354.77405 * 14 12 11 16 16 C 1.53240 * 108.68273 * 231.11641 * 15 14 12 17 17 N 1.46957 * 109.31055 * 306.46663 * 16 15 14 18 18 C 1.46893 * 110.96924 * 187.06611 * 17 16 15 19 19 C 1.50702 * 109.47268 * 65.49712 * 18 17 16 20 20 H 1.08001 * 119.99990 * 359.97438 * 19 18 17 21 21 C 1.37875 * 120.00164 * 179.97438 * 19 18 17 22 22 N 1.31514 * 119.99982 * 359.97438 * 21 19 18 23 23 Si 1.86795 * 120.00096 * 179.97438 * 21 19 18 24 24 H 1.48504 * 109.47278 * 90.00005 * 23 21 19 25 25 H 1.48503 * 109.47203 * 209.99923 * 23 21 19 26 26 H 1.48497 * 109.47323 * 329.99922 * 23 21 19 27 27 C 1.46904 * 111.44347 * 62.80995 * 17 16 15 28 28 C 1.47078 * 120.88449 * 174.80498 * 14 12 11 29 29 H 1.09009 * 109.47153 * 299.89104 * 1 2 3 30 30 H 1.08998 * 109.47301 * 59.88838 * 1 2 3 31 31 H 1.08999 * 109.47556 * 179.89055 * 1 2 3 32 32 H 1.09000 * 109.50679 * 63.28344 * 4 2 1 33 33 H 1.08997 * 109.50433 * 303.14903 * 4 2 1 34 34 H 1.09011 * 109.50515 * 181.39389 * 6 5 4 35 35 H 1.09002 * 109.54443 * 61.28828 * 6 5 4 36 36 H 1.08998 * 109.51079 * 296.91784 * 7 6 5 37 37 H 1.09000 * 109.54398 * 176.80486 * 7 6 5 38 38 H 1.08998 * 109.47292 * 182.41009 * 9 8 2 39 39 H 1.09001 * 109.47120 * 62.41181 * 9 8 2 40 40 H 1.08993 * 109.46873 * 60.00222 * 10 9 8 41 41 H 1.08994 * 109.46992 * 299.99694 * 10 9 8 42 42 H 0.97002 * 120.00495 * 0.02562 * 11 10 9 43 43 H 1.09001 * 109.62899 * 350.91585 * 15 14 12 44 44 H 1.08992 * 109.62950 * 111.30261 * 15 14 12 45 45 H 1.09000 * 109.49641 * 66.22525 * 16 15 14 46 46 H 1.09007 * 109.49012 * 186.31119 * 16 15 14 47 47 H 1.09001 * 109.47035 * 185.49268 * 18 17 16 48 48 H 1.09002 * 109.47364 * 305.49704 * 18 17 16 49 49 H 0.97003 * 119.99894 * 180.02562 * 22 21 19 50 50 H 1.08992 * 109.49508 * 177.02434 * 27 17 16 51 51 H 1.09008 * 109.48931 * 56.94244 * 27 17 16 52 52 H 1.09002 * 109.62371 * 248.50474 * 28 14 12 53 53 H 1.08999 * 109.62697 * 8.88051 * 28 14 12 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.8946 1.0272 0.0000 4 6 2.0411 -0.7239 -1.2484 5 8 3.4693 -0.7867 -1.2057 6 6 3.9817 -1.4579 -0.0514 7 6 3.4903 -0.7414 1.2091 8 7 2.0209 -0.6952 1.1988 9 6 1.5054 -0.0644 2.4212 10 6 1.9158 1.4093 2.4468 11 7 1.4013 2.0384 3.6656 12 6 1.6420 3.3429 3.9039 13 8 2.2073 4.0196 3.0667 14 7 1.2562 3.8974 5.0702 15 6 0.4716 3.1288 6.0484 16 6 -0.7591 3.9510 6.4453 17 7 -0.3295 5.2822 6.8958 18 6 -1.4611 6.0428 7.4425 19 6 -1.9596 5.3705 8.6958 20 1 -1.4836 4.4683 9.0505 21 6 -3.0232 5.9073 9.3896 22 7 -3.6025 7.0063 8.9580 23 14 -3.6415 5.0737 10.9427 24 1 -2.9190 5.6204 12.1193 25 1 -5.0970 5.3264 11.0943 26 1 -3.4002 3.6114 10.8493 27 6 0.3328 6.0267 5.8164 28 6 1.6081 5.2872 5.3983 29 1 -0.3634 0.5122 0.8910 30 1 -0.3634 0.5155 -0.8889 31 1 -0.3634 -1.0276 -0.0020 32 1 1.7272 -0.1794 -2.1390 33 1 1.6330 -1.7342 -1.2759 34 1 5.0715 -1.4443 -0.0764 35 1 3.6308 -2.4899 -0.0439 36 1 3.8862 0.2739 1.2284 37 1 3.8312 -1.2819 2.0921 38 1 1.9179 -0.5716 3.2934 39 1 0.4181 -0.1392 2.4379 40 1 1.5033 1.9164 1.5746 41 1 3.0031 1.4840 2.4301 42 1 0.8879 1.5160 4.3016 43 1 0.1538 2.1859 5.6032 44 1 1.0797 2.9304 6.9309 45 1 -1.4191 4.0569 5.5843 46 1 -1.2903 3.4447 7.2513 47 1 -1.1375 7.0566 7.6782 48 1 -2.2639 6.0795 6.7061 49 1 -4.3506 7.3841 9.4464 50 1 0.5903 7.0256 6.1682 51 1 -0.3397 6.1037 4.9619 52 1 2.3245 5.2993 6.2197 53 1 2.0451 5.7703 4.5243 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 ZINC000601583179.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:13:24 Heat of formation + Delta-G solvation = -45.470734 kcal Electronic energy + Delta-G solvation = -30594.931855 eV Core-core repulsion = 26502.290588 eV Total energy + Delta-G solvation = -4092.641267 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 0.77 seconds Orbital eigenvalues (eV) -40.60709 -40.03405 -37.79254 -36.89752 -34.62882 -34.34148 -33.02648 -32.42691 -30.89675 -29.91203 -28.30399 -27.62017 -26.25740 -25.21488 -23.50670 -23.23648 -22.95812 -22.08472 -20.80413 -20.04971 -19.28767 -17.91144 -17.02858 -16.91317 -16.44203 -16.20824 -15.98191 -15.70299 -15.23492 -15.16083 -14.83319 -14.58174 -14.43803 -14.22091 -13.95540 -13.89718 -13.71609 -13.58974 -13.10373 -12.74306 -12.66594 -12.47082 -12.30859 -12.21898 -12.04422 -11.90562 -11.62793 -11.31046 -11.22534 -11.07367 -11.02595 -10.93139 -10.78126 -10.68009 -10.54530 -10.04099 -10.02195 -9.88762 -9.31888 -9.21569 -8.82675 -8.51260 -8.36622 -7.34686 -6.09074 -4.11339 2.62182 3.24967 3.27342 3.27876 3.33744 3.35859 3.47056 3.91745 4.30043 4.35652 4.40546 4.47795 4.58294 4.69530 4.75670 4.80409 4.82405 4.96387 5.09637 5.16390 5.22080 5.22518 5.27167 5.29768 5.36563 5.38334 5.41355 5.44051 5.53255 5.57308 5.61503 5.63022 5.73253 5.82716 5.86903 6.04533 6.27594 6.29170 6.41265 6.43240 6.47449 6.50390 6.61018 6.72685 6.99666 7.19776 7.25804 7.70741 7.94884 7.99102 8.09709 8.44498 8.62625 8.95712 9.53210 11.01667 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.021226 B = 0.001967 C = 0.001902 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1318.833667 B =14228.161455 C =14715.665378 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.147 4.147 2 C 0.035 3.965 3 H 0.059 0.941 4 C 0.069 3.931 5 O -0.385 6.385 6 C 0.064 3.936 7 C 0.009 3.991 8 N -0.521 5.521 9 C 0.040 3.960 10 C 0.103 3.897 11 N -0.747 5.747 12 C 0.708 3.292 13 O -0.598 6.598 14 N -0.624 5.624 15 C 0.115 3.885 16 C 0.052 3.948 17 N -0.535 5.535 18 C 0.193 3.807 19 C -0.510 4.510 20 H 0.064 0.936 21 C 0.055 3.945 22 N -0.897 5.897 23 Si 0.899 3.101 24 H -0.280 1.280 25 H -0.288 1.288 26 H -0.274 1.274 27 C 0.041 3.959 28 C 0.136 3.864 29 H 0.063 0.937 30 H 0.064 0.936 31 H 0.066 0.934 32 H 0.098 0.902 33 H 0.060 0.940 34 H 0.098 0.902 35 H 0.059 0.941 36 H 0.050 0.950 37 H 0.089 0.911 38 H 0.084 0.916 39 H 0.089 0.911 40 H 0.067 0.933 41 H 0.063 0.937 42 H 0.397 0.603 43 H 0.068 0.932 44 H 0.075 0.925 45 H 0.024 0.976 46 H 0.097 0.903 47 H 0.032 0.968 48 H -0.009 1.009 49 H 0.262 0.738 50 H 0.082 0.918 51 H 0.026 0.974 52 H 0.070 0.930 53 H 0.088 0.912 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 11.525 -18.001 -16.581 27.052 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.073 4.073 3 H 0.077 0.923 4 C -0.008 4.008 5 O -0.304 6.304 6 C -0.013 4.013 7 C -0.119 4.119 8 N -0.246 5.246 9 C -0.088 4.088 10 C -0.023 4.023 11 N -0.406 5.406 12 C 0.409 3.591 13 O -0.477 6.477 14 N -0.361 5.361 15 C -0.008 4.008 16 C -0.078 4.078 17 N -0.261 5.261 18 C 0.068 3.932 19 C -0.548 4.548 20 H 0.082 0.918 21 C -0.146 4.146 22 N -0.536 5.536 23 Si 0.727 3.273 24 H -0.206 1.206 25 H -0.213 1.213 26 H -0.199 1.199 27 C -0.089 4.089 28 C 0.013 3.987 29 H 0.082 0.918 30 H 0.083 0.917 31 H 0.085 0.915 32 H 0.116 0.884 33 H 0.078 0.922 34 H 0.116 0.884 35 H 0.078 0.922 36 H 0.068 0.932 37 H 0.107 0.893 38 H 0.103 0.897 39 H 0.108 0.892 40 H 0.085 0.915 41 H 0.081 0.919 42 H 0.233 0.767 43 H 0.087 0.913 44 H 0.094 0.906 45 H 0.042 0.958 46 H 0.116 0.884 47 H 0.050 0.950 48 H 0.009 0.991 49 H 0.071 0.929 50 H 0.100 0.900 51 H 0.045 0.955 52 H 0.088 0.912 53 H 0.106 0.894 Dipole moment (debyes) X Y Z Total from point charges 12.093 -17.862 -16.954 27.436 hybrid contribution -1.247 0.334 0.927 1.589 sum 10.846 -17.528 -16.027 26.110 Atomic orbital electron populations 1.21831 0.93142 1.02875 1.02640 1.22438 0.96714 0.97202 0.90975 0.92279 1.22730 0.80812 1.00315 0.96915 1.88006 1.20753 1.78445 1.43153 1.22666 0.99071 0.95583 0.83987 1.22900 0.88014 1.01002 0.99941 1.62563 1.08557 1.51821 1.01619 1.22472 1.01899 0.91303 0.93172 1.21858 0.99307 0.93686 0.87484 1.45113 1.60518 1.09748 1.25221 1.15817 0.79461 0.81661 0.82160 1.90886 1.49406 1.61565 1.45853 1.47388 1.56930 1.14609 1.17201 1.21856 0.93526 0.93548 0.91874 1.22321 0.94034 0.92021 0.99400 1.62030 1.31648 1.02619 1.29775 1.19730 0.85577 0.93646 0.94200 1.21127 1.14356 1.10815 1.08540 0.91822 1.25005 0.96574 0.95109 0.97956 1.69959 1.23438 1.20417 1.39803 0.86485 0.77801 0.79713 0.83297 1.20605 1.21341 1.19870 1.22420 0.94998 0.97758 0.93696 1.21639 0.96989 0.79553 1.00563 0.91806 0.91681 0.91479 0.88390 0.92189 0.88360 0.92225 0.93152 0.89284 0.89726 0.89244 0.91515 0.91908 0.76684 0.91337 0.90642 0.95802 0.88444 0.94977 0.99128 0.92852 0.89998 0.95549 0.91156 0.89352 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 -0.98 8.89 37.16 0.33 -0.65 16 2 C 0.04 0.27 2.05 -67.60 -0.14 0.13 16 3 H 0.06 0.56 6.41 -51.93 -0.33 0.23 16 4 C 0.07 0.51 6.29 37.21 0.23 0.75 16 5 O -0.39 -3.55 10.83 -35.23 -0.38 -3.93 16 6 C 0.06 0.44 7.01 37.21 0.26 0.70 16 7 C 0.01 0.07 5.33 -3.71 -0.02 0.05 16 8 N -0.52 -3.92 4.75 -226.34 -1.07 -5.00 16 9 C 0.04 0.29 4.97 -3.83 -0.02 0.27 16 10 C 0.10 1.04 4.97 -4.04 -0.02 1.02 16 11 N -0.75 -8.53 5.46 -65.23 -0.36 -8.88 16 12 C 0.71 11.01 8.31 -86.92 -0.72 10.29 16 13 O -0.60 -11.42 15.80 5.28 0.08 -11.34 16 14 N -0.62 -9.51 3.01 -178.20 -0.54 -10.05 16 15 C 0.12 1.59 6.53 -3.60 -0.02 1.56 16 16 C 0.05 0.92 5.46 -3.67 -0.02 0.90 16 17 N -0.54 -10.99 4.62 -235.71 -1.09 -12.08 16 18 C 0.19 4.87 4.73 -4.98 -0.02 4.85 16 19 C -0.51 -13.71 8.41 -34.44 -0.29 -14.00 16 20 H 0.06 1.67 7.11 -52.49 -0.37 1.29 16 21 C 0.06 1.52 5.46 -17.82 -0.10 1.42 16 22 N -0.90 -26.10 13.29 55.43 0.74 -25.36 16 23 Si 0.90 20.95 29.54 -169.99 -5.02 15.93 16 24 H -0.28 -6.49 7.11 56.52 0.40 -6.09 16 25 H -0.29 -6.68 7.11 56.52 0.40 -6.28 16 26 H -0.27 -6.19 7.11 56.52 0.40 -5.79 16 27 C 0.04 0.73 5.62 -3.70 -0.02 0.71 16 28 C 0.14 2.15 6.61 -3.61 -0.02 2.13 16 29 H 0.06 0.46 6.43 -51.92 -0.33 0.13 16 30 H 0.06 0.38 8.14 -51.93 -0.42 -0.04 16 31 H 0.07 0.42 8.14 -51.93 -0.42 0.00 16 32 H 0.10 0.61 8.14 -51.93 -0.42 0.18 16 33 H 0.06 0.41 8.14 -51.93 -0.42 -0.01 16 34 H 0.10 0.53 8.14 -51.92 -0.42 0.11 16 35 H 0.06 0.38 8.14 -51.93 -0.42 -0.05 16 36 H 0.05 0.43 6.59 -51.93 -0.34 0.08 16 37 H 0.09 0.52 8.12 -51.93 -0.42 0.10 16 38 H 0.08 0.54 8.11 -51.93 -0.42 0.12 16 39 H 0.09 0.61 6.40 -51.93 -0.33 0.28 16 40 H 0.07 0.78 6.42 -51.93 -0.33 0.44 16 41 H 0.06 0.70 6.61 -51.93 -0.34 0.36 16 42 H 0.40 3.62 6.01 -40.82 -0.25 3.37 16 43 H 0.07 0.72 5.74 -51.93 -0.30 0.42 16 44 H 0.08 1.02 8.14 -51.93 -0.42 0.60 16 45 H 0.02 0.44 8.14 -51.93 -0.42 0.02 16 46 H 0.10 1.83 6.67 -51.92 -0.35 1.48 16 47 H 0.03 0.87 7.97 -51.93 -0.41 0.46 16 48 H -0.01 -0.25 8.13 -51.93 -0.42 -0.67 16 49 H 0.26 7.22 8.31 -40.82 -0.34 6.88 16 50 H 0.08 1.43 7.99 -51.93 -0.41 1.01 16 51 H 0.03 0.50 8.14 -51.92 -0.42 0.08 16 52 H 0.07 1.07 8.14 -51.93 -0.42 0.64 16 53 H 0.09 1.42 7.04 -51.93 -0.37 1.05 16 LS Contribution 400.72 15.07 6.04 6.04 Total: -1.00 -32.84 400.72 -11.29 -44.14 By element: Atomic # 1 Polarization: 9.51 SS G_CDS: -9.10 Total: 0.40 kcal Atomic # 6 Polarization: 10.72 SS G_CDS: -0.59 Total: 10.13 kcal Atomic # 7 Polarization: -59.05 SS G_CDS: -2.32 Total: -61.37 kcal Atomic # 8 Polarization: -14.97 SS G_CDS: -0.30 Total: -15.27 kcal Atomic # 14 Polarization: 20.95 SS G_CDS: -5.02 Total: 15.93 kcal Total LS contribution 6.04 Total: 6.04 kcal Total: -32.84 -11.29 -44.14 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. ZINC000601583179.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 -1.333 kcal (2) G-P(sol) polarization free energy of solvation -32.844 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -34.176 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.294 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.138 kcal (6) G-S(sol) free energy of system = (1) + (5) -45.471 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.77 seconds