Wall clock time and date at job start Wed Apr 1 2020 00:24: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.53125 * 1 3 3 H 1.08865 * 109.97722 * 2 1 4 4 O 1.43668 * 109.31539 * 121.17532 * 2 1 3 5 5 C 1.43788 * 110.17575 * 114.26931 * 4 2 1 6 6 C 1.55009 * 105.52708 * 340.83335 * 5 4 2 7 7 C 1.53664 * 110.20039 * 238.48553 * 2 1 3 8 8 C 1.52314 * 110.62874 * 152.35323 * 7 2 1 9 9 N 1.47001 * 110.19037 * 283.02131 * 8 7 2 10 10 C 1.46915 * 111.05298 * 178.94046 * 9 8 7 11 11 C 1.52997 * 109.47044 * 188.75679 * 10 9 8 12 12 C 1.50697 * 109.47466 * 179.97438 * 11 10 9 13 13 O 1.21288 * 120.00038 * 0.02562 * 12 11 10 14 14 N 1.34777 * 120.00369 * 180.02562 * 12 11 10 15 15 C 1.46504 * 119.99736 * 359.97438 * 14 12 11 16 16 C 1.46505 * 120.00115 * 180.02562 * 14 12 11 17 17 C 1.50702 * 109.46836 * 89.99960 * 16 14 12 18 18 H 1.07994 * 120.00445 * 0.02562 * 17 16 14 19 19 C 1.37877 * 119.99368 * 179.97438 * 17 16 14 20 20 N 1.31515 * 120.00175 * 0.02562 * 19 17 16 21 21 Si 1.86804 * 119.99773 * 179.97438 * 19 17 16 22 22 H 1.48500 * 109.47075 * 0.02562 * 21 19 17 23 23 H 1.48507 * 109.46791 * 120.00371 * 21 19 17 24 24 H 1.48497 * 109.47353 * 240.00199 * 21 19 17 25 25 C 1.47171 * 110.05991 * 302.13546 * 9 8 7 26 26 C 1.53145 * 108.56767 * 62.23409 * 25 9 8 27 27 H 1.09085 * 109.10161 * 63.40730 * 26 25 9 28 28 C 1.53192 * 109.24070 * 182.51726 * 26 25 9 29 29 O 1.41267 * 112.03432 * 25.71056 * 7 2 1 30 30 H 1.09001 * 109.47155 * 181.50953 * 1 2 3 31 31 H 1.08995 * 109.47335 * 301.51253 * 1 2 3 32 32 H 1.09005 * 109.46863 * 61.51612 * 1 2 3 33 33 H 1.09326 * 110.60783 * 100.08700 * 5 4 2 34 34 H 1.09326 * 110.60318 * 221.61095 * 5 4 2 35 35 H 1.08700 * 111.40315 * 277.17410 * 6 5 4 36 36 H 1.08792 * 111.26554 * 152.07108 * 6 5 4 37 37 H 1.09145 * 109.43071 * 162.72789 * 8 7 2 38 38 H 1.09147 * 109.45791 * 43.27335 * 8 7 2 39 39 H 1.09002 * 109.47236 * 308.75739 * 10 9 8 40 40 H 1.09001 * 109.47372 * 68.75929 * 10 9 8 41 41 H 1.09002 * 109.47153 * 300.00111 * 11 10 9 42 42 H 1.08998 * 109.47533 * 59.99740 * 11 10 9 43 43 H 1.08996 * 109.46784 * 90.00263 * 15 14 12 44 44 H 1.09003 * 109.47146 * 209.99886 * 15 14 12 45 45 H 1.08994 * 109.47096 * 330.00470 * 15 14 12 46 46 H 1.08998 * 109.46420 * 209.99986 * 16 14 12 47 47 H 1.08995 * 109.46904 * 329.99992 * 16 14 12 48 48 H 0.97004 * 119.99914 * 180.02562 * 20 19 17 49 49 H 1.08852 * 109.58387 * 181.97313 * 25 9 8 50 50 H 1.08884 * 109.54646 * 302.57948 * 25 9 8 51 51 H 1.09001 * 109.46982 * 60.00511 * 28 26 25 52 52 H 1.09003 * 109.46867 * 180.02562 * 28 26 25 53 53 H 1.09002 * 109.47341 * 300.00771 * 28 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.5313 0.0000 0.0000 3 1 1.9032 1.0231 0.0000 4 8 2.0065 -0.7018 1.1600 5 6 2.6940 -1.9018 0.7665 6 6 2.2213 -2.1955 -0.6802 7 6 2.0619 -0.7538 -1.2294 8 6 3.4322 -0.2253 -1.6328 9 7 3.2939 1.0569 -2.3384 10 6 4.6075 1.5967 -2.7145 11 6 4.4401 3.0215 -3.2462 12 6 5.7876 3.5756 -3.6314 13 8 6.7838 2.8974 -3.4946 14 7 5.8857 4.8247 -4.1279 15 6 4.6824 5.6437 -4.2937 16 6 7.1957 5.3631 -4.5029 17 6 7.4777 5.0378 -5.9472 18 1 6.7535 4.4857 -6.5277 19 6 8.6598 5.4460 -6.5275 20 7 9.5415 6.1188 -5.8206 21 14 9.0098 5.0423 -8.3175 22 1 7.8676 4.2825 -8.8860 23 1 9.1954 6.3019 -9.0818 24 1 10.2444 4.2221 -8.4077 25 6 2.4357 0.8914 -3.5225 26 6 1.0416 0.4601 -3.0580 27 1 0.6059 1.2548 -2.4509 28 6 0.1526 0.2223 -4.2826 29 8 1.1376 -0.7425 -2.2977 30 1 -0.3633 -1.0273 0.0271 31 1 -0.3634 0.5371 0.8761 32 1 -0.3633 0.4901 -0.9033 33 1 3.7764 -1.7497 0.7879 34 1 2.4285 -2.7315 1.4271 35 1 1.2692 -2.7199 -0.6877 36 1 2.9798 -2.7379 -1.2407 37 1 3.9217 -0.9480 -2.2880 38 1 4.0439 -0.0821 -0.7402 39 1 5.2589 1.6092 -1.8407 40 1 5.0500 0.9697 -3.4886 41 1 3.7887 3.0091 -4.1201 42 1 3.9976 3.6486 -2.4722 43 1 4.5040 6.2157 -3.3832 44 1 4.8208 6.3276 -5.1311 45 1 3.8271 4.9972 -4.4901 46 1 7.1962 6.4446 -4.3672 47 1 7.9655 4.9173 -3.8731 48 1 10.3731 6.4063 -6.2290 49 1 2.3670 1.8365 -4.0582 50 1 2.8566 0.1268 -4.1735 51 1 0.5896 -0.5607 -4.9024 52 1 -0.8411 -0.0846 -3.9561 53 1 0.0765 1.1432 -4.8608 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 ZINC000363939292.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Wed Apr 1 2020 00:24:59 Heat of formation + Delta-G solvation = -116.508082 kcal Electronic energy + Delta-G solvation = -31448.770100 eV Core-core repulsion = 27320.755672 eV Total energy + Delta-G solvation = -4128.014428 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.97 seconds Orbital eigenvalues (eV) -41.21457 -39.15917 -38.39977 -36.33298 -34.81074 -33.58961 -32.68443 -31.68289 -31.38304 -30.45736 -27.98553 -27.46721 -26.33152 -25.64768 -24.43670 -22.99602 -22.36912 -21.98124 -21.12823 -19.87983 -18.40635 -17.23158 -16.80528 -16.67498 -16.56481 -15.96682 -15.88022 -15.65339 -15.31424 -14.90066 -14.60739 -14.51796 -14.12731 -14.09211 -13.64902 -13.46226 -13.31089 -13.06074 -12.75461 -12.65636 -12.48358 -12.30405 -12.27389 -12.17894 -12.05106 -12.03743 -11.82168 -11.68569 -11.51624 -11.35940 -11.18616 -11.15937 -10.92254 -10.77312 -10.59879 -10.39062 -10.18303 -9.79259 -9.76653 -9.63217 -8.92932 -8.60428 -8.46121 -7.65958 -6.06554 -4.11403 3.04534 3.27385 3.34189 3.36921 3.42126 3.43272 3.64183 4.00722 4.11320 4.31761 4.42962 4.43796 4.61599 4.68809 4.75343 4.85465 4.97785 5.02123 5.09966 5.12432 5.16260 5.21275 5.27022 5.28690 5.32248 5.34186 5.37908 5.43794 5.45812 5.55370 5.56776 5.59113 5.65560 5.70484 5.78946 5.85906 5.88158 5.90734 6.00955 6.02069 6.04920 6.19587 6.68150 6.74290 6.78054 7.04184 7.07421 7.42760 7.52266 7.89857 8.21290 8.40098 8.88017 9.07376 9.53786 11.01569 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.018390 B = 0.002353 C = 0.002285 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1522.198658 B =11896.740885 C =12250.324762 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.153 4.153 2 C 0.062 3.938 3 H 0.096 0.904 4 O -0.375 6.375 5 C 0.050 3.950 6 C -0.145 4.145 7 C 0.105 3.895 8 C 0.032 3.968 9 N -0.533 5.533 10 C 0.072 3.928 11 C -0.134 4.134 12 C 0.504 3.496 13 O -0.546 6.546 14 N -0.603 5.603 15 C 0.084 3.916 16 C 0.244 3.756 17 C -0.525 4.525 18 H 0.057 0.943 19 C 0.060 3.940 20 N -0.894 5.894 21 Si 0.898 3.102 22 H -0.274 1.274 23 H -0.283 1.283 24 H -0.283 1.283 25 C 0.013 3.987 26 C 0.092 3.908 27 H 0.069 0.931 28 C -0.145 4.145 29 O -0.360 6.360 30 H 0.058 0.942 31 H 0.058 0.942 32 H 0.068 0.932 33 H 0.066 0.934 34 H 0.085 0.915 35 H 0.090 0.910 36 H 0.089 0.911 37 H 0.047 0.953 38 H 0.091 0.909 39 H 0.085 0.915 40 H 0.046 0.954 41 H 0.085 0.915 42 H 0.094 0.906 43 H 0.045 0.955 44 H 0.073 0.927 45 H 0.056 0.944 46 H 0.032 0.968 47 H 0.044 0.956 48 H 0.263 0.737 49 H 0.092 0.908 50 H 0.046 0.954 51 H 0.063 0.937 52 H 0.065 0.935 53 H 0.075 0.925 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -21.305 -13.107 9.453 26.741 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.210 4.210 2 C 0.004 3.996 3 H 0.114 0.886 4 O -0.295 6.295 5 C -0.026 4.026 6 C -0.183 4.183 7 C 0.064 3.936 8 C -0.097 4.097 9 N -0.258 5.258 10 C -0.056 4.056 11 C -0.175 4.175 12 C 0.290 3.710 13 O -0.424 6.424 14 N -0.335 5.335 15 C -0.060 4.060 16 C 0.125 3.875 17 C -0.564 4.564 18 H 0.075 0.925 19 C -0.142 4.142 20 N -0.534 5.534 21 Si 0.726 3.274 22 H -0.198 1.198 23 H -0.209 1.209 24 H -0.209 1.209 25 C -0.115 4.115 26 C 0.035 3.965 27 H 0.087 0.913 28 C -0.203 4.203 29 O -0.278 6.278 30 H 0.077 0.923 31 H 0.078 0.922 32 H 0.087 0.913 33 H 0.084 0.916 34 H 0.103 0.897 35 H 0.109 0.891 36 H 0.107 0.893 37 H 0.065 0.935 38 H 0.109 0.891 39 H 0.103 0.897 40 H 0.064 0.936 41 H 0.103 0.897 42 H 0.112 0.888 43 H 0.063 0.937 44 H 0.092 0.908 45 H 0.075 0.925 46 H 0.050 0.950 47 H 0.062 0.938 48 H 0.073 0.927 49 H 0.110 0.890 50 H 0.064 0.936 51 H 0.082 0.918 52 H 0.084 0.916 53 H 0.094 0.906 Dipole moment (debyes) X Y Z Total from point charges -20.938 -13.089 9.991 26.637 hybrid contribution 0.889 -0.449 -1.855 2.105 sum -20.048 -13.537 8.137 25.522 Atomic orbital electron populations 1.22027 0.94523 1.02518 1.01970 1.23158 0.95021 0.96982 0.84456 0.88625 1.88543 1.74621 1.36415 1.29882 1.23411 0.96916 0.88041 0.94240 1.23246 1.03507 0.95453 0.96063 1.22234 0.88539 0.95390 0.87429 1.22643 0.96230 0.90831 0.99974 1.62464 1.14557 1.24999 1.23763 1.21844 0.88350 0.94391 1.00982 1.21633 0.96716 0.96904 1.02239 1.20261 0.90302 0.82088 0.78303 1.90670 1.38780 1.57992 1.54938 1.48160 1.09866 1.13765 1.61681 1.21633 0.87400 0.95692 1.01279 1.19296 0.78365 0.94000 0.95859 1.21459 1.03299 1.35642 0.96010 0.92505 1.24947 0.95833 0.95060 0.98342 1.69945 1.08942 1.38119 1.36355 0.86445 0.78337 0.78040 0.84601 1.19849 1.20925 1.20861 1.22992 0.93862 1.02103 0.92558 1.22180 0.94624 0.86970 0.92733 0.91271 1.21817 0.99563 1.02816 0.96095 1.87310 1.59825 1.36736 1.43972 0.92330 0.92248 0.91324 0.91632 0.89688 0.89116 0.89274 0.93492 0.89098 0.89665 0.93555 0.89666 0.88778 0.93681 0.90829 0.92503 0.94999 0.93841 0.92749 0.88968 0.93556 0.91833 0.91559 0.90614 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 -1.39 8.95 37.22 0.33 -1.05 16 2 C 0.06 0.62 3.44 -26.32 -0.09 0.53 16 3 H 0.10 1.06 7.86 -52.01 -0.41 0.65 16 4 O -0.37 -4.10 11.01 -35.23 -0.39 -4.49 16 5 C 0.05 0.41 7.31 38.22 0.28 0.69 16 6 C -0.14 -1.16 6.27 -24.56 -0.15 -1.32 16 7 C 0.10 1.01 0.78 -89.52 -0.07 0.94 16 8 C 0.03 0.33 3.96 -4.12 -0.02 0.32 16 9 N -0.53 -6.04 4.70 -234.63 -1.10 -7.14 16 10 C 0.07 0.99 4.24 -3.82 -0.02 0.97 16 11 C -0.13 -1.91 4.19 -27.89 -0.12 -2.03 16 12 C 0.50 10.28 7.67 -10.99 -0.08 10.20 16 13 O -0.55 -13.56 15.56 5.55 0.09 -13.47 16 14 N -0.60 -12.69 2.86 -182.95 -0.52 -13.21 16 15 C 0.08 1.30 10.24 59.85 0.61 1.91 16 16 C 0.24 6.41 5.16 -5.19 -0.03 6.38 16 17 C -0.53 -14.92 9.39 -34.44 -0.32 -15.24 16 18 H 0.06 1.61 7.81 -52.49 -0.41 1.20 16 19 C 0.06 1.72 5.46 -17.82 -0.10 1.63 16 20 N -0.89 -26.67 13.29 55.43 0.74 -25.93 16 21 Si 0.90 21.85 29.54 -169.99 -5.02 16.83 16 22 H -0.27 -6.56 7.11 56.52 0.40 -6.15 16 23 H -0.28 -6.70 7.11 56.52 0.40 -6.30 16 24 H -0.28 -6.88 7.11 56.52 0.40 -6.48 16 25 C 0.01 0.13 4.73 -3.60 -0.02 0.11 16 26 C 0.09 0.82 2.93 -26.56 -0.08 0.75 16 27 H 0.07 0.61 6.96 -51.88 -0.36 0.25 16 28 C -0.15 -1.04 9.57 37.26 0.36 -0.68 16 29 O -0.36 -3.61 8.59 -35.03 -0.30 -3.91 16 30 H 0.06 0.53 8.14 -51.93 -0.42 0.11 16 31 H 0.06 0.52 8.14 -51.93 -0.42 0.09 16 32 H 0.07 0.59 5.95 -51.93 -0.31 0.28 16 33 H 0.07 0.50 7.84 -51.74 -0.41 0.09 16 34 H 0.08 0.56 8.17 -51.74 -0.42 0.14 16 35 H 0.09 0.71 8.12 -52.10 -0.42 0.29 16 36 H 0.09 0.63 7.84 -52.05 -0.41 0.23 16 37 H 0.05 0.51 7.96 -51.85 -0.41 0.10 16 38 H 0.09 0.96 7.40 -51.85 -0.38 0.58 16 39 H 0.08 1.26 8.03 -51.93 -0.42 0.84 16 40 H 0.05 0.71 8.08 -51.93 -0.42 0.29 16 41 H 0.09 1.07 5.71 -51.93 -0.30 0.78 16 42 H 0.09 1.19 8.11 -51.93 -0.42 0.77 16 43 H 0.04 0.60 8.14 -51.93 -0.42 0.18 16 44 H 0.07 1.18 8.07 -51.93 -0.42 0.76 16 45 H 0.06 0.70 6.44 -51.93 -0.33 0.37 16 46 H 0.03 0.84 8.07 -51.93 -0.42 0.43 16 47 H 0.04 1.29 7.42 -51.93 -0.39 0.90 16 48 H 0.26 7.36 8.31 -40.82 -0.34 7.02 16 49 H 0.09 0.87 6.24 -52.01 -0.32 0.55 16 50 H 0.05 0.51 8.13 -51.99 -0.42 0.09 16 51 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 52 H 0.07 0.44 8.14 -51.93 -0.42 0.02 16 53 H 0.08 0.45 8.14 -51.93 -0.42 0.03 16 LS Contribution 408.53 15.07 6.16 6.16 Total: -1.00 -33.57 408.53 -9.34 -42.91 By element: Atomic # 1 Polarization: 7.62 SS G_CDS: -9.48 Total: -1.86 kcal Atomic # 6 Polarization: 3.62 SS G_CDS: 0.49 Total: 4.11 kcal Atomic # 7 Polarization: -45.39 SS G_CDS: -0.89 Total: -46.28 kcal Atomic # 8 Polarization: -21.27 SS G_CDS: -0.60 Total: -21.87 kcal Atomic # 14 Polarization: 21.85 SS G_CDS: -5.02 Total: 16.83 kcal Total LS contribution 6.16 Total: 6.16 kcal Total: -33.57 -9.34 -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. ZINC000363939292.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 -73.595 kcal (2) G-P(sol) polarization free energy of solvation -33.574 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -107.169 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.339 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) -116.508 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.97 seconds