Wall clock time and date at job start Sat Apr 11 2020 06:10:50 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.53002 * 1 3 3 H 1.08997 * 109.46640 * 2 1 4 4 C 1.50701 * 109.46921 * 240.00046 * 2 1 3 5 5 H 1.08002 * 120.00258 * 274.99625 * 4 2 1 6 6 C 1.37884 * 120.00120 * 95.00144 * 4 2 1 7 7 N 1.31505 * 119.99965 * 184.91245 * 6 4 2 8 8 Si 1.86802 * 119.99826 * 4.91268 * 6 4 2 9 9 H 1.48504 * 109.46936 * 84.83366 * 8 6 4 10 10 H 1.48499 * 109.46887 * 204.83265 * 8 6 4 11 11 H 1.48496 * 109.47224 * 324.83342 * 8 6 4 12 12 N 1.46494 * 109.47363 * 119.99762 * 2 1 3 13 13 C 1.34775 * 120.00474 * 275.00142 * 12 2 1 14 14 O 1.21285 * 119.99515 * 0.02562 * 13 12 2 15 15 C 1.50696 * 120.00472 * 179.97438 * 13 12 2 16 16 H 1.09010 * 109.58437 * 319.71700 * 15 13 12 17 17 C 1.53184 * 109.58549 * 80.00029 * 15 13 12 18 18 C 1.53036 * 109.31404 * 65.15527 * 17 15 13 19 19 H 1.09006 * 109.46006 * 301.35661 * 18 17 15 20 20 C 1.52998 * 109.46202 * 181.38867 * 18 17 15 21 21 C 1.53036 * 109.53743 * 61.36353 * 18 17 15 22 22 C 1.53187 * 109.31316 * 298.63193 * 21 18 17 23 23 N 1.46922 * 109.58742 * 199.29428 * 15 13 12 24 24 C 1.34773 * 120.63693 * 113.82120 * 23 15 13 25 25 O 1.21518 * 119.99680 * 0.02562 * 24 23 15 26 26 O 1.34628 * 120.00374 * 179.97438 * 24 23 15 27 27 C 1.45197 * 117.00381 * 179.97438 * 26 24 23 28 28 C 1.53001 * 109.46753 * 60.00286 * 27 26 24 29 29 C 1.52998 * 109.47586 * 180.02562 * 27 26 24 30 30 C 1.53001 * 109.47371 * 300.00910 * 27 26 24 31 31 H 1.09006 * 109.46962 * 175.03937 * 1 2 3 32 32 H 1.08996 * 109.47356 * 295.03097 * 1 2 3 33 33 H 1.08998 * 109.47235 * 55.04312 * 1 2 3 34 34 H 0.96999 * 120.00415 * 179.97438 * 7 6 4 35 35 H 0.97007 * 119.99983 * 95.27136 * 12 2 1 36 36 H 1.09004 * 109.49853 * 185.10496 * 17 15 13 37 37 H 1.09007 * 109.49982 * 305.16967 * 17 15 13 38 38 H 1.08999 * 109.47790 * 59.99795 * 20 18 17 39 39 H 1.09003 * 109.47533 * 180.02562 * 20 18 17 40 40 H 1.09007 * 109.46923 * 299.99382 * 20 18 17 41 41 H 1.09001 * 109.49922 * 58.59270 * 21 18 17 42 42 H 1.09005 * 109.49654 * 178.68324 * 21 18 17 43 43 H 1.08998 * 109.58646 * 174.42134 * 22 21 18 44 44 H 1.08999 * 109.59411 * 294.70844 * 22 21 18 45 45 H 1.09001 * 109.47292 * 60.00909 * 28 27 26 46 46 H 1.09009 * 109.47038 * 180.02562 * 28 27 26 47 47 H 1.08995 * 109.47767 * 299.99983 * 28 27 26 48 48 H 1.08997 * 109.47242 * 60.00058 * 29 27 26 49 49 H 1.09000 * 109.47356 * 180.02562 * 29 27 26 50 50 H 1.09003 * 109.46731 * 300.00009 * 29 27 26 51 51 H 1.09002 * 109.47582 * 59.99440 * 30 27 26 52 52 H 1.09000 * 109.47124 * 179.97438 * 30 27 26 53 53 H 1.08997 * 109.46628 * 299.99434 * 30 27 26 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 1 1.8933 1.0277 0.0000 4 6 2.0323 -0.7104 -1.2305 5 1 2.1355 -1.7855 -1.2291 6 6 2.3603 0.0121 -2.3581 7 7 2.7059 -0.6131 -3.4622 8 14 2.3135 1.8791 -2.3175 9 1 3.5991 2.3949 -1.7822 10 1 2.1057 2.3991 -3.6929 11 1 1.1986 2.3275 -1.4451 12 7 2.0184 -0.6905 1.1961 13 6 2.1471 -0.0182 2.3571 14 8 1.8582 1.1584 2.4117 15 6 2.6490 -0.7287 3.5877 16 1 2.2111 -1.7257 3.6381 17 6 4.1756 -0.8406 3.5293 18 6 4.7887 0.5600 3.5956 19 1 4.4150 1.1597 2.7656 20 6 6.3122 0.4545 3.5028 21 6 4.4034 1.2253 4.9189 22 6 2.8793 1.3607 4.9929 23 7 2.2763 0.0370 4.7849 24 6 1.3981 -0.4651 5.6755 25 8 0.9067 -1.5613 5.4923 26 8 1.0660 0.2473 6.7685 27 6 0.1138 -0.3554 7.6841 28 6 0.6785 -1.6764 8.2103 29 6 -0.1385 0.5950 8.8562 30 6 -1.2024 -0.6198 6.9501 31 1 -0.3633 -1.0239 -0.0889 32 1 -0.3634 0.4348 0.9311 33 1 -0.3633 0.5888 -0.8422 34 1 2.9363 -0.1049 -4.2556 35 1 2.2529 -1.6308 1.1516 36 1 4.5295 -1.4329 4.3732 37 1 4.4697 -1.3232 2.5973 38 1 6.5870 -0.0181 2.5598 39 1 6.7491 1.4520 3.5504 40 1 6.6861 -0.1454 4.3325 41 1 4.7551 0.6135 5.7497 42 1 4.8596 2.2136 4.9756 43 1 2.5935 1.7424 5.9731 44 1 2.5353 2.0469 4.2191 45 1 1.6161 -1.4881 8.7334 46 1 -0.0366 -2.1292 8.8974 47 1 0.8579 -2.3534 7.3752 48 1 -0.5409 1.5359 8.4810 49 1 -0.8530 0.1423 9.5438 50 1 0.7991 0.7836 9.3793 51 1 -1.0230 -1.2967 6.1148 52 1 -1.9170 -1.0727 7.6373 53 1 -1.6046 0.3213 6.5751 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 ZINC000495493352.mol2 53 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:10:50 Heat of formation + Delta-G solvation = -133.059781 kcal Electronic energy + Delta-G solvation = -32449.800277 eV Core-core repulsion = 28321.068114 eV Total energy + Delta-G solvation = -4128.732164 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.60 seconds Orbital eigenvalues (eV) -41.08818 -39.77446 -37.74089 -36.68114 -34.87773 -34.46009 -33.50214 -32.50778 -31.51613 -29.50564 -28.28109 -27.57938 -27.50009 -25.34255 -24.53573 -23.40035 -22.46106 -20.99529 -20.35347 -19.60174 -19.11943 -17.92442 -17.26750 -16.67255 -16.26909 -15.82049 -15.58719 -15.09736 -15.02609 -14.59671 -14.43729 -14.08139 -13.96925 -13.84255 -13.52056 -13.37638 -13.28165 -12.74316 -12.55585 -12.52173 -12.40119 -12.34615 -12.13445 -12.04177 -11.96062 -11.84941 -11.68086 -11.46933 -11.44563 -11.30467 -11.11587 -10.89940 -10.84694 -10.66316 -10.46934 -10.45615 -10.40689 -9.96020 -9.90366 -9.53193 -9.21951 -8.74072 -8.37809 -8.10373 -6.25267 -4.28892 2.05110 2.84697 3.18051 3.24832 3.31939 3.39029 3.76504 3.99677 4.30343 4.56477 4.60736 4.70454 4.74944 4.76278 4.79690 5.01690 5.09325 5.23239 5.28443 5.33001 5.37321 5.45470 5.47136 5.51193 5.58036 5.60865 5.65253 5.66262 5.71969 5.75693 5.80609 5.85931 5.89587 6.04703 6.11635 6.15339 6.19942 6.23172 6.33045 6.42238 6.42646 6.47744 6.54268 6.68099 6.92016 7.18389 7.35243 7.41250 7.53462 8.10132 8.35247 8.54212 8.75426 8.92451 9.36623 10.86446 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.013968 B = 0.003470 C = 0.003086 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2004.033388 B = 8067.078106 C = 9071.575290 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.147 4.147 2 C 0.233 3.767 3 H 0.050 0.950 4 C -0.530 4.530 5 H 0.079 0.921 6 C 0.082 3.918 7 N -0.888 5.888 8 Si 0.874 3.126 9 H -0.273 1.273 10 H -0.286 1.286 11 H -0.261 1.261 12 N -0.698 5.698 13 C 0.497 3.503 14 O -0.550 6.550 15 C 0.137 3.863 16 H 0.120 0.880 17 C -0.116 4.116 18 C -0.099 4.099 19 H 0.082 0.918 20 C -0.144 4.144 21 C -0.124 4.124 22 C 0.105 3.895 23 N -0.570 5.570 24 C 0.657 3.343 25 O -0.573 6.573 26 O -0.360 6.360 27 C 0.135 3.865 28 C -0.182 4.182 29 C -0.140 4.140 30 C -0.184 4.184 31 H 0.047 0.953 32 H 0.028 0.972 33 H 0.051 0.949 34 H 0.268 0.732 35 H 0.394 0.606 36 H 0.076 0.924 37 H 0.088 0.912 38 H 0.058 0.942 39 H 0.054 0.946 40 H 0.052 0.948 41 H 0.062 0.938 42 H 0.075 0.925 43 H 0.085 0.915 44 H 0.096 0.904 45 H 0.058 0.942 46 H 0.068 0.932 47 H 0.088 0.912 48 H 0.066 0.934 49 H 0.079 0.921 50 H 0.065 0.935 51 H 0.089 0.911 52 H 0.068 0.932 53 H 0.059 0.941 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 1.229 1.070 23.177 23.234 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.206 4.206 2 C 0.131 3.869 3 H 0.068 0.932 4 C -0.568 4.568 5 H 0.097 0.903 6 C -0.123 4.123 7 N -0.524 5.524 8 Si 0.701 3.299 9 H -0.199 1.199 10 H -0.212 1.212 11 H -0.185 1.185 12 N -0.347 5.347 13 C 0.280 3.720 14 O -0.429 6.429 15 C 0.030 3.970 16 H 0.138 0.862 17 C -0.154 4.154 18 C -0.117 4.117 19 H 0.100 0.900 20 C -0.201 4.201 21 C -0.162 4.162 22 C -0.017 4.017 23 N -0.309 5.309 24 C 0.412 3.588 25 O -0.463 6.463 26 O -0.275 6.275 27 C 0.098 3.902 28 C -0.239 4.239 29 C -0.197 4.197 30 C -0.241 4.241 31 H 0.066 0.934 32 H 0.047 0.953 33 H 0.070 0.930 34 H 0.078 0.922 35 H 0.227 0.773 36 H 0.094 0.906 37 H 0.106 0.894 38 H 0.077 0.923 39 H 0.074 0.926 40 H 0.071 0.929 41 H 0.081 0.919 42 H 0.093 0.907 43 H 0.103 0.897 44 H 0.114 0.886 45 H 0.077 0.923 46 H 0.087 0.913 47 H 0.107 0.893 48 H 0.085 0.915 49 H 0.097 0.903 50 H 0.084 0.916 51 H 0.107 0.893 52 H 0.087 0.913 53 H 0.079 0.921 Dipole moment (debyes) X Y Z Total from point charges 0.862 1.606 23.005 23.078 hybrid contribution 0.040 0.455 -0.195 0.496 sum 0.902 2.060 22.811 22.921 Atomic orbital electron populations 1.22212 0.98608 1.00625 0.99155 1.19118 0.90882 0.94250 0.82665 0.93244 1.21014 1.45193 0.95667 0.94918 0.90263 1.24935 0.93630 0.98565 0.95155 1.69877 1.44704 1.36186 1.01679 0.86901 0.79734 0.85930 0.77371 1.19851 1.21203 1.18543 1.46653 1.64357 1.13989 1.09682 1.20697 0.78610 0.87807 0.84853 1.90587 1.51474 1.15937 1.84890 1.21160 0.96521 0.96515 0.82766 0.86196 1.21936 0.91671 0.98104 1.03719 1.21169 0.95874 0.96131 0.98573 0.89956 1.21794 0.94791 1.01683 1.01837 1.21939 0.95434 1.01269 0.97573 1.21942 0.92664 0.83806 1.03285 1.46498 1.43628 1.19983 1.20749 1.17855 0.78957 0.82328 0.79669 1.90876 1.56687 1.21008 1.77705 1.86289 1.60304 1.48473 1.32447 1.22046 0.87266 0.93345 0.87530 1.22509 1.01202 0.97292 1.02929 1.21846 1.02574 0.99431 0.95820 1.22516 0.95696 1.03356 1.02546 0.93377 0.95333 0.92951 0.92228 0.77271 0.90584 0.89379 0.92267 0.92649 0.92907 0.91880 0.90664 0.89712 0.88586 0.92318 0.91311 0.89282 0.91463 0.90266 0.91580 0.89270 0.91300 0.92143 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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 -3.29 9.47 37.16 0.35 -2.93 16 2 C 0.23 5.49 2.22 -69.08 -0.15 5.33 16 3 H 0.05 1.26 4.73 -51.93 -0.25 1.02 16 4 C -0.53 -14.04 9.05 -34.44 -0.31 -14.35 16 5 H 0.08 2.21 8.06 -52.49 -0.42 1.79 16 6 C 0.08 2.26 5.47 -17.82 -0.10 2.16 16 7 N -0.89 -25.87 13.96 55.43 0.77 -25.10 16 8 Si 0.87 21.34 26.99 -169.99 -4.59 16.75 16 9 H -0.27 -6.62 7.11 56.52 0.40 -6.22 16 10 H -0.29 -6.87 7.11 56.52 0.40 -6.47 16 11 H -0.26 -6.32 6.82 56.52 0.39 -5.93 16 12 N -0.70 -13.86 5.07 -54.94 -0.28 -14.14 16 13 C 0.50 9.37 6.10 -10.99 -0.07 9.30 16 14 O -0.55 -11.91 13.63 5.55 0.08 -11.84 16 15 C 0.14 2.07 3.00 -68.76 -0.21 1.86 16 16 H 0.12 1.80 7.01 -51.92 -0.36 1.43 16 17 C -0.12 -1.45 4.82 -26.62 -0.13 -1.58 16 18 C -0.10 -1.24 2.52 -90.58 -0.23 -1.46 16 19 H 0.08 1.28 7.10 -51.93 -0.37 0.91 16 20 C -0.14 -1.51 9.19 37.16 0.34 -1.17 16 21 C -0.12 -1.38 5.36 -26.62 -0.14 -1.53 16 22 C 0.10 1.40 5.94 -3.71 -0.02 1.38 16 23 N -0.57 -8.50 2.54 -172.30 -0.44 -8.94 16 24 C 0.66 10.15 7.74 54.05 0.42 10.57 16 25 O -0.57 -9.68 11.30 -19.29 -0.22 -9.90 16 26 O -0.36 -4.88 8.52 -40.29 -0.34 -5.23 16 27 C 0.13 1.48 1.13 -90.62 -0.10 1.38 16 28 C -0.18 -1.87 8.37 37.16 0.31 -1.56 16 29 C -0.14 -1.11 8.85 37.16 0.33 -0.78 16 30 C -0.18 -2.10 8.37 37.16 0.31 -1.79 16 31 H 0.05 1.02 8.14 -51.93 -0.42 0.60 16 32 H 0.03 0.61 8.14 -51.93 -0.42 0.19 16 33 H 0.05 1.17 8.14 -51.93 -0.42 0.75 16 34 H 0.27 7.34 8.31 -40.82 -0.34 7.00 16 35 H 0.39 6.95 8.59 -40.82 -0.35 6.60 16 36 H 0.08 0.84 8.14 -51.93 -0.42 0.42 16 37 H 0.09 1.08 8.14 -51.93 -0.42 0.66 16 38 H 0.06 0.65 8.14 -51.93 -0.42 0.22 16 39 H 0.05 0.55 8.14 -51.93 -0.42 0.12 16 40 H 0.05 0.49 8.14 -51.93 -0.42 0.06 16 41 H 0.06 0.65 8.14 -51.93 -0.42 0.23 16 42 H 0.07 0.74 8.14 -51.93 -0.42 0.32 16 43 H 0.08 1.02 6.93 -51.93 -0.36 0.66 16 44 H 0.10 1.51 6.17 -51.93 -0.32 1.19 16 45 H 0.06 0.56 8.14 -51.93 -0.42 0.14 16 46 H 0.07 0.56 8.14 -51.92 -0.42 0.14 16 47 H 0.09 1.13 5.88 -51.93 -0.31 0.83 16 48 H 0.07 0.54 8.14 -51.93 -0.42 0.12 16 49 H 0.08 0.48 8.14 -51.93 -0.42 0.06 16 50 H 0.07 0.49 8.14 -51.93 -0.42 0.07 16 51 H 0.09 1.26 5.88 -51.93 -0.31 0.96 16 52 H 0.07 0.62 8.14 -51.93 -0.42 0.20 16 53 H 0.06 0.70 8.14 -51.93 -0.42 0.27 16 LS Contribution 407.72 15.07 6.14 6.14 Total: -1.00 -31.46 407.72 -7.65 -39.11 By element: Atomic # 1 Polarization: 17.68 SS G_CDS: -9.38 Total: 8.30 kcal Atomic # 6 Polarization: 4.23 SS G_CDS: 0.60 Total: 4.84 kcal Atomic # 7 Polarization: -48.24 SS G_CDS: 0.06 Total: -48.18 kcal Atomic # 8 Polarization: -26.47 SS G_CDS: -0.49 Total: -26.96 kcal Atomic # 14 Polarization: 21.34 SS G_CDS: -4.59 Total: 16.75 kcal Total LS contribution 6.14 Total: 6.14 kcal Total: -31.46 -7.65 -39.11 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. ZINC000495493352.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 -93.947 kcal (2) G-P(sol) polarization free energy of solvation -31.462 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -125.409 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.651 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.113 kcal (6) G-S(sol) free energy of system = (1) + (5) -133.060 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.60 seconds