Wall clock time and date at job start Sat Apr 11 2020 06:10:35 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.52999 * 1 3 3 H 1.09004 * 109.47340 * 2 1 4 4 C 1.53005 * 109.47140 * 119.99403 * 2 1 3 5 5 N 1.46501 * 109.46922 * 295.00363 * 4 2 1 6 6 C 1.34782 * 119.99351 * 179.97438 * 5 4 2 7 7 O 1.21508 * 120.00306 * 359.97438 * 6 5 4 8 8 O 1.34641 * 119.99271 * 179.97438 * 6 5 4 9 9 C 1.45205 * 116.99710 * 180.02562 * 8 6 5 10 10 C 1.53000 * 109.47350 * 59.99638 * 9 8 6 11 11 C 1.52994 * 109.47239 * 179.97438 * 9 8 6 12 12 C 1.53006 * 109.46749 * 300.00045 * 9 8 6 13 13 C 1.50694 * 109.47093 * 239.99506 * 2 1 3 14 14 O 1.21290 * 120.00197 * 258.18592 * 13 2 1 15 15 N 1.34775 * 120.00297 * 78.19218 * 13 2 1 16 16 C 1.37094 * 120.00725 * 184.08519 * 15 13 2 17 17 H 1.08001 * 120.00237 * 226.54778 * 16 15 13 18 18 C 1.38953 * 119.99744 * 46.53412 * 16 15 13 19 19 N 1.31414 * 119.99776 * 13.75452 * 18 16 15 20 20 Si 1.86797 * 120.00102 * 193.75372 * 18 16 15 21 21 H 1.48497 * 109.47384 * 59.99989 * 20 18 16 22 22 H 1.48500 * 109.47220 * 180.02562 * 20 18 16 23 23 H 1.48506 * 109.46979 * 300.00013 * 20 18 16 24 24 C 1.46509 * 119.99589 * 4.08228 * 15 13 2 25 25 C 1.55234 * 110.57892 * 228.56299 * 24 15 13 26 26 C 1.54893 * 102.58509 * 203.34042 * 25 24 15 27 27 C 1.54265 * 104.31383 * 37.93397 * 26 25 24 28 28 C 1.53875 * 106.59577 * 336.31995 * 27 26 25 29 29 H 1.09008 * 109.47099 * 303.95522 * 1 2 3 30 30 H 1.09006 * 109.46791 * 63.95525 * 1 2 3 31 31 H 1.08993 * 109.47569 * 183.95308 * 1 2 3 32 32 H 1.09008 * 109.46833 * 54.99806 * 4 2 1 33 33 H 1.08997 * 109.47295 * 174.99684 * 4 2 1 34 34 H 0.96991 * 120.00474 * 359.97438 * 5 4 2 35 35 H 1.09006 * 109.46859 * 60.00026 * 10 9 8 36 36 H 1.08998 * 109.47389 * 179.97438 * 10 9 8 37 37 H 1.09001 * 109.47403 * 300.00710 * 10 9 8 38 38 H 1.09003 * 109.47109 * 60.00332 * 11 9 8 39 39 H 1.09003 * 109.47126 * 179.97438 * 11 9 8 40 40 H 1.08997 * 109.47127 * 299.99918 * 11 9 8 41 41 H 1.08997 * 109.47405 * 60.00504 * 12 9 8 42 42 H 1.09006 * 109.46447 * 179.97438 * 12 9 8 43 43 H 1.08999 * 109.47067 * 299.99799 * 12 9 8 44 44 H 0.96999 * 119.99344 * 179.97438 * 19 18 16 45 45 H 1.08996 * 110.43537 * 351.14516 * 24 15 13 46 46 H 1.08999 * 110.90329 * 85.04936 * 25 24 15 47 47 H 1.09005 * 110.76891 * 321.56677 * 25 24 15 48 48 H 1.09000 * 110.48352 * 279.18250 * 26 25 24 49 49 H 1.09002 * 110.48362 * 156.67277 * 26 25 24 50 50 H 1.08995 * 110.03337 * 95.61067 * 27 26 25 51 51 H 1.09004 * 110.02707 * 216.99514 * 27 26 25 52 52 H 1.09006 * 110.02922 * 119.28837 * 28 27 26 53 53 H 1.08998 * 110.03463 * 240.66729 * 28 27 26 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.8934 1.0277 0.0000 4 6 2.0400 -0.7211 1.2494 5 7 1.6523 0.0356 2.4424 6 6 1.9955 -0.4134 3.6660 7 8 2.6260 -1.4459 3.7797 8 8 1.6397 0.2824 4.7624 9 6 2.0437 -0.2602 6.0473 10 6 3.5694 -0.3609 6.1015 11 6 1.5538 0.6607 7.1666 12 6 1.4340 -1.6519 6.2276 13 6 2.0323 -0.7105 -1.2304 14 8 2.4372 -1.8509 -1.1493 15 7 2.0318 -0.0787 -2.4208 16 6 2.4095 -0.7513 -3.5542 17 1 3.0993 -0.2947 -4.2485 18 6 1.9051 -2.0207 -3.8093 19 7 0.8754 -2.4680 -3.1263 20 14 2.6906 -3.0911 -5.1232 21 1 4.1144 -3.3336 -4.7783 22 1 1.9716 -4.3882 -5.2001 23 1 2.6098 -2.4023 -6.4364 24 6 1.6274 1.3270 -2.5039 25 6 0.5968 1.5318 -3.6465 26 6 0.7566 3.0360 -3.9796 27 6 2.2714 3.2968 -3.8489 28 6 2.8250 2.2102 -2.9105 29 1 -0.3634 0.5740 0.8525 30 1 -0.3633 0.4513 -0.9234 31 1 -0.3634 -1.0251 0.0708 32 1 1.6053 -1.7197 1.2953 33 1 3.1263 -0.7994 1.2057 34 1 1.1488 0.8596 2.3517 35 1 4.0037 0.6306 5.9729 36 1 3.8727 -0.7678 7.0661 37 1 3.9188 -1.0167 5.3040 38 1 0.4669 0.7327 7.1282 39 1 1.8567 0.2531 8.1311 40 1 1.9885 1.6519 7.0379 41 1 1.7835 -2.3078 5.4302 42 1 1.7377 -2.0594 7.1919 43 1 0.3471 -1.5802 6.1890 44 1 0.5231 -3.3540 -3.3047 45 1 1.2130 1.6606 -1.5525 46 1 -0.4146 1.3180 -3.3010 47 1 0.8495 0.9161 -4.5098 48 1 0.2015 3.6456 -3.2666 49 1 0.4209 3.2393 -4.9965 50 1 2.7484 3.2247 -4.8262 51 1 2.4446 4.2843 -3.4210 52 1 3.5709 1.6107 -3.4325 53 1 3.2678 2.6700 -2.0270 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 ZINC000495617541.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:10:35 Heat of formation + Delta-G solvation = -125.744126 kcal Electronic energy + Delta-G solvation = -31603.992597 eV Core-core repulsion = 27475.577664 eV Total energy + Delta-G solvation = -4128.414933 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.77 seconds Orbital eigenvalues (eV) -40.65883 -39.15425 -37.55139 -36.76770 -35.12261 -34.08272 -33.35383 -32.17992 -31.04494 -28.10460 -27.54236 -27.46141 -26.82173 -26.29289 -25.14719 -22.80842 -21.42339 -20.17325 -19.97388 -19.76750 -18.94794 -17.98177 -16.72925 -16.28302 -16.13126 -16.00444 -15.08508 -15.03083 -14.67515 -14.53833 -14.46352 -14.17684 -13.79683 -13.47051 -13.43055 -13.10107 -12.75929 -12.65146 -12.53658 -12.30667 -12.16717 -11.99551 -11.96752 -11.74604 -11.54994 -11.47555 -11.35764 -11.08661 -10.96033 -10.82202 -10.76729 -10.64955 -10.54432 -10.50122 -10.44969 -10.37142 -10.18434 -9.99559 -9.95735 -9.73910 -8.92102 -8.67356 -8.03017 -7.05205 -6.08676 -3.85954 2.28094 3.01135 3.31782 3.38997 3.42339 4.14462 4.22813 4.39972 4.61148 4.64726 4.79544 4.87722 5.16092 5.17390 5.21144 5.29562 5.31589 5.38788 5.40430 5.50350 5.55161 5.60830 5.63347 5.65164 5.75897 5.82580 5.84281 5.93922 6.02923 6.09600 6.12727 6.15242 6.17926 6.21446 6.23730 6.25748 6.28401 6.39068 6.39922 6.44375 6.57284 6.60493 6.64604 6.74935 6.82877 6.92673 7.16672 7.60317 7.66173 8.31936 8.54823 9.11128 9.19092 9.68590 10.31261 11.00533 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.013448 B = 0.002944 C = 0.002582 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2081.592520 B = 9508.767913 C =10842.283978 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.134 4.134 2 C -0.121 4.121 3 H 0.081 0.919 4 C 0.160 3.840 5 N -0.707 5.707 6 C 0.647 3.353 7 O -0.566 6.566 8 O -0.369 6.369 9 C 0.136 3.864 10 C -0.184 4.184 11 C -0.140 4.140 12 C -0.184 4.184 13 C 0.495 3.505 14 O -0.534 6.534 15 N -0.446 5.446 16 C -0.337 4.337 17 H 0.075 0.925 18 C 0.071 3.929 19 N -0.838 5.838 20 Si 0.909 3.091 21 H -0.275 1.275 22 H -0.287 1.287 23 H -0.279 1.279 24 C 0.127 3.873 25 C -0.118 4.118 26 C -0.125 4.125 27 C -0.118 4.118 28 C -0.135 4.135 29 H 0.047 0.953 30 H 0.064 0.936 31 H 0.064 0.936 32 H 0.081 0.919 33 H 0.073 0.927 34 H 0.402 0.598 35 H 0.058 0.942 36 H 0.068 0.932 37 H 0.089 0.911 38 H 0.065 0.935 39 H 0.078 0.922 40 H 0.065 0.935 41 H 0.090 0.910 42 H 0.068 0.932 43 H 0.059 0.941 44 H 0.263 0.737 45 H 0.072 0.928 46 H 0.056 0.944 47 H 0.088 0.912 48 H 0.056 0.944 49 H 0.057 0.943 50 H 0.058 0.942 51 H 0.053 0.947 52 H 0.074 0.926 53 H 0.051 0.949 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -1.821 11.189 12.920 17.189 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.192 4.192 2 C -0.144 4.144 3 H 0.099 0.901 4 C 0.038 3.962 5 N -0.366 5.366 6 C 0.399 3.601 7 O -0.455 6.455 8 O -0.284 6.284 9 C 0.099 3.901 10 C -0.241 4.241 11 C -0.197 4.197 12 C -0.241 4.241 13 C 0.281 3.719 14 O -0.412 6.412 15 N -0.167 5.167 16 C -0.464 4.464 17 H 0.092 0.908 18 C -0.139 4.139 19 N -0.471 5.471 20 Si 0.736 3.264 21 H -0.200 1.200 22 H -0.213 1.213 23 H -0.204 1.204 24 C 0.021 3.979 25 C -0.157 4.157 26 C -0.163 4.163 27 C -0.156 4.156 28 C -0.174 4.174 29 H 0.066 0.934 30 H 0.083 0.917 31 H 0.083 0.917 32 H 0.099 0.901 33 H 0.091 0.909 34 H 0.236 0.764 35 H 0.077 0.923 36 H 0.087 0.913 37 H 0.108 0.892 38 H 0.084 0.916 39 H 0.097 0.903 40 H 0.084 0.916 41 H 0.108 0.892 42 H 0.087 0.913 43 H 0.078 0.922 44 H 0.073 0.927 45 H 0.090 0.910 46 H 0.075 0.925 47 H 0.106 0.894 48 H 0.075 0.925 49 H 0.076 0.924 50 H 0.077 0.923 51 H 0.071 0.929 52 H 0.093 0.907 53 H 0.070 0.930 Dipole moment (debyes) X Y Z Total from point charges -1.410 9.883 13.114 16.481 hybrid contribution 0.016 -0.127 -0.663 0.676 sum -1.394 9.756 12.450 15.879 Atomic orbital electron populations 1.21697 0.93273 1.02884 1.01317 1.21504 0.97745 1.00214 0.94931 0.90071 1.21164 1.00497 0.94851 0.79714 1.44904 1.59616 1.28841 1.03213 1.18078 0.78647 0.81963 0.81385 1.90910 1.43835 1.24732 1.86064 1.86282 1.74689 1.52813 1.14639 1.22009 0.95013 0.93996 0.79053 1.22522 0.94119 1.03298 1.04172 1.21848 1.02248 0.99546 0.96034 1.22523 1.01052 0.96551 1.04008 1.19462 0.82095 0.87354 0.83004 1.90544 1.48123 1.17884 1.84657 1.43168 1.60162 1.07814 1.05570 1.18724 1.26735 1.00140 1.00773 0.90762 1.25245 0.93893 0.98326 0.96460 1.69799 1.20295 1.17837 1.39187 0.86334 0.79840 0.79532 0.80649 1.20046 1.21296 1.20438 1.22034 0.93351 0.84931 0.97554 1.22726 0.98303 0.94259 1.00386 1.22330 0.95902 0.97594 1.00502 1.22038 0.96033 0.97949 0.99560 1.22497 0.99443 0.96778 0.98687 0.93388 0.91662 0.91726 0.90107 0.90905 0.76359 0.92262 0.91334 0.89190 0.91550 0.90311 0.91557 0.89170 0.91340 0.92227 0.92669 0.91050 0.92539 0.89373 0.92524 0.92403 0.92322 0.92861 0.90710 0.92992 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.13 -2.19 8.69 37.16 0.32 -1.87 16 2 C -0.12 -2.18 1.89 -91.77 -0.17 -2.35 16 3 H 0.08 1.26 6.72 -51.93 -0.35 0.91 16 4 C 0.16 2.86 4.42 -4.04 -0.02 2.84 16 5 N -0.71 -10.35 5.32 -65.27 -0.35 -10.69 16 6 C 0.65 10.21 8.06 54.06 0.44 10.65 16 7 O -0.57 -10.46 11.88 -19.28 -0.23 -10.68 16 8 O -0.37 -4.99 9.94 -40.33 -0.40 -5.39 16 9 C 0.14 1.53 1.13 -90.62 -0.10 1.43 16 10 C -0.18 -2.05 8.37 37.16 0.31 -1.74 16 11 C -0.14 -1.12 8.85 37.15 0.33 -0.79 16 12 C -0.18 -2.08 8.37 37.16 0.31 -1.77 16 13 C 0.49 11.81 6.01 -10.99 -0.07 11.75 16 14 O -0.53 -15.25 14.17 5.55 0.08 -15.17 16 15 N -0.45 -10.37 2.52 -105.78 -0.27 -10.64 16 16 C -0.34 -8.76 7.24 -14.93 -0.11 -8.87 16 17 H 0.07 1.90 7.14 -52.49 -0.37 1.52 16 18 C 0.07 1.92 4.54 -17.47 -0.08 1.84 16 19 N -0.84 -23.86 11.77 55.50 0.65 -23.21 16 20 Si 0.91 21.60 29.57 -169.99 -5.03 16.57 16 21 H -0.27 -6.67 7.11 56.52 0.40 -6.27 16 22 H -0.29 -6.86 7.11 56.52 0.40 -6.46 16 23 H -0.28 -6.32 7.11 56.52 0.40 -5.92 16 24 C 0.13 2.33 2.54 -66.08 -0.17 2.16 16 25 C -0.12 -2.09 6.38 -24.55 -0.16 -2.25 16 26 C -0.13 -1.74 6.82 -25.07 -0.17 -1.91 16 27 C -0.12 -1.63 6.80 -25.61 -0.17 -1.80 16 28 C -0.14 -2.22 6.44 -25.62 -0.16 -2.38 16 29 H 0.05 0.61 7.28 -51.92 -0.38 0.23 16 30 H 0.06 1.11 6.99 -51.93 -0.36 0.75 16 31 H 0.06 1.20 8.14 -51.93 -0.42 0.78 16 32 H 0.08 1.65 8.07 -51.92 -0.42 1.23 16 33 H 0.07 1.42 8.14 -51.93 -0.42 1.00 16 34 H 0.40 4.78 7.78 -40.82 -0.32 4.47 16 35 H 0.06 0.62 8.14 -51.93 -0.42 0.20 16 36 H 0.07 0.61 8.14 -51.93 -0.42 0.19 16 37 H 0.09 1.27 5.88 -51.93 -0.31 0.96 16 38 H 0.07 0.51 8.14 -51.93 -0.42 0.09 16 39 H 0.08 0.48 8.14 -51.93 -0.42 0.06 16 40 H 0.07 0.51 8.14 -51.93 -0.42 0.08 16 41 H 0.09 1.29 5.88 -51.93 -0.31 0.98 16 42 H 0.07 0.62 8.14 -51.93 -0.42 0.19 16 43 H 0.06 0.64 8.14 -51.93 -0.42 0.22 16 44 H 0.26 7.34 8.31 -40.82 -0.34 7.00 16 45 H 0.07 1.18 4.88 -51.93 -0.25 0.93 16 46 H 0.06 0.99 8.14 -51.93 -0.42 0.57 16 47 H 0.09 1.81 6.98 -51.93 -0.36 1.44 16 48 H 0.06 0.70 8.14 -51.93 -0.42 0.28 16 49 H 0.06 0.75 8.14 -51.93 -0.42 0.32 16 50 H 0.06 0.80 8.14 -51.93 -0.42 0.38 16 51 H 0.05 0.63 8.14 -51.93 -0.42 0.21 16 52 H 0.07 1.39 6.91 -51.93 -0.36 1.03 16 53 H 0.05 0.77 8.14 -51.93 -0.42 0.35 16 LS Contribution 407.99 15.07 6.15 6.15 Total: -1.00 -32.08 407.99 -8.33 -40.41 By element: Atomic # 1 Polarization: 16.98 SS G_CDS: -9.26 Total: 7.72 kcal Atomic # 6 Polarization: 4.61 SS G_CDS: 0.33 Total: 4.94 kcal Atomic # 7 Polarization: -44.58 SS G_CDS: 0.04 Total: -44.54 kcal Atomic # 8 Polarization: -30.70 SS G_CDS: -0.55 Total: -31.25 kcal Atomic # 14 Polarization: 21.60 SS G_CDS: -5.03 Total: 16.57 kcal Total LS contribution 6.15 Total: 6.15 kcal Total: -32.08 -8.33 -40.41 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. ZINC000495617541.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 -85.334 kcal (2) G-P(sol) polarization free energy of solvation -32.085 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -117.419 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.325 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.410 kcal (6) G-S(sol) free energy of system = (1) + (5) -125.744 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.77 seconds