Wall clock time and date at job start Sat Apr 11 2020 05:30:05 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.46901 * 1 3 3 C 1.46892 * 111.00306 * 2 1 4 4 C 1.50701 * 109.46993 * 66.04420 * 3 2 1 5 5 H 1.07994 * 120.00328 * 0.02562 * 4 3 2 6 6 C 1.37887 * 119.99865 * 179.97438 * 4 3 2 7 7 N 1.31513 * 120.00027 * 0.02562 * 6 4 3 8 8 Si 1.86798 * 119.99584 * 180.02562 * 6 4 3 9 9 H 1.48498 * 109.47098 * 60.00416 * 8 6 4 10 10 H 1.48504 * 109.46649 * 180.02562 * 8 6 4 11 11 H 1.48506 * 109.47435 * 300.00203 * 8 6 4 12 12 C 1.46906 * 110.99925 * 236.04735 * 2 1 3 13 13 C 1.52998 * 109.47197 * 199.76710 * 12 2 1 14 14 H 1.09003 * 109.47068 * 302.16128 * 13 12 2 15 15 O 1.42905 * 109.47031 * 62.15935 * 13 12 2 16 16 C 1.52999 * 109.47318 * 182.15693 * 13 12 2 17 17 N 1.46901 * 109.47409 * 179.97438 * 16 13 12 18 18 C 1.46893 * 111.00024 * 180.02562 * 17 16 13 19 19 C 1.53005 * 109.47093 * 179.97438 * 18 17 16 20 20 C 1.53000 * 109.46854 * 300.00154 * 19 18 17 21 21 C 1.52999 * 109.47054 * 59.99907 * 19 18 17 22 22 C 1.50697 * 109.47039 * 179.97438 * 19 18 17 23 23 C 1.38269 * 119.98742 * 359.73254 * 22 19 18 24 24 C 1.38231 * 120.05006 * 179.72687 * 23 22 19 25 25 C 1.38219 * 120.02346 * 0.54128 * 24 23 22 26 26 C 1.38447 * 119.97644 * 359.72593 * 25 24 23 27 27 F 1.35101 * 120.02494 * 180.02562 * 26 25 24 28 28 C 1.38174 * 119.98576 * 179.97438 * 22 19 18 29 29 H 1.08999 * 109.47031 * 54.96694 * 1 2 3 30 30 H 1.09002 * 109.47006 * 174.96219 * 1 2 3 31 31 H 1.08993 * 109.46940 * 294.96483 * 1 2 3 32 32 H 1.08994 * 109.47530 * 186.04788 * 3 2 1 33 33 H 1.09008 * 109.47253 * 306.04661 * 3 2 1 34 34 H 0.97000 * 119.99606 * 179.97438 * 7 6 4 35 35 H 1.09000 * 109.46958 * 319.76806 * 12 2 1 36 36 H 1.09002 * 109.47247 * 79.76447 * 12 2 1 37 37 H 0.96701 * 114.00201 * 59.99801 * 15 13 12 38 38 H 1.08995 * 109.47044 * 300.00138 * 16 13 12 39 39 H 1.09002 * 109.46836 * 59.99632 * 16 13 12 40 40 H 1.00901 * 110.99676 * 56.04861 * 17 16 13 41 41 H 1.09006 * 109.47427 * 300.00241 * 18 17 16 42 42 H 1.08997 * 109.47301 * 59.99976 * 18 17 16 43 43 H 1.08997 * 109.47385 * 60.00219 * 20 19 18 44 44 H 1.09005 * 109.46246 * 179.97438 * 20 19 18 45 45 H 1.08997 * 109.47025 * 299.99704 * 20 19 18 46 46 H 1.09006 * 109.46727 * 60.00206 * 21 19 18 47 47 H 1.08995 * 109.47368 * 179.97438 * 21 19 18 48 48 H 1.09001 * 109.47120 * 300.00336 * 21 19 18 49 49 H 1.07991 * 119.97873 * 359.97014 * 23 22 19 50 50 H 1.07997 * 119.98918 * 180.26947 * 24 23 22 51 51 H 1.08002 * 120.01381 * 179.75260 * 25 24 23 52 52 H 1.08000 * 120.01467 * 359.97438 * 28 22 19 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.9955 1.3713 0.0000 4 6 1.6370 2.0470 1.2984 5 1 1.0746 1.5112 2.0488 6 6 2.0278 3.3491 1.5289 7 7 2.7130 4.0015 0.6154 8 14 1.5840 4.1863 3.1387 9 1 2.1757 3.4285 4.2704 10 1 2.1161 5.5728 3.1386 11 1 0.1067 4.2236 3.2865 12 6 1.9955 -0.7660 -1.1377 13 6 3.4455 -1.1649 -0.8566 14 1 4.0400 -0.2695 -0.6748 15 8 3.4925 -2.0098 0.2950 16 6 4.0109 -1.9161 -2.0636 17 7 5.4034 -2.2986 -1.7941 18 6 5.9805 -3.0231 -2.9342 19 6 7.4282 -3.4033 -2.6169 20 6 8.2425 -2.1350 -2.3535 21 6 7.4605 -4.2970 -1.3756 22 6 8.0200 -4.1471 -3.7864 23 6 7.2540 -4.3863 -4.9124 24 6 7.7973 -5.0635 -5.9880 25 6 9.1036 -5.5120 -5.9346 26 6 9.8698 -5.2779 -4.8055 27 9 11.1469 -5.7154 -4.7525 28 6 9.3263 -4.5940 -3.7309 29 1 -0.3633 0.5899 -0.8415 30 1 -0.3633 -1.0237 -0.0902 31 1 -0.3633 0.4337 0.9316 32 1 3.0797 1.3443 -0.1083 33 1 1.5612 1.9273 -0.8310 34 1 2.9876 4.9176 0.7775 35 1 1.3932 -1.6630 -1.2815 36 1 1.9552 -0.1531 -2.0382 37 1 2.9842 -2.8273 0.2040 38 1 3.4165 -2.8115 -2.2453 39 1 3.9750 -1.2717 -2.9420 40 1 5.4714 -2.8432 -0.9473 41 1 5.4001 -3.9264 -3.1222 42 1 5.9579 -2.3862 -3.8184 43 1 7.8143 -1.5974 -1.5075 44 1 9.2738 -2.4060 -2.1271 45 1 8.2195 -1.4983 -3.2379 46 1 6.8803 -5.2006 -1.5633 47 1 8.4917 -4.5683 -1.1498 48 1 7.0322 -3.7594 -0.5296 49 1 6.2318 -4.0405 -4.9531 50 1 7.1995 -5.2461 -6.8687 51 1 9.5267 -6.0446 -6.7736 52 1 9.9233 -4.4101 -2.8499 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) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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 ZINC000566114724.mol2 52 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 05:30:05 Heat of formation + Delta-G solvation = -73.785166 kcal Electronic energy + Delta-G solvation = -29181.327974 eV Core-core repulsion = 25094.850418 eV Total energy + Delta-G solvation = -4086.477556 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 338.211 amu Computer time = 0.69 seconds Orbital eigenvalues (eV) -49.19823 -40.24614 -38.29966 -37.05715 -35.55374 -33.84515 -32.45525 -31.09791 -30.30536 -29.51625 -27.42950 -27.31866 -26.58166 -24.88560 -23.53602 -22.93904 -22.45709 -21.47370 -20.07942 -19.52853 -17.71766 -16.75510 -16.53215 -16.41177 -16.11576 -15.94349 -15.88995 -15.19393 -14.92228 -14.81599 -14.56039 -14.32513 -14.06798 -13.82743 -13.76663 -13.46597 -13.33481 -13.04600 -12.76726 -12.50405 -12.33361 -12.13877 -12.08551 -11.94319 -11.91180 -11.74293 -11.38248 -11.16587 -11.13451 -10.94557 -10.91653 -10.78935 -10.69623 -10.60850 -10.34713 -9.73356 -9.60876 -9.28309 -9.11705 -9.00620 -8.42829 -8.38002 -7.18594 -5.98256 -3.98760 0.81674 0.85961 3.23143 3.36918 3.40612 3.44146 3.75263 4.03784 4.28913 4.36161 4.47027 4.65955 4.68112 4.72939 4.79593 4.82955 4.95276 4.96094 5.03317 5.26506 5.34337 5.38267 5.38723 5.45344 5.51030 5.55825 5.57321 5.67032 5.70968 5.73127 5.75471 5.78277 5.80852 5.84167 5.96683 6.03543 6.05362 6.08657 6.20489 6.31051 6.35731 6.54471 6.55522 6.71932 6.94166 7.30689 7.34113 7.37640 7.83196 8.06821 8.13547 8.23754 8.58916 9.08042 9.64834 11.12751 Molecular weight = 338.21amu Principal moments of inertia in cm(-1) A = 0.020238 B = 0.001879 C = 0.001806 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1383.169847 B =14898.312962 C =15496.310899 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.041 3.959 2 N -0.531 5.531 3 C 0.193 3.807 4 C -0.505 4.505 5 H 0.062 0.938 6 C 0.054 3.946 7 N -0.899 5.899 8 Si 0.898 3.102 9 H -0.276 1.276 10 H -0.288 1.288 11 H -0.278 1.278 12 C 0.026 3.974 13 C 0.116 3.884 14 H 0.112 0.888 15 O -0.541 6.541 16 C 0.030 3.970 17 N -0.680 5.680 18 C 0.068 3.932 19 C -0.024 4.024 20 C -0.132 4.132 21 C -0.144 4.144 22 C -0.035 4.035 23 C -0.133 4.133 24 C -0.088 4.088 25 C -0.163 4.163 26 C 0.099 3.901 27 F -0.142 7.142 28 C -0.149 4.149 29 H 0.011 0.989 30 H 0.041 0.959 31 H 0.072 0.928 32 H 0.032 0.968 33 H -0.012 1.012 34 H 0.261 0.739 35 H 0.058 0.942 36 H 0.032 0.968 37 H 0.358 0.642 38 H 0.027 0.973 39 H 0.078 0.922 40 H 0.349 0.651 41 H 0.030 0.970 42 H 0.075 0.925 43 H 0.079 0.921 44 H 0.049 0.951 45 H 0.055 0.945 46 H 0.054 0.946 47 H 0.059 0.941 48 H 0.071 0.929 49 H 0.132 0.868 50 H 0.132 0.868 51 H 0.132 0.868 52 H 0.135 0.865 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 10.169 -20.898 -13.292 26.773 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.108 4.108 2 N -0.255 5.255 3 C 0.068 3.932 4 C -0.543 4.543 5 H 0.080 0.920 6 C -0.147 4.147 7 N -0.539 5.539 8 Si 0.726 3.274 9 H -0.201 1.201 10 H -0.214 1.214 11 H -0.204 1.204 12 C -0.104 4.104 13 C 0.057 3.943 14 H 0.129 0.871 15 O -0.343 6.343 16 C -0.100 4.100 17 N -0.313 5.313 18 C -0.062 4.062 19 C -0.024 4.024 20 C -0.189 4.189 21 C -0.201 4.201 22 C -0.035 4.035 23 C -0.151 4.151 24 C -0.106 4.106 25 C -0.182 4.182 26 C 0.080 3.920 27 F -0.120 7.120 28 C -0.167 4.167 29 H 0.029 0.971 30 H 0.060 0.940 31 H 0.090 0.910 32 H 0.050 0.950 33 H 0.006 0.994 34 H 0.071 0.929 35 H 0.076 0.924 36 H 0.050 0.950 37 H 0.203 0.797 38 H 0.045 0.955 39 H 0.096 0.904 40 H 0.169 0.831 41 H 0.048 0.952 42 H 0.093 0.907 43 H 0.098 0.902 44 H 0.068 0.932 45 H 0.074 0.926 46 H 0.073 0.927 47 H 0.078 0.922 48 H 0.090 0.910 49 H 0.150 0.850 50 H 0.149 0.851 51 H 0.150 0.850 52 H 0.153 0.847 Dipole moment (debyes) X Y Z Total from point charges 11.092 -19.856 -12.845 26.121 hybrid contribution -1.025 -0.406 -0.326 1.149 sum 10.067 -20.263 -13.171 26.180 Atomic orbital electron populations 1.22191 0.88921 0.98495 1.01145 1.62044 1.06446 1.17902 1.39130 1.19778 0.95537 0.81374 0.96492 1.21111 1.34703 0.96812 1.01696 0.91988 1.24991 0.95360 0.97392 0.96963 1.69946 1.38418 1.11465 1.34090 0.86560 0.78580 0.78295 0.83953 1.20120 1.21409 1.20391 1.22677 0.95964 0.98013 0.93736 1.21984 0.95140 0.93576 0.83651 0.87066 1.86705 1.73803 1.34336 1.39505 1.22668 0.88378 0.98906 1.00047 1.60332 1.18081 1.47582 1.05334 1.22078 0.93127 0.96206 0.94754 1.19543 0.95025 0.95397 0.92472 1.21800 0.98337 0.95622 1.03126 1.21945 1.02367 0.98970 0.96855 1.19729 0.92715 0.95833 0.95255 1.21284 1.00731 0.99950 0.93183 1.21021 0.94252 0.96649 0.98655 1.21136 0.95312 1.03879 0.97845 1.17701 0.79486 0.98939 0.95903 1.91501 1.36088 1.87644 1.96807 1.20967 0.96345 1.00781 0.98647 0.97059 0.94016 0.90955 0.95009 0.99412 0.92934 0.92400 0.94968 0.79731 0.95493 0.90373 0.83062 0.95176 0.90688 0.90181 0.93176 0.92635 0.92724 0.92200 0.91026 0.84985 0.85050 0.84975 0.84727 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.04 0.80 9.31 60.07 0.56 1.35 16 2 N -0.53 -11.38 4.59 -244.37 -1.12 -12.50 16 3 C 0.19 5.06 4.55 -4.98 -0.02 5.04 16 4 C -0.50 -14.40 8.58 -34.44 -0.30 -14.70 16 5 H 0.06 1.76 7.10 -52.49 -0.37 1.39 16 6 C 0.05 1.55 5.46 -17.82 -0.10 1.46 16 7 N -0.90 -26.99 13.29 55.43 0.74 -26.25 16 8 Si 0.90 21.99 29.54 -169.99 -5.02 16.97 16 9 H -0.28 -6.69 7.11 56.52 0.40 -6.29 16 10 H -0.29 -6.97 7.11 56.52 0.40 -6.57 16 11 H -0.28 -6.66 7.11 56.52 0.40 -6.26 16 12 C 0.03 0.44 4.65 -3.82 -0.02 0.42 16 13 C 0.12 1.82 3.31 -26.73 -0.09 1.73 16 14 H 0.11 2.17 6.73 -51.93 -0.35 1.82 16 15 O -0.54 -8.51 12.58 -35.23 -0.44 -8.96 16 16 C 0.03 0.36 5.39 -3.82 -0.02 0.33 16 17 N -0.68 -7.54 5.71 -113.78 -0.65 -8.19 16 18 C 0.07 0.61 3.31 -3.83 -0.01 0.60 16 19 C -0.02 -0.21 0.67 -155.66 -0.10 -0.31 16 20 C -0.13 -1.25 8.13 37.16 0.30 -0.95 16 21 C -0.14 -1.16 7.96 37.16 0.30 -0.87 16 22 C -0.03 -0.29 4.76 -104.64 -0.50 -0.79 16 23 C -0.13 -1.03 8.97 -39.57 -0.36 -1.39 16 24 C -0.09 -0.61 10.04 -39.59 -0.40 -1.01 16 25 C -0.16 -1.26 10.03 -39.51 -0.40 -1.66 16 26 C 0.10 0.91 7.30 -39.41 -0.29 0.62 16 27 F -0.14 -1.43 17.27 2.25 0.04 -1.39 16 28 C -0.15 -1.31 8.69 -39.51 -0.34 -1.65 16 29 H 0.01 0.21 8.08 -51.93 -0.42 -0.21 16 30 H 0.04 0.67 7.85 -51.93 -0.41 0.26 16 31 H 0.07 1.49 7.01 -51.93 -0.36 1.13 16 32 H 0.03 0.89 6.60 -51.93 -0.34 0.55 16 33 H -0.01 -0.33 8.07 -51.92 -0.42 -0.75 16 34 H 0.26 7.43 8.31 -40.82 -0.34 7.09 16 35 H 0.06 0.83 7.75 -51.93 -0.40 0.42 16 36 H 0.03 0.56 8.00 -51.93 -0.42 0.14 16 37 H 0.36 4.18 9.05 45.56 0.41 4.60 16 38 H 0.03 0.28 8.14 -51.93 -0.42 -0.15 16 39 H 0.08 0.93 8.14 -51.93 -0.42 0.51 16 40 H 0.35 3.85 6.16 -39.29 -0.24 3.60 16 41 H 0.03 0.25 7.50 -51.93 -0.39 -0.14 16 42 H 0.07 0.70 7.53 -51.93 -0.39 0.31 16 43 H 0.08 0.88 7.60 -51.93 -0.39 0.48 16 44 H 0.05 0.43 7.39 -51.93 -0.38 0.05 16 45 H 0.05 0.55 8.14 -51.93 -0.42 0.13 16 46 H 0.05 0.42 8.14 -51.93 -0.42 0.00 16 47 H 0.06 0.45 7.39 -51.93 -0.38 0.07 16 48 H 0.07 0.61 6.57 -51.93 -0.34 0.27 16 49 H 0.13 0.88 5.89 -52.49 -0.31 0.57 16 50 H 0.13 0.61 8.06 -52.49 -0.42 0.19 16 51 H 0.13 0.79 8.06 -52.49 -0.42 0.37 16 52 H 0.13 1.07 7.10 -52.49 -0.37 0.70 16 LS Contribution 411.77 15.07 6.21 6.21 Total: -1.00 -31.62 411.77 -10.00 -41.62 By element: Atomic # 1 Polarization: 12.24 SS G_CDS: -7.96 Total: 4.28 kcal Atomic # 6 Polarization: -9.99 SS G_CDS: -1.78 Total: -11.77 kcal Atomic # 7 Polarization: -45.91 SS G_CDS: -1.03 Total: -46.95 kcal Atomic # 8 Polarization: -8.51 SS G_CDS: -0.44 Total: -8.96 kcal Atomic # 9 Polarization: -1.43 SS G_CDS: 0.04 Total: -1.39 kcal Atomic # 14 Polarization: 21.99 SS G_CDS: -5.02 Total: 16.97 kcal Total LS contribution 6.21 Total: 6.21 kcal Total: -31.62 -10.00 -41.62 kcal The number of atoms in the molecule is 52 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. ZINC000566114724.mol2 52 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 -32.167 kcal (2) G-P(sol) polarization free energy of solvation -31.622 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -63.789 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.996 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.618 kcal (6) G-S(sol) free energy of system = (1) + (5) -73.785 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.69 seconds