Wall clock time and date at job start Wed Apr 1 2020 05:36:12
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.53000 * 1
3 3 C 1.53009 * 109.46828 * 2 1
4 4 C 1.52991 * 109.47619 * 239.99935 * 2 1 3
5 5 C 1.50696 * 109.46800 * 119.99358 * 2 1 3
6 6 C 1.39680 * 126.06736 * 75.00474 * 5 2 1
7 7 C 1.38401 * 105.35964 * 180.02562 * 6 5 2
8 8 C 1.47166 * 126.88391 * 179.97438 * 7 6 5
9 9 O 1.21629 * 119.99705 * 332.79703 * 8 7 6
10 10 N 1.34774 * 120.00046 * 152.80108 * 8 7 6
11 11 C 1.46598 * 119.69024 * 0.06208 * 10 8 7
12 12 C 1.54781 * 110.71556 * 228.03732 * 11 10 8
13 13 C 1.51520 * 106.80080 * 308.79697 * 12 11 10
14 14 H 1.09002 * 109.28321 * 306.63012 * 13 12 11
15 15 C 1.55365 * 107.78560 * 65.19025 * 13 12 11
16 16 H 1.09003 * 111.51004 * 49.83358 * 15 13 12
17 17 C 1.53734 * 100.22423 * 168.51391 * 15 13 12
18 18 C 1.55291 * 101.84216 * 320.90373 * 17 15 13
19 19 N 1.45468 * 116.86673 * 181.93021 * 13 12 11
20 20 C 1.36935 * 126.40751 * 30.89182 * 19 13 12
21 21 H 1.07997 * 119.99997 * 359.97438 * 20 19 13
22 22 C 1.38812 * 119.99963 * 179.97438 * 20 19 13
23 23 N 1.31613 * 120.00348 * 359.97438 * 22 20 19
24 24 Si 1.86802 * 119.99872 * 179.97438 * 22 20 19
25 25 H 1.48498 * 109.47325 * 89.99795 * 24 22 20
26 26 H 1.48498 * 109.46976 * 210.00073 * 24 22 20
27 27 H 1.48508 * 109.46860 * 329.99988 * 24 22 20
28 28 C 1.47722 * 119.68585 * 180.34178 * 10 8 7
29 29 N 1.36598 * 106.23411 * 359.97438 * 7 6 5
30 30 H 0.96996 * 125.19485 * 179.97438 * 29 7 6
31 31 N 1.28449 * 109.61949 * 0.22928 * 29 7 6
32 32 H 1.08998 * 109.47388 * 300.00638 * 1 2 3
33 33 H 1.09004 * 109.46948 * 60.00537 * 1 2 3
34 34 H 1.08996 * 109.46928 * 179.97438 * 1 2 3
35 35 H 1.09005 * 109.46900 * 179.97438 * 3 2 1
36 36 H 1.08997 * 109.47267 * 299.99360 * 3 2 1
37 37 H 1.08997 * 109.47170 * 59.99722 * 3 2 1
38 38 H 1.08996 * 109.47075 * 300.00038 * 4 2 1
39 39 H 1.09004 * 109.47346 * 59.99771 * 4 2 1
40 40 H 1.08998 * 109.47530 * 180.02562 * 4 2 1
41 41 H 1.07997 * 127.31610 * 359.97438 * 6 5 2
42 42 H 1.09004 * 109.25378 * 348.13880 * 11 10 8
43 43 H 1.08999 * 109.25344 * 107.56820 * 11 10 8
44 44 H 1.08996 * 109.87145 * 189.66555 * 12 11 10
45 45 H 1.09005 * 110.49486 * 68.27263 * 12 11 10
46 46 H 1.09002 * 110.95738 * 79.03729 * 17 15 13
47 47 H 1.08998 * 110.96030 * 202.76592 * 17 15 13
48 48 H 1.08997 * 109.91884 * 263.13835 * 18 17 15
49 49 H 1.08997 * 110.09049 * 141.91964 * 18 17 15
50 50 H 0.96998 * 120.00329 * 179.97438 * 23 22 20
51 51 H 1.08993 * 110.07351 * 247.12497 * 28 10 8
52 52 H 1.09003 * 110.06543 * 8.59609 * 28 10 8
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.4426 0.0000
4 6 2.0401 -0.7212 -1.2491
5 6 2.0322 -0.7103 1.2305
6 6 2.0309 -2.0911 1.4412
7 6 2.5812 -2.2707 2.6983
8 6 2.8022 -3.5501 3.3912
9 8 2.0947 -4.5073 3.1409
10 7 3.7884 -3.6611 4.3030
11 6 4.6361 -2.5026 4.6004
12 6 6.1322 -2.8987 4.5779
13 6 6.2770 -4.1138 5.4714
14 1 5.8496 -3.8941 6.4498
15 6 5.4854 -5.2803 4.8183
16 1 5.7439 -5.3988 3.7660
17 6 5.9843 -6.4590 5.6698
18 6 7.4864 -6.1086 5.8502
19 7 7.6194 -4.6476 5.6427
20 6 8.7874 -3.9334 5.6167
21 1 8.7641 -2.8657 5.4561
22 6 10.0015 -4.5818 5.7968
23 7 10.0300 -5.8829 5.9930
24 14 11.5950 -3.6076 5.7606
25 1 11.9225 -3.1432 7.1325
26 1 12.6920 -4.4717 5.2554
27 1 11.4372 -2.4336 4.8650
28 6 4.0089 -4.9507 4.9889
29 7 2.8781 -1.0225 3.1671
30 1 3.2774 -0.8297 4.0297
31 7 2.5430 -0.1427 2.2933
32 1 -0.3634 0.5139 0.8899
33 1 -0.3633 0.5138 -0.8901
34 1 -0.3633 -1.0276 -0.0005
35 1 3.1300 1.4426 0.0005
36 1 1.6767 1.9564 -0.8900
37 1 1.6766 1.9564 0.8899
38 1 1.6768 -1.7488 -1.2491
39 1 1.6768 -0.2074 -2.1392
40 1 3.1301 -0.7216 -1.2489
41 1 1.6752 -2.8559 0.7667
42 1 4.4619 -1.7284 3.8531
43 1 4.3805 -2.1138 5.5861
44 1 6.7382 -2.0796 4.9650
45 1 6.4465 -3.1449 3.5636
46 1 5.4728 -6.4890 6.6319
47 1 5.8655 -7.4026 5.1372
48 1 8.0839 -6.6442 5.1126
49 1 7.8142 -6.3722 6.8558
50 1 10.8784 -6.3361 6.1185
51 1 3.7663 -4.8578 6.0474
52 1 3.3955 -5.7271 4.5317
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
ZINC001036756224.mol2 52
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
SCF FIELD WAS ACHIEVED
Wall clock time and date at job start Wed Apr 1 2020 05:36:12
Heat of formation + Delta-G solvation = 74.207772 kcal
Electronic energy + Delta-G solvation = -31469.514693 eV
Core-core repulsion = 27471.604466 eV
Total energy + Delta-G solvation = -3997.910227 eV
No. of doubly occupied orbitals = 66
Molecular weight
(most abundant/longest-lived isotopes) = 346.211 amu
Computer time = 1.43 seconds
Orbital eigenvalues (eV)
-42.90032 -40.16427 -38.29158 -37.18887 -35.59083 -32.90378 -31.72986 -31.24355
-31.17913 -29.53476 -28.16743 -27.47961 -27.44519 -25.31261 -24.19111 -23.17919
-22.41799 -22.14155 -20.68928 -20.12454 -19.61019 -18.93593 -17.75744 -16.70489
-16.31694 -16.01804 -15.81918 -15.20086 -15.02549 -14.86032 -14.45716 -14.18926
-14.13165 -13.92165 -13.88252 -13.51509 -13.30338 -13.12196 -12.70929 -12.59393
-12.48244 -12.36332 -11.86879 -11.76064 -11.60273 -11.55196 -11.49822 -11.30535
-11.23010 -10.94999 -10.86548 -10.84466 -10.71902 -10.29108 -10.17497 -10.03022
-9.94261 -9.71175 -9.59163 -9.27545 -9.26067 -8.71798 -8.68716 -6.91339
-5.78084 -3.11826 1.17570 2.36673 2.55518 2.72358 3.12811 3.34959
3.39973 3.42139 4.07560 4.29716 4.42606 4.53508 4.65370 4.77685
4.88019 4.95633 5.03234 5.10492 5.20907 5.22866 5.33253 5.39975
5.43269 5.46397 5.48301 5.52885 5.57813 5.67733 5.68652 5.69893
5.76192 5.83408 5.88989 5.91128 6.00230 6.04096 6.09029 6.11023
6.12995 6.20966 6.30666 6.34757 6.50893 6.54147 6.68403 6.75791
6.95645 7.20868 7.48312 7.54702 7.74856 7.79581 7.96936 8.17513
9.12367 9.75934 10.43410 11.34322
Molecular weight = 346.21amu
Principal moments of inertia in cm(-1)
A = 0.014932 B = 0.002642 C = 0.002366
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 1874.656372 B =10593.659758 C =11833.363037
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C -0.137 4.137
2 C 0.017 3.983
3 C -0.135 4.135
4 C -0.140 4.140
5 C 0.041 3.959
6 C -0.153 4.153
7 C -0.027 4.027
8 C 0.586 3.414
9 O -0.526 6.526
10 N -0.600 5.600
11 C 0.089 3.911
12 C -0.126 4.126
13 C 0.167 3.833
14 H 0.021 0.979
15 C -0.125 4.125
16 H 0.081 0.919
17 C -0.134 4.134
18 C 0.176 3.824
19 N -0.594 5.594
20 C -0.229 4.229
21 H 0.078 0.922
22 C -0.006 4.006
23 N -0.920 5.920
24 Si 0.910 3.090
25 H -0.289 1.289
26 H -0.290 1.290
27 H -0.281 1.281
28 C 0.127 3.873
29 N -0.383 5.383
30 H 0.448 0.552
31 N -0.271 5.271
32 H 0.056 0.944
33 H 0.062 0.938
34 H 0.059 0.941
35 H 0.057 0.943
36 H 0.056 0.944
37 H 0.061 0.939
38 H 0.060 0.940
39 H 0.062 0.938
40 H 0.058 0.942
41 H 0.161 0.839
42 H 0.078 0.922
43 H 0.073 0.927
44 H 0.079 0.921
45 H 0.071 0.929
46 H 0.059 0.941
47 H 0.062 0.938
48 H 0.055 0.945
49 H 0.040 0.960
50 H 0.253 0.747
51 H 0.073 0.927
52 H 0.094 0.906
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges -17.864 9.996 -8.731 22.255
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.195 4.195
2 C 0.015 3.985
3 C -0.193 4.193
4 C -0.198 4.198
5 C -0.118 4.118
6 C -0.182 4.182
7 C -0.155 4.155
8 C 0.372 3.628
9 O -0.402 6.402
10 N -0.333 5.333
11 C -0.035 4.035
12 C -0.165 4.165
13 C 0.058 3.942
14 H 0.039 0.961
15 C -0.145 4.145
16 H 0.099 0.901
17 C -0.173 4.173
18 C 0.052 3.948
19 N -0.319 5.319
20 C -0.357 4.357
21 H 0.096 0.904
22 C -0.203 4.203
23 N -0.567 5.567
24 Si 0.740 3.260
25 H -0.216 1.216
26 H -0.217 1.217
27 H -0.207 1.207
28 C 0.006 3.994
29 N -0.104 5.104
30 H 0.298 0.702
31 N -0.103 5.103
32 H 0.075 0.925
33 H 0.081 0.919
34 H 0.078 0.922
35 H 0.076 0.924
36 H 0.075 0.925
37 H 0.080 0.920
38 H 0.079 0.921
39 H 0.081 0.919
40 H 0.077 0.923
41 H 0.178 0.822
42 H 0.096 0.904
43 H 0.091 0.909
44 H 0.098 0.902
45 H 0.090 0.910
46 H 0.078 0.922
47 H 0.081 0.919
48 H 0.073 0.927
49 H 0.058 0.942
50 H 0.063 0.937
51 H 0.091 0.909
52 H 0.112 0.888
Dipole moment (debyes) X Y Z Total
from point charges -17.693 10.442 -8.329 22.168
hybrid contribution 1.139 -0.184 0.278 1.187
sum -16.554 10.258 -8.051 21.073
Atomic orbital electron populations
1.21870 0.93346 1.02341 1.01992 1.18906 0.95360 0.94844 0.89385
1.21774 1.01408 0.94247 1.01836 1.21887 1.01526 1.00476 0.95890
1.21889 1.02940 0.91532 0.95435 1.21938 1.04272 0.96565 0.95458
1.21713 1.09022 0.86747 0.98054 1.17126 0.80806 0.85960 0.78880
1.90820 1.47061 1.38223 1.64051 1.48251 1.31442 1.13390 1.40181
1.22402 0.92768 0.87540 1.00833 1.22023 0.92121 1.00425 1.01971
1.20665 0.87745 0.91270 0.94513 0.96146 1.22655 0.92976 0.98247
1.00635 0.90061 1.22697 0.99703 0.96584 0.98342 1.21549 0.88737
0.84053 1.00485 1.44205 0.99795 1.03682 1.84187 1.18967 0.82809
0.91095 1.42796 0.90422 1.24795 0.97649 0.95882 1.01979 1.69930
1.24615 1.07013 1.55163 0.86157 0.83423 0.78737 0.77656 1.21593
1.21653 1.20708 1.21719 0.97255 0.83764 0.96642 1.44503 1.45999
1.02881 1.16997 0.70171 1.73074 1.16761 1.20122 1.00322 0.92509
0.91938 0.92189 0.92405 0.92472 0.92034 0.92082 0.91886 0.92277
0.82154 0.90385 0.90900 0.90238 0.91001 0.92227 0.91904 0.92686
0.94229 0.93677 0.90886 0.88809
Number of geometries 1
Number of calculations of the screened coulomb radii 1
The total number of SCF iterations 15.
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 -1.03 8.46 37.16 0.31 -0.71 16
2 C 0.02 0.15 0.96 -155.66 -0.15 0.00 16
3 C -0.14 -1.10 8.04 37.16 0.30 -0.81 16
4 C -0.14 -1.06 8.20 37.15 0.30 -0.76 16
5 C 0.04 0.45 5.65 -82.47 -0.47 -0.01 16
6 C -0.15 -1.85 10.22 -38.98 -0.40 -2.25 16
7 C -0.03 -0.35 6.43 -82.06 -0.53 -0.88 16
8 C 0.59 8.54 7.63 -12.64 -0.10 8.44 16
9 O -0.53 -9.11 16.37 5.25 0.09 -9.02 16
10 N -0.60 -8.19 3.00 -173.33 -0.52 -8.71 16
11 C 0.09 1.07 5.27 -3.05 -0.02 1.05 16
12 C -0.13 -2.01 5.52 -26.54 -0.15 -2.15 16
13 C 0.17 3.13 2.57 -67.98 -0.17 2.96 16
14 H 0.02 0.40 7.92 -51.93 -0.41 -0.01 16
15 C -0.12 -2.11 2.86 -89.37 -0.26 -2.37 16
16 H 0.08 1.51 8.11 -51.93 -0.42 1.09 16
17 C -0.13 -2.38 6.54 -25.14 -0.16 -2.55 16
18 C 0.18 4.10 5.45 -1.91 -0.01 4.09 16
19 N -0.59 -14.28 3.35 -163.09 -0.55 -14.83 16
20 C -0.23 -6.04 9.86 -15.06 -0.15 -6.19 16
21 H 0.08 1.96 7.18 -52.49 -0.38 1.59 16
22 C -0.01 -0.18 5.49 -17.43 -0.10 -0.27 16
23 N -0.92 -26.42 10.43 55.49 0.58 -25.84 16
24 Si 0.91 21.57 29.55 -169.99 -5.02 16.54 16
25 H -0.29 -6.84 7.11 56.52 0.40 -6.44 16
26 H -0.29 -7.00 7.11 56.52 0.40 -6.59 16
27 H -0.28 -6.49 7.11 56.52 0.40 -6.08 16
28 C 0.13 1.82 5.96 -3.76 -0.02 1.80 16
29 N -0.38 -4.23 5.61 33.48 0.19 -4.05 16
30 H 0.45 3.94 6.35 -190.46 -1.21 2.73 16
31 N -0.27 -3.19 11.48 33.64 0.39 -2.80 16
32 H 0.06 0.45 8.14 -51.93 -0.42 0.03 16
33 H 0.06 0.37 8.14 -51.93 -0.42 -0.05 16
34 H 0.06 0.45 8.14 -51.93 -0.42 0.03 16
35 H 0.06 0.52 8.14 -51.93 -0.42 0.09 16
36 H 0.06 0.37 8.14 -51.93 -0.42 -0.05 16
37 H 0.06 0.54 8.08 -51.93 -0.42 0.12 16
38 H 0.06 0.45 7.60 -51.93 -0.39 0.06 16
39 H 0.06 0.38 8.14 -51.93 -0.42 -0.04 16
40 H 0.06 0.49 8.14 -51.93 -0.42 0.07 16
41 H 0.16 1.77 7.79 -52.49 -0.41 1.36 16
42 H 0.08 0.83 4.02 -51.93 -0.21 0.63 16
43 H 0.07 0.74 8.14 -51.93 -0.42 0.32 16
44 H 0.08 1.27 7.59 -51.93 -0.39 0.87 16
45 H 0.07 1.25 8.11 -51.93 -0.42 0.83 16
46 H 0.06 0.97 8.12 -51.93 -0.42 0.54 16
47 H 0.06 1.02 8.14 -51.93 -0.42 0.60 16
48 H 0.05 1.44 7.21 -51.93 -0.37 1.06 16
49 H 0.04 1.01 7.76 -51.93 -0.40 0.61 16
50 H 0.25 7.09 8.31 -40.82 -0.34 6.75 16
51 H 0.07 0.96 8.08 -51.93 -0.42 0.54 16
52 H 0.09 1.38 6.99 -51.93 -0.36 1.01 16
LS Contribution 398.71 15.07 6.01 6.01
Total: -1.00 -31.45 398.71 -10.18 -41.63
By element:
Atomic # 1 Polarization: 11.23 SS G_CDS: -9.59 Total: 1.65 kcal
Atomic # 6 Polarization: 1.16 SS G_CDS: -1.75 Total: -0.59 kcal
Atomic # 7 Polarization: -56.30 SS G_CDS: 0.09 Total: -56.22 kcal
Atomic # 8 Polarization: -9.11 SS G_CDS: 0.09 Total: -9.02 kcal
Atomic # 14 Polarization: 21.57 SS G_CDS: -5.02 Total: 16.54 kcal
Total LS contribution 6.01 Total: 6.01 kcal
Total: -31.45 -10.18 -41.63 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.
ZINC001036756224.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 115.840 kcal
(2) G-P(sol) polarization free energy of solvation -31.448 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system 84.392 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.184 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.632 kcal
(6) G-S(sol) free energy of system = (1) + (5) 74.208 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 1.44 seconds