Wall clock time and date at job start Tue Mar 31 2020 23:18:16 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.52996 * 1 3 3 C 1.53008 * 109.47276 * 2 1 4 4 O 1.42891 * 109.47169 * 64.04126 * 3 2 1 5 5 C 1.53000 * 109.47044 * 184.04257 * 3 2 1 6 6 C 1.52997 * 109.47165 * 304.99970 * 5 3 2 7 7 C 1.53009 * 109.46647 * 304.03950 * 3 2 1 8 8 C 1.53001 * 109.46608 * 184.12831 * 7 3 2 9 9 N 1.46497 * 109.46896 * 175.17281 * 8 7 3 10 10 C 1.34781 * 120.00114 * 180.27540 * 9 8 7 11 11 O 1.21350 * 119.99661 * 359.72657 * 10 9 8 12 12 C 1.49449 * 120.00055 * 179.72563 * 10 9 8 13 13 O 1.21345 * 120.00120 * 6.64420 * 12 10 9 14 14 N 1.34776 * 119.99665 * 186.64662 * 12 10 9 15 15 C 1.46924 * 120.63181 * 6.54976 * 14 12 10 16 16 C 1.52755 * 109.00735 * 233.58331 * 15 14 12 17 17 C 1.53046 * 109.38679 * 305.59087 * 16 15 14 18 18 C 1.50699 * 109.45013 * 181.25815 * 17 16 15 19 19 H 1.08003 * 119.99593 * 239.94244 * 18 17 16 20 20 C 1.37867 * 120.00421 * 59.93461 * 18 17 16 21 21 N 1.31520 * 120.00256 * 0.02562 * 20 18 17 22 22 Si 1.86803 * 120.00136 * 180.02562 * 20 18 17 23 23 H 1.48498 * 109.47369 * 59.99996 * 22 20 18 24 24 H 1.48497 * 109.47374 * 180.02562 * 22 20 18 25 25 H 1.48504 * 109.47063 * 299.99858 * 22 20 18 26 26 C 1.53037 * 109.53715 * 61.24014 * 17 16 15 27 27 C 1.46935 * 120.62709 * 186.83117 * 14 12 10 28 28 H 1.09006 * 109.47002 * 55.15992 * 1 2 3 29 29 H 1.08994 * 109.47189 * 175.16284 * 1 2 3 30 30 H 1.09002 * 109.46635 * 295.16575 * 1 2 3 31 31 H 1.08996 * 109.47330 * 240.00086 * 2 1 3 32 32 H 1.09002 * 109.47537 * 119.99917 * 2 1 3 33 33 H 0.96705 * 114.00087 * 180.02562 * 4 3 2 34 34 H 1.09001 * 109.47074 * 64.99959 * 5 3 2 35 35 H 1.09001 * 109.46851 * 184.99511 * 5 3 2 36 36 H 1.08999 * 109.46671 * 299.99690 * 6 5 3 37 37 H 1.08994 * 109.47262 * 60.00105 * 6 5 3 38 38 H 1.09002 * 109.47118 * 180.02562 * 6 5 3 39 39 H 1.09000 * 109.46837 * 64.13112 * 7 3 2 40 40 H 1.08991 * 109.47177 * 304.13011 * 7 3 2 41 41 H 1.09000 * 109.47112 * 55.17158 * 8 7 3 42 42 H 1.09004 * 109.47474 * 295.17425 * 8 7 3 43 43 H 0.97006 * 120.00221 * 0.02562 * 9 8 7 44 44 H 1.08993 * 109.58887 * 353.42017 * 15 14 12 45 45 H 1.09006 * 109.58382 * 113.84681 * 15 14 12 46 46 H 1.09004 * 109.56819 * 185.55554 * 16 15 14 47 47 H 1.09001 * 109.39087 * 65.51666 * 16 15 14 48 48 H 1.09004 * 109.45667 * 301.21652 * 17 16 15 49 49 H 0.96999 * 120.00232 * 180.02562 * 21 20 18 50 50 H 1.08997 * 109.49715 * 178.68893 * 26 17 16 51 51 H 1.09000 * 109.52488 * 58.60098 * 26 17 16 52 52 H 1.08997 * 109.58312 * 246.15329 * 27 14 12 53 53 H 1.08998 * 109.58517 * 6.43849 * 27 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 6 2.0400 1.4426 0.0000 4 8 1.6428 2.0882 1.2113 5 6 3.5666 1.4437 -0.1017 6 6 4.1534 0.5892 1.0236 7 6 1.4487 2.1925 -1.1954 8 6 2.0380 3.6032 -1.2546 9 7 1.5610 4.2814 -2.4623 10 6 1.9572 5.5419 -2.7285 11 8 2.7051 6.1167 -1.9651 12 6 1.4755 6.2315 -3.9637 13 8 0.6356 5.7055 -4.6639 14 7 1.9739 7.4377 -4.2999 15 6 3.0905 8.0231 -3.5453 16 6 2.6962 9.4234 -3.0795 17 6 2.2487 10.2528 -4.2855 18 6 1.8888 11.6441 -3.8319 19 1 0.8895 12.0220 -3.9901 20 6 2.8350 12.4345 -3.2150 21 7 4.0521 11.9745 -3.0226 22 14 2.3886 14.1589 -2.6522 23 1 1.2824 14.0863 -1.6642 24 1 3.5732 14.7968 -2.0237 25 1 1.9554 14.9648 -3.8219 26 6 1.0280 9.5960 -4.9338 27 6 1.4070 8.1950 -5.4242 28 1 -0.3633 0.5871 -0.8435 29 1 -0.3633 -1.0239 -0.0867 30 1 -0.3633 0.4370 0.9302 31 1 1.8933 -0.5138 -0.8899 32 1 1.8934 -0.5138 0.8900 33 1 1.9324 3.0084 1.2798 34 1 3.8669 1.0316 -1.0651 35 1 3.9359 2.4654 -0.0125 36 1 3.8531 1.0013 1.9870 37 1 3.7841 -0.4324 0.9344 38 1 5.2410 0.5904 0.9514 39 1 0.3661 2.2553 -1.0855 40 1 1.6902 1.6590 -2.1146 41 1 3.1260 3.5417 -1.2781 42 1 1.7248 4.1651 -0.3746 43 1 0.9629 3.8220 -3.0724 44 1 3.3123 7.4002 -2.6789 45 1 3.9700 8.0883 -4.1859 46 1 3.5512 9.9044 -2.6044 47 1 1.8766 9.3496 -2.3646 48 1 3.0606 10.3060 -5.0109 49 1 4.7177 12.5305 -2.5882 50 1 0.6957 10.1995 -5.7784 51 1 0.2236 9.5200 -4.2021 52 1 2.1465 8.2745 -6.2210 53 1 0.5199 7.6842 -5.7985 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 ZINC001325514432.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 23:18:16 Heat of formation + Delta-G solvation = -126.252955 kcal Electronic energy + Delta-G solvation = -30671.123264 eV Core-core repulsion = 26542.686266 eV Total energy + Delta-G solvation = -4128.436998 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.75 seconds Orbital eigenvalues (eV) -40.76323 -39.70275 -38.42428 -37.03875 -35.30914 -33.97507 -33.43461 -32.21968 -31.21938 -30.81778 -30.06005 -26.92388 -26.10404 -24.53753 -23.35814 -23.09168 -22.38368 -21.83430 -20.50578 -20.04888 -18.86984 -17.43039 -17.09212 -16.57677 -16.28563 -16.06287 -15.80895 -15.55838 -15.07475 -14.94851 -14.66859 -14.52209 -14.16720 -14.10194 -14.06933 -13.54229 -13.05544 -12.93238 -12.83402 -12.70246 -12.49429 -12.33537 -12.26375 -12.17656 -12.02511 -11.73249 -11.57822 -11.40898 -11.28523 -11.17665 -10.93467 -10.86387 -10.63353 -10.61017 -10.35075 -10.21185 -10.05116 -9.69596 -9.38249 -9.23729 -9.16841 -8.52480 -8.30976 -8.09605 -5.99470 -4.06004 1.57890 2.46404 3.31654 3.36453 3.38508 3.39489 3.69612 3.84380 4.18883 4.30479 4.40910 4.43135 4.52977 4.66387 4.75433 4.87886 5.04725 5.07497 5.11908 5.24790 5.30157 5.35757 5.38674 5.55046 5.57088 5.63006 5.63627 5.68864 5.79379 5.81521 5.84570 5.88378 5.94813 6.01888 6.04079 6.09560 6.24275 6.46766 6.62201 6.71887 6.79900 6.97381 7.09044 7.18343 7.31095 7.32492 7.67143 7.79071 7.88874 7.97283 8.09178 8.23174 8.43398 8.96889 9.58693 11.05961 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.017576 B = 0.002179 C = 0.002047 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1592.671026 B =12848.349203 C =13675.204562 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.154 4.154 2 C -0.110 4.110 3 C 0.141 3.859 4 O -0.555 6.555 5 C -0.156 4.156 6 C -0.151 4.151 7 C -0.163 4.163 8 C 0.110 3.890 9 N -0.704 5.704 10 C 0.510 3.490 11 O -0.499 6.499 12 C 0.499 3.501 13 O -0.546 6.546 14 N -0.559 5.559 15 C 0.044 3.956 16 C -0.118 4.118 17 C 0.045 3.955 18 C -0.513 4.513 19 H 0.054 0.946 20 C 0.061 3.939 21 N -0.890 5.890 22 Si 0.891 3.109 23 H -0.278 1.278 24 H -0.287 1.287 25 H -0.278 1.278 26 C -0.123 4.123 27 C 0.107 3.893 28 H 0.057 0.943 29 H 0.051 0.949 30 H 0.062 0.938 31 H 0.073 0.927 32 H 0.074 0.926 33 H 0.377 0.623 34 H 0.072 0.928 35 H 0.072 0.928 36 H 0.061 0.939 37 H 0.056 0.944 38 H 0.052 0.948 39 H 0.086 0.914 40 H 0.081 0.919 41 H 0.077 0.923 42 H 0.072 0.928 43 H 0.412 0.588 44 H 0.140 0.860 45 H 0.056 0.944 46 H 0.105 0.895 47 H 0.036 0.964 48 H 0.027 0.973 49 H 0.261 0.739 50 H 0.063 0.937 51 H 0.052 0.948 52 H 0.066 0.934 53 H 0.080 0.920 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -1.965 -24.425 4.669 24.945 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.211 4.211 2 C -0.147 4.147 3 C 0.100 3.900 4 O -0.362 6.362 5 C -0.194 4.194 6 C -0.208 4.208 7 C -0.201 4.201 8 C -0.013 4.013 9 N -0.358 5.358 10 C 0.292 3.708 11 O -0.375 6.375 12 C 0.281 3.719 13 O -0.426 6.426 14 N -0.290 5.290 15 C -0.075 4.075 16 C -0.156 4.156 17 C 0.027 3.973 18 C -0.551 4.551 19 H 0.072 0.928 20 C -0.140 4.140 21 N -0.529 5.529 22 Si 0.719 3.281 23 H -0.203 1.203 24 H -0.213 1.213 25 H -0.204 1.204 26 C -0.161 4.161 27 C -0.015 4.015 28 H 0.076 0.924 29 H 0.070 0.930 30 H 0.081 0.919 31 H 0.092 0.908 32 H 0.093 0.907 33 H 0.225 0.775 34 H 0.091 0.909 35 H 0.090 0.910 36 H 0.080 0.920 37 H 0.075 0.925 38 H 0.071 0.929 39 H 0.104 0.896 40 H 0.100 0.900 41 H 0.095 0.905 42 H 0.091 0.909 43 H 0.248 0.752 44 H 0.156 0.844 45 H 0.075 0.925 46 H 0.123 0.877 47 H 0.055 0.945 48 H 0.045 0.955 49 H 0.071 0.929 50 H 0.082 0.918 51 H 0.071 0.929 52 H 0.085 0.915 53 H 0.098 0.902 Dipole moment (debyes) X Y Z Total from point charges -1.930 -24.782 5.042 25.363 hybrid contribution -0.152 1.514 -0.185 1.532 sum -2.082 -23.268 4.857 23.861 Atomic orbital electron populations 1.21810 0.95252 1.01409 1.02651 1.21551 0.96740 0.93696 1.02742 1.21131 0.92511 0.92898 0.83484 1.86543 1.80251 1.29988 1.39401 1.22331 0.95335 1.01555 1.00188 1.21760 1.00245 0.99874 0.98900 1.22614 1.02425 0.94935 1.00138 1.21578 0.99853 0.92261 0.87613 1.45760 1.50055 1.14180 1.25760 1.20463 0.81026 0.82135 0.87162 1.90406 1.34762 1.63076 1.49261 1.19825 0.83345 0.82948 0.85766 1.90692 1.33861 1.65381 1.52655 1.48247 1.33364 1.15961 1.31451 1.23395 0.86335 0.99543 0.98270 1.21852 1.02774 0.91107 0.99862 1.18690 0.91906 0.93417 0.93291 1.21146 0.97195 0.98391 1.38400 0.92800 1.24961 0.95004 0.98165 0.95840 1.69837 1.03046 1.36729 1.43307 0.86676 0.77789 0.84687 0.78936 1.20348 1.21310 1.20350 1.22019 1.00929 0.92792 1.00346 1.21895 0.98642 0.93311 0.87619 0.92436 0.92957 0.91918 0.90794 0.90706 0.77517 0.90947 0.90985 0.92027 0.92490 0.92861 0.89598 0.90012 0.90475 0.90947 0.75208 0.84358 0.92541 0.87737 0.94469 0.95500 0.92917 0.91849 0.92899 0.91532 0.90199 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.15 -0.99 9.44 37.16 0.35 -0.64 16 2 C -0.11 -0.67 4.59 -26.73 -0.12 -0.80 16 3 C 0.14 0.91 0.53 -90.61 -0.05 0.86 16 4 O -0.56 -4.26 12.09 -35.23 -0.43 -4.68 16 5 C -0.16 -0.97 4.57 -26.73 -0.12 -1.09 16 6 C -0.15 -0.94 9.32 37.16 0.35 -0.59 16 7 C -0.16 -1.05 4.03 -26.73 -0.11 -1.15 16 8 C 0.11 0.94 5.03 -4.04 -0.02 0.92 16 9 N -0.70 -7.78 5.55 -60.73 -0.34 -8.11 16 10 C 0.51 7.73 7.78 -11.58 -0.09 7.64 16 11 O -0.50 -8.56 13.10 5.50 0.07 -8.48 16 12 C 0.50 8.68 7.74 -11.59 -0.09 8.59 16 13 O -0.55 -9.91 15.94 5.50 0.09 -9.83 16 14 N -0.56 -10.45 2.98 -173.54 -0.52 -10.96 16 15 C 0.04 0.88 6.24 -3.94 -0.02 0.86 16 16 C -0.12 -2.75 5.14 -26.83 -0.14 -2.89 16 17 C 0.05 1.11 2.07 -91.73 -0.19 0.92 16 18 C -0.51 -13.90 8.56 -34.44 -0.29 -14.20 16 19 H 0.05 1.45 7.74 -52.48 -0.41 1.04 16 20 C 0.06 1.68 5.30 -17.82 -0.09 1.59 16 21 N -0.89 -24.98 11.31 55.43 0.63 -24.35 16 22 Si 0.89 21.18 29.54 -169.99 -5.02 16.16 16 23 H -0.28 -6.58 7.11 56.52 0.40 -6.17 16 24 H -0.29 -6.78 7.11 56.52 0.40 -6.38 16 25 H -0.28 -6.54 7.11 56.52 0.40 -6.14 16 26 C -0.12 -2.53 5.30 -26.61 -0.14 -2.67 16 27 C 0.11 1.95 6.42 -3.71 -0.02 1.93 16 28 H 0.06 0.36 6.18 -51.93 -0.32 0.03 16 29 H 0.05 0.29 8.14 -51.93 -0.42 -0.13 16 30 H 0.06 0.48 7.93 -51.93 -0.41 0.07 16 31 H 0.07 0.37 8.14 -51.93 -0.42 -0.05 16 32 H 0.07 0.48 6.59 -51.93 -0.34 0.14 16 33 H 0.38 2.18 7.60 45.56 0.35 2.52 16 34 H 0.07 0.42 8.14 -51.93 -0.42 -0.01 16 35 H 0.07 0.47 6.55 -51.93 -0.34 0.13 16 36 H 0.06 0.43 7.82 -51.93 -0.41 0.02 16 37 H 0.06 0.33 6.44 -51.93 -0.33 -0.01 16 38 H 0.05 0.31 8.14 -51.93 -0.42 -0.11 16 39 H 0.09 0.58 6.49 -51.93 -0.34 0.24 16 40 H 0.08 0.46 8.14 -51.93 -0.42 0.04 16 41 H 0.08 0.67 6.20 -51.93 -0.32 0.34 16 42 H 0.07 0.64 7.09 -51.93 -0.37 0.27 16 43 H 0.41 4.05 7.96 -40.82 -0.33 3.72 16 44 H 0.14 2.70 3.77 -56.81 -0.21 2.49 16 45 H 0.06 1.11 8.14 -51.93 -0.42 0.68 16 46 H 0.10 2.72 6.03 -51.93 -0.31 2.41 16 47 H 0.04 0.90 8.14 -51.93 -0.42 0.47 16 48 H 0.03 0.71 8.14 -51.93 -0.42 0.28 16 49 H 0.26 7.03 8.31 -40.82 -0.34 6.69 16 50 H 0.06 1.24 8.14 -51.93 -0.42 0.82 16 51 H 0.05 1.12 8.14 -51.93 -0.42 0.70 16 52 H 0.07 1.15 8.14 -51.93 -0.42 0.72 16 53 H 0.08 1.39 7.05 -51.93 -0.37 1.02 16 LS Contribution 403.18 15.07 6.08 6.08 Total: -1.00 -30.57 403.18 -8.49 -39.06 By element: Atomic # 1 Polarization: 14.11 SS G_CDS: -8.25 Total: 5.86 kcal Atomic # 6 Polarization: 0.08 SS G_CDS: -0.81 Total: -0.73 kcal Atomic # 7 Polarization: -43.20 SS G_CDS: -0.23 Total: -43.43 kcal Atomic # 8 Polarization: -22.73 SS G_CDS: -0.27 Total: -22.99 kcal Atomic # 14 Polarization: 21.18 SS G_CDS: -5.02 Total: 16.16 kcal Total LS contribution 6.08 Total: 6.08 kcal Total: -30.57 -8.49 -39.06 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. ZINC001325514432.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 -87.191 kcal (2) G-P(sol) polarization free energy of solvation -30.567 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -117.759 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.494 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.062 kcal (6) G-S(sol) free energy of system = (1) + (5) -126.253 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.75 seconds