Wall clock time and date at job start Sat Apr 11 2020 06:10:25 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.52991 * 1 3 3 C 1.52999 * 109.47744 * 2 1 4 4 C 1.53001 * 109.47165 * 239.99826 * 2 1 3 5 5 O 1.45209 * 109.47276 * 119.99943 * 2 1 3 6 6 C 1.34636 * 116.99514 * 59.99799 * 5 2 1 7 7 O 1.21509 * 120.00675 * 359.97438 * 6 5 2 8 8 N 1.34778 * 119.99556 * 180.02562 * 6 5 2 9 9 C 1.46500 * 119.99850 * 179.97438 * 8 6 5 10 10 C 1.53004 * 109.47150 * 179.97438 * 9 8 6 11 11 C 1.53346 * 109.41904 * 295.13157 * 10 9 8 12 12 C 1.53195 * 109.14901 * 178.67276 * 11 10 9 13 13 N 1.46919 * 108.89273 * 54.57640 * 12 11 10 14 14 C 1.34778 * 120.63255 * 126.34238 * 13 12 11 15 15 O 1.21279 * 120.00169 * 359.97438 * 14 13 12 16 16 C 1.50698 * 120.00023 * 179.97438 * 14 13 12 17 17 C 1.53003 * 109.46931 * 180.02562 * 16 14 13 18 18 C 1.50706 * 109.47067 * 179.97438 * 17 16 14 19 19 H 1.07997 * 119.99699 * 0.02562 * 18 17 16 20 20 C 1.37874 * 119.99846 * 179.97438 * 18 17 16 21 21 N 1.31517 * 120.00051 * 0.02562 * 20 18 17 22 22 Si 1.86808 * 119.99888 * 180.02562 * 20 18 17 23 23 H 1.48498 * 109.46775 * 60.00383 * 22 20 18 24 24 H 1.48500 * 109.47178 * 180.02562 * 22 20 18 25 25 H 1.48502 * 109.46977 * 300.00371 * 22 20 18 26 26 C 1.46933 * 118.74003 * 306.32349 * 13 12 11 27 27 C 1.53034 * 109.46014 * 175.00023 * 10 9 8 28 28 H 1.09007 * 109.46853 * 300.00619 * 1 2 3 29 29 H 1.09001 * 109.47339 * 59.99878 * 1 2 3 30 30 H 1.08996 * 109.47849 * 180.02562 * 1 2 3 31 31 H 1.08999 * 109.47499 * 180.02562 * 3 2 1 32 32 H 1.09003 * 109.46461 * 300.00045 * 3 2 1 33 33 H 1.08997 * 109.47332 * 59.99745 * 3 2 1 34 34 H 1.09008 * 109.46851 * 300.00322 * 4 2 1 35 35 H 1.09003 * 109.47253 * 59.99825 * 4 2 1 36 36 H 1.08999 * 109.47279 * 180.02562 * 4 2 1 37 37 H 0.96999 * 120.00446 * 359.97438 * 8 6 5 38 38 H 1.09001 * 109.47353 * 60.00073 * 9 8 6 39 39 H 1.08997 * 109.47478 * 299.99831 * 9 8 6 40 40 H 1.09001 * 109.45799 * 55.03185 * 10 9 8 41 41 H 1.08997 * 109.68637 * 58.70485 * 11 10 9 42 42 H 1.09002 * 109.47277 * 298.50741 * 11 10 9 43 43 H 1.09002 * 109.58720 * 294.71670 * 12 11 10 44 44 H 1.08996 * 109.58651 * 174.43567 * 12 11 10 45 45 H 1.09001 * 109.47488 * 299.99559 * 16 14 13 46 46 H 1.09001 * 109.47269 * 59.99982 * 16 14 13 47 47 H 1.08996 * 109.47074 * 299.99864 * 17 16 14 48 48 H 1.09001 * 109.47188 * 59.99897 * 17 16 14 49 49 H 0.97008 * 119.99666 * 180.02562 * 21 20 18 50 50 H 1.08994 * 109.58580 * 173.39910 * 26 13 12 51 51 H 1.09006 * 109.57831 * 293.67403 * 26 13 12 52 52 H 1.08994 * 109.50232 * 61.43679 * 27 10 9 53 53 H 1.09007 * 109.52628 * 301.35340 * 27 10 9 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.5299 0.0000 0.0000 3 6 2.0401 1.4424 0.0000 4 6 2.0399 -0.7213 -1.2492 5 8 2.0140 -0.6845 1.1856 6 6 1.6521 -1.9723 1.3383 7 8 0.9507 -2.5095 0.5041 8 7 2.0671 -2.6654 2.4172 9 6 1.6730 -4.0665 2.5837 10 6 2.2732 -4.6146 3.8800 11 6 3.8024 -4.6304 3.7670 12 6 4.3980 -5.1488 5.0798 13 7 3.8007 -6.4547 5.3903 14 6 4.5810 -7.5269 5.6312 15 8 5.7884 -7.4193 5.5923 16 6 3.9532 -8.8595 5.9488 17 6 5.0520 -9.8990 6.1791 18 6 4.4242 -11.2316 6.4974 19 1 3.3490 -11.3275 6.5315 20 6 5.2224 -12.3284 6.7439 21 7 6.5317 -12.2116 6.7018 22 14 4.4442 -13.9803 7.1377 23 1 3.5884 -14.4076 6.0018 24 1 5.5107 -14.9892 7.3613 25 1 3.6150 -13.8572 8.3635 26 6 2.3365 -6.5695 5.4336 27 6 1.7610 -6.0370 4.1176 28 1 -0.3633 0.5140 0.8900 29 1 -0.3634 0.5138 -0.8900 30 1 -0.3635 -1.0276 0.0005 31 1 3.1301 1.4424 -0.0005 32 1 1.6767 1.9562 -0.8900 33 1 1.6768 1.9563 0.8899 34 1 1.6765 -1.7490 -1.2492 35 1 1.6766 -0.2075 -2.1393 36 1 3.1299 -0.7217 -1.2490 37 1 2.6268 -2.2366 3.0834 38 1 2.0394 -4.6499 1.7390 39 1 0.5861 -4.1354 2.6286 40 1 1.9780 -3.9783 4.7143 41 1 4.1669 -3.6212 3.5756 42 1 4.0996 -5.2863 2.9487 43 1 4.1781 -4.4460 5.8835 44 1 5.4774 -5.2566 4.9739 45 1 3.3446 -8.7693 6.8486 46 1 3.3248 -9.1731 5.1151 47 1 5.6609 -9.9894 5.2796 48 1 5.6800 -9.5852 7.0129 49 1 7.0933 -12.9834 6.8749 50 1 2.0556 -7.6152 5.5587 51 1 1.9488 -5.9846 6.2677 52 1 2.0772 -6.6795 3.2959 53 1 0.6725 -6.0276 4.1740 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 ZINC000497196055.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:25 Heat of formation + Delta-G solvation = -141.920197 kcal Electronic energy + Delta-G solvation = -29813.967721 eV Core-core repulsion = 25684.851341 eV Total energy + Delta-G solvation = -4129.116380 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.76 seconds Orbital eigenvalues (eV) -41.09993 -39.98942 -38.27437 -37.19302 -35.25013 -34.37258 -33.52482 -32.59926 -31.62413 -30.53891 -27.81735 -27.76000 -27.75325 -25.76525 -23.82757 -23.44721 -22.19965 -21.55692 -20.60413 -20.17052 -18.68783 -18.34355 -17.02587 -16.96556 -16.57274 -15.92782 -15.66381 -15.44803 -15.23239 -15.05784 -14.91648 -14.58393 -14.41348 -13.97470 -13.86081 -13.64137 -13.46547 -13.07532 -12.90418 -12.86135 -12.56320 -12.37445 -12.30288 -12.29116 -12.11823 -12.05127 -11.78141 -11.69753 -11.40455 -11.26034 -11.14564 -10.90962 -10.86684 -10.74689 -10.67960 -10.58651 -10.36825 -9.95641 -9.61046 -9.50957 -9.18139 -8.99232 -8.45494 -8.30750 -5.95569 -4.00162 1.83821 2.54046 3.19963 3.26816 3.34588 3.39529 3.41207 3.78901 4.31134 4.31241 4.35720 4.40095 4.47458 4.49949 4.51212 4.91217 4.99274 5.06323 5.11886 5.13926 5.16282 5.19045 5.28655 5.29772 5.35539 5.39252 5.40574 5.44957 5.46036 5.54171 5.58393 5.66910 5.67374 5.69386 5.82419 5.87674 5.90009 5.97530 6.03577 6.13996 6.20208 6.20358 6.21494 6.72584 6.90472 7.07741 7.73367 7.74659 7.86393 8.06587 8.07458 8.38578 8.73853 9.08045 9.63330 11.15272 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.028617 B = 0.001649 C = 0.001621 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 978.209614 B =16978.453674 C =17267.088774 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.183 4.183 2 C 0.136 3.864 3 C -0.141 4.141 4 C -0.184 4.184 5 O -0.366 6.366 6 C 0.649 3.351 7 O -0.568 6.568 8 N -0.709 5.709 9 C 0.142 3.858 10 C -0.098 4.098 11 C -0.130 4.130 12 C 0.110 3.890 13 N -0.626 5.626 14 C 0.513 3.487 15 O -0.534 6.534 16 C -0.129 4.129 17 C 0.037 3.963 18 C -0.528 4.528 19 H 0.055 0.945 20 C 0.058 3.942 21 N -0.889 5.889 22 Si 0.893 3.107 23 H -0.277 1.277 24 H -0.287 1.287 25 H -0.277 1.277 26 C 0.110 3.890 27 C -0.127 4.127 28 H 0.059 0.941 29 H 0.070 0.930 30 H 0.088 0.912 31 H 0.067 0.933 32 H 0.080 0.920 33 H 0.066 0.934 34 H 0.088 0.912 35 H 0.070 0.930 36 H 0.060 0.940 37 H 0.406 0.594 38 H 0.073 0.927 39 H 0.072 0.928 40 H 0.081 0.919 41 H 0.071 0.929 42 H 0.071 0.929 43 H 0.060 0.940 44 H 0.092 0.908 45 H 0.079 0.921 46 H 0.079 0.921 47 H 0.023 0.977 48 H 0.023 0.977 49 H 0.261 0.739 50 H 0.095 0.905 51 H 0.065 0.935 52 H 0.074 0.926 53 H 0.079 0.921 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -13.706 28.027 -11.354 33.201 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.241 4.241 2 C 0.100 3.900 3 C -0.198 4.198 4 C -0.241 4.241 5 O -0.281 6.281 6 C 0.402 3.598 7 O -0.458 6.458 8 N -0.369 5.369 9 C 0.020 3.980 10 C -0.117 4.117 11 C -0.168 4.168 12 C -0.013 4.013 13 N -0.362 5.362 14 C 0.303 3.697 15 O -0.409 6.409 16 C -0.168 4.168 17 C -0.001 4.001 18 C -0.566 4.566 19 H 0.073 0.927 20 C -0.143 4.143 21 N -0.528 5.528 22 Si 0.721 3.279 23 H -0.203 1.203 24 H -0.213 1.213 25 H -0.203 1.203 26 C -0.013 4.013 27 C -0.166 4.166 28 H 0.078 0.922 29 H 0.089 0.911 30 H 0.106 0.894 31 H 0.086 0.914 32 H 0.099 0.901 33 H 0.085 0.915 34 H 0.106 0.894 35 H 0.089 0.911 36 H 0.079 0.921 37 H 0.241 0.759 38 H 0.091 0.909 39 H 0.091 0.909 40 H 0.099 0.901 41 H 0.089 0.911 42 H 0.090 0.910 43 H 0.079 0.921 44 H 0.110 0.890 45 H 0.097 0.903 46 H 0.097 0.903 47 H 0.041 0.959 48 H 0.041 0.959 49 H 0.071 0.929 50 H 0.113 0.887 51 H 0.084 0.916 52 H 0.093 0.907 53 H 0.097 0.903 Dipole moment (debyes) X Y Z Total from point charges -13.216 27.756 -11.816 32.934 hybrid contribution -0.243 -0.636 0.423 0.802 sum -13.458 27.120 -11.393 32.348 Atomic orbital electron populations 1.22516 0.93974 1.05087 1.02487 1.22072 0.94362 0.91510 0.82076 1.21867 1.02330 0.92499 1.03136 1.22523 1.01638 1.03533 0.96438 1.86302 1.71598 1.23676 1.46529 1.18098 0.78931 0.80509 0.82284 1.90910 1.41105 1.70461 1.43276 1.44781 1.55638 1.12829 1.23635 1.21456 0.99926 0.78490 0.98176 1.21470 0.95931 0.97584 0.96747 1.22063 0.94286 1.01562 0.98904 1.21445 0.99874 0.84510 0.95423 1.48300 1.06410 1.09489 1.72024 1.21592 0.87219 0.85973 0.74884 1.90657 1.14745 1.84911 1.50618 1.21881 1.01939 0.91400 1.01613 1.18907 0.92178 0.90805 0.98195 1.21178 0.92840 0.94172 1.48419 0.92669 1.25022 0.94802 0.99016 0.95497 1.69992 1.03240 1.30485 1.49082 0.86655 0.78829 0.83625 0.78823 1.20257 1.21292 1.20259 1.21665 0.80022 1.03015 0.96624 1.21898 1.00735 0.94995 0.98980 0.92197 0.91116 0.89362 0.91412 0.90110 0.91464 0.89351 0.91112 0.92114 0.75898 0.90906 0.90931 0.90061 0.91062 0.91011 0.92133 0.88994 0.90294 0.90264 0.95862 0.95854 0.92905 0.88691 0.91640 0.90711 0.90260 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.18 -1.37 8.37 37.15 0.31 -1.05 16 2 C 0.14 1.05 1.13 -90.62 -0.10 0.95 16 3 C -0.14 -0.74 8.85 37.16 0.33 -0.41 16 4 C -0.18 -1.43 8.37 37.16 0.31 -1.12 16 5 O -0.37 -3.41 9.94 -40.31 -0.40 -3.81 16 6 C 0.65 6.62 8.06 54.06 0.44 7.05 16 7 O -0.57 -7.10 11.88 -19.28 -0.23 -7.33 16 8 N -0.71 -5.55 5.28 -65.28 -0.34 -5.90 16 9 C 0.14 1.08 5.07 -4.04 -0.02 1.06 16 10 C -0.10 -0.70 2.19 -90.42 -0.20 -0.89 16 11 C -0.13 -1.13 5.27 -26.45 -0.14 -1.27 16 12 C 0.11 1.28 6.42 -3.71 -0.02 1.26 16 13 N -0.63 -8.68 2.97 -172.77 -0.51 -9.20 16 14 C 0.51 9.85 7.70 -10.99 -0.08 9.76 16 15 O -0.53 -12.58 15.39 5.56 0.09 -12.49 16 16 C -0.13 -2.58 3.99 -27.88 -0.11 -2.70 16 17 C 0.04 0.99 4.26 -27.88 -0.12 0.87 16 18 C -0.53 -14.91 9.11 -34.44 -0.31 -15.22 16 19 H 0.06 1.46 7.74 -52.49 -0.41 1.05 16 20 C 0.06 1.69 5.46 -17.82 -0.10 1.60 16 21 N -0.89 -27.41 13.29 55.43 0.74 -26.67 16 22 Si 0.89 21.60 29.54 -169.99 -5.02 16.58 16 23 H -0.28 -6.54 7.11 56.52 0.40 -6.14 16 24 H -0.29 -6.91 7.11 56.52 0.40 -6.51 16 25 H -0.28 -6.53 7.11 56.52 0.40 -6.13 16 26 C 0.11 1.11 6.34 -3.71 -0.02 1.09 16 27 C -0.13 -1.01 5.36 -26.62 -0.14 -1.15 16 28 H 0.06 0.42 8.14 -51.92 -0.42 -0.01 16 29 H 0.07 0.40 8.14 -51.93 -0.42 -0.02 16 30 H 0.09 0.86 5.88 -51.93 -0.31 0.56 16 31 H 0.07 0.35 8.14 -51.93 -0.42 -0.07 16 32 H 0.08 0.29 8.14 -51.93 -0.42 -0.13 16 33 H 0.07 0.33 8.14 -51.93 -0.42 -0.09 16 34 H 0.09 0.90 5.88 -51.92 -0.31 0.59 16 35 H 0.07 0.42 8.14 -51.93 -0.42 0.00 16 36 H 0.06 0.46 8.14 -51.93 -0.42 0.04 16 37 H 0.41 2.43 7.76 -40.82 -0.32 2.11 16 38 H 0.07 0.65 8.14 -51.93 -0.42 0.23 16 39 H 0.07 0.56 8.14 -51.93 -0.42 0.14 16 40 H 0.08 0.50 8.14 -51.93 -0.42 0.08 16 41 H 0.07 0.50 7.16 -51.93 -0.37 0.13 16 42 H 0.07 0.75 8.14 -51.93 -0.42 0.33 16 43 H 0.06 0.60 8.14 -51.93 -0.42 0.18 16 44 H 0.09 1.30 7.00 -51.93 -0.36 0.93 16 45 H 0.08 1.45 7.94 -51.93 -0.41 1.03 16 46 H 0.08 1.45 7.69 -51.93 -0.40 1.05 16 47 H 0.02 0.68 8.10 -51.93 -0.42 0.26 16 48 H 0.02 0.68 8.10 -51.93 -0.42 0.26 16 49 H 0.26 7.58 8.31 -40.82 -0.34 7.24 16 50 H 0.09 1.00 5.93 -51.93 -0.31 0.69 16 51 H 0.07 0.56 8.14 -51.93 -0.42 0.14 16 52 H 0.07 0.68 8.14 -51.93 -0.42 0.26 16 53 H 0.08 0.49 8.14 -51.92 -0.42 0.06 16 LS Contribution 415.20 15.07 6.26 6.26 Total: -1.00 -35.57 415.20 -8.93 -44.49 By element: Atomic # 1 Polarization: 7.76 SS G_CDS: -9.51 Total: -1.74 kcal Atomic # 6 Polarization: -0.19 SS G_CDS: 0.01 Total: -0.18 kcal Atomic # 7 Polarization: -41.64 SS G_CDS: -0.12 Total: -41.76 kcal Atomic # 8 Polarization: -23.09 SS G_CDS: -0.54 Total: -23.63 kcal Atomic # 14 Polarization: 21.60 SS G_CDS: -5.02 Total: 16.58 kcal Total LS contribution 6.26 Total: 6.26 kcal Total: -35.57 -8.93 -44.49 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. ZINC000497196055.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 -97.428 kcal (2) G-P(sol) polarization free energy of solvation -35.565 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -132.994 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.926 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.492 kcal (6) G-S(sol) free energy of system = (1) + (5) -141.920 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.76 seconds