Wall clock time and date at job start Tue Mar 31 2020 20:26:03 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.52996 * 1 3 3 C 1.52995 * 109.47364 * 2 1 4 4 H 1.09007 * 109.47338 * 54.57008 * 3 2 1 5 5 C 1.53004 * 109.47045 * 174.56786 * 3 2 1 6 6 N 1.46496 * 109.47396 * 294.56759 * 3 2 1 7 7 C 1.46947 * 124.35593 * 304.42867 * 6 3 2 8 8 C 1.53681 * 105.33129 * 196.25343 * 7 6 3 9 9 C 1.53845 * 102.59640 * 335.42523 * 8 7 6 10 10 C 1.34527 * 124.35341 * 124.36910 * 6 3 2 11 11 O 1.21284 * 124.94527 * 359.68694 * 10 6 3 12 12 C 1.51861 * 110.78905 * 217.39782 * 8 7 6 13 13 C 1.53204 * 109.44195 * 184.15155 * 12 8 7 14 14 C 1.53057 * 109.46579 * 59.72638 * 13 12 8 15 15 C 1.53187 * 109.48279 * 298.05306 * 14 13 12 16 16 N 1.45966 * 110.84496 * 93.65429 * 8 7 6 17 17 C 1.34771 * 120.64190 * 354.79388 * 16 8 7 18 18 O 1.21278 * 120.00184 * 351.81717 * 17 16 8 19 19 C 1.50704 * 120.00245 * 171.54007 * 17 16 8 20 20 C 1.52996 * 109.47637 * 180.27686 * 19 17 16 21 21 C 1.50695 * 109.47462 * 180.02562 * 20 19 17 22 22 H 1.08000 * 120.00048 * 0.02562 * 21 20 19 23 23 C 1.37877 * 120.00293 * 180.02562 * 21 20 19 24 24 N 1.31516 * 119.99590 * 359.97438 * 23 21 20 25 25 Si 1.86795 * 120.00066 * 180.02562 * 23 21 20 26 26 H 1.48501 * 109.47413 * 359.97438 * 25 23 21 27 27 H 1.48507 * 109.47156 * 119.99816 * 25 23 21 28 28 H 1.48500 * 109.47010 * 239.99743 * 25 23 21 29 29 H 1.08997 * 109.47400 * 175.14971 * 1 2 3 30 30 H 1.09000 * 109.47493 * 295.15676 * 1 2 3 31 31 H 1.09009 * 109.47168 * 55.14487 * 1 2 3 32 32 H 1.09003 * 109.47259 * 239.99754 * 2 1 3 33 33 H 1.09000 * 109.47039 * 120.00242 * 2 1 3 34 34 H 1.09003 * 109.46834 * 180.02562 * 5 3 2 35 35 H 1.08999 * 109.47349 * 300.00080 * 5 3 2 36 36 H 1.09000 * 109.46928 * 60.00598 * 5 3 2 37 37 H 1.08995 * 110.27468 * 77.28472 * 7 6 3 38 38 H 1.09000 * 110.26988 * 315.21675 * 7 6 3 39 39 H 1.09003 * 110.51359 * 142.63927 * 9 8 7 40 40 H 1.09004 * 110.51979 * 265.35138 * 9 8 7 41 41 H 1.09003 * 109.56323 * 64.09299 * 12 8 7 42 42 H 1.09002 * 109.35907 * 304.08430 * 12 8 7 43 43 H 1.09004 * 109.47491 * 299.73481 * 13 12 8 44 44 H 1.09000 * 109.47540 * 179.72970 * 13 12 8 45 45 H 1.08994 * 109.46745 * 178.04001 * 14 13 12 46 46 H 1.09003 * 109.46200 * 58.05157 * 14 13 12 47 47 H 1.09003 * 109.61227 * 295.52648 * 15 14 13 48 48 H 1.09002 * 109.75441 * 175.10190 * 15 14 13 49 49 H 1.09002 * 109.47038 * 300.27953 * 19 17 16 50 50 H 1.09005 * 109.46757 * 60.27335 * 19 17 16 51 51 H 1.09002 * 109.47124 * 300.00214 * 20 19 17 52 52 H 1.09005 * 109.46918 * 59.99662 * 20 19 17 53 53 H 0.97003 * 119.99433 * 180.02562 * 24 23 21 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 6 2.0400 1.4424 0.0000 4 1 1.5995 1.9836 0.8374 5 6 3.5639 1.4445 0.1366 6 7 1.6614 2.0943 -1.2561 7 6 1.9999 1.5820 -2.5912 8 6 1.7774 2.7636 -3.5483 9 6 0.7557 3.6250 -2.7861 10 6 0.9582 3.2384 -1.3358 11 8 0.5464 3.8684 -0.3847 12 6 1.1948 2.2867 -4.8671 13 6 1.0825 3.4707 -5.8328 14 6 2.4709 4.0702 -6.0685 15 6 3.0321 4.5958 -4.7436 16 7 3.0190 3.5014 -3.7601 17 6 4.1332 3.1827 -3.0721 18 8 4.0771 2.3793 -2.1654 19 6 5.4501 3.8201 -3.4337 20 6 6.5464 3.2846 -2.5105 21 6 7.8632 3.9224 -2.8716 22 1 7.9130 4.6384 -3.6786 23 6 9.0042 3.5919 -2.1719 24 7 8.9434 2.7203 -1.1889 25 14 10.6362 4.3830 -2.6191 26 1 10.4371 5.3222 -3.7521 27 1 11.1640 5.1260 -1.4465 28 1 11.6051 3.3294 -3.0145 29 1 -0.3634 -1.0239 -0.0869 30 1 -0.3634 0.4368 0.9302 31 1 -0.3634 0.5874 -0.8434 32 1 1.8933 -0.5139 -0.8900 33 1 1.8933 -0.5139 0.8900 34 1 3.9272 2.4722 0.1370 35 1 3.8457 0.9595 1.0712 36 1 4.0045 0.9034 -0.7008 37 1 1.3432 0.7531 -2.8553 38 1 3.0419 1.2633 -2.6203 39 1 0.9652 4.6847 -2.9318 40 1 -0.2588 3.3877 -3.1066 41 1 0.2058 1.8613 -4.6967 42 1 1.8481 1.5274 -5.2971 43 1 0.4276 4.2289 -5.4036 44 1 0.6689 3.1280 -6.7813 45 1 2.3961 4.8904 -6.7824 46 1 3.1354 3.3026 -6.4654 47 1 2.4128 5.4178 -4.3844 48 1 4.0545 4.9447 -4.8890 49 1 5.6975 3.5813 -4.4681 50 1 5.3742 4.9014 -3.3185 51 1 6.2991 3.5233 -1.4761 52 1 6.6223 2.2033 -2.6257 53 1 9.7462 2.4880 -0.6963 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 ZINC001269990418.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 20:26:03 Heat of formation + Delta-G solvation = -57.906300 kcal Electronic energy + Delta-G solvation = -32305.191507 eV Core-core repulsion = 28370.174294 eV Total energy + Delta-G solvation = -3935.017213 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.83 seconds Orbital eigenvalues (eV) -41.54173 -39.05792 -37.54501 -36.45044 -34.71508 -33.46748 -32.94117 -31.44050 -30.81277 -29.63699 -28.36101 -26.82595 -25.62354 -24.44489 -23.54360 -22.54912 -21.64818 -21.21064 -19.99026 -18.83075 -18.55800 -17.04133 -16.64292 -16.23088 -15.97634 -15.43309 -15.11450 -14.90954 -14.51671 -14.44252 -14.13799 -14.00302 -13.90757 -13.78247 -13.44049 -13.20451 -13.06740 -12.87414 -12.61880 -12.42278 -12.18527 -12.06774 -11.88299 -11.69336 -11.53258 -11.41449 -11.20071 -11.07735 -10.99132 -10.96825 -10.80921 -10.76815 -10.64939 -10.59269 -10.40975 -9.99370 -9.60119 -9.57918 -9.12760 -8.80460 -8.54747 -8.35007 -7.96447 -5.99762 -4.03652 2.80742 3.03733 3.08241 3.30028 3.30898 3.37846 3.40481 4.22887 4.36086 4.44676 4.56157 4.70441 4.82043 4.89709 4.91260 5.03193 5.12817 5.18617 5.24061 5.30270 5.35957 5.37144 5.45899 5.49108 5.60108 5.62613 5.68419 5.74792 5.76922 5.79083 5.80629 5.83221 5.87653 5.96529 5.99994 6.08446 6.12254 6.15011 6.18638 6.25919 6.32660 6.53228 6.75654 6.84135 6.94327 7.03595 7.48640 7.65604 7.70570 7.82036 8.02518 8.35552 8.44969 8.70387 9.04463 9.60279 11.10989 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.011136 B = 0.004145 C = 0.003282 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2513.870454 B = 6754.154473 C = 8529.391693 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.147 4.147 2 C -0.132 4.132 3 C 0.155 3.845 4 H 0.084 0.916 5 C -0.144 4.144 6 N -0.580 5.580 7 C 0.072 3.928 8 C 0.166 3.834 9 C -0.164 4.164 10 C 0.515 3.485 11 O -0.549 6.549 12 C -0.114 4.114 13 C -0.118 4.118 14 C -0.135 4.135 15 C 0.108 3.892 16 N -0.611 5.611 17 C 0.509 3.491 18 O -0.504 6.504 19 C -0.129 4.129 20 C 0.041 3.959 21 C -0.527 4.527 22 H 0.055 0.945 23 C 0.060 3.940 24 N -0.888 5.888 25 Si 0.890 3.110 26 H -0.273 1.273 27 H -0.284 1.284 28 H -0.284 1.284 29 H 0.050 0.950 30 H 0.051 0.949 31 H 0.051 0.949 32 H 0.066 0.934 33 H 0.072 0.928 34 H 0.063 0.937 35 H 0.049 0.951 36 H 0.076 0.924 37 H 0.055 0.945 38 H 0.150 0.850 39 H 0.099 0.901 40 H 0.104 0.896 41 H 0.076 0.924 42 H 0.074 0.926 43 H 0.065 0.935 44 H 0.066 0.934 45 H 0.074 0.926 46 H 0.070 0.930 47 H 0.064 0.936 48 H 0.095 0.905 49 H 0.078 0.922 50 H 0.076 0.924 51 H 0.023 0.977 52 H 0.023 0.977 53 H 0.260 0.740 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -18.716 -2.686 -8.806 20.858 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.205 4.205 2 C -0.170 4.170 3 C 0.053 3.947 4 H 0.102 0.898 5 C -0.202 4.202 6 N -0.312 5.312 7 C -0.050 4.050 8 C 0.079 3.921 9 C -0.205 4.205 10 C 0.303 3.697 11 O -0.428 6.428 12 C -0.153 4.153 13 C -0.156 4.156 14 C -0.174 4.174 15 C -0.015 4.015 16 N -0.346 5.346 17 C 0.301 3.699 18 O -0.381 6.381 19 C -0.169 4.169 20 C 0.004 3.996 21 C -0.565 4.565 22 H 0.073 0.927 23 C -0.142 4.142 24 N -0.527 5.527 25 Si 0.718 3.282 26 H -0.198 1.198 27 H -0.210 1.210 28 H -0.210 1.210 29 H 0.070 0.930 30 H 0.070 0.930 31 H 0.070 0.930 32 H 0.085 0.915 33 H 0.091 0.909 34 H 0.082 0.918 35 H 0.068 0.932 36 H 0.095 0.905 37 H 0.073 0.927 38 H 0.167 0.833 39 H 0.118 0.882 40 H 0.122 0.878 41 H 0.094 0.906 42 H 0.092 0.908 43 H 0.083 0.917 44 H 0.085 0.915 45 H 0.093 0.907 46 H 0.089 0.911 47 H 0.083 0.917 48 H 0.113 0.887 49 H 0.096 0.904 50 H 0.094 0.906 51 H 0.041 0.959 52 H 0.041 0.959 53 H 0.070 0.930 Dipole moment (debyes) X Y Z Total from point charges -18.789 -2.926 -7.594 20.475 hybrid contribution 0.876 0.574 -0.364 1.108 sum -17.913 -2.352 -7.957 19.742 Atomic orbital electron populations 1.21722 0.95974 1.00958 1.01814 1.21878 0.95907 0.96565 1.02659 1.20971 0.96419 0.92291 0.85059 0.89763 1.22026 0.93166 1.03664 1.01358 1.48045 1.52598 1.24043 1.06490 1.23113 1.07041 0.94434 0.80370 1.20489 0.84878 0.91859 0.94842 1.22738 1.00893 1.02701 0.94160 1.20772 0.77438 0.80629 0.90871 1.90679 1.53167 1.56056 1.42878 1.21633 1.01181 0.98529 0.93955 1.21535 0.97247 0.97590 0.99207 1.21913 0.96942 1.01340 0.97166 1.21770 1.03019 0.88899 0.87812 1.47948 1.09556 1.35069 1.42032 1.21686 0.85308 0.82083 0.80795 1.90524 1.85534 1.33670 1.28375 1.21916 0.91752 1.00532 1.02683 1.18869 0.90053 0.96671 0.94040 1.21241 0.93050 1.24973 1.17237 0.92717 1.24994 0.99036 0.95147 0.95036 1.69998 1.28636 1.32839 1.21232 0.86576 0.84060 0.79237 0.78331 1.19776 1.21027 1.21033 0.93050 0.92960 0.92997 0.91512 0.90912 0.91797 0.93217 0.90504 0.92739 0.83343 0.88229 0.87805 0.90561 0.90763 0.91656 0.91503 0.90745 0.91122 0.91744 0.88719 0.90354 0.90559 0.95888 0.95894 0.93043 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.15 -1.61 9.86 37.16 0.37 -1.25 16 2 C -0.13 -1.57 5.36 -26.74 -0.14 -1.71 16 3 C 0.15 2.31 2.88 -67.93 -0.20 2.11 16 4 H 0.08 1.37 7.54 -51.92 -0.39 0.98 16 5 C -0.14 -2.44 7.75 37.16 0.29 -2.15 16 6 N -0.58 -8.63 1.62 -165.09 -0.27 -8.90 16 7 C 0.07 0.94 4.67 -3.44 -0.02 0.92 16 8 C 0.17 1.87 0.87 -131.89 -0.11 1.76 16 9 C -0.16 -1.95 5.40 -27.07 -0.15 -2.10 16 10 C 0.52 8.34 8.04 -10.73 -0.09 8.25 16 11 O -0.55 -10.96 17.36 5.55 0.10 -10.86 16 12 C -0.11 -0.94 4.65 -27.20 -0.13 -1.07 16 13 C -0.12 -0.83 5.79 -26.60 -0.15 -0.98 16 14 C -0.14 -1.05 6.05 -26.61 -0.16 -1.21 16 15 C 0.11 1.08 6.10 -3.59 -0.02 1.06 16 16 N -0.61 -7.85 2.69 -158.33 -0.43 -8.28 16 17 C 0.51 8.64 7.47 -10.99 -0.08 8.56 16 18 O -0.50 -9.56 7.11 5.56 0.04 -9.52 16 19 C -0.13 -2.43 4.01 -27.88 -0.11 -2.54 16 20 C 0.04 1.04 4.26 -27.89 -0.12 0.92 16 21 C -0.53 -14.44 9.11 -34.44 -0.31 -14.75 16 22 H 0.05 1.43 7.74 -52.49 -0.41 1.02 16 23 C 0.06 1.69 5.46 -17.82 -0.10 1.59 16 24 N -0.89 -26.63 13.29 55.43 0.74 -25.90 16 25 Si 0.89 21.30 29.54 -169.99 -5.02 16.28 16 26 H -0.27 -6.33 7.11 56.52 0.40 -5.93 16 27 H -0.28 -6.79 7.11 56.52 0.40 -6.39 16 28 H -0.28 -6.74 7.11 56.52 0.40 -6.34 16 29 H 0.05 0.48 8.14 -51.93 -0.42 0.06 16 30 H 0.05 0.59 8.14 -51.93 -0.42 0.16 16 31 H 0.05 0.61 7.64 -51.92 -0.40 0.21 16 32 H 0.07 0.78 7.63 -51.93 -0.40 0.39 16 33 H 0.07 0.80 8.14 -51.93 -0.42 0.38 16 34 H 0.06 1.29 7.47 -51.93 -0.39 0.90 16 35 H 0.05 0.76 8.14 -51.93 -0.42 0.34 16 36 H 0.08 1.39 6.17 -51.93 -0.32 1.07 16 37 H 0.05 0.60 7.86 -51.93 -0.41 0.19 16 38 H 0.15 2.32 4.27 -56.81 -0.24 2.08 16 39 H 0.10 1.20 7.16 -51.93 -0.37 0.83 16 40 H 0.10 0.99 7.72 -51.93 -0.40 0.59 16 41 H 0.08 0.53 7.93 -51.93 -0.41 0.11 16 42 H 0.07 0.64 8.14 -51.93 -0.42 0.21 16 43 H 0.06 0.44 7.65 -51.93 -0.40 0.05 16 44 H 0.07 0.40 8.14 -51.93 -0.42 -0.03 16 45 H 0.07 0.47 8.14 -51.93 -0.42 0.05 16 46 H 0.07 0.60 8.14 -51.93 -0.42 0.18 16 47 H 0.06 0.60 7.44 -51.93 -0.39 0.21 16 48 H 0.09 0.98 5.96 -51.93 -0.31 0.67 16 49 H 0.08 1.36 7.99 -51.93 -0.41 0.95 16 50 H 0.08 1.37 7.65 -51.93 -0.40 0.97 16 51 H 0.02 0.64 8.10 -51.93 -0.42 0.22 16 52 H 0.02 0.63 8.10 -51.93 -0.42 0.21 16 53 H 0.26 7.35 8.31 -40.82 -0.34 7.01 16 LS Contribution 396.11 15.07 5.97 5.97 Total: -1.00 -32.92 396.11 -9.50 -42.42 By element: Atomic # 1 Polarization: 10.76 SS G_CDS: -9.40 Total: 1.36 kcal Atomic # 6 Polarization: -1.35 SS G_CDS: -1.23 Total: -2.59 kcal Atomic # 7 Polarization: -43.11 SS G_CDS: 0.04 Total: -43.07 kcal Atomic # 8 Polarization: -20.51 SS G_CDS: 0.14 Total: -20.38 kcal Atomic # 14 Polarization: 21.30 SS G_CDS: -5.02 Total: 16.28 kcal Total LS contribution 5.97 Total: 5.97 kcal Total: -32.92 -9.50 -42.42 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. ZINC001269990418.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 -15.481 kcal (2) G-P(sol) polarization free energy of solvation -32.920 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -48.402 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.505 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.425 kcal (6) G-S(sol) free energy of system = (1) + (5) -57.906 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.83 seconds