Wall clock time and date at job start Sat Apr 11 2020 06:44:52
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.53003 * 1
3 3 C 1.50696 * 109.47175 * 2 1
4 4 C 1.38234 * 120.00045 * 90.00065 * 3 2 1
5 5 C 1.38237 * 120.00290 * 180.02562 * 4 3 2
6 6 C 1.38240 * 119.99721 * 359.94747 * 5 4 3
7 7 C 1.50695 * 120.00026 * 180.02562 * 6 5 4
8 8 C 1.50698 * 109.47112 * 90.00241 * 7 6 5
9 9 O 1.21271 * 120.00107 * 359.97438 * 8 7 6
10 10 N 1.34780 * 119.99605 * 180.02562 * 8 7 6
11 11 C 1.46951 * 120.43847 * 179.97438 * 10 8 7
12 12 C 1.53593 * 108.52798 * 125.72437 * 11 10 8
13 13 H 1.08998 * 110.59777 * 172.77123 * 12 11 10
14 14 C 1.53369 * 108.68595 * 294.36273 * 12 11 10
15 15 C 1.53007 * 109.27201 * 60.30010 * 14 12 11
16 16 C 1.53007 * 109.80054 * 60.56533 * 15 14 12
17 17 C 1.53370 * 109.25252 * 299.42421 * 16 15 14
18 18 H 1.08998 * 110.59526 * 178.14430 * 17 16 15
19 19 C 1.46941 * 120.44855 * 359.97438 * 10 8 7
20 20 N 1.47566 * 107.98751 * 56.37731 * 17 16 15
21 21 C 1.36936 * 121.60911 * 122.25605 * 20 17 16
22 22 H 1.08006 * 119.99698 * 206.85727 * 21 20 17
23 23 C 1.38809 * 120.00342 * 26.85503 * 21 20 17
24 24 N 1.31619 * 119.99678 * 18.53643 * 23 21 20
25 25 Si 1.86805 * 120.00262 * 198.53564 * 23 21 20
26 26 H 1.48492 * 109.46901 * 60.00133 * 25 23 21
27 27 H 1.48498 * 109.47028 * 180.02562 * 25 23 21
28 28 H 1.48500 * 109.47100 * 300.00042 * 25 23 21
29 29 C 1.38231 * 119.99730 * 0.28073 * 6 5 4
30 30 C 1.38235 * 119.99864 * 269.97452 * 3 2 1
31 31 H 1.08992 * 109.47352 * 300.00571 * 1 2 3
32 32 H 1.09000 * 109.46831 * 60.00471 * 1 2 3
33 33 H 1.09002 * 109.46957 * 179.97438 * 1 2 3
34 34 H 1.09003 * 109.46949 * 239.99671 * 2 1 3
35 35 H 1.08999 * 109.46624 * 120.00032 * 2 1 3
36 36 H 1.08003 * 120.00237 * 359.96828 * 4 3 2
37 37 H 1.07993 * 120.00047 * 179.97438 * 5 4 3
38 38 H 1.09008 * 109.47078 * 329.99847 * 7 6 5
39 39 H 1.09002 * 109.47190 * 210.00181 * 7 6 5
40 40 H 1.09001 * 109.62855 * 6.01378 * 11 10 8
41 41 H 1.09001 * 109.78526 * 245.53275 * 11 10 8
42 42 H 1.09006 * 109.50216 * 180.23662 * 14 12 11
43 43 H 1.08994 * 109.50468 * 300.35414 * 14 12 11
44 44 H 1.09004 * 109.41161 * 300.46412 * 15 14 12
45 45 H 1.09001 * 109.36151 * 180.63776 * 15 14 12
46 46 H 1.08994 * 109.51032 * 179.48505 * 16 15 14
47 47 H 1.09001 * 109.53907 * 59.37277 * 16 15 14
48 48 H 1.09010 * 109.62041 * 353.98033 * 19 10 8
49 49 H 1.09000 * 109.62481 * 114.54117 * 19 10 8
50 50 H 0.97007 * 120.00025 * 179.97438 * 24 23 21
51 51 H 1.08003 * 120.00005 * 179.72633 * 29 6 5
52 52 H 1.07993 * 120.00334 * 0.02562 * 30 3 2
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.0324 1.4208 0.0000
4 6 2.2628 2.0724 1.1971
5 6 2.7231 3.3759 1.1972
6 6 2.9540 4.0274 0.0000
7 6 3.4563 5.4481 0.0000
8 6 4.9633 5.4481 0.0006
9 8 5.5697 4.3979 0.0014
10 7 5.6371 6.6154 0.0012
11 6 7.1066 6.6267 0.0023
12 6 7.5890 7.4809 -1.1795
13 1 8.6695 7.6179 -1.1378
14 6 7.1833 6.7866 -2.4855
15 6 5.6599 6.6522 -2.5338
16 6 5.0154 8.0391 -2.4875
17 6 5.4125 8.7379 -1.1813
18 1 4.9911 9.7423 -1.1404
19 6 4.9123 7.8936 0.0007
20 7 6.8862 8.7768 -1.1150
21 6 7.5685 9.9584 -0.9993
22 1 8.5783 10.0378 -1.3741
23 6 6.9623 11.0542 -0.4003
24 7 5.9086 10.8839 0.3697
25 14 7.6422 12.7703 -0.6870
26 1 7.5977 13.0843 -2.1377
27 1 6.8226 13.7585 0.0592
28 1 9.0471 12.8322 -0.2099
29 6 2.7191 3.3773 -1.1971
30 6 2.2632 2.0722 -1.1972
31 1 -0.3633 0.5139 0.8898
32 1 -0.3633 0.5138 -0.8900
33 1 -0.3633 -1.0277 -0.0005
34 1 1.8933 -0.5139 -0.8900
35 1 1.8933 -0.5138 0.8900
36 1 2.0823 1.5635 2.1325
37 1 2.9027 3.8851 2.1325
38 1 3.0926 5.9622 0.8898
39 1 3.0933 5.9619 -0.8902
40 1 7.4807 5.6086 -0.1050
41 1 7.4715 7.0543 0.9362
42 1 7.5233 7.3803 -3.3341
43 1 7.6376 5.7968 -2.5278
44 1 5.3206 6.0662 -1.6795
45 1 5.3706 6.1487 -3.4563
46 1 3.9310 7.9383 -2.5314
47 1 5.3607 8.6300 -3.3359
48 1 3.8436 7.7083 -0.1078
49 1 5.1020 8.4223 0.9348
50 1 5.4851 11.6496 0.7886
51 1 2.8953 3.8877 -2.1324
52 1 2.0833 1.5632 -2.1324
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
ZINC001039934872.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:44:52
Heat of formation + Delta-G solvation = 26.153867 kcal
Electronic energy + Delta-G solvation = -31346.085944 eV
Core-core repulsion = 27504.255243 eV
Total energy + Delta-G solvation = -3841.830701 eV
No. of doubly occupied orbitals = 65
Molecular weight
(most abundant/longest-lived isotopes) = 342.205 amu
Computer time = 1.12 seconds
Orbital eigenvalues (eV)
-40.68109 -39.37712 -37.61404 -34.66194 -33.83098 -32.97025 -32.12009 -31.05595
-30.86555 -30.24383 -27.72823 -26.40089 -25.47125 -23.72551 -23.13396 -22.69381
-21.64224 -21.25150 -20.60689 -19.26905 -18.25486 -17.11018 -16.33179 -16.04629
-15.67074 -15.36849 -15.03435 -14.75662 -14.57823 -14.09626 -13.85931 -13.81335
-13.55159 -13.41569 -13.23076 -13.14321 -12.87706 -12.75025 -12.56847 -12.19601
-12.09982 -11.96526 -11.81396 -11.71861 -11.35108 -11.28126 -11.06017 -10.96974
-10.82942 -10.80709 -10.72499 -10.37013 -10.18700 -10.11759 -10.04596 -9.60251
-9.51274 -9.20832 -8.83090 -8.68167 -8.42308 -7.75060 -6.65709 -5.90304
-3.29418 1.37510 1.49558 3.06017 3.33510 3.39548 3.41223 3.51341
3.91062 4.41746 4.63153 4.73592 4.85675 4.99561 5.03826 5.07029
5.17907 5.19700 5.27122 5.35620 5.37633 5.50655 5.62763 5.63416
5.67024 5.69049 5.72126 5.74760 5.82051 5.84974 5.93541 5.96676
6.03256 6.12489 6.13712 6.24913 6.35614 6.40088 6.45795 6.46953
6.51941 6.60349 6.63466 6.66994 6.69768 6.72887 6.85931 7.01637
7.06804 7.11289 7.23401 7.48916 7.70648 8.03616 8.30821 8.66505
9.08530 9.70622 10.32162 11.18381
Molecular weight = 342.20amu
Principal moments of inertia in cm(-1)
A = 0.025872 B = 0.002656 C = 0.002595
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 1081.980583 B =10541.545010 C =10785.537677
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C -0.146 4.146
2 C -0.084 4.084
3 C -0.089 4.089
4 C -0.118 4.118
5 C -0.103 4.103
6 C -0.059 4.059
7 C -0.106 4.106
8 C 0.514 3.486
9 O -0.552 6.552
10 N -0.590 5.590
11 C 0.102 3.898
12 C 0.153 3.847
13 H 0.077 0.923
14 C -0.121 4.121
15 C -0.139 4.139
16 C -0.133 4.133
17 C 0.178 3.822
18 H 0.116 0.884
19 C 0.099 3.901
20 N -0.579 5.579
21 C -0.246 4.246
22 H 0.067 0.933
23 C 0.003 3.997
24 N -0.903 5.903
25 Si 0.904 3.096
26 H -0.285 1.285
27 H -0.291 1.291
28 H -0.284 1.284
29 C -0.102 4.102
30 C -0.117 4.117
31 H 0.055 0.945
32 H 0.055 0.945
33 H 0.055 0.945
34 H 0.069 0.931
35 H 0.069 0.931
36 H 0.118 0.882
37 H 0.119 0.881
38 H 0.100 0.900
39 H 0.101 0.899
40 H 0.076 0.924
41 H 0.063 0.937
42 H 0.054 0.946
43 H 0.046 0.954
44 H 0.085 0.915
45 H 0.031 0.969
46 H 0.043 0.957
47 H 0.054 0.946
48 H 0.062 0.938
49 H 0.087 0.913
50 H 0.255 0.745
51 H 0.119 0.881
52 H 0.118 0.882
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges -9.908 -14.759 -0.902 17.799
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.203 4.203
2 C -0.121 4.121
3 C -0.090 4.090
4 C -0.136 4.136
5 C -0.121 4.121
6 C -0.059 4.059
7 C -0.146 4.146
8 C 0.303 3.697
9 O -0.432 6.432
10 N -0.323 5.323
11 C -0.023 4.023
12 C 0.045 3.955
13 H 0.095 0.905
14 C -0.160 4.160
15 C -0.177 4.177
16 C -0.172 4.172
17 C 0.071 3.929
18 H 0.133 0.867
19 C -0.024 4.024
20 N -0.300 5.300
21 C -0.375 4.375
22 H 0.084 0.916
23 C -0.195 4.195
24 N -0.547 5.547
25 Si 0.734 3.266
26 H -0.211 1.211
27 H -0.217 1.217
28 H -0.210 1.210
29 C -0.120 4.120
30 C -0.135 4.135
31 H 0.074 0.926
32 H 0.074 0.926
33 H 0.074 0.926
34 H 0.088 0.912
35 H 0.088 0.912
36 H 0.136 0.864
37 H 0.137 0.863
38 H 0.118 0.882
39 H 0.119 0.881
40 H 0.095 0.905
41 H 0.081 0.919
42 H 0.073 0.927
43 H 0.065 0.935
44 H 0.103 0.897
45 H 0.050 0.950
46 H 0.062 0.938
47 H 0.072 0.928
48 H 0.080 0.920
49 H 0.105 0.895
50 H 0.065 0.935
51 H 0.137 0.863
52 H 0.136 0.864
Dipole moment (debyes) X Y Z Total
from point charges -9.645 -15.119 -0.787 17.950
hybrid contribution 0.371 1.157 -0.141 1.223
sum -9.274 -13.962 -0.928 16.786
Atomic orbital electron populations
1.21746 0.94902 1.01446 1.02218 1.20517 0.95999 0.92935 1.02653
1.19819 0.99808 0.96030 0.93295 1.20910 1.00565 0.94775 0.97326
1.20958 0.98851 0.94842 0.97443 1.19333 0.98266 0.94856 0.93471
1.20779 0.93608 0.95015 1.05172 1.20108 0.89395 0.83659 0.76513
1.90671 1.70081 1.31080 1.51330 1.48263 1.07373 1.06128 1.70511
1.22189 0.79276 1.00603 1.00199 1.21225 0.96199 0.86589 0.91445
0.90494 1.21901 0.96139 0.98802 0.99177 1.22101 0.96288 0.98006
1.01302 1.22137 0.97799 0.96473 1.00820 1.21093 0.84630 0.98640
0.88585 0.86683 1.22483 0.93788 0.82962 1.03186 1.44350 1.00704
1.00353 1.84554 1.19123 0.99109 0.87647 1.31656 0.91577 1.24937
0.98880 0.98825 0.96863 1.69932 1.23857 1.29565 1.31374 0.86348
0.78284 0.84232 0.77778 1.21149 1.21680 1.20992 1.20966 0.98749
0.94814 0.97471 1.20916 1.00530 0.94752 0.97332 0.92584 0.92581
0.92617 0.91216 0.91230 0.86430 0.86335 0.88174 0.88079 0.90519
0.91904 0.92736 0.93492 0.89675 0.95013 0.93807 0.92756 0.91987
0.89488 0.93506 0.86292 0.86414
Number of geometries 1
Number of calculations of the screened coulomb radii 1
The total number of SCF iterations 19.
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.08 9.90 37.16 0.37 -0.72 16
2 C -0.08 -0.72 6.06 -27.88 -0.17 -0.89 16
3 C -0.09 -1.01 5.13 -104.63 -0.54 -1.55 16
4 C -0.12 -1.39 9.70 -39.58 -0.38 -1.78 16
5 C -0.10 -1.35 9.70 -39.58 -0.38 -1.74 16
6 C -0.06 -0.82 4.37 -104.63 -0.46 -1.28 16
7 C -0.11 -1.48 4.14 -29.04 -0.12 -1.60 16
8 C 0.51 9.04 5.01 -10.98 -0.06 8.99 16
9 O -0.55 -11.30 15.43 5.57 0.09 -11.21 16
10 N -0.59 -10.87 1.65 -172.22 -0.28 -11.15 16
11 C 0.10 1.97 5.89 -3.49 -0.02 1.95 16
12 C 0.15 3.09 3.61 -66.85 -0.24 2.85 16
13 H 0.08 1.53 8.03 -51.93 -0.42 1.11 16
14 C -0.12 -2.24 5.58 -26.54 -0.15 -2.39 16
15 C -0.14 -2.49 5.55 -26.72 -0.15 -2.64 16
16 C -0.13 -2.46 5.59 -26.54 -0.15 -2.61 16
17 C 0.18 3.74 2.63 -66.83 -0.18 3.56 16
18 H 0.12 2.75 5.93 -51.93 -0.31 2.44 16
19 C 0.10 1.87 5.75 -3.47 -0.02 1.85 16
20 N -0.58 -13.22 2.54 -155.41 -0.39 -13.61 16
21 C -0.25 -6.26 9.73 -15.06 -0.15 -6.41 16
22 H 0.07 1.68 7.90 -52.48 -0.41 1.27 16
23 C 0.00 0.08 4.98 -17.43 -0.09 0.00 16
24 N -0.90 -23.48 10.73 55.49 0.60 -22.89 16
25 Si 0.90 20.62 29.59 -169.99 -5.03 15.59 16
26 H -0.28 -6.53 7.11 56.52 0.40 -6.13 16
27 H -0.29 -6.72 7.11 56.52 0.40 -6.32 16
28 H -0.28 -6.49 7.11 56.52 0.40 -6.09 16
29 C -0.10 -1.35 9.70 -39.58 -0.38 -1.74 16
30 C -0.12 -1.40 9.70 -39.58 -0.38 -1.78 16
31 H 0.06 0.42 8.14 -51.93 -0.42 0.00 16
32 H 0.06 0.43 8.14 -51.93 -0.42 0.00 16
33 H 0.05 0.34 8.14 -51.93 -0.42 -0.08 16
34 H 0.07 0.56 8.09 -51.93 -0.42 0.14 16
35 H 0.07 0.55 8.09 -51.93 -0.42 0.13 16
36 H 0.12 1.14 8.06 -52.48 -0.42 0.72 16
37 H 0.12 1.42 8.06 -52.49 -0.42 1.00 16
38 H 0.10 1.25 7.87 -51.92 -0.41 0.84 16
39 H 0.10 1.30 7.62 -51.93 -0.40 0.90 16
40 H 0.08 1.47 6.99 -51.93 -0.36 1.11 16
41 H 0.06 1.28 8.14 -51.93 -0.42 0.86 16
42 H 0.05 1.03 8.14 -51.93 -0.42 0.61 16
43 H 0.05 0.82 8.14 -51.93 -0.42 0.40 16
44 H 0.09 1.54 2.49 -52.44 -0.13 1.41 16
45 H 0.03 0.53 8.14 -51.93 -0.42 0.11 16
46 H 0.04 0.71 8.14 -51.93 -0.42 0.29 16
47 H 0.05 1.06 8.14 -51.93 -0.42 0.64 16
48 H 0.06 0.99 5.90 -51.92 -0.31 0.68 16
49 H 0.09 1.88 7.75 -51.93 -0.40 1.48 16
50 H 0.25 6.53 8.31 -40.82 -0.34 6.19 16
51 H 0.12 1.47 8.06 -52.48 -0.42 1.04 16
52 H 0.12 1.16 8.06 -52.49 -0.42 0.74 16
LS Contribution 394.44 15.07 5.94 5.94
Total: -1.00 -28.44 394.44 -11.34 -39.78
By element:
Atomic # 1 Polarization: 14.08 SS G_CDS: -8.62 Total: 5.47 kcal
Atomic # 6 Polarization: -4.27 SS G_CDS: -3.64 Total: -7.91 kcal
Atomic # 7 Polarization: -47.57 SS G_CDS: -0.08 Total: -47.65 kcal
Atomic # 8 Polarization: -11.30 SS G_CDS: 0.09 Total: -11.21 kcal
Atomic # 14 Polarization: 20.62 SS G_CDS: -5.03 Total: 15.59 kcal
Total LS contribution 5.94 Total: 5.94 kcal
Total: -28.44 -11.34 -39.78 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.
ZINC001039934872.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 65.933 kcal
(2) G-P(sol) polarization free energy of solvation -28.435 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system 37.497 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.344 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.779 kcal
(6) G-S(sol) free energy of system = (1) + (5) 26.154 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 1.12 seconds