Wall clock time and date at job start Sat Apr 11 2020 06:17:27 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.53009 * 1 3 3 C 1.53000 * 109.47081 * 2 1 4 4 C 1.52997 * 109.47075 * 120.00458 * 2 1 3 5 5 O 1.42900 * 109.46867 * 174.99712 * 4 2 1 6 6 C 1.42906 * 113.99748 * 179.97438 * 5 4 2 7 7 C 1.52999 * 109.46720 * 179.97438 * 6 5 4 8 8 C 1.52998 * 109.46994 * 180.02562 * 7 6 5 9 9 N 1.46502 * 109.47360 * 180.02562 * 8 7 6 10 10 C 1.34782 * 120.00166 * 180.02562 * 9 8 7 11 11 O 1.21350 * 119.99528 * 359.97438 * 10 9 8 12 12 C 1.49446 * 120.00084 * 179.97438 * 10 9 8 13 13 O 1.21356 * 119.99803 * 353.37756 * 12 10 9 14 14 N 1.34768 * 120.00057 * 173.37736 * 12 10 9 15 15 C 1.46923 * 120.63140 * 353.06161 * 14 12 10 16 16 C 1.53197 * 108.77603 * 126.41596 * 15 14 12 17 17 C 1.53036 * 109.34037 * 54.62030 * 16 15 14 18 18 C 1.53043 * 109.53781 * 298.68543 * 17 16 15 19 19 H 1.09000 * 109.49772 * 301.34791 * 18 17 16 20 20 C 1.50699 * 109.49861 * 181.28187 * 18 17 16 21 21 H 1.08009 * 119.99826 * 0.06421 * 20 18 17 22 22 C 1.37874 * 120.00391 * 180.34183 * 20 18 17 23 23 N 1.31511 * 119.99935 * 359.72505 * 22 20 18 24 24 Si 1.86803 * 120.00371 * 179.72468 * 22 20 18 25 25 H 1.48501 * 109.46889 * 59.99781 * 24 22 20 26 26 H 1.48493 * 109.47092 * 179.97438 * 24 22 20 27 27 H 1.48501 * 109.46907 * 300.00065 * 24 22 20 28 28 C 1.46927 * 120.63303 * 173.08900 * 14 12 10 29 29 H 1.09005 * 109.46752 * 299.99900 * 1 2 3 30 30 H 1.08995 * 109.46893 * 59.99303 * 1 2 3 31 31 H 1.09000 * 109.47088 * 180.02562 * 1 2 3 32 32 H 1.09005 * 109.46752 * 240.00100 * 2 1 3 33 33 H 1.08996 * 109.47213 * 180.02562 * 3 2 1 34 34 H 1.08993 * 109.47145 * 300.00418 * 3 2 1 35 35 H 1.08998 * 109.46895 * 60.00371 * 3 2 1 36 36 H 1.09002 * 109.47378 * 54.99717 * 4 2 1 37 37 H 1.09001 * 109.47250 * 294.99793 * 4 2 1 38 38 H 1.08999 * 109.47003 * 60.00559 * 6 5 4 39 39 H 1.08993 * 109.47283 * 299.99924 * 6 5 4 40 40 H 1.09004 * 109.47307 * 60.00310 * 7 6 5 41 41 H 1.08999 * 109.47615 * 300.00205 * 7 6 5 42 42 H 1.09007 * 109.47352 * 60.00069 * 8 7 6 43 43 H 1.09010 * 109.47442 * 300.00155 * 8 7 6 44 44 H 0.96996 * 120.00059 * 359.97438 * 9 8 7 45 45 H 1.08998 * 109.59171 * 246.19792 * 15 14 12 46 46 H 1.09003 * 109.58922 * 6.48999 * 15 14 12 47 47 H 1.08995 * 109.50150 * 174.58993 * 16 15 14 48 48 H 1.09009 * 109.49617 * 294.68452 * 16 15 14 49 49 H 1.08997 * 109.49188 * 178.64856 * 17 16 15 50 50 H 1.09001 * 109.43006 * 58.68768 * 17 16 15 51 51 H 0.97000 * 120.00590 * 180.02562 * 23 22 20 52 52 H 1.09003 * 109.59192 * 113.79627 * 28 14 12 53 53 H 1.08998 * 109.70190 * 353.44355 * 28 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.5301 0.0000 0.0000 3 6 2.0401 1.4425 0.0000 4 6 2.0401 -0.7213 1.2492 5 8 3.4642 -0.8239 1.1918 6 6 4.0451 -1.4848 2.3179 7 6 5.5650 -1.5376 2.1500 8 6 6.1869 -2.2457 3.3552 9 7 7.6422 -2.2957 3.1949 10 6 8.4049 -2.8781 4.1413 11 8 7.8861 -3.3607 5.1264 12 6 9.8895 -2.9296 3.9774 13 8 10.3988 -2.5567 2.9409 14 7 10.6626 -3.3902 4.9806 15 6 10.0793 -3.7022 6.2925 16 6 10.4605 -5.1363 6.6731 17 6 11.9818 -5.2904 6.6090 18 6 12.4634 -5.0247 5.1809 19 1 11.9795 -5.7229 4.4979 20 6 13.9577 -5.2074 5.1117 21 1 14.5092 -5.4863 5.9976 22 6 14.6216 -5.0129 3.9191 23 7 13.9502 -4.6732 2.8406 24 14 16.4749 -5.2319 3.8353 25 1 16.8268 -6.6225 4.2192 26 1 16.9423 -4.9677 2.4508 27 1 17.1285 -4.2800 4.7691 28 6 12.1049 -3.5893 4.7834 29 1 -0.3633 0.5139 0.8901 30 1 -0.3633 0.5139 -0.8899 31 1 -0.3633 -1.0277 0.0005 32 1 1.8934 -0.5139 -0.8901 33 1 3.1300 1.4425 -0.0005 34 1 1.6767 1.9563 -0.8899 35 1 1.6768 1.9563 0.8900 36 1 1.6054 -1.7199 1.2951 37 1 1.7517 -0.1585 2.1369 38 1 3.6521 -2.4992 2.3863 39 1 3.7985 -0.9374 3.2275 40 1 5.9584 -0.5233 2.0813 41 1 5.8113 -2.0852 1.2403 42 1 5.7935 -3.2600 3.4236 43 1 5.9406 -1.6984 4.2652 44 1 8.0569 -1.9103 2.4072 45 1 10.4685 -3.0103 7.0394 46 1 8.9943 -3.6118 6.2408 47 1 10.1158 -5.3480 7.6852 48 1 9.9950 -5.8336 5.9764 49 1 12.2575 -6.3032 6.9026 50 1 12.4475 -4.5762 7.2881 51 1 14.4172 -4.5368 2.0014 52 1 12.6593 -2.8881 5.4072 53 1 12.3586 -3.4260 3.7360 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 ZINC001347416877.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:17:27 Heat of formation + Delta-G solvation = -127.016012 kcal Electronic energy + Delta-G solvation = -29515.728395 eV Core-core repulsion = 25387.258309 eV Total energy + Delta-G solvation = -4128.470086 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 2.52 seconds Orbital eigenvalues (eV) -40.51529 -38.85378 -38.30511 -36.87383 -36.12144 -35.13320 -33.27077 -31.98851 -31.66182 -29.73213 -27.58652 -27.27374 -26.05560 -25.09080 -24.74688 -22.58075 -21.89049 -21.77783 -20.64582 -20.08637 -18.76156 -17.36835 -17.07477 -16.83187 -16.56304 -16.10039 -15.85051 -15.62295 -15.09818 -14.98773 -14.61304 -14.55049 -14.35560 -13.80943 -13.66855 -13.43016 -13.27218 -13.18548 -12.85960 -12.68281 -12.58113 -12.50328 -12.22705 -12.08079 -11.80122 -11.67882 -11.55752 -11.44913 -11.33971 -11.19148 -11.00433 -10.96033 -10.71115 -10.46629 -10.26961 -10.05597 -9.97639 -9.88210 -9.30062 -9.21358 -9.10815 -8.52908 -8.41077 -7.82138 -6.01529 -4.09837 1.78331 2.59636 3.26290 3.32245 3.37571 3.39210 3.64359 3.97711 4.38913 4.42070 4.44046 4.56307 4.59008 4.64862 4.78105 4.86178 5.00675 5.03013 5.09734 5.13237 5.17047 5.26135 5.28945 5.35203 5.43526 5.50864 5.63009 5.64745 5.69991 5.72306 5.74471 5.85047 5.89103 5.93004 6.07040 6.26424 6.39967 6.60355 6.62440 6.75339 6.82089 6.87335 7.00101 7.05130 7.18468 7.23207 7.59565 7.83388 7.95718 7.98839 8.16821 8.35665 8.57895 8.93980 9.55678 11.04083 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.022247 B = 0.001703 C = 0.001610 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1258.303631 B =16433.568539 C =17386.560699 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.147 4.147 2 C -0.096 4.096 3 C -0.145 4.145 4 C 0.066 3.934 5 O -0.386 6.386 6 C 0.068 3.932 7 C -0.133 4.133 8 C 0.121 3.879 9 N -0.704 5.704 10 C 0.513 3.487 11 O -0.515 6.515 12 C 0.495 3.505 13 O -0.524 6.524 14 N -0.567 5.567 15 C 0.050 3.950 16 C -0.128 4.128 17 C -0.113 4.113 18 C 0.041 3.959 19 H 0.029 0.971 20 C -0.520 4.520 21 H 0.059 0.941 22 C 0.062 3.938 23 N -0.883 5.883 24 Si 0.894 3.106 25 H -0.277 1.277 26 H -0.287 1.287 27 H -0.277 1.277 28 C 0.108 3.892 29 H 0.056 0.944 30 H 0.061 0.939 31 H 0.056 0.944 32 H 0.078 0.922 33 H 0.063 0.937 34 H 0.053 0.947 35 H 0.052 0.948 36 H 0.048 0.952 37 H 0.049 0.951 38 H 0.053 0.947 39 H 0.052 0.948 40 H 0.081 0.919 41 H 0.082 0.918 42 H 0.074 0.926 43 H 0.074 0.926 44 H 0.415 0.585 45 H 0.054 0.946 46 H 0.143 0.857 47 H 0.059 0.941 48 H 0.058 0.942 49 H 0.050 0.950 50 H 0.049 0.951 51 H 0.264 0.736 52 H 0.048 0.952 53 H 0.124 0.876 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -25.331 9.148 -1.534 26.976 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.204 4.204 2 C -0.115 4.115 3 C -0.202 4.202 4 C -0.010 4.010 5 O -0.305 6.305 6 C -0.010 4.010 7 C -0.170 4.170 8 C -0.003 4.003 9 N -0.357 5.357 10 C 0.295 3.705 11 O -0.392 6.392 12 C 0.275 3.725 13 O -0.402 6.402 14 N -0.298 5.298 15 C -0.072 4.072 16 C -0.166 4.166 17 C -0.151 4.151 18 C 0.021 3.979 19 H 0.047 0.953 20 C -0.558 4.558 21 H 0.077 0.923 22 C -0.140 4.140 23 N -0.521 5.521 24 Si 0.722 3.278 25 H -0.203 1.203 26 H -0.213 1.213 27 H -0.202 1.202 28 C -0.012 4.012 29 H 0.075 0.925 30 H 0.080 0.920 31 H 0.075 0.925 32 H 0.097 0.903 33 H 0.082 0.918 34 H 0.072 0.928 35 H 0.071 0.929 36 H 0.067 0.933 37 H 0.067 0.933 38 H 0.071 0.929 39 H 0.071 0.929 40 H 0.099 0.901 41 H 0.100 0.900 42 H 0.092 0.908 43 H 0.092 0.908 44 H 0.251 0.749 45 H 0.073 0.927 46 H 0.160 0.840 47 H 0.077 0.923 48 H 0.077 0.923 49 H 0.069 0.931 50 H 0.067 0.933 51 H 0.074 0.926 52 H 0.066 0.934 53 H 0.141 0.859 Dipole moment (debyes) X Y Z Total from point charges -25.275 9.121 -1.727 26.926 hybrid contribution 1.131 -0.121 -0.135 1.146 sum -24.144 8.999 -1.861 25.833 Atomic orbital electron populations 1.21795 0.94307 1.02018 1.02288 1.21110 0.95463 0.97020 0.97901 1.21790 1.02031 0.94998 1.01417 1.22585 0.81060 0.99225 0.98150 1.87858 1.20072 1.79142 1.43477 1.22304 0.93119 0.96622 0.88909 1.21684 0.93939 1.02565 0.98860 1.21259 0.79399 1.01472 0.98174 1.45778 1.06985 1.58573 1.24413 1.20793 0.89114 0.77879 0.82664 1.90302 1.74307 1.46023 1.28603 1.19732 0.86578 0.81008 0.85139 1.90762 1.74263 1.47431 1.27720 1.48438 1.10693 1.59782 1.10917 1.23078 1.04294 0.96312 0.83492 1.21856 0.98206 0.96309 1.00200 1.21665 0.92941 1.00399 1.00080 1.18972 0.95955 0.91969 0.90959 0.95263 1.21061 0.91676 1.47674 0.95420 0.92296 1.25056 0.98450 0.94818 0.95708 1.69856 1.35605 1.46409 1.00196 0.86621 0.85478 0.78640 0.77058 1.20260 1.21252 1.20221 1.22491 0.73700 0.99101 1.05886 0.92500 0.91953 0.92490 0.90311 0.91752 0.92777 0.92927 0.93335 0.93309 0.92896 0.92932 0.90054 0.89970 0.90775 0.90820 0.74876 0.92743 0.84028 0.92283 0.92345 0.93073 0.93254 0.92559 0.93381 0.85860 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.15 -0.67 9.19 37.16 0.34 -0.33 16 2 C -0.10 -0.59 3.16 -90.62 -0.29 -0.88 16 3 C -0.15 -0.91 9.09 37.16 0.34 -0.57 16 4 C 0.07 0.45 5.32 37.16 0.20 0.64 16 5 O -0.39 -3.38 9.88 -35.23 -0.35 -3.73 16 6 C 0.07 0.49 5.94 37.16 0.22 0.71 16 7 C -0.13 -1.03 5.70 -26.73 -0.15 -1.18 16 8 C 0.12 1.14 5.67 -4.04 -0.02 1.11 16 9 N -0.70 -8.88 5.55 -60.72 -0.34 -9.22 16 10 C 0.51 8.41 7.78 -11.58 -0.09 8.32 16 11 O -0.51 -8.73 13.10 5.50 0.07 -8.66 16 12 C 0.49 10.03 7.74 -11.59 -0.09 9.94 16 13 O -0.52 -11.97 15.95 5.49 0.09 -11.88 16 14 N -0.57 -11.94 2.97 -173.27 -0.51 -12.45 16 15 C 0.05 0.90 6.27 -3.71 -0.02 0.87 16 16 C -0.13 -2.29 6.08 -26.61 -0.16 -2.46 16 17 C -0.11 -2.35 5.11 -26.69 -0.14 -2.49 16 18 C 0.04 1.02 2.25 -91.65 -0.21 0.81 16 19 H 0.03 0.82 8.14 -51.93 -0.42 0.39 16 20 C -0.52 -14.15 8.57 -34.44 -0.30 -14.44 16 21 H 0.06 1.57 7.74 -52.48 -0.41 1.16 16 22 C 0.06 1.72 5.30 -17.82 -0.09 1.63 16 23 N -0.88 -25.55 11.38 55.43 0.63 -24.92 16 24 Si 0.89 21.23 29.54 -169.99 -5.02 16.21 16 25 H -0.28 -6.50 7.11 56.52 0.40 -6.10 16 26 H -0.29 -6.80 7.11 56.52 0.40 -6.40 16 27 H -0.28 -6.47 7.11 56.52 0.40 -6.07 16 28 C 0.11 2.65 5.38 -3.71 -0.02 2.63 16 29 H 0.06 0.24 8.14 -51.93 -0.42 -0.18 16 30 H 0.06 0.26 8.14 -51.93 -0.42 -0.16 16 31 H 0.06 0.24 8.14 -51.93 -0.42 -0.18 16 32 H 0.08 0.53 8.14 -51.93 -0.42 0.11 16 33 H 0.06 0.49 7.82 -51.93 -0.41 0.08 16 34 H 0.05 0.30 8.14 -51.93 -0.42 -0.12 16 35 H 0.05 0.31 8.14 -51.93 -0.42 -0.11 16 36 H 0.05 0.29 8.14 -51.93 -0.42 -0.13 16 37 H 0.05 0.29 8.14 -51.93 -0.42 -0.13 16 38 H 0.05 0.34 8.14 -51.93 -0.42 -0.08 16 39 H 0.05 0.33 8.14 -51.93 -0.42 -0.09 16 40 H 0.08 0.63 8.14 -51.93 -0.42 0.21 16 41 H 0.08 0.67 8.14 -51.93 -0.42 0.25 16 42 H 0.07 0.71 8.14 -51.93 -0.42 0.28 16 43 H 0.07 0.67 8.14 -51.92 -0.42 0.25 16 44 H 0.41 4.96 7.96 -40.82 -0.33 4.63 16 45 H 0.05 0.92 8.14 -51.93 -0.42 0.49 16 46 H 0.14 2.46 3.78 -56.75 -0.21 2.25 16 47 H 0.06 0.90 8.14 -51.93 -0.42 0.47 16 48 H 0.06 1.10 8.14 -51.92 -0.42 0.67 16 49 H 0.05 1.04 8.14 -51.93 -0.42 0.62 16 50 H 0.05 1.00 8.14 -51.93 -0.42 0.58 16 51 H 0.26 7.36 8.31 -40.82 -0.34 7.02 16 52 H 0.05 1.17 8.14 -51.93 -0.42 0.75 16 53 H 0.12 3.45 5.17 -51.93 -0.27 3.18 16 LS Contribution 420.00 15.07 6.33 6.33 Total: -1.00 -31.16 420.00 -9.22 -40.38 By element: Atomic # 1 Polarization: 13.27 SS G_CDS: -9.63 Total: 3.63 kcal Atomic # 6 Polarization: 4.79 SS G_CDS: -0.48 Total: 4.31 kcal Atomic # 7 Polarization: -46.37 SS G_CDS: -0.22 Total: -46.59 kcal Atomic # 8 Polarization: -24.08 SS G_CDS: -0.19 Total: -24.27 kcal Atomic # 14 Polarization: 21.23 SS G_CDS: -5.02 Total: 16.21 kcal Total LS contribution 6.33 Total: 6.33 kcal Total: -31.16 -9.22 -40.38 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. ZINC001347416877.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 -86.641 kcal (2) G-P(sol) polarization free energy of solvation -31.160 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -117.800 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.216 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.375 kcal (6) G-S(sol) free energy of system = (1) + (5) -127.016 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 2.52 seconds