Wall clock time and date at job start Mon Mar 8 2021 04:40:21 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= 0 * AM1 - THE AM1 HAMILTONIAN TO BE USED * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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.46901 * 1 3 3 C 1.46906 * 110.99793 * 2 1 4 4 C 1.50696 * 109.47351 * 293.95574 * 3 2 1 5 5 O 1.21285 * 119.99962 * 359.97438 * 4 3 2 6 6 N 1.34775 * 120.00393 * 180.02562 * 4 3 2 7 7 C 1.46497 * 120.00052 * 180.02562 * 6 4 3 8 8 H 1.08994 * 112.85368 * 23.80473 * 7 6 4 9 9 C 1.53790 * 113.61321 * 155.00233 * 7 6 4 10 10 C 1.53775 * 87.08086 * 89.11176 * 9 7 6 11 11 H 1.08996 * 113.61999 * 270.87940 * 10 9 7 12 12 C 1.50707 * 113.61179 * 140.06207 * 10 9 7 13 13 C 1.38200 * 119.98490 * 299.73046 * 12 10 9 14 14 C 1.38342 * 120.00315 * 179.72640 * 13 12 10 15 15 Cl 1.73596 * 120.02374 * 180.27033 * 14 13 12 16 16 C 1.38448 * 119.95989 * 0.54398 * 14 13 12 17 17 Cl 1.73597 * 120.02127 * 179.75212 * 16 14 13 18 18 C 1.38355 * 119.96089 * 359.72596 * 16 14 13 19 19 C 1.38202 * 119.99718 * 0.02562 * 18 16 14 20 20 C 1.53780 * 87.08364 * 25.43301 * 10 9 7 21 21 C 1.39673 * 111.00219 * 123.95459 * 2 1 3 22 22 C 1.38890 * 120.07192 * 115.76903 * 21 2 1 23 23 C 1.50699 * 120.03969 * 0.02562 * 22 21 2 24 24 O 1.42895 * 109.46882 * 284.99754 * 23 22 21 25 25 Si 1.86300 * 109.46948 * 164.99833 * 23 22 21 26 26 C 1.38129 * 119.92191 * 179.78839 * 22 21 2 27 27 C 1.38262 * 120.07326 * 0.49339 * 26 22 21 28 28 C 1.50701 * 119.92475 * 179.72367 * 27 26 22 29 29 F 1.39903 * 109.47073 * 180.02562 * 28 27 26 30 30 F 1.39899 * 109.46925 * 300.01416 * 28 27 26 31 31 F 1.39892 * 109.47561 * 60.01626 * 28 27 26 32 32 C 1.38313 * 120.15150 * 359.73818 * 27 26 22 33 33 C 1.38077 * 120.07164 * 0.02562 * 32 27 26 34 34 H 1.09002 * 109.47401 * 59.99808 * 1 2 3 35 35 H 1.09001 * 109.46874 * 299.99821 * 1 2 3 36 36 H 1.09004 * 109.46970 * 179.97438 * 1 2 3 37 37 H 1.08997 * 109.46995 * 53.95554 * 3 2 1 38 38 H 1.09001 * 109.46848 * 173.95560 * 3 2 1 39 39 H 0.96996 * 119.99654 * 0.02562 * 6 4 3 40 40 H 1.08996 * 113.60920 * 203.66542 * 9 7 6 41 41 H 1.08998 * 113.61162 * 334.56091 * 9 7 6 42 42 H 1.08009 * 119.99859 * 359.97438 * 13 12 10 43 43 H 1.07990 * 119.99802 * 179.97438 * 18 16 14 44 44 H 1.08003 * 119.98285 * 179.97438 * 19 18 16 45 45 H 1.09000 * 113.61449 * 220.01651 * 20 10 9 46 46 H 1.09003 * 113.61883 * 89.11282 * 20 10 9 47 47 H 1.09001 * 109.47365 * 45.00015 * 23 22 21 48 48 H 0.96697 * 114.00408 * 300.00097 * 24 23 22 49 49 H 1.48507 * 109.99784 * 179.97438 * 25 23 22 50 50 H 1.48503 * 110.00100 * 301.31695 * 25 23 22 51 51 H 1.07999 * 119.96080 * 180.23130 * 26 22 21 52 52 H 1.08000 * 119.96216 * 180.02562 * 32 27 26 53 53 H 1.07997 * 120.03526 * 179.97438 * 33 32 27 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.4690 0.0000 0.0000 3 6 1.9954 1.3715 0.0000 4 6 1.6369 2.0472 -1.2984 5 8 1.0054 1.4454 -2.1410 6 7 2.0188 3.3200 -1.5237 7 6 1.6697 3.9770 -2.7857 8 1 0.7960 3.5356 -3.2652 9 6 1.6087 5.5104 -2.6856 10 6 3.1060 5.5143 -3.0360 11 1 3.7583 5.4458 -2.1655 12 6 3.5024 6.6213 -3.9786 13 6 2.9335 6.6946 -5.2360 14 6 3.2929 7.7139 -6.0995 15 17 2.5734 7.8106 -7.6764 16 6 4.2311 8.6536 -5.7075 17 17 4.6880 9.9296 -6.7923 18 6 4.8038 8.5755 -4.4505 19 6 4.4397 7.5595 -3.5873 20 6 2.8659 4.1651 -3.7336 21 6 1.9696 -0.7283 1.0816 22 6 2.7069 -1.8844 0.8606 23 6 2.9685 -2.3548 -0.5470 24 8 1.7783 -2.9350 -1.0842 25 14 4.3286 -3.6279 -0.5335 26 6 3.2042 -2.5999 1.9324 27 6 2.9601 -2.1719 3.2243 28 6 3.4990 -2.9585 4.3913 29 9 3.1251 -2.3348 5.5865 30 9 2.9773 -4.2562 4.3596 31 9 4.8945 -3.0148 4.3112 32 6 2.2222 -1.0236 3.4480 33 6 1.7263 -0.3011 2.3810 34 1 -0.3634 0.5139 -0.8900 35 1 -0.3633 0.5138 0.8900 36 1 -0.3633 -1.0277 -0.0005 37 1 1.5611 1.9274 0.8309 38 1 3.0797 1.3445 0.1082 39 1 2.5247 3.8009 -0.8502 40 1 0.9808 5.9729 -3.4470 41 1 1.3944 5.8741 -1.6806 42 1 2.2052 5.9585 -5.5429 43 1 5.5359 9.3079 -4.1443 44 1 4.8872 7.4985 -2.6062 45 1 3.6666 3.4428 -3.5749 46 1 2.5817 4.2598 -4.7817 47 1 3.2699 -1.5077 -1.1631 48 1 1.4517 -3.6961 -0.5851 49 1 4.5957 -4.1003 -1.9159 50 1 5.5629 -3.0463 0.0524 51 1 3.7807 -3.4969 1.7612 52 1 2.0337 -0.6919 4.4583 53 1 1.1498 0.5952 2.5563 There are 83 doubly filled levels RHF calculation, no. of doubly occupied orbitals= 83 No. of singly occupied orbitals= 1 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) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). Si: (AM1): M.J.S.DEWAR, C. JIE, ORGANOMETALLICS, 6, 1486-1490 (1987). Cl: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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=0 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=WATER s_22_9750034_6406254.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Mon Mar 8 2021 04:40:21 Heat of formation + Delta-G solvation = -119.345866 kcal Electronic energy + Delta-G solvation = -44546.417513 eV Core-core repulsion = 38216.145015 eV Total energy + Delta-G solvation = -6330.272497 eV No. of doubly occupied orbitals = 83 and no. of open levels = 1 Molecular weight (most abundant/longest-lived isotopes) = 489.087 amu Computer time = 1.98 seconds Orbital eigenvalues (eV) -51.16775 -48.69549 -48.62748 -41.85547 -41.42464 -40.45132 -39.18606 -37.70036 -37.59481 -37.56069 -37.33893 -34.56004 -32.96495 -32.16072 -31.97136 -31.70971 -29.35659 -27.58068 -26.96690 -26.72071 -26.22593 -25.02649 -24.09639 -23.37594 -22.72766 -22.13108 -21.54907 -20.21507 -19.58895 -19.48645 -19.38233 -18.74872 -18.35593 -17.71997 -17.49560 -17.22925 -16.99657 -16.87246 -16.59268 -16.39797 -16.26339 -15.87276 -15.71118 -15.53446 -15.51848 -15.35488 -15.23811 -15.15797 -15.06841 -15.00997 -14.84168 -14.78685 -14.77575 -14.61156 -14.47598 -14.22502 -13.93749 -13.80351 -13.65038 -13.58680 -13.45356 -13.23228 -12.97905 -12.94206 -12.88285 -12.77686 -12.57748 -12.37111 -12.29210 -12.26871 -12.06348 -11.95399 -11.84002 -11.57671 -11.26636 -10.90397 -10.38479 -10.19280 -10.10153 -10.06682 -9.99157 -9.62411 -9.52835 -5.06375 -0.31546 -0.14258 0.03247 0.04433 0.77763 1.17726 1.26654 1.57910 1.59637 1.62709 1.92471 1.94097 2.21212 2.33700 2.66281 2.78456 2.97163 3.39305 3.41220 3.51098 3.54276 3.61762 3.64047 3.68843 3.72672 3.76118 3.90304 4.04035 4.06517 4.07867 4.09263 4.12527 4.23349 4.24253 4.25409 4.30018 4.43952 4.51168 4.55595 4.59564 4.67185 4.71131 4.76058 4.83615 4.86583 4.90759 4.96506 5.00547 5.00765 5.02781 5.07553 5.16522 5.21833 5.39263 5.43493 5.56532 5.63431 5.68663 6.18903 6.56066 6.85381 7.28308 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.034 3.966 2 N -0.502 5.502 3 C 0.037 3.963 4 C 0.515 3.485 5 O -0.560 6.560 6 N -0.701 5.701 7 C 0.119 3.881 8 H 0.082 0.918 9 C -0.124 4.124 10 C -0.062 4.062 11 H 0.128 0.872 12 C -0.052 4.052 13 C -0.103 4.103 14 C -0.045 4.045 15 Cl -0.063 7.063 16 C -0.053 4.053 17 Cl -0.063 7.063 18 C -0.100 4.100 19 C -0.093 4.093 20 C -0.147 4.147 21 C 0.141 3.859 22 C -0.104 4.104 23 C -0.045 4.045 24 O -0.543 6.543 25 Si 0.746 3.254 26 C -0.060 4.060 27 C -0.186 4.186 28 C 0.509 3.491 29 F -0.178 7.178 30 F -0.177 7.177 31 F -0.182 7.182 32 C -0.053 4.053 33 C -0.143 4.143 34 H 0.073 0.927 35 H 0.045 0.955 36 H 0.066 0.934 37 H 0.101 0.899 38 H 0.121 0.879 39 H 0.418 0.582 40 H 0.095 0.905 41 H 0.095 0.905 42 H 0.148 0.852 43 H 0.158 0.842 44 H 0.169 0.831 45 H 0.069 0.931 46 H 0.090 0.910 47 H 0.102 0.898 48 H 0.389 0.611 49 H -0.298 1.298 50 H -0.280 1.280 51 H 0.152 0.848 52 H 0.157 0.843 53 H 0.157 0.843 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 0.255 2.278 8.347 8.656 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.111 4.111 2 N -0.223 5.223 3 C -0.092 4.092 4 C 0.301 3.699 5 O -0.437 6.437 6 N -0.354 5.354 7 C 0.013 3.987 8 H 0.100 0.900 9 C -0.162 4.162 10 C -0.080 4.080 11 H 0.146 0.854 12 C -0.053 4.053 13 C -0.121 4.121 14 C -0.072 4.072 15 Cl -0.034 7.034 16 C -0.081 4.081 17 Cl -0.034 7.034 18 C -0.118 4.118 19 C -0.111 4.111 20 C -0.186 4.186 21 C 0.047 3.953 22 C -0.105 4.105 23 C -0.155 4.155 24 O -0.350 6.350 25 Si 0.652 3.348 26 C -0.078 4.078 27 C -0.188 4.188 28 C 0.455 3.545 29 F -0.160 7.160 30 F -0.159 7.159 31 F -0.163 7.163 32 C -0.072 4.072 33 C -0.163 4.163 34 H 0.091 0.909 35 H 0.064 0.936 36 H 0.085 0.915 37 H 0.119 0.881 38 H 0.139 0.861 39 H 0.255 0.745 40 H 0.114 0.886 41 H 0.114 0.886 42 H 0.166 0.834 43 H 0.175 0.825 44 H 0.186 0.814 45 H 0.088 0.912 46 H 0.108 0.892 47 H 0.120 0.880 48 H 0.238 0.762 49 H -0.226 1.226 50 H -0.206 1.206 51 H 0.170 0.830 52 H 0.174 0.826 53 H 0.174 0.826 Dipole moment (debyes) X Y Z Total from point charges 0.794 1.995 6.366 6.719 hybrid contribution -0.036 1.172 0.169 1.184 sum 0.758 3.167 6.535 7.301 Atomic orbital electron populations 1.22468 0.84139 1.02872 1.01665 1.59088 1.09850 1.12717 1.40637 1.22172 1.03360 0.88466 0.95196 1.20796 0.78495 0.81673 0.88910 1.90734 1.40583 1.64676 1.47727 1.45503 1.57532 1.13072 1.19289 1.21751 0.97944 0.93739 0.85263 0.90012 1.23371 0.95139 0.93462 1.04230 1.21677 0.97675 0.91060 0.97592 0.85434 1.19942 0.95300 0.96560 0.93467 1.21062 1.01047 0.96856 0.93167 1.20953 0.98739 0.98364 0.89140 1.98404 1.82236 1.97964 1.24792 1.21052 1.00004 0.92932 0.94063 1.98402 1.91443 1.49574 1.64003 1.21354 0.99861 1.00014 0.90615 1.21500 0.96160 0.92546 1.00845 1.23573 0.98676 0.96693 0.99620 1.18870 0.94346 0.93129 0.88996 1.20239 0.99011 0.95556 0.95705 1.24832 0.96238 1.02933 0.91539 1.86588 1.42690 1.50453 1.55275 0.94457 0.76105 0.83507 0.80745 1.20771 0.97419 0.99193 0.90408 1.19943 1.06678 0.99757 0.92374 1.18353 0.75536 0.77941 0.82705 1.93176 1.91188 1.82313 1.49290 1.93235 1.86467 1.39557 1.96595 1.93228 1.30350 1.96369 1.96351 1.21142 0.93006 0.91682 1.01320 1.21896 1.01056 1.01792 0.91555 0.90852 0.93626 0.91470 0.88141 0.86140 0.74494 0.88606 0.88606 0.83396 0.82479 0.81385 0.91203 0.89198 0.88018 0.76219 1.22608 1.20626 0.83046 0.82569 0.82568 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.03 0.00 8.94 127.69 1.14 1.14 16 2 N -0.50 0.47 3.71 -724.42 -2.69 -2.22 16 3 C 0.04 -0.23 5.03 85.55 0.43 0.20 16 4 C 0.51 -0.54 7.55 87.66 0.66 0.12 16 5 O -0.56 -5.21 13.91 -3.04 -0.04 -5.25 16 6 N -0.70 5.12 5.26 -440.04 -2.32 2.81 16 7 C 0.12 -0.61 4.32 45.52 0.20 -0.41 16 8 H 0.08 0.05 7.46 -2.39 -0.02 0.04 16 9 C -0.12 1.15 7.29 31.12 0.23 1.38 16 10 C -0.06 0.59 3.85 -11.01 -0.04 0.55 16 11 H 0.13 -1.76 8.09 -2.39 -0.02 -1.78 16 12 C -0.05 0.29 5.29 -19.86 -0.11 0.18 16 13 C -0.10 0.18 8.41 22.29 0.19 0.36 16 14 C -0.04 -0.11 6.32 22.34 0.14 0.04 16 15 Cl -0.06 -0.30 27.96 -2.72 -0.08 -0.37 16 16 C -0.05 -0.11 6.32 22.35 0.14 0.04 16 17 Cl -0.06 -0.26 27.96 -2.72 -0.08 -0.34 16 18 C -0.10 0.29 9.83 22.29 0.22 0.51 16 19 C -0.09 0.61 9.71 22.27 0.22 0.82 16 20 C -0.15 0.79 7.25 31.12 0.23 1.02 16 21 C 0.14 -0.40 4.65 38.21 0.18 -0.22 16 22 C -0.10 0.20 5.28 -19.73 -0.10 0.09 16 23 C -0.04 -0.02 2.34 71.24 0.17 0.15 16 24 O -0.54 -1.06 12.88 -148.98 -1.92 -2.98 16 25 Si 0.75 -0.70 29.92 68.60 2.05 1.35 16 26 C -0.06 0.20 8.20 22.25 0.18 0.38 16 27 C -0.19 0.54 5.37 -19.84 -0.11 0.43 16 28 C 0.51 -0.29 3.31 71.24 0.24 -0.05 16 29 F -0.18 -0.01 15.97 44.97 0.72 0.70 16 30 F -0.18 -0.12 17.23 44.97 0.78 0.65 16 31 F -0.18 -0.28 17.24 44.97 0.78 0.49 16 32 C -0.05 0.30 9.42 22.25 0.21 0.51 16 33 C -0.14 0.94 8.70 22.40 0.19 1.14 16 34 H 0.07 0.37 5.88 -2.39 -0.01 0.35 16 35 H 0.05 -0.21 7.48 -2.39 -0.02 -0.22 16 36 H 0.07 0.11 8.13 -2.38 -0.02 0.09 16 37 H 0.10 -1.20 7.47 -2.39 -0.02 -1.22 16 38 H 0.12 -1.08 8.06 -2.39 -0.02 -1.10 16 39 H 0.42 -5.93 8.70 -92.71 -0.81 -6.74 16 40 H 0.10 -0.75 7.73 -2.39 -0.02 -0.77 16 41 H 0.10 -1.13 8.14 -2.39 -0.02 -1.15 16 42 H 0.15 -0.51 6.65 -2.91 -0.02 -0.52 16 43 H 0.16 -0.98 8.06 -2.92 -0.02 -1.00 16 44 H 0.17 -2.14 8.05 -2.91 -0.02 -2.16 16 45 H 0.07 -0.24 8.14 -2.39 -0.02 -0.26 16 46 H 0.09 -0.36 6.31 -2.39 -0.02 -0.38 16 47 H 0.10 0.42 7.78 -2.39 -0.02 0.40 16 48 H 0.39 -2.54 8.70 -74.06 -0.64 -3.18 16 49 H -0.30 -2.17 7.11 99.48 0.71 -1.46 16 50 H -0.28 -0.92 7.11 99.48 0.71 -0.21 16 51 H 0.15 -0.95 5.12 -2.91 -0.01 -0.96 16 52 H 0.16 -1.33 7.15 -2.91 -0.02 -1.35 16 53 H 0.16 -1.73 7.16 -2.91 -0.02 -1.75 16 Total: 0.00 -23.56 473.93 1.42 -22.14 By element: Atomic # 1 Polarization: -24.97 SS G_CDS: -0.38 Total: -25.34 kcal Atomic # 6 Polarization: 3.77 SS G_CDS: 4.60 Total: 8.37 kcal Atomic # 7 Polarization: 5.59 SS G_CDS: -5.01 Total: 0.58 kcal Atomic # 8 Polarization: -6.27 SS G_CDS: -1.96 Total: -8.23 kcal Atomic # 9 Polarization: -0.42 SS G_CDS: 2.27 Total: 1.85 kcal Atomic # 14 Polarization: -0.70 SS G_CDS: 2.05 Total: 1.35 kcal Atomic # 17 Polarization: -0.56 SS G_CDS: -0.15 Total: -0.71 kcal Total: -23.56 1.42 -22.14 kcal The number of atoms in the molecule is 53 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. s_22_9750034_6406254.mol2 53 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 -97.210 kcal (2) G-P(sol) polarization free energy of solvation -23.559 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -120.768 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 1.423 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -22.136 kcal (6) G-S(sol) free energy of system = (1) + (5) -119.346 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.98 seconds