Wall clock time and date at job start Sun Mar 7 2021 03:56:20 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= 0 * 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 N 1.46507 * 1 3 3 C 1.34771 * 120.00173 * 2 1 4 4 O 1.21286 * 120.00396 * 180.02562 * 3 2 1 5 5 C 1.50701 * 119.99957 * 0.02562 * 3 2 1 6 6 N 1.46506 * 109.47436 * 179.97438 * 5 3 2 7 7 C 1.36832 * 125.06219 * 265.03384 * 6 5 3 8 8 C 1.34242 * 109.93768 * 179.97438 * 7 6 5 9 9 C 1.50699 * 126.49906 * 180.02562 * 8 7 6 10 10 O 1.42902 * 109.46646 * 104.97731 * 9 8 7 11 11 Si 1.86299 * 109.47458 * 344.97436 * 9 8 7 12 12 C 1.46490 * 106.99769 * 359.97438 * 8 7 6 13 13 C 1.39591 * 134.03949 * 180.02562 * 12 8 7 14 14 C 1.36662 * 119.76515 * 179.97438 * 13 12 8 15 15 C 1.38792 * 120.53033 * 0.02562 * 14 13 12 16 16 C 1.37787 * 120.66396 * 359.97438 * 15 14 13 17 17 C 1.37681 * 125.06064 * 84.72013 * 6 5 3 18 18 C 1.46501 * 119.99744 * 180.02562 * 2 1 3 19 19 H 1.09008 * 112.84595 * 245.18787 * 18 2 1 20 20 C 1.54050 * 113.61370 * 16.49901 * 18 2 1 21 21 C 1.53774 * 86.98479 * 89.11837 * 20 18 2 22 22 H 1.09003 * 113.61561 * 270.83846 * 21 20 18 23 23 C 1.50699 * 113.61494 * 140.01960 * 21 20 18 24 24 C 1.38242 * 119.99874 * 320.58649 * 23 21 20 25 25 C 1.38234 * 119.99705 * 179.78594 * 24 23 21 26 26 C 1.38237 * 120.00130 * 0.44619 * 25 24 23 27 27 C 1.38234 * 119.99819 * 359.80966 * 26 25 24 28 28 C 1.38226 * 120.00210 * 359.97438 * 27 26 25 29 29 C 1.53781 * 87.17688 * 25.42897 * 21 20 18 30 30 H 1.08997 * 109.46801 * 264.66129 * 1 2 3 31 31 H 1.09000 * 109.46898 * 24.65943 * 1 2 3 32 32 H 1.08995 * 109.47117 * 144.65551 * 1 2 3 33 33 H 1.09000 * 109.47661 * 300.00111 * 5 3 2 34 34 H 1.09003 * 109.47232 * 60.00173 * 5 3 2 35 35 H 1.07999 * 125.03281 * 359.95205 * 7 6 5 36 36 H 1.08996 * 109.46843 * 224.97269 * 9 8 7 37 37 H 0.96704 * 113.99694 * 179.97438 * 10 9 8 38 38 H 1.48496 * 109.99968 * 298.68057 * 11 9 8 39 39 H 1.48503 * 109.99861 * 60.00060 * 11 9 8 40 40 H 1.08000 * 120.11858 * 0.02562 * 13 12 8 41 41 H 1.08003 * 119.73439 * 180.02562 * 14 13 12 42 42 H 1.07995 * 119.66451 * 179.97438 * 15 14 13 43 43 H 1.08000 * 120.11239 * 180.02562 * 16 15 14 44 44 H 1.08998 * 113.69395 * 334.59810 * 20 18 2 45 45 H 1.08997 * 113.61075 * 203.62440 * 20 18 2 46 46 H 1.07993 * 119.99967 * 0.02562 * 24 23 21 47 47 H 1.08003 * 119.99586 * 180.27455 * 25 24 23 48 48 H 1.08000 * 119.99592 * 179.82889 * 26 25 24 49 49 H 1.07999 * 119.99637 * 180.02562 * 27 26 25 50 50 H 1.08002 * 120.00080 * 180.02562 * 28 27 26 51 51 H 1.08999 * 113.61422 * 89.06843 * 29 21 20 52 52 H 1.08997 * 113.61358 * 219.89754 * 29 21 20 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 7 1.4651 0.0000 0.0000 3 6 2.1390 1.1671 0.0000 4 8 3.3518 1.1672 -0.0005 5 6 1.3855 2.4723 0.0006 6 7 2.3375 3.5858 0.0011 7 6 2.7746 4.2458 1.1172 8 6 3.6401 5.2110 0.7687 9 6 4.3468 6.1607 1.7013 10 8 5.7138 5.7667 1.8362 11 14 3.5213 6.1210 3.3709 12 6 3.7684 5.1650 -0.6898 13 6 4.5034 5.8982 -1.6229 14 6 4.4014 5.5999 -2.9527 15 6 3.5722 4.5736 -3.3834 16 6 2.8379 3.8372 -2.4796 17 6 2.9302 4.1201 -1.1209 18 6 2.1975 -1.2688 -0.0006 19 1 2.7742 -1.4249 0.9112 20 6 1.3339 -2.4796 -0.4020 21 6 1.7382 -2.1832 -1.8557 22 1 1.0677 -1.4892 -2.3625 23 6 2.0225 -3.4206 -2.6676 24 6 2.6938 -4.4843 -2.0942 25 6 2.9506 -5.6211 -2.8376 26 6 2.5442 -5.6905 -4.1571 27 6 1.8776 -4.6247 -4.7319 28 6 1.6169 -3.4897 -3.9873 29 6 3.0017 -1.5153 -1.2880 30 1 -0.3633 -0.0956 1.0232 31 1 -0.3633 0.9340 -0.4288 32 1 -0.3633 -0.8382 -0.5945 33 1 0.7584 2.5303 -0.8891 34 1 0.7592 2.5298 0.8909 35 1 2.4683 4.0235 2.1287 36 1 4.2983 7.1713 1.2959 37 1 6.2308 6.3351 2.4233 38 1 2.1078 6.5592 3.2492 39 1 3.5809 4.7483 3.9344 40 1 5.1509 6.6980 -1.2951 41 1 4.9706 6.1673 -3.6741 42 1 3.5020 4.3500 -4.4376 43 1 2.1958 3.0403 -2.8245 44 1 0.2702 -2.3425 -0.2079 45 1 1.7143 -3.4302 -0.0282 46 1 3.0142 -4.4288 -1.0644 47 1 3.4712 -6.4539 -2.3883 48 1 2.7474 -6.5776 -4.7387 49 1 1.5598 -4.6791 -5.7627 50 1 1.0957 -2.6571 -4.4363 51 1 3.8372 -2.2027 -1.1559 52 1 3.2842 -0.6009 -1.8097 There are 69 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 69 No. of singly occupied orbitals= 1 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=0 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 s_22_14797354_9818708.mol2 52 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sun Mar 7 2021 03:56:20 Heat of formation + Delta-G solvation = 76.234267 kcal Electronic energy + Delta-G solvation = -34191.638562 eV Core-core repulsion = 29896.893905 eV Total energy + Delta-G solvation = -4294.744657 eV No. of doubly occupied orbitals = 69 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 377.178 amu Computer time = 1.43 seconds Orbital eigenvalues (eV) -42.54356 -41.78893 -40.27263 -38.79668 -37.91082 -37.44965 -35.79170 -34.07838 -32.08287 -31.78728 -31.23884 -30.17174 -29.03685 -28.67204 -26.71987 -26.05159 -25.35239 -23.80080 -23.54256 -22.88686 -22.25065 -21.18543 -20.41396 -19.55534 -18.55298 -18.23282 -18.11170 -18.04229 -17.67358 -17.29814 -16.76428 -16.35188 -16.19635 -16.10126 -15.55252 -15.43755 -15.28827 -14.94511 -14.79794 -14.64826 -14.58975 -14.32790 -14.30644 -14.18628 -14.08745 -14.01741 -13.88051 -13.48491 -13.46957 -13.10592 -12.84751 -12.80042 -12.62227 -12.43389 -12.34866 -12.27408 -11.89734 -11.80951 -11.63955 -11.62126 -11.47118 -10.96290 -10.40822 -10.08805 -9.84542 -9.81984 -9.57905 -9.01930 -8.52864 -4.96386 0.27242 0.40433 0.47846 0.79706 1.32921 1.37397 1.39792 1.66028 1.76599 1.91648 2.03268 2.79899 3.04031 3.15073 3.34968 3.45060 3.54427 3.59485 3.66673 3.73431 3.76203 3.85114 3.85974 3.86444 3.94537 4.08033 4.10511 4.12539 4.15104 4.21814 4.32228 4.36144 4.37360 4.46294 4.48176 4.55401 4.62026 4.63244 4.70428 4.70983 4.78627 4.86753 4.90840 4.97893 4.99420 5.07253 5.07862 5.11975 5.25625 5.43499 5.55435 5.61663 5.64842 5.69807 5.80545 5.86041 5.89697 6.04886 6.07628 6.25866 6.36315 6.84313 6.97335 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.064 3.936 2 N -0.600 5.600 3 C 0.492 3.508 4 O -0.504 6.504 5 C 0.094 3.906 6 N -0.466 5.466 7 C 0.065 3.935 8 C -0.151 4.151 9 C -0.032 4.032 10 O -0.545 6.545 11 Si 0.682 3.318 12 C -0.085 4.085 13 C -0.063 4.063 14 C -0.145 4.145 15 C -0.096 4.096 16 C -0.127 4.127 17 C 0.110 3.890 18 C 0.123 3.877 19 H 0.087 0.913 20 C -0.137 4.137 21 C -0.070 4.070 22 H 0.095 0.905 23 C -0.072 4.072 24 C -0.117 4.117 25 C -0.117 4.117 26 C -0.124 4.124 27 C -0.117 4.117 28 C -0.116 4.116 29 C -0.129 4.129 30 H 0.064 0.936 31 H 0.074 0.926 32 H 0.083 0.917 33 H 0.118 0.882 34 H 0.120 0.880 35 H 0.162 0.838 36 H 0.070 0.930 37 H 0.377 0.623 38 H -0.256 1.256 39 H -0.254 1.254 40 H 0.123 0.877 41 H 0.124 0.876 42 H 0.126 0.874 43 H 0.114 0.886 44 H 0.083 0.917 45 H 0.101 0.899 46 H 0.121 0.879 47 H 0.122 0.878 48 H 0.122 0.878 49 H 0.122 0.878 50 H 0.119 0.881 51 H 0.094 0.906 52 H 0.076 0.924 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -5.003 0.205 1.576 5.250 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.079 4.079 2 N -0.336 5.336 3 C 0.277 3.723 4 O -0.377 6.377 5 C -0.032 4.032 6 N -0.149 5.149 7 C -0.060 4.060 8 C -0.158 4.158 9 C -0.144 4.144 10 O -0.351 6.351 11 Si 0.583 3.417 12 C -0.089 4.089 13 C -0.082 4.082 14 C -0.165 4.165 15 C -0.115 4.115 16 C -0.148 4.148 17 C 0.006 3.994 18 C 0.020 3.980 19 H 0.105 0.895 20 C -0.175 4.175 21 C -0.089 4.089 22 H 0.113 0.887 23 C -0.072 4.072 24 C -0.135 4.135 25 C -0.135 4.135 26 C -0.142 4.142 27 C -0.135 4.135 28 C -0.134 4.134 29 C -0.167 4.167 30 H 0.083 0.917 31 H 0.093 0.907 32 H 0.101 0.899 33 H 0.136 0.864 34 H 0.138 0.862 35 H 0.179 0.821 36 H 0.088 0.912 37 H 0.224 0.776 38 H -0.181 1.181 39 H -0.179 1.179 40 H 0.141 0.859 41 H 0.142 0.858 42 H 0.143 0.857 43 H 0.132 0.868 44 H 0.101 0.899 45 H 0.119 0.881 46 H 0.139 0.861 47 H 0.140 0.860 48 H 0.140 0.860 49 H 0.140 0.860 50 H 0.137 0.863 51 H 0.112 0.888 52 H 0.095 0.905 Dipole moment (debyes) X Y Z Total from point charges -5.419 -0.854 1.687 5.740 hybrid contribution -0.215 -0.391 -0.587 0.738 sum -5.634 -1.246 1.100 5.874 Atomic orbital electron populations 1.22021 0.79184 1.03825 1.02851 1.47593 1.06869 1.05012 1.74145 1.20921 0.88664 0.85159 0.77508 1.90830 1.13257 1.86668 1.46988 1.21653 0.91075 0.83989 1.06522 1.43202 1.38740 1.27939 1.05044 1.21395 0.99535 0.94420 0.90665 1.19529 1.03032 1.01266 0.91993 1.23934 0.87007 0.97281 1.06152 1.86579 1.20384 1.57059 1.71031 0.95301 0.81587 0.84688 0.80102 1.18925 0.99864 0.97862 0.92280 1.20281 0.97347 0.98569 0.91959 1.20669 1.00950 0.99283 0.95550 1.20823 0.96201 0.95485 0.98964 1.19939 1.03512 1.02141 0.89182 1.18164 0.97170 0.95943 0.88157 1.21948 0.94260 0.82737 0.99045 0.89485 1.23529 1.00793 0.97771 0.95364 1.21556 0.98550 0.94629 0.94136 0.88672 1.19689 0.98241 0.95079 0.94238 1.20996 0.99146 0.93362 0.99999 1.21180 0.99801 0.98007 0.94511 1.21204 0.98712 0.98837 0.95415 1.21183 0.98762 0.93650 0.99946 1.21057 1.00421 0.98160 0.93809 1.23298 0.96414 1.01284 0.95751 0.91719 0.90744 0.89851 0.86422 0.86238 0.82136 0.91243 0.77576 1.18059 1.17904 0.85874 0.85805 0.85654 0.86786 0.89895 0.88053 0.86068 0.85964 0.86030 0.86003 0.86266 0.88765 0.90502 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.06 -0.27 9.74 59.85 0.58 0.31 16 2 N -0.60 0.36 3.01 -172.83 -0.52 -0.17 16 3 C 0.49 1.05 7.46 -10.99 -0.08 0.97 16 4 O -0.50 -3.59 15.60 5.55 0.09 -3.51 16 5 C 0.09 -0.02 4.96 -5.19 -0.03 -0.05 16 6 N -0.47 -1.68 2.27 -115.70 -0.26 -1.94 16 7 C 0.06 0.27 9.89 -16.73 -0.17 0.10 16 8 C -0.15 -0.75 5.61 -102.71 -0.58 -1.32 16 9 C -0.03 -0.06 3.05 36.01 0.11 0.04 16 10 O -0.54 -1.46 12.86 -35.23 -0.45 -1.91 16 11 Si 0.68 -0.08 30.79 -169.99 -5.23 -5.31 16 12 C -0.09 -0.50 6.35 -105.14 -0.67 -1.16 16 13 C -0.06 -0.28 10.05 -39.65 -0.40 -0.68 16 14 C -0.15 -0.47 9.97 -39.95 -0.40 -0.87 16 15 C -0.10 -0.28 10.01 -39.54 -0.40 -0.68 16 16 C -0.13 -0.47 10.00 -39.43 -0.39 -0.86 16 17 C 0.11 0.56 6.79 -83.99 -0.57 -0.01 16 18 C 0.12 0.00 4.10 -66.97 -0.27 -0.27 16 19 H 0.09 0.09 8.14 -51.92 -0.42 -0.33 16 20 C -0.14 0.26 6.59 -25.78 -0.17 0.09 16 21 C -0.07 0.04 3.75 -90.96 -0.34 -0.30 16 22 H 0.10 -0.09 7.98 -51.93 -0.41 -0.51 16 23 C -0.07 -0.05 5.29 -104.62 -0.55 -0.61 16 24 C -0.12 -0.08 8.53 -39.58 -0.34 -0.41 16 25 C -0.12 -0.09 10.05 -39.58 -0.40 -0.49 16 26 C -0.12 -0.11 10.05 -39.58 -0.40 -0.51 16 27 C -0.12 -0.11 10.05 -39.59 -0.40 -0.50 16 28 C -0.12 -0.11 9.71 -39.59 -0.38 -0.49 16 29 C -0.13 -0.14 6.89 -25.92 -0.18 -0.31 16 30 H 0.06 -0.32 8.14 -51.93 -0.42 -0.74 16 31 H 0.07 -0.43 6.25 -51.93 -0.32 -0.76 16 32 H 0.08 -0.40 5.63 -51.93 -0.29 -0.69 16 33 H 0.12 -0.32 7.42 -51.93 -0.39 -0.71 16 34 H 0.12 -0.37 8.08 -51.93 -0.42 -0.79 16 35 H 0.16 0.37 6.30 -52.49 -0.33 0.03 16 36 H 0.07 0.03 7.96 -51.93 -0.41 -0.38 16 37 H 0.38 -1.08 8.74 45.56 0.40 -0.68 16 38 H -0.26 -0.50 7.11 56.52 0.40 -0.10 16 39 H -0.25 -0.72 7.04 56.52 0.40 -0.32 16 40 H 0.12 0.36 8.06 -52.49 -0.42 -0.06 16 41 H 0.12 0.12 8.06 -52.48 -0.42 -0.30 16 42 H 0.13 0.07 8.06 -52.49 -0.42 -0.35 16 43 H 0.11 0.19 8.06 -52.49 -0.42 -0.23 16 44 H 0.08 -0.26 5.82 -51.93 -0.30 -0.56 16 45 H 0.10 -0.22 6.41 -51.93 -0.33 -0.56 16 46 H 0.12 -0.10 6.73 -52.49 -0.35 -0.45 16 47 H 0.12 -0.10 8.06 -52.48 -0.42 -0.52 16 48 H 0.12 -0.09 8.06 -52.49 -0.42 -0.51 16 49 H 0.12 -0.09 8.06 -52.49 -0.42 -0.51 16 50 H 0.12 -0.06 8.05 -52.49 -0.42 -0.48 16 51 H 0.09 0.11 7.62 -51.93 -0.40 -0.28 16 52 H 0.08 0.27 7.42 -51.93 -0.39 -0.12 16 LS Contribution 420.66 15.07 6.34 6.34 Total: 0.00 -11.58 420.66 -13.84 -25.42 By element: Atomic # 1 Polarization: -3.54 SS G_CDS: -7.38 Total: -10.92 kcal Atomic # 6 Polarization: -1.59 SS G_CDS: -6.41 Total: -8.00 kcal Atomic # 7 Polarization: -1.33 SS G_CDS: -0.78 Total: -2.11 kcal Atomic # 8 Polarization: -5.05 SS G_CDS: -0.37 Total: -5.41 kcal Atomic # 14 Polarization: -0.08 SS G_CDS: -5.23 Total: -5.31 kcal Total LS contribution 6.34 Total: 6.34 kcal Total: -11.58 -13.84 -25.42 kcal The number of atoms in the molecule is 52 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. s_22_14797354_9818708.mol2 52 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 101.652 kcal (2) G-P(sol) polarization free energy of solvation -11.579 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 90.073 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -13.839 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -25.418 kcal (6) G-S(sol) free energy of system = (1) + (5) 76.234 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.43 seconds