Wall clock time and date at job start Tue Mar 9 2021 13:08:13 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 * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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 O 2 2 C 1.42895 * 1 3 3 Si 1.86299 * 109.47599 * 2 1 4 4 C 1.50710 * 109.47121 * 120.00009 * 2 1 3 5 5 C 1.38236 * 119.93898 * 244.97827 * 4 2 1 6 6 C 1.38192 * 120.06682 * 179.69814 * 5 4 2 7 7 C 1.38678 * 119.95252 * 0.60002 * 6 5 4 8 8 O 1.35922 * 120.05505 * 179.70058 * 7 6 5 9 9 C 1.42891 * 117.00122 * 359.71940 * 8 7 6 10 10 C 1.50706 * 109.47034 * 179.97438 * 9 8 7 11 11 O 1.21282 * 119.99788 * 0.02562 * 10 9 8 12 12 N 1.34777 * 119.99710 * 180.02562 * 10 9 8 13 13 C 1.46507 * 119.99689 * 180.02562 * 12 10 9 14 14 H 1.09001 * 109.47073 * 34.99816 * 13 12 10 15 15 C 1.52994 * 109.47066 * 155.00169 * 13 12 10 16 16 C 1.53002 * 109.46917 * 180.02562 * 15 13 12 17 17 C 1.52996 * 109.47124 * 299.99452 * 16 15 13 18 18 H 1.08997 * 109.47566 * 300.00558 * 17 16 15 19 19 N 1.46498 * 109.47627 * 180.02562 * 17 16 15 20 20 C 1.46451 * 126.08163 * 300.00727 * 19 17 16 21 21 C 1.53200 * 113.57016 * 181.76496 * 20 19 17 22 22 C 1.55352 * 106.34440 * 332.23815 * 21 20 19 23 23 S 1.67282 * 126.08456 * 119.95662 * 19 17 16 24 24 O 1.42096 * 105.41147 * 314.82146 * 23 19 17 25 25 O 1.42097 * 105.40736 * 85.08275 * 23 19 17 26 26 C 1.53005 * 109.47055 * 60.00046 * 17 16 15 27 27 C 1.52996 * 109.47137 * 300.00452 * 26 17 16 28 28 C 1.38829 * 119.88279 * 359.42350 * 7 6 5 29 29 Br 1.89105 * 120.03660 * 180.25897 * 28 7 6 30 30 C 1.38192 * 119.92778 * 0.27491 * 28 7 6 31 31 H 0.96698 * 113.99669 * 300.00088 * 1 2 3 32 32 H 1.09002 * 109.47071 * 239.99741 * 2 1 3 33 33 H 1.48499 * 109.99666 * 300.00129 * 3 2 1 34 34 H 1.48503 * 109.99846 * 61.31542 * 3 2 1 35 35 H 1.07995 * 119.96822 * 359.97438 * 5 4 2 36 36 H 1.07993 * 120.02273 * 180.31767 * 6 5 4 37 37 H 1.09004 * 109.47121 * 60.00611 * 9 8 7 38 38 H 1.08997 * 109.47377 * 300.00302 * 9 8 7 39 39 H 0.96994 * 120.00255 * 359.97438 * 12 10 9 40 40 H 1.09001 * 109.47463 * 59.99886 * 15 13 12 41 41 H 1.09007 * 109.47269 * 299.99371 * 15 13 12 42 42 H 1.09003 * 109.46833 * 179.97438 * 16 15 13 43 43 H 1.08997 * 109.47354 * 60.00140 * 16 15 13 44 44 H 1.09002 * 108.64642 * 302.78736 * 20 19 17 45 45 H 1.08999 * 108.64561 * 60.74530 * 20 19 17 46 46 H 1.08995 * 110.10365 * 91.55325 * 21 20 19 47 47 H 1.08993 * 110.10394 * 213.16883 * 21 20 19 48 48 H 1.08998 * 111.49741 * 279.58252 * 22 21 20 49 49 H 1.09001 * 111.33660 * 154.71785 * 22 21 20 50 50 H 1.09001 * 109.46984 * 59.99904 * 26 17 16 51 51 H 1.08999 * 109.47593 * 179.97438 * 26 17 16 52 52 H 1.09008 * 109.47142 * 299.99614 * 27 26 17 53 53 H 1.09004 * 109.47351 * 179.97438 * 27 26 17 54 54 H 1.08002 * 119.97564 * 179.97438 * 30 28 7 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 8 0.0000 0.0000 0.0000 2 6 1.4289 0.0000 0.0000 3 14 2.0501 1.7564 0.0000 4 6 1.9313 -0.7105 1.2305 5 6 2.6390 -0.0112 2.1903 6 6 3.0954 -0.6586 3.3226 7 6 2.8540 -2.0139 3.4901 8 8 3.3073 -2.6541 4.6001 9 6 4.0308 -1.8704 5.5510 10 6 4.4498 -2.7430 6.7061 11 8 4.1580 -3.9201 6.7139 12 7 5.1489 -2.2151 7.7304 13 6 5.5567 -3.0634 8.8531 14 1 5.7928 -4.0624 8.4865 15 6 6.7925 -2.4613 9.5246 16 6 7.2188 -3.3476 10.6966 17 6 6.0777 -3.4339 11.7122 18 1 5.8415 -2.4351 12.0791 19 7 6.4860 -4.2821 12.8347 20 6 6.8822 -5.6896 12.7528 21 6 7.2141 -6.3127 14.1124 22 6 7.7075 -5.1503 15.0173 23 16 6.5837 -3.8148 14.4380 24 8 7.2384 -2.5537 14.4246 25 8 5.2717 -3.9768 14.9591 26 6 4.8418 -4.0360 11.0406 27 6 4.4156 -3.1498 9.8685 28 6 2.1435 -2.7144 2.5247 29 35 1.8054 -4.5610 2.7521 30 6 1.6839 -2.0607 1.3972 31 1 -0.3933 0.4417 0.7650 32 1 1.7923 -0.5139 -0.8900 33 1 1.5616 2.4678 -1.2085 34 1 1.5879 2.4585 1.2242 35 1 2.8314 1.0431 2.0573 36 1 3.6438 -0.1105 4.0744 37 1 3.3945 -1.0646 5.9170 38 1 4.9154 -1.4472 5.0752 39 1 5.3825 -1.2737 7.7239 40 1 6.5564 -1.4623 9.8912 41 1 7.6057 -2.4002 8.8013 42 1 8.0995 -2.9188 11.1748 43 1 7.4551 -4.3464 10.3297 44 1 6.0679 -6.2550 12.2997 45 1 7.7591 -5.7679 12.1102 46 1 6.3233 -6.7708 14.5420 47 1 8.0022 -7.0573 14.0009 48 1 8.7509 -4.9051 14.8191 49 1 7.5481 -5.3738 16.0722 50 1 5.0780 -5.0348 10.6737 51 1 4.0289 -4.0980 11.7641 52 1 4.1794 -2.1509 10.2354 53 1 3.5353 -3.5789 9.3898 54 1 1.1314 -2.6049 0.6455 There are 77 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 77 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). S: (AM1): M.J.S. DEWAR, Y.-C. YUAN, INORGANIC CHEM., 29, 589 (1990) Br: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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=WATER s_22_16729364_6201210.mol2 54 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 9 2021 13:08:13 Heat of formation + Delta-G solvation = -132.254360 kcal Electronic energy + Delta-G solvation = -38941.059912 eV Core-core repulsion = 33632.391501 eV Total energy + Delta-G solvation = -5308.668411 eV No. of doubly occupied orbitals = 77 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 489.075 amu Computer time = 2.02 seconds Orbital eigenvalues (eV) -41.71052 -40.98946 -39.82863 -39.06029 -38.28588 -37.61965 -37.30415 -36.97891 -35.39654 -33.63184 -33.16246 -31.94450 -31.83586 -30.66127 -30.09913 -28.57343 -27.20816 -26.31357 -24.26063 -23.80978 -23.57029 -22.77091 -22.34153 -21.57388 -20.71613 -19.92647 -19.57867 -18.60383 -18.36721 -18.18868 -17.81478 -17.48070 -17.27021 -17.09629 -16.93663 -16.43146 -16.25769 -16.02895 -15.80855 -15.69726 -15.64101 -15.23517 -15.20238 -14.93843 -14.70900 -14.56979 -14.48344 -14.45659 -14.03664 -13.89247 -13.78951 -13.63261 -13.57757 -13.44837 -13.26572 -13.21557 -13.13712 -12.68695 -12.43192 -12.40026 -12.35159 -12.20947 -12.17483 -12.15240 -12.00661 -11.88716 -11.69363 -11.54779 -11.48509 -11.26788 -11.23490 -10.92939 -10.60024 -10.14984 -9.96640 -9.84905 -9.32905 -4.89195 -0.80461 -0.05181 0.12233 0.64352 1.33551 1.40351 1.48644 1.49461 1.61831 2.01316 2.39258 2.52117 2.58822 2.93509 3.28437 3.41084 3.49209 3.57855 3.60595 3.62321 3.70259 3.76326 3.82482 3.91439 3.95727 4.03797 4.07337 4.09298 4.12271 4.17815 4.19467 4.33141 4.35048 4.41513 4.44419 4.51789 4.58312 4.62098 4.65191 4.68439 4.71555 4.85615 4.90809 4.94797 4.98885 5.00663 5.07219 5.11077 5.12892 5.16991 5.23711 5.32109 5.37544 5.46488 5.50524 5.68259 6.07409 6.49587 6.50502 7.18050 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 O -0.551 6.551 2 C -0.084 4.084 3 Si 0.778 3.222 4 C -0.149 4.149 5 C -0.056 4.056 6 C -0.189 4.189 7 C 0.149 3.851 8 O -0.305 6.305 9 C 0.031 3.969 10 C 0.512 3.488 11 O -0.567 6.567 12 N -0.706 5.706 13 C 0.149 3.851 14 H 0.072 0.928 15 C -0.117 4.117 16 C -0.121 4.121 17 C 0.183 3.817 18 H 0.083 0.917 19 N -0.984 5.984 20 C 0.069 3.931 21 C -0.126 4.126 22 C -0.557 4.557 23 S 2.542 3.458 24 O -0.999 6.999 25 O -1.000 7.000 26 C -0.124 4.124 27 C -0.125 4.125 28 C -0.144 4.144 29 Br -0.055 7.055 30 C -0.035 4.035 31 H 0.387 0.613 32 H 0.122 0.878 33 H -0.308 1.308 34 H -0.270 1.270 35 H 0.164 0.836 36 H 0.179 0.821 37 H 0.128 0.872 38 H 0.134 0.866 39 H 0.421 0.579 40 H 0.087 0.913 41 H 0.090 0.910 42 H 0.081 0.919 43 H 0.078 0.922 44 H 0.081 0.919 45 H 0.094 0.906 46 H 0.096 0.904 47 H 0.140 0.860 48 H 0.165 0.835 49 H 0.159 0.841 50 H 0.073 0.927 51 H 0.059 0.941 52 H 0.078 0.922 53 H 0.050 0.950 54 H 0.136 0.864 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 9.423 4.036 0.826 10.284 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 O -0.359 6.359 2 C -0.193 4.193 3 Si 0.682 3.318 4 C -0.151 4.151 5 C -0.074 4.074 6 C -0.208 4.208 7 C 0.101 3.899 8 O -0.217 6.217 9 C -0.046 4.046 10 C 0.298 3.702 11 O -0.445 6.445 12 N -0.360 5.360 13 C 0.045 3.955 14 H 0.090 0.910 15 C -0.156 4.156 16 C -0.161 4.161 17 C 0.079 3.921 18 H 0.101 0.899 19 N -0.825 5.825 20 C -0.056 4.056 21 C -0.165 4.165 22 C -0.692 4.692 23 S 2.645 3.355 24 O -0.989 6.989 25 O -0.990 6.990 26 C -0.164 4.164 27 C -0.164 4.164 28 C -0.224 4.224 29 Br 0.030 6.970 30 C -0.055 4.055 31 H 0.236 0.764 32 H 0.140 0.860 33 H -0.237 1.237 34 H -0.197 1.197 35 H 0.181 0.819 36 H 0.197 0.803 37 H 0.146 0.854 38 H 0.152 0.848 39 H 0.259 0.741 40 H 0.105 0.895 41 H 0.109 0.891 42 H 0.100 0.900 43 H 0.097 0.903 44 H 0.099 0.901 45 H 0.113 0.887 46 H 0.115 0.885 47 H 0.159 0.841 48 H 0.183 0.817 49 H 0.177 0.823 50 H 0.091 0.909 51 H 0.078 0.922 52 H 0.097 0.903 53 H 0.069 0.931 54 H 0.154 0.846 Dipole moment (debyes) X Y Z Total from point charges 8.447 3.684 0.854 9.255 hybrid contribution -1.187 -0.227 0.347 1.257 sum 7.260 3.457 1.201 8.130 Atomic orbital electron populations 1.86594 1.18684 1.80440 1.50188 1.25321 0.83376 1.12402 0.98151 0.94095 0.89868 0.71095 0.76746 1.20280 1.02563 0.93598 0.98634 1.21291 0.93592 1.00859 0.91623 1.21560 1.06089 0.91522 1.01660 1.18812 0.93230 0.91737 0.86094 1.86248 1.68658 1.37920 1.28922 1.22877 0.99675 0.93689 0.88370 1.20497 0.77931 0.87995 0.83758 1.90789 1.55345 1.18702 1.79638 1.45736 1.55985 1.14506 1.19789 1.20981 0.94226 0.95089 0.85243 0.91022 1.21754 0.96350 1.00064 0.97385 1.21821 0.97829 1.01529 0.94896 1.20705 0.94061 0.91987 0.85358 0.89925 1.57048 1.83496 1.18674 1.23302 1.21785 1.03075 0.81193 0.99563 1.22538 1.05152 0.97661 0.91179 1.30558 1.14908 1.12893 1.10820 1.07015 0.76613 0.76507 0.75382 1.93547 1.76916 1.42578 1.85851 1.93601 1.38561 1.87381 1.79431 1.21829 0.97454 1.00120 0.96971 1.21858 0.98210 1.00900 0.95464 1.21277 1.09772 0.91907 0.99485 1.96760 1.94964 1.07538 1.97722 1.20890 0.95727 0.94457 0.94406 0.76417 0.85988 1.23667 1.19656 0.81898 0.80335 0.85416 0.84839 0.74063 0.89461 0.89139 0.90037 0.90316 0.90079 0.88731 0.88491 0.84145 0.81679 0.82271 0.90871 0.92174 0.90342 0.93100 0.84585 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 O -0.55 -0.77 12.83 -148.98 -1.91 -2.68 16 2 C -0.08 0.23 2.80 71.24 0.20 0.43 16 3 Si 0.78 -4.34 29.88 68.60 2.05 -2.29 16 4 C -0.15 0.30 5.35 -19.85 -0.11 0.20 16 5 C -0.06 0.34 8.66 22.26 0.19 0.54 16 6 C -0.19 1.11 9.04 22.36 0.20 1.31 16 7 C 0.15 0.27 6.68 22.52 0.15 0.43 16 8 O -0.31 -1.25 9.01 -73.64 -0.66 -1.92 16 9 C 0.03 -0.18 5.22 71.24 0.37 0.19 16 10 C 0.51 0.28 7.87 87.66 0.69 0.97 16 11 O -0.57 -7.14 16.16 16.00 0.26 -6.88 16 12 N -0.71 4.58 5.08 -442.52 -2.25 2.33 16 13 C 0.15 -0.57 2.45 44.99 0.11 -0.46 16 14 H 0.07 0.08 7.58 -2.39 -0.02 0.06 16 15 C -0.12 0.97 5.39 30.59 0.16 1.14 16 16 C -0.12 0.73 5.25 30.59 0.16 0.89 16 17 C 0.18 -0.33 2.86 44.99 0.13 -0.20 16 18 H 0.08 -0.03 8.10 -2.39 -0.02 -0.05 16 19 N -0.98 1.89 3.14 -493.15 -1.55 0.34 16 20 C 0.07 -0.54 6.20 86.46 0.54 -0.01 16 21 C -0.13 1.36 6.59 31.47 0.21 1.57 16 22 C -0.56 5.06 7.65 72.79 0.56 5.61 16 23 S 2.54 4.55 6.15 -56.49 -0.35 4.20 16 24 O -1.00 -8.87 18.04 -127.81 -2.31 -11.18 16 25 O -1.00 -10.53 18.14 -127.81 -2.32 -12.85 16 26 C -0.12 0.08 5.24 30.59 0.16 0.24 16 27 C -0.12 0.15 5.50 30.59 0.17 0.32 16 28 C -0.14 -0.78 6.24 22.40 0.14 -0.64 16 29 Br -0.06 -0.52 33.16 -20.37 -0.68 -1.19 16 30 C -0.03 -0.09 9.42 22.27 0.21 0.12 16 31 H 0.39 -1.73 8.70 -74.06 -0.64 -2.38 16 32 H 0.12 -0.43 7.96 -2.39 -0.02 -0.45 16 33 H -0.31 -1.05 7.11 99.48 0.71 -0.35 16 34 H -0.27 0.98 7.01 99.48 0.70 1.67 16 35 H 0.16 -1.90 4.83 -2.91 -0.01 -1.91 16 36 H 0.18 -2.53 6.30 -2.91 -0.02 -2.55 16 37 H 0.13 -1.58 7.65 -2.38 -0.02 -1.60 16 38 H 0.13 -1.77 7.67 -2.39 -0.02 -1.79 16 39 H 0.42 -6.04 8.60 -92.71 -0.80 -6.84 16 40 H 0.09 -0.90 8.14 -2.39 -0.02 -0.92 16 41 H 0.09 -0.92 8.14 -2.38 -0.02 -0.94 16 42 H 0.08 -0.47 8.14 -2.39 -0.02 -0.49 16 43 H 0.08 -0.54 7.69 -2.39 -0.02 -0.56 16 44 H 0.08 -0.60 7.60 -2.39 -0.02 -0.62 16 45 H 0.09 -1.06 7.74 -2.39 -0.02 -1.08 16 46 H 0.10 -0.99 8.14 -2.39 -0.02 -1.01 16 47 H 0.14 -2.20 8.14 -2.39 -0.02 -2.22 16 48 H 0.17 -2.12 8.14 -2.39 -0.02 -2.14 16 49 H 0.16 -1.90 8.14 -2.39 -0.02 -1.92 16 50 H 0.07 -0.06 7.60 -2.39 -0.02 -0.08 16 51 H 0.06 0.15 8.14 -2.39 -0.02 0.13 16 52 H 0.08 -0.26 8.14 -2.38 -0.02 -0.28 16 53 H 0.05 0.21 8.14 -2.38 -0.02 0.19 16 54 H 0.14 0.41 8.06 -2.91 -0.02 0.38 16 Total: 0.00 -41.30 461.60 -5.92 -47.22 By element: Atomic # 1 Polarization: -27.28 SS G_CDS: -0.45 Total: -27.73 kcal Atomic # 6 Polarization: 8.39 SS G_CDS: 4.24 Total: 12.63 kcal Atomic # 7 Polarization: 6.47 SS G_CDS: -3.79 Total: 2.68 kcal Atomic # 8 Polarization: -28.57 SS G_CDS: -6.94 Total: -35.51 kcal Atomic # 14 Polarization: -4.34 SS G_CDS: 2.05 Total: -2.29 kcal Atomic # 16 Polarization: 4.55 SS G_CDS: -0.35 Total: 4.20 kcal Atomic # 35 Polarization: -0.52 SS G_CDS: -0.68 Total: -1.19 kcal Total: -41.30 -5.92 -47.22 kcal The number of atoms in the molecule is 54 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_16729364_6201210.mol2 54 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 -85.035 kcal (2) G-P(sol) polarization free energy of solvation -41.301 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -126.336 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -5.918 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -47.220 kcal (6) G-S(sol) free energy of system = (1) + (5) -132.254 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 2.02 seconds