Wall clock time and date at job start Sat Apr 11 2020 06:09:07 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.53006 * 1 3 3 C 1.52998 * 109.46713 * 2 1 4 4 H 1.08998 * 109.47514 * 304.99989 * 3 2 1 5 5 C 1.50704 * 109.47041 * 185.00043 * 3 2 1 6 6 H 1.08003 * 119.99791 * 240.00161 * 5 3 2 7 7 C 1.37878 * 119.99855 * 59.99944 * 5 3 2 8 8 N 1.31509 * 120.00035 * 0.02562 * 7 5 3 9 9 Si 1.86803 * 119.99571 * 179.97438 * 7 5 3 10 10 H 1.48496 * 109.47359 * 0.02562 * 9 7 5 11 11 H 1.48502 * 109.46817 * 120.00079 * 9 7 5 12 12 H 1.48493 * 109.47195 * 239.99964 * 9 7 5 13 13 N 1.46503 * 109.46984 * 65.00125 * 3 2 1 14 14 C 1.34774 * 120.00062 * 274.99237 * 13 3 2 15 15 O 1.21287 * 119.99796 * 0.02562 * 14 13 3 16 16 C 1.50695 * 120.00178 * 180.02562 * 14 13 3 17 17 H 1.09002 * 110.64209 * 336.42091 * 16 14 13 18 18 C 1.55162 * 110.72325 * 213.28307 * 16 14 13 19 19 C 1.54914 * 101.57837 * 276.96272 * 18 16 14 20 20 N 1.47423 * 104.82804 * 322.98687 * 19 18 16 21 21 C 1.34779 * 125.64537 * 204.09127 * 20 19 18 22 22 O 1.21515 * 120.00075 * 179.75332 * 21 20 19 23 23 O 1.34641 * 119.99839 * 0.03141 * 21 20 19 24 24 C 1.45203 * 116.99908 * 179.97438 * 23 21 20 25 25 C 1.53008 * 109.47115 * 59.99604 * 24 23 21 26 26 C 1.52997 * 109.47118 * 179.97438 * 24 23 21 27 27 C 1.53003 * 109.47190 * 300.00131 * 24 23 21 28 28 C 1.47019 * 108.70782 * 24.35978 * 20 19 18 29 29 H 1.08996 * 109.88108 * 239.29656 * 28 20 19 30 30 C 1.53000 * 109.88258 * 118.12421 * 28 20 19 31 31 H 1.08998 * 109.47350 * 299.99597 * 1 2 3 32 32 H 1.09005 * 109.47019 * 60.00163 * 1 2 3 33 33 H 1.09003 * 109.46738 * 179.97438 * 1 2 3 34 34 H 1.08998 * 109.47350 * 240.00403 * 2 1 3 35 35 H 1.08998 * 109.47302 * 120.00136 * 2 1 3 36 36 H 0.97000 * 119.99648 * 179.97438 * 8 7 5 37 37 H 0.96994 * 119.99931 * 95.00208 * 13 3 2 38 38 H 1.08991 * 111.00518 * 158.88766 * 18 16 14 39 39 H 1.08998 * 110.99941 * 35.02882 * 18 16 14 40 40 H 1.09001 * 110.36439 * 204.15028 * 19 18 16 41 41 H 1.08992 * 110.36824 * 81.82567 * 19 18 16 42 42 H 1.08999 * 109.46844 * 59.99765 * 25 24 23 43 43 H 1.08995 * 109.47093 * 179.97438 * 25 24 23 44 44 H 1.08996 * 109.47088 * 299.99997 * 25 24 23 45 45 H 1.09003 * 109.46689 * 60.00325 * 26 24 23 46 46 H 1.09004 * 109.46718 * 179.97438 * 26 24 23 47 47 H 1.08995 * 109.47393 * 299.99972 * 26 24 23 48 48 H 1.08999 * 109.47286 * 59.99218 * 27 24 23 49 49 H 1.08995 * 109.47094 * 180.02562 * 27 24 23 50 50 H 1.09007 * 109.46685 * 299.99588 * 27 24 23 51 51 H 1.09003 * 109.46622 * 65.88198 * 30 28 20 52 52 H 1.09000 * 109.47045 * 185.87627 * 30 28 20 53 53 H 1.08997 * 109.47380 * 305.87807 * 30 28 20 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.5301 0.0000 0.0000 3 6 2.0400 1.4425 0.0000 4 1 1.6053 1.9816 -0.8418 5 6 3.5419 1.4444 -0.1238 6 1 4.1471 1.8890 0.6523 7 6 4.1433 0.8785 -1.2280 8 7 3.4062 0.3376 -2.1733 9 14 6.0049 0.8815 -1.3818 10 1 6.5987 1.5466 -0.1942 11 1 6.4971 -0.5172 -1.4629 12 1 6.3990 1.6170 -2.6101 13 7 1.6523 2.0974 1.2519 14 6 0.4351 2.6644 1.3670 15 8 -0.3396 2.6320 0.4344 16 6 0.0366 3.3386 2.6545 17 1 0.6159 2.9449 3.4897 18 6 -1.4840 3.1713 2.9140 19 6 -2.0738 4.3195 2.0576 20 7 -1.1115 5.4274 2.1978 21 6 -1.3901 6.7356 2.0317 22 8 -0.5154 7.5664 2.1774 23 8 -2.6390 7.1223 1.7097 24 6 -2.8666 8.5471 1.5470 25 6 -1.9693 9.0826 0.4293 26 6 -4.3328 8.7912 1.1844 27 6 -2.5372 9.2688 2.8553 28 6 0.2039 4.8687 2.5427 29 1 0.5401 5.2737 3.4971 30 6 1.2160 5.2053 1.4457 31 1 -0.3634 0.5137 0.8900 32 1 -0.3633 0.5138 -0.8900 33 1 -0.3633 -1.0277 -0.0005 34 1 1.8934 -0.5137 -0.8900 35 1 1.8934 -0.5138 0.8899 36 1 3.8293 -0.0602 -2.9502 37 1 2.2718 2.1232 1.9978 38 1 -1.7179 3.3111 3.9693 39 1 -1.8359 2.2013 2.5626 40 1 -3.0508 4.6148 2.4403 41 1 -2.1507 4.0136 1.0143 42 1 -2.2039 8.5684 -0.5027 43 1 -2.1405 10.1520 0.3069 44 1 -0.9247 8.9086 0.6873 45 1 -4.9722 8.4096 1.9805 46 1 -4.5037 9.8608 1.0628 47 1 -4.5670 8.2771 0.2523 48 1 -1.4926 9.0948 3.1132 49 1 -2.7085 10.3383 2.7333 50 1 -3.1765 8.8873 3.6515 51 1 1.3775 6.2830 1.4188 52 1 2.1597 4.7016 1.6551 53 1 0.8319 4.8712 0.4819 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 ZINC000492695021.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:09:07 Heat of formation + Delta-G solvation = -124.443498 kcal Electronic energy + Delta-G solvation = -32129.901287 eV Core-core repulsion = 28001.542754 eV Total energy + Delta-G solvation = -4128.358534 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.69 seconds Orbital eigenvalues (eV) -40.98411 -39.90178 -37.87851 -36.35889 -34.94868 -34.16983 -32.52818 -31.73374 -31.28750 -30.58026 -27.56902 -27.49762 -27.18430 -25.91858 -23.88451 -23.18959 -22.23064 -21.25671 -19.70764 -19.37545 -19.06863 -17.79036 -17.37157 -16.52117 -16.32575 -15.88984 -15.30122 -15.23286 -14.90795 -14.66461 -14.35343 -13.98222 -13.86171 -13.47041 -13.26255 -13.10481 -12.86498 -12.67796 -12.62480 -12.38056 -12.21081 -12.03641 -12.00936 -11.87833 -11.59552 -11.50319 -11.45391 -11.33698 -11.32882 -11.28099 -11.13384 -10.95613 -10.85208 -10.75817 -10.41639 -10.31087 -10.18042 -9.95637 -9.64513 -9.06007 -8.89981 -8.65555 -8.51162 -7.83164 -6.13857 -4.22992 2.08879 2.91005 3.22536 3.28237 3.31220 3.33395 3.86779 4.00930 4.35502 4.48021 4.57687 4.61800 4.63633 4.76278 4.98613 5.17953 5.22428 5.29240 5.31096 5.37947 5.42137 5.46406 5.51554 5.52058 5.61882 5.70609 5.72356 5.77159 5.77962 5.82658 5.97646 6.14844 6.16497 6.19905 6.22970 6.34072 6.36876 6.43231 6.45965 6.58988 6.64122 7.01336 7.17698 7.34766 7.43952 7.55725 7.60686 7.71662 7.88949 8.16882 8.39550 8.65402 8.76812 9.11863 9.44712 10.91403 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.016851 B = 0.002831 C = 0.002648 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1661.258355 B = 9889.070918 C =10570.108187 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.136 4.136 2 C -0.128 4.128 3 C 0.272 3.728 4 H 0.048 0.952 5 C -0.533 4.533 6 H 0.061 0.939 7 C 0.067 3.933 8 N -0.882 5.882 9 Si 0.899 3.101 10 H -0.273 1.273 11 H -0.283 1.283 12 H -0.283 1.283 13 N -0.706 5.706 14 C 0.516 3.484 15 O -0.543 6.543 16 C -0.127 4.127 17 H 0.105 0.895 18 C -0.111 4.111 19 C 0.105 3.895 20 N -0.596 5.596 21 C 0.670 3.330 22 O -0.570 6.570 23 O -0.362 6.362 24 C 0.136 3.864 25 C -0.184 4.184 26 C -0.140 4.140 27 C -0.182 4.182 28 C 0.175 3.825 29 H 0.084 0.916 30 C -0.172 4.172 31 H 0.038 0.962 32 H 0.049 0.951 33 H 0.021 0.979 34 H 0.098 0.902 35 H 0.030 0.970 36 H 0.263 0.737 37 H 0.387 0.613 38 H 0.079 0.921 39 H 0.092 0.908 40 H 0.074 0.926 41 H 0.089 0.911 42 H 0.060 0.940 43 H 0.068 0.932 44 H 0.088 0.912 45 H 0.065 0.935 46 H 0.078 0.922 47 H 0.067 0.933 48 H 0.088 0.912 49 H 0.068 0.932 50 H 0.057 0.943 51 H 0.069 0.931 52 H 0.074 0.926 53 H 0.069 0.931 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -12.535 11.669 11.367 20.555 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.193 4.193 2 C -0.167 4.167 3 C 0.170 3.830 4 H 0.066 0.934 5 C -0.571 4.571 6 H 0.079 0.921 7 C -0.137 4.137 8 N -0.519 5.519 9 Si 0.727 3.273 10 H -0.197 1.197 11 H -0.209 1.209 12 H -0.209 1.209 13 N -0.353 5.353 14 C 0.300 3.700 15 O -0.422 6.422 16 C -0.149 4.149 17 H 0.123 0.877 18 C -0.149 4.149 19 C -0.017 4.017 20 N -0.336 5.336 21 C 0.425 3.575 22 O -0.460 6.460 23 O -0.277 6.277 24 C 0.099 3.901 25 C -0.241 4.241 26 C -0.197 4.197 27 C -0.239 4.239 28 C 0.072 3.928 29 H 0.102 0.898 30 C -0.230 4.230 31 H 0.057 0.943 32 H 0.068 0.932 33 H 0.040 0.960 34 H 0.116 0.884 35 H 0.049 0.951 36 H 0.074 0.926 37 H 0.220 0.780 38 H 0.097 0.903 39 H 0.110 0.890 40 H 0.093 0.907 41 H 0.107 0.893 42 H 0.079 0.921 43 H 0.087 0.913 44 H 0.107 0.893 45 H 0.084 0.916 46 H 0.097 0.903 47 H 0.086 0.914 48 H 0.107 0.893 49 H 0.087 0.913 50 H 0.076 0.924 51 H 0.088 0.912 52 H 0.093 0.907 53 H 0.088 0.912 Dipole moment (debyes) X Y Z Total from point charges -12.850 11.941 10.499 20.444 hybrid contribution 1.359 -0.052 0.200 1.375 sum -11.491 11.890 10.699 19.694 Atomic orbital electron populations 1.21831 0.98361 0.97957 1.01162 1.21934 0.92600 0.99883 1.02235 1.18761 0.94141 0.89088 0.81003 0.93374 1.21266 0.90764 1.38855 1.06231 0.92070 1.24877 0.98515 0.94840 0.95426 1.69899 1.35001 1.37713 1.09319 0.86431 0.85013 0.77746 0.78143 1.19730 1.20918 1.20920 1.46551 1.14464 1.55755 1.18536 1.19824 0.82892 0.78899 0.88400 1.90660 1.51497 1.59813 1.40205 1.22576 0.98856 0.94738 0.98732 0.87658 1.22629 0.92721 0.99690 0.99886 1.22433 0.91846 0.85370 1.02095 1.47224 1.11946 1.02714 1.71752 1.17406 0.81307 0.82209 0.76555 1.90872 1.47508 1.50422 1.57170 1.86313 1.29829 1.28642 1.82895 1.22036 0.96172 0.75750 0.96168 1.22515 1.02593 1.01234 0.97755 1.21848 0.92832 1.02294 1.02729 1.22511 1.04816 1.00379 0.96215 1.21308 0.81931 0.90579 0.98964 0.89844 1.22289 0.98707 1.03064 0.98969 0.94284 0.93193 0.95969 0.88377 0.95092 0.92648 0.78034 0.90263 0.88962 0.90734 0.89283 0.92099 0.91326 0.89295 0.91591 0.90276 0.91413 0.89300 0.91333 0.92356 0.91197 0.90737 0.91234 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.14 -2.88 8.51 37.16 0.32 -2.56 16 2 C -0.13 -2.98 5.45 -26.73 -0.15 -3.13 16 3 C 0.27 6.49 2.24 -69.08 -0.15 6.34 16 4 H 0.05 1.36 7.58 -51.93 -0.39 0.96 16 5 C -0.53 -13.59 8.73 -34.44 -0.30 -13.89 16 6 H 0.06 1.50 7.74 -52.48 -0.41 1.10 16 7 C 0.07 1.75 5.29 -17.82 -0.09 1.66 16 8 N -0.88 -24.29 11.29 55.43 0.63 -23.67 16 9 Si 0.90 20.36 29.54 -169.99 -5.02 15.34 16 10 H -0.27 -6.05 7.11 56.52 0.40 -5.65 16 11 H -0.28 -6.39 7.11 56.52 0.40 -5.99 16 12 H -0.28 -6.48 7.11 56.52 0.40 -6.08 16 13 N -0.71 -13.40 4.82 -54.93 -0.26 -13.67 16 14 C 0.52 9.07 3.29 -10.99 -0.04 9.04 16 15 O -0.54 -11.34 11.85 5.55 0.07 -11.27 16 16 C -0.13 -1.65 3.36 -89.93 -0.30 -1.96 16 17 H 0.11 1.09 8.13 -51.93 -0.42 0.67 16 18 C -0.11 -1.27 6.34 -24.58 -0.16 -1.43 16 19 C 0.10 1.37 6.52 -2.53 -0.02 1.35 16 20 N -0.60 -8.61 3.06 -168.75 -0.52 -9.13 16 21 C 0.67 10.26 7.71 54.06 0.42 10.67 16 22 O -0.57 -9.48 11.01 -19.28 -0.21 -9.69 16 23 O -0.36 -5.00 9.42 -40.33 -0.38 -5.38 16 24 C 0.14 1.51 1.13 -90.62 -0.10 1.40 16 25 C -0.18 -2.09 8.37 37.16 0.31 -1.78 16 26 C -0.14 -1.15 8.85 37.16 0.33 -0.82 16 27 C -0.18 -1.87 8.37 37.16 0.31 -1.56 16 28 C 0.17 2.42 3.38 -66.95 -0.23 2.19 16 29 H 0.08 1.04 8.14 -51.93 -0.42 0.62 16 30 C -0.17 -2.76 8.12 37.16 0.30 -2.46 16 31 H 0.04 0.78 5.64 -51.93 -0.29 0.49 16 32 H 0.05 1.14 7.85 -51.93 -0.41 0.73 16 33 H 0.02 0.42 8.14 -51.93 -0.42 -0.01 16 34 H 0.10 2.57 5.95 -51.93 -0.31 2.27 16 35 H 0.03 0.69 8.14 -51.93 -0.42 0.27 16 36 H 0.26 6.92 8.31 -40.82 -0.34 6.58 16 37 H 0.39 6.42 8.53 -40.82 -0.35 6.07 16 38 H 0.08 0.72 8.14 -51.93 -0.42 0.29 16 39 H 0.09 1.10 7.99 -51.93 -0.42 0.68 16 40 H 0.07 0.84 8.04 -51.93 -0.42 0.42 16 41 H 0.09 1.38 6.96 -51.93 -0.36 1.02 16 42 H 0.06 0.70 8.14 -51.93 -0.42 0.28 16 43 H 0.07 0.62 8.14 -51.93 -0.42 0.19 16 44 H 0.09 1.24 5.89 -51.93 -0.31 0.93 16 45 H 0.07 0.51 8.14 -51.93 -0.42 0.08 16 46 H 0.08 0.49 8.14 -51.93 -0.42 0.07 16 47 H 0.07 0.56 8.14 -51.93 -0.42 0.14 16 48 H 0.09 1.13 5.87 -51.93 -0.30 0.82 16 49 H 0.07 0.56 8.14 -51.93 -0.42 0.14 16 50 H 0.06 0.56 8.14 -51.93 -0.42 0.14 16 51 H 0.07 1.15 7.04 -51.93 -0.37 0.79 16 52 H 0.07 1.15 7.87 -51.93 -0.41 0.74 16 53 H 0.07 1.33 7.02 -51.93 -0.36 0.96 16 LS Contribution 403.95 15.07 6.09 6.09 Total: -1.00 -30.09 403.95 -8.47 -38.57 By element: Atomic # 1 Polarization: 19.03 SS G_CDS: -9.31 Total: 9.73 kcal Atomic # 6 Polarization: 2.63 SS G_CDS: 0.45 Total: 3.08 kcal Atomic # 7 Polarization: -46.31 SS G_CDS: -0.16 Total: -46.46 kcal Atomic # 8 Polarization: -25.81 SS G_CDS: -0.53 Total: -26.34 kcal Atomic # 14 Polarization: 20.36 SS G_CDS: -5.02 Total: 15.34 kcal Total LS contribution 6.09 Total: 6.09 kcal Total: -30.09 -8.47 -38.57 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. ZINC000492695021.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 -85.878 kcal (2) G-P(sol) polarization free energy of solvation -30.094 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -115.972 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.471 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -38.565 kcal (6) G-S(sol) free energy of system = (1) + (5) -124.443 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.69 seconds