Wall clock time and date at job start Sat Apr 11 2020 06:14:17 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 * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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.53004 * 1 3 3 C 1.53004 * 109.47072 * 2 1 4 4 C 1.52993 * 109.47301 * 179.97438 * 3 2 1 5 5 O 1.42898 * 109.47476 * 179.97438 * 4 3 2 6 6 C 1.42895 * 114.00192 * 179.97438 * 5 4 3 7 7 C 1.50700 * 109.47075 * 180.02562 * 6 5 4 8 8 O 1.21282 * 120.00058 * 0.02562 * 7 6 5 9 9 N 1.34777 * 119.99997 * 179.97438 * 7 6 5 10 10 C 1.46925 * 120.62980 * 180.02562 * 9 7 6 11 11 C 1.53187 * 108.77688 * 126.41364 * 10 9 7 12 12 C 1.53046 * 109.31043 * 54.63722 * 11 10 9 13 13 C 1.53001 * 109.46008 * 178.63085 * 12 11 10 14 14 C 1.53000 * 109.46981 * 300.02190 * 13 12 11 15 15 C 1.52997 * 109.47440 * 60.02370 * 13 12 11 16 16 C 1.53000 * 109.47072 * 180.02562 * 13 12 11 17 17 N 1.46499 * 109.47170 * 179.97438 * 16 13 12 18 18 C 1.36932 * 120.00350 * 179.97438 * 17 16 13 19 19 H 1.08005 * 119.99868 * 359.97438 * 18 17 16 20 20 C 1.38810 * 120.00619 * 180.02562 * 18 17 16 21 21 N 1.31620 * 120.00233 * 180.02562 * 20 18 17 22 22 Si 1.86805 * 119.99577 * 359.97438 * 20 18 17 23 23 H 1.48498 * 109.47355 * 90.00244 * 22 20 18 24 24 H 1.48499 * 109.47094 * 210.00226 * 22 20 18 25 25 H 1.48503 * 109.47252 * 329.99924 * 22 20 18 26 26 C 1.53043 * 109.53475 * 298.63386 * 12 11 10 27 27 C 1.46921 * 120.63225 * 0.02846 * 9 7 6 28 28 H 1.08991 * 109.47344 * 300.00192 * 1 2 3 29 29 H 1.09001 * 109.47136 * 60.00437 * 1 2 3 30 30 H 1.09005 * 109.46987 * 179.97438 * 1 2 3 31 31 H 1.08991 * 109.47344 * 239.99808 * 2 1 3 32 32 H 1.09001 * 109.47136 * 119.99563 * 2 1 3 33 33 H 1.09003 * 109.47091 * 59.99283 * 3 2 1 34 34 H 1.09001 * 109.46784 * 299.99779 * 3 2 1 35 35 H 1.09003 * 109.47433 * 59.99738 * 4 3 2 36 36 H 1.09001 * 109.47222 * 299.99756 * 4 3 2 37 37 H 1.08997 * 109.47036 * 60.00303 * 6 5 4 38 38 H 1.08999 * 109.47491 * 300.00329 * 6 5 4 39 39 H 1.09003 * 109.59222 * 6.62634 * 10 9 7 40 40 H 1.09006 * 109.70276 * 246.27263 * 10 9 7 41 41 H 1.08996 * 109.50293 * 294.68191 * 11 10 9 42 42 H 1.09001 * 109.50064 * 174.58451 * 11 10 9 43 43 H 1.09002 * 109.46068 * 58.65418 * 12 11 10 44 44 H 1.09005 * 109.47243 * 60.00407 * 14 13 12 45 45 H 1.09000 * 109.47231 * 180.02562 * 14 13 12 46 46 H 1.09002 * 109.47199 * 300.00391 * 14 13 12 47 47 H 1.08995 * 109.47302 * 60.00023 * 15 13 12 48 48 H 1.09000 * 109.47089 * 180.02562 * 15 13 12 49 49 H 1.09005 * 109.47083 * 300.00337 * 15 13 12 50 50 H 1.09008 * 109.46946 * 299.99481 * 16 13 12 51 51 H 1.08992 * 109.47337 * 59.99460 * 16 13 12 52 52 H 0.97007 * 120.00299 * 359.97438 * 17 16 13 53 53 H 0.97004 * 119.99279 * 179.97438 * 21 20 18 54 54 H 1.09004 * 109.49679 * 181.31796 * 26 12 11 55 55 H 1.08994 * 109.49744 * 301.41312 * 26 12 11 56 56 H 1.09002 * 109.58756 * 353.36576 * 27 9 7 57 57 H 1.09002 * 109.59282 * 113.79441 * 27 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 6 1.5300 0.0000 0.0000 3 6 2.0400 1.4425 0.0000 4 6 3.5700 1.4426 0.0006 5 8 4.0463 2.7898 0.0000 6 6 5.4708 2.9027 0.0012 7 6 5.8594 4.3587 0.0010 8 8 5.0009 5.2155 0.0006 9 7 7.1609 4.7089 0.0025 10 6 7.5552 6.1242 0.0030 11 6 8.5156 6.3681 1.1712 12 6 9.6880 5.3890 1.0755 13 6 10.6669 5.6522 2.2215 14 6 11.1905 7.0866 2.1256 15 6 9.9498 5.4594 3.5592 16 6 11.8393 4.6739 2.1257 17 7 12.7769 4.9263 3.2226 18 6 13.9104 4.1675 3.3429 19 1 14.1130 3.3825 2.6293 20 6 14.7985 4.4063 4.3826 21 7 15.8883 3.6772 4.4979 22 14 14.4478 5.7635 5.6174 23 1 13.6591 5.2121 6.7483 24 1 15.7300 6.3110 6.1287 25 1 13.6756 6.8463 4.9568 26 6 9.1652 3.9538 1.1714 27 6 8.2122 3.6826 0.0036 28 1 -0.3633 0.5138 0.8899 29 1 -0.3633 0.5138 -0.8900 30 1 -0.3633 -1.0277 -0.0005 31 1 1.8934 -0.5138 -0.8899 32 1 1.8934 -0.5138 0.8900 33 1 1.6766 1.9564 0.8899 34 1 1.6767 1.9563 -0.8900 35 1 3.9338 0.9286 -0.8891 36 1 3.9330 0.9289 0.8909 37 1 5.8743 2.4186 -0.8881 38 1 5.8729 2.4198 0.8918 39 1 6.6706 6.7501 0.1205 40 1 8.0527 6.3674 -0.9360 41 1 7.9892 6.2144 2.1132 42 1 8.8910 7.3904 1.1264 43 1 10.1990 5.5264 0.1225 44 1 10.3552 7.7837 2.1940 45 1 11.8882 7.2740 2.9418 46 1 11.7013 7.2240 1.1726 47 1 9.5768 4.4375 3.6276 48 1 10.6471 5.6474 4.3756 49 1 9.1145 6.1563 3.6275 50 1 11.4663 3.6519 2.1940 51 1 12.3501 4.8112 1.1727 52 1 12.5948 5.6311 3.8639 53 1 16.5090 3.8444 5.2243 54 1 10.0029 3.2577 1.1266 55 1 8.6327 3.8227 2.1133 56 1 7.7610 2.6974 0.1217 57 1 8.7642 3.7238 -0.9354 RHF calculation, no. of doubly occupied orbitals= 67 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=WATER ZINC001101524129.mol2 57 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:14:17 Heat of formation + Delta-G solvation = -111.590759 kcal Electronic energy + Delta-G solvation = -30632.508074 eV Core-core repulsion = 26640.539972 eV Total energy + Delta-G solvation = -3991.968103 eV No. of doubly occupied orbitals = 67 Molecular weight (most abundant/longest-lived isotopes) = 340.251 amu Computer time = 1.02 seconds Orbital eigenvalues (eV) -41.49640 -40.09055 -38.79634 -37.41758 -35.81563 -35.57168 -33.41923 -32.35478 -31.88182 -30.50863 -28.50960 -28.18202 -27.72224 -26.09012 -24.40195 -24.09818 -23.62918 -22.26045 -21.93204 -20.36634 -19.86977 -19.06686 -18.30335 -17.57907 -17.47301 -16.82768 -16.56785 -16.44358 -16.07324 -15.68932 -15.53178 -15.39814 -15.23408 -15.02404 -14.73150 -14.62580 -14.31033 -14.26185 -14.02873 -13.79610 -13.57674 -13.45900 -13.19831 -13.15400 -12.90173 -12.81996 -12.76317 -12.65349 -12.43048 -12.30627 -12.16193 -11.97540 -11.94254 -11.86654 -11.82300 -11.73018 -11.53324 -11.28805 -11.17294 -10.97813 -10.85479 -10.70455 -10.68477 -9.66226 -9.21903 -8.01193 -5.31803 1.35857 1.38294 1.45663 1.52324 1.72220 2.48195 2.80583 3.30810 3.50163 3.58182 3.76040 3.81688 3.98174 3.98566 4.05855 4.12729 4.17353 4.20178 4.22888 4.24015 4.29240 4.36610 4.40366 4.42235 4.46969 4.49791 4.62655 4.68735 4.79992 4.83006 4.86814 4.90928 4.92684 4.94450 4.95525 4.99647 5.02289 5.04733 5.10234 5.17116 5.21393 5.23891 5.25052 5.34573 5.34611 5.40229 5.44251 5.50130 5.62198 5.68123 5.79307 6.09088 6.23852 6.37957 7.11567 7.11895 7.90548 8.07928 9.20638 Molecular weight = 340.25amu Principal moments of inertia in cm(-1) A = 0.017359 B = 0.001933 C = 0.001824 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1612.568644 B =14482.815390 C =15344.042048 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.154 4.154 2 C -0.117 4.117 3 C -0.119 4.119 4 C 0.056 3.944 5 O -0.400 6.400 6 C 0.046 3.954 7 C 0.509 3.491 8 O -0.578 6.578 9 N -0.604 5.604 10 C 0.113 3.887 11 C -0.123 4.123 12 C -0.064 4.064 13 C -0.060 4.060 14 C -0.139 4.139 15 C -0.138 4.138 16 C 0.163 3.837 17 N -0.755 5.755 18 C -0.227 4.227 19 H 0.058 0.942 20 C -0.064 4.064 21 N -0.984 5.984 22 Si 0.947 3.053 23 H -0.280 1.280 24 H -0.297 1.297 25 H -0.254 1.254 26 C -0.136 4.136 27 C 0.106 3.894 28 H 0.051 0.949 29 H 0.053 0.947 30 H 0.064 0.936 31 H 0.072 0.928 32 H 0.070 0.930 33 H 0.059 0.941 34 H 0.061 0.939 35 H 0.067 0.933 36 H 0.064 0.936 37 H 0.108 0.892 38 H 0.102 0.898 39 H 0.073 0.927 40 H 0.083 0.917 41 H 0.062 0.938 42 H 0.088 0.912 43 H 0.081 0.919 44 H 0.074 0.926 45 H 0.032 0.968 46 H 0.056 0.944 47 H 0.058 0.942 48 H 0.032 0.968 49 H 0.068 0.932 50 H 0.032 0.968 51 H 0.020 0.980 52 H 0.364 0.636 53 H 0.256 0.744 54 H 0.083 0.917 55 H 0.067 0.933 56 H 0.106 0.894 57 H 0.090 0.910 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -27.359 -3.311 -12.318 30.187 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.212 4.212 2 C -0.155 4.155 3 C -0.157 4.157 4 C -0.020 4.020 5 O -0.320 6.320 6 C -0.034 4.034 7 C 0.298 3.702 8 O -0.458 6.458 9 N -0.338 5.338 10 C -0.009 4.009 11 C -0.161 4.161 12 C -0.083 4.083 13 C -0.061 4.061 14 C -0.196 4.196 15 C -0.196 4.196 16 C 0.037 3.963 17 N -0.401 5.401 18 C -0.358 4.358 19 H 0.076 0.924 20 C -0.255 4.255 21 N -0.636 5.636 22 Si 0.775 3.225 23 H -0.206 1.206 24 H -0.224 1.224 25 H -0.179 1.179 26 C -0.174 4.174 27 C -0.016 4.016 28 H 0.071 0.929 29 H 0.072 0.928 30 H 0.083 0.917 31 H 0.091 0.909 32 H 0.088 0.912 33 H 0.078 0.922 34 H 0.080 0.920 35 H 0.086 0.914 36 H 0.082 0.918 37 H 0.126 0.874 38 H 0.120 0.880 39 H 0.092 0.908 40 H 0.101 0.899 41 H 0.080 0.920 42 H 0.106 0.894 43 H 0.099 0.901 44 H 0.093 0.907 45 H 0.051 0.949 46 H 0.075 0.925 47 H 0.077 0.923 48 H 0.052 0.948 49 H 0.087 0.913 50 H 0.049 0.951 51 H 0.038 0.962 52 H 0.199 0.801 53 H 0.065 0.935 54 H 0.101 0.899 55 H 0.085 0.915 56 H 0.124 0.876 57 H 0.108 0.892 Dipole moment (debyes) X Y Z Total from point charges -27.808 -2.868 -12.397 30.581 hybrid contribution 0.817 0.536 1.135 1.498 sum -26.991 -2.332 -11.262 29.339 Atomic orbital electron populations 1.21774 0.95283 1.02237 1.01857 1.21423 0.96389 0.94704 1.02966 1.21555 0.93290 0.98270 1.02577 1.22788 0.95462 0.82049 1.01707 1.87959 1.19783 1.28472 1.95754 1.22377 0.83290 0.92860 1.04827 1.20825 0.82233 0.91545 0.75616 1.90687 1.50403 1.53396 1.51291 1.48172 1.04985 1.08286 1.72341 1.21628 0.98683 0.81186 0.99421 1.22008 0.95473 1.00441 0.98130 1.21016 0.94439 0.95044 0.97812 1.20662 0.94196 0.95175 0.96028 1.21877 1.00723 0.95684 1.01362 1.21898 0.99657 1.02035 0.95966 1.20919 0.89024 0.96826 0.89494 1.42809 1.23811 1.40730 1.32786 1.19074 0.97473 1.10689 1.08558 0.92389 1.25420 0.96908 1.01940 1.01229 1.69877 1.14188 1.44475 1.35011 0.85533 0.76935 0.80355 0.79721 1.20586 1.22357 1.17876 1.22102 0.99413 0.97049 0.98821 1.22023 0.88622 0.91175 0.99750 0.92939 0.92779 0.91696 0.90916 0.91153 0.92249 0.91985 0.91445 0.91818 0.87413 0.87999 0.90841 0.89916 0.91954 0.89369 0.90080 0.90689 0.94872 0.92490 0.92253 0.94841 0.91324 0.95051 0.96238 0.80128 0.93512 0.89852 0.91471 0.87587 0.89210 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.15 -1.39 9.91 71.98 0.71 -0.67 16 2 C -0.12 -1.07 5.93 30.60 0.18 -0.88 16 3 C -0.12 -1.59 5.58 30.59 0.17 -1.41 16 4 C 0.06 0.74 5.94 71.98 0.43 1.16 16 5 O -0.40 -8.40 9.88 -129.94 -1.28 -9.69 16 6 C 0.05 0.70 4.76 71.24 0.34 1.04 16 7 C 0.51 11.41 7.78 87.66 0.68 12.09 16 8 O -0.58 -18.77 15.91 16.00 0.25 -18.51 16 9 N -0.60 -10.61 2.97 -818.99 -2.43 -13.05 16 10 C 0.11 2.02 6.43 86.36 0.56 2.57 16 11 C -0.12 -2.12 4.65 30.67 0.14 -1.98 16 12 C -0.06 -1.20 1.51 -10.76 -0.02 -1.22 16 13 C -0.06 -1.54 0.33 -52.18 -0.02 -1.56 16 14 C -0.14 -3.59 7.95 71.98 0.57 -3.02 16 15 C -0.14 -3.60 7.60 71.98 0.55 -3.05 16 16 C 0.16 5.67 4.80 86.38 0.41 6.09 16 17 N -0.76 -34.12 4.81 -343.46 -1.65 -35.77 16 18 C -0.23 -12.75 10.47 84.03 0.88 -11.87 16 19 H 0.06 3.55 8.06 -2.91 -0.02 3.53 16 20 C -0.06 -3.74 5.50 84.86 0.47 -3.27 16 21 N -0.98 -62.91 13.97 21.65 0.30 -62.61 16 22 Si 0.95 45.53 27.74 68.60 1.90 47.43 16 23 H -0.28 -13.14 7.11 99.48 0.71 -12.43 16 24 H -0.30 -14.48 7.11 99.48 0.71 -13.78 16 25 H -0.25 -11.08 6.52 99.48 0.65 -10.43 16 26 C -0.14 -2.26 4.65 30.66 0.14 -2.11 16 27 C 0.11 1.30 6.34 86.35 0.55 1.84 16 28 H 0.05 0.50 8.14 -2.39 -0.02 0.48 16 29 H 0.05 0.50 8.14 -2.39 -0.02 0.48 16 30 H 0.06 0.47 8.14 -2.38 -0.02 0.45 16 31 H 0.07 0.54 8.14 -2.39 -0.02 0.52 16 32 H 0.07 0.56 8.14 -2.39 -0.02 0.54 16 33 H 0.06 0.92 8.14 -2.39 -0.02 0.90 16 34 H 0.06 0.93 8.14 -2.39 -0.02 0.91 16 35 H 0.07 0.64 8.14 -2.39 -0.02 0.63 16 36 H 0.06 0.67 8.14 -2.39 -0.02 0.65 16 37 H 0.11 1.01 7.94 -2.39 -0.02 0.99 16 38 H 0.10 1.11 7.69 -2.39 -0.02 1.09 16 39 H 0.07 1.61 7.01 -2.39 -0.02 1.59 16 40 H 0.08 1.24 8.14 -2.38 -0.02 1.22 16 41 H 0.06 1.21 6.51 -2.39 -0.02 1.19 16 42 H 0.09 1.39 6.36 -2.39 -0.02 1.38 16 43 H 0.08 1.46 8.14 -2.39 -0.02 1.44 16 44 H 0.07 1.54 6.40 -2.38 -0.02 1.53 16 45 H 0.03 1.02 6.86 -2.39 -0.02 1.01 16 46 H 0.06 1.43 8.14 -2.39 -0.02 1.41 16 47 H 0.06 1.49 6.46 -2.39 -0.02 1.48 16 48 H 0.03 1.02 6.86 -2.39 -0.02 1.00 16 49 H 0.07 1.51 6.46 -2.38 -0.02 1.50 16 50 H 0.03 1.09 6.40 -2.38 -0.02 1.08 16 51 H 0.02 0.71 8.14 -2.39 -0.02 0.69 16 52 H 0.36 16.02 3.22 -92.70 -0.30 15.72 16 53 H 0.26 15.65 8.31 -92.71 -0.77 14.88 16 54 H 0.08 1.50 6.36 -2.39 -0.02 1.48 16 55 H 0.07 1.21 6.51 -2.39 -0.02 1.19 16 56 H 0.11 0.80 5.93 -2.39 -0.01 0.78 16 57 H 0.09 0.89 8.14 -2.39 -0.02 0.87 16 Total: -1.00 -76.80 423.47 4.31 -72.49 By element: Atomic # 1 Polarization: 25.49 SS G_CDS: 0.47 Total: 25.96 kcal Atomic # 6 Polarization: -13.00 SS G_CDS: 6.75 Total: -6.25 kcal Atomic # 7 Polarization: -107.64 SS G_CDS: -3.78 Total: -111.43 kcal Atomic # 8 Polarization: -27.17 SS G_CDS: -1.03 Total: -28.20 kcal Atomic # 14 Polarization: 45.53 SS G_CDS: 1.90 Total: 47.43 kcal Total: -76.80 4.31 -72.49 kcal The number of atoms in the molecule is 57 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. ZINC001101524129.mol2 57 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 -39.100 kcal (2) G-P(sol) polarization free energy of solvation -76.804 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -115.904 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 4.313 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -72.491 kcal (6) G-S(sol) free energy of system = (1) + (5) -111.591 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.02 seconds