4.4: Hückel Calculations Using Hulis (Part II)
- Page ID
- 371632
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)The ELN template for this part of the experiment is here: ELN01_2_Huckle2022.docx
Huckle Calculation Tasks
Perform the tasks listed here using the instructions for using Hulis that are given below.
- Perform Huckle Molecular Orbital (HMO) calculations for each of the four dyes. In the HMO analysis, omit the ethyl groups since they are not part of the conjugated network. Save a screen-shot image of the entire program window showing the structure and the orbital diagram before moving on to the next one.
- Calculate the energy difference, ∆E, between the HOMO and LUMO levels (in terms of \( \beta \)) for each dye.
- Plot the experimentally observed ∆E's (kJ/mol) (from UV-vis in Part I, plot as y or ordinate) vs the calcuated HMO ∆E's (abscissa or x).
- From the least-squares slope of the resultant line, determine a semiempirical value for \( \beta \) (\(\beta\) is the spectroscopic exchange integral). Report the uncertainty in \( \beta \).
- Compare your derived value of \(\beta\) with the accepted range for this constant (a.k.a spectroscopic exchange integral), which is – 250 to – 350 kJ/mole.9
- Use your value of \( \beta \) to calculate the energy difference, ∆E, between the HOMO and LUMO levels in units of kJ/mol and eV for each dye. Using these energies, calculate the wavelength (nm) at which you would expect to observe this HOMO \( \rightarrow \) LUMO transition for each dye, using the HMO model.
- Calculate the percent error between your calculated HMO results your experimental values. Tabulate your data in the format shown below.
WAVELENGTH (nm)
Dye #
\( \lambda_{max} \)
Literature*
\( \lambda_{max} \) Experimental
\( \lambda_{max} \) Particle-in-box
\( \lambda_{max} \) HMO
Value (nm)
Absolute Error (nm)
Value (nm)
% Error
(from exp.)
Value (nm)
% Error
(from exp.)
±
±
±
±
* reference for your literature values goes here
- Use this table as the basis for discussion of your results. Which theoretical model gives the best fit to the experimental data, and why? Are the results for each dye equally satisfactory? If your answer is no, explain why you think this is so. How are these models useful and what might be done to improve the results?
- Calculate the wavelength (nm) and minimum excitation energy (eV) of the molecule CH2=CHCH=CHCH=CHCH=CH2 using the HMO model program. In calculation use your value of \( \beta \) to convert to electron volts. What color does the molecule absorb from white light? What color does it then appear?
Instructions For Using The Program "Hulis" ver. 3.0.x
Hulis is a Java applet. Java must be installed to run the applet. You are welcome to run it from your own laptop if Java is installed. But for your convenience in the event that you do not want Java on your computer, the program is available on all of the lab computers.
Hulis and Java can be downloaded here: Huckel Theory and HuLis
- Double-click on the Hulis icon to Start the Program
- Running the Program:
The main screen shows a molecule builder. Highlight Build on the left side of the screen if it is not already highlighted when the program starts. To build a molecule, for example, Dye 1, begin by clicking anywhere in the blank window in the center of the screen. A CH3 group should appear. Click on one of the available hydrogens to add a conjugated carbon to your molecule. As you add carbons to form the first ring, you will see the orbital diagram for the molecule appear in the right-hand panel of the program.
You can follow the electron count from this pane of the program window. After putting 6 carbons to form the ring, you will need to form a bond between the two open and adjacent carbons. To do this, simply click on one of the open carbons and drag a bond over to the adjacent open carbon. The hydrogens will disappear and a delocalized C-C bond will form.
- Continue adding carbon atoms as appropriate to form the template molecule. When you are done your structure will look like this:
You can rescale the display on the screen using Zoom out from the Display menu. You can also move, rotate, and center the molecule in the display window with the appropriate control buttons in the left-hand (blue colored) pane of the program window. In the image above the Move button is currently highlighted.
- You now need to change two of the carbon atoms to nitrogen to complete the structure. Click on the Change button in the left-hand pane, and then double-click on the appropriate carbon in your structure. A Change Atom box appears. Select =NR to change this to the pyridinium group. Check to be sure the value in the Hx box = 0.51. Click on OK. Now change the other carbon to a nitrogen in the same way, but this time select the –NR2 option. Be Sure the value in the Hx box = 1.37. Click OK. Finally, make the molecule a cation by clickingh the + in the left-hand pane to increase the Total charge to 1. You final molecule will look as shown below.
- Click on the Results button, and a Results window appears with significant information about the Hückel calculation. Please note the following things:
- As you scan down the Hamiltonian matrix check to see that the Coulomb and Exchange integrals for heteroatoms and carbon-to-heteroatom bonds are as shown in the table below. If you need to see how your atoms are numbered in your structure to make it easier to find them in the Result table, put a check in the Numbering box in the left hand pane of the program window.
Atom
Bond
Number of
p-Electrons
Coulomb
Integral
Exchange
Integral
C
1
0.00
1.00
N
Pyridine
2
1.37
0.89
N
Pyridinium
1
0.51
1.02
- In the Orbitals section of the Results box count over the appropriate number of orbitals (11 in the case of Dye 1), and record the results for the HOMO and LUMO in this molecule.
Notice that in this case the energy of the HOMO is \( \alpha + 0.2737 \beta \); and, the energy of the LUMO is \( \alpha -0.4390 \beta \).
NOTE: you can change the number of digits in the display from Preferences, under the File menu.
Close the Results Box.
- To get a visual representation of the pi-molecular orbitals generated in the Hückel program, simply click on one of the orbitals in the right-hand pane of the program window. The orbital will appear on top of the molecule in the main screen. The image below is of the HOMO.
You can turn off the orbital display from your molecule by clicking below the orbitals in the right-hand pane of the program window. Save a screen-shot image of the entire program window showing the structure and the orbital diagram.
- As you scan down the Hamiltonian matrix check to see that the Coulomb and Exchange integrals for heteroatoms and carbon-to-heteroatom bonds are as shown in the table below. If you need to see how your atoms are numbered in your structure to make it easier to find them in the Result table, put a check in the Numbering box in the left hand pane of the program window.