Wall clock time and date at job start Sat Apr 11 2020 06:13:08 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 N 1.46902 * 1 3 3 C 1.46898 * 110.99543 * 2 1 4 4 C 1.46894 * 110.99947 * 123.95688 * 2 1 3 5 5 C 1.52997 * 109.47515 * 66.04655 * 4 2 1 6 6 C 1.53000 * 109.47060 * 180.02562 * 5 4 2 7 7 C 1.50706 * 109.46735 * 180.02562 * 6 5 4 8 8 O 1.21273 * 120.00140 * 359.97438 * 7 6 5 9 9 N 1.34779 * 119.99767 * 180.02562 * 7 6 5 10 10 C 1.47029 * 125.64556 * 359.68841 * 9 7 6 11 11 C 1.54335 * 107.25970 * 178.85288 * 10 9 7 12 12 C 1.55157 * 103.02283 * 338.13530 * 11 10 9 13 13 H 1.08996 * 111.00938 * 277.23391 * 12 11 10 14 14 C 1.53006 * 111.03463 * 153.39519 * 12 11 10 15 15 N 1.46902 * 109.46761 * 72.79037 * 14 12 11 16 16 C 1.46897 * 110.99705 * 179.97438 * 15 14 12 17 17 C 1.50708 * 109.47031 * 179.97438 * 16 15 14 18 18 O 1.21286 * 119.99781 * 0.02562 * 17 16 15 19 19 N 1.34780 * 119.99835 * 180.02562 * 17 16 15 20 20 C 1.37096 * 119.99689 * 180.02562 * 19 17 16 21 21 H 1.07998 * 120.00110 * 359.97438 * 20 19 17 22 22 C 1.38945 * 120.00059 * 180.02562 * 20 19 17 23 23 N 1.31416 * 120.00034 * 359.97438 * 22 20 19 24 24 Si 1.86799 * 120.00143 * 179.97438 * 22 20 19 25 25 H 1.48502 * 109.47381 * 89.99833 * 24 22 20 26 26 H 1.48506 * 109.46950 * 209.99833 * 24 22 20 27 27 H 1.48505 * 109.47258 * 329.99373 * 24 22 20 28 28 C 1.47421 * 125.64901 * 179.97438 * 9 7 6 29 29 H 1.08997 * 109.47044 * 296.04083 * 1 2 3 30 30 H 1.09001 * 109.47499 * 56.04429 * 1 2 3 31 31 H 1.09003 * 109.47124 * 176.04281 * 1 2 3 32 32 H 1.09002 * 109.46576 * 183.95470 * 3 2 1 33 33 H 1.09002 * 109.47680 * 303.94979 * 3 2 1 34 34 H 1.08997 * 109.47173 * 63.95741 * 3 2 1 35 35 H 1.09002 * 109.47104 * 306.04387 * 4 2 1 36 36 H 1.08998 * 109.46828 * 186.05219 * 4 2 1 37 37 H 1.09002 * 109.47153 * 59.99762 * 5 4 2 38 38 H 1.09001 * 109.46832 * 299.99953 * 5 4 2 39 39 H 1.09001 * 109.47194 * 299.99650 * 6 5 4 40 40 H 1.08992 * 109.47399 * 59.99823 * 6 5 4 41 41 H 1.08994 * 109.88052 * 298.39407 * 10 9 7 42 42 H 1.08993 * 109.88202 * 59.56327 * 10 9 7 43 43 H 1.08999 * 110.84691 * 219.67348 * 11 10 9 44 44 H 1.08995 * 110.67211 * 96.52282 * 11 10 9 45 45 H 1.08998 * 109.47052 * 192.78312 * 14 12 11 46 46 H 1.08993 * 109.47480 * 312.78957 * 14 12 11 47 47 H 1.00901 * 110.99731 * 303.94924 * 15 14 12 48 48 H 1.08998 * 109.47327 * 299.99863 * 16 15 14 49 49 H 1.09000 * 109.47213 * 59.99839 * 16 15 14 50 50 H 0.96994 * 119.99991 * 359.97438 * 19 17 16 51 51 H 0.97000 * 120.00174 * 180.02562 * 23 22 20 52 52 H 1.09001 * 110.36259 * 322.92716 * 28 9 7 53 53 H 1.08995 * 110.36906 * 85.14076 * 28 9 7 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 7 1.4690 0.0000 0.0000 3 6 1.9953 1.3714 0.0000 4 6 1.9954 -0.7660 1.1375 5 6 1.6315 -2.2426 0.9703 6 6 2.1793 -3.0403 2.1554 7 6 1.8213 -4.4949 1.9903 8 8 1.1899 -4.8576 1.0206 9 7 2.2021 -5.3924 2.9209 10 6 2.9722 -5.0917 4.1367 11 6 3.1457 -6.4170 4.9084 12 6 1.9856 -7.2924 4.3650 13 1 1.0526 -7.0701 4.8829 14 6 2.3256 -8.7806 4.4683 15 7 2.2609 -9.2011 5.8743 16 6 2.5818 -10.6280 6.0122 17 6 2.5060 -11.0218 7.4649 18 8 2.2130 -10.1953 8.3028 19 7 2.7638 -12.2925 7.8330 20 6 2.6943 -12.6508 9.1545 21 1 2.4334 -11.9148 9.9006 22 6 2.9606 -13.9606 9.5341 23 7 3.2785 -14.8561 8.6263 24 14 2.8666 -14.4487 11.3347 25 1 1.4880 -14.8989 11.6543 26 1 3.8221 -15.5546 11.5983 27 1 3.2152 -13.2808 12.1831 28 6 1.9204 -6.8390 2.8854 29 1 -0.3633 0.4511 0.9233 30 1 -0.3634 0.5740 -0.8524 31 1 -0.3633 -1.0252 -0.0709 32 1 3.0827 1.3432 -0.0709 33 1 1.5898 1.9164 -0.8525 34 1 1.7043 1.8723 0.9233 35 1 1.5611 -0.3873 2.0627 36 1 3.0797 -0.6610 1.1756 37 1 2.0658 -2.6214 0.0450 38 1 0.5472 -2.3475 0.9322 39 1 3.2636 -2.9355 2.1939 40 1 1.7446 -2.6615 3.0803 41 1 3.9477 -4.6871 3.8670 42 1 2.4285 -4.3750 4.7521 43 1 3.0388 -6.2580 5.9814 44 1 4.1103 -6.8709 4.6814 45 1 1.6101 -9.3582 3.8830 46 1 3.3312 -8.9511 4.0839 47 1 2.8684 -8.6364 6.4489 48 1 1.8678 -11.2182 5.4377 49 1 3.5891 -10.8117 5.6384 50 1 2.9977 -12.9535 7.1629 51 1 3.4640 -15.7706 8.8912 52 1 0.9269 -7.0220 2.4759 53 1 2.6750 -7.3580 2.2944 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 ZINC001690733062.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:13:08 Heat of formation + Delta-G solvation = -73.643768 kcal Electronic energy + Delta-G solvation = -28503.962690 eV Core-core repulsion = 24410.099749 eV Total energy + Delta-G solvation = -4093.862941 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 1.04 seconds Orbital eigenvalues (eV) -40.41399 -38.58380 -38.18396 -37.72049 -35.39795 -33.88939 -33.65643 -31.96048 -31.19880 -29.75066 -28.35414 -28.01014 -26.95960 -24.99425 -23.91731 -23.08917 -22.83466 -21.19860 -20.52325 -19.38820 -19.23068 -17.90338 -17.40636 -16.94870 -16.44480 -16.04962 -15.76732 -15.54885 -15.26172 -15.15083 -14.81026 -14.62370 -14.46728 -14.03213 -13.93541 -13.65222 -13.40269 -13.21246 -13.00105 -12.86598 -12.72089 -12.65217 -12.44116 -12.24384 -12.13743 -12.05401 -11.85191 -11.66737 -11.57980 -11.52858 -11.22004 -11.06526 -11.01279 -10.91040 -10.83758 -10.40456 -10.23038 -10.08207 -9.52137 -8.92387 -8.83687 -8.48820 -8.09244 -7.97215 -6.42861 -3.83802 2.38174 2.65385 3.17085 3.25809 3.31362 3.48596 3.89278 3.97433 4.13160 4.28040 4.28873 4.48078 4.49828 4.61179 4.70209 4.75789 4.85486 4.92785 4.97520 4.98762 5.05236 5.08023 5.13059 5.16414 5.21695 5.22951 5.24292 5.25075 5.38318 5.39542 5.46712 5.50699 5.54984 5.61690 5.67135 5.74333 5.77803 5.88363 5.92107 5.96738 6.01601 6.13707 6.25653 6.53788 6.55932 7.16993 7.23470 7.54058 7.62667 7.90390 8.05066 8.59978 9.16654 9.53094 10.04377 10.74860 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.043275 B = 0.001304 C = 0.001285 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 646.862170 B =21475.353343 C =21778.449405 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.019 3.981 2 N -0.549 5.549 3 C 0.024 3.976 4 C 0.054 3.946 5 C -0.102 4.102 6 C -0.135 4.135 7 C 0.528 3.472 8 O -0.540 6.540 9 N -0.620 5.620 10 C 0.101 3.899 11 C -0.133 4.133 12 C -0.097 4.097 13 H 0.097 0.903 14 C 0.070 3.930 15 N -0.668 5.668 16 C 0.050 3.950 17 C 0.440 3.560 18 O -0.584 6.584 19 N -0.563 5.563 20 C -0.300 4.300 21 H 0.103 0.897 22 C 0.023 3.977 23 N -0.902 5.902 24 Si 0.930 3.070 25 H -0.278 1.278 26 H -0.284 1.284 27 H -0.270 1.270 28 C 0.111 3.889 29 H 0.025 0.975 30 H 0.067 0.933 31 H 0.074 0.926 32 H 0.068 0.932 33 H 0.068 0.932 34 H 0.025 0.975 35 H 0.035 0.965 36 H 0.075 0.925 37 H 0.079 0.921 38 H 0.071 0.929 39 H 0.093 0.907 40 H 0.092 0.908 41 H 0.069 0.931 42 H 0.070 0.930 43 H 0.090 0.910 44 H 0.076 0.924 45 H 0.075 0.925 46 H 0.030 0.970 47 H 0.349 0.651 48 H 0.087 0.913 49 H 0.043 0.957 50 H 0.392 0.608 51 H 0.272 0.728 52 H 0.080 0.920 53 H 0.068 0.932 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 0.822 26.440 -16.267 31.055 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.129 4.129 2 N -0.272 5.272 3 C -0.124 4.124 4 C -0.074 4.074 5 C -0.140 4.140 6 C -0.175 4.175 7 C 0.317 3.683 8 O -0.417 6.417 9 N -0.355 5.355 10 C -0.022 4.022 11 C -0.171 4.171 12 C -0.117 4.117 13 H 0.115 0.885 14 C -0.060 4.060 15 N -0.300 5.300 16 C -0.082 4.082 17 C 0.226 3.774 18 O -0.468 6.468 19 N -0.200 5.200 20 C -0.430 4.430 21 H 0.121 0.879 22 C -0.179 4.179 23 N -0.543 5.543 24 Si 0.757 3.243 25 H -0.204 1.204 26 H -0.209 1.209 27 H -0.194 1.194 28 C -0.012 4.012 29 H 0.044 0.956 30 H 0.086 0.914 31 H 0.093 0.907 32 H 0.087 0.913 33 H 0.087 0.913 34 H 0.043 0.957 35 H 0.053 0.947 36 H 0.093 0.907 37 H 0.098 0.902 38 H 0.090 0.910 39 H 0.111 0.889 40 H 0.110 0.890 41 H 0.087 0.913 42 H 0.088 0.912 43 H 0.108 0.892 44 H 0.095 0.905 45 H 0.094 0.906 46 H 0.049 0.951 47 H 0.168 0.832 48 H 0.105 0.895 49 H 0.061 0.939 50 H 0.226 0.774 51 H 0.081 0.919 52 H 0.098 0.902 53 H 0.087 0.913 Dipole moment (debyes) X Y Z Total from point charges -0.203 26.337 -16.697 31.184 hybrid contribution 0.889 -0.407 1.165 1.521 sum 0.686 25.929 -15.532 30.233 Atomic orbital electron populations 1.22470 0.87193 1.02811 1.00409 1.61942 1.08281 1.13778 1.43226 1.22388 1.00807 0.88725 1.00512 1.22169 1.00204 0.93758 0.91278 1.21343 1.02100 0.92633 0.97925 1.21605 1.03149 0.91609 1.01132 1.20345 0.78924 0.87561 0.81484 1.90719 1.40184 1.80823 1.30003 1.48547 1.55385 1.08500 1.23069 1.22071 0.95722 0.98555 0.85805 1.22771 0.99192 0.94227 1.00879 1.22384 0.99945 0.94767 0.94575 0.88464 1.21852 0.99632 0.93772 0.90736 1.60039 1.60099 1.03225 1.06648 1.21674 0.99670 0.89916 0.96960 1.19808 0.84544 0.83674 0.89399 1.90516 1.53104 1.50607 1.52597 1.41862 1.60617 1.08323 1.09182 1.19706 1.44185 0.95053 0.84033 0.87919 1.25257 0.97106 0.95514 0.99998 1.69613 1.47173 1.03752 1.33759 0.85862 0.78067 0.78267 0.82110 1.20375 1.20932 1.19423 1.22223 1.02054 0.80135 0.96774 0.95634 0.91399 0.90730 0.91305 0.91306 0.95688 0.94725 0.90657 0.90219 0.91037 0.88867 0.89004 0.91306 0.91174 0.89165 0.90504 0.90622 0.95145 0.83197 0.89528 0.93865 0.77449 0.91861 0.90165 0.91348 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.02 0.12 9.32 60.07 0.56 0.68 16 2 N -0.55 -3.93 5.02 -251.27 -1.26 -5.19 16 3 C 0.02 0.13 9.57 60.06 0.57 0.70 16 4 C 0.05 0.34 4.96 -3.83 -0.02 0.32 16 5 C -0.10 -0.77 4.43 -26.73 -0.12 -0.89 16 6 C -0.13 -0.91 4.89 -27.88 -0.14 -1.04 16 7 C 0.53 5.55 7.85 -10.98 -0.09 5.46 16 8 O -0.54 -7.80 16.24 5.56 0.09 -7.71 16 9 N -0.62 -5.83 3.33 -170.53 -0.57 -6.40 16 10 C 0.10 0.70 6.49 -3.05 -0.02 0.68 16 11 C -0.13 -1.05 5.90 -24.89 -0.15 -1.20 16 12 C -0.10 -0.94 3.25 -88.48 -0.29 -1.22 16 13 H 0.10 1.06 8.14 -51.93 -0.42 0.64 16 14 C 0.07 0.78 5.19 -3.82 -0.02 0.76 16 15 N -0.67 -10.13 6.13 -115.05 -0.70 -10.83 16 16 C 0.05 0.86 6.17 -4.97 -0.03 0.83 16 17 C 0.44 10.23 7.82 -10.98 -0.09 10.15 16 18 O -0.58 -15.62 15.78 5.55 0.09 -15.53 16 19 N -0.56 -14.13 5.29 -10.96 -0.06 -14.19 16 20 C -0.30 -8.32 9.68 -14.94 -0.14 -8.46 16 21 H 0.10 2.98 7.16 -52.49 -0.38 2.61 16 22 C 0.02 0.61 5.50 -17.47 -0.10 0.51 16 23 N -0.90 -24.38 13.05 55.50 0.72 -23.65 16 24 Si 0.93 20.89 29.55 -169.99 -5.02 15.86 16 25 H -0.28 -6.19 7.11 56.52 0.40 -5.78 16 26 H -0.28 -6.23 7.11 56.52 0.40 -5.83 16 27 H -0.27 -6.03 7.11 56.52 0.40 -5.63 16 28 C 0.11 1.09 6.37 -2.54 -0.02 1.07 16 29 H 0.03 0.15 8.14 -51.93 -0.42 -0.27 16 30 H 0.07 0.38 8.12 -51.93 -0.42 -0.04 16 31 H 0.07 0.52 6.47 -51.93 -0.34 0.18 16 32 H 0.07 0.33 8.12 -51.93 -0.42 -0.09 16 33 H 0.07 0.33 8.12 -51.93 -0.42 -0.09 16 34 H 0.02 0.12 8.14 -51.93 -0.42 -0.30 16 35 H 0.03 0.20 8.12 -51.93 -0.42 -0.22 16 36 H 0.07 0.41 8.06 -51.93 -0.42 -0.01 16 37 H 0.08 0.74 8.10 -51.93 -0.42 0.32 16 38 H 0.07 0.63 6.51 -51.93 -0.34 0.29 16 39 H 0.09 0.47 8.14 -51.93 -0.42 0.05 16 40 H 0.09 0.49 8.14 -51.93 -0.42 0.06 16 41 H 0.07 0.38 8.14 -51.93 -0.42 -0.05 16 42 H 0.07 0.41 8.14 -51.93 -0.42 -0.01 16 43 H 0.09 0.72 8.14 -51.93 -0.42 0.30 16 44 H 0.08 0.59 7.91 -51.93 -0.41 0.18 16 45 H 0.08 0.84 8.14 -51.93 -0.42 0.41 16 46 H 0.03 0.31 7.97 -51.93 -0.41 -0.10 16 47 H 0.35 5.36 7.94 -39.29 -0.31 5.05 16 48 H 0.09 1.38 8.14 -51.93 -0.42 0.96 16 49 H 0.04 0.66 8.14 -51.93 -0.42 0.24 16 50 H 0.39 9.40 7.90 -40.82 -0.32 9.08 16 51 H 0.27 6.95 8.31 -40.82 -0.34 6.61 16 52 H 0.08 0.88 8.04 -51.93 -0.42 0.46 16 53 H 0.07 0.65 8.04 -51.93 -0.42 0.23 16 LS Contribution 427.52 15.07 6.44 6.44 Total: -1.00 -33.61 427.52 -10.00 -43.60 By element: Atomic # 1 Polarization: 18.89 SS G_CDS: -9.65 Total: 9.24 kcal Atomic # 6 Polarization: 8.43 SS G_CDS: -0.07 Total: 8.36 kcal Atomic # 7 Polarization: -58.40 SS G_CDS: -1.87 Total: -60.27 kcal Atomic # 8 Polarization: -23.42 SS G_CDS: 0.18 Total: -23.24 kcal Atomic # 14 Polarization: 20.89 SS G_CDS: -5.02 Total: 15.86 kcal Total LS contribution 6.44 Total: 6.44 kcal Total: -33.61 -10.00 -43.60 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. ZINC001690733062.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 -30.039 kcal (2) G-P(sol) polarization free energy of solvation -33.608 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -63.647 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.997 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.605 kcal (6) G-S(sol) free energy of system = (1) + (5) -73.644 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.05 seconds