Wall clock time and date at job start Tue Mar 31 2020 16:51: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.52997 * 1 3 3 C 1.53006 * 109.47401 * 2 1 4 4 C 1.52996 * 109.46898 * 240.00356 * 2 1 3 5 5 C 1.52995 * 109.47292 * 120.00190 * 2 1 3 6 6 C 1.54419 * 129.01550 * 59.99829 * 5 2 1 7 7 C 1.54413 * 78.89400 * 180.64455 * 6 5 2 8 8 C 1.50699 * 129.00845 * 180.64375 * 7 6 5 9 9 O 1.21281 * 120.00345 * 359.97438 * 8 7 6 10 10 N 1.34778 * 119.99885 * 180.02562 * 8 7 6 11 11 C 1.47851 * 134.36377 * 180.02562 * 10 8 7 12 12 C 1.52564 * 85.66371 * 156.34532 * 11 10 8 13 13 C 1.47856 * 134.36190 * 359.96979 * 10 8 7 14 14 C 1.53852 * 115.93570 * 265.03424 * 12 11 10 15 15 C 1.55397 * 101.30364 * 196.14659 * 14 12 11 16 16 N 1.47538 * 104.88714 * 35.97382 * 15 14 12 17 17 C 1.36939 * 125.66309 * 155.98205 * 16 15 14 18 18 H 1.08004 * 120.00116 * 192.96608 * 17 16 15 19 19 C 1.38812 * 119.99588 * 12.95699 * 17 16 15 20 20 N 1.31623 * 119.99875 * 204.39639 * 19 17 16 21 21 Si 1.86800 * 120.00059 * 24.40410 * 19 17 16 22 22 H 1.48497 * 109.46967 * 269.99382 * 21 19 17 23 23 H 1.48499 * 109.46957 * 29.99382 * 21 19 17 24 24 H 1.48501 * 109.46777 * 149.99591 * 21 19 17 25 25 C 1.47361 * 108.67965 * 335.71831 * 16 15 14 26 26 C 1.54891 * 84.49440 * 43.15348 * 7 6 5 27 27 C 1.54886 * 84.49553 * 318.13385 * 7 6 5 28 28 H 1.09005 * 109.46669 * 299.99947 * 1 2 3 29 29 H 1.08998 * 109.47177 * 59.99109 * 1 2 3 30 30 H 1.08998 * 109.47422 * 180.02562 * 1 2 3 31 31 H 1.08998 * 109.47320 * 180.02562 * 3 2 1 32 32 H 1.08992 * 109.47007 * 300.00631 * 3 2 1 33 33 H 1.09004 * 109.46773 * 59.99988 * 3 2 1 34 34 H 1.08999 * 109.47423 * 299.99344 * 4 2 1 35 35 H 1.08996 * 109.47504 * 60.00243 * 4 2 1 36 36 H 1.09009 * 109.47241 * 179.97438 * 4 2 1 37 37 H 1.09002 * 114.92666 * 68.15452 * 6 5 2 38 38 H 1.08998 * 114.92706 * 293.13718 * 6 5 2 39 39 H 1.08994 * 113.84247 * 42.16055 * 11 10 8 40 40 H 1.09005 * 113.95506 * 270.54484 * 11 10 8 41 41 H 1.09005 * 113.83772 * 89.46771 * 13 10 8 42 42 H 1.08991 * 113.84006 * 317.83695 * 13 10 8 43 43 H 1.09005 * 111.03576 * 314.38077 * 14 12 11 44 44 H 1.09007 * 111.03594 * 78.26563 * 14 12 11 45 45 H 1.08995 * 110.48631 * 277.08089 * 15 14 12 46 46 H 1.09007 * 110.30390 * 154.77005 * 15 14 12 47 47 H 0.96996 * 120.00134 * 354.58443 * 20 19 17 48 48 H 1.09005 * 109.98347 * 121.64103 * 25 16 15 49 49 H 1.08996 * 109.98663 * 242.94002 * 25 16 15 50 50 H 1.08998 * 114.97425 * 69.39129 * 26 7 6 51 51 H 1.08997 * 114.97673 * 204.56286 * 26 7 6 52 52 H 1.09003 * 114.97842 * 155.43668 * 27 7 6 53 53 H 1.08993 * 114.98177 * 290.60853 * 27 7 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.0401 1.4425 0.0000 4 6 2.0399 -0.7212 -1.2493 5 6 2.0400 -0.7213 1.2492 6 6 1.7984 -2.1794 1.6965 7 6 2.6731 -1.6537 2.8554 8 6 3.1648 -2.3831 4.0790 9 8 2.8743 -3.5484 4.2481 10 7 3.9267 -1.7402 4.9860 11 6 4.5647 -2.1290 6.2619 12 6 5.5311 -0.9789 5.9954 13 6 4.4582 -0.3652 5.1011 14 6 6.8506 -1.3562 5.2999 15 6 7.7131 -0.0937 5.5776 16 7 7.3185 0.3450 6.9298 17 6 8.0848 1.0957 7.7810 18 1 7.8030 1.1914 8.8192 19 6 9.2240 1.7334 7.3093 20 7 10.1996 2.0369 8.1392 21 14 9.3724 2.1612 5.4970 22 1 10.0066 1.0311 4.7720 23 1 8.0209 2.4133 4.9357 24 1 10.2076 3.3794 5.3428 25 6 5.9603 -0.1532 7.2101 26 6 1.7757 -0.3946 2.7637 27 6 3.4908 -0.9714 1.7307 28 1 -0.3633 0.5139 0.8900 29 1 -0.3633 0.5140 -0.8899 30 1 -0.3634 -1.0276 0.0005 31 1 3.1300 1.4425 -0.0005 32 1 1.6767 1.9564 -0.8899 33 1 1.6767 1.9563 0.8900 34 1 1.6767 -1.7489 -1.2493 35 1 1.6766 -0.2073 -2.1392 36 1 3.1300 -0.7208 -1.2496 37 1 2.2713 -2.9362 1.0706 38 1 0.7686 -2.4122 1.9671 39 1 5.0334 -3.1127 6.2383 40 1 3.9321 -1.9791 7.1368 41 1 3.7916 0.3227 5.6213 42 1 4.8431 0.0381 4.1646 43 1 7.2919 -2.2439 5.7532 44 1 6.7019 -1.4991 4.2296 45 1 7.4945 0.6856 4.8476 46 1 8.7730 -0.3475 5.5560 47 1 10.0908 1.8867 9.0913 48 1 5.9700 -0.7799 8.1019 49 1 5.2777 0.6849 7.3508 50 1 0.7450 -0.5396 3.0873 51 1 2.2303 0.5228 3.1374 52 1 4.0525 -0.0899 2.0400 53 1 4.0550 -1.6525 1.0938 RHF calculation, no. of doubly occupied orbitals= 64 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 ZINC001073476120.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 16:51:21 Heat of formation + Delta-G solvation = 115.624199 kcal Electronic energy + Delta-G solvation = -29868.682715 eV Core-core repulsion = 26131.646401 eV Total energy + Delta-G solvation = -3737.036314 eV No. of doubly occupied orbitals = 64 Molecular weight (most abundant/longest-lived isotopes) = 332.221 amu Computer time = 2.18 seconds Orbital eigenvalues (eV) -42.43659 -40.51255 -38.11710 -36.24709 -35.89412 -32.39714 -30.11326 -29.01212 -27.80078 -27.77094 -26.82699 -26.30278 -25.72804 -25.61760 -24.61925 -23.17123 -21.46626 -20.68455 -18.74392 -18.49200 -17.85891 -17.30405 -16.77371 -16.20969 -15.77900 -15.09632 -14.80621 -14.64522 -14.54407 -14.43099 -14.37187 -14.19617 -13.71926 -13.40453 -13.34091 -13.05247 -12.66332 -12.63935 -12.61838 -12.26451 -12.01280 -11.86853 -11.64149 -11.57125 -11.53915 -11.39125 -11.24129 -11.15591 -11.08364 -11.06188 -10.98603 -10.62849 -10.42757 -10.30870 -10.26972 -10.06171 -9.63100 -9.59210 -9.26895 -8.68824 -8.52559 -7.10899 -6.33763 -3.22801 2.43738 2.70116 3.38957 3.48300 3.64014 3.85434 3.97223 4.05683 4.13873 4.27249 4.36446 4.42926 4.54716 4.57108 4.77157 4.78707 4.79987 4.88245 5.06786 5.19853 5.29742 5.30572 5.32855 5.35923 5.36347 5.40280 5.42302 5.46827 5.50743 5.53887 5.58316 5.59240 5.63860 5.65027 5.66188 5.71571 5.72961 5.77285 5.84013 5.84772 5.95181 6.03549 6.04628 6.07564 6.14841 6.24949 6.44233 6.52302 6.57769 6.71605 7.64646 7.70024 8.13938 8.32081 9.01582 9.56329 10.23547 11.14165 Molecular weight = 332.22amu Principal moments of inertia in cm(-1) A = 0.014955 B = 0.003111 C = 0.002737 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1871.776269 B = 8998.843325 C =10227.801784 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.140 4.140 2 C -0.005 4.005 3 C -0.141 4.141 4 C -0.141 4.141 5 C -0.104 4.104 6 C -0.087 4.087 7 C -0.158 4.158 8 C 0.594 3.406 9 O -0.542 6.542 10 N -0.634 5.634 11 C 0.125 3.875 12 C -0.105 4.105 13 C 0.120 3.880 14 C -0.104 4.104 15 C 0.106 3.894 16 N -0.608 5.608 17 C -0.327 4.327 18 H 0.075 0.925 19 C -0.020 4.020 20 N -0.922 5.922 21 Si 1.065 2.935 22 H -0.294 1.294 23 H -0.287 1.287 24 H -0.299 1.299 25 C 0.151 3.849 26 C -0.099 4.099 27 C -0.097 4.097 28 H 0.056 0.944 29 H 0.057 0.943 30 H 0.056 0.944 31 H 0.058 0.942 32 H 0.057 0.943 33 H 0.057 0.943 34 H 0.056 0.944 35 H 0.058 0.942 36 H 0.057 0.943 37 H 0.075 0.925 38 H 0.075 0.925 39 H 0.083 0.917 40 H 0.087 0.913 41 H 0.089 0.911 42 H 0.083 0.917 43 H 0.076 0.924 44 H 0.067 0.933 45 H 0.031 0.969 46 H 0.061 0.939 47 H 0.263 0.737 48 H 0.037 0.963 49 H 0.032 0.968 50 H 0.075 0.925 51 H 0.079 0.921 52 H 0.080 0.920 53 H 0.075 0.925 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -19.204 -4.348 -16.987 26.005 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.197 4.197 2 C -0.005 4.005 3 C -0.198 4.198 4 C -0.198 4.198 5 C -0.104 4.104 6 C -0.125 4.125 7 C -0.161 4.161 8 C 0.383 3.617 9 O -0.420 6.420 10 N -0.369 5.369 11 C 0.003 3.997 12 C -0.107 4.107 13 C -0.003 4.003 14 C -0.142 4.142 15 C -0.019 4.019 16 N -0.335 5.335 17 C -0.454 4.454 18 H 0.093 0.907 19 C -0.221 4.221 20 N -0.563 5.563 21 Si 0.897 3.103 22 H -0.222 1.222 23 H -0.214 1.214 24 H -0.227 1.227 25 C 0.025 3.975 26 C -0.137 4.137 27 C -0.134 4.134 28 H 0.075 0.925 29 H 0.077 0.923 30 H 0.075 0.925 31 H 0.077 0.923 32 H 0.076 0.924 33 H 0.076 0.924 34 H 0.075 0.925 35 H 0.077 0.923 36 H 0.076 0.924 37 H 0.093 0.907 38 H 0.093 0.907 39 H 0.102 0.898 40 H 0.105 0.895 41 H 0.107 0.893 42 H 0.101 0.899 43 H 0.094 0.906 44 H 0.085 0.915 45 H 0.049 0.951 46 H 0.079 0.921 47 H 0.073 0.927 48 H 0.055 0.945 49 H 0.050 0.950 50 H 0.093 0.907 51 H 0.097 0.903 52 H 0.098 0.902 53 H 0.093 0.907 Dipole moment (debyes) X Y Z Total from point charges -18.408 -4.385 -18.079 26.171 hybrid contribution -1.076 -0.607 1.353 1.832 sum -19.483 -4.992 -16.726 26.158 Atomic orbital electron populations 1.21824 0.94010 1.02003 1.01891 1.19535 0.95221 0.94573 0.91215 1.21834 1.01425 0.94653 1.01921 1.21828 1.01353 1.00313 0.96296 1.24331 0.96360 0.96056 0.93611 1.23765 0.99020 0.95719 0.93994 1.23715 0.99860 0.98338 0.94151 1.18297 0.76234 0.85479 0.81687 1.90694 1.57373 1.16521 1.77386 1.50044 1.50608 1.14061 1.22144 1.23046 0.91773 0.99320 0.85569 1.22470 0.94575 0.95943 0.97677 1.23173 0.92848 0.84429 0.99838 1.22446 0.93447 0.99205 0.99134 1.22125 0.98316 0.93269 0.88185 1.44875 1.12437 1.61734 1.14484 1.20027 1.07185 1.27131 0.91098 0.90710 1.25875 0.91973 0.97503 1.06749 1.70904 1.28617 1.52311 1.04424 0.83836 0.74822 0.76181 0.75425 1.22165 1.21374 1.22712 1.21028 0.87941 0.94967 0.93580 1.24531 0.98744 0.95552 0.94833 1.24074 0.95415 0.99551 0.94370 0.92469 0.92349 0.92505 0.92315 0.92375 0.92398 0.92491 0.92344 0.92392 0.90667 0.90681 0.89812 0.89457 0.89286 0.89886 0.90575 0.91472 0.95106 0.92113 0.92693 0.94549 0.95033 0.90669 0.90283 0.90162 0.90655 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 16. 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.14 -0.88 8.47 37.16 0.31 -0.57 16 2 C -0.01 -0.04 1.02 -154.51 -0.16 -0.20 16 3 C -0.14 -0.95 8.47 37.16 0.31 -0.63 16 4 C -0.14 -0.91 8.47 37.16 0.31 -0.60 16 5 C -0.10 -0.87 1.67 -169.83 -0.28 -1.15 16 6 C -0.09 -0.88 7.78 -26.71 -0.21 -1.08 16 7 C -0.16 -1.81 3.10 -172.26 -0.53 -2.34 16 8 C 0.59 8.78 8.17 -10.99 -0.09 8.69 16 9 O -0.54 -9.93 17.82 5.56 0.10 -9.83 16 10 N -0.63 -8.65 3.62 -185.64 -0.67 -9.32 16 11 C 0.12 1.68 8.16 -3.53 -0.03 1.65 16 12 C -0.11 -1.57 0.97 -154.72 -0.15 -1.72 16 13 C 0.12 1.52 7.85 -3.53 -0.03 1.49 16 14 C -0.10 -1.67 6.43 -25.02 -0.16 -1.84 16 15 C 0.11 2.19 4.72 -2.20 -0.01 2.17 16 16 N -0.61 -14.19 3.28 -170.20 -0.56 -14.75 16 17 C -0.33 -9.10 10.45 -15.06 -0.16 -9.25 16 18 H 0.08 2.20 8.06 -52.48 -0.42 1.78 16 19 C -0.02 -0.59 5.42 -17.43 -0.09 -0.69 16 20 N -0.92 -28.47 13.70 55.49 0.76 -27.71 16 21 Si 1.07 25.97 24.69 -169.99 -4.20 21.77 16 22 H -0.29 -6.79 6.62 56.52 0.37 -6.41 16 23 H -0.29 -6.52 6.21 56.52 0.35 -6.17 16 24 H -0.30 -7.43 7.11 56.52 0.40 -7.03 16 25 C 0.15 2.79 6.75 -3.56 -0.02 2.77 16 26 C -0.10 -0.93 7.95 -24.71 -0.20 -1.12 16 27 C -0.10 -0.94 7.84 -25.98 -0.20 -1.14 16 28 H 0.06 0.37 8.14 -51.93 -0.42 -0.06 16 29 H 0.06 0.32 8.14 -51.93 -0.42 -0.10 16 30 H 0.06 0.36 8.14 -51.93 -0.42 -0.06 16 31 H 0.06 0.42 8.14 -51.93 -0.42 -0.01 16 32 H 0.06 0.34 8.14 -51.93 -0.42 -0.08 16 33 H 0.06 0.40 8.14 -51.93 -0.42 -0.02 16 34 H 0.06 0.37 8.14 -51.93 -0.42 -0.05 16 35 H 0.06 0.33 8.14 -51.93 -0.42 -0.10 16 36 H 0.06 0.39 8.14 -51.92 -0.42 -0.03 16 37 H 0.07 0.76 7.94 -51.93 -0.41 0.35 16 38 H 0.07 0.74 7.91 -51.93 -0.41 0.33 16 39 H 0.08 1.14 8.12 -51.93 -0.42 0.72 16 40 H 0.09 1.06 8.14 -51.93 -0.42 0.64 16 41 H 0.09 1.03 8.00 -51.93 -0.42 0.61 16 42 H 0.08 1.01 7.63 -51.93 -0.40 0.61 16 43 H 0.08 1.21 8.14 -51.93 -0.42 0.79 16 44 H 0.07 0.97 8.06 -51.93 -0.42 0.55 16 45 H 0.03 0.64 5.55 -51.93 -0.29 0.35 16 46 H 0.06 1.36 6.25 -51.92 -0.32 1.03 16 47 H 0.26 8.03 8.73 -40.82 -0.36 7.68 16 48 H 0.04 0.68 8.09 -51.93 -0.42 0.26 16 49 H 0.03 0.58 7.86 -51.93 -0.41 0.17 16 50 H 0.07 0.69 7.91 -51.93 -0.41 0.27 16 51 H 0.08 0.71 7.97 -51.93 -0.41 0.30 16 52 H 0.08 0.76 7.48 -51.93 -0.39 0.38 16 53 H 0.07 0.73 7.93 -51.93 -0.41 0.32 16 LS Contribution 409.78 15.07 6.18 6.18 Total: -1.00 -32.59 409.78 -9.62 -42.20 By element: Atomic # 1 Polarization: 6.85 SS G_CDS: -9.84 Total: -2.99 kcal Atomic # 6 Polarization: -4.17 SS G_CDS: -1.38 Total: -5.55 kcal Atomic # 7 Polarization: -51.30 SS G_CDS: -0.47 Total: -51.77 kcal Atomic # 8 Polarization: -9.93 SS G_CDS: 0.10 Total: -9.83 kcal Atomic # 14 Polarization: 25.97 SS G_CDS: -4.20 Total: 21.77 kcal Total LS contribution 6.18 Total: 6.18 kcal Total: -32.59 -9.62 -42.20 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. ZINC001073476120.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 157.827 kcal (2) G-P(sol) polarization free energy of solvation -32.586 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system 125.241 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.617 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.203 kcal (6) G-S(sol) free energy of system = (1) + (5) 115.624 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 2.18 seconds