Wall clock time and date at job start Wed Apr 1 2020 05:21:49
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 N 1.46506 * 1
3 3 C 1.34944 * 126.04069 * 2 1
4 4 C 1.35309 * 107.74803 * 180.14789 * 3 2 1
5 5 C 1.50698 * 126.04044 * 179.87155 * 4 3 2
6 6 C 1.52993 * 109.47189 * 52.69573 * 5 4 3
7 7 C 1.53007 * 109.47167 * 292.69352 * 5 4 3
8 8 C 1.50703 * 109.47283 * 172.69722 * 5 4 3
9 9 O 1.21277 * 120.00085 * 124.74726 * 8 5 4
10 10 N 1.34778 * 119.99833 * 304.74373 * 8 5 4
11 11 C 1.46803 * 119.57634 * 356.06682 * 10 8 5
12 12 C 1.53461 * 107.60200 * 127.97452 * 11 10 8
13 13 C 1.52733 * 108.43963 * 58.06930 * 12 11 10
14 14 C 1.54524 * 110.41434 * 295.69237 * 13 12 11
15 15 H 1.08997 * 113.49166 * 291.36050 * 14 13 12
16 16 C 1.45117 * 119.57708 * 176.08213 * 10 8 5
17 17 H 1.09001 * 109.97151 * 348.96342 * 16 10 8
18 18 C 1.54425 * 111.60619 * 111.32954 * 16 10 8
19 19 C 1.55464 * 102.03261 * 155.95219 * 18 16 10
20 20 N 1.47120 * 107.60199 * 339.51701 * 19 18 16
21 21 C 1.36937 * 125.76677 * 176.63333 * 20 19 18
22 22 H 1.08007 * 119.99965 * 180.02562 * 21 20 19
23 23 C 1.38811 * 120.00077 * 0.04466 * 21 20 19
24 24 N 1.31614 * 120.00064 * 359.97438 * 23 21 20
25 25 Si 1.86804 * 120.00023 * 179.97438 * 23 21 20
26 26 H 1.48500 * 109.46900 * 59.99989 * 25 23 21
27 27 H 1.48498 * 109.47105 * 179.97438 * 25 23 21
28 28 H 1.48490 * 109.47131 * 299.99760 * 25 23 21
29 29 C 1.39868 * 107.91881 * 359.60350 * 4 3 2
30 30 N 1.30838 * 108.20241 * 0.24594 * 29 4 3
31 31 H 1.08999 * 109.46929 * 270.26818 * 1 2 3
32 32 H 1.08996 * 109.46718 * 30.26607 * 1 2 3
33 33 H 1.08995 * 109.47055 * 150.26805 * 1 2 3
34 34 H 1.08002 * 126.12584 * 359.97438 * 3 2 1
35 35 H 1.09001 * 109.47627 * 184.39169 * 6 5 4
36 36 H 1.09006 * 109.47646 * 304.39050 * 6 5 4
37 37 H 1.09001 * 109.47689 * 64.39522 * 6 5 4
38 38 H 1.08996 * 109.46689 * 292.58532 * 7 5 4
39 39 H 1.08998 * 109.46593 * 52.58249 * 7 5 4
40 40 H 1.08999 * 109.47104 * 172.58024 * 7 5 4
41 41 H 1.08997 * 109.83486 * 8.44440 * 11 10 8
42 42 H 1.09000 * 109.82965 * 247.49707 * 11 10 8
43 43 H 1.08991 * 109.64575 * 177.75937 * 12 11 10
44 44 H 1.09007 * 109.63635 * 298.39008 * 12 11 10
45 45 H 1.08998 * 109.21097 * 56.57922 * 13 12 11
46 46 H 1.09007 * 109.20747 * 175.89123 * 13 12 11
47 47 H 1.08997 * 110.80815 * 274.81671 * 18 16 10
48 48 H 1.08997 * 110.81159 * 38.33668 * 18 16 10
49 49 H 1.09002 * 110.31946 * 219.64936 * 19 18 16
50 50 H 1.08997 * 109.23997 * 98.77227 * 19 18 16
51 51 H 0.96999 * 120.00205 * 180.02562 * 24 23 21
52 52 H 1.08004 * 125.89781 * 180.22520 * 29 4 3
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 7 1.4651 0.0000 0.0000
3 6 2.2590 1.0912 0.0000
4 6 3.5437 0.6665 -0.0033
5 6 4.7680 1.5451 -0.0073
6 6 4.6758 2.5534 1.1397
7 6 4.8542 2.2939 -1.3389
8 6 5.9984 0.6932 0.1704
9 8 6.8981 0.7507 -0.6408
10 7 6.0989 -0.1324 1.2309
11 6 5.0447 -0.1313 2.2526
12 6 4.5895 -1.5833 2.4512
13 6 5.8009 -2.4175 2.8625
14 6 6.8441 -2.4512 1.7231
15 1 6.5149 -3.0130 0.8490
16 6 7.2480 -1.0109 1.3484
17 1 7.8122 -1.0126 0.4158
18 6 8.1646 -0.5947 2.5194
19 6 8.8927 -1.9277 2.8508
20 7 8.0799 -3.0260 2.3054
21 6 8.4097 -4.3548 2.3307
22 1 7.7408 -5.0856 1.9005
23 6 9.6036 -4.7627 2.9097
24 7 10.4188 -3.8722 3.4338
25 14 10.0532 -6.5755 2.9449
26 1 10.1333 -7.0927 1.5552
27 1 11.3685 -6.7464 3.6127
28 1 9.0156 -7.3283 3.6943
29 6 3.5346 -0.7322 0.0038
30 7 2.2891 -1.1328 0.0058
31 1 -0.3633 0.0048 1.0277
32 1 -0.3632 0.8876 -0.5180
33 1 -0.3633 -0.8923 -0.5096
34 1 1.9283 2.1193 -0.0004
35 1 5.5221 3.2387 1.0901
36 1 3.7466 3.1169 1.0545
37 1 4.6932 2.0232 2.0919
38 1 5.0518 1.5844 -2.1424
39 1 3.9108 2.8054 -1.5296
40 1 5.6612 3.0251 -1.2930
41 1 4.2049 0.4770 1.9168
42 1 5.4378 0.2644 3.1891
43 1 3.8323 -1.6261 3.2340
44 1 4.1786 -1.9692 1.5182
45 1 6.2459 -1.9859 3.7590
46 1 5.4784 -3.4362 3.0779
47 1 8.8668 0.1778 2.2059
48 1 7.5748 -0.2601 3.3729
49 1 8.9971 -2.0440 3.9295
50 1 9.8764 -1.9260 2.3813
51 1 11.2532 -4.1571 3.8380
52 1 4.4054 -1.3711 0.0067
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
ZINC001049829313.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:21:49
Heat of formation + Delta-G solvation = 99.077131 kcal
Electronic energy + Delta-G solvation = -33283.396973 eV
Core-core repulsion = 29286.565162 eV
Total energy + Delta-G solvation = -3996.831811 eV
No. of doubly occupied orbitals = 66
Molecular weight
(most abundant/longest-lived isotopes) = 346.211 amu
Computer time = 1.20 seconds
Orbital eigenvalues (eV)
-41.98755 -40.16709 -38.13211 -35.77055 -34.90829 -33.16161 -32.21620 -32.06531
-29.57281 -28.97089 -28.27776 -27.31064 -26.82341 -25.12720 -24.20835 -23.05293
-22.01756 -21.38169 -20.64110 -19.97557 -19.07200 -18.33416 -17.40809 -16.80833
-16.28557 -15.74674 -15.48220 -15.03635 -14.72046 -14.54606 -14.25655 -14.19484
-14.03117 -13.70558 -13.38426 -13.26080 -13.00800 -12.74685 -12.66034 -12.58128
-12.40324 -12.29490 -11.94487 -11.89066 -11.78572 -11.73380 -11.32155 -11.27230
-10.91147 -10.77352 -10.69566 -10.62748 -10.49293 -10.35971 -10.06626 -9.98470
-9.60792 -9.53315 -9.40944 -9.30191 -8.57089 -8.45191 -8.16165 -6.80385
-5.66290 -3.00140 1.70236 2.32363 2.62306 2.86907 3.26624 3.44420
3.48682 3.49697 4.25595 4.42340 4.53390 4.63807 4.75779 4.84251
4.94780 4.99574 5.07886 5.12073 5.19397 5.26981 5.37363 5.45487
5.48258 5.51387 5.57392 5.63534 5.66618 5.69021 5.72359 5.78857
5.86804 5.93137 5.97217 6.04476 6.14777 6.20445 6.22654 6.27760
6.39236 6.44121 6.47810 6.58395 6.68948 6.79375 6.82551 7.07095
7.13312 7.54300 7.60145 7.66120 7.73904 7.87653 8.28491 8.46760
9.25469 9.85196 10.51507 11.44741
Molecular weight = 346.21amu
Principal moments of inertia in cm(-1)
A = 0.016651 B = 0.003408 C = 0.003062
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 1681.144097 B = 8214.428279 C = 9141.982205
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C 0.103 3.897
2 N -0.379 5.379
3 C 0.063 3.937
4 C -0.227 4.227
5 C 0.000 4.000
6 C -0.149 4.149
7 C -0.121 4.121
8 C 0.521 3.479
9 O -0.539 6.539
10 N -0.588 5.588
11 C 0.093 3.907
12 C -0.145 4.145
13 C -0.126 4.126
14 C 0.156 3.844
15 H 0.062 0.938
16 C 0.128 3.872
17 H 0.099 0.901
18 C -0.151 4.151
19 C 0.182 3.818
20 N -0.598 5.598
21 C -0.225 4.225
22 H 0.077 0.923
23 C -0.007 4.007
24 N -0.922 5.922
25 Si 0.911 3.089
26 H -0.285 1.285
27 H -0.291 1.291
28 H -0.286 1.286
29 C 0.064 3.936
30 N -0.291 5.291
31 H 0.073 0.927
32 H 0.085 0.915
33 H 0.088 0.912
34 H 0.173 0.827
35 H 0.064 0.936
36 H 0.064 0.936
37 H 0.072 0.928
38 H 0.063 0.937
39 H 0.056 0.944
40 H 0.062 0.938
41 H 0.101 0.899
42 H 0.064 0.936
43 H 0.062 0.938
44 H 0.077 0.923
45 H 0.062 0.938
46 H 0.057 0.943
47 H 0.061 0.939
48 H 0.057 0.943
49 H 0.041 0.959
50 H 0.049 0.951
51 H 0.252 0.748
52 H 0.188 0.812
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges -20.664 15.335 -4.417 26.109
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.038 4.038
2 N -0.161 5.161
3 C -0.083 4.083
4 C -0.239 4.239
5 C -0.005 4.005
6 C -0.207 4.207
7 C -0.178 4.178
8 C 0.309 3.691
9 O -0.416 6.416
10 N -0.319 5.319
11 C -0.029 4.029
12 C -0.184 4.184
13 C -0.166 4.166
14 C 0.047 3.953
15 H 0.080 0.920
16 C 0.023 3.977
17 H 0.117 0.883
18 C -0.191 4.191
19 C 0.056 3.944
20 N -0.322 5.322
21 C -0.354 4.354
22 H 0.095 0.905
23 C -0.203 4.203
24 N -0.569 5.569
25 Si 0.742 3.258
26 H -0.212 1.212
27 H -0.217 1.217
28 H -0.212 1.212
29 C -0.122 4.122
30 N -0.112 5.112
31 H 0.092 0.908
32 H 0.103 0.897
33 H 0.107 0.893
34 H 0.190 0.810
35 H 0.083 0.917
36 H 0.083 0.917
37 H 0.090 0.910
38 H 0.082 0.918
39 H 0.075 0.925
40 H 0.081 0.919
41 H 0.119 0.881
42 H 0.082 0.918
43 H 0.081 0.919
44 H 0.095 0.905
45 H 0.081 0.919
46 H 0.076 0.924
47 H 0.080 0.920
48 H 0.076 0.924
49 H 0.059 0.941
50 H 0.068 0.932
51 H 0.063 0.937
52 H 0.205 0.795
Dipole moment (debyes) X Y Z Total
from point charges -20.921 14.185 -4.749 25.719
hybrid contribution 0.777 -0.305 0.325 0.896
sum -20.144 13.880 -4.424 24.860
Atomic orbital electron populations
1.21823 0.74314 1.04365 1.03304 1.48451 1.06546 1.03847 1.57210
1.22096 0.90427 0.93328 1.02494 1.19632 0.93394 0.93652 1.17256
1.19040 0.89638 0.93785 0.98035 1.22165 1.02571 0.98330 0.97591
1.21659 1.01299 0.99900 0.94977 1.20394 0.87843 0.80023 0.80882
1.90682 1.40852 1.66568 1.43502 1.47597 1.22736 1.39005 1.22578
1.22300 0.90709 0.98383 0.91511 1.22267 0.98689 0.95114 1.02299
1.22349 0.94836 0.99326 1.00112 1.21694 0.88571 0.90968 0.94039
0.91973 1.21817 0.85751 0.89724 1.00367 0.88277 1.23241 0.98579
1.00087 0.97146 1.21301 0.93250 0.81890 0.97911 1.44691 1.17522
1.02869 1.67137 1.18991 1.02442 0.81489 1.32457 0.90497 1.24782
0.97051 0.97698 1.00805 1.69947 1.11894 1.31428 1.43628 0.86213
0.77404 0.84620 0.77601 1.21187 1.21677 1.21250 1.24117 0.95532
0.94532 0.98017 1.77477 1.01366 1.07739 1.24621 0.90797 0.89669
0.89322 0.81019 0.91714 0.91696 0.90963 0.91789 0.92475 0.91875
0.88133 0.91787 0.91948 0.90516 0.91940 0.92390 0.92007 0.92417
0.94094 0.93231 0.93745 0.79519
Number of geometries 1
Number of calculations of the screened coulomb radii 1
The total number of SCF iterations 14.
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.10 0.52 11.47 59.85 0.69 1.21 16
2 N -0.38 -3.53 4.01 -68.16 -0.27 -3.80 16
3 C 0.06 0.52 10.62 -17.23 -0.18 0.34 16
4 C -0.23 -2.46 2.90 -105.06 -0.30 -2.76 16
5 C 0.00 -0.01 0.67 -156.81 -0.11 -0.11 16
6 C -0.15 -1.32 8.02 37.15 0.30 -1.03 16
7 C -0.12 -1.28 7.82 37.16 0.29 -0.99 16
8 C 0.52 8.05 5.73 -10.98 -0.06 7.99 16
9 O -0.54 -10.36 15.26 5.56 0.08 -10.27 16
10 N -0.59 -9.03 1.94 -165.02 -0.32 -9.35 16
11 C 0.09 1.21 3.56 -3.64 -0.01 1.19 16
12 C -0.14 -2.04 5.46 -26.63 -0.15 -2.18 16
13 C -0.13 -2.26 5.11 -26.07 -0.13 -2.39 16
14 C 0.16 3.19 3.74 -65.56 -0.25 2.95 16
15 H 0.06 1.30 8.14 -51.93 -0.42 0.88 16
16 C 0.13 2.42 3.03 -67.29 -0.20 2.21 16
17 H 0.10 2.13 6.98 -51.93 -0.36 1.77 16
18 C -0.15 -2.96 6.16 -24.67 -0.15 -3.11 16
19 C 0.18 4.44 5.17 -2.39 -0.01 4.43 16
20 N -0.60 -15.08 2.89 -160.74 -0.46 -15.54 16
21 C -0.23 -6.26 10.29 -15.06 -0.16 -6.42 16
22 H 0.08 2.07 7.83 -52.48 -0.41 1.66 16
23 C -0.01 -0.20 5.49 -17.43 -0.10 -0.30 16
24 N -0.92 -27.22 10.33 55.49 0.57 -26.64 16
25 Si 0.91 22.36 29.55 -169.99 -5.02 17.34 16
26 H -0.29 -6.94 7.11 56.52 0.40 -6.53 16
27 H -0.29 -7.21 7.11 56.52 0.40 -6.81 16
28 H -0.29 -6.91 7.11 56.52 0.40 -6.51 16
29 C 0.06 0.83 7.07 -17.37 -0.12 0.71 16
30 N -0.29 -3.73 12.58 30.63 0.39 -3.34 16
31 H 0.07 0.30 8.14 -51.93 -0.42 -0.13 16
32 H 0.08 0.21 8.14 -51.93 -0.42 -0.22 16
33 H 0.09 0.45 8.14 -51.93 -0.42 0.03 16
34 H 0.17 0.83 7.88 -52.49 -0.41 0.41 16
35 H 0.06 0.59 8.14 -51.93 -0.42 0.17 16
36 H 0.06 0.41 7.60 -51.93 -0.39 0.02 16
37 H 0.07 0.66 5.09 -51.93 -0.26 0.40 16
38 H 0.06 0.80 7.94 -51.93 -0.41 0.39 16
39 H 0.06 0.45 7.97 -51.93 -0.41 0.04 16
40 H 0.06 0.68 8.14 -51.93 -0.42 0.26 16
41 H 0.10 1.13 3.01 -51.93 -0.16 0.97 16
42 H 0.06 0.82 7.51 -51.93 -0.39 0.43 16
43 H 0.06 0.80 8.14 -51.93 -0.42 0.38 16
44 H 0.08 1.11 4.51 -51.93 -0.23 0.87 16
45 H 0.06 1.20 7.10 -51.93 -0.37 0.83 16
46 H 0.06 1.10 8.14 -51.92 -0.42 0.68 16
47 H 0.06 1.18 8.14 -51.93 -0.42 0.75 16
48 H 0.06 1.05 6.80 -51.93 -0.35 0.69 16
49 H 0.04 1.10 7.60 -51.93 -0.39 0.70 16
50 H 0.05 1.36 7.29 -51.93 -0.38 0.98 16
51 H 0.25 7.28 8.31 -40.82 -0.34 6.95 16
52 H 0.19 2.75 5.48 -52.48 -0.29 2.46 16
LS Contribution 382.35 15.07 5.76 5.76
Total: -1.00 -33.47 382.35 -8.11 -41.58
By element:
Atomic # 1 Polarization: 10.71 SS G_CDS: -8.17 Total: 2.53 kcal
Atomic # 6 Polarization: 2.40 SS G_CDS: -0.66 Total: 1.74 kcal
Atomic # 7 Polarization: -58.58 SS G_CDS: -0.10 Total: -58.68 kcal
Atomic # 8 Polarization: -10.36 SS G_CDS: 0.08 Total: -10.27 kcal
Atomic # 14 Polarization: 22.36 SS G_CDS: -5.02 Total: 17.34 kcal
Total LS contribution 5.76 Total: 5.76 kcal
Total: -33.47 -8.11 -41.58 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.
ZINC001049829313.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 140.658 kcal
(2) G-P(sol) polarization free energy of solvation -33.469 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system 107.188 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.111 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.580 kcal
(6) G-S(sol) free energy of system = (1) + (5) 99.077 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 1.20 seconds