Wall clock time and date at job start Sat Apr 11 2020 06:11:06 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.30997 * 1 3 3 C 1.50698 * 120.00248 * 2 1 4 4 C 1.52999 * 109.47589 * 119.99851 * 3 2 1 5 5 C 1.53006 * 109.46758 * 240.00058 * 3 2 1 6 6 C 1.53005 * 109.47064 * 299.99461 * 5 3 2 7 7 O 1.42901 * 109.47106 * 180.02562 * 6 5 3 8 8 C 1.42898 * 113.99883 * 180.02562 * 7 6 5 9 9 C 1.50693 * 109.47077 * 359.97438 * 3 2 1 10 10 O 1.21284 * 120.00376 * 239.99898 * 9 3 2 11 11 N 1.34780 * 120.00056 * 60.00560 * 9 3 2 12 12 C 1.46498 * 120.00013 * 179.97438 * 11 9 3 13 13 C 1.53006 * 109.46791 * 180.02562 * 12 11 9 14 14 H 1.08999 * 108.45350 * 57.87439 * 13 12 11 15 15 C 1.46919 * 108.61886 * 175.75769 * 13 12 11 16 16 N 1.49896 * 106.62338 * 177.93437 * 15 13 12 17 17 C 1.36937 * 120.49372 * 111.40039 * 16 15 13 18 18 H 1.08001 * 120.00051 * 155.30993 * 17 16 15 19 19 C 1.38811 * 120.00022 * 335.30609 * 17 16 15 20 20 N 1.31614 * 120.00004 * 340.86507 * 19 17 16 21 21 Si 1.86799 * 119.99944 * 160.86761 * 19 17 16 22 22 H 1.48504 * 109.47185 * 90.00319 * 21 19 17 23 23 H 1.48499 * 109.47251 * 209.99595 * 21 19 17 24 24 H 1.48495 * 109.47289 * 330.00004 * 21 19 17 25 25 C 1.43428 * 118.87357 * 291.09372 * 16 15 13 26 26 C 1.55323 * 116.08867 * 89.55027 * 25 16 15 27 27 C 1.55192 * 108.87677 * 319.57310 * 26 25 16 28 28 O 1.39640 * 110.56596 * 322.14578 * 27 26 25 29 29 H 1.07997 * 120.00403 * 180.02562 * 1 2 3 30 30 H 1.08004 * 119.99709 * 0.02562 * 1 2 3 31 31 H 1.07998 * 120.00028 * 179.97438 * 2 1 3 32 32 H 1.08995 * 109.47051 * 299.99964 * 4 3 2 33 33 H 1.09007 * 109.46633 * 60.00044 * 4 3 2 34 34 H 1.08994 * 109.47497 * 179.97438 * 4 3 2 35 35 H 1.08995 * 109.47834 * 59.99740 * 5 3 2 36 36 H 1.09010 * 109.46978 * 180.02562 * 5 3 2 37 37 H 1.08999 * 109.47415 * 300.00499 * 6 5 3 38 38 H 1.08994 * 109.46851 * 60.00483 * 6 5 3 39 39 H 1.09005 * 109.46865 * 299.99574 * 8 7 6 40 40 H 1.08996 * 109.47553 * 59.99356 * 8 7 6 41 41 H 1.08996 * 109.46939 * 180.02562 * 8 7 6 42 42 H 0.97001 * 119.99717 * 359.97438 * 11 9 3 43 43 H 1.08997 * 109.47145 * 60.00300 * 12 11 9 44 44 H 1.08994 * 109.47171 * 299.99664 * 12 11 9 45 45 H 1.09003 * 109.98705 * 297.23377 * 15 13 12 46 46 H 1.08999 * 109.99256 * 58.60484 * 15 13 12 47 47 H 0.97005 * 120.00063 * 179.97438 * 20 19 17 48 48 H 1.09005 * 108.11937 * 327.90086 * 25 16 15 49 49 H 1.09000 * 108.12953 * 211.21042 * 25 16 15 50 50 H 1.08999 * 109.56721 * 79.38848 * 26 25 16 51 51 H 1.08994 * 109.57265 * 199.75347 * 26 25 16 52 52 H 1.09000 * 109.25615 * 82.41831 * 27 26 25 53 53 H 1.08995 * 109.25769 * 201.89876 * 27 26 25 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.3100 0.0000 0.0000 3 6 2.0635 1.3051 0.0000 4 6 2.9431 1.3862 1.2492 5 6 2.9431 1.3860 -1.2493 6 6 3.9374 0.2230 -1.2492 7 8 4.7593 0.2989 -2.4157 8 6 5.7275 -0.7475 -2.5136 9 6 1.0843 2.4505 -0.0006 10 8 1.0897 3.2525 -0.9104 11 7 0.2025 2.5838 1.0099 12 6 -0.7491 3.6976 1.0096 13 6 -1.6288 3.6172 2.2589 14 1 -0.9835 3.6381 3.1371 15 6 -2.5050 4.7960 2.2923 16 7 -3.3790 4.6466 3.5008 17 6 -3.2329 5.4868 4.5722 18 1 -4.0634 5.6585 5.2410 19 6 -2.0174 6.1176 4.7988 20 7 -0.9192 5.6482 4.2458 21 14 -1.9403 7.6327 5.8888 22 1 -1.7052 7.2190 7.2955 23 1 -0.8300 8.5102 5.4388 24 1 -3.2243 8.3732 5.7986 25 6 -4.3741 3.6137 3.5125 26 6 -3.9212 2.2351 4.0667 27 6 -2.4870 1.9351 3.5556 28 8 -2.3318 2.3792 2.2408 29 1 -0.5401 -0.9352 0.0004 30 1 -0.5400 0.9354 -0.0004 31 1 1.8500 -0.9353 0.0004 32 1 2.3165 1.3285 2.1392 33 1 3.6515 0.5576 1.2492 34 1 3.4885 2.3299 1.2490 35 1 2.3166 1.3282 -2.1393 36 1 3.4886 2.3298 -1.2491 37 1 4.5640 0.2808 -0.3592 38 1 3.3924 -0.7209 -1.2492 39 1 6.3911 -0.7104 -1.6496 40 1 5.2194 -1.7114 -2.5403 41 1 6.3110 -0.6185 -3.4252 42 1 0.1979 1.9422 1.7374 43 1 -1.3755 3.6401 0.1194 44 1 -0.2038 4.6413 1.0098 45 1 -1.9062 5.7032 2.3739 46 1 -3.1163 4.8316 1.3905 47 1 -0.0698 6.0893 4.4038 48 1 -4.7252 3.4696 2.4906 49 1 -5.2158 3.9603 4.1120 50 1 -3.9224 2.2615 5.1564 51 1 -4.6029 1.4595 3.7178 52 1 -1.7642 2.4430 4.1941 53 1 -2.3088 0.8605 3.5949 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 ZINC001150047440.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:11:06 Heat of formation + Delta-G solvation = -74.898455 kcal Electronic energy + Delta-G solvation = -31827.109233 eV Core-core repulsion = 27700.899133 eV Total energy + Delta-G solvation = -4126.210099 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.80 seconds Orbital eigenvalues (eV) -40.45732 -38.98911 -38.09241 -36.07926 -35.05750 -34.61296 -32.29791 -31.70161 -31.08579 -29.34594 -28.51587 -26.56111 -25.80364 -25.34131 -24.63024 -22.70807 -21.95695 -20.96131 -20.80048 -19.62962 -18.75081 -17.14671 -16.89602 -16.53793 -16.24257 -15.98631 -15.60389 -15.47863 -15.09234 -14.85464 -14.46237 -14.30269 -13.96192 -13.87352 -13.77441 -13.41472 -13.04662 -13.02643 -12.82227 -12.48890 -12.22446 -12.12633 -12.03265 -11.93184 -11.82495 -11.43036 -11.40911 -11.23309 -11.07137 -10.86226 -10.77653 -10.69243 -10.51124 -10.17559 -9.98881 -9.82789 -9.71236 -9.55825 -9.47575 -9.46667 -8.80989 -8.61591 -8.51365 -6.83211 -5.94215 -3.32770 2.17056 2.85826 3.27455 3.31189 3.39864 3.42426 3.70885 4.41809 4.55432 4.61057 4.70280 4.82046 4.93824 5.01663 5.09924 5.11546 5.20436 5.21908 5.23860 5.26946 5.33160 5.39165 5.62776 5.63355 5.72289 5.80466 5.85168 5.89696 5.91929 5.99698 6.06310 6.13449 6.19422 6.31805 6.39730 6.45406 6.52245 6.56659 6.60733 6.65674 6.69463 6.80119 6.91299 6.94428 6.94664 7.10933 7.31174 7.70177 7.80527 8.12304 8.30687 8.47624 9.05592 9.73676 10.27583 11.14791 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.016833 B = 0.002779 C = 0.002547 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1663.022552 B =10072.912844 C =10990.893066 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.184 4.184 2 C -0.160 4.160 3 C -0.046 4.046 4 C -0.142 4.142 5 C -0.090 4.090 6 C 0.060 3.940 7 O -0.391 6.391 8 C 0.032 3.968 9 C 0.522 3.478 10 O -0.542 6.542 11 N -0.711 5.711 12 C 0.119 3.881 13 C 0.061 3.939 14 H 0.107 0.893 15 C 0.144 3.856 16 N -0.587 5.587 17 C -0.248 4.248 18 H 0.068 0.932 19 C 0.002 3.998 20 N -0.907 5.907 21 Si 0.903 3.097 22 H -0.289 1.289 23 H -0.289 1.289 24 H -0.280 1.280 25 C 0.161 3.839 26 C -0.172 4.172 27 C 0.068 3.932 28 O -0.369 6.369 29 H 0.101 0.899 30 H 0.117 0.883 31 H 0.117 0.883 32 H 0.065 0.935 33 H 0.068 0.932 34 H 0.066 0.934 35 H 0.079 0.921 36 H 0.080 0.920 37 H 0.053 0.947 38 H 0.051 0.949 39 H 0.040 0.960 40 H 0.040 0.960 41 H 0.084 0.916 42 H 0.402 0.598 43 H 0.062 0.938 44 H 0.078 0.922 45 H 0.104 0.896 46 H 0.019 0.981 47 H 0.255 0.745 48 H 0.031 0.969 49 H 0.032 0.968 50 H 0.056 0.944 51 H 0.049 0.951 52 H 0.062 0.938 53 H 0.050 0.950 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 8.419 -14.039 -6.775 17.716 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.221 4.221 2 C -0.178 4.178 3 C -0.049 4.049 4 C -0.200 4.200 5 C -0.128 4.128 6 C -0.016 4.016 7 O -0.310 6.310 8 C -0.065 4.065 9 C 0.308 3.692 10 O -0.419 6.419 11 N -0.362 5.362 12 C -0.006 4.006 13 C 0.002 3.998 14 H 0.125 0.875 15 C 0.019 3.981 16 N -0.308 5.308 17 C -0.378 4.378 18 H 0.086 0.914 19 C -0.196 4.196 20 N -0.552 5.552 21 Si 0.733 3.267 22 H -0.216 1.216 23 H -0.215 1.215 24 H -0.206 1.206 25 C 0.032 3.968 26 C -0.213 4.213 27 C -0.010 4.010 28 O -0.287 6.287 29 H 0.119 0.881 30 H 0.135 0.865 31 H 0.135 0.865 32 H 0.084 0.916 33 H 0.087 0.913 34 H 0.085 0.915 35 H 0.097 0.903 36 H 0.099 0.901 37 H 0.071 0.929 38 H 0.069 0.931 39 H 0.058 0.942 40 H 0.058 0.942 41 H 0.103 0.897 42 H 0.237 0.763 43 H 0.080 0.920 44 H 0.096 0.904 45 H 0.122 0.878 46 H 0.037 0.963 47 H 0.065 0.935 48 H 0.050 0.950 49 H 0.051 0.949 50 H 0.075 0.925 51 H 0.068 0.932 52 H 0.080 0.920 53 H 0.068 0.932 Dipole moment (debyes) X Y Z Total from point charges 8.456 -13.273 -7.708 17.524 hybrid contribution 0.232 0.124 1.510 1.533 sum 8.688 -13.149 -6.198 16.935 Atomic orbital electron populations 1.24467 0.96101 1.03700 0.97824 1.22650 0.95839 0.96484 1.02866 1.19854 0.95113 0.92481 0.97454 1.21915 0.99446 1.02946 0.95650 1.21240 0.96183 0.98493 0.96842 1.22696 0.92536 0.98768 0.87642 1.87889 1.49284 1.49887 1.43967 1.22726 0.91566 0.90316 1.01870 1.20188 0.81093 0.85844 0.82114 1.90694 1.71146 1.44112 1.35975 1.46131 1.35076 1.30260 1.24769 1.21730 0.90316 0.89264 0.99241 1.22200 0.93767 0.85951 0.97858 0.87505 1.20705 0.92437 0.98024 0.86970 1.44031 1.39974 1.32211 1.14555 1.19244 0.95204 1.18787 1.04551 0.91387 1.24894 0.94866 1.00480 0.99319 1.69858 0.97664 1.39759 1.47969 0.86281 0.77261 0.83167 0.79956 1.21563 1.21529 1.20594 1.20367 0.91656 0.87136 0.97600 1.23579 0.97278 1.00172 1.00275 1.22329 0.97270 0.96891 0.84484 1.87509 1.77951 1.37358 1.25902 0.88067 0.86528 0.86522 0.91566 0.91332 0.91471 0.90253 0.90113 0.92934 0.93087 0.94152 0.94180 0.89712 0.76336 0.91986 0.90391 0.87846 0.96308 0.93489 0.95014 0.94919 0.92539 0.93215 0.91956 0.93159 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 18. 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.18 -2.13 13.17 29.02 0.38 -1.75 16 2 C -0.16 -1.69 7.99 -36.03 -0.29 -1.97 16 3 C -0.05 -0.55 0.62 -156.81 -0.10 -0.64 16 4 C -0.14 -1.45 8.00 37.16 0.30 -1.15 16 5 C -0.09 -1.10 3.92 -26.73 -0.10 -1.20 16 6 C 0.06 0.59 5.33 37.16 0.20 0.79 16 7 O -0.39 -4.42 10.70 -35.23 -0.38 -4.80 16 8 C 0.03 0.24 11.01 101.05 1.11 1.35 16 9 C 0.52 8.23 5.34 -10.99 -0.06 8.17 16 10 O -0.54 -10.52 15.98 5.55 0.09 -10.43 16 11 N -0.71 -11.11 3.60 -60.30 -0.22 -11.33 16 12 C 0.12 2.15 5.24 -4.04 -0.02 2.13 16 13 C 0.06 1.20 1.94 -29.64 -0.06 1.15 16 14 H 0.11 2.38 5.21 -51.93 -0.27 2.11 16 15 C 0.14 3.20 4.71 -5.06 -0.02 3.17 16 16 N -0.59 -13.84 2.59 -174.70 -0.45 -14.29 16 17 C -0.25 -6.38 9.61 -15.06 -0.14 -6.53 16 18 H 0.07 1.75 7.88 -52.49 -0.41 1.34 16 19 C 0.00 0.05 4.91 -17.43 -0.09 -0.03 16 20 N -0.91 -23.30 9.54 55.49 0.53 -22.77 16 21 Si 0.90 20.69 29.59 -169.99 -5.03 15.66 16 22 H -0.29 -6.66 7.11 56.52 0.40 -6.26 16 23 H -0.29 -6.73 7.11 56.52 0.40 -6.33 16 24 H -0.28 -6.43 7.11 56.52 0.40 -6.03 16 25 C 0.16 3.35 5.34 -4.43 -0.02 3.32 16 26 C -0.17 -3.13 6.75 -24.34 -0.16 -3.30 16 27 C 0.07 1.23 6.38 38.32 0.24 1.48 16 28 O -0.37 -6.94 9.54 -31.94 -0.30 -7.24 16 29 H 0.10 1.00 8.06 -52.49 -0.42 0.58 16 30 H 0.12 1.64 5.18 -52.48 -0.27 1.37 16 31 H 0.12 0.97 7.00 -52.49 -0.37 0.60 16 32 H 0.07 0.67 7.57 -51.93 -0.39 0.28 16 33 H 0.07 0.54 6.43 -51.93 -0.33 0.21 16 34 H 0.07 0.72 8.14 -51.93 -0.42 0.30 16 35 H 0.08 1.09 8.10 -51.93 -0.42 0.67 16 36 H 0.08 1.08 8.10 -51.92 -0.42 0.66 16 37 H 0.05 0.45 6.43 -51.93 -0.33 0.12 16 38 H 0.05 0.44 6.62 -51.93 -0.34 0.10 16 39 H 0.04 0.26 8.14 -51.93 -0.42 -0.16 16 40 H 0.04 0.26 8.14 -51.93 -0.42 -0.17 16 41 H 0.08 0.53 8.14 -51.93 -0.42 0.11 16 42 H 0.40 5.78 7.08 -40.82 -0.29 5.49 16 43 H 0.06 1.15 8.14 -51.93 -0.42 0.72 16 44 H 0.08 1.51 8.14 -51.93 -0.42 1.09 16 45 H 0.10 2.54 6.03 -51.93 -0.31 2.22 16 46 H 0.02 0.41 8.14 -51.93 -0.42 -0.02 16 47 H 0.26 6.49 8.31 -40.82 -0.34 6.15 16 48 H 0.03 0.64 7.99 -51.93 -0.41 0.23 16 49 H 0.03 0.67 7.96 -51.93 -0.41 0.26 16 50 H 0.06 1.04 8.14 -51.93 -0.42 0.62 16 51 H 0.05 0.80 8.14 -51.93 -0.42 0.38 16 52 H 0.06 1.28 6.66 -51.93 -0.35 0.93 16 53 H 0.05 0.78 8.14 -51.93 -0.42 0.36 16 LS Contribution 405.14 15.07 6.11 6.11 Total: -1.00 -28.53 405.14 -7.62 -36.16 By element: Atomic # 1 Polarization: 17.08 SS G_CDS: -9.13 Total: 7.95 kcal Atomic # 6 Polarization: 3.82 SS G_CDS: 1.16 Total: 4.98 kcal Atomic # 7 Polarization: -48.25 SS G_CDS: -0.14 Total: -48.39 kcal Atomic # 8 Polarization: -21.88 SS G_CDS: -0.59 Total: -22.47 kcal Atomic # 14 Polarization: 20.69 SS G_CDS: -5.03 Total: 15.66 kcal Total LS contribution 6.11 Total: 6.11 kcal Total: -28.53 -7.62 -36.16 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. ZINC001150047440.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 -38.740 kcal (2) G-P(sol) polarization free energy of solvation -28.534 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -67.274 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.624 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -36.158 kcal (6) G-S(sol) free energy of system = (1) + (5) -74.898 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.80 seconds