Wall clock time and date at job start Wed Apr 1 2020 04:13:41 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.53007 * 1 3 3 C 1.53000 * 109.47393 * 2 1 4 4 H 1.08994 * 109.47107 * 185.59800 * 3 2 1 5 5 C 1.50705 * 109.47053 * 305.60458 * 3 2 1 6 6 H 1.07999 * 120.00329 * 119.99874 * 5 3 2 7 7 C 1.37881 * 119.99606 * 299.99890 * 5 3 2 8 8 N 1.31515 * 119.99715 * 0.02562 * 7 5 3 9 9 Si 1.86797 * 120.00040 * 180.02562 * 7 5 3 10 10 H 1.48497 * 109.47478 * 89.99818 * 9 7 5 11 11 H 1.48503 * 109.47188 * 210.00205 * 9 7 5 12 12 H 1.48507 * 109.46808 * 329.99613 * 9 7 5 13 13 N 1.46901 * 109.47012 * 65.59949 * 3 2 1 14 14 C 1.46892 * 110.99904 * 87.52460 * 13 3 2 15 15 C 1.50706 * 109.47370 * 184.28361 * 14 13 3 16 16 O 1.21279 * 119.99692 * 359.97438 * 15 14 13 17 17 N 1.34774 * 120.00449 * 179.97438 * 15 14 13 18 18 C 1.46496 * 120.00146 * 179.97438 * 17 15 14 19 19 C 1.53002 * 109.47180 * 180.02562 * 18 17 15 20 20 H 1.09000 * 110.87555 * 57.39054 * 19 18 17 21 21 C 1.55141 * 110.98980 * 293.55066 * 19 18 17 22 22 C 1.54320 * 103.01697 * 206.21962 * 21 19 18 23 23 N 1.47438 * 107.23305 * 22.27897 * 22 21 19 24 24 C 1.39677 * 110.99889 * 123.26487 * 23 22 21 25 25 C 1.38913 * 120.07734 * 36.48354 * 24 23 22 26 26 C 1.38101 * 119.92640 * 179.97438 * 25 24 23 27 27 C 1.38283 * 120.07456 * 0.02562 * 26 25 24 28 28 C 1.38286 * 120.14917 * 0.27475 * 27 26 25 29 29 C 1.38104 * 120.07511 * 359.45072 * 28 27 26 30 30 C 1.47838 * 108.57061 * 0.99335 * 23 22 21 31 31 H 1.08993 * 109.47479 * 183.51109 * 1 2 3 32 32 H 1.09001 * 109.47076 * 303.51327 * 1 2 3 33 33 H 1.09004 * 109.46621 * 63.50877 * 1 2 3 34 34 H 1.08997 * 109.46996 * 239.99512 * 2 1 3 35 35 H 1.08999 * 109.46844 * 119.99683 * 2 1 3 36 36 H 0.97011 * 119.99583 * 180.02562 * 8 7 5 37 37 H 1.00905 * 111.00239 * 323.57038 * 13 3 2 38 38 H 1.08995 * 109.47484 * 64.28546 * 14 13 3 39 39 H 1.09006 * 109.47312 * 304.28036 * 14 13 3 40 40 H 0.97004 * 120.00032 * 0.02562 * 17 15 14 41 41 H 1.09003 * 109.47005 * 60.00326 * 18 17 15 42 42 H 1.08997 * 109.47490 * 300.00276 * 18 17 15 43 43 H 1.08993 * 110.70974 * 87.83128 * 21 19 18 44 44 H 1.09001 * 110.85102 * 324.66493 * 21 19 18 45 45 H 1.09004 * 109.88167 * 262.87708 * 22 21 19 46 46 H 1.09002 * 109.89125 * 141.68924 * 22 21 19 47 47 H 1.08003 * 120.03756 * 359.96088 * 25 24 23 48 48 H 1.08003 * 119.96410 * 180.02562 * 26 25 24 49 49 H 1.08006 * 119.92545 * 179.97438 * 27 26 25 50 50 H 1.08004 * 119.95829 * 179.69997 * 28 27 26 51 51 H 1.07999 * 120.03611 * 180.27276 * 29 28 27 52 52 H 1.09007 * 110.37452 * 94.73091 * 30 23 22 53 53 H 1.09000 * 110.37746 * 217.07537 * 30 23 22 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 1 3.1255 1.4441 -0.1002 5 6 1.4277 2.1918 -1.1552 6 1 0.8199 3.0653 -0.9709 7 6 1.6434 1.7622 -2.4475 8 7 2.3832 0.6982 -2.6719 9 14 0.8850 2.6914 -3.8796 10 1 1.8238 3.7467 -4.3378 11 1 0.6173 1.7486 -4.9952 12 1 -0.3893 3.3174 -3.4441 13 7 1.6639 2.0949 1.2613 14 6 2.6792 1.8627 2.2972 15 6 2.2108 2.4553 3.6014 16 8 1.1372 3.0162 3.6617 17 7 2.9850 2.3619 4.7006 18 6 2.5299 2.9385 5.9682 19 6 3.5873 2.6966 7.0473 20 1 4.5477 3.1189 6.7516 21 6 3.7199 1.1841 7.3662 22 6 4.2378 1.1796 8.8200 23 7 3.8682 2.4750 9.4193 24 6 3.1034 2.2923 10.5738 25 6 3.3753 1.2314 11.4283 26 6 2.6175 1.0541 12.5692 27 6 1.5888 1.9308 12.8610 28 6 1.3115 2.9844 12.0093 29 6 2.0689 3.1709 10.8697 30 6 3.1279 3.2653 8.4129 31 1 -0.3634 -1.0256 0.0629 32 1 -0.3633 0.5674 0.8568 33 1 -0.3633 0.4584 -0.9198 34 1 1.8934 -0.5139 -0.8899 35 1 1.8933 -0.5138 0.8900 36 1 2.5352 0.3962 -3.5811 37 1 0.7544 1.7864 1.5708 38 1 3.6160 2.3338 1.9999 39 1 2.8344 0.7908 2.4202 40 1 3.8439 1.9136 4.6523 41 1 1.5927 2.4681 6.2655 42 1 2.3755 4.0104 5.8452 43 1 4.4404 0.7112 6.6990 44 1 2.7519 0.6874 7.3004 45 1 5.3216 1.0634 8.8271 46 1 3.7718 0.3668 9.3771 47 1 4.1782 0.5459 11.2007 48 1 2.8286 0.2297 13.2342 49 1 0.9974 1.7896 13.7537 50 1 0.5042 3.6645 12.2377 51 1 1.8546 3.9965 10.2072 52 1 2.0537 3.1292 8.5387 53 1 3.3893 4.3206 8.4918 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 ZINC000579855194.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 04:13:41 Heat of formation + Delta-G solvation = -1.145392 kcal Electronic energy + Delta-G solvation = -29384.772305 eV Core-core repulsion = 25449.020431 eV Total energy + Delta-G solvation = -3935.751875 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 345.216 amu Computer time = 1.03 seconds Orbital eigenvalues (eV) -40.94352 -39.47389 -38.20959 -36.31615 -35.12230 -32.75987 -32.44535 -31.70666 -30.71335 -29.87068 -28.89335 -27.32582 -25.89766 -24.30377 -23.41101 -23.16483 -22.56757 -21.08484 -20.55871 -20.13317 -19.30363 -18.05508 -17.04664 -16.85172 -16.09171 -15.84504 -15.56653 -15.40564 -15.33373 -14.77972 -14.63397 -14.61891 -14.17670 -13.79316 -13.64129 -13.55625 -13.38950 -13.22047 -12.91008 -12.72284 -12.58662 -12.27533 -12.13959 -11.94236 -11.81453 -11.66687 -11.56020 -11.27867 -11.24457 -11.09504 -10.89939 -10.76991 -10.51060 -10.18242 -10.00573 -9.89725 -9.56882 -9.49644 -9.28558 -9.04299 -8.96788 -8.33080 -8.20063 -7.23021 -5.96753 -3.97331 1.04104 1.04310 2.93905 3.11226 3.25550 3.34979 3.41255 3.46695 3.60804 4.03480 4.15292 4.34994 4.47043 4.52500 4.60518 4.75682 4.80274 4.86030 4.89588 5.06652 5.10085 5.15178 5.17606 5.25552 5.33328 5.35095 5.38419 5.39141 5.44013 5.57180 5.60255 5.66615 5.78761 5.84726 5.87355 5.89412 6.05692 6.21854 6.33233 6.37592 6.39480 6.45313 6.52620 6.56419 6.89041 6.89366 7.20428 7.25088 7.49126 7.57121 7.69401 7.91160 8.26409 8.44822 8.55555 8.61684 9.09101 9.66249 11.10513 Molecular weight = 345.22amu Principal moments of inertia in cm(-1) A = 0.027218 B = 0.001584 C = 0.001575 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1028.500247 B =17669.279179 C =17772.382226 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.128 4.128 2 C -0.145 4.145 3 C 0.215 3.785 4 H 0.033 0.967 5 C -0.486 4.486 6 H 0.063 0.937 7 C 0.056 3.944 8 N -0.898 5.898 9 Si 0.895 3.105 10 H -0.282 1.282 11 H -0.289 1.289 12 H -0.274 1.274 13 N -0.669 5.669 14 C 0.043 3.957 15 C 0.513 3.487 16 O -0.526 6.526 17 N -0.739 5.739 18 C 0.133 3.867 19 C -0.110 4.110 20 H 0.095 0.905 21 C -0.138 4.138 22 C 0.037 3.963 23 N -0.477 5.477 24 C 0.113 3.887 25 C -0.159 4.159 26 C -0.093 4.093 27 C -0.154 4.154 28 C -0.091 4.091 29 C -0.165 4.165 30 C 0.047 3.953 31 H 0.022 0.978 32 H 0.032 0.968 33 H 0.066 0.934 34 H 0.094 0.906 35 H 0.009 0.991 36 H 0.259 0.741 37 H 0.337 0.663 38 H 0.094 0.906 39 H 0.047 0.953 40 H 0.399 0.601 41 H 0.065 0.935 42 H 0.069 0.931 43 H 0.083 0.917 44 H 0.075 0.925 45 H 0.085 0.915 46 H 0.071 0.929 47 H 0.124 0.876 48 H 0.122 0.878 49 H 0.119 0.881 50 H 0.122 0.878 51 H 0.126 0.874 52 H 0.068 0.932 53 H 0.091 0.909 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 4.563 -0.813 29.570 29.931 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.186 4.186 2 C -0.185 4.185 3 C 0.107 3.893 4 H 0.051 0.949 5 C -0.525 4.525 6 H 0.081 0.919 7 C -0.144 4.144 8 N -0.538 5.538 9 Si 0.723 3.277 10 H -0.208 1.208 11 H -0.215 1.215 12 H -0.199 1.199 13 N -0.302 5.302 14 C -0.090 4.090 15 C 0.299 3.701 16 O -0.400 6.400 17 N -0.394 5.394 18 C 0.008 3.992 19 C -0.129 4.129 20 H 0.113 0.887 21 C -0.176 4.176 22 C -0.088 4.088 23 N -0.190 5.190 24 C 0.014 3.986 25 C -0.179 4.179 26 C -0.111 4.111 27 C -0.173 4.173 28 C -0.109 4.109 29 C -0.186 4.186 30 C -0.080 4.080 31 H 0.041 0.959 32 H 0.051 0.949 33 H 0.085 0.915 34 H 0.112 0.888 35 H 0.028 0.972 36 H 0.068 0.932 37 H 0.157 0.843 38 H 0.112 0.888 39 H 0.065 0.935 40 H 0.232 0.768 41 H 0.084 0.916 42 H 0.088 0.912 43 H 0.101 0.899 44 H 0.094 0.906 45 H 0.103 0.897 46 H 0.090 0.910 47 H 0.142 0.858 48 H 0.140 0.860 49 H 0.137 0.863 50 H 0.140 0.860 51 H 0.144 0.856 52 H 0.086 0.914 53 H 0.109 0.891 Dipole moment (debyes) X Y Z Total from point charges 4.575 0.387 29.325 29.683 hybrid contribution -0.538 -1.082 -0.691 1.392 sum 4.037 -0.695 28.635 28.926 Atomic orbital electron populations 1.21568 0.97540 0.97873 1.01613 1.22115 0.95222 1.00954 1.00164 1.19267 0.96054 0.88703 0.85254 0.94876 1.21068 1.30385 1.11500 0.89514 0.91853 1.24907 0.95614 0.95195 0.98709 1.69882 1.39369 1.24087 1.20445 0.86547 0.79712 0.79512 0.81942 1.20831 1.21488 1.19897 1.60679 1.03405 1.63415 1.02706 1.21826 0.98970 0.98324 0.89876 1.21105 0.84445 0.79286 0.85215 1.90681 1.22160 1.42089 1.85068 1.45998 1.23204 1.60748 1.09407 1.21140 0.96265 0.99260 0.82534 1.22568 1.00102 0.95989 0.94228 0.88656 1.22890 1.02514 0.96181 0.96030 1.23125 1.02023 0.89250 0.94407 1.60531 1.37812 1.16084 1.04602 1.19174 0.91316 0.94712 0.93411 1.20882 1.02292 0.98512 0.96259 1.21004 0.95100 0.98765 0.96252 1.21068 0.99315 0.96442 1.00465 1.21028 0.98705 0.97386 0.93826 1.21062 0.96420 1.01446 0.99663 1.23244 0.97259 0.98720 0.88758 0.95900 0.94906 0.91539 0.88768 0.97237 0.93161 0.84302 0.88848 0.93519 0.76762 0.91622 0.91230 0.89881 0.90622 0.89666 0.91023 0.85765 0.86027 0.86258 0.85976 0.85568 0.91384 0.89091 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.13 -2.79 9.38 37.16 0.35 -2.44 16 2 C -0.15 -3.34 5.51 -26.73 -0.15 -3.49 16 3 C 0.21 5.25 2.34 -68.86 -0.16 5.09 16 4 H 0.03 0.87 7.81 -51.93 -0.41 0.47 16 5 C -0.49 -13.23 7.40 -34.44 -0.25 -13.49 16 6 H 0.06 1.79 7.81 -52.49 -0.41 1.38 16 7 C 0.06 1.53 5.24 -17.82 -0.09 1.43 16 8 N -0.90 -25.03 11.38 55.43 0.63 -24.40 16 9 Si 0.89 21.48 29.54 -169.99 -5.02 16.46 16 10 H -0.28 -6.76 7.11 56.52 0.40 -6.36 16 11 H -0.29 -6.88 7.11 56.52 0.40 -6.48 16 12 H -0.27 -6.57 7.11 56.52 0.40 -6.17 16 13 N -0.67 -14.57 5.94 -109.22 -0.65 -15.22 16 14 C 0.04 0.74 5.95 -4.98 -0.03 0.71 16 15 C 0.51 8.17 7.86 -10.98 -0.09 8.09 16 16 O -0.53 -10.54 16.10 5.56 0.09 -10.45 16 17 N -0.74 -8.14 5.41 -60.29 -0.33 -8.47 16 18 C 0.13 1.28 5.12 -4.04 -0.02 1.26 16 19 C -0.11 -0.76 3.10 -88.48 -0.27 -1.04 16 20 H 0.10 0.63 8.14 -51.93 -0.42 0.21 16 21 C -0.14 -0.89 6.11 -24.91 -0.15 -1.04 16 22 C 0.04 0.25 6.18 -2.84 -0.02 0.23 16 23 N -0.48 -3.83 5.27 -167.92 -0.89 -4.71 16 24 C 0.11 0.97 5.78 -83.78 -0.48 0.49 16 25 C -0.16 -1.26 8.59 -39.38 -0.34 -1.59 16 26 C -0.09 -0.67 10.04 -39.62 -0.40 -1.06 16 27 C -0.15 -1.10 10.04 -39.55 -0.40 -1.50 16 28 C -0.09 -0.68 10.04 -39.61 -0.40 -1.08 16 29 C -0.17 -1.36 8.66 -39.38 -0.34 -1.70 16 30 C 0.05 0.34 4.94 -2.30 -0.01 0.33 16 31 H 0.02 0.43 8.14 -51.93 -0.42 0.01 16 32 H 0.03 0.67 6.34 -51.93 -0.33 0.34 16 33 H 0.07 1.58 7.15 -51.93 -0.37 1.20 16 34 H 0.09 2.42 6.24 -51.93 -0.32 2.10 16 35 H 0.01 0.19 7.48 -51.93 -0.39 -0.20 16 36 H 0.26 6.95 8.31 -40.81 -0.34 6.61 16 37 H 0.34 7.45 6.05 -39.29 -0.24 7.21 16 38 H 0.09 1.51 8.13 -51.93 -0.42 1.09 16 39 H 0.05 0.74 7.47 -51.93 -0.39 0.35 16 40 H 0.40 3.30 8.34 -40.82 -0.34 2.96 16 41 H 0.07 0.70 8.11 -51.93 -0.42 0.28 16 42 H 0.07 0.74 8.14 -51.93 -0.42 0.32 16 43 H 0.08 0.45 8.14 -51.93 -0.42 0.03 16 44 H 0.08 0.56 8.02 -51.93 -0.42 0.14 16 45 H 0.09 0.51 8.14 -51.93 -0.42 0.09 16 46 H 0.07 0.47 5.95 -51.93 -0.31 0.16 16 47 H 0.12 0.82 6.78 -52.48 -0.36 0.47 16 48 H 0.12 0.62 8.06 -52.48 -0.42 0.20 16 49 H 0.12 0.63 8.06 -52.48 -0.42 0.21 16 50 H 0.12 0.69 8.06 -52.48 -0.42 0.26 16 51 H 0.13 0.91 6.40 -52.49 -0.34 0.57 16 52 H 0.07 0.51 6.38 -51.93 -0.33 0.18 16 53 H 0.09 0.62 8.14 -51.93 -0.42 0.19 16 LS Contribution 413.05 15.07 6.22 6.22 Total: -1.00 -31.64 413.05 -11.92 -43.56 By element: Atomic # 1 Polarization: 16.54 SS G_CDS: -8.73 Total: 7.82 kcal Atomic # 6 Polarization: -7.56 SS G_CDS: -3.26 Total: -10.81 kcal Atomic # 7 Polarization: -51.57 SS G_CDS: -1.23 Total: -52.80 kcal Atomic # 8 Polarization: -10.54 SS G_CDS: 0.09 Total: -10.45 kcal Atomic # 14 Polarization: 21.48 SS G_CDS: -5.02 Total: 16.46 kcal Total LS contribution 6.22 Total: 6.22 kcal Total: -31.64 -11.92 -43.56 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. ZINC000579855194.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 42.414 kcal (2) G-P(sol) polarization free energy of solvation -31.641 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 10.773 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.919 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.559 kcal (6) G-S(sol) free energy of system = (1) + (5) -1.145 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.03 seconds