Wall clock time and date at job start Sat Apr 11 2020 06:46:30
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.50703 * 1
3 3 C 1.38237 * 119.99416 * 2 1
4 4 C 1.38234 * 119.99412 * 179.97438 * 3 2 1
5 5 C 1.38236 * 120.00190 * 0.02562 * 4 3 2
6 6 C 1.50711 * 119.99825 * 179.97438 * 5 4 3
7 7 C 1.50698 * 109.46978 * 89.99959 * 6 5 4
8 8 O 1.21283 * 120.00287 * 359.97438 * 7 6 5
9 9 N 1.34770 * 119.99828 * 180.02562 * 7 6 5
10 10 C 1.46501 * 120.00044 * 179.97438 * 9 7 6
11 11 H 1.08996 * 108.96871 * 33.92067 * 10 9 7
12 12 C 1.52525 * 109.04189 * 274.99559 * 10 9 7
13 13 C 1.53817 * 111.70096 * 185.38963 * 12 10 9
14 14 C 1.51626 * 107.02221 * 57.54435 * 13 12 10
15 15 H 1.08997 * 110.50424 * 55.19860 * 14 13 12
16 16 C 1.53059 * 116.67684 * 183.74058 * 14 13 12
17 17 N 1.48660 * 103.44270 * 155.88509 * 16 14 13
18 18 C 1.36941 * 125.78119 * 155.62542 * 17 16 14
19 19 H 1.08003 * 119.99664 * 179.97438 * 18 17 16
20 20 C 1.38810 * 119.99881 * 359.97438 * 18 17 16
21 21 N 1.31620 * 120.00114 * 0.02562 * 20 18 17
22 22 Si 1.86806 * 119.99853 * 180.02562 * 20 18 17
23 23 H 1.48498 * 109.47128 * 59.99373 * 22 20 18
24 24 H 1.48495 * 109.47120 * 179.97438 * 22 20 18
25 25 H 1.48506 * 109.46689 * 300.00031 * 22 20 18
26 26 C 1.48895 * 108.52288 * 335.22832 * 17 16 14
27 27 C 1.52328 * 105.53885 * 357.42989 * 26 17 16
28 28 H 1.09004 * 109.87734 * 271.58291 * 27 26 17
29 29 C 1.51439 * 117.45352 * 145.05752 * 27 26 17
30 30 C 1.38232 * 120.00110 * 0.24933 * 5 4 3
31 31 C 1.38231 * 120.00059 * 180.02562 * 2 1 3
32 32 H 1.08997 * 109.47006 * 90.00557 * 1 2 3
33 33 H 1.09002 * 109.46903 * 209.99666 * 1 2 3
34 34 H 1.08996 * 109.47493 * 329.99694 * 1 2 3
35 35 H 1.07998 * 120.00611 * 0.02562 * 3 2 1
36 36 H 1.08000 * 120.00037 * 179.97438 * 4 3 2
37 37 H 1.09002 * 109.47087 * 329.99851 * 6 5 4
38 38 H 1.09001 * 109.46749 * 209.99619 * 6 5 4
39 39 H 0.97005 * 120.00063 * 359.97438 * 9 7 6
40 40 H 1.09003 * 109.03179 * 64.82598 * 12 10 9
41 41 H 1.08997 * 108.98114 * 305.92589 * 12 10 9
42 42 H 1.08997 * 109.94324 * 176.91863 * 13 12 10
43 43 H 1.08998 * 109.94516 * 298.16537 * 13 12 10
44 44 H 1.08998 * 110.63008 * 274.37883 * 16 14 13
45 45 H 1.09002 * 110.62692 * 37.39655 * 16 14 13
46 46 H 0.97001 * 119.99435 * 180.02562 * 21 20 18
47 47 H 1.09001 * 110.55437 * 238.31562 * 26 17 16
48 48 H 1.08999 * 110.10436 * 116.30390 * 26 17 16
49 49 H 1.08996 * 109.94185 * 61.33813 * 29 27 26
50 50 H 1.09005 * 110.00685 * 300.13988 * 29 27 26
51 51 H 1.07996 * 120.00054 * 179.72775 * 30 5 4
52 52 H 1.07997 * 120.00287 * 359.97438 * 31 2 1
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.5070 0.0000 0.0000
3 6 2.1981 1.1972 0.0000
4 6 3.5804 1.1972 0.0005
5 6 4.2716 0.0001 0.0005
6 6 5.7788 0.0002 0.0005
7 6 6.2811 0.0002 1.4213
8 8 5.4929 0.0007 2.3432
9 7 7.6061 0.0008 1.6676
10 6 8.0944 0.0015 3.0488
11 1 7.4063 -0.5739 3.6682
12 6 8.1414 1.4379 3.5594
13 6 8.7716 1.5234 4.9599
14 6 10.1591 0.9204 4.8599
15 1 10.7364 1.4179 4.0806
16 6 10.9569 0.8123 6.1616
17 7 11.9163 -0.2900 5.8885
18 6 13.1226 -0.4739 6.5100
19 1 13.7483 -1.3110 6.2375
20 6 13.5412 0.4155 7.4900
21 7 12.7784 1.4354 7.8223
22 14 15.1871 0.1650 8.3373
23 1 16.2750 0.2058 7.3273
24 1 15.3988 1.2432 9.3361
25 1 15.1960 -1.1530 9.0214
26 6 11.3408 -1.1669 4.8317
27 6 9.9664 -0.5868 4.5238
28 1 9.2221 -1.0233 5.1899
29 6 9.4820 -0.6617 3.0909
30 6 3.5806 -1.1970 -0.0052
31 6 2.1982 -1.1971 -0.0005
32 1 -0.3633 -0.0001 -1.0276
33 1 -0.3633 -0.8900 0.5138
34 1 -0.3634 0.8899 0.5138
35 1 1.6581 2.1325 0.0004
36 1 4.1204 2.1325 0.0001
37 1 6.1420 0.8902 -0.5133
38 1 6.1421 -0.8898 -0.5133
39 1 8.2364 0.0009 0.9303
40 1 7.1265 1.8335 3.6016
41 1 8.7293 2.0411 2.8676
42 1 8.8407 2.5650 5.2736
43 1 8.1712 0.9585 5.6730
44 1 11.4881 1.7422 6.3648
45 1 10.3016 0.5537 6.9933
46 1 13.0712 2.0571 8.5069
47 1 11.2442 -2.1903 5.1943
48 1 11.9680 -1.1408 3.9406
49 1 9.4058 -1.7031 2.7782
50 1 10.1698 -0.1265 2.4362
51 1 4.1206 -2.1323 -0.0097
52 1 1.6583 -2.1324 -0.0005
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)
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
ZINC000980152759.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:46:30
Heat of formation + Delta-G solvation = -3.164197 kcal
Electronic energy + Delta-G solvation = -29262.609189 eV
Core-core repulsion = 25419.507162 eV
Total energy + Delta-G solvation = -3843.102027 eV
No. of doubly occupied orbitals = 65
Molecular weight
(most abundant/longest-lived isotopes) = 342.205 amu
Computer time = 0.85 seconds
Orbital eigenvalues (eV)
-40.41763 -39.53075 -38.59299 -35.66223 -34.01848 -33.50993 -31.84017 -31.13192
-30.83127 -29.86128 -28.56312 -26.35332 -25.60260 -24.72693 -22.95129 -22.01351
-21.67194 -20.64062 -20.39891 -19.75132 -19.29880 -17.22428 -16.68906 -16.60525
-16.00837 -15.72304 -15.19839 -14.92411 -14.86882 -14.41439 -14.22519 -14.04590
-13.80233 -13.71445 -13.63381 -13.32075 -13.04872 -12.76491 -12.46051 -12.29473
-12.12099 -12.07105 -11.93315 -11.84746 -11.60182 -11.36045 -11.30942 -11.17245
-10.91411 -10.82451 -10.63139 -10.45950 -10.06935 -10.04068 -9.86848 -9.80331
-9.55300 -9.23951 -9.10376 -8.94598 -8.70633 -8.61948 -6.79431 -5.60658
-3.01488 1.08645 1.22553 2.57124 3.06538 3.42723 3.46730 3.48255
3.61249 4.37647 4.51193 4.52937 4.66941 4.73223 4.79046 4.83898
4.89508 4.93916 5.14398 5.18612 5.22935 5.26609 5.28544 5.35472
5.41635 5.48967 5.51745 5.54583 5.66641 5.71222 5.71593 5.82103
5.85822 5.86998 6.08802 6.18332 6.18611 6.22252 6.28159 6.30789
6.36701 6.38138 6.42636 6.51988 6.59438 6.66947 6.71062 6.80090
6.89512 7.07402 7.44429 7.58721 7.64613 7.90086 8.17596 8.37421
9.24424 9.85269 10.52653 11.43690
Molecular weight = 342.20amu
Principal moments of inertia in cm(-1)
A = 0.035868 B = 0.001765 C = 0.001735
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 780.446890 B =15856.125305 C =16130.066544
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C -0.114 4.114
2 C -0.089 4.089
3 C -0.120 4.120
4 C -0.101 4.101
5 C -0.070 4.070
6 C -0.103 4.103
7 C 0.512 3.488
8 O -0.530 6.530
9 N -0.728 5.728
10 C 0.149 3.851
11 H 0.084 0.916
12 C -0.123 4.123
13 C -0.099 4.099
14 C -0.114 4.114
15 H 0.068 0.932
16 C 0.164 3.836
17 N -0.608 5.608
18 C -0.221 4.221
19 H 0.076 0.924
20 C -0.013 4.013
21 N -0.923 5.923
22 Si 0.913 3.087
23 H -0.286 1.286
24 H -0.290 1.290
25 H -0.286 1.286
26 C 0.144 3.856
27 C -0.105 4.105
28 H 0.069 0.931
29 C -0.105 4.105
30 C -0.101 4.101
31 C -0.120 4.120
32 H 0.068 0.932
33 H 0.065 0.935
34 H 0.065 0.935
35 H 0.120 0.880
36 H 0.119 0.881
37 H 0.102 0.898
38 H 0.101 0.899
39 H 0.398 0.602
40 H 0.060 0.940
41 H 0.059 0.941
42 H 0.065 0.935
43 H 0.062 0.938
44 H 0.104 0.896
45 H 0.013 0.987
46 H 0.252 0.748
47 H 0.025 0.975
48 H 0.023 0.977
49 H 0.064 0.936
50 H 0.068 0.932
51 H 0.119 0.881
52 H 0.120 0.880
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges -20.031 -1.566 -20.520 28.719
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.171 4.171
2 C -0.089 4.089
3 C -0.138 4.138
4 C -0.119 4.119
5 C -0.070 4.070
6 C -0.143 4.143
7 C 0.298 3.702
8 O -0.406 6.406
9 N -0.379 5.379
10 C 0.045 3.955
11 H 0.102 0.898
12 C -0.162 4.162
13 C -0.137 4.137
14 C -0.134 4.134
15 H 0.087 0.913
16 C 0.040 3.960
17 N -0.335 5.335
18 C -0.349 4.349
19 H 0.094 0.906
20 C -0.208 4.208
21 N -0.572 5.572
22 Si 0.744 3.256
23 H -0.213 1.213
24 H -0.216 1.216
25 H -0.213 1.213
26 C 0.018 3.982
27 C -0.124 4.124
28 H 0.087 0.913
29 C -0.144 4.144
30 C -0.119 4.119
31 C -0.138 4.138
32 H 0.087 0.913
33 H 0.084 0.916
34 H 0.084 0.916
35 H 0.138 0.862
36 H 0.137 0.863
37 H 0.120 0.880
38 H 0.120 0.880
39 H 0.232 0.768
40 H 0.078 0.922
41 H 0.077 0.923
42 H 0.084 0.916
43 H 0.081 0.919
44 H 0.122 0.878
45 H 0.031 0.969
46 H 0.062 0.938
47 H 0.043 0.957
48 H 0.041 0.959
49 H 0.082 0.918
50 H 0.087 0.913
51 H 0.137 0.863
52 H 0.138 0.862
Dipole moment (debyes) X Y Z Total
from point charges -20.474 -1.477 -19.673 28.432
hybrid contribution 1.208 -0.203 0.152 1.235
sum -19.266 -1.680 -19.520 27.478
Atomic orbital electron populations
1.20788 0.91330 1.02808 1.02143 1.19915 0.96012 0.93277 0.99744
1.20971 0.94067 0.97505 1.01246 1.20956 0.94249 0.97356 0.99339
1.19247 0.93574 0.93606 1.00569 1.20695 0.93376 1.05180 0.95087
1.20512 0.82401 0.76222 0.91113 1.90741 1.55181 1.48525 1.46172
1.46208 1.06773 1.75597 1.09317 1.21101 0.95669 0.97024 0.81734
0.89789 1.21890 1.00859 0.95776 0.97665 1.21256 0.96035 1.01201
0.95255 1.22391 0.95249 0.94112 1.01613 0.91329 1.21624 0.89729
0.94812 0.89844 1.45285 1.18427 1.27047 1.42753 1.19054 0.93083
1.08483 1.14287 0.90647 1.24769 0.98789 0.97825 0.99447 1.69809
1.36779 1.16232 1.34356 0.86170 0.83329 0.77030 0.79068 1.21252
1.21608 1.21302 1.21069 0.91237 0.93642 0.92210 1.21717 0.97913
0.98208 0.94605 0.91267 1.21534 0.94700 1.00091 0.98051 1.20953
0.94249 0.97338 0.99402 1.20967 0.94064 0.97499 1.01268 0.91295
0.91586 0.91578 0.86165 0.86254 0.88024 0.88037 0.76797 0.92157
0.92272 0.91648 0.91915 0.87800 0.96933 0.93773 0.95735 0.95899
0.91778 0.91343 0.86263 0.86176
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.11 -0.73 10.01 36.01 0.36 -0.37 16
2 C -0.09 -0.80 5.88 -104.62 -0.62 -1.42 16
3 C -0.12 -1.11 9.70 -39.58 -0.38 -1.49 16
4 C -0.10 -0.98 9.70 -39.58 -0.38 -1.37 16
5 C -0.07 -0.68 4.37 -104.62 -0.46 -1.14 16
6 C -0.10 -0.81 5.83 -29.03 -0.17 -0.98 16
7 C 0.51 5.85 7.78 -10.99 -0.09 5.77 16
8 O -0.53 -8.50 15.66 5.55 0.09 -8.41 16
9 N -0.73 -7.39 5.03 -53.61 -0.27 -7.66 16
10 C 0.15 1.91 2.44 -67.73 -0.17 1.74 16
11 H 0.08 1.26 7.51 -51.93 -0.39 0.87 16
12 C -0.12 -1.68 5.37 -26.55 -0.14 -1.83 16
13 C -0.10 -1.57 5.71 -27.00 -0.15 -1.72 16
14 C -0.11 -2.12 2.51 -89.89 -0.23 -2.35 16
15 H 0.07 1.30 8.01 -51.93 -0.42 0.88 16
16 C 0.16 3.99 5.53 -2.82 -0.02 3.97 16
17 N -0.61 -15.73 3.46 -169.87 -0.59 -16.32 16
18 C -0.22 -6.22 10.30 -15.06 -0.16 -6.38 16
19 H 0.08 2.05 7.83 -52.48 -0.41 1.64 16
20 C -0.01 -0.37 5.49 -17.43 -0.10 -0.47 16
21 N -0.92 -27.16 10.11 55.49 0.56 -26.60 16
22 Si 0.91 22.26 29.55 -169.99 -5.02 17.24 16
23 H -0.29 -6.84 7.11 56.52 0.40 -6.44 16
24 H -0.29 -7.12 7.11 56.52 0.40 -6.71 16
25 H -0.29 -6.90 7.11 56.52 0.40 -6.50 16
26 C 0.14 2.99 7.04 -3.06 -0.02 2.97 16
27 C -0.10 -1.80 2.28 -90.36 -0.21 -2.01 16
28 H 0.07 1.23 7.92 -51.93 -0.41 0.82 16
29 C -0.11 -1.38 5.16 -27.08 -0.14 -1.52 16
30 C -0.10 -0.98 9.70 -39.58 -0.38 -1.37 16
31 C -0.12 -1.10 9.70 -39.58 -0.38 -1.49 16
32 H 0.07 0.36 8.14 -51.93 -0.42 -0.06 16
33 H 0.07 0.38 8.09 -51.93 -0.42 -0.04 16
34 H 0.07 0.38 8.09 -51.93 -0.42 -0.04 16
35 H 0.12 0.88 8.06 -52.49 -0.42 0.45 16
36 H 0.12 0.99 8.06 -52.49 -0.42 0.57 16
37 H 0.10 0.58 8.09 -51.93 -0.42 0.16 16
38 H 0.10 0.57 8.08 -51.93 -0.42 0.15 16
39 H 0.40 2.85 8.61 -40.82 -0.35 2.50 16
40 H 0.06 0.87 8.14 -51.93 -0.42 0.44 16
41 H 0.06 0.76 8.14 -51.93 -0.42 0.34 16
42 H 0.06 1.02 8.14 -51.93 -0.42 0.60 16
43 H 0.06 1.02 7.98 -51.93 -0.41 0.61 16
44 H 0.10 2.85 5.93 -51.93 -0.31 2.54 16
45 H 0.01 0.33 8.14 -51.93 -0.42 -0.09 16
46 H 0.25 7.20 8.31 -40.82 -0.34 6.86 16
47 H 0.02 0.51 8.14 -51.93 -0.42 0.09 16
48 H 0.02 0.47 8.08 -51.93 -0.42 0.05 16
49 H 0.06 0.77 8.14 -51.93 -0.42 0.35 16
50 H 0.07 0.85 8.04 -51.93 -0.42 0.44 16
51 H 0.12 0.98 8.06 -52.49 -0.42 0.56 16
52 H 0.12 0.87 8.06 -52.49 -0.42 0.45 16
LS Contribution 409.44 15.07 6.17 6.17
Total: -1.00 -33.65 409.44 -11.89 -45.54
By element:
Atomic # 1 Polarization: 10.47 SS G_CDS: -9.00 Total: 1.47 kcal
Atomic # 6 Polarization: -7.61 SS G_CDS: -3.82 Total: -11.44 kcal
Atomic # 7 Polarization: -50.28 SS G_CDS: -0.30 Total: -50.57 kcal
Atomic # 8 Polarization: -8.50 SS G_CDS: 0.09 Total: -8.41 kcal
Atomic # 14 Polarization: 22.26 SS G_CDS: -5.02 Total: 17.24 kcal
Total LS contribution 6.17 Total: 6.17 kcal
Total: -33.65 -11.89 -45.54 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.
ZINC000980152759.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 42.374 kcal
(2) G-P(sol) polarization free energy of solvation -33.646 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system 8.728 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.892 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -45.538 kcal
(6) G-S(sol) free energy of system = (1) + (5) -3.164 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 0.86 seconds