Wall clock time and date at job start Tue Mar 31 2020 21:00:44 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.09006 * 109.48186 * 2 1 4 4 C 1.53092 * 109.48104 * 120.02236 * 2 1 3 5 5 C 1.53183 * 109.16395 * 177.60578 * 4 2 1 6 6 H 1.08998 * 109.54550 * 62.97120 * 5 4 2 7 7 C 1.50695 * 109.54659 * 183.15469 * 5 4 2 8 8 O 1.21283 * 119.99979 * 93.17315 * 7 5 4 9 9 N 1.34770 * 119.99981 * 273.17722 * 7 5 4 10 10 C 1.47082 * 125.62881 * 355.33735 * 9 7 5 11 11 C 1.54464 * 107.29944 * 179.04056 * 10 9 7 12 12 C 1.54367 * 102.87363 * 338.49936 * 11 10 9 13 13 C 1.47569 * 125.63142 * 175.37486 * 9 7 5 14 14 C 1.52111 * 110.56734 * 152.60108 * 12 11 10 15 15 C 1.53348 * 109.44075 * 182.79562 * 14 12 11 16 16 N 1.46926 * 108.68533 * 306.69825 * 15 14 12 17 17 C 1.36941 * 120.55012 * 234.36383 * 16 15 14 18 18 H 1.08000 * 119.99775 * 39.04004 * 17 16 15 19 19 C 1.38816 * 119.99861 * 219.04244 * 17 16 15 20 20 N 1.31614 * 120.00402 * 3.06386 * 19 17 16 21 21 Si 1.86799 * 119.99625 * 183.06479 * 19 17 16 22 22 H 1.48501 * 109.47352 * 90.00117 * 21 19 17 23 23 H 1.48501 * 109.47085 * 209.99807 * 21 19 17 24 24 H 1.48495 * 109.47196 * 329.99731 * 21 19 17 25 25 C 1.46936 * 118.89791 * 54.35103 * 16 15 14 26 26 C 1.52116 * 110.55057 * 276.68957 * 12 11 10 27 27 C 1.53180 * 109.10601 * 303.06799 * 5 4 2 28 28 C 1.53104 * 109.16187 * 56.84445 * 27 5 4 29 29 O 1.42897 * 109.49154 * 239.97472 * 2 1 3 30 30 H 1.09004 * 109.46778 * 299.97195 * 1 2 3 31 31 H 1.08994 * 109.47660 * 59.97368 * 1 2 3 32 32 H 1.08995 * 109.47108 * 179.97438 * 1 2 3 33 33 H 1.08995 * 109.52526 * 57.70071 * 4 2 1 34 34 H 1.08997 * 109.57504 * 297.53869 * 4 2 1 35 35 H 1.09004 * 109.88190 * 59.61770 * 10 9 7 36 36 H 1.09000 * 109.87988 * 298.46353 * 10 9 7 37 37 H 1.09000 * 110.73553 * 220.15129 * 11 10 9 38 38 H 1.09001 * 110.88912 * 96.91407 * 11 10 9 39 39 H 1.09008 * 110.37284 * 84.46377 * 13 9 7 40 40 H 1.09002 * 110.38173 * 322.11875 * 13 9 7 41 41 H 1.08994 * 109.47379 * 302.78401 * 14 12 11 42 42 H 1.08994 * 109.47500 * 62.80343 * 14 12 11 43 43 H 1.09000 * 109.59472 * 186.94123 * 15 14 12 44 44 H 1.08996 * 109.73304 * 66.54192 * 15 14 12 45 45 H 0.96995 * 120.00834 * 179.97438 * 20 19 17 46 46 H 1.09003 * 109.60242 * 185.88289 * 25 16 15 47 47 H 1.08998 * 109.73700 * 65.48427 * 25 16 15 48 48 H 1.08995 * 109.47678 * 297.18465 * 26 12 11 49 49 H 1.09007 * 109.46882 * 57.19784 * 26 12 11 50 50 H 1.09001 * 109.52624 * 176.84469 * 27 5 4 51 51 H 1.09009 * 109.51992 * 297.02880 * 27 5 4 52 52 H 1.09003 * 109.48244 * 182.49955 * 28 27 5 53 53 H 1.08998 * 109.41167 * 62.51459 * 28 27 5 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.8935 1.0276 0.0000 4 6 2.0405 -0.7221 1.2496 5 6 3.5711 -0.7704 1.2124 6 1 3.9688 0.2440 1.2426 7 6 4.0801 -1.5470 2.3993 8 8 4.2980 -2.7358 2.2980 9 7 4.2930 -0.9206 3.5734 10 6 4.1639 0.5255 3.8088 11 6 4.4792 0.7783 5.2996 12 6 4.2341 -0.6038 5.9421 13 6 4.6950 -1.5722 4.8350 14 6 5.0867 -0.7714 7.1906 15 6 4.8760 -2.1765 7.7676 16 7 3.4356 -2.3965 7.9569 17 6 2.9438 -2.7752 9.1775 18 1 3.3559 -2.3507 10.0810 19 6 1.9157 -3.7048 9.2535 20 7 1.4544 -4.2671 8.1565 21 14 1.1867 -4.1575 10.9127 22 1 0.0534 -3.2485 11.2202 23 1 0.6997 -5.5599 10.8734 24 1 2.2273 -4.0242 11.9636 25 6 2.5317 -2.2032 6.8147 26 6 2.7558 -0.7980 6.2433 27 6 4.0198 -1.4583 -0.0806 28 6 3.4309 -0.7089 -1.2789 29 8 2.0068 -0.6741 -1.1663 30 1 -0.3633 0.5134 0.8903 31 1 -0.3634 0.5142 -0.8897 32 1 -0.3633 -1.0276 -0.0005 33 1 1.6445 -1.7374 1.2710 34 1 1.7147 -0.1852 2.1404 35 1 3.1470 0.8481 3.5852 36 1 4.8728 1.0677 3.1830 37 1 3.8015 1.5233 5.7165 38 1 5.5166 1.0862 5.4307 39 1 5.7771 -1.6994 4.8684 40 1 4.1954 -2.5352 4.9409 41 1 6.1375 -0.6379 6.9337 42 1 4.7961 -0.0269 7.9317 43 1 5.3875 -2.2583 8.7267 44 1 5.2721 -2.9211 7.0772 45 1 0.7363 -4.9170 8.2096 46 1 1.4979 -2.3028 7.1456 47 1 2.7434 -2.9479 6.0475 48 1 2.4324 -0.0538 6.9711 49 1 2.1797 -0.6828 5.3252 50 1 5.1081 -1.4465 -0.1409 51 1 3.6649 -2.4890 -0.0870 52 1 3.7104 -1.2199 -2.2003 53 1 3.8193 0.3094 -1.2958 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 ZINC001040895803.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:00:44 Heat of formation + Delta-G solvation = -75.627269 kcal Electronic energy + Delta-G solvation = -31318.936171 eV Core-core repulsion = 27383.150520 eV Total energy + Delta-G solvation = -3935.785652 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.11 seconds Orbital eigenvalues (eV) -40.70752 -39.56807 -38.14118 -35.82171 -33.88553 -33.00020 -32.94803 -31.06270 -30.02087 -29.32383 -28.30581 -26.73476 -25.22414 -24.29262 -23.00748 -22.47876 -22.23114 -20.74748 -20.70331 -18.94451 -18.19684 -17.14200 -16.58712 -16.43753 -16.16142 -15.97012 -15.17885 -15.00507 -14.88631 -14.38992 -14.12615 -13.78799 -13.59731 -13.55386 -13.20217 -13.15768 -12.91514 -12.76380 -12.58745 -12.55488 -12.29432 -12.09889 -11.91003 -11.68422 -11.58972 -11.33621 -11.25799 -11.14344 -10.95134 -10.85312 -10.79027 -10.64953 -10.51083 -10.21597 -10.16311 -9.99341 -9.83951 -9.53893 -9.43901 -9.35815 -8.48718 -8.43722 -6.46844 -5.78065 -3.34349 2.69911 3.06520 3.36531 3.43485 3.44918 3.53420 4.32621 4.45690 4.48408 4.52034 4.65453 4.66695 4.83350 4.90452 4.93963 5.03296 5.09525 5.22140 5.23225 5.29999 5.38134 5.40348 5.43799 5.44073 5.48511 5.58012 5.66515 5.67848 5.74509 5.76323 5.89256 5.92286 5.95937 5.98460 6.07089 6.09306 6.11495 6.25768 6.32886 6.37476 6.48317 6.63488 6.66661 6.76766 6.88229 7.00504 7.01026 7.37915 7.65804 7.70892 7.81185 8.28434 8.54590 9.25731 9.83758 10.46966 11.29509 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.015647 B = 0.003014 C = 0.002784 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1789.005683 B = 9288.958024 C =10054.383658 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.146 4.146 2 C 0.091 3.909 3 H 0.058 0.942 4 C -0.139 4.139 5 C -0.103 4.103 6 H 0.087 0.913 7 C 0.531 3.469 8 O -0.536 6.536 9 N -0.614 5.614 10 C 0.103 3.897 11 C -0.127 4.127 12 C -0.055 4.055 13 C 0.115 3.885 14 C -0.117 4.117 15 C 0.144 3.856 16 N -0.590 5.590 17 C -0.260 4.260 18 H 0.060 0.940 19 C 0.013 3.987 20 N -0.900 5.900 21 Si 0.902 3.098 22 H -0.287 1.287 23 H -0.290 1.290 24 H -0.278 1.278 25 C 0.175 3.825 26 C -0.125 4.125 27 C -0.140 4.140 28 C 0.065 3.935 29 O -0.379 6.379 30 H 0.073 0.927 31 H 0.063 0.937 32 H 0.066 0.934 33 H 0.081 0.919 34 H 0.079 0.921 35 H 0.065 0.935 36 H 0.062 0.938 37 H 0.080 0.920 38 H 0.072 0.928 39 H 0.061 0.939 40 H 0.089 0.911 41 H 0.051 0.949 42 H 0.063 0.937 43 H 0.054 0.946 44 H 0.018 0.982 45 H 0.255 0.745 46 H 0.081 0.919 47 H 0.020 0.980 48 H 0.057 0.943 49 H 0.046 0.954 50 H 0.075 0.925 51 H 0.086 0.914 52 H 0.097 0.903 53 H 0.049 0.951 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 4.024 14.224 -15.925 21.728 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.203 4.203 2 C 0.033 3.967 3 H 0.076 0.924 4 C -0.177 4.177 5 C -0.124 4.124 6 H 0.106 0.894 7 C 0.319 3.681 8 O -0.413 6.413 9 N -0.348 5.348 10 C -0.020 4.020 11 C -0.165 4.165 12 C -0.056 4.056 13 C -0.007 4.007 14 C -0.157 4.157 15 C 0.017 3.983 16 N -0.313 5.313 17 C -0.390 4.390 18 H 0.078 0.922 19 C -0.187 4.187 20 N -0.543 5.543 21 Si 0.731 3.269 22 H -0.214 1.214 23 H -0.216 1.216 24 H -0.203 1.203 25 C 0.050 3.950 26 C -0.164 4.164 27 C -0.178 4.178 28 C -0.012 4.012 29 O -0.299 6.299 30 H 0.091 0.909 31 H 0.082 0.918 32 H 0.085 0.915 33 H 0.100 0.900 34 H 0.098 0.902 35 H 0.083 0.917 36 H 0.080 0.920 37 H 0.099 0.901 38 H 0.091 0.909 39 H 0.079 0.921 40 H 0.107 0.893 41 H 0.070 0.930 42 H 0.081 0.919 43 H 0.072 0.928 44 H 0.036 0.964 45 H 0.065 0.935 46 H 0.099 0.901 47 H 0.038 0.962 48 H 0.075 0.925 49 H 0.064 0.936 50 H 0.094 0.906 51 H 0.105 0.895 52 H 0.115 0.885 53 H 0.067 0.933 Dipole moment (debyes) X Y Z Total from point charges 3.866 13.308 -16.197 21.317 hybrid contribution -0.141 -0.074 1.111 1.123 sum 3.726 13.234 -15.086 20.411 Atomic orbital electron populations 1.21773 0.92729 1.03331 1.02490 1.22068 0.95315 0.94665 0.84659 0.92361 1.22178 0.93017 1.03361 0.99113 1.21313 0.97281 1.00034 0.93799 0.89436 1.20122 0.77585 0.87508 0.82848 1.90717 1.51775 1.15447 1.83358 1.48589 1.67997 1.08373 1.09809 1.22049 1.00836 0.81490 0.97625 1.22724 1.02035 0.98528 0.93175 1.21000 0.97823 0.92551 0.94251 1.22377 0.99155 0.98022 0.81106 1.21857 0.97007 1.00953 0.95854 1.20971 0.87374 0.92015 0.97932 1.43898 0.99573 1.79573 1.08236 1.19086 1.17446 1.17783 0.84671 0.92212 1.25402 1.00531 0.96030 0.96719 1.69940 1.28566 1.29551 1.26212 0.86316 0.79062 0.77599 0.83883 1.21383 1.21596 1.20344 1.20761 0.96511 0.89139 0.88635 1.22019 0.92603 1.02378 0.99402 1.22088 1.01343 1.01554 0.92826 1.22565 0.81010 0.99345 0.98294 1.87941 1.20953 1.78056 1.42916 0.90853 0.91789 0.91537 0.89985 0.90226 0.91655 0.91988 0.90132 0.90895 0.92119 0.89296 0.92986 0.91868 0.92759 0.96424 0.93512 0.90094 0.96245 0.92451 0.93559 0.90606 0.89500 0.88515 0.93298 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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 -1.22 9.57 37.16 0.36 -0.87 16 2 C 0.09 0.85 3.26 -26.68 -0.09 0.76 16 3 H 0.06 0.47 8.14 -51.93 -0.42 0.04 16 4 C -0.14 -1.50 4.69 -26.59 -0.12 -1.63 16 5 C -0.10 -1.16 2.18 -91.59 -0.20 -1.36 16 6 H 0.09 0.79 7.25 -51.93 -0.38 0.42 16 7 C 0.53 8.18 6.95 -10.99 -0.08 8.10 16 8 O -0.54 -10.58 16.37 5.56 0.09 -10.49 16 9 N -0.61 -8.89 3.32 -170.75 -0.57 -9.46 16 10 C 0.10 1.12 5.72 -2.96 -0.02 1.10 16 11 C -0.13 -1.46 5.89 -25.25 -0.15 -1.61 16 12 C -0.06 -0.86 0.54 -154.10 -0.08 -0.94 16 13 C 0.11 1.92 4.87 -2.86 -0.01 1.90 16 14 C -0.12 -1.91 4.83 -27.00 -0.13 -2.04 16 15 C 0.14 2.97 6.24 -3.63 -0.02 2.94 16 16 N -0.59 -14.61 2.94 -172.13 -0.51 -15.11 16 17 C -0.26 -7.09 9.49 -15.06 -0.14 -7.23 16 18 H 0.06 1.61 7.83 -52.49 -0.41 1.20 16 19 C 0.01 0.36 5.17 -17.43 -0.09 0.27 16 20 N -0.90 -26.31 11.21 55.49 0.62 -25.68 16 21 Si 0.90 21.64 29.55 -169.99 -5.02 16.62 16 22 H -0.29 -6.83 7.11 56.52 0.40 -6.43 16 23 H -0.29 -7.05 7.11 56.52 0.40 -6.65 16 24 H -0.28 -6.53 7.11 56.52 0.40 -6.13 16 25 C 0.18 4.14 4.90 -3.63 -0.02 4.13 16 26 C -0.12 -2.27 4.83 -27.01 -0.13 -2.40 16 27 C -0.14 -1.51 5.04 -26.58 -0.13 -1.64 16 28 C 0.06 0.57 6.84 37.21 0.25 0.83 16 29 O -0.38 -4.15 10.41 -35.23 -0.37 -4.52 16 30 H 0.07 0.52 8.14 -51.93 -0.42 0.10 16 31 H 0.06 0.49 8.14 -51.93 -0.42 0.07 16 32 H 0.07 0.61 8.14 -51.93 -0.42 0.19 16 33 H 0.08 1.07 8.14 -51.93 -0.42 0.64 16 34 H 0.08 0.81 7.70 -51.93 -0.40 0.41 16 35 H 0.07 0.64 7.70 -51.93 -0.40 0.24 16 36 H 0.06 0.55 7.98 -51.93 -0.41 0.13 16 37 H 0.08 0.86 8.01 -51.93 -0.42 0.44 16 38 H 0.07 0.76 7.99 -51.93 -0.42 0.34 16 39 H 0.06 0.98 7.97 -51.92 -0.41 0.57 16 40 H 0.09 1.83 6.15 -51.93 -0.32 1.52 16 41 H 0.05 0.73 8.09 -51.93 -0.42 0.31 16 42 H 0.06 1.04 8.14 -51.93 -0.42 0.62 16 43 H 0.05 1.11 7.93 -51.93 -0.41 0.70 16 44 H 0.02 0.38 7.64 -51.93 -0.40 -0.02 16 45 H 0.25 7.23 8.31 -40.82 -0.34 6.89 16 46 H 0.08 2.17 6.43 -51.93 -0.33 1.83 16 47 H 0.02 0.51 6.58 -51.93 -0.34 0.16 16 48 H 0.06 1.05 8.14 -51.93 -0.42 0.63 16 49 H 0.05 0.80 7.84 -51.93 -0.41 0.39 16 50 H 0.08 0.76 8.14 -51.93 -0.42 0.34 16 51 H 0.09 1.13 7.76 -51.92 -0.40 0.73 16 52 H 0.10 0.73 8.14 -51.93 -0.42 0.30 16 53 H 0.05 0.36 8.14 -51.93 -0.42 -0.07 16 LS Contribution 396.74 15.07 5.98 5.98 Total: -1.00 -32.20 396.74 -10.22 -42.43 By element: Atomic # 1 Polarization: 9.55 SS G_CDS: -9.64 Total: -0.09 kcal Atomic # 6 Polarization: 1.13 SS G_CDS: -0.81 Total: 0.32 kcal Atomic # 7 Polarization: -49.80 SS G_CDS: -0.45 Total: -50.25 kcal Atomic # 8 Polarization: -14.73 SS G_CDS: -0.28 Total: -15.00 kcal Atomic # 14 Polarization: 21.64 SS G_CDS: -5.02 Total: 16.62 kcal Total LS contribution 5.98 Total: 5.98 kcal Total: -32.20 -10.22 -42.43 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. ZINC001040895803.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 -33.201 kcal (2) G-P(sol) polarization free energy of solvation -32.203 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -65.404 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.223 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.427 kcal (6) G-S(sol) free energy of system = (1) + (5) -75.627 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.11 seconds