Wall clock time and date at job start Tue Mar 31 2020 19:41:43 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.53003 * 1 3 3 C 1.53000 * 109.47372 * 2 1 4 4 C 1.53001 * 109.47138 * 240.00404 * 2 1 3 5 5 C 1.50699 * 109.47240 * 59.99719 * 4 2 1 6 6 H 1.07997 * 119.99965 * 0.02562 * 5 4 2 7 7 C 1.37883 * 120.00460 * 180.02562 * 5 4 2 8 8 N 1.31509 * 120.00301 * 180.02562 * 7 5 4 9 9 Si 1.86801 * 119.99441 * 359.97438 * 7 5 4 10 10 H 1.48506 * 109.46632 * 90.00181 * 9 7 5 11 11 H 1.48496 * 109.47096 * 209.99990 * 9 7 5 12 12 H 1.48502 * 109.47428 * 330.00541 * 9 7 5 13 13 C 1.52994 * 109.47053 * 120.00227 * 2 1 3 14 14 N 1.46492 * 109.47266 * 299.99877 * 13 2 1 15 15 C 1.34776 * 120.00086 * 180.02562 * 14 13 2 16 16 O 1.21282 * 120.00274 * 359.97438 * 15 14 13 17 17 C 1.50699 * 120.00181 * 180.02562 * 15 14 13 18 18 H 1.09000 * 109.49910 * 60.06007 * 17 15 14 19 19 C 1.53038 * 109.49914 * 299.99782 * 17 15 14 20 20 C 1.53035 * 109.54175 * 178.68279 * 19 17 15 21 21 C 1.53191 * 109.31327 * 61.36680 * 20 19 17 22 22 N 1.46926 * 108.89109 * 305.29292 * 21 20 19 23 23 C 1.34778 * 120.63336 * 233.68039 * 22 21 20 24 24 O 1.21273 * 120.00022 * 180.02562 * 23 22 21 25 25 C 1.50706 * 119.99581 * 0.02562 * 23 22 21 26 26 C 1.53001 * 117.49617 * 248.74453 * 25 23 22 27 27 C 1.52998 * 117.49369 * 180.08940 * 25 23 22 28 28 C 1.47207 * 118.55531 * 53.54119 * 22 21 20 29 29 H 1.08996 * 109.47043 * 299.99747 * 1 2 3 30 30 H 1.09003 * 109.46745 * 59.99734 * 1 2 3 31 31 H 1.08997 * 109.47135 * 179.97438 * 1 2 3 32 32 H 1.08991 * 109.47354 * 179.97438 * 3 2 1 33 33 H 1.09002 * 109.47435 * 300.00647 * 3 2 1 34 34 H 1.09006 * 109.46884 * 60.00258 * 3 2 1 35 35 H 1.08992 * 109.47646 * 180.02562 * 4 2 1 36 36 H 1.09009 * 109.46943 * 300.00284 * 4 2 1 37 37 H 0.96996 * 119.99895 * 179.97438 * 8 7 5 38 38 H 1.09002 * 109.46809 * 179.97438 * 13 2 1 39 39 H 1.09005 * 109.47129 * 60.00241 * 13 2 1 40 40 H 0.96999 * 120.00102 * 359.97438 * 14 13 2 41 41 H 1.09003 * 109.46001 * 58.65612 * 19 17 15 42 42 H 1.09004 * 109.46331 * 298.70671 * 19 17 15 43 43 H 1.09003 * 109.49576 * 301.41349 * 20 19 17 44 44 H 1.09000 * 109.52700 * 181.33397 * 20 19 17 45 45 H 1.09002 * 109.58786 * 65.15336 * 21 20 19 46 46 H 1.08996 * 109.58870 * 185.43616 * 21 20 19 47 47 H 1.08999 * 115.55287 * 34.41336 * 25 23 22 48 48 H 1.09008 * 117.49649 * 0.02562 * 26 25 23 49 49 H 1.08997 * 117.49980 * 145.02365 * 26 25 23 50 50 H 1.09001 * 117.49584 * 214.98027 * 27 25 23 51 51 H 1.09002 * 117.50282 * 359.97039 * 27 25 23 52 52 H 1.08996 * 109.53610 * 186.67231 * 28 22 21 53 53 H 1.09002 * 109.68003 * 66.36699 * 28 22 21 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.0401 1.4425 0.0000 4 6 2.0400 -0.7212 -1.2493 5 6 1.5377 -0.0107 -2.4797 6 1 0.9171 0.8675 -2.3806 7 6 1.8703 -0.4820 -3.7321 8 7 1.4323 0.1383 -4.8058 9 14 2.9443 -2.0007 -3.9035 10 1 2.0857 -3.2116 -3.9484 11 1 3.7357 -1.9084 -5.1566 12 1 3.8651 -2.0896 -2.7418 13 6 2.0400 -0.7213 1.2492 14 7 1.5517 -0.0308 2.4453 15 6 1.8777 -0.4910 3.6694 16 8 2.5750 -1.4771 3.7807 17 6 1.3750 0.2190 4.8999 18 1 0.2850 0.2175 4.9010 19 6 1.8837 1.6624 4.9006 20 6 1.3468 2.3938 6.1329 21 6 1.8607 1.6995 7.3981 22 7 1.5000 0.2764 7.3397 23 6 0.8248 -0.3030 8.3521 24 8 0.5333 -1.4786 8.2921 25 6 0.4322 0.5102 9.5587 26 6 1.1322 0.1939 10.8819 27 6 -0.3309 -0.2056 10.6750 28 6 1.8892 -0.5000 6.1511 29 1 -0.3633 0.5138 0.8900 30 1 -0.3633 0.5139 -0.8900 31 1 -0.3633 -1.0276 -0.0005 32 1 3.1300 1.4425 0.0005 33 1 1.6767 1.9564 -0.8899 34 1 1.6768 1.9563 0.8901 35 1 3.1300 -0.7216 -1.2492 36 1 1.6766 -1.7489 -1.2494 37 1 1.6666 -0.1930 -5.6868 38 1 3.1300 -0.7216 1.2489 39 1 1.6767 -1.7490 1.2491 40 1 0.9944 0.7581 2.3564 41 1 1.5397 2.1691 3.9989 42 1 2.9735 1.6633 4.9239 43 1 0.2570 2.3718 6.1229 44 1 1.6899 3.4283 6.1208 45 1 2.9446 1.7997 7.4552 46 1 1.4048 2.1563 8.2764 47 1 0.1760 1.5510 9.3611 48 1 1.8770 -0.6020 10.8796 49 1 1.3361 1.0266 11.5550 50 1 -1.0890 0.3644 11.2120 51 1 -0.5493 -1.2645 10.5365 52 1 1.4501 -1.4960 6.2066 53 1 2.9752 -0.5814 6.1053 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC000359767610.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 19:41:43 Heat of formation + Delta-G solvation = -67.031258 kcal Electronic energy + Delta-G solvation = -29696.659678 eV Core-core repulsion = 25761.246778 eV Total energy + Delta-G solvation = -3935.412901 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.86 seconds Orbital eigenvalues (eV) -40.89766 -39.80594 -37.94119 -37.08478 -36.36483 -33.89281 -33.17149 -32.07149 -31.61527 -27.91031 -26.92455 -26.59776 -26.13448 -24.79770 -23.37771 -22.85435 -22.42344 -21.45390 -20.02956 -19.18892 -18.37490 -17.92668 -17.04360 -16.78824 -16.16214 -16.02561 -15.54901 -15.26760 -14.94995 -14.77266 -14.50632 -14.22172 -14.14034 -13.78298 -13.51077 -13.42709 -13.23952 -12.97831 -12.81753 -12.74191 -12.30859 -12.27281 -12.01334 -11.85660 -11.57127 -11.39348 -11.27702 -11.15698 -11.03030 -11.01654 -11.00888 -10.89527 -10.64297 -10.52064 -10.11590 -10.00117 -9.89646 -9.75687 -9.46070 -9.03280 -8.98692 -8.84242 -8.55303 -6.06116 -4.03483 2.14530 2.50501 2.67189 3.30498 3.33475 3.38630 3.42560 3.43939 4.01120 4.21258 4.36557 4.38141 4.44138 4.44946 4.45431 4.58096 4.61218 4.81229 5.01027 5.01571 5.14242 5.16620 5.26367 5.27031 5.34774 5.37362 5.41573 5.47747 5.49724 5.53761 5.58673 5.60989 5.69947 5.77813 5.86747 6.06686 6.08488 6.10576 6.16499 6.27730 6.44381 6.55953 6.65261 6.86853 6.98223 7.01951 7.03530 7.17777 7.34138 7.60532 7.77258 7.85286 8.29631 8.50748 8.97659 9.56596 11.06518 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.028998 B = 0.002028 C = 0.001974 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 965.347999 B =13801.383063 C =14182.521527 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.129 4.129 2 C -0.063 4.063 3 C -0.128 4.128 4 C -0.002 4.002 5 C -0.511 4.511 6 H 0.078 0.922 7 C 0.066 3.934 8 N -0.903 5.903 9 Si 0.875 3.125 10 H -0.278 1.278 11 H -0.289 1.289 12 H -0.265 1.265 13 C 0.125 3.875 14 N -0.718 5.718 15 C 0.516 3.484 16 O -0.543 6.543 17 C -0.122 4.122 18 H 0.100 0.900 19 C -0.113 4.113 20 C -0.134 4.134 21 C 0.104 3.896 22 N -0.609 5.609 23 C 0.544 3.456 24 O -0.538 6.538 25 C -0.197 4.197 26 C -0.152 4.152 27 C -0.148 4.148 28 C 0.133 3.867 29 H 0.030 0.970 30 H 0.073 0.927 31 H 0.047 0.953 32 H 0.047 0.953 33 H 0.073 0.927 34 H 0.030 0.970 35 H 0.018 0.982 36 H 0.014 0.986 37 H 0.262 0.738 38 H 0.064 0.936 39 H 0.064 0.936 40 H 0.399 0.601 41 H 0.075 0.925 42 H 0.073 0.927 43 H 0.071 0.929 44 H 0.083 0.917 45 H 0.072 0.928 46 H 0.089 0.911 47 H 0.126 0.874 48 H 0.099 0.901 49 H 0.105 0.895 50 H 0.103 0.897 51 H 0.108 0.892 52 H 0.099 0.901 53 H 0.079 0.921 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -1.313 7.945 29.949 31.013 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.186 4.186 2 C -0.063 4.063 3 C -0.185 4.185 4 C -0.040 4.040 5 C -0.549 4.549 6 H 0.096 0.904 7 C -0.135 4.135 8 N -0.542 5.542 9 Si 0.702 3.298 10 H -0.204 1.204 11 H -0.215 1.215 12 H -0.190 1.190 13 C 0.002 3.998 14 N -0.369 5.369 15 C 0.301 3.699 16 O -0.421 6.421 17 C -0.144 4.144 18 H 0.118 0.882 19 C -0.151 4.151 20 C -0.173 4.173 21 C -0.019 4.019 22 N -0.342 5.342 23 C 0.333 3.667 24 O -0.414 6.414 25 C -0.217 4.217 26 C -0.190 4.190 27 C -0.185 4.185 28 C 0.011 3.989 29 H 0.049 0.951 30 H 0.092 0.908 31 H 0.066 0.934 32 H 0.066 0.934 33 H 0.091 0.909 34 H 0.049 0.951 35 H 0.036 0.964 36 H 0.032 0.968 37 H 0.072 0.928 38 H 0.082 0.918 39 H 0.082 0.918 40 H 0.235 0.765 41 H 0.094 0.906 42 H 0.092 0.908 43 H 0.089 0.911 44 H 0.101 0.899 45 H 0.091 0.909 46 H 0.107 0.893 47 H 0.143 0.857 48 H 0.118 0.882 49 H 0.123 0.877 50 H 0.122 0.878 51 H 0.126 0.874 52 H 0.117 0.883 53 H 0.097 0.903 Dipole moment (debyes) X Y Z Total from point charges -0.951 6.584 29.779 30.513 hybrid contribution 0.287 -0.260 -0.560 0.681 sum -0.664 6.325 29.218 29.902 Atomic orbital electron populations 1.21781 0.94805 1.01309 1.00746 1.20691 0.95811 0.95906 0.93939 1.21768 1.00440 0.95584 1.00728 1.19591 0.96282 0.95511 0.92602 1.21070 1.31892 1.13343 0.88611 0.90406 1.25003 0.95706 0.96759 0.95998 1.69896 1.45907 1.40544 0.97885 0.86857 0.81333 0.83918 0.77651 1.20365 1.21496 1.19002 1.21653 0.98908 0.96090 0.83145 1.46156 1.53125 1.31936 1.05720 1.20300 0.80550 0.84541 0.84549 1.90716 1.38114 1.26374 1.86868 1.21662 1.02597 0.96500 0.93660 0.88206 1.21538 1.01676 0.92989 0.98865 1.21937 1.01634 0.99597 0.94086 1.21838 1.00785 0.80296 0.98996 1.48094 1.56417 1.10917 1.18804 1.19677 0.76821 0.87152 0.83081 1.90713 1.54605 1.16399 1.79733 1.22233 1.09119 1.00908 0.89488 1.22954 0.96498 1.03771 0.95729 1.22982 0.94602 1.02272 0.98688 1.21386 1.00671 0.95590 0.81300 0.95062 0.90845 0.93376 0.93392 0.90854 0.95082 0.96365 0.96814 0.92783 0.91777 0.91794 0.76518 0.90635 0.90814 0.91061 0.89895 0.90936 0.89326 0.85668 0.88215 0.87714 0.87842 0.87358 0.88262 0.90268 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.13 -2.19 8.05 37.16 0.30 -1.89 16 2 C -0.06 -1.12 0.54 -154.51 -0.08 -1.21 16 3 C -0.13 -2.12 8.05 37.16 0.30 -1.82 16 4 C 0.00 -0.05 3.11 -27.88 -0.09 -0.13 16 5 C -0.51 -13.55 8.20 -34.44 -0.28 -13.83 16 6 H 0.08 2.18 5.92 -52.49 -0.31 1.86 16 7 C 0.07 1.87 5.47 -17.82 -0.10 1.77 16 8 N -0.90 -27.04 13.96 55.43 0.77 -26.27 16 9 Si 0.87 21.47 27.47 -169.99 -4.67 16.80 16 10 H -0.28 -6.77 7.11 56.52 0.40 -6.37 16 11 H -0.29 -7.08 7.11 56.52 0.40 -6.68 16 12 H -0.27 -6.31 6.54 56.52 0.37 -5.94 16 13 C 0.13 2.04 4.56 -4.05 -0.02 2.03 16 14 N -0.72 -9.48 4.85 -60.30 -0.29 -9.77 16 15 C 0.52 7.05 7.16 -10.99 -0.08 6.97 16 16 O -0.54 -9.63 15.92 5.56 0.09 -9.54 16 17 C -0.12 -1.21 2.68 -91.65 -0.25 -1.46 16 18 H 0.10 0.94 8.14 -51.93 -0.42 0.52 16 19 C -0.11 -0.82 5.07 -26.69 -0.14 -0.95 16 20 C -0.13 -0.68 6.08 -26.61 -0.16 -0.85 16 21 C 0.10 0.57 6.34 -3.71 -0.02 0.55 16 22 N -0.61 -5.54 2.98 -172.54 -0.51 -6.05 16 23 C 0.54 5.78 7.49 -10.98 -0.08 5.70 16 24 O -0.54 -7.99 15.69 5.56 0.09 -7.90 16 25 C -0.20 -1.35 5.12 -91.77 -0.47 -1.82 16 26 C -0.15 -0.86 10.24 -26.69 -0.27 -1.13 16 27 C -0.15 -0.97 10.05 -26.70 -0.27 -1.23 16 28 C 0.13 1.49 5.08 -3.56 -0.02 1.47 16 29 H 0.03 0.42 6.86 -51.93 -0.36 0.07 16 30 H 0.07 1.40 6.62 -51.93 -0.34 1.06 16 31 H 0.05 0.83 8.14 -51.93 -0.42 0.40 16 32 H 0.05 0.80 8.14 -51.93 -0.42 0.38 16 33 H 0.07 1.38 6.62 -51.93 -0.34 1.04 16 34 H 0.03 0.41 6.86 -51.93 -0.36 0.05 16 35 H 0.02 0.42 7.20 -51.93 -0.37 0.05 16 36 H 0.01 0.31 7.57 -51.92 -0.39 -0.08 16 37 H 0.26 7.43 8.31 -40.82 -0.34 7.09 16 38 H 0.06 1.11 8.14 -51.93 -0.42 0.68 16 39 H 0.06 1.12 8.14 -51.93 -0.42 0.70 16 40 H 0.40 4.40 5.56 -40.82 -0.23 4.18 16 41 H 0.07 0.51 7.53 -51.93 -0.39 0.11 16 42 H 0.07 0.57 8.14 -51.93 -0.42 0.15 16 43 H 0.07 0.37 8.14 -51.93 -0.42 -0.05 16 44 H 0.08 0.29 8.14 -51.93 -0.42 -0.13 16 45 H 0.07 0.37 8.14 -51.93 -0.42 -0.06 16 46 H 0.09 0.32 5.51 -51.93 -0.29 0.03 16 47 H 0.13 0.56 6.57 -51.93 -0.34 0.22 16 48 H 0.10 0.67 8.14 -51.92 -0.42 0.24 16 49 H 0.10 0.36 8.14 -51.93 -0.42 -0.06 16 50 H 0.10 0.46 8.14 -51.93 -0.42 0.04 16 51 H 0.11 0.91 7.67 -51.93 -0.40 0.51 16 52 H 0.10 1.38 7.02 -51.93 -0.36 1.02 16 53 H 0.08 0.93 7.95 -51.93 -0.41 0.51 16 LS Contribution 406.38 15.07 6.12 6.12 Total: -1.00 -33.65 406.38 -9.27 -42.91 By element: Atomic # 1 Polarization: 10.68 SS G_CDS: -9.14 Total: 1.54 kcal Atomic # 6 Polarization: -6.12 SS G_CDS: -1.73 Total: -7.85 kcal Atomic # 7 Polarization: -42.06 SS G_CDS: -0.03 Total: -42.09 kcal Atomic # 8 Polarization: -17.62 SS G_CDS: 0.18 Total: -17.44 kcal Atomic # 14 Polarization: 21.47 SS G_CDS: -4.67 Total: 16.80 kcal Total LS contribution 6.12 Total: 6.12 kcal Total: -33.65 -9.27 -42.91 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. ZINC000359767610.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 -24.117 kcal (2) G-P(sol) polarization free energy of solvation -33.645 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -57.763 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.269 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.914 kcal (6) G-S(sol) free energy of system = (1) + (5) -67.031 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.87 seconds