Wall clock time and date at job start Tue Mar 31 2020 20:37:59 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.53002 * 1 3 3 C 1.53000 * 109.47482 * 2 1 4 4 H 1.09000 * 109.46713 * 54.78001 * 3 2 1 5 5 O 1.42901 * 109.46750 * 294.78670 * 3 2 1 6 6 C 1.42910 * 113.99820 * 151.08169 * 5 3 2 7 7 C 1.54266 * 110.02995 * 269.52739 * 6 5 3 8 8 C 1.54888 * 104.19590 * 217.10742 * 7 6 5 9 9 C 1.54889 * 102.74960 * 37.94058 * 8 7 6 10 10 C 1.53873 * 110.02957 * 152.36060 * 6 5 3 11 11 C 1.50696 * 109.47317 * 174.78185 * 3 2 1 12 12 O 1.21288 * 119.99549 * 120.00141 * 11 3 2 13 13 N 1.34770 * 120.00308 * 300.00285 * 11 3 2 14 14 C 1.46495 * 120.00040 * 180.02562 * 13 11 3 15 15 C 1.50699 * 109.47492 * 179.97438 * 14 13 11 16 16 C 1.29449 * 117.70055 * 94.99292 * 15 14 13 17 17 C 1.50321 * 123.65528 * 180.58271 * 16 15 14 18 18 N 1.47005 * 109.16453 * 349.28073 * 17 16 15 19 19 C 1.36940 * 121.37487 * 222.97107 * 18 17 16 20 20 H 1.07998 * 120.00041 * 159.55642 * 19 18 17 21 21 C 1.38807 * 119.99931 * 339.55799 * 19 18 17 22 22 N 1.31624 * 120.00366 * 343.49006 * 21 19 18 23 23 Si 1.86807 * 120.00023 * 163.49084 * 21 19 18 24 24 H 1.48500 * 109.46685 * 89.99821 * 23 21 19 25 25 H 1.48493 * 109.47016 * 209.99972 * 23 21 19 26 26 H 1.48498 * 109.46994 * 330.00425 * 23 21 19 27 27 C 1.47337 * 117.24770 * 42.94861 * 18 17 16 28 28 C 1.49966 * 117.68945 * 275.22239 * 15 14 13 29 29 H 1.09000 * 109.46931 * 65.29227 * 1 2 3 30 30 H 1.09000 * 109.47238 * 185.30058 * 1 2 3 31 31 H 1.09000 * 109.46998 * 305.30007 * 1 2 3 32 32 H 1.09000 * 109.47345 * 239.99614 * 2 1 3 33 33 H 1.09005 * 109.47002 * 119.99744 * 2 1 3 34 34 H 1.08996 * 110.01466 * 30.96934 * 6 5 3 35 35 H 1.09006 * 110.48783 * 335.78800 * 7 6 5 36 36 H 1.08998 * 110.49091 * 98.42680 * 7 6 5 37 37 H 1.08996 * 110.75771 * 279.62830 * 8 7 6 38 38 H 1.08995 * 110.91524 * 156.33158 * 8 7 6 39 39 H 1.08998 * 110.48963 * 203.38358 * 9 8 7 40 40 H 1.09005 * 110.59827 * 80.79413 * 9 8 7 41 41 H 1.08996 * 110.03163 * 238.56504 * 10 6 5 42 42 H 1.08999 * 110.03291 * 359.95150 * 10 6 5 43 43 H 0.96997 * 120.00019 * 359.97438 * 13 11 3 44 44 H 1.08998 * 109.47691 * 59.99728 * 14 13 11 45 45 H 1.09000 * 109.47426 * 299.99781 * 14 13 11 46 46 H 1.08002 * 118.17125 * 0.86474 * 16 15 14 47 47 H 1.08998 * 109.52550 * 109.19478 * 17 16 15 48 48 H 1.09005 * 109.51888 * 229.37984 * 17 16 15 49 49 H 0.96998 * 119.99838 * 180.02562 * 22 21 19 50 50 H 1.08999 * 109.77740 * 57.81650 * 27 18 17 51 51 H 1.09000 * 109.77626 * 178.61643 * 27 18 17 52 52 H 1.09003 * 109.43865 * 40.81564 * 28 15 14 53 53 H 1.09003 * 109.12738 * 281.42566 * 28 15 14 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 1 1.6024 1.9825 0.8396 5 8 1.6664 2.0798 -1.2232 6 6 1.4837 3.4927 -1.1103 7 6 0.0150 3.8093 -0.7603 8 6 -0.2945 5.1024 -1.5548 9 6 0.4647 4.8758 -2.8857 10 6 1.7697 4.1707 -2.4616 11 6 3.5415 1.4444 0.1292 12 8 4.2243 1.9426 -0.7406 13 7 4.1255 0.8926 1.2112 14 6 5.5850 0.8950 1.3370 15 6 5.9792 0.2232 2.6272 16 6 6.2809 -1.0348 2.5829 17 6 6.6794 -1.8389 3.7888 18 7 6.9138 -0.9292 4.9195 19 6 7.9861 -1.0923 5.7554 20 1 8.3338 -0.2645 6.3556 21 6 8.6263 -2.3216 5.8312 22 7 8.0329 -3.4023 5.3702 23 14 10.3301 -2.4423 6.5876 24 1 10.2075 -2.6957 8.0457 25 1 11.0758 -3.5583 5.9522 26 1 11.0613 -1.1689 6.3662 27 6 5.9475 0.1656 5.1160 28 6 5.9893 1.0582 3.8728 29 1 -0.3633 0.4296 -0.9336 30 1 -0.3634 -1.0233 0.0949 31 1 -0.3633 0.5938 0.8387 32 1 1.8934 -0.5139 -0.8899 33 1 1.8933 -0.5138 0.8901 34 1 2.1460 3.8943 -0.3434 35 1 -0.6375 2.9979 -1.0830 36 1 -0.0930 3.9841 0.3101 37 1 0.0925 5.9779 -1.0336 38 1 -1.3651 5.2027 -1.7330 39 1 0.6850 5.8290 -3.3663 40 1 -0.1141 4.2377 -3.5535 41 1 2.5701 4.9021 -2.3504 42 1 2.0486 3.4223 -3.2033 43 1 3.5794 0.4945 1.9070 44 1 5.9481 1.9227 1.3373 45 1 6.0230 0.3552 0.4975 46 1 6.2486 -1.5335 1.6254 47 1 5.8808 -2.5363 4.0418 48 1 7.5927 -2.3929 3.5716 49 1 8.4802 -4.2613 5.4235 50 1 4.9457 -0.2448 5.2419 51 1 6.2223 0.7464 5.9965 52 1 5.1221 1.7186 3.8739 53 1 6.8976 1.6605 3.8918 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 ZINC001755894372.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 20:37:59 Heat of formation + Delta-G solvation = -64.755441 kcal Electronic energy + Delta-G solvation = -29636.236408 eV Core-core repulsion = 25700.922195 eV Total energy + Delta-G solvation = -3935.314214 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.21 seconds Orbital eigenvalues (eV) -40.30075 -38.80929 -37.76159 -36.34763 -34.97957 -33.66285 -32.71858 -30.93246 -29.75945 -29.00886 -28.47422 -27.58208 -26.18626 -24.04878 -22.77510 -22.30654 -21.20317 -21.05951 -20.59251 -20.34927 -17.92287 -17.54829 -17.39294 -16.16865 -15.88984 -15.72023 -15.22005 -15.06832 -14.92712 -14.44004 -14.24203 -14.09290 -13.89097 -13.62812 -13.57367 -13.45001 -13.03071 -12.76055 -12.51597 -12.28563 -12.16897 -12.06798 -11.99638 -11.60405 -11.45788 -11.39972 -11.39187 -11.31776 -11.14072 -10.94604 -10.92996 -10.75162 -10.74434 -10.46793 -10.39522 -10.33151 -9.85881 -9.66053 -9.31942 -8.97741 -8.72502 -8.06698 -6.83083 -5.82974 -3.19136 2.52965 2.59711 3.27157 3.35943 3.43165 3.45526 3.60969 4.26594 4.30539 4.46551 4.56257 4.59926 4.66518 4.75222 4.94618 5.03674 5.11547 5.16791 5.19383 5.27532 5.32315 5.33941 5.36174 5.37369 5.48998 5.54493 5.55728 5.59856 5.69595 5.72839 5.78306 5.79490 5.89872 5.92147 5.96289 6.01732 6.02714 6.02757 6.12642 6.28630 6.42550 6.50237 6.68156 6.81322 6.95819 7.13580 7.24439 7.29886 7.56817 7.68736 7.80143 8.23492 8.32777 9.15645 9.78870 10.34696 11.27210 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.025592 B = 0.002042 C = 0.001996 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1093.811378 B =13705.649966 C =14026.216418 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.148 4.148 2 C -0.117 4.117 3 C 0.076 3.924 4 H 0.085 0.915 5 O -0.346 6.346 6 C 0.086 3.914 7 C -0.149 4.149 8 C -0.123 4.123 9 C -0.123 4.123 10 C -0.119 4.119 11 C 0.484 3.516 12 O -0.531 6.531 13 N -0.719 5.719 14 C 0.157 3.843 15 C -0.168 4.168 16 C -0.123 4.123 17 C 0.191 3.809 18 N -0.596 5.596 19 C -0.243 4.243 20 H 0.070 0.930 21 C -0.005 4.005 22 N -0.911 5.911 23 Si 0.905 3.095 24 H -0.288 1.288 25 H -0.289 1.289 26 H -0.280 1.280 27 C 0.144 3.856 28 C -0.090 4.090 29 H 0.062 0.938 30 H 0.059 0.941 31 H 0.054 0.946 32 H 0.078 0.922 33 H 0.082 0.918 34 H 0.062 0.938 35 H 0.078 0.922 36 H 0.065 0.935 37 H 0.067 0.933 38 H 0.068 0.932 39 H 0.066 0.934 40 H 0.067 0.933 41 H 0.071 0.929 42 H 0.079 0.921 43 H 0.400 0.600 44 H 0.064 0.936 45 H 0.065 0.935 46 H 0.094 0.906 47 H 0.040 0.960 48 H 0.101 0.899 49 H 0.254 0.746 50 H 0.023 0.977 51 H 0.057 0.943 52 H 0.052 0.948 53 H 0.064 0.936 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -18.322 12.443 -12.499 25.431 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.155 4.155 3 C 0.014 3.986 4 H 0.102 0.898 5 O -0.263 6.263 6 C 0.029 3.971 7 C -0.187 4.187 8 C -0.161 4.161 9 C -0.161 4.161 10 C -0.157 4.157 11 C 0.269 3.731 12 O -0.408 6.408 13 N -0.372 5.372 14 C 0.035 3.965 15 C -0.169 4.169 16 C -0.144 4.144 17 C 0.066 3.934 18 N -0.319 5.319 19 C -0.372 4.372 20 H 0.088 0.912 21 C -0.202 4.202 22 N -0.557 5.557 23 Si 0.735 3.265 24 H -0.215 1.215 25 H -0.215 1.215 26 H -0.206 1.206 27 C 0.018 3.982 28 C -0.129 4.129 29 H 0.081 0.919 30 H 0.078 0.922 31 H 0.073 0.927 32 H 0.096 0.904 33 H 0.100 0.900 34 H 0.080 0.920 35 H 0.097 0.903 36 H 0.084 0.916 37 H 0.086 0.914 38 H 0.086 0.914 39 H 0.084 0.916 40 H 0.086 0.914 41 H 0.090 0.910 42 H 0.098 0.902 43 H 0.235 0.765 44 H 0.083 0.917 45 H 0.083 0.917 46 H 0.112 0.888 47 H 0.058 0.942 48 H 0.119 0.881 49 H 0.064 0.936 50 H 0.041 0.959 51 H 0.075 0.925 52 H 0.071 0.929 53 H 0.082 0.918 Dipole moment (debyes) X Y Z Total from point charges -17.750 12.552 -13.307 25.489 hybrid contribution 0.812 -0.218 0.922 1.248 sum -16.939 12.334 -12.385 24.340 Atomic orbital electron populations 1.21721 0.94431 1.01984 1.02351 1.21601 0.97303 0.92850 1.03746 1.22118 0.90800 0.96334 0.89344 0.89763 1.88046 1.87860 1.20077 1.30317 1.22834 0.95029 0.81093 0.98181 1.22954 0.97057 0.98559 1.00102 1.22423 1.00242 0.97329 0.96098 1.22353 0.96188 1.00280 0.97285 1.22151 0.97668 0.99388 0.96463 1.21435 0.91225 0.78461 0.81978 1.90701 1.65099 1.44947 1.40013 1.46265 1.08339 1.59898 1.22696 1.20445 0.79204 1.00658 0.96195 1.20479 1.06429 0.95588 0.94354 1.23298 0.95108 0.95407 1.00577 1.19854 1.00932 0.90178 0.82455 1.44166 1.26894 1.31011 1.29864 1.19202 1.04377 0.94816 1.18837 0.91203 1.24915 0.99575 0.95018 1.00689 1.70005 1.36626 0.97880 1.51193 0.86248 0.84376 0.77264 0.78640 1.21475 1.21546 1.20580 1.20937 0.92718 0.90112 0.94475 1.20749 1.00314 0.97469 0.94342 0.91898 0.92205 0.92729 0.90369 0.89968 0.92020 0.90304 0.91610 0.91440 0.91362 0.91551 0.91417 0.91046 0.90190 0.76518 0.91748 0.91693 0.88777 0.94238 0.88146 0.93590 0.95925 0.92462 0.92944 0.91793 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 19. 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 -1.09 9.76 37.16 0.36 -0.73 16 2 C -0.12 -1.04 5.59 -26.73 -0.15 -1.19 16 3 C 0.08 0.79 2.35 -27.89 -0.07 0.73 16 4 H 0.08 0.72 7.83 -51.93 -0.41 0.32 16 5 O -0.35 -4.06 8.73 -55.14 -0.48 -4.54 16 6 C 0.09 0.87 3.08 -25.61 -0.08 0.79 16 7 C -0.15 -1.15 6.54 -25.07 -0.16 -1.31 16 8 C -0.12 -0.83 6.91 -24.74 -0.17 -1.00 16 9 C -0.12 -0.91 6.83 -25.07 -0.17 -1.08 16 10 C -0.12 -1.14 6.61 -25.61 -0.17 -1.31 16 11 C 0.48 6.60 7.14 -10.99 -0.08 6.52 16 12 O -0.53 -9.38 16.44 -14.36 -0.24 -9.62 16 13 N -0.72 -9.15 5.41 -61.48 -0.33 -9.48 16 14 C 0.16 2.35 5.66 -5.20 -0.03 2.32 16 15 C -0.17 -2.81 5.23 -100.99 -0.53 -3.34 16 16 C -0.12 -2.29 9.47 -35.78 -0.34 -2.63 16 17 C 0.19 4.41 5.28 -5.10 -0.03 4.39 16 18 N -0.60 -14.56 3.02 -173.68 -0.52 -15.08 16 19 C -0.24 -6.55 9.74 -15.06 -0.15 -6.70 16 20 H 0.07 1.84 7.88 -52.49 -0.41 1.43 16 21 C 0.00 -0.13 5.02 -17.43 -0.09 -0.22 16 22 N -0.91 -25.82 10.88 55.49 0.60 -25.22 16 23 Si 0.90 21.52 29.59 -169.99 -5.03 16.49 16 24 H -0.29 -6.83 7.11 56.52 0.40 -6.43 16 25 H -0.29 -6.94 7.11 56.52 0.40 -6.53 16 26 H -0.28 -6.58 7.11 56.52 0.40 -6.17 16 27 C 0.14 2.87 6.58 -3.54 -0.02 2.85 16 28 C -0.09 -1.51 5.94 -28.19 -0.17 -1.68 16 29 H 0.06 0.50 7.97 -51.93 -0.41 0.08 16 30 H 0.06 0.40 8.14 -51.93 -0.42 -0.02 16 31 H 0.05 0.35 8.14 -51.93 -0.42 -0.07 16 32 H 0.08 0.79 8.14 -51.93 -0.42 0.36 16 33 H 0.08 0.65 7.42 -51.93 -0.39 0.27 16 34 H 0.06 0.67 7.83 -51.93 -0.41 0.26 16 35 H 0.08 0.62 8.08 -51.93 -0.42 0.20 16 36 H 0.07 0.45 8.14 -51.93 -0.42 0.03 16 37 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 38 H 0.07 0.39 8.14 -51.93 -0.42 -0.03 16 39 H 0.07 0.42 8.14 -51.93 -0.42 0.00 16 40 H 0.07 0.49 8.14 -51.93 -0.42 0.06 16 41 H 0.07 0.68 8.14 -51.93 -0.42 0.25 16 42 H 0.08 0.87 7.95 -51.93 -0.41 0.46 16 43 H 0.40 4.34 7.82 -40.82 -0.32 4.02 16 44 H 0.06 0.99 8.11 -51.93 -0.42 0.57 16 45 H 0.06 1.03 7.78 -51.93 -0.40 0.62 16 46 H 0.09 1.59 7.73 -52.49 -0.41 1.18 16 47 H 0.04 0.97 8.14 -51.93 -0.42 0.55 16 48 H 0.10 2.62 6.16 -51.93 -0.32 2.30 16 49 H 0.25 6.98 8.31 -40.82 -0.34 6.65 16 50 H 0.02 0.44 8.14 -51.93 -0.42 0.02 16 51 H 0.06 1.11 7.96 -51.93 -0.41 0.70 16 52 H 0.05 0.75 8.14 -51.93 -0.42 0.33 16 53 H 0.06 1.13 8.14 -51.93 -0.42 0.71 16 LS Contribution 417.80 15.07 6.30 6.30 Total: -1.00 -31.12 417.80 -11.51 -42.63 By element: Atomic # 1 Polarization: 11.88 SS G_CDS: -9.77 Total: 2.10 kcal Atomic # 6 Polarization: -1.55 SS G_CDS: -2.03 Total: -3.59 kcal Atomic # 7 Polarization: -49.52 SS G_CDS: -0.25 Total: -49.78 kcal Atomic # 8 Polarization: -13.44 SS G_CDS: -0.72 Total: -14.16 kcal Atomic # 14 Polarization: 21.52 SS G_CDS: -5.03 Total: 16.49 kcal Total LS contribution 6.30 Total: 6.30 kcal Total: -31.12 -11.51 -42.63 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. ZINC001755894372.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 -22.122 kcal (2) G-P(sol) polarization free energy of solvation -31.123 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -53.245 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.510 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.634 kcal (6) G-S(sol) free energy of system = (1) + (5) -64.755 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.21 seconds