Wall clock time and date at job start Tue Mar 31 2020 20:30:21 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.53001 * 1 3 3 C 1.50702 * 109.46975 * 2 1 4 4 C 1.32988 * 119.99877 * 125.26569 * 3 2 1 5 5 C 1.50707 * 119.99997 * 175.69777 * 4 3 2 6 6 C 1.47100 * 120.00242 * 355.70484 * 4 3 2 7 7 O 1.21621 * 119.99568 * 263.32139 * 6 4 3 8 8 N 1.34770 * 120.00839 * 83.32020 * 6 4 3 9 9 C 1.46929 * 120.57706 * 174.84018 * 8 6 4 10 10 C 1.53302 * 108.71061 * 125.86649 * 9 8 6 11 11 C 1.52424 * 109.38595 * 53.70864 * 10 9 8 12 12 C 1.53971 * 110.25032 * 62.55551 * 11 10 9 13 13 O 1.44552 * 104.69163 * 94.88878 * 12 11 10 14 14 C 1.44467 * 105.35992 * 40.02080 * 13 12 11 15 15 H 1.09008 * 110.36988 * 78.77912 * 14 13 12 16 16 C 1.52995 * 110.36727 * 200.98910 * 14 13 12 17 17 N 1.46507 * 109.47327 * 298.54838 * 16 14 13 18 18 C 1.36935 * 119.99783 * 179.97438 * 17 16 14 19 19 H 1.08000 * 120.00717 * 0.02562 * 18 17 16 20 20 C 1.38813 * 119.99576 * 180.02562 * 18 17 16 21 21 N 1.31619 * 120.00391 * 359.72548 * 20 18 17 22 22 Si 1.86803 * 119.99989 * 180.02562 * 20 18 17 23 23 H 1.48498 * 109.47273 * 60.00004 * 22 20 18 24 24 H 1.48497 * 109.47258 * 180.02562 * 22 20 18 25 25 H 1.48507 * 109.47184 * 299.99984 * 22 20 18 26 26 C 1.53985 * 110.26886 * 177.53776 * 11 10 9 27 27 C 1.52406 * 111.00537 * 300.07326 * 11 10 9 28 28 C 1.46935 * 120.56380 * 354.86608 * 8 6 4 29 29 H 1.08998 * 109.47028 * 299.99761 * 1 2 3 30 30 H 1.09003 * 109.47166 * 60.00202 * 1 2 3 31 31 H 1.09004 * 109.46806 * 179.97438 * 1 2 3 32 32 H 1.08999 * 109.46895 * 239.99701 * 2 1 3 33 33 H 1.09003 * 109.47166 * 119.99798 * 2 1 3 34 34 H 1.07997 * 119.99699 * 305.26493 * 3 2 1 35 35 H 1.09001 * 109.46973 * 275.04584 * 5 4 3 36 36 H 1.08995 * 109.46557 * 35.04089 * 5 4 3 37 37 H 1.08993 * 109.47015 * 155.04394 * 5 4 3 38 38 H 1.08997 * 109.59550 * 6.10757 * 9 8 6 39 39 H 1.08996 * 109.72974 * 245.71434 * 9 8 6 40 40 H 1.08993 * 109.48986 * 293.73234 * 10 9 8 41 41 H 1.09000 * 109.48283 * 173.67784 * 10 9 8 42 42 H 1.09000 * 110.39499 * 336.08793 * 12 11 10 43 43 H 1.08997 * 110.48963 * 213.74099 * 12 11 10 44 44 H 1.08998 * 109.47372 * 178.54804 * 16 14 13 45 45 H 1.09000 * 109.47289 * 58.54832 * 16 14 13 46 46 H 0.96999 * 119.99600 * 0.02562 * 17 16 14 47 47 H 0.96997 * 120.00269 * 179.97438 * 21 20 18 48 48 H 1.09000 * 110.54477 * 122.76650 * 26 11 10 49 49 H 1.08997 * 110.54695 * 359.96343 * 26 11 10 50 50 H 1.08993 * 109.48575 * 179.88978 * 27 11 10 51 51 H 1.09003 * 109.48053 * 299.93896 * 27 11 10 52 52 H 1.08997 * 109.59820 * 353.89070 * 28 8 6 53 53 H 1.09001 * 109.59000 * 114.19659 * 28 8 6 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.0323 1.4208 0.0000 4 6 2.8809 1.8261 0.9404 5 6 3.4585 3.2170 0.8854 6 6 3.2558 0.9168 2.0341 7 8 2.6813 0.9831 3.1040 8 7 4.2359 0.0101 1.8508 9 6 4.6979 -0.8270 2.9665 10 6 4.6091 -2.2988 2.5466 11 6 5.3633 -2.4978 1.2371 12 6 6.8561 -2.1722 1.4272 13 8 7.4776 -3.4482 1.7012 14 6 6.8311 -4.3872 0.8138 15 1 7.2399 -4.3010 -0.1930 16 6 7.0000 -5.8156 1.3354 17 7 8.4250 -6.1524 1.3843 18 6 8.8219 -7.3876 1.8224 19 1 8.0845 -8.1136 2.1315 20 6 10.1721 -7.7064 1.8693 21 7 11.0705 -6.8236 1.4875 22 14 10.7136 -9.3916 2.4663 23 1 10.1279 -10.4386 1.5911 24 1 12.1952 -9.4802 2.4213 25 1 10.2482 -9.5956 3.8617 26 6 5.3392 -3.9821 0.8279 27 6 4.7717 -1.6344 0.1293 28 6 4.8613 -0.1579 0.5318 29 1 -0.3633 0.5138 0.8900 30 1 -0.3634 0.5138 -0.8900 31 1 -0.3633 -1.0277 -0.0005 32 1 1.8933 -0.5139 -0.8900 33 1 1.8934 -0.5138 0.8900 34 1 1.7031 2.1099 -0.7637 35 1 4.3357 3.2207 0.2384 36 1 2.7118 3.9050 0.4890 37 1 3.7458 3.5310 1.8888 38 1 4.0655 -0.6567 3.8377 39 1 5.7307 -0.5781 3.2101 40 1 3.5637 -2.5754 2.4098 41 1 5.0528 -2.9250 3.3206 42 1 6.9928 -1.4945 2.2698 43 1 7.2701 -1.7355 0.5184 44 1 6.4855 -6.5086 0.6697 45 1 6.5746 -5.8897 2.3362 46 1 9.0872 -5.5001 1.1070 47 1 12.0139 -7.0463 1.5206 48 1 4.9019 -4.0984 -0.1637 49 1 4.7895 -4.5722 1.5611 50 1 5.3296 -1.7941 -0.7933 51 1 3.7275 -1.9063 -0.0252 52 1 4.3365 0.4536 -0.2021 53 1 5.9066 0.1471 0.5815 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 ZINC001276912641.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:30:21 Heat of formation + Delta-G solvation = -63.800298 kcal Electronic energy + Delta-G solvation = -30461.924017 eV Core-core repulsion = 26526.651221 eV Total energy + Delta-G solvation = -3935.272796 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.89 seconds Orbital eigenvalues (eV) -41.12181 -39.67152 -36.93919 -36.21443 -34.45728 -33.39631 -32.60160 -31.61826 -31.45047 -29.58065 -27.49255 -27.17297 -25.30149 -24.57300 -23.75890 -22.99413 -22.45130 -21.62078 -20.64060 -19.41477 -18.21841 -17.08486 -16.94051 -16.27262 -15.96150 -15.83777 -15.41944 -15.26711 -14.86155 -14.64804 -14.37932 -14.17941 -14.00783 -13.62643 -13.45505 -13.27864 -13.19380 -13.08065 -12.96855 -12.79257 -12.75453 -12.56086 -12.36491 -12.03590 -11.79614 -11.57592 -11.48458 -11.38072 -11.34822 -11.15567 -11.01817 -10.76889 -10.59304 -10.45256 -10.35368 -10.06331 -9.98983 -9.61795 -9.34196 -9.32264 -8.87845 -8.82718 -6.84974 -5.69306 -2.97813 1.38503 2.31293 2.74803 3.47899 3.51748 3.53185 4.06848 4.11108 4.32109 4.38752 4.50301 4.59906 4.64503 4.65950 4.77528 4.89907 4.96592 4.99126 4.99848 5.01574 5.06075 5.17417 5.23704 5.26806 5.30138 5.35265 5.42241 5.43730 5.47667 5.53417 5.55384 5.62227 5.65639 5.68479 5.77934 5.80162 5.81799 5.90234 5.94295 5.98319 6.03529 6.07793 6.16444 6.17449 6.29058 6.38056 6.59484 6.78323 7.09053 7.27164 7.73801 7.86625 8.49144 9.60421 10.19234 10.55040 11.53441 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.025937 B = 0.002212 C = 0.002162 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1079.260269 B =12655.923084 C =12945.703499 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.107 4.107 3 C -0.114 4.114 4 C -0.164 4.164 5 C -0.104 4.104 6 C 0.561 3.439 7 O -0.539 6.539 8 N -0.604 5.604 9 C 0.106 3.894 10 C -0.116 4.116 11 C -0.090 4.090 12 C 0.049 3.951 13 O -0.351 6.351 14 C 0.071 3.929 15 H 0.066 0.934 16 C 0.160 3.840 17 N -0.728 5.728 18 C -0.239 4.239 19 H 0.069 0.931 20 C -0.014 4.014 21 N -0.932 5.932 22 Si 0.907 3.093 23 H -0.286 1.286 24 H -0.291 1.291 25 H -0.286 1.286 26 C -0.138 4.138 27 C -0.120 4.120 28 C 0.105 3.895 29 H 0.056 0.944 30 H 0.056 0.944 31 H 0.058 0.942 32 H 0.079 0.921 33 H 0.069 0.931 34 H 0.126 0.874 35 H 0.073 0.927 36 H 0.072 0.928 37 H 0.066 0.934 38 H 0.094 0.906 39 H 0.073 0.927 40 H 0.067 0.933 41 H 0.085 0.915 42 H 0.087 0.913 43 H 0.055 0.945 44 H 0.029 0.971 45 H 0.024 0.976 46 H 0.388 0.612 47 H 0.255 0.745 48 H 0.075 0.925 49 H 0.079 0.921 50 H 0.080 0.920 51 H 0.071 0.929 52 H 0.088 0.912 53 H 0.072 0.928 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -19.763 20.453 -6.482 29.170 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.204 4.204 2 C -0.144 4.144 3 C -0.133 4.133 4 C -0.167 4.167 5 C -0.161 4.161 6 C 0.350 3.650 7 O -0.417 6.417 8 N -0.338 5.338 9 C -0.016 4.016 10 C -0.154 4.154 11 C -0.090 4.090 12 C -0.028 4.028 13 O -0.269 6.269 14 C 0.013 3.987 15 H 0.084 0.916 16 C 0.032 3.968 17 N -0.372 5.372 18 C -0.371 4.371 19 H 0.087 0.913 20 C -0.210 4.210 21 N -0.581 5.581 22 Si 0.738 3.262 23 H -0.213 1.213 24 H -0.217 1.217 25 H -0.213 1.213 26 C -0.176 4.176 27 C -0.158 4.158 28 C -0.018 4.018 29 H 0.075 0.925 30 H 0.075 0.925 31 H 0.077 0.923 32 H 0.098 0.902 33 H 0.088 0.912 34 H 0.143 0.857 35 H 0.092 0.908 36 H 0.091 0.909 37 H 0.085 0.915 38 H 0.112 0.888 39 H 0.091 0.909 40 H 0.086 0.914 41 H 0.103 0.897 42 H 0.105 0.895 43 H 0.074 0.926 44 H 0.047 0.953 45 H 0.042 0.958 46 H 0.222 0.778 47 H 0.065 0.935 48 H 0.094 0.906 49 H 0.097 0.903 50 H 0.099 0.901 51 H 0.089 0.911 52 H 0.106 0.894 53 H 0.090 0.910 Dipole moment (debyes) X Y Z Total from point charges -19.878 20.263 -5.677 28.947 hybrid contribution 0.469 -0.678 -0.003 0.824 sum -19.409 19.585 -5.680 28.152 Atomic orbital electron populations 1.21684 0.94528 1.02074 1.02083 1.20837 0.96913 0.93791 1.02884 1.23368 0.95734 0.97555 0.96642 1.21510 1.01786 0.97410 0.95976 1.20669 1.00937 0.91516 1.03020 1.18531 0.78780 0.79870 0.87827 1.90815 1.52680 1.70931 1.27236 1.48171 1.36909 1.37485 1.11195 1.21764 1.00373 0.89695 0.89732 1.21848 1.00686 0.95770 0.97129 1.21483 0.96483 0.95610 0.95452 1.23500 0.90049 0.88397 1.00882 1.89357 1.60417 1.16728 1.60414 1.23291 0.90041 0.91530 0.93856 0.91605 1.20735 0.88664 0.89660 0.97695 1.42259 1.00411 1.15418 1.79102 1.19183 0.91185 0.90967 1.35739 0.91329 1.24744 0.95827 0.98087 1.02377 1.69666 1.00396 1.34259 1.53736 0.86267 0.77362 0.83948 0.78594 1.21322 1.21707 1.21299 1.23009 0.97981 0.94332 1.02258 1.21921 1.01964 0.93754 0.98167 1.21829 0.98707 0.97224 0.84049 0.92467 0.92457 0.92287 0.90222 0.91249 0.85667 0.90792 0.90947 0.91511 0.88803 0.90878 0.91381 0.89663 0.89534 0.92640 0.95338 0.95802 0.77838 0.93528 0.90631 0.90266 0.90117 0.91060 0.89397 0.90953 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 -0.96 9.89 37.16 0.37 -0.59 16 2 C -0.11 -0.75 6.01 -27.88 -0.17 -0.92 16 3 C -0.11 -0.77 8.42 -36.13 -0.30 -1.07 16 4 C -0.16 -1.35 4.14 -102.85 -0.43 -1.78 16 5 C -0.10 -0.66 10.02 36.01 0.36 -0.29 16 6 C 0.56 6.22 5.14 -12.67 -0.07 6.16 16 7 O -0.54 -8.13 16.52 5.26 0.09 -8.05 16 8 N -0.60 -5.73 2.84 -174.02 -0.49 -6.23 16 9 C 0.11 1.08 6.13 -3.65 -0.02 1.06 16 10 C -0.12 -1.20 4.69 -26.87 -0.13 -1.32 16 11 C -0.09 -0.93 0.68 -154.09 -0.11 -1.03 16 12 C 0.05 0.61 5.63 37.66 0.21 0.82 16 13 O -0.35 -6.32 11.20 -35.23 -0.39 -6.71 16 14 C 0.07 1.17 3.37 -25.89 -0.09 1.09 16 15 H 0.07 1.10 8.14 -51.92 -0.42 0.67 16 16 C 0.16 3.32 5.89 -4.04 -0.02 3.30 16 17 N -0.73 -19.33 5.53 -59.91 -0.33 -19.66 16 18 C -0.24 -6.93 9.99 -15.06 -0.15 -7.08 16 19 H 0.07 1.91 7.83 -52.49 -0.41 1.50 16 20 C -0.01 -0.43 5.50 -17.43 -0.10 -0.52 16 21 N -0.93 -29.05 13.06 55.49 0.72 -28.33 16 22 Si 0.91 22.89 29.55 -169.99 -5.02 17.86 16 23 H -0.29 -7.06 7.11 56.52 0.40 -6.65 16 24 H -0.29 -7.37 7.11 56.52 0.40 -6.97 16 25 H -0.29 -7.11 7.11 56.52 0.40 -6.71 16 26 C -0.14 -1.66 5.90 -25.37 -0.15 -1.81 16 27 C -0.12 -0.94 4.55 -26.88 -0.12 -1.06 16 28 C 0.10 0.77 5.89 -3.65 -0.02 0.75 16 29 H 0.06 0.44 8.14 -51.93 -0.42 0.01 16 30 H 0.06 0.30 8.14 -51.93 -0.42 -0.12 16 31 H 0.06 0.37 8.14 -51.93 -0.42 -0.05 16 32 H 0.08 0.46 8.14 -51.93 -0.42 0.04 16 33 H 0.07 0.64 6.20 -51.93 -0.32 0.31 16 34 H 0.13 0.59 8.06 -52.49 -0.42 0.17 16 35 H 0.07 0.36 8.14 -51.93 -0.42 -0.06 16 36 H 0.07 0.36 8.10 -51.93 -0.42 -0.06 16 37 H 0.07 0.49 8.08 -51.93 -0.42 0.07 16 38 H 0.09 1.07 6.99 -51.93 -0.36 0.71 16 39 H 0.07 0.71 6.44 -51.93 -0.33 0.38 16 40 H 0.07 0.70 8.14 -51.93 -0.42 0.27 16 41 H 0.08 0.98 8.14 -51.93 -0.42 0.55 16 42 H 0.09 1.01 5.99 -51.93 -0.31 0.70 16 43 H 0.06 0.63 7.71 -51.93 -0.40 0.23 16 44 H 0.03 0.58 8.14 -51.93 -0.42 0.15 16 45 H 0.02 0.51 8.13 -51.93 -0.42 0.09 16 46 H 0.39 10.98 7.61 -40.82 -0.31 10.67 16 47 H 0.25 7.66 8.31 -40.82 -0.34 7.32 16 48 H 0.08 0.78 7.71 -51.93 -0.40 0.38 16 49 H 0.08 0.97 7.43 -51.93 -0.39 0.59 16 50 H 0.08 0.58 8.12 -51.93 -0.42 0.16 16 51 H 0.07 0.54 8.14 -51.93 -0.42 0.12 16 52 H 0.09 0.54 5.71 -51.93 -0.30 0.24 16 53 H 0.07 0.51 7.46 -51.93 -0.39 0.12 16 LS Contribution 409.19 15.07 6.17 6.17 Total: -1.00 -34.85 409.19 -9.58 -44.43 By element: Atomic # 1 Polarization: 14.22 SS G_CDS: -9.39 Total: 4.82 kcal Atomic # 6 Polarization: -3.38 SS G_CDS: -0.93 Total: -4.31 kcal Atomic # 7 Polarization: -54.12 SS G_CDS: -0.10 Total: -54.22 kcal Atomic # 8 Polarization: -14.45 SS G_CDS: -0.31 Total: -14.76 kcal Atomic # 14 Polarization: 22.89 SS G_CDS: -5.02 Total: 17.86 kcal Total LS contribution 6.17 Total: 6.17 kcal Total: -34.85 -9.58 -44.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. ZINC001276912641.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 -19.368 kcal (2) G-P(sol) polarization free energy of solvation -34.848 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -54.215 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.585 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.433 kcal (6) G-S(sol) free energy of system = (1) + (5) -63.800 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.89 seconds