Wall clock time and date at job start Sat Apr 11 2020 06:10:00 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.52996 * 1 3 3 H 1.09002 * 115.54951 * 2 1 4 4 C 1.53000 * 117.50107 * 214.32788 * 2 1 3 5 5 C 1.53002 * 117.49888 * 145.67026 * 2 1 3 6 6 H 1.08999 * 117.49664 * 0.02562 * 5 2 1 7 7 C 1.52999 * 117.49968 * 215.00840 * 5 2 1 8 8 N 1.46494 * 109.47048 * 174.99938 * 7 5 2 9 9 C 1.34783 * 119.99885 * 179.97438 * 8 7 5 10 10 O 1.21279 * 120.00134 * 359.97438 * 9 8 7 11 11 C 1.50700 * 120.00020 * 179.97438 * 9 8 7 12 12 H 1.09004 * 109.46838 * 60.00084 * 11 9 8 13 13 C 1.53004 * 109.47457 * 300.00473 * 11 9 8 14 14 C 1.52999 * 109.47067 * 174.99698 * 13 11 9 15 15 C 1.50694 * 109.47127 * 179.97438 * 14 13 11 16 16 H 1.08005 * 119.99900 * 0.02562 * 15 14 13 17 17 C 1.37876 * 120.00240 * 180.02562 * 15 14 13 18 18 N 1.31516 * 120.00119 * 0.02562 * 17 15 14 19 19 Si 1.86800 * 120.00030 * 180.02562 * 17 15 14 20 20 H 1.48505 * 109.47451 * 59.99934 * 19 17 15 21 21 H 1.48499 * 109.47208 * 179.97438 * 19 17 15 22 22 H 1.48495 * 109.47302 * 300.00097 * 19 17 15 23 23 N 1.46492 * 109.47234 * 179.97438 * 11 9 8 24 24 C 1.34777 * 120.00678 * 205.00035 * 23 11 9 25 25 O 1.21512 * 119.99660 * 0.02562 * 24 23 11 26 26 O 1.34634 * 120.00593 * 180.02562 * 24 23 11 27 27 C 1.45202 * 117.00374 * 180.02562 * 26 24 23 28 28 C 1.52999 * 109.47039 * 59.99750 * 27 26 24 29 29 C 1.52994 * 109.47210 * 180.02562 * 27 26 24 30 30 C 1.53002 * 109.46845 * 300.00345 * 27 26 24 31 31 H 1.09003 * 109.47442 * 154.41657 * 1 2 3 32 32 H 1.09000 * 109.47181 * 274.41811 * 1 2 3 33 33 H 1.09005 * 109.47217 * 34.42156 * 1 2 3 34 34 H 1.09008 * 117.49673 * 359.97438 * 4 2 1 35 35 H 1.08999 * 117.49942 * 144.97366 * 4 2 1 36 36 H 1.09002 * 109.46825 * 294.99320 * 7 5 2 37 37 H 1.08997 * 109.47412 * 54.99052 * 7 5 2 38 38 H 0.97004 * 120.00338 * 359.97438 * 8 7 5 39 39 H 1.09001 * 109.46496 * 294.99782 * 13 11 9 40 40 H 1.09001 * 109.46772 * 54.99861 * 13 11 9 41 41 H 1.08999 * 109.47564 * 300.00229 * 14 13 11 42 42 H 1.09009 * 109.46772 * 59.99682 * 14 13 11 43 43 H 0.96998 * 120.00050 * 179.97438 * 18 17 15 44 44 H 0.97000 * 120.00140 * 25.00468 * 23 11 9 45 45 H 1.08997 * 109.47124 * 60.00040 * 28 27 26 46 46 H 1.09000 * 109.47026 * 180.02562 * 28 27 26 47 47 H 1.09004 * 109.47558 * 300.00298 * 28 27 26 48 48 H 1.08995 * 109.47612 * 59.99791 * 29 27 26 49 49 H 1.09004 * 109.47206 * 179.97438 * 29 27 26 50 50 H 1.08999 * 109.47189 * 299.99470 * 29 27 26 51 51 H 1.08992 * 109.47407 * 59.99537 * 30 27 26 52 52 H 1.08995 * 109.47117 * 179.97438 * 30 27 26 53 53 H 1.09008 * 109.46759 * 299.99606 * 30 27 26 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 1 2.0001 0.9834 0.0000 4 6 2.2365 -1.1207 -0.7653 5 6 2.2364 -1.1207 0.7654 6 1 1.6112 -1.8583 1.2685 7 6 3.5486 -0.7753 1.4723 8 7 4.0081 -1.9337 2.2423 9 6 5.1521 -1.8609 2.9513 10 8 5.8013 -0.8366 2.9513 11 6 5.6245 -3.0524 3.7439 12 1 5.7843 -3.8956 3.0718 13 6 4.5671 -3.4239 4.7855 14 6 4.9888 -4.7064 5.5054 15 6 3.9477 -5.0718 6.5318 16 1 3.0725 -4.4516 6.6577 17 6 4.1124 -6.1980 7.3100 18 7 5.1783 -6.9529 7.1571 19 14 2.8215 -6.6514 8.5818 20 1 1.5168 -6.8694 7.9069 21 1 3.2308 -7.8961 9.2806 22 1 2.6902 -5.5494 9.5685 23 7 6.8815 -2.7221 4.4197 24 6 7.7435 -3.7020 4.7562 25 8 7.4781 -4.8595 4.4987 26 8 8.8986 -3.3985 5.3778 27 6 9.7790 -4.5044 5.7098 28 6 10.1818 -5.2379 4.4289 29 6 11.0312 -3.9680 6.4063 30 6 9.0512 -5.4727 6.6445 31 1 -0.3634 -0.9269 -0.4438 32 1 -0.3633 0.0792 1.0246 33 1 -0.3634 0.8478 -0.5809 34 1 1.6112 -1.8584 -1.2685 35 1 3.1711 -0.8746 -1.2692 36 1 4.3020 -0.5102 0.7305 37 1 3.3883 0.0676 2.1445 38 1 3.4890 -2.7532 2.2421 39 1 4.4712 -2.6152 5.5101 40 1 3.6094 -3.5834 4.2901 41 1 5.0848 -5.5150 4.7809 42 1 5.9466 -4.5468 6.0008 43 1 5.2944 -7.7450 7.7049 44 1 7.0936 -1.7981 4.6250 45 1 10.7002 -4.5481 3.7629 46 1 10.8423 -6.0683 4.6782 47 1 9.2895 -5.6201 3.9331 48 1 10.7439 -3.4458 7.3188 49 1 11.6919 -4.7983 6.6559 50 1 11.5498 -3.2781 5.7407 51 1 8.1589 -5.8547 6.1488 52 1 9.7119 -6.3031 6.8934 53 1 8.7643 -4.9501 7.5571 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 ZINC000496961654.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:10:00 Heat of formation + Delta-G solvation = -110.878474 kcal Electronic energy + Delta-G solvation = -30999.247701 eV Core-core repulsion = 26871.477391 eV Total energy + Delta-G solvation = -4127.770310 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.63 seconds Orbital eigenvalues (eV) -40.49949 -39.70030 -38.28600 -36.99428 -34.71450 -34.46697 -33.73663 -32.88551 -30.93088 -29.46354 -27.84211 -27.14117 -26.86333 -26.01369 -23.96741 -22.05103 -21.41942 -21.05923 -20.49382 -19.78037 -19.58067 -18.04314 -16.93508 -16.41700 -16.12404 -15.96489 -15.54051 -15.34239 -15.06661 -14.67341 -14.40585 -14.16419 -14.03673 -13.93447 -13.47754 -13.38263 -13.13433 -13.08842 -12.65919 -12.47226 -12.23524 -12.07293 -11.95413 -11.73356 -11.69412 -11.59681 -11.57811 -11.36730 -11.06510 -10.94012 -10.85470 -10.83793 -10.65862 -10.38467 -10.28484 -9.93427 -9.92439 -9.76256 -9.47761 -9.21567 -9.09477 -8.99850 -8.45135 -8.23256 -5.88177 -3.91860 2.48169 2.68542 3.15652 3.36090 3.41049 3.42918 3.47254 3.57456 4.31892 4.50842 4.52374 4.59235 4.61323 4.83446 4.96572 5.05096 5.07771 5.10274 5.29014 5.30392 5.45281 5.48315 5.52000 5.53025 5.57202 5.58431 5.67602 5.69956 5.73773 5.82846 5.89244 5.94354 6.00907 6.01268 6.08987 6.22516 6.29294 6.39243 6.39879 6.56460 6.61049 6.79695 7.00696 7.06853 7.17453 7.47022 7.82354 7.95830 7.98823 8.16074 8.21471 8.50337 9.00557 9.16687 9.70426 11.20874 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.007866 B = 0.003429 C = 0.002504 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 3558.803860 B = 8163.633660 C =11178.190284 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.124 4.124 2 C -0.145 4.145 3 H 0.106 0.894 4 C -0.172 4.172 5 C -0.167 4.167 6 H 0.114 0.886 7 C 0.153 3.847 8 N -0.723 5.723 9 C 0.517 3.483 10 O -0.544 6.544 11 C 0.148 3.852 12 H 0.114 0.886 13 C -0.120 4.120 14 C 0.033 3.967 15 C -0.527 4.527 16 H 0.050 0.950 17 C 0.061 3.939 18 N -0.882 5.882 19 Si 0.885 3.115 20 H -0.278 1.278 21 H -0.288 1.288 22 H -0.279 1.279 23 N -0.695 5.695 24 C 0.652 3.348 25 O -0.538 6.538 26 O -0.370 6.370 27 C 0.138 3.862 28 C -0.181 4.181 29 C -0.138 4.138 30 C -0.183 4.183 31 H 0.057 0.943 32 H 0.061 0.939 33 H 0.060 0.940 34 H 0.096 0.904 35 H 0.096 0.904 36 H 0.068 0.932 37 H 0.071 0.929 38 H 0.407 0.593 39 H 0.056 0.944 40 H 0.062 0.938 41 H 0.024 0.976 42 H 0.027 0.973 43 H 0.258 0.742 44 H 0.407 0.593 45 H 0.053 0.947 46 H 0.066 0.934 47 H 0.090 0.910 48 H 0.065 0.935 49 H 0.076 0.924 50 H 0.061 0.939 51 H 0.108 0.892 52 H 0.060 0.940 53 H 0.057 0.943 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 1.213 10.363 -12.333 16.155 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.181 4.181 2 C -0.164 4.164 3 H 0.124 0.876 4 C -0.209 4.209 5 C -0.185 4.185 6 H 0.132 0.868 7 C 0.029 3.971 8 N -0.376 5.376 9 C 0.303 3.697 10 O -0.422 6.422 11 C 0.041 3.959 12 H 0.131 0.869 13 C -0.159 4.159 14 C -0.004 4.004 15 C -0.565 4.565 16 H 0.068 0.932 17 C -0.140 4.140 18 N -0.520 5.520 19 Si 0.713 3.287 20 H -0.204 1.204 21 H -0.214 1.214 22 H -0.205 1.205 23 N -0.356 5.356 24 C 0.405 3.595 25 O -0.424 6.424 26 O -0.285 6.285 27 C 0.101 3.899 28 C -0.238 4.238 29 C -0.196 4.196 30 C -0.240 4.240 31 H 0.076 0.924 32 H 0.080 0.920 33 H 0.079 0.921 34 H 0.115 0.885 35 H 0.114 0.886 36 H 0.086 0.914 37 H 0.089 0.911 38 H 0.242 0.758 39 H 0.075 0.925 40 H 0.080 0.920 41 H 0.043 0.957 42 H 0.046 0.954 43 H 0.068 0.932 44 H 0.243 0.757 45 H 0.072 0.928 46 H 0.084 0.916 47 H 0.109 0.891 48 H 0.084 0.916 49 H 0.095 0.905 50 H 0.080 0.920 51 H 0.126 0.874 52 H 0.079 0.921 53 H 0.076 0.924 Dipole moment (debyes) X Y Z Total from point charges 1.721 10.312 -12.372 16.197 hybrid contribution -0.905 -0.349 0.932 1.345 sum 0.815 9.962 -11.440 15.192 Atomic orbital electron populations 1.21379 0.92575 1.02022 1.02134 1.22377 0.95438 0.94607 1.03957 0.87611 1.23275 1.02733 1.01349 0.93568 1.22606 0.98550 1.03019 0.94333 0.86837 1.20845 0.93244 0.88552 0.94448 1.45876 1.25048 1.17005 1.49699 1.20857 0.80407 0.89658 0.78793 1.90663 1.55727 1.32036 1.63735 1.20929 0.84955 0.94760 0.95273 0.86865 1.22194 0.95074 1.02308 0.96352 1.19132 0.96777 0.89698 0.94800 1.21217 1.06498 1.09204 1.19613 0.93223 1.24964 0.95382 0.96530 0.97168 1.69999 1.26897 1.15332 1.39799 0.86793 0.82224 0.77695 0.82002 1.20424 1.21397 1.20485 1.44139 1.18705 1.09557 1.63210 1.17799 0.79562 0.84513 0.77626 1.91013 1.76296 1.17657 1.57417 1.86227 1.29472 1.40354 1.72489 1.22011 0.87934 0.84667 0.95268 1.22496 1.03643 1.00377 0.97281 1.21837 0.95191 1.01778 1.00752 1.22515 1.03391 0.98474 0.99578 0.92430 0.91959 0.92090 0.88550 0.88605 0.91410 0.91053 0.75755 0.92489 0.91951 0.95721 0.95446 0.93237 0.75668 0.92819 0.91552 0.89103 0.91632 0.90507 0.91996 0.87371 0.92082 0.92423 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.12 -0.81 10.19 37.16 0.38 -0.43 16 2 C -0.15 -1.07 6.12 -90.62 -0.56 -1.63 16 3 H 0.11 0.75 8.14 -51.93 -0.42 0.33 16 4 C -0.17 -1.28 9.83 -26.70 -0.26 -1.54 16 5 C -0.17 -1.39 5.57 -90.58 -0.50 -1.89 16 6 H 0.11 0.93 8.04 -51.93 -0.42 0.51 16 7 C 0.15 1.59 5.96 -4.05 -0.02 1.57 16 8 N -0.72 -8.68 5.42 -60.30 -0.33 -9.01 16 9 C 0.52 7.92 7.22 -10.98 -0.08 7.84 16 10 O -0.54 -9.59 15.84 5.56 0.09 -9.50 16 11 C 0.15 2.55 2.66 -69.08 -0.18 2.36 16 12 H 0.11 1.97 7.58 -51.93 -0.39 1.58 16 13 C -0.12 -2.36 4.77 -26.73 -0.13 -2.49 16 14 C 0.03 0.84 3.57 -27.89 -0.10 0.74 16 15 C -0.53 -14.80 9.11 -34.44 -0.31 -15.11 16 16 H 0.05 1.35 7.74 -52.48 -0.41 0.94 16 17 C 0.06 1.76 5.46 -17.82 -0.10 1.66 16 18 N -0.88 -26.38 13.29 55.43 0.74 -25.64 16 19 Si 0.88 21.71 29.54 -169.99 -5.02 16.69 16 20 H -0.28 -6.69 7.11 56.52 0.40 -6.29 16 21 H -0.29 -7.03 7.11 56.52 0.40 -6.63 16 22 H -0.28 -6.78 7.11 56.52 0.40 -6.38 16 23 N -0.70 -12.58 4.91 -59.87 -0.29 -12.87 16 24 C 0.65 12.95 6.62 54.06 0.36 13.31 16 25 O -0.54 -11.96 9.02 -19.28 -0.17 -12.14 16 26 O -0.37 -6.76 9.94 -40.31 -0.40 -7.16 16 27 C 0.14 2.23 1.13 -90.62 -0.10 2.13 16 28 C -0.18 -2.71 8.37 37.16 0.31 -2.40 16 29 C -0.14 -1.72 8.85 37.15 0.33 -1.39 16 30 C -0.18 -3.37 8.37 37.16 0.31 -3.06 16 31 H 0.06 0.36 7.90 -51.93 -0.41 -0.05 16 32 H 0.06 0.42 8.14 -51.93 -0.42 -0.01 16 33 H 0.06 0.35 8.14 -51.93 -0.42 -0.07 16 34 H 0.10 0.66 8.05 -51.92 -0.42 0.24 16 35 H 0.10 0.70 8.04 -51.93 -0.42 0.28 16 36 H 0.07 0.73 7.88 -51.93 -0.41 0.32 16 37 H 0.07 0.80 8.14 -51.93 -0.42 0.37 16 38 H 0.41 4.26 7.76 -40.82 -0.32 3.94 16 39 H 0.06 1.15 8.14 -51.93 -0.42 0.73 16 40 H 0.06 1.10 7.43 -51.93 -0.39 0.71 16 41 H 0.02 0.68 7.84 -51.93 -0.41 0.28 16 42 H 0.03 0.76 5.65 -51.92 -0.29 0.47 16 43 H 0.26 7.37 8.31 -40.82 -0.34 7.03 16 44 H 0.41 6.93 7.90 -40.82 -0.32 6.61 16 45 H 0.05 0.73 8.14 -51.93 -0.42 0.30 16 46 H 0.07 0.85 8.14 -51.93 -0.42 0.43 16 47 H 0.09 1.61 5.88 -51.93 -0.31 1.30 16 48 H 0.06 0.81 8.14 -51.93 -0.42 0.39 16 49 H 0.08 0.80 8.14 -51.93 -0.42 0.37 16 50 H 0.06 0.70 8.14 -51.93 -0.42 0.28 16 51 H 0.11 2.49 5.88 -51.93 -0.31 2.19 16 52 H 0.06 0.97 8.14 -51.93 -0.42 0.55 16 53 H 0.06 1.09 8.14 -51.92 -0.42 0.67 16 LS Contribution 422.65 15.07 6.37 6.37 Total: -1.00 -33.10 422.65 -9.10 -42.20 By element: Atomic # 1 Polarization: 20.80 SS G_CDS: -9.42 Total: 11.39 kcal Atomic # 6 Polarization: 0.33 SS G_CDS: -0.66 Total: -0.33 kcal Atomic # 7 Polarization: -47.64 SS G_CDS: 0.12 Total: -47.52 kcal Atomic # 8 Polarization: -28.31 SS G_CDS: -0.49 Total: -28.80 kcal Atomic # 14 Polarization: 21.71 SS G_CDS: -5.02 Total: 16.69 kcal Total LS contribution 6.37 Total: 6.37 kcal Total: -33.10 -9.10 -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. ZINC000496961654.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 -68.678 kcal (2) G-P(sol) polarization free energy of solvation -33.100 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -101.778 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.100 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.200 kcal (6) G-S(sol) free energy of system = (1) + (5) -110.878 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.63 seconds