Wall clock time and date at job start Tue Mar 31 2020 21:32:33 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.52998 * 1 3 3 H 1.09000 * 109.50792 * 2 1 4 4 C 1.53280 * 109.55077 * 120.12349 * 2 1 3 5 5 N 1.47170 * 108.52065 * 173.03068 * 4 2 1 6 6 C 1.34785 * 121.07701 * 128.67256 * 5 4 2 7 7 O 1.21281 * 120.00205 * 180.02562 * 6 5 4 8 8 C 1.50701 * 120.00003 * 0.02562 * 6 5 4 9 9 C 1.53004 * 109.47066 * 180.02562 * 8 6 5 10 10 C 1.52995 * 109.46760 * 179.97438 * 9 8 6 11 11 C 1.53000 * 115.55464 * 34.41727 * 10 9 8 12 12 C 1.53000 * 117.49774 * 248.74558 * 10 9 8 13 13 C 1.52998 * 117.49553 * 180.08765 * 10 9 8 14 14 C 1.47320 * 117.84856 * 308.69516 * 5 4 2 15 15 C 1.52362 * 108.46571 * 49.13081 * 14 5 4 16 16 C 1.54363 * 110.41525 * 73.70362 * 15 14 5 17 17 C 1.54469 * 102.87855 * 152.35347 * 16 15 14 18 18 N 1.47086 * 107.29487 * 338.51408 * 17 16 15 19 19 C 1.36940 * 125.62726 * 178.95584 * 18 17 16 20 20 H 1.07998 * 119.99851 * 180.02562 * 19 18 17 21 21 C 1.38811 * 119.99974 * 0.03021 * 19 18 17 22 22 N 1.31616 * 120.00374 * 359.97438 * 21 19 18 23 23 Si 1.86798 * 119.99863 * 179.97438 * 21 19 18 24 24 H 1.48505 * 109.47184 * 59.99576 * 23 21 19 25 25 H 1.48498 * 109.47156 * 179.97438 * 23 21 19 26 26 H 1.48501 * 109.47135 * 299.99535 * 23 21 19 27 27 C 1.47561 * 108.73903 * 358.99559 * 18 17 16 28 28 O 1.41872 * 111.48389 * 310.13891 * 15 14 5 29 29 H 1.08998 * 109.47110 * 180.09879 * 1 2 3 30 30 H 1.09003 * 109.47425 * 300.10302 * 1 2 3 31 31 H 1.09004 * 109.47178 * 60.10000 * 1 2 3 32 32 H 1.09007 * 109.63103 * 292.73905 * 4 2 1 33 33 H 1.08997 * 109.63351 * 53.12570 * 4 2 1 34 34 H 1.09000 * 109.47148 * 300.00075 * 8 6 5 35 35 H 1.09000 * 109.47168 * 59.99859 * 8 6 5 36 36 H 1.08995 * 109.46848 * 299.99618 * 9 8 6 37 37 H 1.09000 * 109.46671 * 60.00237 * 9 8 6 38 38 H 1.08998 * 109.47098 * 85.57564 * 11 10 9 39 39 H 1.08999 * 109.47429 * 205.58156 * 11 10 9 40 40 H 1.09000 * 109.46807 * 325.58470 * 11 10 9 41 41 H 1.08997 * 117.50357 * 0.02562 * 12 10 9 42 42 H 1.09005 * 117.49920 * 145.02795 * 12 10 9 43 43 H 1.09004 * 117.49902 * 214.97835 * 13 10 9 44 44 H 1.08998 * 117.50349 * 359.97438 * 13 10 9 45 45 H 1.08995 * 109.63370 * 289.44345 * 14 5 4 46 46 H 1.09000 * 109.63800 * 168.82477 * 14 5 4 47 47 H 1.09007 * 110.72978 * 34.00869 * 16 15 14 48 48 H 1.08995 * 110.73164 * 270.70183 * 16 15 14 49 49 H 1.08999 * 109.87633 * 219.09195 * 17 16 15 50 50 H 1.09000 * 109.87602 * 97.93611 * 17 16 15 51 51 H 0.97001 * 120.00283 * 180.02562 * 22 21 19 52 52 H 1.09003 * 110.37908 * 142.15389 * 27 18 17 53 53 H 1.09008 * 110.37665 * 264.38340 * 27 18 17 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 1 1.8940 1.0274 0.0000 4 6 2.0429 -0.7249 1.2494 5 7 3.5047 -0.8598 1.1445 6 6 4.3133 -0.4691 2.1495 7 8 5.5154 -0.5888 2.0421 8 6 3.7238 0.1166 3.4067 9 6 4.8495 0.4746 4.3791 10 6 4.2509 1.0697 5.6552 11 6 2.9036 0.4869 6.0863 12 6 4.4963 2.5557 5.9242 13 6 5.2214 1.4966 6.7582 14 6 4.0555 -1.4353 -0.0948 15 6 3.4214 -0.7212 -1.2819 16 6 4.0036 0.7009 -1.4286 17 6 3.8558 0.9792 -2.9408 18 7 3.7467 -0.3228 -3.6164 19 6 3.5816 -0.5115 -4.9626 20 1 3.5115 -1.5121 -5.3630 21 6 3.5043 0.5832 -5.8126 22 7 3.5893 1.8026 -5.3247 23 14 3.2798 0.3256 -7.6491 24 1 2.0130 -0.4094 -7.8946 25 1 3.2269 1.6437 -8.3309 26 1 4.4204 -0.4622 -8.1816 27 6 3.8436 -1.3919 -2.6039 28 8 2.0081 -0.6765 -1.1670 29 1 -0.3633 -1.0276 0.0018 30 1 -0.3634 0.5154 0.8891 31 1 -0.3634 0.5123 -0.8909 32 1 1.7915 -0.1457 2.1379 33 1 1.5850 -1.7119 1.3140 34 1 3.1588 1.0152 3.1587 35 1 3.0604 -0.6131 3.8710 36 1 5.4147 -0.4238 4.6272 37 1 5.5125 1.2044 3.9145 38 1 3.0699 -0.4091 6.6844 39 1 2.3617 1.2239 6.6790 40 1 2.3192 0.2304 5.2027 41 1 5.0838 3.1095 5.1920 42 1 3.6981 3.1181 6.4088 43 1 4.8999 1.3626 7.7911 44 1 6.2863 1.3531 6.5751 45 1 3.8221 -2.4990 -0.1399 46 1 5.1360 -1.2930 -0.1161 47 1 5.0527 0.7206 -1.1333 48 1 3.4241 1.4170 -0.8460 49 1 4.7323 1.5142 -3.3065 50 1 2.9568 1.5678 -3.1231 51 1 3.5350 2.5675 -5.9186 52 1 3.1649 -2.2090 -2.8485 53 1 4.8680 -1.7580 -2.5354 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 ZINC001131692542.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 21:32:33 Heat of formation + Delta-G solvation = -25.882027 kcal Electronic energy + Delta-G solvation = -30892.550125 eV Core-core repulsion = 26958.921589 eV Total energy + Delta-G solvation = -3933.628536 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.15 seconds Orbital eigenvalues (eV) -40.65898 -39.16299 -38.22884 -36.35699 -34.59521 -33.22166 -32.45602 -31.32043 -29.37964 -28.72882 -28.38172 -26.49561 -26.15092 -24.55615 -23.26144 -22.68550 -21.88830 -20.19297 -19.70092 -19.08285 -18.77306 -17.30009 -17.00607 -16.54206 -15.99083 -15.82197 -15.51860 -14.79071 -14.66106 -14.51342 -14.23927 -14.06343 -13.84481 -13.43286 -13.29376 -13.19625 -12.80708 -12.63298 -12.54277 -12.36713 -12.27474 -12.03113 -11.95888 -11.79077 -11.77719 -11.59255 -11.41983 -11.23982 -10.97710 -10.89770 -10.83002 -10.39794 -10.34072 -10.20643 -10.09241 -10.04571 -9.87351 -9.74451 -9.58986 -9.49260 -8.54233 -8.35555 -6.78448 -5.65695 -2.98849 2.59979 2.86015 3.44529 3.44931 3.48185 3.49639 3.86489 4.22810 4.40285 4.50931 4.53682 4.63138 4.68474 4.82313 4.87098 4.92340 5.04867 5.05314 5.14951 5.19162 5.26939 5.36221 5.38051 5.43976 5.46640 5.49475 5.56022 5.59878 5.60808 5.62236 5.72020 5.77763 5.81755 5.83014 5.85921 5.92779 5.94421 6.00511 6.09771 6.16287 6.18682 6.37045 6.50564 6.53084 6.56876 6.67174 7.34865 7.58979 7.66783 7.70126 7.87337 8.13658 8.29963 9.26554 9.86593 10.52441 11.46169 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.020664 B = 0.002600 C = 0.002525 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1354.715400 B =10768.100973 C =11088.250560 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.141 4.141 2 C 0.068 3.932 3 H 0.065 0.935 4 C 0.079 3.921 5 N -0.612 5.612 6 C 0.514 3.486 7 O -0.541 6.541 8 C -0.138 4.138 9 C -0.065 4.065 10 C -0.129 4.129 11 C -0.122 4.122 12 C -0.176 4.176 13 C -0.176 4.176 14 C 0.097 3.903 15 C 0.093 3.907 16 C -0.165 4.165 17 C 0.186 3.814 18 N -0.599 5.599 19 C -0.225 4.225 20 H 0.078 0.922 21 C -0.009 4.009 22 N -0.923 5.923 23 Si 0.913 3.087 24 H -0.285 1.285 25 H -0.291 1.291 26 H -0.285 1.285 27 C 0.153 3.847 28 O -0.344 6.344 29 H 0.063 0.937 30 H 0.068 0.932 31 H 0.071 0.929 32 H 0.089 0.911 33 H 0.077 0.923 34 H 0.094 0.906 35 H 0.096 0.904 36 H 0.073 0.927 37 H 0.073 0.927 38 H 0.059 0.941 39 H 0.056 0.944 40 H 0.066 0.934 41 H 0.094 0.906 42 H 0.092 0.908 43 H 0.092 0.908 44 H 0.094 0.906 45 H 0.071 0.929 46 H 0.101 0.899 47 H 0.064 0.936 48 H 0.060 0.940 49 H 0.042 0.958 50 H 0.046 0.954 51 H 0.252 0.748 52 H 0.050 0.950 53 H 0.037 0.963 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -4.254 -2.212 24.135 24.607 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.199 4.199 2 C 0.009 3.991 3 H 0.083 0.917 4 C -0.045 4.045 5 N -0.346 5.346 6 C 0.303 3.697 7 O -0.418 6.418 8 C -0.178 4.178 9 C -0.103 4.103 10 C -0.130 4.130 11 C -0.179 4.179 12 C -0.213 4.213 13 C -0.213 4.213 14 C -0.026 4.026 15 C 0.052 3.948 16 C -0.203 4.203 17 C 0.062 3.938 18 N -0.325 5.325 19 C -0.353 4.353 20 H 0.096 0.904 21 C -0.205 4.205 22 N -0.571 5.571 23 Si 0.743 3.257 24 H -0.212 1.212 25 H -0.217 1.217 26 H -0.212 1.212 27 C 0.026 3.974 28 O -0.262 6.262 29 H 0.082 0.918 30 H 0.087 0.913 31 H 0.090 0.910 32 H 0.108 0.892 33 H 0.095 0.905 34 H 0.112 0.888 35 H 0.114 0.886 36 H 0.092 0.908 37 H 0.091 0.909 38 H 0.078 0.922 39 H 0.075 0.925 40 H 0.085 0.915 41 H 0.113 0.887 42 H 0.111 0.889 43 H 0.111 0.889 44 H 0.112 0.888 45 H 0.089 0.911 46 H 0.119 0.881 47 H 0.083 0.917 48 H 0.079 0.921 49 H 0.060 0.940 50 H 0.065 0.935 51 H 0.062 0.938 52 H 0.069 0.931 53 H 0.054 0.946 Dipole moment (debyes) X Y Z Total from point charges -3.742 -2.214 23.738 24.133 hybrid contribution -0.069 0.153 -0.755 0.773 sum -3.812 -2.061 22.983 23.388 Atomic orbital electron populations 1.21691 0.92481 1.02964 1.02738 1.22442 0.95306 0.96075 0.85260 0.91729 1.22519 0.80201 1.01961 0.99790 1.48218 1.06731 1.62689 1.16965 1.20724 0.86780 0.77647 0.84558 1.90696 1.14697 1.55434 1.80979 1.21780 1.00193 1.02201 0.93578 1.20854 0.96274 1.01550 0.91621 1.21535 0.97318 0.97040 0.97056 1.21483 0.94312 1.00664 1.01439 1.23197 1.03331 0.91379 1.03376 1.23031 0.98162 1.03969 0.96117 1.22140 1.00564 0.97287 0.82581 1.22390 0.81028 0.96483 0.94911 1.23354 0.99931 0.95733 1.01301 1.21095 1.00709 0.87786 0.84252 1.44411 1.85828 1.00081 1.02218 1.19001 1.42745 0.92866 0.80670 0.90440 1.24788 1.02300 0.95272 0.98129 1.69957 1.56849 0.96942 1.33320 0.86186 0.77984 0.76583 0.84921 1.21197 1.21676 1.21196 1.21336 0.97508 0.92598 0.85993 1.87912 1.20183 1.76551 1.41570 0.91763 0.91314 0.91042 0.89240 0.90515 0.88802 0.88598 0.90820 0.90881 0.92206 0.92494 0.91469 0.88731 0.88927 0.88921 0.88758 0.91053 0.88061 0.91712 0.92114 0.94006 0.93546 0.93794 0.93144 0.94595 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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.14 -1.51 9.57 37.16 0.36 -1.16 16 2 C 0.07 0.90 3.20 -26.59 -0.08 0.82 16 3 H 0.06 0.92 6.25 -51.93 -0.32 0.60 16 4 C 0.08 0.83 5.82 -3.53 -0.02 0.81 16 5 N -0.61 -8.13 2.95 -172.70 -0.51 -8.64 16 6 C 0.51 6.89 7.71 -10.98 -0.08 6.81 16 7 O -0.54 -9.27 15.47 5.56 0.09 -9.18 16 8 C -0.14 -1.27 3.08 -27.88 -0.09 -1.35 16 9 C -0.07 -0.63 4.69 -26.73 -0.13 -0.75 16 10 C -0.13 -1.00 2.46 -154.51 -0.38 -1.39 16 11 C -0.12 -0.76 9.49 37.16 0.35 -0.41 16 12 C -0.18 -1.32 9.82 -26.69 -0.26 -1.58 16 13 C -0.18 -1.30 9.97 -26.69 -0.27 -1.57 16 14 C 0.10 1.51 5.30 -3.92 -0.02 1.48 16 15 C 0.09 1.70 0.84 -89.64 -0.08 1.63 16 16 C -0.16 -3.28 5.88 -25.25 -0.15 -3.43 16 17 C 0.19 4.70 5.38 -2.95 -0.02 4.68 16 18 N -0.60 -15.74 3.40 -170.31 -0.58 -16.32 16 19 C -0.22 -6.40 10.28 -15.06 -0.15 -6.55 16 20 H 0.08 2.11 7.83 -52.49 -0.41 1.70 16 21 C -0.01 -0.26 5.49 -17.43 -0.10 -0.35 16 22 N -0.92 -27.68 10.31 55.49 0.57 -27.11 16 23 Si 0.91 22.51 29.55 -169.99 -5.02 17.49 16 24 H -0.29 -6.93 7.11 56.52 0.40 -6.53 16 25 H -0.29 -7.22 7.11 56.52 0.40 -6.82 16 26 H -0.29 -6.92 7.11 56.52 0.40 -6.52 16 27 C 0.15 3.25 6.41 -2.87 -0.02 3.23 16 28 O -0.34 -6.20 9.92 -35.22 -0.35 -6.55 16 29 H 0.06 0.65 8.14 -51.93 -0.42 0.23 16 30 H 0.07 0.53 8.14 -51.93 -0.42 0.10 16 31 H 0.07 0.88 8.14 -51.93 -0.42 0.46 16 32 H 0.09 0.67 6.05 -51.93 -0.31 0.36 16 33 H 0.08 0.75 8.14 -51.93 -0.42 0.33 16 34 H 0.09 0.76 7.67 -51.93 -0.40 0.36 16 35 H 0.10 0.69 6.31 -51.93 -0.33 0.36 16 36 H 0.07 0.78 8.10 -51.93 -0.42 0.36 16 37 H 0.07 0.81 7.86 -51.93 -0.41 0.40 16 38 H 0.06 0.36 8.14 -51.93 -0.42 -0.07 16 39 H 0.06 0.32 7.90 -51.93 -0.41 -0.09 16 40 H 0.07 0.37 5.59 -51.93 -0.29 0.08 16 41 H 0.09 0.74 8.05 -51.93 -0.42 0.33 16 42 H 0.09 0.61 8.05 -51.93 -0.42 0.19 16 43 H 0.09 0.60 8.14 -51.93 -0.42 0.17 16 44 H 0.09 0.71 8.14 -51.93 -0.42 0.29 16 45 H 0.07 1.03 8.14 -51.93 -0.42 0.60 16 46 H 0.10 1.70 6.85 -51.93 -0.36 1.34 16 47 H 0.06 1.27 7.82 -51.93 -0.41 0.87 16 48 H 0.06 1.13 6.10 -51.93 -0.32 0.82 16 49 H 0.04 1.16 7.50 -51.93 -0.39 0.78 16 50 H 0.05 1.30 7.38 -51.93 -0.38 0.92 16 51 H 0.25 7.34 8.31 -40.82 -0.34 7.01 16 52 H 0.05 1.06 8.14 -51.93 -0.42 0.64 16 53 H 0.04 0.74 8.00 -51.92 -0.42 0.32 16 LS Contribution 403.23 15.07 6.08 6.08 Total: -1.00 -33.49 403.23 -10.21 -43.70 By element: Atomic # 1 Polarization: 8.95 SS G_CDS: -9.35 Total: -0.40 kcal Atomic # 6 Polarization: 2.07 SS G_CDS: -1.13 Total: 0.94 kcal Atomic # 7 Polarization: -51.55 SS G_CDS: -0.52 Total: -52.07 kcal Atomic # 8 Polarization: -15.47 SS G_CDS: -0.26 Total: -15.73 kcal Atomic # 14 Polarization: 22.51 SS G_CDS: -5.02 Total: 17.49 kcal Total LS contribution 6.08 Total: 6.08 kcal Total: -33.49 -10.21 -43.70 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. ZINC001131692542.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 17.817 kcal (2) G-P(sol) polarization free energy of solvation -33.492 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -15.676 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.206 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.699 kcal (6) G-S(sol) free energy of system = (1) + (5) -25.882 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.15 seconds