Wall clock time and date at job start Sat Apr 11 2020 06:26:58
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.53008 * 1
3 3 H 1.09001 * 109.46796 * 2 1
4 4 C 1.53001 * 109.46508 * 119.99286 * 2 1 3
5 5 N 1.46902 * 109.46534 * 175.00072 * 4 2 1
6 6 C 1.46907 * 110.99776 * 180.02562 * 5 4 2
7 7 C 1.50699 * 109.46914 * 179.97438 * 6 5 4
8 8 O 1.21291 * 120.00563 * 0.02562 * 7 6 5
9 9 N 1.34778 * 119.99610 * 179.97438 * 7 6 5
10 10 C 1.37107 * 119.99906 * 179.97438 * 9 7 6
11 11 H 1.08004 * 119.99254 * 0.02562 * 10 9 7
12 12 C 1.38948 * 120.00293 * 180.02562 * 10 9 7
13 13 N 1.31408 * 119.99998 * 180.02562 * 12 10 9
14 14 Si 1.86800 * 119.99599 * 0.02562 * 12 10 9
15 15 H 1.48500 * 109.47280 * 90.00005 * 14 12 10
16 16 H 1.48505 * 109.47000 * 210.00048 * 14 12 10
17 17 H 1.48496 * 109.47304 * 329.99765 * 14 12 10
18 18 N 1.46497 * 109.47335 * 239.99136 * 2 1 3
19 19 C 1.34777 * 119.99525 * 84.99977 * 18 2 1
20 20 O 1.21283 * 120.00463 * 0.02562 * 19 18 2
21 21 C 1.50700 * 119.99675 * 179.72583 * 19 18 2
22 22 O 1.42899 * 109.47143 * 180.27800 * 21 19 18
23 23 C 1.42899 * 114.00265 * 179.97438 * 22 21 19
24 24 C 1.52993 * 109.47343 * 270.00118 * 23 22 21
25 25 C 1.53006 * 109.47236 * 179.97438 * 24 23 22
26 26 C 1.53007 * 109.46948 * 60.00052 * 25 24 23
27 27 C 1.52996 * 109.46807 * 300.00230 * 26 25 24
28 28 C 1.52997 * 109.47558 * 60.00317 * 27 26 25
29 29 H 1.08998 * 109.46811 * 60.00625 * 1 2 3
30 30 H 1.08994 * 109.47201 * 180.02562 * 1 2 3
31 31 H 1.08995 * 109.47273 * 300.00977 * 1 2 3
32 32 H 1.08995 * 109.47553 * 294.99926 * 4 2 1
33 33 H 1.08998 * 109.47075 * 55.00372 * 4 2 1
34 34 H 1.00896 * 111.00427 * 56.05228 * 5 4 2
35 35 H 1.08998 * 109.46749 * 299.99469 * 6 5 4
36 36 H 1.08998 * 109.47510 * 59.99863 * 6 5 4
37 37 H 0.96996 * 120.00640 * 0.02562 * 9 7 6
38 38 H 0.97004 * 119.99617 * 179.97438 * 13 12 10
39 39 H 0.97000 * 120.00580 * 265.00723 * 18 2 1
40 40 H 1.09009 * 109.47352 * 300.27556 * 21 19 18
41 41 H 1.08998 * 109.47585 * 60.27844 * 21 19 18
42 42 H 1.08998 * 109.46957 * 30.00235 * 23 22 21
43 43 H 1.08998 * 109.47340 * 299.99441 * 24 23 22
44 44 H 1.09000 * 109.47181 * 60.00247 * 24 23 22
45 45 H 1.08998 * 109.47366 * 300.00301 * 25 24 23
46 46 H 1.09006 * 109.47367 * 179.97438 * 25 24 23
47 47 H 1.08994 * 109.46966 * 179.97438 * 26 25 24
48 48 H 1.08998 * 109.46883 * 59.99954 * 26 25 24
49 49 H 1.08998 * 109.47173 * 299.99546 * 27 26 25
50 50 H 1.09007 * 109.46594 * 179.97438 * 27 26 25
51 51 H 1.08997 * 109.47543 * 180.02562 * 28 27 26
52 52 H 1.08998 * 109.47508 * 59.99434 * 28 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.5301 0.0000 0.0000
3 1 1.8934 1.0277 0.0000
4 6 2.0399 -0.7211 1.2494
5 7 3.5040 -0.8265 1.1906
6 6 4.0302 -1.5157 2.3765
7 6 5.5315 -1.6052 2.2805
8 8 6.1066 -1.1404 1.3191
9 7 6.2351 -2.2011 3.2635
10 6 7.6010 -2.2831 3.1760
11 1 8.1129 -1.8698 2.3195
12 6 8.3264 -2.8978 4.1891
13 7 9.6355 -2.9759 4.1054
14 14 7.4408 -3.6131 5.6701
15 1 7.0905 -5.0308 5.4002
16 1 8.3268 -3.5400 6.8596
17 1 6.2011 -2.8373 5.9279
18 7 2.0185 -0.6908 -1.1960
19 6 2.1472 -0.0185 -2.3570
20 8 1.8591 1.1584 -2.4116
21 6 2.6554 -0.7275 -3.5859
22 8 2.7084 0.1918 -4.6786
23 6 3.1804 -0.3824 -5.8990
24 6 4.7016 -0.2387 -5.9755
25 6 5.2070 -0.8542 -7.2820
26 6 4.5660 -0.1306 -8.4681
27 6 3.0448 -0.2743 -8.3914
28 6 2.5394 0.3411 -7.0851
29 1 -0.3633 0.5137 -0.8900
30 1 -0.3633 -1.0276 0.0005
31 1 -0.3633 0.5140 0.8898
32 1 1.7516 -0.1581 2.1370
33 1 1.6053 -1.7197 1.2953
34 1 3.9280 0.0826 1.0826
35 1 3.7556 -0.9585 3.2721
36 1 3.6094 -2.5197 2.4299
37 1 5.7754 -2.5726 4.0325
38 1 10.1420 -3.4047 4.8129
39 1 2.2491 -1.6319 -1.1524
40 1 3.6539 -1.1197 -3.3923
41 1 1.9842 -1.5494 -3.8349
42 1 2.9132 -1.4387 -5.9307
43 1 4.9689 0.8175 -5.9439
44 1 5.1583 -0.7546 -5.1309
45 1 4.9398 -1.9104 -7.3134
46 1 6.2908 -0.7514 -7.3368
47 1 4.9264 -0.5687 -9.3988
48 1 4.8332 0.9257 -8.4367
49 1 2.7776 -1.3306 -8.4228
50 1 2.5881 0.2416 -9.2362
51 1 1.4557 0.2383 -7.0303
52 1 2.8066 1.3974 -7.0537
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
ZINC001275215036.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 06:26:58
Heat of formation + Delta-G solvation = -135.939560 kcal
Electronic energy + Delta-G solvation = -29565.540498 eV
Core-core repulsion = 25371.257547 eV
Total energy + Delta-G solvation = -4194.282951 eV
No. of doubly occupied orbitals = 66
Molecular weight
(most abundant/longest-lived isotopes) = 341.215 amu
Computer time = 0.64 seconds
Orbital eigenvalues (eV)
-40.03156 -39.17066 -38.49914 -36.94357 -36.22273 -34.34742 -33.73570 -31.71356
-31.45246 -31.07124 -30.12965 -27.79067 -27.45160 -25.04170 -23.19016 -22.90514
-22.48625 -21.58050 -20.57179 -20.09446 -19.40949 -18.39509 -17.44342 -17.09661
-16.39506 -16.10477 -15.89989 -15.64131 -15.15547 -14.92573 -14.72040 -14.63005
-14.38824 -14.06992 -13.96175 -13.84296 -13.39150 -13.24328 -13.21297 -12.76045
-12.69663 -12.61373 -12.48668 -12.29726 -12.25056 -12.09863 -11.74996 -11.71134
-11.43136 -11.32156 -11.18856 -11.00955 -10.90437 -10.80702 -10.62476 -10.58243
-10.00400 -9.74963 -9.64486 -9.60242 -9.07348 -8.64066 -8.25175 -8.19730
-6.33716 -3.85339 2.62482 2.75218 2.99611 3.03996 3.12554 3.70380
3.76208 4.00536 4.07803 4.40946 4.43670 4.63728 4.67895 4.74621
4.80407 4.87940 4.91376 4.96864 5.01299 5.03951 5.12754 5.17288
5.19514 5.20958 5.23404 5.31366 5.36080 5.49559 5.50145 5.52706
5.58022 5.68091 5.72873 5.75314 5.81678 5.86203 5.88108 5.88187
5.92742 6.01172 6.02372 6.11615 6.45944 6.94355 6.99295 7.46972
7.53896 7.55381 8.11623 8.26038 8.56762 8.92229 9.39484 9.79505
10.79366
Molecular weight = 341.21amu
Principal moments of inertia in cm(-1)
A = 0.010843 B = 0.002093 C = 0.001789
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 2581.707522 B =13372.307691 C =15648.131370
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C -0.161 4.161
2 C 0.147 3.853
3 H 0.087 0.913
4 C 0.048 3.952
5 N -0.667 5.667
6 C 0.048 3.952
7 C 0.449 3.551
8 O -0.564 6.564
9 N -0.590 5.590
10 C -0.305 4.305
11 H 0.119 0.881
12 C 0.041 3.959
13 N -0.881 5.881
14 Si 0.881 3.119
15 H -0.271 1.271
16 H -0.277 1.277
17 H -0.270 1.270
18 N -0.708 5.708
19 C 0.518 3.482
20 O -0.528 6.528
21 C 0.044 3.956
22 O -0.357 6.357
23 C 0.090 3.910
24 C -0.145 4.145
25 C -0.117 4.117
26 C -0.123 4.123
27 C -0.116 4.116
28 C -0.115 4.115
29 H 0.063 0.937
30 H 0.058 0.942
31 H 0.065 0.935
32 H 0.044 0.956
33 H 0.078 0.922
34 H 0.353 0.647
35 H 0.045 0.955
36 H 0.087 0.913
37 H 0.370 0.630
38 H 0.273 0.727
39 H 0.405 0.595
40 H 0.083 0.917
41 H 0.076 0.924
42 H 0.060 0.940
43 H 0.073 0.927
44 H 0.067 0.933
45 H 0.062 0.938
46 H 0.066 0.934
47 H 0.062 0.938
48 H 0.064 0.936
49 H 0.062 0.938
50 H 0.067 0.933
51 H 0.068 0.932
52 H 0.072 0.928
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges -19.323 0.899 -16.632 25.511
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.220 4.220
2 C 0.042 3.958
3 H 0.105 0.895
4 C -0.082 4.082
5 N -0.300 5.300
6 C -0.085 4.085
7 C 0.234 3.766
8 O -0.445 6.445
9 N -0.230 5.230
10 C -0.434 4.434
11 H 0.137 0.863
12 C -0.162 4.162
13 N -0.521 5.521
14 Si 0.707 3.293
15 H -0.196 1.196
16 H -0.202 1.202
17 H -0.196 1.196
18 N -0.361 5.361
19 C 0.303 3.697
20 O -0.404 6.404
21 C -0.036 4.036
22 O -0.275 6.275
23 C 0.032 3.968
24 C -0.182 4.182
25 C -0.155 4.155
26 C -0.161 4.161
27 C -0.153 4.153
28 C -0.153 4.153
29 H 0.082 0.918
30 H 0.077 0.923
31 H 0.084 0.916
32 H 0.062 0.938
33 H 0.096 0.904
34 H 0.173 0.827
35 H 0.063 0.937
36 H 0.105 0.895
37 H 0.203 0.797
38 H 0.084 0.916
39 H 0.240 0.760
40 H 0.101 0.899
41 H 0.094 0.906
42 H 0.078 0.922
43 H 0.092 0.908
44 H 0.086 0.914
45 H 0.081 0.919
46 H 0.085 0.915
47 H 0.081 0.919
48 H 0.083 0.917
49 H 0.081 0.919
50 H 0.085 0.915
51 H 0.087 0.913
52 H 0.091 0.909
Dipole moment (debyes) X Y Z Total
from point charges -18.973 1.642 -17.603 25.933
hybrid contribution -0.340 -0.796 2.149 2.317
sum -19.313 0.845 -15.454 24.749
Atomic orbital electron populations
1.22143 0.94956 1.02213 1.02702 1.20924 0.93584 0.97540 0.83741
0.89500 1.22319 0.89348 1.00743 0.95782 1.60055 1.03201 1.27105
1.39660 1.21712 0.95059 1.00022 0.91695 1.19704 0.87944 0.84147
0.84768 1.90576 1.70982 1.50403 1.32584 1.42210 1.07880 1.53082
1.19818 1.19429 0.79570 1.35990 1.08402 0.86301 1.25049 0.95453
0.97665 0.98057 1.69623 1.05077 1.43241 1.34120 0.86659 0.79607
0.79977 0.83040 1.19606 1.20185 1.19562 1.45692 1.67300 1.14437
1.08652 1.20173 0.77893 0.87397 0.84257 1.90709 1.49126 1.15688
1.84840 1.22239 1.01875 0.92414 0.87118 1.88031 1.86490 1.35416
1.17604 1.22176 0.93461 0.96866 0.84297 1.22070 0.95572 1.02568
0.98038 1.21501 1.00333 0.99633 0.94024 1.21587 0.95238 1.01081
0.98167 1.21507 0.95906 1.01051 0.96886 1.21575 1.00334 1.00813
0.92581 0.91818 0.92311 0.91591 0.93792 0.90393 0.82728 0.93704
0.89462 0.79710 0.91595 0.75973 0.89874 0.90630 0.92161 0.90846
0.91447 0.91933 0.91491 0.91893 0.91727 0.91895 0.91452 0.91340
0.90902
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.16 -1.42 9.48 37.16 0.35 -1.06 16
2 C 0.15 1.56 2.81 -67.93 -0.19 1.37 16
3 H 0.09 1.11 7.58 -51.93 -0.39 0.71 16
4 C 0.05 0.49 5.24 -3.82 -0.02 0.47 16
5 N -0.67 -9.45 5.85 -115.05 -0.67 -10.13 16
6 C 0.05 0.74 6.17 -4.97 -0.03 0.71 16
7 C 0.45 9.88 7.82 -10.99 -0.09 9.80 16
8 O -0.56 -14.86 15.78 5.55 0.09 -14.77 16
9 N -0.59 -13.76 5.09 -10.96 -0.06 -13.82 16
10 C -0.31 -8.54 10.16 -14.93 -0.15 -8.70 16
11 H 0.12 3.82 7.39 -52.48 -0.39 3.44 16
12 C 0.04 1.11 5.50 -17.47 -0.10 1.01 16
13 N -0.88 -25.56 13.96 55.50 0.77 -24.79 16
14 Si 0.88 19.09 27.74 -169.99 -4.71 14.37 16
15 H -0.27 -5.73 7.11 56.52 0.40 -5.33 16
16 H -0.28 -6.03 7.11 56.52 0.40 -5.63 16
17 H -0.27 -5.31 6.52 56.52 0.37 -4.95 16
18 N -0.71 -7.73 4.80 -53.80 -0.26 -7.99 16
19 C 0.52 6.87 7.87 -10.98 -0.09 6.78 16
20 O -0.53 -9.20 16.16 -14.36 -0.23 -9.43 16
21 C 0.04 0.47 6.23 36.01 0.22 0.69 16
22 O -0.36 -4.54 9.52 -55.14 -0.53 -5.07 16
23 C 0.09 0.86 2.61 -26.73 -0.07 0.79 16
24 C -0.14 -1.37 5.54 -26.73 -0.15 -1.51 16
25 C -0.12 -0.89 5.92 -26.72 -0.16 -1.05 16
26 C -0.12 -0.83 5.92 -26.73 -0.16 -0.99 16
27 C -0.12 -0.76 5.92 -26.73 -0.16 -0.92 16
28 C -0.12 -0.98 5.57 -26.73 -0.15 -1.13 16
29 H 0.06 0.66 8.14 -51.93 -0.42 0.24 16
30 H 0.06 0.44 8.14 -51.93 -0.42 0.02 16
31 H 0.07 0.50 8.14 -51.93 -0.42 0.07 16
32 H 0.04 0.41 8.14 -51.93 -0.42 -0.02 16
33 H 0.08 0.69 8.14 -51.93 -0.42 0.27 16
34 H 0.35 5.60 8.43 -39.29 -0.33 5.27 16
35 H 0.04 0.61 8.14 -51.93 -0.42 0.19 16
36 H 0.09 1.17 8.14 -51.93 -0.42 0.74 16
37 H 0.37 7.54 5.76 -40.82 -0.24 7.30 16
38 H 0.27 7.32 8.31 -40.82 -0.34 6.98 16
39 H 0.41 3.78 8.60 -40.82 -0.35 3.43 16
40 H 0.08 0.84 7.87 -51.92 -0.41 0.43 16
41 H 0.08 0.59 8.09 -51.93 -0.42 0.18 16
42 H 0.06 0.49 7.79 -51.93 -0.40 0.08 16
43 H 0.07 0.77 8.14 -51.93 -0.42 0.35 16
44 H 0.07 0.67 7.87 -51.93 -0.41 0.26 16
45 H 0.06 0.44 8.14 -51.93 -0.42 0.01 16
46 H 0.07 0.50 8.14 -51.93 -0.42 0.07 16
47 H 0.06 0.37 8.14 -51.93 -0.42 -0.06 16
48 H 0.06 0.46 8.14 -51.93 -0.42 0.03 16
49 H 0.06 0.38 8.14 -51.93 -0.42 -0.04 16
50 H 0.07 0.39 8.14 -51.93 -0.42 -0.04 16
51 H 0.07 0.57 8.14 -51.93 -0.42 0.15 16
52 H 0.07 0.68 8.14 -51.93 -0.42 0.26 16
LS Contribution 420.36 15.07 6.33 6.33
Total: -1.00 -35.14 420.36 -9.46 -44.60
By element:
Atomic # 1 Polarization: 23.69 SS G_CDS: -9.27 Total: 14.42 kcal
Atomic # 6 Polarization: 7.19 SS G_CDS: -0.93 Total: 6.26 kcal
Atomic # 7 Polarization: -56.51 SS G_CDS: -0.21 Total: -56.72 kcal
Atomic # 8 Polarization: -28.60 SS G_CDS: -0.67 Total: -29.27 kcal
Atomic # 14 Polarization: 19.09 SS G_CDS: -4.71 Total: 14.37 kcal
Total LS contribution 6.33 Total: 6.33 kcal
Total: -35.14 -9.46 -44.60 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.
ZINC001275215036.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 -91.335 kcal
(2) G-P(sol) polarization free energy of solvation -35.143 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system -126.477 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.462 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.605 kcal
(6) G-S(sol) free energy of system = (1) + (5) -135.940 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 0.64 seconds