Wall clock time and date at job start Sat Apr 11 2020 06:25:41
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.52995 * 1
3 3 C 1.52997 * 109.47256 * 2 1
4 4 H 1.09000 * 109.46915 * 304.11282 * 3 2 1
5 5 O 1.42901 * 109.47581 * 64.11256 * 3 2 1
6 6 C 1.42895 * 114.00071 * 254.23127 * 5 3 2
7 7 C 1.50704 * 109.47009 * 184.11551 * 3 2 1
8 8 O 1.21275 * 119.99887 * 252.80743 * 7 3 2
9 9 N 1.34776 * 119.99924 * 72.81048 * 7 3 2
10 10 C 1.47016 * 125.64949 * 355.86023 * 9 7 3
11 11 C 1.54325 * 107.27416 * 181.02103 * 10 9 7
12 12 C 1.55160 * 102.93690 * 22.19100 * 11 10 9
13 13 H 1.08997 * 111.03710 * 82.59243 * 12 11 10
14 14 C 1.53002 * 111.00455 * 206.42256 * 12 11 10
15 15 C 1.53004 * 109.47382 * 287.19789 * 14 12 11
16 16 N 1.46900 * 109.47459 * 179.97438 * 15 14 12
17 17 C 1.46897 * 111.00003 * 180.02562 * 16 15 14
18 18 C 1.50699 * 109.47016 * 180.02562 * 17 16 15
19 19 O 1.21296 * 120.00115 * 359.97438 * 18 17 16
20 20 N 1.34783 * 120.00058 * 180.02562 * 18 17 16
21 21 C 1.37094 * 120.00131 * 179.97438 * 20 18 17
22 22 H 1.07999 * 119.99654 * 359.97438 * 21 20 18
23 23 C 1.38942 * 120.00035 * 179.97438 * 21 20 18
24 24 N 1.31413 * 120.00225 * 0.02562 * 23 21 20
25 25 Si 1.86800 * 120.00041 * 179.97438 * 23 21 20
26 26 H 1.48498 * 109.47392 * 0.02562 * 25 23 21
27 27 H 1.48504 * 109.47196 * 120.00036 * 25 23 21
28 28 H 1.48500 * 109.47419 * 240.00442 * 25 23 21
29 29 C 1.47425 * 125.64725 * 175.81447 * 9 7 3
30 30 H 1.09008 * 109.46751 * 59.81863 * 1 2 3
31 31 H 1.09000 * 109.47167 * 179.81525 * 1 2 3
32 32 H 1.09005 * 109.47512 * 299.82052 * 1 2 3
33 33 H 1.08999 * 109.47415 * 239.99477 * 2 1 3
34 34 H 1.09000 * 109.47003 * 119.99328 * 2 1 3
35 35 H 1.09003 * 109.46787 * 304.85659 * 6 5 3
36 36 H 1.08995 * 109.47266 * 64.86176 * 6 5 3
37 37 H 1.08999 * 109.47251 * 184.86632 * 6 5 3
38 38 H 1.08997 * 109.88350 * 300.44380 * 10 9 7
39 39 H 1.09000 * 109.88492 * 61.45211 * 10 9 7
40 40 H 1.09000 * 110.72283 * 263.83361 * 11 10 9
41 41 H 1.08993 * 110.72119 * 140.55209 * 11 10 9
42 42 H 1.08999 * 109.47602 * 167.19688 * 14 12 11
43 43 H 1.09006 * 109.47265 * 47.19833 * 14 12 11
44 44 H 1.08997 * 109.46986 * 300.00470 * 15 14 12
45 45 H 1.09006 * 109.46571 * 59.99730 * 15 14 12
46 46 H 1.00896 * 110.99875 * 303.95628 * 16 15 14
47 47 H 1.09005 * 109.47067 * 300.00378 * 17 16 15
48 48 H 1.08997 * 109.47039 * 60.00066 * 17 16 15
49 49 H 0.96993 * 119.99695 * 0.02562 * 20 18 17
50 50 H 0.96999 * 120.00102 * 179.97438 * 24 23 21
51 51 H 1.09004 * 110.36663 * 37.07140 * 29 9 7
52 52 H 1.09002 * 110.36895 * 274.74879 * 29 9 7
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.5299 0.0000 0.0000
3 6 2.0400 1.4425 0.0000
4 1 1.6177 1.9771 -0.8508
5 8 1.6442 2.0877 1.2121
6 6 0.5779 3.0251 1.0501
7 6 3.5436 1.4436 -0.1020
8 8 4.2167 1.7304 0.8652
9 7 4.1401 1.1260 -1.2680
10 6 3.4529 0.6778 -2.4879
11 6 4.5375 0.4080 -3.5521
12 6 5.8120 0.1764 -2.6981
13 1 5.8575 -0.8497 -2.3332
14 6 7.0734 0.5306 -3.4882
15 6 7.3280 -0.5413 -4.5500
16 7 8.5388 -0.2009 -5.3091
17 6 8.8146 -1.2100 -6.3403
18 6 10.0612 -0.8281 -7.0960
19 8 10.6607 0.1853 -6.8046
20 7 10.5095 -1.6121 -8.0965
21 6 11.6433 -1.2644 -8.7844
22 1 12.1770 -0.3621 -8.5248
23 6 12.1058 -2.0729 -9.8154
24 7 11.4561 -3.1706 -10.1316
25 14 13.6502 -1.5988 -10.7531
26 1 14.1901 -0.3277 -10.2074
27 1 14.6670 -2.6712 -10.6067
28 1 13.3218 -1.4209 -12.1904
29 6 5.5899 1.1712 -1.5317
30 1 -0.3633 0.5167 -0.8884
31 1 -0.3633 -1.0277 -0.0033
32 1 -0.3634 0.5111 0.8916
33 1 1.8933 -0.5139 -0.8899
34 1 1.8933 -0.5137 0.8900
35 1 0.8505 3.7578 0.2906
36 1 -0.3246 2.4986 0.7398
37 1 0.3945 3.5339 1.9965
38 1 2.8966 -0.2373 -2.2852
39 1 2.7730 1.4550 -2.8372
40 1 4.6619 1.2727 -4.2039
41 1 4.2952 -0.4822 -4.1325
42 1 7.9253 0.5796 -2.8100
43 1 6.9390 1.4976 -3.9731
44 1 7.4624 -1.5082 -4.0652
45 1 6.4761 -0.5902 -5.2283
46 1 8.4640 0.7197 -5.7152
47 1 8.9596 -2.1818 -5.8682
48 1 7.9728 -1.2635 -7.0307
49 1 10.0296 -2.4219 -8.3304
50 1 11.7787 -3.7348 -10.8516
51 1 6.1460 0.8524 -0.6500
52 1 5.8905 2.1763 -1.8278
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
ZINC001501541182.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:25:41
Heat of formation + Delta-G solvation = -123.288199 kcal
Electronic energy + Delta-G solvation = -29135.137270 eV
Core-core repulsion = 24941.402923 eV
Total energy + Delta-G solvation = -4193.734347 eV
No. of doubly occupied orbitals = 66
Molecular weight
(most abundant/longest-lived isotopes) = 341.215 amu
Computer time = 0.60 seconds
Orbital eigenvalues (eV)
-40.76202 -38.91655 -38.28180 -37.55804 -35.86802 -34.14440 -33.85413 -32.85869
-31.22625 -30.21597 -29.37115 -28.23755 -26.53476 -25.66122 -24.10451 -23.53984
-22.37112 -22.07200 -20.62024 -19.88198 -19.20421 -17.92224 -17.85874 -17.30280
-16.33314 -16.21887 -16.11714 -15.83864 -15.65657 -14.94115 -14.72377 -14.48416
-14.41694 -14.30491 -13.96324 -13.87316 -13.62240 -13.43696 -13.30788 -13.07612
-12.93379 -12.65061 -12.46054 -12.36626 -12.12559 -12.07727 -11.89728 -11.58704
-11.53242 -11.44712 -11.33306 -11.07832 -11.06154 -10.96148 -10.72754 -10.55587
-10.34593 -10.22837 -9.74952 -9.52505 -9.02550 -8.46629 -8.05953 -7.96213
-6.42365 -3.83121 2.22474 2.33034 3.16041 3.26696 3.34095 3.37856
3.88663 3.90721 4.03418 4.16035 4.23170 4.27624 4.33115 4.35802
4.50721 4.65871 4.76497 4.80043 4.86070 4.93065 4.95514 5.00331
5.09608 5.16408 5.19819 5.22039 5.23860 5.26319 5.36736 5.41833
5.44233 5.49453 5.52862 5.55180 5.63381 5.69665 5.75027 5.78800
5.89790 5.97398 6.09091 6.13110 6.27117 6.53842 7.09312 7.21972
7.55006 7.66208 7.80450 8.06259 8.60891 9.17613 9.53954 10.05432
10.75474
Molecular weight = 341.21amu
Principal moments of inertia in cm(-1)
A = 0.037529 B = 0.001388 C = 0.001359
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 745.910386 B =20169.861035 C =20597.592404
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C -0.152 4.152
2 C -0.138 4.138
3 C 0.078 3.922
4 H 0.086 0.914
5 O -0.357 6.357
6 C 0.025 3.975
7 C 0.522 3.478
8 O -0.518 6.518
9 N -0.619 5.619
10 C 0.094 3.906
11 C -0.129 4.129
12 C -0.097 4.097
13 H 0.088 0.912
14 C -0.107 4.107
15 C 0.061 3.939
16 N -0.671 5.671
17 C 0.051 3.949
18 C 0.440 3.560
19 O -0.583 6.583
20 N -0.564 5.564
21 C -0.299 4.299
22 H 0.104 0.896
23 C 0.022 3.978
24 N -0.902 5.902
25 Si 0.929 3.071
26 H -0.268 1.268
27 H -0.281 1.281
28 H -0.282 1.282
29 C 0.110 3.890
30 H 0.054 0.946
31 H 0.062 0.938
32 H 0.062 0.938
33 H 0.080 0.920
34 H 0.079 0.921
35 H 0.043 0.957
36 H 0.047 0.953
37 H 0.091 0.909
38 H 0.072 0.928
39 H 0.068 0.932
40 H 0.080 0.920
41 H 0.085 0.915
42 H 0.079 0.921
43 H 0.066 0.934
44 H 0.070 0.930
45 H 0.026 0.974
46 H 0.349 0.651
47 H 0.085 0.915
48 H 0.042 0.958
49 H 0.392 0.608
50 H 0.271 0.729
51 H 0.081 0.919
52 H 0.072 0.928
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges -26.824 5.333 16.304 31.840
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.210 4.210
2 C -0.176 4.176
3 C 0.017 3.983
4 H 0.103 0.897
5 O -0.274 6.274
6 C -0.071 4.071
7 C 0.309 3.691
8 O -0.393 6.393
9 N -0.356 5.356
10 C -0.029 4.029
11 C -0.167 4.167
12 C -0.117 4.117
13 H 0.106 0.894
14 C -0.145 4.145
15 C -0.068 4.068
16 N -0.303 5.303
17 C -0.081 4.081
18 C 0.226 3.774
19 O -0.467 6.467
20 N -0.201 5.201
21 C -0.429 4.429
22 H 0.121 0.879
23 C -0.179 4.179
24 N -0.544 5.544
25 Si 0.756 3.244
26 H -0.192 1.192
27 H -0.207 1.207
28 H -0.208 1.208
29 C -0.012 4.012
30 H 0.074 0.926
31 H 0.081 0.919
32 H 0.081 0.919
33 H 0.099 0.901
34 H 0.097 0.903
35 H 0.062 0.938
36 H 0.065 0.935
37 H 0.109 0.891
38 H 0.090 0.910
39 H 0.086 0.914
40 H 0.099 0.901
41 H 0.104 0.896
42 H 0.097 0.903
43 H 0.085 0.915
44 H 0.088 0.912
45 H 0.044 0.956
46 H 0.168 0.832
47 H 0.103 0.897
48 H 0.060 0.940
49 H 0.225 0.775
50 H 0.081 0.919
51 H 0.099 0.901
52 H 0.090 0.910
Dipole moment (debyes) X Y Z Total
from point charges -25.544 5.753 17.612 31.557
hybrid contribution -0.457 -0.147 -1.701 1.767
sum -26.002 5.607 15.912 30.995
Atomic orbital electron populations
1.21805 0.94534 1.02534 1.02088 1.21963 0.97256 0.94521 1.03889
1.21992 0.91413 0.95620 0.89290 0.89655 1.88039 1.53319 1.57228
1.28856 1.23009 0.89440 0.91775 1.02926 1.20353 0.89995 0.77189
0.81524 1.90735 1.65688 1.47097 1.35742 1.48479 1.08418 1.67307
1.11376 1.22221 0.96062 0.99677 0.84891 1.22769 0.93881 1.03177
0.96840 1.22370 0.94812 0.99800 0.94723 0.89369 1.21480 0.96422
0.98925 0.97657 1.21945 0.90797 0.98842 0.95260 1.60051 1.36277
1.05755 1.28218 1.21660 0.96703 0.96546 0.93226 1.19837 0.87751
0.86346 0.83493 1.90525 1.60056 1.27335 1.68799 1.41872 1.25207
1.23283 1.29746 1.19762 1.02140 1.08798 1.12160 0.87858 1.25257
0.99495 0.95494 0.97696 1.69620 1.40276 1.13583 1.30883 0.85801
0.80904 0.78718 0.78933 1.19202 1.20739 1.20785 1.22325 0.80484
0.99761 0.98599 0.92649 0.91875 0.91931 0.90119 0.90260 0.93846
0.93461 0.89098 0.90955 0.91376 0.90091 0.89590 0.90265 0.91505
0.91160 0.95558 0.83222 0.89653 0.93960 0.77476 0.91940 0.90055
0.91022
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.15 -0.67 9.17 37.15 0.34 -0.33 16
2 C -0.14 -0.84 5.56 -26.74 -0.15 -0.99 16
3 C 0.08 0.62 2.07 -27.88 -0.06 0.56 16
4 H 0.09 0.47 7.33 -51.93 -0.38 0.09 16
5 O -0.36 -3.46 9.56 -54.81 -0.52 -3.98 16
6 C 0.02 0.14 10.43 101.05 1.05 1.20 16
7 C 0.52 5.72 7.17 -10.98 -0.08 5.64 16
8 O -0.52 -8.08 16.78 -14.01 -0.24 -8.31 16
9 N -0.62 -5.34 3.33 -170.53 -0.57 -5.91 16
10 C 0.09 0.51 6.11 -3.07 -0.02 0.49 16
11 C -0.13 -0.72 5.90 -24.89 -0.15 -0.87 16
12 C -0.10 -0.71 3.00 -88.47 -0.27 -0.97 16
13 H 0.09 0.65 8.14 -51.93 -0.42 0.22 16
14 C -0.11 -0.96 4.98 -26.73 -0.13 -1.10 16
15 C 0.06 0.67 5.59 -3.82 -0.02 0.65 16
16 N -0.67 -10.43 6.13 -115.06 -0.70 -11.14 16
17 C 0.05 0.89 6.17 -4.98 -0.03 0.86 16
18 C 0.44 10.36 7.82 -10.98 -0.09 10.27 16
19 O -0.58 -15.85 15.78 5.54 0.09 -15.77 16
20 N -0.56 -14.26 5.29 -10.96 -0.06 -14.32 16
21 C -0.30 -8.34 9.68 -14.94 -0.14 -8.49 16
22 H 0.10 3.03 7.16 -52.49 -0.38 2.66 16
23 C 0.02 0.59 5.50 -17.47 -0.10 0.49 16
24 N -0.90 -24.43 13.05 55.50 0.72 -23.71 16
25 Si 0.93 20.99 29.55 -169.99 -5.02 15.96 16
26 H -0.27 -6.06 7.11 56.52 0.40 -5.66 16
27 H -0.28 -6.26 7.11 56.52 0.40 -5.86 16
28 H -0.28 -6.25 7.11 56.52 0.40 -5.85 16
29 C 0.11 0.97 6.40 -2.53 -0.02 0.96 16
30 H 0.05 0.18 8.14 -51.92 -0.42 -0.24 16
31 H 0.06 0.24 8.14 -51.93 -0.42 -0.18 16
32 H 0.06 0.31 6.51 -51.93 -0.34 -0.03 16
33 H 0.08 0.38 5.54 -51.93 -0.29 0.10 16
34 H 0.08 0.60 8.14 -51.93 -0.42 0.18 16
35 H 0.04 0.21 8.04 -51.93 -0.42 -0.20 16
36 H 0.05 0.18 6.43 -51.93 -0.33 -0.15 16
37 H 0.09 0.43 8.14 -51.93 -0.42 0.01 16
38 H 0.07 0.30 5.95 -51.93 -0.31 -0.01 16
39 H 0.07 0.28 8.01 -51.93 -0.42 -0.13 16
40 H 0.08 0.46 7.91 -51.93 -0.41 0.04 16
41 H 0.09 0.41 8.14 -51.93 -0.42 -0.01 16
42 H 0.08 0.81 8.14 -51.93 -0.42 0.39 16
43 H 0.07 0.59 7.97 -51.93 -0.41 0.18 16
44 H 0.07 0.76 8.14 -51.93 -0.42 0.34 16
45 H 0.03 0.26 7.75 -51.93 -0.40 -0.14 16
46 H 0.35 5.57 8.43 -39.29 -0.33 5.24 16
47 H 0.09 1.37 8.14 -51.93 -0.42 0.95 16
48 H 0.04 0.65 8.14 -51.93 -0.42 0.23 16
49 H 0.39 9.43 7.90 -40.82 -0.32 9.11 16
50 H 0.27 6.92 8.31 -40.82 -0.34 6.58 16
51 H 0.08 0.83 8.00 -51.93 -0.42 0.42 16
52 H 0.07 0.64 8.04 -51.93 -0.42 0.22 16
LS Contribution 417.04 15.07 6.28 6.28
Total: -1.00 -35.24 417.04 -8.80 -44.04
By element:
Atomic # 1 Polarization: 17.40 SS G_CDS: -8.94 Total: 8.47 kcal
Atomic # 6 Polarization: 8.22 SS G_CDS: 0.15 Total: 8.37 kcal
Atomic # 7 Polarization: -54.46 SS G_CDS: -0.61 Total: -55.07 kcal
Atomic # 8 Polarization: -27.39 SS G_CDS: -0.67 Total: -28.06 kcal
Atomic # 14 Polarization: 20.99 SS G_CDS: -5.02 Total: 15.96 kcal
Total LS contribution 6.28 Total: 6.28 kcal
Total: -35.24 -8.80 -44.04 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.
ZINC001501541182.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 -79.246 kcal
(2) G-P(sol) polarization free energy of solvation -35.242 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system -114.487 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.801 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.043 kcal
(6) G-S(sol) free energy of system = (1) + (5) -123.288 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 0.60 seconds