Wall clock time and date at job start Sat Apr 11 2020 06:09:18 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.52997 * 1 3 3 C 1.52999 * 109.47049 * 2 1 4 4 H 1.08992 * 109.47365 * 54.73683 * 3 2 1 5 5 O 1.42902 * 109.47275 * 294.73127 * 3 2 1 6 6 C 1.42898 * 113.99945 * 151.12957 * 5 3 2 7 7 C 1.54271 * 110.03156 * 270.07265 * 6 5 3 8 8 C 1.54955 * 104.15857 * 217.13665 * 7 6 5 9 9 C 1.54887 * 102.78287 * 37.89057 * 8 7 6 10 10 C 1.53877 * 110.03021 * 152.85613 * 6 5 3 11 11 C 1.50703 * 109.46975 * 174.73527 * 3 2 1 12 12 O 1.21284 * 120.00278 * 114.99933 * 11 3 2 13 13 N 1.34782 * 119.99824 * 294.99775 * 11 3 2 14 14 C 1.46496 * 120.00126 * 179.72363 * 13 11 3 15 15 C 1.52994 * 109.47560 * 180.02562 * 14 13 11 16 16 O 1.42903 * 110.88486 * 297.57516 * 15 14 13 17 17 C 1.55151 * 111.00850 * 173.65333 * 15 14 13 18 18 C 1.54321 * 102.94088 * 206.42286 * 17 15 14 19 19 N 1.47020 * 107.27541 * 22.19256 * 18 17 15 20 20 C 1.36935 * 125.65079 * 181.09242 * 19 18 17 21 21 H 1.07999 * 119.99932 * 179.97438 * 20 19 18 22 22 C 1.38810 * 120.00015 * 359.95798 * 20 19 18 23 23 N 1.31625 * 120.00059 * 0.02562 * 22 20 19 24 24 Si 1.86801 * 119.99938 * 179.97438 * 22 20 19 25 25 H 1.48499 * 109.47002 * 0.02562 * 24 22 20 26 26 H 1.48500 * 109.46675 * 120.00256 * 24 22 20 27 27 H 1.48495 * 109.47385 * 239.99880 * 24 22 20 28 28 C 1.47423 * 108.69914 * 1.05936 * 19 18 17 29 29 H 1.09003 * 109.47204 * 185.28901 * 1 2 3 30 30 H 1.09003 * 109.46644 * 305.28204 * 1 2 3 31 31 H 1.08995 * 109.47469 * 65.27858 * 1 2 3 32 32 H 1.09007 * 109.46904 * 240.00107 * 2 1 3 33 33 H 1.08996 * 109.47273 * 120.00671 * 2 1 3 34 34 H 1.08995 * 110.01616 * 31.46378 * 6 5 3 35 35 H 1.09001 * 110.48241 * 335.78708 * 7 6 5 36 36 H 1.08995 * 110.48265 * 98.42695 * 7 6 5 37 37 H 1.09006 * 110.72437 * 279.57386 * 8 7 6 38 38 H 1.09004 * 110.76530 * 156.22204 * 8 7 6 39 39 H 1.08996 * 110.49045 * 203.40706 * 9 8 7 40 40 H 1.09000 * 110.60066 * 80.81280 * 9 8 7 41 41 H 1.08998 * 110.03124 * 238.56673 * 10 6 5 42 42 H 1.08997 * 110.03200 * 359.95252 * 10 6 5 43 43 H 0.97001 * 119.99538 * 0.02562 * 13 11 3 44 44 H 1.08999 * 109.47057 * 300.00761 * 14 13 11 45 45 H 1.09003 * 109.47282 * 60.00071 * 14 13 11 46 46 H 0.96700 * 113.99876 * 59.29991 * 16 15 14 47 47 H 1.08999 * 110.72007 * 324.78535 * 17 15 14 48 48 H 1.09001 * 110.72166 * 88.06202 * 17 15 14 49 49 H 1.09003 * 109.88533 * 262.77186 * 18 17 15 50 50 H 1.09011 * 109.88154 * 141.61163 * 18 17 15 51 51 H 0.96997 * 119.99740 * 179.97438 * 23 22 20 52 52 H 1.09003 * 110.36576 * 217.03374 * 28 19 18 53 53 H 1.08996 * 110.37269 * 94.71383 * 28 19 18 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.0399 1.4425 0.0000 4 1 1.6017 1.9828 0.8390 5 8 1.6673 2.0795 -1.2237 6 6 1.4856 3.4924 -1.1116 7 6 0.0199 3.8096 -0.7492 8 6 -0.2961 5.1015 -1.5444 9 6 0.4524 4.8729 -2.8810 10 6 1.7607 4.1684 -2.4663 11 6 3.5413 1.4445 0.1304 12 8 4.2285 1.8598 -0.7786 13 7 4.1205 0.9847 1.2572 14 6 5.5802 0.9809 1.3816 15 6 5.9743 0.4130 2.7465 16 8 5.4720 1.2327 3.8037 17 6 7.5128 0.2443 2.8555 18 6 7.6635 -0.8938 3.8867 19 7 6.4014 -1.6479 3.8922 20 6 6.1195 -2.7264 4.6876 21 1 5.1582 -3.2134 4.6161 22 6 7.0694 -3.1937 5.5854 23 7 8.2409 -2.6000 5.6726 24 14 6.6845 -4.6644 6.6709 25 1 5.3100 -5.1466 6.3822 26 1 7.6581 -5.7514 6.3951 27 1 6.7799 -4.2649 8.0979 28 6 5.4894 -1.0502 2.9000 29 1 -0.3634 -1.0233 0.0947 30 1 -0.3633 0.5936 0.8389 31 1 -0.3634 0.4297 -0.9334 32 1 1.8933 -0.5139 -0.8901 33 1 1.8933 -0.5139 0.8899 34 1 2.1539 3.8953 -0.3507 35 1 -0.6352 2.9976 -1.0647 36 1 -0.0785 3.9862 0.3218 37 1 0.0951 5.9779 -1.0277 38 1 -1.3685 5.1983 -1.7139 39 1 0.6688 5.8253 -3.3649 40 1 -0.1316 4.2337 -3.5432 41 1 2.5620 4.9000 -2.3626 42 1 2.0337 3.4188 -3.2090 43 1 3.5709 0.6522 1.9841 44 1 6.0104 0.3632 0.5934 45 1 5.9553 2.0003 1.2898 46 1 4.5095 1.3263 3.8054 47 1 7.9379 -0.0468 1.8949 48 1 7.9762 1.1604 3.2218 49 1 7.8497 -0.4756 4.8759 50 1 8.4856 -1.5490 3.5982 51 1 8.9044 -2.9264 6.3003 52 1 4.4629 -1.0732 3.2661 53 1 5.5656 -1.5769 1.9488 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 ZINC001278811308.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Apr 11 2020 06:09:18 Heat of formation + Delta-G solvation = -122.296411 kcal Electronic energy + Delta-G solvation = -31069.566822 eV Core-core repulsion = 26941.301393 eV Total energy + Delta-G solvation = -4128.265429 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.75 seconds Orbital eigenvalues (eV) -40.41714 -39.10909 -38.25865 -36.36800 -35.13449 -34.31185 -33.37103 -32.11151 -29.76152 -29.02827 -28.76121 -28.05342 -26.87592 -25.46579 -23.55343 -22.30082 -21.17177 -20.73419 -20.58895 -20.10035 -19.85202 -17.98448 -17.46828 -16.91659 -16.30065 -16.14300 -15.61464 -15.27591 -15.07420 -14.73227 -14.48682 -14.37633 -14.18603 -13.95932 -13.76692 -13.63956 -13.32838 -12.90820 -12.80945 -12.59580 -12.29889 -12.14631 -12.03240 -11.94968 -11.73004 -11.46672 -11.40978 -11.34127 -11.16906 -11.01613 -10.94580 -10.91786 -10.77811 -10.74460 -10.49777 -10.35002 -10.11983 -10.10229 -9.78961 -9.48703 -9.27812 -9.11547 -8.61413 -6.85529 -5.73541 -3.07313 2.50626 2.93340 3.38075 3.42841 3.44613 3.45889 4.24290 4.27800 4.43695 4.46343 4.49976 4.57710 4.64280 4.67725 4.87275 4.92788 5.01993 5.11859 5.14969 5.16791 5.24579 5.32077 5.35167 5.35870 5.47461 5.50082 5.53598 5.58280 5.63538 5.68612 5.71856 5.77045 5.81866 5.88333 5.94047 5.94368 6.00124 6.02094 6.15148 6.29603 6.33265 6.41189 6.52431 6.58570 6.91294 7.38488 7.51570 7.62225 7.81588 7.95185 8.13402 8.23084 9.18628 9.78813 10.43912 11.38217 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.020720 B = 0.002293 C = 0.002156 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1351.027518 B =12206.488652 C =12983.084842 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.117 4.117 3 C 0.076 3.924 4 H 0.085 0.915 5 O -0.345 6.345 6 C 0.085 3.915 7 C -0.149 4.149 8 C -0.123 4.123 9 C -0.123 4.123 10 C -0.119 4.119 11 C 0.487 3.513 12 O -0.525 6.525 13 N -0.725 5.725 14 C 0.104 3.896 15 C 0.102 3.898 16 O -0.536 6.536 17 C -0.128 4.128 18 C 0.181 3.819 19 N -0.598 5.598 20 C -0.230 4.230 21 H 0.077 0.923 22 C -0.008 4.008 23 N -0.919 5.919 24 Si 0.913 3.087 25 H -0.280 1.280 26 H -0.290 1.290 27 H -0.290 1.290 28 C 0.118 3.882 29 H 0.059 0.941 30 H 0.054 0.946 31 H 0.062 0.938 32 H 0.079 0.921 33 H 0.082 0.918 34 H 0.061 0.939 35 H 0.079 0.921 36 H 0.065 0.935 37 H 0.067 0.933 38 H 0.068 0.932 39 H 0.066 0.934 40 H 0.067 0.933 41 H 0.071 0.929 42 H 0.080 0.920 43 H 0.398 0.602 44 H 0.078 0.922 45 H 0.073 0.927 46 H 0.362 0.638 47 H 0.072 0.928 48 H 0.072 0.928 49 H 0.051 0.949 50 H 0.044 0.956 51 H 0.252 0.748 52 H 0.037 0.963 53 H 0.027 0.973 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -17.176 13.507 -13.442 25.654 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.155 4.155 3 C 0.014 3.986 4 H 0.102 0.898 5 O -0.262 6.262 6 C 0.027 3.973 7 C -0.187 4.187 8 C -0.161 4.161 9 C -0.161 4.161 10 C -0.157 4.157 11 C 0.273 3.727 12 O -0.401 6.401 13 N -0.380 5.380 14 C -0.020 4.020 15 C 0.060 3.940 16 O -0.340 6.340 17 C -0.167 4.167 18 C 0.057 3.943 19 N -0.325 5.325 20 C -0.358 4.358 21 H 0.095 0.905 22 C -0.204 4.204 23 N -0.567 5.567 24 Si 0.743 3.257 25 H -0.205 1.205 26 H -0.217 1.217 27 H -0.217 1.217 28 C -0.009 4.009 29 H 0.079 0.921 30 H 0.073 0.927 31 H 0.081 0.919 32 H 0.097 0.903 33 H 0.101 0.899 34 H 0.079 0.921 35 H 0.097 0.903 36 H 0.084 0.916 37 H 0.086 0.914 38 H 0.087 0.913 39 H 0.085 0.915 40 H 0.086 0.914 41 H 0.090 0.910 42 H 0.099 0.901 43 H 0.232 0.768 44 H 0.096 0.904 45 H 0.091 0.909 46 H 0.208 0.792 47 H 0.091 0.909 48 H 0.091 0.909 49 H 0.069 0.931 50 H 0.062 0.938 51 H 0.063 0.937 52 H 0.055 0.945 53 H 0.044 0.956 Dipole moment (debyes) X Y Z Total from point charges -15.930 13.625 -13.981 25.197 hybrid contribution -0.825 -0.868 1.124 1.642 sum -16.755 12.758 -12.857 24.674 Atomic orbital electron populations 1.21719 0.94362 1.02052 1.02348 1.21598 0.97350 0.92689 1.03892 1.22119 0.90548 0.96489 0.89403 0.89771 1.88063 1.87860 1.20045 1.30218 1.22852 0.95155 0.81088 0.98159 1.22970 0.97020 0.98609 1.00095 1.22440 1.00286 0.97327 0.96040 1.22350 0.96180 1.00292 0.97276 1.22155 0.97707 0.99406 0.96424 1.21429 0.91348 0.77758 0.82214 1.90711 1.64814 1.45659 1.38873 1.46067 1.08213 1.65077 1.18598 1.22036 0.79540 1.02140 0.98276 1.22721 0.96857 0.90559 0.83910 1.86649 1.24989 1.70770 1.51601 1.22910 0.91055 1.01359 1.01367 1.21142 0.88626 0.89102 0.95448 1.44352 1.11820 1.33741 1.42539 1.19060 0.99395 1.03905 1.13418 0.90466 1.24816 0.96240 0.99595 0.99750 1.69973 1.05665 1.41880 1.39149 0.86036 0.77378 0.81985 0.80274 1.20541 1.21737 1.21723 1.22124 0.94371 0.90744 0.93693 0.92150 0.92719 0.91886 0.90265 0.89909 0.92068 0.90276 0.91576 0.91414 0.91345 0.91531 0.91396 0.91038 0.90141 0.76768 0.90418 0.90870 0.79185 0.90904 0.90944 0.93117 0.93762 0.93738 0.94452 0.95584 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.15 -1.07 9.76 37.16 0.36 -0.71 16 2 C -0.12 -1.03 5.61 -26.73 -0.15 -1.18 16 3 C 0.08 0.76 2.34 -27.88 -0.07 0.70 16 4 H 0.08 0.66 7.84 -51.93 -0.41 0.25 16 5 O -0.34 -3.91 8.73 -55.09 -0.48 -4.39 16 6 C 0.09 0.82 3.14 -25.61 -0.08 0.74 16 7 C -0.15 -1.11 6.55 -25.03 -0.16 -1.27 16 8 C -0.12 -0.80 6.91 -24.70 -0.17 -0.97 16 9 C -0.12 -0.88 6.83 -25.07 -0.17 -1.05 16 10 C -0.12 -1.10 6.61 -25.61 -0.17 -1.27 16 11 C 0.49 6.32 7.08 -10.98 -0.08 6.24 16 12 O -0.53 -8.95 16.51 -14.31 -0.24 -9.19 16 13 N -0.73 -8.55 5.06 -60.30 -0.31 -8.85 16 14 C 0.10 1.48 4.90 -4.05 -0.02 1.46 16 15 C 0.10 1.70 0.95 -88.48 -0.08 1.61 16 16 O -0.54 -9.01 13.19 -35.23 -0.46 -9.48 16 17 C -0.13 -2.40 6.16 -24.90 -0.15 -2.55 16 18 C 0.18 4.42 5.37 -3.07 -0.02 4.41 16 19 N -0.60 -15.03 3.39 -170.07 -0.58 -15.61 16 20 C -0.23 -6.33 10.28 -15.06 -0.15 -6.48 16 21 H 0.08 2.04 7.83 -52.49 -0.41 1.63 16 22 C -0.01 -0.21 5.49 -17.43 -0.10 -0.31 16 23 N -0.92 -27.05 10.31 55.49 0.57 -26.48 16 24 Si 0.91 21.98 29.55 -169.99 -5.02 16.95 16 25 H -0.28 -6.50 7.11 56.52 0.40 -6.10 16 26 H -0.29 -7.04 7.11 56.52 0.40 -6.64 16 27 H -0.29 -7.04 7.11 56.52 0.40 -6.64 16 28 C 0.12 2.28 5.92 -2.53 -0.01 2.26 16 29 H 0.06 0.41 8.14 -51.93 -0.42 -0.02 16 30 H 0.05 0.34 8.14 -51.93 -0.42 -0.08 16 31 H 0.06 0.49 7.97 -51.93 -0.41 0.08 16 32 H 0.08 0.80 8.14 -51.93 -0.42 0.37 16 33 H 0.08 0.65 7.76 -51.93 -0.40 0.24 16 34 H 0.06 0.63 7.83 -51.93 -0.41 0.22 16 35 H 0.08 0.60 8.08 -51.93 -0.42 0.18 16 36 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 37 H 0.07 0.42 8.14 -51.93 -0.42 -0.01 16 38 H 0.07 0.38 8.14 -51.93 -0.42 -0.04 16 39 H 0.07 0.41 8.14 -51.93 -0.42 -0.01 16 40 H 0.07 0.48 8.14 -51.93 -0.42 0.05 16 41 H 0.07 0.65 8.14 -51.93 -0.42 0.23 16 42 H 0.08 0.85 7.95 -51.93 -0.41 0.44 16 43 H 0.40 3.81 7.00 -40.82 -0.29 3.52 16 44 H 0.08 1.17 8.11 -51.93 -0.42 0.74 16 45 H 0.07 1.07 8.14 -51.93 -0.42 0.65 16 46 H 0.36 4.51 8.17 45.56 0.37 4.88 16 47 H 0.07 1.23 8.01 -51.93 -0.42 0.82 16 48 H 0.07 1.29 8.14 -51.93 -0.42 0.87 16 49 H 0.05 1.40 7.38 -51.93 -0.38 1.01 16 50 H 0.04 1.18 7.51 -51.92 -0.39 0.79 16 51 H 0.25 7.19 8.31 -40.82 -0.34 6.85 16 52 H 0.04 0.67 7.96 -51.93 -0.41 0.25 16 53 H 0.03 0.51 8.04 -51.93 -0.42 0.09 16 LS Contribution 417.25 15.07 6.29 6.29 Total: -1.00 -34.01 417.25 -10.47 -44.48 By element: Atomic # 1 Polarization: 13.67 SS G_CDS: -9.01 Total: 4.65 kcal Atomic # 6 Polarization: 2.85 SS G_CDS: -1.23 Total: 1.62 kcal Atomic # 7 Polarization: -50.62 SS G_CDS: -0.31 Total: -50.93 kcal Atomic # 8 Polarization: -21.88 SS G_CDS: -1.18 Total: -23.06 kcal Atomic # 14 Polarization: 21.98 SS G_CDS: -5.02 Total: 16.95 kcal Total LS contribution 6.29 Total: 6.29 kcal Total: -34.01 -10.47 -44.48 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. ZINC001278811308.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 -77.818 kcal (2) G-P(sol) polarization free energy of solvation -34.012 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -111.830 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.467 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.479 kcal (6) G-S(sol) free energy of system = (1) + (5) -122.296 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.75 seconds