Full

Download as pdf or txt
Download as pdf or txt
You are on page 1of 194

CHEM310L - PHYSICAL

CHEMISTRY I LAB
MANUAL
CHEM310L
Physical Chemistry Laboratory
This text is disseminated via the Open Education Resource (OER) LibreTexts Project (https://LibreTexts.org) and like the hundreds
of other texts available within this powerful platform, it is freely available for reading, printing and "consuming." Most, but not all,
pages in the library have licenses that may allow individuals to make changes, save, and print this book. Carefully
consult the applicable license(s) before pursuing such effects.
Instructors can adopt existing LibreTexts texts or Remix them to quickly build course-specific resources to meet the needs of their
students. Unlike traditional textbooks, LibreTexts’ web based origins allow powerful integration of advanced features and new
technologies to support learning.

The LibreTexts mission is to unite students, faculty and scholars in a cooperative effort to develop an easy-to-use online platform
for the construction, customization, and dissemination of OER content to reduce the burdens of unreasonable textbook costs to our
students and society. The LibreTexts project is a multi-institutional collaborative venture to develop the next generation of open-
access texts to improve postsecondary education at all levels of higher learning by developing an Open Access Resource
environment. The project currently consists of 14 independently operating and interconnected libraries that are constantly being
optimized by students, faculty, and outside experts to supplant conventional paper-based books. These free textbook alternatives are
organized within a central environment that is both vertically (from advance to basic level) and horizontally (across different fields)
integrated.
The LibreTexts libraries are Powered by NICE CXOne and are supported by the Department of Education Open Textbook Pilot
Project, the UC Davis Office of the Provost, the UC Davis Library, the California State University Affordable Learning Solutions
Program, and Merlot. This material is based upon work supported by the National Science Foundation under Grant No. 1246120,
1525057, and 1413739.
Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not
necessarily reflect the views of the National Science Foundation nor the US Department of Education.
Have questions or comments? For information about adoptions or adaptions contact info@LibreTexts.org. More information on our
activities can be found via Facebook (https://facebook.com/Libretexts), Twitter (https://twitter.com/libretexts), or our blog
(http://Blog.Libretexts.org).
This text was compiled on 08/15/2023
TABLE OF CONTENTS
Licensing

1: Orientation to this course


1.1: Pre-lab orientation assignment
1.2: Introductory Details
1.3: Rubrics
1.4: Statements

2: Introduction to Matlab
2.1: Pre-lab Assignment
2.2: Matlab Practice (In Class Activity)
2.2.1: Practice using a custom script
2.2.2: FTIR spectrum of HCl
2.2.3: Emission Spectrum of Iodine Vapor
2.3: How to Access Matlab
2.3.1: Invoking Matlab From the Lab Computers
2.3.2: Install FastX and access software with your own computer
2.3.3: Access from off-campus (The Duke VPN)
2.3.4: Install Matlab on your own computer
2.4: Using Matlab (Supplemental)
2.4.1: Entering Matricies
2.4.2: Matrix Computation
2.4.3: Vectors
2.4.4: Subscripting
2.4.5: Finding, Deleting, Adding or Replacing Data
2.4.6: Output Format
2.4.7: The Help Facility
2.4.8: Array Operations
2.4.8.1: Array Addition and Subtraction
2.4.8.2: Array Multiplication and Division
2.4.8.3: Array Powers
2.4.8.4: Math Functions
2.4.9: Graphics
2.4.9.1: Basic Plots
2.4.9.2: Line Types
2.4.9.3: Multiple Plots
2.4.9.4: Manual Axis Scaling
2.4.9.5: Overlaying Plots
2.4.10: Linear Least Squares Analysis
2.4.11: Last Line Editing and Recall
2.4.12: Saving Data
2.4.13: Deleting Matrices
2.4.14: Quitting Matlab
2.4.15: Transferring Files Between Your Computer and Your Duke NetID Account

1 https://chem.libretexts.org/@go/page/369891
3: The Treatment of Experimental Error
3.1: Pre-lab Assignment
3.2: Learning Objectives, Important Info
3.3: Characterizing Experimental Errors
3.4: Propagation of Uncertainty
3.5: Least Squares Linear Regression
3.5.1: Rejection of Discordant Data
3.5.2: Weighted Least Squares Analysis
3.5.3: References
3.6: Summary - General Guidelines For The Treatment of Experimental Errors
3.6.1: Error Analysis in a Nutshell
3.7: Homework Problems

4: Absorption Spectrum of Conjugated Dyes


4.1: Prelaboratory Assignment
4.2: Introduction
4.3: Experimental (Part I)
4.4: Hückel Calculations Using Hulis (Part II)
4.5: References
4.6: Appendix A - Use of the Agilent 100 Series UV-Vis Spectrophotometer
4.7: Appendix B - The Hückel Approximation

5: Molecular Spectroscopy of Iodine


5.1: Prelaboratory Assignment
5.2: Introduction
5.3: Potential Energy Curves
5.4: Part I - Absorption Spectrum of Iodine Vapor - Experimental
5.5: Part I - Absorption Spectrum of Iodine Vapor - Data Analysis
5.6: Part I - Absorption Spectrum of Iodine Vapor - Questions
5.7: Part II - Emission Spectrum of Iodine Vapor - Apparatus
5.8: Part II - Emission Spectrum of Iodine Vapor - Experimental
5.9: Part II - Emission Spectrum of Iodine Vapor - Data Analysis
5.10: Part II - Emission Spectrum of Iodine Vapor - Questions
5.11: Part II - Emission Spectrum of Iodine Vapor - References

6: Rotation-Vibration Spectrum of HCl AND DCl; Vibrational Spectrum of


SO2
6.1: Prelaboratory Exercise
6.2: Transition Intensities
6.3: Part I - Experimental
6.4: Part I - Data Analysis
6.5: Part I - References
6.6: Part II - Experimental
6.7: Part II - Data Analysis
6.8: Part II - References

2 https://chem.libretexts.org/@go/page/369891
7: Molecular Electronic Structure Calculations
7.1: Pre-Lab Assignment
7.2: Introduction to Electronic Structure Calculations
7.3: Accessing Spartan
7.4: Exercise 1 - Estimation of the IR Spectrum of HCl and DCl
7.5: Exercise 2 - Molecular Orbitals of Formaldehyde
7.6: Exercise 3 - Modeling Sulfur Dioxide - Comparing the Results of Different Basis Sets and Calculation Models
7.7: Exercise 4 - Modeling the Absorption Spectrum of Cyanine Dyes
7.8: Exercise 5 - Modeling an Intramolecular Rearrangement

8: Independent Project
8.1: Introduction
8.2: Literature
8.3: Computational

Index
Glossary

Detailed Licensing

3 https://chem.libretexts.org/@go/page/369891
Licensing
A detailed breakdown of this resource's licensing can be found in Back Matter/Detailed Licensing.

1 https://chem.libretexts.org/@go/page/417347
CHAPTER OVERVIEW

1: Orientation to this course


1.1: Pre-lab orientation assignment
1.2: Introductory Details
1.3: Rubrics
1.4: Statements

This page titled 1: Orientation to this course is shared under a not declared license and was authored, remixed, and/or curated by Kathryn Haas.

1
1.1: Pre-lab orientation assignment
Prior to our first laboratory meeting (it is the first week of classes!), please complete the following.
Plan at least three hours in the week prior to classes to complete this assignment.
1. Log into Duke's Sakai platform and familiarize yourself with this course's Sakai site.
2. Read the entire orientation module in this manual.
3. Find the syllabus, download it, and read the entire document carefully.
4. Complete the "Orientation Assignment" that is found at the end of the syllabus. Complete the assignment and submit it on
Gradescope prior to the first course meeting.
5. Review the Duke Chemistry Safety Manual.

Things to bring to the first lab meeting:


Please bring your working computer with you to our first meeting.

1.1: Pre-lab orientation assignment is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

1.1.1 https://chem.libretexts.org/@go/page/453019
1.2: Introductory Details
Course Information
This laboratory course compliments the undergraduate Physical Chemistry lecture course to give undergraduate students practical
laboratory experience with instruments and methods relevant to physical chemistry.
The online site for this course is on Sakai. Details about the course, (announcements, grading, instructor contact, assignment
submission) are available only by visiting our Sakai site. Links to this manual will be available on Sakai, although the manual
"lives" on LibreTexts. The manual is a work that is constantly revised, and you will find the latest version of all experiments on the
"live" LibreTexts site for this course.

Preparation for Laboratory


Pre-Lab Assignment:
You should prepare for lab meetings by reading the appropriate module in this manual, and completing the specified pre-laboratory
assignments. The pre-laboratory assignments are designed to prepare you to operate efficiently and safely during the experimental
module. You will not be allowed to conduct experiments if the pre-lab assignment is not completed before your lab meeting. If
there is an extenuating circumstance that prevents your preparation, please inform your instructors as soon as possible.
For a given experiment, you should complete the assigned reading and write out the specified pre-lab assignments. The pre-lab
assignment must be submitted as a single pdf document on Sakai at least 24 hours before the scheduled lab meeting. Late pre-
laboratory work will not be accepted and will receive a score of zero. Your TA will evaluate the work and return it no later than the
beginning of your lab section. Your preparedness will also be evaluated by your answers to oral questions pertaining to the
laboratory from the TA and other lab instructors. The pre-lab and preparation accounts for 10 points of each experiment requiring a
pre-lab.
For every experimental module there is a general pre-lab assignment. It is described here:
General Pre-lab assignment:
Write a brief introduction to the experiment and the experimental plan using the following outline:
I. Experiment Introduction: In your own words, briefly describe the underlying theory and explain the principle upon which the
experiment isbased. Keep it concise; lengthy discussions or derivations of equations are unnecessary.
a. Background (importance/relevance): A brief description and justification of the importance of the topic, with references to
source material. What is the topic and why is it interesting and important? This is meant to "hook" the reader so they want
to read about the topic.
b. Theory: A concise paragraph or more describing the general theory for the method/technique used; this section should
contain a reference to the appropriate pages of the textbook, or other source material from the literature.
c. Goals/Purpose: Brief statement of purpose, which should indicate what is analyzed and the technique used. Limit to three to
five sentences.
II. Experimental Plan: a short summary of the specific things that you would need to know in order to do the experiment without
having access to your lab manual, with the understanding that you will always have a TA available to teach you how to use the
instrument. The best way to approach this task is to thoroughly read the procedure for the experiment. While reading, take notes
on the specifics of the work. What chemicals will you use? If you have to make solutions, what solutions do you need to make
and of what concentration? How will you prepare them? What instrument will you use, and what specific instrument settings
will you need to input (sometimes you will find this information in the Appendices that discuss the use of the instruments)?
How many runs of each sample will you do; and over what wavelength, or temperature, etc. range? You do not need to write
down how to use the instrument; we will teach you that. Bottom Line: we want you to write down only that pertinent
information about samples, sample preparation, instrument and instrument parameters, numbers of runs, and other things, that
you will need to know in order to do your work.
Pre-lab assignments should be typed. If you prefer to write out equations by hand, they should be neat and inserted as an image in
the appropriate place. If you need advice for how to insert equations or images into your documents, please ask for help at least 48
hours before the assignment is due; an in-person (or zoom) meeting is most appropriate for this issue.

1.2.1 https://chem.libretexts.org/@go/page/408416
Always Wear Protective Clothing:
Please dress appropriately for a laboratory environment - this is considered an important part of preparation for the laboratory
period (see also Laboratory Safety below). You must wear close-toes shoes and clothing must cover your body from your upper
arms to your knees. You should also wear protective eyewear. If lab attire is not jiving with your wardrobe, please purchase a lab
coat and bring it with you so that you can cover up in lab. No student should be in lab without wearing protective eyewear and
appropriate protective clothing.

Laboratory Work and Reports


You will complete each experiment as a group, and there will be a dedicated TA to answer questions. You will record observations
using an electronic laboratory notebook (ELN) template. The ELN will contain a written record of your group's activities during
lab. Each student should collect a copy of the group data immediately after the lab meeting; Every lab computer has internet access,
and you may either email data to yourself or collect it using a usb drive. After the lab meeting, you will make the ELN your own by
filling in the Background, Theory, and Purpose statements that you wrote for the pre-lab assignment. Edit these statements as
necessary using the feedback from your TA. Then, complete the questions and data analysis as specified in the Lab manual and in
the ELN. Your completed ELN must be submitted for grading before the next schedule experiment (within 1 week). The
ELN accounts for 30 points for each experimental module.

Using the ELN


Note: before beginning any electronic laboratory notebook entry, select Highlight Changes from the Track Changes option under
the Tools menu in Microsoft Word. This way if you accidentally erase good data it can be retrieved afterwards. SAVE OFTEN!
There are three sections to each ELN: The "Observations and Data" section should be edited during lab while your group is
working together. The other sections must include your own work.
Pre-Lab section:

(Your own work) This is where you will copy/paste your own pre-lab work.
Observations and Data:
(Group data collection) Record all observations and primary data as a group. This section should not evaluate the data/observations
in any way. If data are being collected by hand it should be in clear, easy to read format (tables work well). Data collected at a
station away from your computer may be written on paper and immediately entered into the ELN file. For data that are collected by
a computer interface to your equipment, note the specific run and the name of the data file in case you have to retrieve the file later!
Take screen shots as necessary to capture important data to insert into the ELN. Write observations for each run, including "bad" or
aborted runs, in a narrative style and state where the actual data file can be found (including filename) when needed.
Post-Lab work:
(your own work)
Calculations of results and uncertainties: Unless otherwise instructed all calculations of results and uncertainties should be
included as part of the laboratory notebook. In this course doing so will require the typing of a lot of mathematics. Please remember
that you always have the option of writing your calculations out neatly on plain paper and scanning it to make a digital image that
can be incorporated into an electronic notebook entry. See instructions provided earlier in this section to learn how to do that.
Discussion of results: In this section you truly “think in the notebook”. Reflect on what you have done. Write this section in your
own words after the observations are completed. Include charts, graphs (see below), tables or rearranged data and your rambling
thoughts. Remember this is not a formal report. You want to be neat and complete, not necessarily a literary star. This section need
not be long, but it should contain all pertinent information.
Graphs: Often you will generate graphs from your primary data files. They must be represented in this section of the laboratory
notebook. It is not difficult to include digital copies of graphs within the electronic laboratory notebook, but you should remember
that such images can easily make the file sizes large and thus difficult to work with. Below is a procedure that will work to include
copies of your graphs in the notebook while minimizing the overall size of the resulting file. A reliable way to capture graphs/plots
to insert into your lab notebook is to use WinSnap for PC (or screenshot on Mac). Each lab computer has WinSnap installed. You
can take a screen shot by pressing SHIFT + COMMAND + S on the keyboard. This will allow you to select a region for capture

1.2.2 https://chem.libretexts.org/@go/page/408416
and the image will be saved to the clipboard. Press CONTROL + V to paste. You can also open the WinSnap program from the
desktop shortcut instead of using keyboard shortcuts. (SAVE OFTEN!)

Here are a few more pointers for lab reports


BE BRIEF AND CONCISE! This is one of the most difficult aspects of good writing to master. Read over your first draft and be
ruthless in chopping out superfluous words, phrases, or sentences. Learn to be economical in expressing yourself. The total length
of the lab report should be less than 15 pages (but please don't aim for 15! That's really long!). Short is sweet!

WATCH YOUR GRAMMAR! You may be a brilliant scientist in the lab, but your credibility will be diminished if you can't
communicate your results. Science has two parts: doing novel and reproducible experiments, and being able to tell other people
what you did!

WATCH YOUR SPELLING! You will not be marked down for the occasional typographical error, but a report full of spelling
errors will be penalized. Take the time to proofread the draft of your text at least once.

Some particular comments on writing a scientific report:


All tables should be numbered and titled. The table headings should include the units.
All figures should be numbered and titled. If the figure is a graph, the axes should be clearly labeled, including the units. If
more than one set of data is included on the graph, use different symbols and a legend to identify them.
Tables and figures should be presented on different pages, not included with the text. More than one table or figure can be
placed on a page, however. Be sure to refer to any table or figure in the text, e.g. "The current and potential values from which
E and n were calculated are reported in Table 2, and plotted in Figure 3."
All equations presented in the text should be numbered and all symbols defined, as shown in the example below:
"The diffusion-limited current, i , is related to the bulk analyte concentration, c , by the following expression:
d a

1/2 1/2
id = nF S Da ca /(πΔt) (1.2.1)

where n is the number of electrons transferred, F is the Faraday constant, S is the electrode area, Da is the diffusion
coefficient of the analyte ion, and t is the time during which the current is sampled."
Remember, the all-important criteria in determining the quality of the report are that it be clear, concise, and provide enough
information so that anyone else with your level of experience could repeat your work without consulting you.
Students without experience in writing reports are sometimes concerned about the length of the report and the time they take to
write it. The first report may take longer, but with experience you will find that a good report can be written in a few hours. The
most important aspect is that it presents the data and their analysis clearly and concisely. This can usually be done in 4 or 5 pages.
If you have access to a word processor, you should use it to prepare reports. Editing is then very easy, and a neat copy is easily
produced.

Work well with your team of instructors


There is a team of instructors to help you, including the course instructor (a Ph.D. chemist) and severalTeaching Assistants (TA's)
who are graduate students earning their PhD's in chemistry. One or more TA's will be in the lab with you each day, and one of them
will be dedicated to helping you and at least one of your labmates to run a given experiment. That dedicated TA will be present
during the laboratory to assist you, offer tips, and advice.
You are strongly advised to consult with your TA about a given report before writing it. This consultation should happen during
scheduled lab time when possible - do not leave lab early just because you have finished the experiment and cleaned up. Use lab
time wisely to work on your data and consult with your TA so they can help you during your time with them. The TA can show you
the important points to be emphasized in evaluating your data and will discuss the grading scheme.
You may also email your TA at least 48 hours before your lab report is due to ask questions about an assignment; do NOT wait till
the last moment to do this because you will not get help in time. Please be polite in your email correspondences. Please also
recognize that questions about chemistry are often difficult to answer in an email, and a meeting may be much more effective.
Please take advantage of lab time and office hours for this reason.

1.2.3 https://chem.libretexts.org/@go/page/408416
The office hours schedule with your team of instructors is posted is posted on Sakai in the course syllabus.
A list of TAs with their experiment assignments and office hours is posted on Sakai in the course syllabus.

LABORATORY SAFETY
The educational objectives of this course include learning about the hazards of materials, equipment and procedures in the chemical
laboratory. This is very important because it determines the level of risk to which you and other students are exposed by your
actions and the actions of others while you are in the laboratory. Through this type of training, you will gain the knowledge and
experience to determine what precautions are necessary in order to limit risk to yourself and others from chemical hazards. We
hope that you will make a serious effort at all times to ensure that your work in the laboratory is conducted as safely as possible.
The following rules have been developed with your safety in mind:
1. SAFETY GLASSES MUST BE WORN AT ALL TIMES IN THE LABORATORY-NO EXCEPTIONS. Contact lenses
may be worn in the laboratory. Of course, it is understood that you will be wearing safety glasses too.
2. A copy of the Duke Chemistry Department Safety Manual is available on the Chemistry web site. All students should read this
manual prior to beginning laboratory work.
3. No unauthorized experimentation will be allowed in the laboratory.
4. Proper shoes are required in the laboratory - no bare feet, sandals or open-toe footwear.
5. Never eat (including chewing gum), drink, or smoke in the laboratory.
6. To minimize skin exposure to chemical hazards, skin must be covered from neck to knee. Long sleeved shirt and pants should
be worn. We know it can get hot in Durham, so bring extra clothes with you to lab on hot days. Clothing that are not appropriate
for lab include shorts, mini-skirts, crop tops, low cut shirts, etc.
7. Long hair should always be tied back while in the laboratory.
8. Safety Data Sheets (SDS) and other safety notes for every chemical to be used in this laboratory course are kept in a ring binder
in the lab and also online. Students should read these notes prior to beginning each experiment.
9. Failure to observe safety regulations will result in substantial deduction in the lab grade and may result in ejection from the lab
session. Persistent violation of safety regulations will result in permanent expulsion from the lab course.
10. Disposal of Solvents and Chemicals: Never pour solvents or chemicals down a drain. Departmental waste disposal procedures
are described in the safety manual. However, you should ALWAYS consult with a TA before attempting to dispose of any
material.
11. Precautions for individual experiments will be discussed with you by the TA's and are described in the lab manual.

This page titled 1.2: Introductory Details is shared under a not declared license and was authored, remixed, and/or curated by Kathryn Haas.
Introductory Details by W. R. Fawcett, John Berg, P. B. Kelley, Carlito B. Lebrilla, Gang-yu Liu, Delmar Larsen, Paul Hrvatin, David
Goodin, and Brooke McMahon is licensed CC BY-NC-SA 4.0.

1.2.4 https://chem.libretexts.org/@go/page/408416
1.3: Rubrics
Each of the modules and their components are allotted points based on the information below. Teaching assistants will modify the
point breakdown as necessary. These rubrics, and all changes to these rubrics, will be applied equally across all student's work in a
given semester.

Overall points
35 pts Matlab participation and completion.
35 pts Treatment of experimental error problem set: correct completion of all problems.
160 pts Four laboratory "wet-bench" experiments (40 pts each including pre-lab and ELN)
30 pts Spartan calculations
50 pts Written report

Points for each module


MatLab
Introduction to MatLab

Preparation for meeting by installing MatLab on computer, completing


10 points Matlab basics tutorials, bringing working computer with MatLab
installed to lab meeting.
Participation in the MatLab Tutorial and correct completion of
25 points
assignment.

Error
Treatment of Experimental Error

5 points Preparation for meeting by submission of attempt of all problems.

5 points Participation in the Error study hall session

25 points Correct response on all questions

Dyes with Huckle


Absorption spectra of conjugated Dyes

10 points preparation for lab (pre-lab assignment)

Spectra of the following: Validation of instrument, 4 dyes.


Data: Uncertainty in peak wavelengths, comparison with literature
3 points
values
(See Q1-3)

2 points Calculation of experimental Energy in eV from Λ m ax (See Q4)

Selection of Γ using Matlab "dye" script (3 pts) and explanation of how


1 points
dye script works (1 pt). (See Q5)
Calculate Λ from particle in the box model and compare with
3 points max

experimental value. (See Q6-7)

3 points Discussion of results. (See Q8)

1 points Calculation of Λ max and Energy in eV for octatetraene. (See Q9)

1.3.1 https://chem.libretexts.org/@go/page/408417
Absorption spectra of conjugated Dyes

HMO calculations for each dye; correct structures; # of π elecrons and


3 points
connectivity matrix

2 points Δ E (beta units) for each dye

2 points MO diagram for each dye

3 points Plot HMO ΔE vs observed ΔE , determine beta and uncertainty

2 points Calculate energy difference in kL/mole, eV, and nm

1 point completed table

Discussion: Which theoretical model gives the best fit to the


experimental data, and why? Are the results for each dye equally
3 points
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
octatetraene using the HMO model program. In calculation use your
1 points
value of β to convert to electron volts. What color does the molecule
absorb from white light? What color does it then appear?

Iodine Spectroscopy (Electronic and Vibrational Coupling)


Students will participate in only one calorimetry experiment.

Iodine spectroscopy

10 points preparation for lab (pre-lab assignment)

PART I: Absorption spectrum


Validation spectrum and calibrated absorption spectrum - peaks
1 point
assigned

1 point Plot delta G's vs (v'+1) with linear regression

2 point Calculate w'e, X'e, D'e and a'

1 point Comparison with literature

2 point Uncertainty in w'e, X'e, and D'e

1 point Q1: Determine value of v00

1 point Q2: Why do the vibrational spacings get narrower as V(R) increases?

1 point Q3: Why is the Morse potential asymmetric?

1 point Q4: Why are some transitions called "Hot bands"?

1 point Q5: How would you calculate X"e and D"e using the hot bands?

1 point Q6: What are the units of the Morse parameter a?

Part II: Emission Spectrum

1 point Emission spectrum - peaks assigned

1 point Plot deltaG" vs (v"+1) w/linear regression

2 point Calculate w"e, X"e, and D"e

2 point Uncertainty in w"e, X"e, and D"e

2 point Calculate D"e using equation 8/ comparison with D'e

1 point Calculate a"

1.3.2 https://chem.libretexts.org/@go/page/408417
Iodine spectroscopy

1 point Plot of V"(R) and V'(R) (morse in MATLAB)

2 point Estimate R"e - R'e and R'e

1 point Plot both V"(R) and V'(R) on the same graph (morse 2 in MATLAB)

1 point Comparison with literature

Q1: Why are the vibrational spacings different in the ground and excited
1 point
states?
Q2: Why is the shape of the Morse potential different in the ground and
1 point
excited states?
Q3: How its it possible to observe emission peaks at higher energies
1 point
than the pump energy?

FTIR
Rotovibrational spectroscopy of HCl/DCl using FTIR

10 points Preparation for meeting by submission of attempt of all problems.

17 Points for Part I Part I: HCL/DCl

1 point Spectra of HCl & DCl - assignments labeled; list of peaks (see Step 1)

2 points Plot of ν̃ (m) vs m for 4 isotopes (see Steps 2 & 3)

5 points Be , αe , ν̃ o , Ie , re , and k for each isotope (see Steps 4-8)

Comparison with literature value (1 points); parts a-c (3 points) (see


1 points
Step/Q 9)
∗ ∗
ν̃ o Be
2 points Ratios and for H, D, and 35
Cl and 37
Cl (see Step 10)
ν̃ o Be

Explanation of why the 35Cl /37Cl isotope effect so much larger in DCl
1 points
than it is in HCl (see Step/Q 11)
Application of Boltzmann populations of the ground-state levels to
1 points
explain relative band intensities (see Step/Q 12)

1 points Determination of relative abundance (see Step/Q 13)

2 points Uncertainty in B e , αe , ν̃ o , Ie , re , and k for one isotope (see Step/Q 14)

Explanation of why we purge the FTIR spectrometer, and why we are


1 point
careful to purge background and sample similarly (see Step/Q 15)

13 points for Part II Part II SO2

1 point Q1 Spec. of SO , vibrations labeled


2

2 points Assign Bands

1 point Compare to literature

2 points Q2 other bands

2 points Q3 k and
1

2
l

2 points Q4 C (vib) at 298 K and 500 K


v

1 point Q5 Comp. of Spec. and exptl C v$

2 points Q6 Uncertainty in k and 1



2
(Q7)
l

1.3.3 https://chem.libretexts.org/@go/page/408417
Computational
Electronic Structure Calculations with Spartan

2 point Exercise I: Table summarizing HCl and DCl calculations

Exercise I: Comparison of computational results to experiment and


1 point
literature data
Exercise I: Conclusions about the accuracy of the Hartree-Fock
1 point
calculation

1 point Exercise I: Analysis of bond lengths for HCl and DCl

2 point Exercise II: C-O and C-H bond lengths; CHC and HCO bond angles

1 points Exercise II: Comparison with literature

Exercise II: Calculation of energy difference between HOMO and


2 points
LUMO
Exercise II: Images and descriptions of HOMO, LUMO, and other
2 points
orbitals

3 points Exercise III: Completed table

2 points Exercise III: Comparison of computational results to literature

1 point Exercise III: Evaluation of different theories/basis sets

2 points Exercise IV: Completed table

Exercise IV: Comparison of computational λ to experimental and


2 points max

theoretical values

3 points Exercise IV: Questions from Step 9

Exercise V: Equilibrium geometry energies of initial, transition, and


2 points
final states

1 points Exercise V: Reaction profile plots

2 points Exercise V: Calculation of activation energy and ΔH

This page titled 1.3: Rubrics is shared under a not declared license and was authored, remixed, and/or curated by Kathryn Haas.

1.3.4 https://chem.libretexts.org/@go/page/408417
1.4: Statements

Promoting Diversity, Inclusion, and Avoiding Harassment


All members of this course will foster an environment or listening and considering new ideas from a diverse group, with respect for
all participants without regard to gender, race, ethnicity, color, sexual orientation, political orientation, age, religion, national
origin, disability, or any other aspect of how we identify ourselves other than as fellow scientists.
There is zero tolerance in this course for any form of harassment or hate speech, during class meetings or outside of scheduled
class time, which could include verbal or physical conduct that has the purpose or effect of interfering with anyone else’s
participation or performance, or of creating an intimidating, hostile, or offensive environment; any such harassment will result in
dismissal from the course and will be reported to the college.
Duke University's Institutional Statement of Commitment to Diversity and Inclusion is found here (click):
"Duke aspires to create a community built on collaboration, innovation, creativity, and belonging. Our collective success
depends on the robust exchange of ideas—an exchange that is best when the rich diversity of our perspectives, backgrounds,
and experiences flourishes. To achieve this exchange, it is essential that all members of the community feel secure and
welcome, that the contributions of all individuals are respected, and that all voices are heard. All members of our community
have a responsibility to uphold these values."

This page titled 1.4: Statements is shared under a CC BY-SA license and was authored, remixed, and/or curated by Kathryn Haas.

1.4.1 https://chem.libretexts.org/@go/page/408418
CHAPTER OVERVIEW

2: Introduction to Matlab
 Learning Objectives

After completing the activities in this module, you will be able to


Access and open MatLab
Enter data sets on Matlab in the form of Matrices
Perform simple mathematical operations using Matlab
Plot data using Matlab
Run a custom scripts using Matlab

Matlab is an interactive software package for scientific and engineering numeric computation. The name Matlab stands for matrix
laboratory. Matlab works using essentially only one kind of object, a rectangular numerical matrix. In some situations, special
meaning is attached to 1-by-1 matrices, which are scalars, and to matrices with only one row or column, which are vectors.
Problem solutions are expressed in Matlab almost exactly as they are written mathematically - with no programming required.
In this laboratory course you will use Matlab to perform data analysis. The purpose of this activity is to review the basics of using
Matlab and to show you a little of the power and graphics capability of this package.
2.1: Pre-lab Assignment
2.2: Matlab Practice (In Class Activity)
2.2.1: Practice using a custom script
2.2.2: FTIR spectrum of HCl
2.2.3: Emission Spectrum of Iodine Vapor
2.3: How to Access Matlab
2.3.1: Invoking Matlab From the Lab Computers
2.3.2: Install FastX and access software with your own computer
2.3.3: Access from off-campus (The Duke VPN)
2.3.4: Install Matlab on your own computer
2.4: Using Matlab (Supplemental)
2.4.1: Entering Matricies
2.4.2: Matrix Computation
2.4.3: Vectors
2.4.4: Subscripting
2.4.5: Finding, Deleting, Adding or Replacing Data
2.4.6: Output Format
2.4.7: The Help Facility
2.4.8: Array Operations
2.4.8.1: Array Addition and Subtraction
2.4.8.2: Array Multiplication and Division
2.4.8.3: Array Powers
2.4.8.4: Math Functions
2.4.9: Graphics
2.4.9.1: Basic Plots
2.4.9.2: Line Types
2.4.9.3: Multiple Plots

1
2.4.9.4: Manual Axis Scaling
2.4.9.5: Overlaying Plots
2.4.10: Linear Least Squares Analysis
2.4.11: Last Line Editing and Recall
2.4.12: Saving Data
2.4.13: Deleting Matrices
2.4.14: Quitting Matlab
2.4.15: Transferring Files Between Your Computer and Your Duke NetID Account

2: Introduction to Matlab is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2
2.1: Pre-lab Assignment
 Pre-Lab assignment for MatLab

To prepare for the laboratory session on MatLab, please do the following prior to the laboratory period:
Download and install MatLab on your personal laptop computer.
Install the custom scripts and binary files (Duke students click here).
Bring your computer, with MatLab installed, to the lab meeting.
Proceed through the first two Matlab Getting Started tutorials that are provided when you open the software (Desktop
Basics and Matrices and Arrays).
We will meet in a classroom, thus lab attire is not necessary for this session.

This week is an exception in that there is no written work to submit prior to the meeting. However, for nearly the entirety of
the semester, you should expect to submit a pre-lab assignment 24-hours prior to the start of each new module. The general
instructions for preparation for laboratory are copied below for reference.

Reminder about general expectation for pre-lab preparation


Pre-Lab Assignment:
You should prepare for lab meetings by reading the appropriate module in this manual, and completing the specified pre-laboratory
assignments. The pre-laboratory assignments are designed to prepare you to operate efficiently and safely during the experimental
module. You will not be allowed to conduct experiments if the pre-lab assignment is not completed before your lab meeting. If
there is an extenuating circumstance that prevents your preparation, please inform your instructors as soon as possible.
For a given experiment, you should complete the assigned reading and write out the specified pre-lab assignments. The pre-lab
assignment must be submitted as a single pdf document on Sakai at least 24 hours before the scheduled lab meeting. Late pre-
laboratory work will not be accepted and will receive a score of zero. Your TA will evaluate the work and return it no later than the
beginning of your lab section. Your preparedness will also be evaluated by your answers to oral questions pertaining to the
laboratory from the TA and other lab instructors. The pre-lab and preparation accounts for 10 points of each experiment requiring a
pre-lab.
For every experimental module there is a general pre-lab assignment. It is described here:
General Pre-lab assignment:
Write a brief introduction to the experiment and the experimental plan using the following outline:
I. Experiment Introduction: In your own words, briefly describe the underlying theory and explain the principle upon which the
experiment isbased. Keep it concise; lengthy discussions or derivations of equations are unnecessary.
a. Background (importance/relevance): A brief description and justification of the importance of the topic, with references to
source material. What is the topic and why is it interesting and important? This is meant to "hook" the reader so they want
to read about the topic.
b. Theory: A concise paragraph or more describing the general theory for the method/technique used; this section should
contain a reference to the appropriate pages of the textbook, or other source material from the literature.
c. Goals/Purpose: Brief statement of purpose, which should indicate what is analyzed and the technique used. Limit to three to
five sentences.
II. Experimental Plan: a short summary of the specific things that you would need to know in order to do the experiment without
having access to your lab manual, with the understanding that you will always have a TA available to teach you how to use the
instrument. The best way to approach this task is to thoroughly read the procedure for the experiment. While reading, take notes
on the specifics of the work. What chemicals will you use? If you have to make solutions, what solutions do you need to make
and of what concentration? How will you prepare them? What instrument will you use, and what specific instrument settings
will you need to input (sometimes you will find this information in the Appendices that discuss the use of the instruments)?
How many runs of each sample will you do; and over what wavelength, or temperature, etc. range? You do not need to write
down how to use the instrument; we will teach you that. Bottom Line: we want you to write down only that pertinent

2.1.1 https://chem.libretexts.org/@go/page/453023
information about samples, sample preparation, instrument and instrument parameters, numbers of runs, and other things, that
you will need to know in order to do your work.
Pre-lab assignments should be typed. If you prefer to write out equations by hand, they should be neat and inserted as an image in
the appropriate place. If you need advice for how to insert equations or images into your documents, please ask for help at least 48
hours before the assignment is due; an in-person (or zoom) meeting is most appropriate for this issue.

Always Wear Protective Clothing:


Please dress appropriately for a laboratory environment - this is considered an important part of preparation for the laboratory
period (see also Laboratory Safety below). You must wear close-toes shoes and clothing must cover your body from your upper
arms to your knees. You should also wear protective eyewear. If lab attire is not jiving with your wardrobe, please purchase a lab
coat and bring it with you so that you can cover up in lab. No student should be in lab without wearing protective eyewear and
appropriate protective clothing.

2.1: Pre-lab Assignment is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.
Introductory Details by W. R. Fawcett, John Berg, P. B. Kelley, Carlito B. Lebrilla, Gang-yu Liu, Delmar Larsen, Paul Hrvatin, David
Goodin, and Brooke McMahon is licensed CC BY-NC-SA 4.0.

2.1.2 https://chem.libretexts.org/@go/page/453023
SECTION OVERVIEW

2.2: Matlab Practice (In Class Activity)


 Note

Before attempting these exercises, please make sure you have access to a working copy of the MatLab software and you have
the custom scripts and binary files installed in the matlab directory. If you are using FastX to access MatLab on campus, you
are good to go. Otherwise, you should install MatLab and the custom files on your computer before proceeding.

This module walks you through the steps for installing and using a custom script file. Then, there are two excercises which are
intended to be completed during our lab meeting. These in-class Matlab practice activities are are designed to help you through the
data analysis for two of the experiments later in this lab semester. You are essentially building a template that you can use for data
analysis later. All of the activities below can be completed during the first laboratory meeting, and may be performed individually
or with a partner. You should complete both during the lab meeting, and your work should be inspected by your teaching assistant
(TA) before you leave for the day. If your TA determines that you have done the activities correctly, you will receive full credit for
this assignment.

2.2.1: Practice using a custom script

2.2.2: FTIR spectrum of HCl

2.2.3: Emission Spectrum of Iodine Vapor

2.2: Matlab Practice (In Class Activity) is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.2.1 https://chem.libretexts.org/@go/page/370515
2.2.1: Practice using a custom script
Install a custom script
The custom scripts (files with .m extension) and binary files (files with .mat extension) for this course are available here (click to
go to OneDrive) using your Duke NetID. If you are using the Chemistry Department bohr server to access Matlab, the scripts are
already installed for you.
If you are using your own computer, the scripts and binary files should be added to the correct folder on your computer. Follow the
instructions in this video to install the lsq script, which performs least squares analysis.

Matlab Install custom scripts

Practice using a custom script and a binary file


This is an exercise to give you some practice in using Matlab to analyze data provided in a binary file, and plot the results using a
custom script for least squares analysis. Before you start, make sure the custom files lsq.m and expt-1.mat are in your working
MatLab directory (see video above). The binary file "expt-1" simply contains pre-defined 80 x 1 vector matrixes "x" and "y". This
is just to save you the hassle of typing in the values for each vector.
Follow the steps below (there is a video below to walk you through these steps):
1. Open the MatLab program and be sure that you are working in the directory containing the lsq script and expt-1 binary file.
2. Type the following: load expt-1
Your workspace now contains two variables, x and y which obey the following relation:
−C x
y = Ae (2.2.1.1)

3. Plot y (ordinate) versus x (abscissa). Using Matlab, recalculate your data in linear form (ln y versus x) and plot the result. Using
least squares analysis of the linearized data, calculate the values of A and C and report the standard deviation in these values.
Using your calculated values of A and C, use Matlab to draw (a) the least-squares fitted line through your linear plot and (b) the
best fit curve through your exponential plot. (NOTE: data should appear as points using the symbol 'o' and least squares fitted
lines should be solid lines.)

Expectation for a full-credit plot:


Data is plotted in scatter format.
Axes are correctly labeled, with correct units.
When appropriate, a curve fit to the data is shown
equation of the line is shown on the plot, but not obscuring the data
uncertainty in slope and intercept provided
correlation coefficient (R) provided
Approrpate units are given for every value and every axis label.
Plot is centered on graph.

2.2.1.1 https://chem.libretexts.org/@go/page/408552
Text is easily readable (usually this means it is large enough to read easily.

4. You should show the plots described in (3) together with the values of A and C to your instructor. Each graph should have a
title, and the x- and y-axes should be labeled.

Matlab Using a Custom Script and Plotti…


Plotti…

2.2.1: Practice using a custom script is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.2.1.2 https://chem.libretexts.org/@go/page/408552
2.2.2: FTIR spectrum of HCl
Using MatLab, open the "FTIR_StudentGuide...mlx" file, and follow the instructions in the document.

FTIR Spectrum of a mixture of HCl


35
and HCl 37

Figure 2.2.2.1 : FTIR Spectrum of HCl. (CC-NC-BY; DUKE CHEM)

2.2.2.1 https://chem.libretexts.org/@go/page/370517
Table of Peak Energies in wavenumbers:
Table 2.2.2.1 : ν~ (m) ( cm
−1
) for HCl from m = 4 to -4
~
ν (m) (cm
−1
) m isotope HCl 35
or HCl
37

2963.01

2960.84

2944.57

2942.52

2925.65

2923.48

2906.00

2903.95

2864.78

2862.85

2843.33

2841.40

2821.27

2819.34

2798.61

2796.80

1. Determine from Figure 2.2.2.1 which peaks arise from H35Cl and which peaks arise from H37Cl and give each pair of peaks an
appropriate m value.
2. Enter the peak positions ν~ (m) for H35Cl and H37Cl as two separate vectors into Matlab.
3. Create a corresponding vector of m values from +4 to –4.
4. For each isotope, use Matlab to calculate the separation between adjacent peaks Δν~ (m)
~ ~
= ν (m+1) − ν (m) . NOTE:
*Δν~ (m) cannot be calculated for m = –1 since there is no peak for m = 0.
5. For each isotope, do a separate plot of Δν~ (m)
(cm
−1
) (y-axis) versus m (x-axis) with each data point indicated by a small circle.
CAREFUL!: When plotting an x-y graph, the vectors x and y must be the same length. If Δν~ (m)
cannot be calculated for a
particular m value (see step (4)), then that value of m must be deleted from the vector of m values.
6. Perform a linear least-squares regression analysis and determine the slope, y-intercept, and standard deviation in the slope and
intercept.
7. For each isotope, show your TA your plot showing the raw data and the least-squares line.

2.2.2: FTIR spectrum of HCl is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.2.2.2 https://chem.libretexts.org/@go/page/370517
2.2.3: Emission Spectrum of Iodine Vapor
This semester, you use a laser spectrometer to investigate the Molecular Spectroscopy of Iodine Vapor. In this assignment you will
plot the electronic potential energy curves (Morse potentials) for the iodine molecule. A discussion of the terms used in this
assignment can be found in your textbook (D. A. McQuarrie, and J. D. Simon "Physical Chemistry: A Molecular Approach",
University Science Books, CA, 1997, Chapter 13, Section 13.1-13.7.
The functional form of the ground and excited state Morse potential wells for iodine molecule I2 are given by the expression
2
−a(R−Re )
V (R) = Te + De [1 − e ]

where V(R) is the potential energy of the system (y-axis), R is the interatomic separation (x-axis) between the iodine atoms, and Te,
and De, a and Re are constants. Literature values for these constants are given below.

Constant Ground State Excited State

Te 0 (cm 1
) (arbitrarily set) 15730 (cm 1
)

De 12244 (cm 1
) 4112 (cm 1
)

−1 −1
∘ ∘

a 1.87 (A) 2.00 (A)

∘ ∘
Re 2.67 (A) 2.97 (A)

1. Generate a vector R that varies between 2 and 6 Å, in increments of 0.01 Å.


2. For both ground and excited states, calculate V(R) for each value of R using the above expression and the constants from the table.
NOTE: These constants should be used in the units given above. No conversions are necessary.
3. On the same graph, plot V(R) (y-axis) versus R (x-axis), for both ground and excited states. Rescale the axes using the following
axis scales:
V(R): 0 – 22500 cm–1
R: 2 – 6 Å
4. Show the graph to your TA for approval.

2.2.3: Emission Spectrum of Iodine Vapor is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.2.3.1 https://chem.libretexts.org/@go/page/370546
SECTION OVERVIEW

2.3: How to Access Matlab


There are two options for accessing Matlab at Duke:
1. You can access the Matlab program using a lab computer. However, Matlab is not installed on the lab computer hard drive.
Rather, you can access the software that is installed on the chemistry department server (bohr) using an emulation software
(FastX).
Alternatively, FastX is already installed on all laboratory computers in our P-Chem lab, and you are permitted to use those
computers any time that the labs are open.
2. You can download and install software onto your own computer to access Matlab.
a. Install FastX onto your own computer to access the Matlab installed on the chemistry department server.
b. Install Matlab on your own computer.
The ability to access Matlab from your own computer will allow you to do your data analysis without having to come in to the lab
computers. There are two ways that you might accomplish this. The first is to access the version of Matlab on the chemistry server,
bohr, either by downloading and installing FastX (and the Duke VPN if you are off campus) onto your computer, or by using a
browser-based version of FastX. The second is to download and install the actual Matlab program on your own computer. Both
methods have advantages and disadvantages. Both the FastX and, if needed, VPN software packages are small and do not take up a
lot of space on your computer. Additionally running the server version of Matlab through bohr gives you instant access to our in-
house scripts that are used throughout the advanced lab classes. However, running a remote connection from your computer to bohr
can result in a slower response since the speed will be related to the speed of your network connection. Installing the Matlab
program on your own computer will likely give you a faster response, but the Matlab software will take up as much as six
gigabytes of space. Additionally, you will have to get copies of our in-house scripts and place them in the Matlab path so that they
can be used when needed. Instructions for both options are given in sections 1.1.2 and 1.1.3 of this Chapter.

2.3.1: Invoking Matlab From the Lab Computers

2.3.2: Install FastX and access software with your own computer

2.3.3: Access from off-campus (The Duke VPN)

2.3.4: Install Matlab on your own computer

2.3: How to Access Matlab is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.3.1 https://chem.libretexts.org/@go/page/370554
2.3.1: Invoking Matlab From the Lab Computers
INVOKING MATLAB FROM COMPUTERS IN THE LABS
1. To start Matlab on any of the computers in the lab:

1. Click on the FastX icon, on the main desktop screen, and then on the bohr.chem.duke.edu session. Enter your NetID
and password and click OK. Once logged in, click on the + in the upper right of the bohr.chem.duke.edu window.

Figure 2.3.1.1 : Connecting to Bohr through FastX. (CC-NC-BY-DUKE CHEM)


2. Click on one of the xterm choices from the new window that appears.

Figure 2.3.1.1 : Choose xterm. (CC-NC-BY-DUKE CHEM)


3. An x-terminal window will open. At the prompt type,

matlab

2.3.1: Invoking Matlab From the Lab Computers is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.3.1.1 https://chem.libretexts.org/@go/page/369914
2.3.2: Install FastX and access software with your own computer
This document describes two ways to access the server versions of Matlab, and later in the semester Spartan, from your own
computer. These programs are housed on a chemistry server called bohr (at bohr.chem.duke.edu). You can access this server either
by downloading and installing FastX; or, you can access bohr directly from your web browser. Both options are described below
Accessing Matlab and Spartan Through FastX Installed on Your Computer:
1. copies of both the Windows and Macintosh Installers are available in the Resources section of you Sakai site or can be obtained
from the Software section of the OIT web site, www.oit.duke.edu/software. Download and install the appropriate one for your
computer. Installation is quick and straightforward, just follow the prompts.

2. Once installed, Launch the program from the icon.

Figure 2.3.2.1 : Configuring FastX. (CC-NC-BY-DUKE CHEM)


3. Click on the + sign at the upper right of the main window. Select SSH from the menu that appears.
4. Complete the New SSH connection box as follows:

Figure 2.3.2.2 : bohr.chem.duke.edu. (CC-NC-BY-DUKE CHEM)


For Name, call it Bohr, or really anything you want. For Host, enter bohr.chem.duke.edu. The Port should be 22. Enter your NetID
username in the User box. The Path box should already be correct. Click on the Save button when done. The program will attempt
to connect to bohr.chem.duke.edu. If you get a message saying the system is not recognized by this computer (or something to that
effect), click Continue. When prompted, enter your NetID (if you didn’t add it to the connections setup box) and password.
5. Click on one of the xterm choices from the new window that appears.

2.3.2.1 https://chem.libretexts.org/@go/page/370644
Figure 2.3.2.3 : Choose xterm. (CC-NC-BY-DUKE CHEM)
6. An x-terminal window will open. At the prompt type,

matlab

Accessing Matlab and Spartan Through FastX Directly From a Web Browser:
1. Open a browser. We suggest using Chrome or Firefox, but likely any browser will work.
2. Navigate to https://bohr.chem.duke.edu:3300
When Prompted enter your NetID Username and Password

Figure 2.3.2.4 : FastX Login From Web Browser. (CC-NC-BY-DUKE CHEM)


3. Click on the blue + at the top left of the next screen.

Figure 2.3.2.5 : The New Session Button. (CC-NC-BY-DUKE CHEM)


4. Select xterm from the options on the next screen. Click Launch to start the session.

2.3.2.2 https://chem.libretexts.org/@go/page/370644
Figure 2.3.2.6 : Launch a new Xterm Session. (CC-NC-BY-DUKE CHEM)
When your xterminal window appears, type,

matlab

at the prompt.

2.3.2: Install FastX and access software with your own computer is shared under a not declared license and was authored, remixed, and/or curated
by LibreTexts.

2.3.2.3 https://chem.libretexts.org/@go/page/370644
2.3.3: Access from off-campus (The Duke VPN)
If you are planning to access Matlab, or later in the semester Spartan, from an off-campus location you will need to install and run
the Duke Virtual Private Network (VPN) software. Complete instructions for downloading, installing, and running the Duke VPN
software can be found here:
For Windows computers.
For Macintosh computers.
Remember to start the VPN before accessing Matlab (or Spartan) from an off-campus location.

2.3.3: Access from off-campus (The Duke VPN) is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.3.3.1 https://chem.libretexts.org/@go/page/370897
2.3.4: Install Matlab on your own computer
Another way to have access to Matlab from your own computer is to download and install your own copy of the software. To do so,
go to oit.duke.edu/software.

Figure 2.3.4.1 : The OIT Software Site. (CC-NC-BY-DUKE CHEM)


Notice that a link to Matlab & Simulink is right there on the main Software page. Click on that link, and then log in with your
NetID and password from the Login button in the upper right of the next screen. Once logged in you will see your choices along
with a statement of the license agreement Duke has with Mathworks.

Figure 2.3.4.2 : DUke Mathworks Account. (CC-NC-BY-DUKE CHEM)

2.3.4.1 https://chem.libretexts.org/@go/page/370891
Add the Duke Mathworks Account to your Shopping Cart. Then Proceed to checkout. Downloading and installing the software is
straightforward, follow the prompts. Please remember that this is a large program and will use approximately six gigabytes of
space. After installing the software you will need to activate the license. Finally, you will need to contact the Instructor for your
laboratory course to get copies of our in-house Matlab scripts that are used to help with your data analysis and calculations.

2.3.4: Install Matlab on your own computer is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.3.4.2 https://chem.libretexts.org/@go/page/370891
2.4: Using Matlab (Supplemental)

If you are a brand new Matlab user, or its been a long time since you've used it, the "hopw-to" sections of this manual may
seem opaque. A great way to start it to run through the "Getting Started" tutorials from Matlab.

We will introduce you to using Matlab to perform numeric calculations. If you wish, you may save a copy of your work for printing
later. This is done using the diary command, which captures a copy of all subsequent keyboard input and most of the resulting
output (with the exception of graphs) to be written to a named file. If the file already exists, the output is appended to the end of the
file. For example, the command,

diary exptA

will save your work in the file exptA. This file can be downloaded to your computer for later use as described in section 1.2.15.
When you have completed your work, the command,

diary off

will suspend the diary.


IMPORTANT NOTE:
The diary command merely "records" your on-screen work. It does not save your data. Data must be saved as described in
Section 1.2.12. Although diary files may be printed, they cannot be retrieved for viewing on the screen.

2.4: Using Matlab (Supplemental) is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.1 https://chem.libretexts.org/@go/page/369939
2.4.1: Entering Matricies
Matlab works with essentially only one kind of object, a rectangular numerical matrix. In some situations, special meaning is
attached to 1-by-1 matrices, which are scalars, and to matrices with only one row or column, which are vectors.
Matrices can be introduced into Matlab in several different ways. The easiest method for small matrices is to use an explicit list.
The explicit list is surrounded by square brackets, and uses the semicolon, ; , to indicate the ends of the rows. For example, to enter
a 3-by-3 matrix, we type:

A = [1 2 3; 4 5 6; 7 8 0]

This results in the output


A=
123
456
780
The matrix A is saved for later use. The individual elements in the input should be separated by commas or blanks and can be any
Matlab expression. For example, typing

x = [-1.3, 4/5, 4*atan(1)]

uses Matlab like a calculator and results in


x=
-1.3000 0.8000 3.1416
Matrices entered, horizontally, as shown above, can sometimes be long enough to require more than one line. Matlab must be told
when a wrap to another line is wanted. To do this, as you near the end of the line, add a space (spacebar) followed by three periods
( …). Then hit enter. You will be allowed to continue entering numbers into the same matrix on the new line. Alternatively, if
entering matrices horizontally is troublesome, you can always enter numbers vertically.
Typing:

y = [1
2
3
4
5]

results in:
y=
1
2
3
4
5

2.4.1: Entering Matricies is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.1.1 https://chem.libretexts.org/@go/page/369941
2.4.2: Matrix Computation
We will run through some basic matrix operations. Operations and commands in Matlab are intended to be natural, in a matrix
sense, not unlike how they might be indicated on paper. For example, our matrix A can be transposed with:

B = A'

which results in:


B=
147
258
360
Transposing x,

x = x'

results in:
x=
-1.3000
0.8000
3.1416
Matrix multiplication is indicated with:

C = A * B

producing:
C=
14 32 23
32 77 68
23 68 113

2.4.2: Matrix Computation is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.2.1 https://chem.libretexts.org/@go/page/369942
2.4.3: Vectors
Vectors with regularly spaced elements may be generated using the colon:operator.
The statement

D = 1:4

results in
D=
1234
Two colons may be used to create vectors with increments between the elements. Typing:

E = 0:0.1:0.5

results in
E=
0.0000 0.1000 0.2000 0.3000 0.4000 0.5000
Vectors are used often in Matlab because of the ease with which they allow data manipulation and function evaluation. When a
function is applied to a vector in Matlab, the function is applied to each element of the vector. For example, consider taking the sine
of the four elements in vector D above. Typing:

Y = sin(D)

results in
Y=
0.8415 0.9093 0.1411 –0.7568
A wide variety of data analysis functions are available that operate on vectors. Typical functions include cumulative sum,

cumsum(Y)

ans =
0.8415 1.7508 1.8919 1.1351
mean and standard deviation,

m = mean(Y), s = std(Y)

m=
0.2838
s=
0.7758
Other useful functions include:

max(x) returns the largest element in the vector x

min(x) returns the smallest element in the vector x

median(x) returns the median value of the vector x

2.4.3.1 https://chem.libretexts.org/@go/page/369943
sort(x) sorts the elements of x in ascending order

sum(x) returns the sum of the elements of x

2.4.3: Vectors is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.3.2 https://chem.libretexts.org/@go/page/369943
2.4.4: Subscripting
Enclosing their subscripts in parentheses may reference individual matrix elements. The statement A(m,n) refers to the mth row,
and nth column, of matrix A. For example, given a matrix A:
A=
123
456
789
the statement

A(3,3) = A(1,3) + A(2,2)

results in
A=
123
456
788
A subscript can be a vector. If X and V are vectors, then X(V) is [X(V(1)), X(V(2)), ..., X(V(n))]. For matrices, vector subscripts
allow access to contiguous and noncontiguous submatrices. For example, suppose that A is a 10-by-10 matrix. Then,

A(1:5,3)

specifies the 5-by-1 submatrix, or column vector, that consists of the first five elements in the third column of A. Similarly,

A(1:5,7:10)

is the 5-by-4 submatrix of elements from the first five rows and the last four columns.
Using the colon by itself in place of a subscript denotes all of the corresponding row or column. For example,

A(:,3)

is the third column and

A(1:5,:)

is the first five rows.


In general, if v and w are vectors with integer components, then

A(v,w)

is the matrix obtained by taking the elements of A with row subscripts from v and column subscripts from w.

2.4.4: Subscripting is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.4.1 https://chem.libretexts.org/@go/page/369983
2.4.5: Finding, Deleting, Adding or Replacing Data
The following examples will illustrate how to find, change, add or delete data in a vector. The same principles apply to
manipulation of data in matrix form.
(a) Finding Data:
Suppose we wish to find the elements of vector V which are greater than or equal to 7:

V = [-2 3 7 -26 48 10 -8]', i = find(V>=7)

results in
V=
-2
3
7
-26
48
10
-8
i=
3
5
6
The statement, i = find(V>=7), returns a vector, i, containing the index of the elements of V which are greater than or equal to 7.
The third, fifth and sixth elements of V are ≥ 7.
There are six relational operators that can be used in Matlab:
Relational Operators
< less than
<= less than or equal to
> greater than
>= greater than or equal to
== equal to
~= not equal to

(b) Deleting Data:


Specific elements in a matrix or vector can be deleted using the empty matrix (or vector), []. The statement

x = []

assigns a matrix of dimension zero-by-zero to x. To delete the elements of V which are greater than or equal to seven, use the
command

V(i) = []

2.4.5.1 https://chem.libretexts.org/@go/page/370473
which returns
V=
-2
3
-26
-8
The elements of V which have indices given by the elements of i, are removed by assigning them to an empty matrix. Single
elements may be removed in similar fashion. The command

V(3) = []

deletes the third element of V and returns


V=
-2
3
-8
An efficient way of removing rows and columns of a matrix is to assign them to an empty matrix. The command, A(:,3)=[], deletes
column 3 of matrix A.
(c) Adding Data:
Suppose we wish to combine the elements of vectors V and U:

U = [3:2:10]', W = [V; U]

This returns
U=
3
5
7
9
W=
-2
3
-8
3
5
7
9
If we wish to add a single data point (-17) to the end of vector W:

W = [W; -17]

This returns

2.4.5.2 https://chem.libretexts.org/@go/page/370473
W=
-2
3
-8
3
5
7
9
-17
(d) Replacing Data:
To replace the fourth element of W with -12, enter the following

W(4) = -12

This returns
W=
-2
3
-8
-12
5
7
9
-17

2.4.5: Finding, Deleting, Adding or Replacing Data is shared under a not declared license and was authored, remixed, and/or curated by
LibreTexts.

2.4.5.3 https://chem.libretexts.org/@go/page/370473
2.4.6: Output Format
The result of any Matlab assignment statement is displayed on the screen, as well as assigned to the specified variable or to ans if
no variable name is given. The numeric display format can be controlled using the format command. format affects only how
matrices are displayed, not how they are computed or saved (Matlab performs all computation in double precision).
If all the elements of a matrix are exact integers, the matrix is displayed in a format without any decimal points. For example,

x = [-1 0 1]

always results in
x=
-1 0 1
If at least one of the elements of a matrix is not an exact integer, there are several possible output formats. The default format,
called the short format, shows approximately 5 significant decimal digits. The other formats show more significant digits or use
scientific notation. As an example, suppose

x = [4/3 1.2345e-6]

The formats, and the resulting output for this vector, are:
format short
1.3333 0.0000
format short e
1.3333E+000 1.2345E-006
format long
1.33333333333333 0.00000123450000
format long e
1.333333333333333e+00 1.234500000000000e-06
For the long formats, the last significant digit may appear to be incorrect, but the output is actually an accurate decimal
representation of the binary number stored in the computer.
With the short and long formats, if the largest element of a matrix is larger than 1000 or smaller than 0.001, a common scale factor
is applied to the entire matrix when it is displayed. For example,

x = 1.e20*x

multiplies x by 1020 and results in the display


x=
1.0E+020 *
1.3333 0.0000
One final command, format compact, suppresses many of the line-feeds that appear between matrix displays and allows more
information to be shown on the screen.

2.4.6: Output Format is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.6.1 https://chem.libretexts.org/@go/page/370478
2.4.7: The Help Facility
A HELP facility is available, providing online information on most Matlab topics. To get a list of HELP topics, type

help

To get HELP on a specific topic, type help topic. For example,

help eig

provides HELP information on the use of the eigenvalue function,

help []

tells how to use brackets to enter matrices, and

help help

is self-referential, but works just fine.

2.4.7: The Help Facility is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.7.1 https://chem.libretexts.org/@go/page/370481
2.4.8: Array Operations
We use the term array operations to refer to element-by-element arithmetic operations, instead of the usual linear algebraic matrix
operations denoted by the symbols * / \ ^ '. Preceding an operator with a period . indicates an array or element-by-element
operation.

2.4.8: Array Operations is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.8.1 https://chem.libretexts.org/@go/page/370482
2.4.8.1: Array Addition and Subtraction
For addition and subtraction, the array operations and the matrix operations are the same, so + and - can be regarded as either
matrix or array operations.

2.4.8.1: Array Addition and Subtraction is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.8.1.1 https://chem.libretexts.org/@go/page/370485
2.4.8.2: Array Multiplication and Division
Array, or element-by-element, multiplication is denoted by .*. If A and B have the same dimensions, then A .* B denotes the array
whose elements are simply the products of the individual elements of A and B. For example, if
x = [1 2 3]; y = [4 5 6];
then,

z = x.*y

results in
z=
4 10 18
X = A./B is a solution to X.*B = A and X = A.\B is a solution to A.*X = B so, the expressions A./B and A.\B give the quotients of
the individual elements.

z = x.\y

results in
z=
4.0000 2.5000 2.0000
To obtain the reciprocal of the elements of x, type

w = x.\1

This results in
w=
1.0000 0.5000 0.3333

2.4.8.2: Array Multiplication and Division is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.8.2.1 https://chem.libretexts.org/@go/page/370486
2.4.8.3: Array Powers
Element-by-element powers are denoted by .^. Here are several examples, using the above vectors x and y. Typing

z = x.^y

results in
z=
1 32 729
The exponent can be a scalar.

z = x.^2

z=
149
Or, the base can be a scalar.

z = 2 .^x

z=
248

 Note
This last example illustrates one of Matlab's syntactic subtleties. Although it is difficult to see, the space between the digit 2
and the period . is important. If it were not there, the period would be interpreted as a decimal point associated with the 2. An
alternative to the space is to use parentheses, forcing the correct precedence.

2.4.8.3: Array Powers is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.8.3.1 https://chem.libretexts.org/@go/page/370487
2.4.8.4: Math Functions
Available math functions include the following:

MATH FUNCTIONS
Use the Help command to obtain further information on any of these functions

abs absolute value or complex magnitude

sqrt square root

real real part

imag imaginary part

round round to nearest integer

exp exponential base e

log natural logarithm

log10 log base 10

sin sine

cos cosine

tan tangent

asin arcsine

acos arccosine

atan arctangent

2.4.8.4: Math Functions is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.8.4.1 https://chem.libretexts.org/@go/page/370488
2.4.9: Graphics
One of the major advantages of Matlab is its graphics capability. The graphics commands are natural and easy-to-use, yet still
succinct, powerful, and extensible.
Scientific and engineering data are examined graphically in Matlab by using a "graph" command to create a plot on the screen.
There are seven different types of "graph" from which to choose:
GRAPH
plot linear X-Y plot
loglog loglog X-Y plot
semilogx semi-log X-Y plot (x-axis logarithmic)
semilogy semi-log X-Y plot (y-axis logarithmic)
polar polar plot
mesh 3-dimensional mesh surface
contour contour plot

The plot command creates linear X-Y plots. Once the plot command is understood, logarithmic and polar plots are created by
substituting the commands loglog, semilogx, semilogy, or polar for plot. All five commands are used the same way; they only
affect how the axis is scaled and the data are displayed.

2.4.9: Graphics is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.9.1 https://chem.libretexts.org/@go/page/370493
2.4.9.1: Basic Plots
If Y is a vector, plot(Y) produces a linear plot of the elements of Y versus the index of the elements of Y. For example, suppose we
would like to plot the numbers {0., .48, .84, 1., .91, .6, .14}. This is quickly accomplished by typing

Y = [0. .48 .84 1. .91 .6 .14];


plot(Y)

which results in a graph on the screen:

Figure 2.4.9.1.1 : A Basic Plot. (CC-NC-BY-DUKE CHEM)


Notice that the data are auto-scaled and x- and y-axes are drawn. Use X Label; Y Label; and, Title from the Insert Menu to add a
title and axes labels to your graph. Or, select Property Editor from the View menu. A Properties Editor box appears from which you
can enter axes labels, a title and etc.

Figure 2.4.9.1.2 : Property Editor. (CC-NC-BY-DUKE CHEM)

2.4.9.1: Basic Plots is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.9.1.1 https://chem.libretexts.org/@go/page/370494
2.4.9.2: Line Types
(b) Line types
The line types used on a graph may be controlled if the defaults are not satisfactory. Typically one does not connect data points
with a line. Point plots using various symbols may be selected directly in the plot command. For example,

plot(X,Y,'x')

draws a point plot using x-mark symbols while

plot(X1,Y1,':',X2,Y2,'+')

uses a dotted line for the first curve and the plus symbol + for the second curve. Create vectors X1, Y1, X2 and Y2 and try this!
Other line and point types are:

Line Type Point Type

solid - point .

dashed -- plus +

imag imaginary part

dotted : star *

dashdot -. circle o

x-mark x

Alternatively, the line type can be modified from within the Properties Editor window after the plot command is issued. Working

with the graph above, with the Properties Editor window open, first select the arrow tool . Then click once somewhere on the
actual data points. The data set will be highlighted and the Properties Editor window will change to display line properties. Under
the Line menu select no line, and choose a marker from the Marker menu.

Figure 2.4.9.2.1 : Line Editor.


Your graph will change accordingly. Finally, selecting More Properties from the Property Editor opens an additional options box.

2.4.9.2.1 https://chem.libretexts.org/@go/page/370495
Figure 2.4.9.2.2 : Property Inspector.

2.4.9.2: Line Types is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.9.2.2 https://chem.libretexts.org/@go/page/370495
2.4.9.3: Multiple Plots
The way to plot multiple lines on a single graph is to use plot with multiple arguments:

plot(X1,Y1,X2,Y2,...,Xn,Yn)

The variables X1,Y1, X2,Y2 etc. are pairs of vectors. Each x-y pair is graphed, generating multiple lines on the plot. Multiple
arguments have the benefit of allowing vectors of differing lengths to be displayed on the same graph. Each pair uses a different
linetype.
We can apply what we learned in the previous sections to plot simple functions. For example, to examine the behavior of the
functions sin(x) and x2/100 on the interval from -15 to 15:

x = -15:.05:15;
y = sin(x);
z = x.^2/100;
plot(x,y,x,z)

which results in:

Figure 2.4.9.3.1 : Multiple Plots.

2.4.9.3: Multiple Plots is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.9.3.1 https://chem.libretexts.org/@go/page/370496
2.4.9.4: Manual Axis Scaling
In certain situations, it may be desirable to override the automatic axis scaling feature of the plot command and manually select the
axis limits. Manual axis scaling is easily accomplished from the Properties Editor window mentioned above.
Alternatively, manual Axes Scaling can be accomplished from the Matlab workspace by creating a vector, V, with the following
configuration
V=[X-min, X-max, Y-min, Y-max]
Working from the example graph above, let’s say you wanted to scale the X axis from –5 to 5 and the Y axis from –1.5 to 1.5.
Enter the following:

V = [-5, 5, -1.5, 1.5]


axis(V)

Retype the plot command plot(x,y,x,z) to observe the rescaled graph. The command axis(V) sets the axis scaling to the prescribed
limits. If you reenter the command:

axis

autoscaling is resumed. Again, retype the command, plot(x,y,x,z), to observe the graph.

2.4.9.4: Manual Axis Scaling is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.9.4.1 https://chem.libretexts.org/@go/page/370501
2.4.9.5: Overlaying Plots
The command,

hold

holds the current graph on the screen. Subsequent plot commands will add to the graph, using the established axis limits and
retaining the previously plotted curves. hold remains in effect until this command is issued again (i.e. the command hold can be
"toggled" on and off).

2.4.9.5: Overlaying Plots is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.9.5.1 https://chem.libretexts.org/@go/page/370502
2.4.10: Linear Least Squares Analysis
To fit a straight line, y = mx + b, to a set of experimentally determined points
(x 1,y 1), (x 2,y 2), ..., (x n,y n) contained in the vectors [x] and [y], use Basic Fitting under the Tools menu in your graph window.
Try this using the following data:

x = [-3.8 -1.2 0.75 3.1 5.6 8];


y = [8.2 6 3.3 1.6 -2 -6];

Selecting Basic Fitting from the Tools menu opens up an options box. Select linear as the fit type and choose how many significant
figure you want in the calculation. A Linear Least-Squares line will be placed on the graph along with the equation for the line.

Figure 2.4.10.1 : Basic Linear Fit. (CC-NC-BY-DUKE CHEM)


NOTE: Often for our purposes in the Physical Chemistry Laboratory it is often important to know the uncertainties (or standard
deviation) in the calculated slope and intercept. To obtain calculated values of the slope, intercept, correlation coefficient along
with the uncertainties in slope and intercept, use the lsq command. Using the current x,y data as an example, typing:

lsq(x,y)

returns the following:

Figure 2.4.10.2 : The lsq Function. (CC-NC-BY-DUKE CHEM)


To perform a weighted linear least-squares analysis, you would use the command wtlsq(x,y). This requires that you enter a
weighting factor e.

2.4.10: Linear Least Squares Analysis is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.10.1 https://chem.libretexts.org/@go/page/370505
2.4.11: Last Line Editing and Recall
The arrow keys on the terminal can be used to edit mistyped commands or to recall previous command lines. For example, suppose
you enter

log(sqt(atan(2*(3+4))))

You have misspelled sqrt. The response from Matlab is the error message,
Undefined function ‘sqt’ for input arguments of type 'double'.
Instead of retyping the entire line, simply hit the Up-Arrow key. The incorrect line will be displayed again and you can move the
cursor over using the Left-Arrow key and insert the missing r.

log(sqrt(atan(2*(3+4))))

ans = 0.2026
The arrow keys on the keyboard work on copies of the previous input lines, which have been saved in a moderately sized input
buffer. Here is a brief description of their function:
Last-line Editing and Recall Keys
Up Arrow Recalls previous lines.
Left Arrow Move left one character.
Right Arrow Move right one character.
Ctrl-U Cancel current line.
Ctrl-E Move to end of line.
Ctrl-A Toggle between insert and overtype mode.
Delete Delete character at cursor.
Backspace Delete character left of cursor.

2.4.11: Last Line Editing and Recall is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.11.1 https://chem.libretexts.org/@go/page/370506
2.4.12: Saving Data
If you wish to save the variables in your workspace for later use, enter:

save filename

This saves all variables in the file filename.mat.. The next time you invoke Matlab, the workspace may be restored using:

load filename

We highly recommend that you save your data frequently as you work. This will prevent a major loss of data in the event of
a system failure.
If you wish to save graphics, work from within the graphics window. Choose Save from the File menu. A dialog box opens up.
Give your file a name, in the File Name box. You can select the file type from the Files of Type menu. Matlab graphics default to a
.fig file type, and if you intend to open the graph again in Matlab that is a good option. If you want a picture of the graph for use in
electronic lab notebooks (ELN’s) or in your formal written reports you might prefer the .jpg format. The file will be stored in the
specified directory (shown in the Save In menu) of your account associated with your NetID when you click on the Save button.

Figure 2.4.12.1 : Saving a Figure.


Saved .fig files can be opened anytime, from the Matlab prompt, by typing

open filename.fig

2.4.12: Saving Data is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.12.1 https://chem.libretexts.org/@go/page/370507
2.4.13: Deleting Matrices
If you wish to delete a variable from your workspace, enter:

clear name

where name is the name of the matrix or variable to be deleted.

 Note

Be careful, because executing,


clear
by itself, clears all variables from the workspace!

2.4.13: Deleting Matrices is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.13.1 https://chem.libretexts.org/@go/page/370508
2.4.14: Quitting Matlab
Before you exit Matlab type,

diary off

if you had a diary running.


Exit Matlab by typing,

exit

Type,

exit

to close any xterm windows you have open. Finally, please remember to exit out of FastX if you are accessing Matlab through a
FastX connection.

2.4.14: Quitting Matlab is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

2.4.14.1 https://chem.libretexts.org/@go/page/370509
2.4.15: Transferring Files Between Your Computer and Your Duke NetID Account
In several experiments this semester, you will use a lab computer for instrument control and data acquisition. There are times when
it is convenient to have your raw data files uploaded to bohr so that you might analyze the data using Matlab. Additionally, there
are times when files and figures (graphs) that are generated by matlab need to downloaded to your own computer. These operations
can be done in the following manner.
1. To upload files from one of the lab computers to your online account:

1.1. Click on the icon found on the computer desktop to start the WinSCP file transfer program. The File protocol
should be SFTP, and the port should be 22. Enter login.oit.duke.edu under Host Name. Enter your NetID (User name) and
password, and click on Login. Click Continue at the Authentication Banner.

Figure 2.4.15.1 : The Main WinSCP Window. (CC-NC-BY-DUKE CHEM)


1.2. You will see a dialog box, part of which is shown above. On the left-hand side of the screen is your local directory (i.e.,
the directory of the lab computer). It should default to the directory containing your data. If not, change the directory to the
correct one. You should see a list of those files found in the local directory. On the right-hand side you will see the contents
of your OIT account.
1.3. To transfer files from the local PC to your OIT account, first highlight the appropriate files you want transferred. Click
on the button at the top of the screen. An Upload dialog box will appear. Click OK to copy your files from the local
PC to your OIT account. Or, click on Transfer Settings if you want to check to see that the transfer is configured properly. (In
general, the Transfer mode should be Text and the File modification should be No change.) Then click on Copy. Check to be
sure all of your files are transferred and then exit your session by clicking the X at the upper right corner of the screen.

Figure 2.4.15.2 : The WinSCP Upload Screen. (CC-NC-BY-DUKE CHEM)


2. Downloading files from your OIT account to your own Windows computer or to one of the lab computers is accomplished in
the same manner, except that you should select the files you want to download, and execute the transfer from your OIT account to
your local computer.
3. Transferring files back and forth between your Macintosh computer and your OIT account can be accomplished in the
following manner.
3.1. From the Go menu (found when you are in your Finder), select Connect to Server.

2.4.15.1 https://chem.libretexts.org/@go/page/370893
Figure 2.4.15.3 : The Macintosh Go Menu. (CC-NC-BY-DUKE CHEM)
3.2. In the text box at the next screen enter,

smb://homedir.win.duke.edu/users/the first letter of your NetID/your full NetID

An example of a correct server address is shown below, where w is the first letter of the NetID, and woerner is the NetID.

Figure 2.4.15.4 : The Connect to Server Screen. (CC-NC-BY-DUKE CHEM)


3.3. Click the Connect button. You will be prompted for your NetID and Password.

Figure 2.4.15.5 : Login Screen. (CC-NC-BY-DUKE CHEM)


3.4. Once connected a new icon will appear on your desktop.

Figure 2.4.15.6 : The Connectred Server Icon. (CC-NC-BY-DUKE CHEM)


If you do not see the icon appear, Check the Finder Preferences,

2.4.15.2 https://chem.libretexts.org/@go/page/370893
Figure 2.4.15.7 : The Finder Menu. (CC-NC-BY-DUKE CHEM)
and make sure Connected Servers is checked.

Figure 2.4.15.8 : Finder Preferences. (CC-NC-BY-DUKE CHEM)


3.5. The connected server on your desktop acts like any other folder on your Mac. Click on it to open it and see the contents
of your online OIT account. Drag and drop to move files back and forth between your Mac and your OIT account.

2.4.15: Transferring Files Between Your Computer and Your Duke NetID Account is shared under a not declared license and was authored,
remixed, and/or curated by LibreTexts.

2.4.15.3 https://chem.libretexts.org/@go/page/370893
CHAPTER OVERVIEW

3: The Treatment of Experimental Error


3.1: Pre-lab Assignment
3.2: Learning Objectives, Important Info
3.3: Characterizing Experimental Errors
3.4: Propagation of Uncertainty
3.5: Least Squares Linear Regression
3.5.1: Rejection of Discordant Data
3.5.2: Weighted Least Squares Analysis
3.5.3: References
3.6: Summary - General Guidelines For The Treatment of Experimental Errors
3.6.1: Error Analysis in a Nutshell
3.7: Homework Problems

3: The Treatment of Experimental Error is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

1
3.1: Pre-lab Assignment
The ability to make reproducible observations is a critical skill for all laboratory workers. The observations may be qualitative
(result not expressed numerically) or quantitative (result expressed numerically). Wherever possible, scientists make quantitative
observations and analyze them for cause and effect relations using well-established techniques of mathematics.
The measurement of physical quantities such as mass, length, volume, etc. always involves some uncertainty. The overall quality of
a measurement depends on the accuracy of the measuring device and the skill of the user. Since no measurement is completely free
of error (i.e. there is always some uncertainty), it is common to repeat the measurement several times. Then the mean (average)
value of the replicate measurements is used as the best measure of the actual value and the standard deviation of the mean is used
as an estimate of the uncertainty in the mean value.
When a quantity is measured only once, it is not possible to have an average result or a standard deviation of the mean. It is still
possible, however, to obtain some information about the estimated uncertainty in the final result. This is accomplished by
considering the known uncertainties of the individual pieces of equipment used (balances, glassware, scientific instruments) and
propagating those uncertainties through the mathematical calculations used to arrive at the final answer. Often, but not always, the
uncertainties determined in this way are smaller than the standard deviation of the mean in a series of repeated trials. This is easy to
understand since the propagated error from a single trial only accounts for the known uncertainty of the equipment and not any
additional errors added by how well a person uses the equipment.
One way to gauge the improvement in your laboratory skills is to see the standard deviation of the mean from a series of repeated
measurements become equal to or less than the propagated uncertainty of any one trial. Then all of the error seen in the final result
is accounted for by the known uncertainties of the equipment and you could safely say that your use of that equipment did not add
additional uncertainty (error) to your results.
The purpose of this exercise is to give you and opportunity to learn and practice the types of statistical analyses and estimates of
uncertainty associated with your experimental results.

 Pre-lab assignment

In preparation for the lab meeting on Error Analysis, please do the following*:
Read D.A. Skoog, D.M. West, F.J. Holler and S.R. Crouch, Analytical Chemistry: An Introduction, 7th Ed., Saunders
College Publishing, New York, 2000, Chapters 5-7.
Read D.P. Shoemaker, C.W. Garland and J.W. Nibler, Experiments in Physical Chemistry, 6th Ed., McGraw Hill Co. Inc.,
New York, 1996, Chapter 2.
Attempt (does not mean you must complete) all seven problems in the last section of this module; work should be done in a
way that you can track your progress and submit your work in the form of a pdf. Submit all attempts in writing 24 hours
prior to your scheduled laboratory meeting.
Note that the textbooks referenced in this module are available in the Perkins library and also in pdf format (Duke students
click here).
We will meet in a classroom, thus lab attire is not necessary for this session.

*In addition to the general pre-lab assignment (copied below as a reminder), please do the above tasks. For this particular
assignment, an experimental plan is unnecessary. However, it is good practice to reflect on the purpose of this assignment, so
please do write a background/purpose statement. You should submit a pdf of your work on Sakai before the deadline (24 hours
prior to the lab meeting).

Reminder about general expectation for Pre-Lab preparation


Pre-Lab Assignment:
You should prepare for lab meetings by reading the appropriate module in this manual, and completing the specified pre-laboratory
assignments. The pre-laboratory assignments are designed to prepare you to operate efficiently and safely during the experimental

3.1.1 https://chem.libretexts.org/@go/page/453055
module. You will not be allowed to conduct experiments if the pre-lab assignment is not completed before your lab meeting. If
there is an extenuating circumstance that prevents your preparation, please inform your instructors as soon as possible.
For a given experiment, you should complete the assigned reading and write out the specified pre-lab assignments. The pre-lab
assignment must be submitted as a single pdf document on Sakai at least 24 hours before the scheduled lab meeting. Late pre-
laboratory work will not be accepted and will receive a score of zero. Your TA will evaluate the work and return it no later than the
beginning of your lab section. Your preparedness will also be evaluated by your answers to oral questions pertaining to the
laboratory from the TA and other lab instructors. The pre-lab and preparation accounts for 10 points of each experiment requiring a
pre-lab.
For every experimental module there is a general pre-lab assignment. It is described here:
General Pre-lab assignment:
Write a brief introduction to the experiment and the experimental plan using the following outline:
I. Experiment Introduction: In your own words, briefly describe the underlying theory and explain the principle upon which the
experiment isbased. Keep it concise; lengthy discussions or derivations of equations are unnecessary.
a. Background (importance/relevance): A brief description and justification of the importance of the topic, with references to
source material. What is the topic and why is it interesting and important? This is meant to "hook" the reader so they want
to read about the topic.
b. Theory: A concise paragraph or more describing the general theory for the method/technique used; this section should
contain a reference to the appropriate pages of the textbook, or other source material from the literature.
c. Goals/Purpose: Brief statement of purpose, which should indicate what is analyzed and the technique used. Limit to three to
five sentences.
II. Experimental Plan: a short summary of the specific things that you would need to know in order to do the experiment without
having access to your lab manual, with the understanding that you will always have a TA available to teach you how to use the
instrument. The best way to approach this task is to thoroughly read the procedure for the experiment. While reading, take notes
on the specifics of the work. What chemicals will you use? If you have to make solutions, what solutions do you need to make
and of what concentration? How will you prepare them? What instrument will you use, and what specific instrument settings
will you need to input (sometimes you will find this information in the Appendices that discuss the use of the instruments)?
How many runs of each sample will you do; and over what wavelength, or temperature, etc. range? You do not need to write
down how to use the instrument; we will teach you that. Bottom Line: we want you to write down only that pertinent
information about samples, sample preparation, instrument and instrument parameters, numbers of runs, and other things, that
you will need to know in order to do your work.
Pre-lab assignments should be typed. If you prefer to write out equations by hand, they should be neat and inserted as an image in
the appropriate place. If you need advice for how to insert equations or images into your documents, please ask for help at least 48
hours before the assignment is due; an in-person (or zoom) meeting is most appropriate for this issue.

Always Wear Protective Clothing:


Please dress appropriately for a laboratory environment - this is considered an important part of preparation for the laboratory
period (see also Laboratory Safety below). You must wear close-toes shoes and clothing must cover your body from your upper
arms to your knees. You should also wear protective eyewear. If lab attire is not jiving with your wardrobe, please purchase a lab
coat and bring it with you so that you can cover up in lab. No student should be in lab without wearing protective eyewear and
appropriate protective clothing.

3.1: Pre-lab Assignment is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.
Introductory Details by W. R. Fawcett, John Berg, P. B. Kelley, Carlito B. Lebrilla, Gang-yu Liu, Delmar Larsen, Paul Hrvatin, David
Goodin, and Brooke McMahon is licensed CC BY-NC-SA 4.0.

3.1.2 https://chem.libretexts.org/@go/page/453055
3.2: Learning Objectives, Important Info
 Learning Objectives

After completing the readings and practice problems recommended in this module, you should be able to:
Describe and give examples of the following types of error that effect physical measurements, and identify the types of
errors that affect precision and those that affect accuracy.
Gross / personal error
Systematic error
Random error
Determinate error
Indeterminate error
Describe strategies for optimizing the accuracy of physical measurements and evaluating the precision of physical
measurements.
Identify the sources of random error in a measurement.
Propagate error using mathematics.
Report values and their precision in the proper form (applying the rules of significant figures to correctly indicate
precision).

A summary of useful equations is below. A major goal of this module is for you to learn when and how to apply these equations to
propagate error. This skill will be required in nearly every future module in this course.

Statistics (on replicate measurements):


Average or mean: The average (or mean) of N measurements on the variable x is:
N
1
x̄ = ∑ xi (3.2.1)
N
i=1

Standard deviation (S ) is the spread around the mean of several replicate measurements. It has the same units as the measurement
and is commonly used to report uncertainty (see Skoog Ch. 6 p 126)
−−−−−−−−−−−−
N 2
∑ (xi − x̄)
i=1
S =√ (3.2.2)
N −1

Variance (S ) is the dispersion between data points. The units are the square of the measured units.
2

The Confidence Limits (CL) are the values that define the confidence interval (which is the bounds about which an experimental
value should be located with a given probability). The confidence limit is useful for small sample sizes. A common confidence
limit is the 95% confidence limit, and its value is given by the following (see Skoog Ch7 pg 150 and Shoemaker Ch 2 eq 34).
S
C L0.95 = t0.95Sm = t0.95 −
− (3.2.3)
√N

or more generally
tS
C L = x̄ ± −
− (3.2.4)
√N

where t is the critical t value (from Student's t test) that depends on the degrees of freedom and the given probability. A table of
critical t values is found online here (click). (See also Shoemaker Ch 2, Table 3)

3.2.1 https://chem.libretexts.org/@go/page/453056
Propagation of error:
Using partial derivatives
To propagate error to determine the error (Δ) in the value derived from of a function (F ), we must propagate error of each part of
the function. This is accomplished by taking the square root of the sum of the squares of the partial derivative of the function with
respect to each variable (eg partial derivative with respect to x is ), multiplied by the error in that variable (eg error in x is
∂F

∂x

Δ(x), and the square of the error is Δ (x). In mathematical form, this is given below (from Shoemaker Ch 2 p 56, eq 54). If you
2

need to brush up on partial derivatives, check out this video (click here). https://youtu.be/SbfRDBmyAMI
2 2 2
∂F ∂F ∂F
2 2 2 2
Δ (F ) = ( ) Δ (x) + ( ) Δ (y) + ( ) Δ (z) + … (3.2.5)
∂x ∂y ∂z

It is convenient to rearrange this equation as follows to give the error in the value derived from the function:
−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
2 2 2
∂F 2
∂F 2
∂F 2
Δ(F ) = √ ( ) Δ (x) + ( ) Δ (y) + ( ) Δ (z) + … (3.2.6)
∂x ∂y ∂z

Simpler shortcuts:
In some cases, there are shortcuts that allow you to approximate errors without taking the partial derivatives. In the examples
below, assume that a, b, c, n are constants and x, y, z are variables with associated errors Δx, Δy, Δz:
For functions with only addition and subtraction of values with assiciated error (eg F = ax ± by ± cz ), we can use the
following "shortcut": (Shoemaker, Ch 2, p 57, eq 55)
−−−−−−−−−−−−−−−−−−−−−−−
2 2 2 2 2 2
Δ(F ) = √ a Δ (x) + b Δ (y) + c Δ (z) (3.2.7)

For functions that involve only multiplication and division (eg F = axyz or F = axy/z or F = a/xyz ), we can use the
following "shortcut": (Shoemaker, Ch 2, p 57, eq 56)
−−−−−−−−−−−−−−−−−−−−
2 2 2
Δ (x) Δ (y) Δ (z)
Δ(F ) = F √ + + (3.2.8)
x2 y2 z2

For functions that involve taking a variable to the power of a constant (eg F = ax
n
), we can use the following "shortcut":
(Shoemaker, Ch 2, p 57, eq 57)
−−−−−−−−
2
Δ (x) Δ(x)
2
Δ(F ) = F √n → Δ(F ) = F n (3.2.9)
2
x x

For functions where a variable is an exponent (eg F = ae


x
), we can use the following "shortcut": (Shoemaker, Ch 2, p 57, eq 58)
−−−−−−−−−
2 2x 2
Δ(F ) = √a e Δ (x) → Δ(F ) = F Δ(x) (3.2.10)

For functions taking the logarithm of a variable (eg F = a ln(x) ), we can use the following "shortcut": (Shoemaker, Ch 2, p 57,
eq 59)
−−−−−−−−
2
a Δ(x)
2
Δ(F ) = √ Δ (x) → Δ(F ) = a (3.2.11)
2
x x

3.2: Learning Objectives, Important Info is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

3.2.2 https://chem.libretexts.org/@go/page/453056
3.3: Characterizing Experimental Errors
In this course, you will need to determine the confidence in your results by doing appropriate error analysis. We assume that you
are familiar with error analysis from the prerequisite courses (P-chem labs CHEM301L or CHEM310L). This and the next section
are just here as a reminder about what you already know. You might be glad to hear that you will rarely need to use differential
equations to perform error analysis in this course. This section and the next are meant to review the types of experimental errors
that affect accuracy and precision and to remind you of the simple methods for propagating error.

Errors That Affect Accuracy


Accuracy is how close a measure of central tendency is to its expected value, μ . We express accuracy either as an absolute error, e
¯¯¯
¯
e = X −μ (3.3.1)

or as a percent relative error, %e


¯¯¯
¯
X −μ
%e = × 100 (3.3.2)
μ

Although Equation 3.3.1 and Equation 3.3.2 use the mean as the measure of central tendency, we also can use the median.

The convention for representing a statistical parameter is to use a Roman letter for a value calculated from experimental data,
¯¯¯
¯
and a Greek letter for its corresponding expected value. For example, the experimentally determined mean is X and its
underlying expected value is μ . Likewise, the experimental standard deviation is s and the underlying expected value is σ.

We identify as determinate an error that affects the accuracy of an analysis. Each source of a determinate error has a specific
magnitude and sign. Some sources of determinate error are positive and others are negative, and some are larger in magnitude and
others are smaller in magnitude. The cumulative effect of these determinate errors is a net positive or negative error in accuracy.

It is possible, although unlikely, that the positive and negative determinate errors will offset each other, producing a result with
no net error in accuracy.

We assign determinate errors into four categories—sampling errors, method errors, measurement errors, and personal errors—each
of which we consider in this section.

Sampling Errors
A determinate sampling error occurs when our sampling strategy does not provide a us with a representative sample. For example,
if we monitor the environmental quality of a lake by sampling from a single site near a point source of pollution, such as an outlet
for industrial effluent, then our results will be misleading. To determine the mass of a U. S. penny, our strategy for selecting
pennies must ensure that we do not include pennies from other countries.

An awareness of potential sampling errors especially is important when we work with heterogeneous materials. Strategies for
obtaining representative samples are covered in Chapter 5 of Harvey's Analytical Chemistry text on LibreTexts.

Method Errors
In any analysis the relationship between the signal, Stotal, and the absolute amount of analyte, nA, or the analyte’s concentration, CA,
is
Stotal = kA nA + Smb (3.3.3)

Stotal = kA CA + Smb (3.3.4)

where kA is the method’s sensitivity for the analyte and Smb is the signal from the method blank. A method error exists when our
value for kA or for Smb is in error. For example, a method in which Stotal is the mass of a precipitate assumes that k is defined by a

3.3.1 https://chem.libretexts.org/@go/page/406453
pure precipitate of known stoichiometry. If this assumption is not true, then the resulting determination of nA or CA is inaccurate.
We can minimize a determinate error in kA by calibrating the method. A method error due to an interferent in the reagents is
minimized by using a proper method blank.

Measurement Errors
The manufacturers of analytical instruments and equipment, such as glassware and balances, usually provide a statement of the
item’s maximum measurement error, or tolerance. For example, a 10-mL volumetric pipet (Figure 3.3.1 ) has a tolerance of ±0.02
mL, which means the pipet delivers an actual volume within the range 9.98–10.02 mL at a temperature of 20 oC. Although we
express this tolerance as a range, the error is determinate; that is, the pipet’s expected volume, μ , is a fixed value within this stated
range.

Figure 3.3.1 : Close-up of a 10-mL volumetric pipet showing that it has a tolerance of ±0.02 mL at 20 oC.
Volumetric glassware is categorized into classes based on its relative accuracy. Class A glassware is manufactured to comply with
tolerances specified by an agency, such as the National Institute of Standards and Technology or the American Society for Testing
and Materials. The tolerance level for Class A glassware is small enough that normally we can use it without calibration. The
tolerance levels for Class B glassware usually are twice that for Class A glassware. In other words, Class B is less accurate. Other
types of volumetric glassware, such as beakers and graduated cylinders, are not used to measure volume accurately. Table 3.3.1
provides a summary of typical measurement errors for Class A volumetric glassware. Tolerances for digital pipets and for balances
are provided in Table 3.3.2 and Table 3.3.3 .
Tolerance Limits for Class A Volumetric Glassware
Table 3.3.1 : Measurement Errors for Class A Volumetric Glassware
Volumetric Pipets Volumetric Flasks Burets

Capacity (mL) Tolerance (mL) Capacity (mL) Tolerance (mL) Capacity (mL) Tolerance (mL)

1 ±0.006 5 ±0.02 10 ±0.02

2 ±0.006 10 ±0.02 25 ±0.03

5 ±0.01 25 ±0.03 50 ±0.05

10 ±0.02 50 ±0.05

20 ±0.03 100 ±0.08

25 ±0.03 250 ±0.12

50 ±0.05 500 ±0.20

100 ±0.08 1000 ±0.30

2000 ±0.50

Tolerance limits for Digital Pipettes


Table 3.3.2 : Measurement Errors for Digital Pipets
Pipet Range Volume (mL or μL) Percent Measurement Error

10–100 μL 10 ±3.0%

50 ±1.0%

100 ±0.8%

100–1000 μL 100 ±3.0%

3.3.2 https://chem.libretexts.org/@go/page/406453
Pipet Range Volume (mL or μL) Percent Measurement Error

500 ±1.0%

1000 ±0.6%

1–10 mL 1 ±3.0%

5 ±0.8%

10 ±0.6%

The tolerance values for the volumetric glassware in Table 3.3.1 are from the ASTM E288, E542, and E694 standards. The
measurement errors for the digital pipets in Table 3.3.2 are from www.eppendorf.com.

Tolerance Limits for Some Common Balances


Table 3.3.3 : Measurement Errors for Selected Balances
Balance Capacity (g) Measurement Error

Precisa 160M 160 ±1 mg

A & D ER 120M 120 ±0.1 mg

Metler H54 160 ±0.01 mg

We can minimize a determinate measurement error by calibrating our equipment. Balances are calibrated using a reference weight
whose mass we can trace back to the SI standard kilogram. Volumetric glassware and digital pipets are calibrated by determining
the mass of water delivered or contained and using the density of water to calculate the actual volume. It is never safe to assume
that a calibration does not change during an analysis or over time. One study, for example, found that repeatedly exposing
volumetric glassware to higher temperatures during machine washing and oven drying, led to small, but significant changes in the
glassware’s calibration [Castanheira, I.; Batista, E.; Valente, A.; Dias, G.; Mora, M.; Pinto, L.; Costa, H. S. Food Control 2006, 17,
719–726]. Many instruments drift out of calibration over time and may require frequent recalibration during an analysis.

Personal Errors
Finally, analytical work is always subject to personal error, examples of which include the ability to see a change in the color of an
indicator that signals the endpoint of a titration, biases, such as consistently overestimating or underestimating the value on an
instrument’s readout scale, failing to calibrate instrumentation, and misinterpreting procedural directions. You can minimize
personal errors by taking proper care.

Identifying Determinate Errors


Determinate errors often are difficult to detect. Without knowing the expected value for an analysis, the usual situation in any
analysis that matters, we often have nothing to which we can compare our experimental result. Nevertheless, there are strategies we
can use to detect determinate errors.
The magnitude of a constant determinate error is the same for all samples and is more significant when we analyze smaller
samples. Analyzing samples of different sizes, therefore, allows us to detect a constant determinate error. For example, consider a
quantitative analysis in which we separate the analyte from its matrix and determine its mass. Let’s assume the sample is 50.0%
w/w analyte. As we see in Table 3.3.4 , the expected amount of analyte in a 0.100 g sample is 0.050 g. If the analysis has a positive
constant determinate error of 0.010 g, then analyzing the sample gives 0.060 g of analyte, or an apparent concentration of 60.0%
w/w. As we increase the size of the sample the experimental results become closer to the expected result. An upward or downward
trend in a graph of the analyte’s experimental concentration versus the sample’s mass (Figure 3.3.2 ) is evidence of a constant
determinate error.
Table 3.3.4 : Effect of a Constant Determinate Error on the Analysis of a Sample That is 50.0% w/w Analyte
Experimental
Expected Mass Experimental
Mass of Sample (g) Constant Error (g) Concentration of Analyte
of Analyte (g) Mass of Analyte (g)
(% w/w)

0.100 0.050 0.010 0.060 60.0

3.3.3 https://chem.libretexts.org/@go/page/406453
Experimental
Expected Mass Experimental
Mass of Sample (g) Constant Error (g) Concentration of Analyte
of Analyte (g) Mass of Analyte (g)
(% w/w)

0.200 0.100 0.010 0.110 55.0

0.400 0.200 0.010 0.210 52.5

0.800 0.400 0.010 0.410 51.2

1.600 0.800 0.010 0.810 50.6

Figure 3.3.2 : Effect of a constant positive determinate error of +0.01 g and a constant negative determinate error of –0.01 g on the
determination of an analyte in samples of varying size. The analyte’s expected concentration of 50% w/w is shown by the dashed
line.
A proportional determinate error, in which the error’s magnitude depends on the amount of sample, is more difficult to detect
because the result of the analysis is independent of the amount of sample. Table 3.3.5 outlines an example that shows the effect of a
positive proportional error of 1.0% on the analysis of a sample that is 50.0% w/w in analyte. Regardless of the sample’s size, each
analysis gives the same result of 50.5% w/w analyte.
Table 3.3.5 : Effect of a Proportional Determinate Error on the Analysis of a Sample That is 50.0% w/w Analyte
Experimental
Expected Mass Proportional Experimental
Mass of Sample (g) Concentration of Analyte
of Analyte (g) Error (%) Mass of Analyte (g)
(% w/w)

0.100 0.050 1.00 0.0505 50.5

0.200 0.100 1.00 0.101 50.5

0.400 0.200 1.00 0.202 50.5

0.800 0.400 1.00 0.404 50.5

1.600 0.800 1.00 0.808 50.5

One approach for detecting a proportional determinate error is to analyze a standard that contains a known amount of analyte in a
matrix similar to our samples. Standards are available from a variety of sources, such as the National Institute of Standards and
Technology (where they are called Standard Reference Materials) or the American Society for Testing and Materials. Table 3.3.6 ,
for example, lists certified values for several analytes in a standard sample of Gingko biloba leaves. Another approach is to
compare our analysis to an analysis carried out using an independent analytical method that is known to give accurate results. If the
two methods give significantly different results, then a determinate error is the likely cause.
Table 3.3.6 : Certified Concentrations for SRM 3246: Gingko bilbo (Leaves)
Class of Analyte Analyte Mass Fraction (mg/g or ng/g)

Flavonoids/Ginkgolide B (mass fraction in


Qurecetin 2.69 ± 0.31
mg/g)

Kaempferol 3.02 ± 0.41

Isorhamnetin 0.517 ± 0.0.99

3.3.4 https://chem.libretexts.org/@go/page/406453
Class of Analyte Analyte Mass Fraction (mg/g or ng/g)

Total Aglycones 6.22 ± 0.77

Selected Terpenes (mass fraction in mg/g) Ginkgolide A 0.57 ± 0.28

Ginkgolide B 0.470 ± 0.090

Ginkgolide C 0.59 ± 0.22

Ginkgolide J 0.18 ± 0.10

Bilobalide 1.52 ± 0.40

Total Terpene Lactones 3.3 ± 1.1

Selected Toxic Elements (mass fraction in


Cadmium 20.8 ± 1.0
ng/g)

Lead 995 ± 30

Mercury 23.08 ± 0.17

The primary purpose of this Standard Reference Material is to validate analytical methods for determining flavonoids, terpene
lactones, and toxic elements in Ginkgo biloba or other materials with a similar matrix. Values are from the official Certificate
of Analysis available at www.nist.gov.

Constant and proportional determinate errors have distinctly different sources, which we can define in terms of the relationship
between the signal and the moles or concentration of analyte (Equation 3.3.3 and Equation 3.3.4). An invalid method blank, Smb, is
a constant determinate error as it adds or subtracts the same value to the signal. A poorly calibrated method, which yields an invalid
sensitivity for the analyte, kA, results in a proportional determinate error.

Errors that Affect Precision


Precision is a measure of the spread of individual measurements or results about a central value, which we express as a range, a
standard deviation, or a variance. Here we draw a distinction between two types of precision: repeatability and reproducibility.
Repeatability is the precision when a single analyst completes an analysis in a single session using the same solutions, equipment,
and instrumentation. Reproducibility, on the other hand, is the precision under any other set of conditions, including between
analysts or between laboratory sessions for a single analyst. Since reproducibility includes additional sources of variability, the
reproducibility of an analysis cannot be better than its repeatability.

The ratio of the standard deviation associated with reproducibility to the standard deviation associated with repeatability is
called the Horowitz ratio. For a wide variety of analytes in foods, for example, the median Horowtiz ratio is 2.0 with larger
values for fatty acids and for trace elements; see Thompson, M.; Wood, R. “The ‘Horowitz Ratio’–A Study of the Ratio
Between Reproducibility and Repeatability in the Analysis of Foodstuffs,” Anal. Methods, 2015, 7, 375–379.

Errors that affect precision are indeterminate and are characterized by random variations in their magnitude and their direction.
Because they are random, positive and negative indeterminate errors tend to cancel, provided that we make a sufficient number of
measurements. In such situations the mean and the median largely are unaffected by the precision of the analysis.

Sources of Indeterminate Error


We can assign indeterminate errors to several sources, including collecting samples, manipulating samples during the analysis, and
making measurements. When we collect a sample, for instance, only a small portion of the available material is taken, which
increases the chance that small-scale inhomogeneities in the sample will affect repeatability. Individual pennies, for example, may
show variations in mass from several sources, including the manufacturing process and the loss of small amounts of metal or the
addition of dirt during circulation. These variations are sources of indeterminate sampling errors.
During an analysis there are many opportunities to introduce indeterminate method errors. If our method for determining the mass
of a penny includes directions for cleaning them of dirt, then we must be careful to treat each penny in the same way. Cleaning

3.3.5 https://chem.libretexts.org/@go/page/406453
some pennies more vigorously than others might introduce an indeterminate method error.
Finally, all measuring devices are subject to indeterminate measurement errors due to limitations in our ability to read its scale. For
example, a buret with scale divisions every 0.1 mL has an inherent indeterminate error of ±0.01–0.03 mL when we estimate the
volume to the hundredth of a milliliter (Figure 3.3.3 ).

Figure 3.3.3 : Close-up of a buret showing the difficulty in estimating volume. With scale divisions every 0.1 mL it is difficult to
read the actual volume to better than ±0.01–0.03 mL.

Evaluating Indeterminate Error


Indeterminate errors associated with our analytical equipment or instrumentation generally are easy to estimate if we measure the
standard deviation for several replicate measurements, or if we monitor the signal’s fluctuations over time in the absence of analyte
(Figure 3.3.4 ) and calculate the standard deviation. Other sources of indeterminate error, such as treating samples inconsistently,
are more difficult to estimate.

Figure 3.3.4 : Background noise in an instrument showing the random fluctuations in the signal.

Error and Uncertainty


Analytical chemists make a distinction between error and uncertainty [Ellison, S.; Wegscheider, W.; Williams, A. Anal. Chem.
1997, 69, 607A–613A]. Error is the difference between a single measurement or result and its expected value. In other words, error
is a measure of bias. As discussed earlier, we divide errors into determinate and indeterminate sources. Although we can find and
correct a source of determinate error, the indeterminate portion of the error remains.
Uncertainty expresses the range of possible values for a measurement or result. Note that this definition of uncertainty is not the
same as our definition of precision. We calculate precision from our experimental data and use it to estimate the magnitude of
indeterminate errors. Uncertainty accounts for all errors—both determinate and indeterminate—that reasonably might affect a
measurement or a result. Although we always try to correct determinate errors before we begin an analysis, the correction itself is
subject to uncertainty.
Here is an example to help illustrate the difference between precision and uncertainty. Suppose you purchase a 10-mL Class A
pipet from a laboratory supply company and use it without any additional calibration. The pipet’s tolerance of ±0.02 mL is its
uncertainty because your best estimate of its expected volume is 10.00 mL ± 0.02 mL. This uncertainty primarily is determinate. If
you use the pipet to dispense several replicate samples of a solution and determine the volume of each sample, the resulting
standard deviation is the pipet’s precision. Table 3.3.8 shows results for ten such trials, with a mean of 9.992 mL and a standard
deviation of ±0.006 mL. This standard deviation is the precision with which we expect to deliver a solution using a Class A 10-mL
pipet. In this case the pipet’s published uncertainty of ±0.02 mL is worse than its experimentally determined precision of ±0.006
ml. Interestingly, the data in Table 3.3.8 allows us to calibrate this specific pipet’s delivery volume as 9.992 mL. If we use this

3.3.6 https://chem.libretexts.org/@go/page/406453
volume as a better estimate of the pipet’s expected volume, then its uncertainty is ±0.006 mL. As expected, calibrating the pipet
allows us to decrease its uncertainty [Kadis, R. Talanta 2004, 64, 167–173].
Table 3.3.8 : Experimental Results for Volume Dispensed by a 10–mL Class A Transfer Pipet
Replicate Volume (ml) Replicate Volume (mL)

1 10.002 6 9.983

2 9.993 7 9.991

3 9.984 8 9.990

4 9.996 9 9.988

5 9.989 10 9.999

This page titled 3.3: Characterizing Experimental Errors is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated
by Kathryn Haas.
4.2: Characterizing Experimental Errors by David Harvey is licensed CC BY-NC-SA 4.0.

3.3.7 https://chem.libretexts.org/@go/page/406453
3.4: Propagation of Uncertainty
Suppose we dispense 20 mL of a reagent using the Class A 10-mL pipet whose calibration information is given in the Table in the
previous section. If the volume and uncertainty for one use of the pipet is 9.992 ± 0.006 mL, what is the volume and uncertainty if
we use the pipet twice?
As a first guess, we might simply add together the volume and the maximum uncertainty for each delivery; thus
(9.992 mL + 9.992 mL) ± (0.006 mL + 0.006 mL) = 19.984 ± 0.012 mL
It is easy to appreciate that combining uncertainties in this way overestimates the total uncertainty. Adding the uncertainty for the
first delivery to that of the second delivery assumes that with each use the indeterminate error is in the same direction and is as
large as possible. At the other extreme, we might assume that the uncertainty for one delivery is positive and the other is negative.
If we subtract the maximum uncertainties for each delivery,
(9.992 mL + 9.992 mL) ± (0.006 mL – 0.006 mL) = 19.984 ± 0.000 mL
we clearly underestimate the total uncertainty.
So what is the total uncertainty? From the discussion above, we reasonably expect that the total uncertainty is greater than ±0.000
mL and that it is less than ±0.012 mL. To estimate the uncertainty we use a mathematical technique known as the propagation of
uncertainty. Our treatment of the propagation of uncertainty is based on a few simple rules.

A Few Symbols
A propagation of uncertainty allows us to estimate the uncertainty in a result from the uncertainties in the measurements used to
calculate that result. For the equations in this section we represent the result with the symbol R, and we represent the measurements
with the symbols A, B, and C. The corresponding uncertainties are uR, uA, uB, and uC. We can define the uncertainties for A, B, and
C using standard deviations, ranges, or tolerances (or any other measure of uncertainty), as long as we use the same form for all
measurements.

The requirement that we express each uncertainty in the same way is a critically important point. Suppose you have a range for
one measurement, such as a pipet’s tolerance, and standard deviations for the other measurements. All is not lost. There are
ways to convert a range to an estimate of the standard deviation. See Appendix 2 of Harvey's Analytical Chemistry Text for
more details.

Uncertainty When Adding or Subtracting


When we add or subtract measurements we propagate their absolute uncertainties. For example, if the result is given by the
equation

R = A+B−C

the the absolute uncertainty in R is


−−−−−−−−−−−
2 2 2
uR = √ u +u +u (3.4.1)
A B C

 Example 3.4.1

If we dispense 20 mL using a 10-mL Class A pipet, what is the total volume dispensed and what is the uncertainty in this
volume? First, complete the calculation using the manufacturer’s tolerance of 10.00 mL±0.02 mL, and then using the
calibration data from the Table 1.3.8 in the previous section.
Solution
To calculate the total volume we add the volumes for each use of the pipet. When using the manufacturer’s values, the total
volume is

V = 10.00 mL + 10.00 mL = 20.00 mL

3.4.1 https://chem.libretexts.org/@go/page/406454
and when using the calibration data, the total volume is

V = 9.992 mL + 9.992 mL = 19.984 mL

Using the pipet’s tolerance as an estimate of its uncertainty gives the uncertainty in the total volume as
2 2
uR = (0.02 ) + (0.02 ) = 0.028 mL = 0.028 mL

and using the standard deviation for the data in Table 1.3.8 gives an uncertainty of
2 2
uR = (0.006 ) + (0.006 ) = 0.0085 mL

Rounding the volumes to four significant figures gives 20.00 mL ± 0.03 mL when we use the tolerance values, and 19.98 ±
0.01 mL when we use the calibration data.

Uncertainty When Multiplying or Dividing


When we multiple or divide measurements we propagate their relative uncertainties. For example, if the result is given by the
equation
A×B
R =
C

then the relative uncertainty in R is


−−−−−−−−−−−−−−−−−−−−− −
2 2 2
uR uA uB uC
= √( ) +( ) +( ) (3.4.2)
R A B C

 Example 3.4.2

The quantity of charge, Q, in coulombs that passes through an electrical circuit is

Q = i ×t

where i is the current in amperes and t is the time in seconds. When a current of 0.15 A ± 0.01 A passes through the circuit for
120 s ± 1 s, what is the total charge and its uncertainty?
Solution
The total charge is

Q = (0.15 A) × (120 s) = 18 C

Since charge is the product of current and time, the relative uncertainty in the charge is
−−−−−−−−−−−−−−−−−
2 2
0.01 1
uR = √( ) +( ) = 0.0672
0.15 120

and the charge’s absolute uncertainty is

uR = R × 0.0672 = (18 C) × (0.0672) = 1.2 C

Thus, we report the total charge as 18 C ± 1 C.

Uncertainty for Mixed Operations


Many chemical calculations involve a combination of adding and subtracting, and of multiply and dividing. As shown in the
following example, we can calculate the uncertainty by separately treating each operation using Equation 3.4.1 and Equation 3.4.2
as needed.

 Example 3.4.3
For a concentration technique, the relationship between the signal and the an analyte’s concentration is

3.4.2 https://chem.libretexts.org/@go/page/406454
Stotal = kA CA + Smb

What is the analyte’s concentration, CA, and its uncertainty if Stotal is 24.37 ± 0.02, Smb is 0.96 ± 0.02, and kA is
0.186 ± 0.003 ppm ?
−1

Solution
Rearranging the equation and solving for CA
Stotal − Smb 24.37 − 0.96 23.41
CA = = = = 125.9 ppm
−1 −1
kA 0.186 ppm 0.186 ppm

gives the analyte’s concentration as 126 ppm. To estimate the uncertainty in CA, we first use Equation 3.4.1 to determine the
uncertainty for the numerator.
−−−−−−−−−−−−−
2 2
uR = √ (0.02 ) + (0.02 ) = 0.028

The numerator, therefore, is 23.41 ± 0.028. To complete the calculation we use Equation 3.4.2 to estimate the relative
uncertainty in CA.
−−−−−−−−−−−−−−−−−−−
2 2
uR 0.028 0.003
= √( ) +( ) = 0.0162
R 23.41 0.186

The absolute uncertainty in the analyte’s concentration is

uR = (125.9 ppm) × (0.0162) = 2.0 ppm

Thus, we report the analyte’s concentration as 126 ppm ± 2 ppm.

 Exercise 3.4.1

To prepare a standard solution of Cu2+ you obtain a piece of copper from a spool of wire. The spool’s initial weight is 74.2991
g and its final weight is 73.3216 g. You place the sample of wire in a 500-mL volumetric flask, dissolve it in 10 mL of HNO3,
and dilute to volume. Next, you pipet a 1 mL portion to a 250-mL volumetric flask and dilute to volume. What is the final
concentration of Cu2+ in mg/L, and its uncertainty? Assume that the uncertainty in the balance is ±0.1 mg and that you are
using Class A glassware.

Answer
The first step is to determine the concentration of Cu2+ in the final solution. The mass of copper is

74.2991 g − 73.3216 g = 0.9775 g Cu

The 10 mL of HNO3 used to dissolve the copper does not factor into our calculation. The concentration of Cu2+ is
0.9775 g Cu 1.000 mL 1000 mg
2 +
× × = 7.820 mg Cu /L
0.5000 L 250.0 mL g

Having found the concentration of Cu2+, we continue with the propagation of uncertainty. The absolute uncertainty in the
mass of Cu wire is
−−−−−−−−−−−−−−−−−
2 2
ug Cu = √ (0.0001 ) + (0.0001 ) = 0.00014 g

The relative uncertainty in the concentration of Cu2+ is


−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
2 2 2 2
umg/L 0.00014 0.20 0.006 0.12
=√ ( ) +( ) +( ) +( ) = 0.00603
7.820 mg/L 0.9775 500.0 1.000 250.0

Solving for umg/L gives the uncertainty as 0.0472. The concentration and uncertainty for Cu2+ is 7.820 mg/L ± 0.047 mg/L.

3.4.3 https://chem.libretexts.org/@go/page/406454
Uncertainty for Other Mathematical Functions
Many other mathematical operations are common in analytical chemistry, including the use of powers, roots, and logarithms. Table
3.4.1 provides equations for propagating uncertainty for some of these function where A and B are independent measurements and
where k is a constant whose value has no uncertainty.
Table 3.4.1 : Propagation of Uncertainty for Selected Mathematical Functions
Function uR Function uR

uA
R = kA uR = kuA R = ln(A) uR =
A

−−−−−−−
uA
2 2
R = A +B uR = √u +u R = log(A) uR = 0.4343 ×
A B A

−−−−−−−
2 2 A uR
R = A −B uR = √u +u R = e = uA
A B R

−−−−−−−−−−−−
2 uR
uA uB 2 A
R = A ×B uR = √( ) +( ) R = 10 = 2.303 × uA
A B R

−−−−−−−−−−−−
2 uR uA
A uA uB 2 k
R = uR = √( ) +( ) R = A = k×
B A B R A

 Example 3.4.4

If the pH of a solution is 3.72 with an absolute uncertainty of ±0.03, what is the [H+] and its uncertainty?
Solution
The concentration of H+ is
+ −pH −3.72 −4
[H ] = 10 = 10 = 1.91 × 10  M

or 1.9 × 10−4
M to two significant figures. From Table 3.4.1 the relative uncertainty in [H+] is
uR
= 2.303 × uA = 2.303 × 0.03 = 0.069
R

The uncertainty in the concentration, therefore, is


−4 −5
(1.91 × 10  M) × (0.069) = 1.3 × 10  M

We report the [H+] as 1.9(±0.1) × 10 −4


M, which is equivalent to 1.9 × 10 −4
 M  ± 0.1 × 10
−4
 M .

 Exercise 3.4.2

A solution of copper ions is blue because it absorbs yellow and orange light. Absorbance, A, is defined as
P
A = − log T = − log( )
Po

where, T is the transmittance, Po is the power of radiation as emitted from the light source and P is its power after it passes
through the solution. What is the absorbance if Po is 3.80 × 10 and P is 1.50 × 10 ? If the uncertainty in measuring Po and P
2 2

is 15, what is the uncertainty in the absorbance?

Answer
The first step is to calculate the absorbance, which is
2
P 1.50 × 10
A = − log T = − log = − log = 0.4037 ≈ 0.404
Po 2
3.80 × 10

Having found the absorbance, we continue with the propagation of uncertainty. First, we find the uncertainty for the ratio
P/Po, which is the transmittance, T.

3.4.4 https://chem.libretexts.org/@go/page/406454
−−−−−−−−−−−−−−−−−−−−−−−−−−
2 2
uT 15 15
= √( ) +( ) = 0.1075
2 2
T 3.80 × 10 1.50 × 10

Finally, from Table 3.4.1 the uncertainty in the absorbance is


uT −2
uA = 0.4343 × = (0.4343) × (0.1075) = 4.669 × 10
T

The absorbance and uncertainty is 0.40 ± 0.05 absorbance units.

Is Calculating Uncertainty Actually Useful?


Given the effort it takes to calculate uncertainty, it is worth asking whether such calculations are useful. The short answer is, yes.
Let’s consider three examples of how we can use a propagation of uncertainty to help guide the development of an analytical
method.
One reason to complete a propagation of uncertainty is that we can compare our estimate of the uncertainty to that obtained
experimentally. For example, to determine the mass of a penny we measure its mass twice—once to tare the balance at 0.000 g and
once to measure the penny’s mass. If the uncertainty in each measurement of mass is ±0.001 g, then we estimate the total
uncertainty in the penny’s mass as
−−−−−−−−−−−−−−−
2 2
uR = √ (0.001 ) + (0.001 ) = 0.0014 g

If we measure a single penny’s mass several times and obtain a standard deviation of ±0.050 g, then we have evidence that the
measurement process is out of control. Knowing this, we can identify and correct the problem.
We also can use a propagation of uncertainty to help us decide how to improve an analytical method’s uncertainty. In Example
3.4.3 , for instance, we calculated an analyte’s concentration as 126 ppm ± 2 ppm, which is a percent uncertainty of 1.6%. Suppose
we want to decrease the percent uncertainty to no more than 0.8%. How might we accomplish this? Looking back at the
calculation, we see that the concentration’s relative uncertainty is determined by the relative uncertainty in the measured signal
(corrected for the reagent blank)
0.028
= 0.0012 or 0.12%
23.41

and the relative uncertainty in the method’s sensitivity, kA,


−1
0.003 ppm
= 0.016 or 1.6%
−1
0.186 ppm

Of these two terms, the uncertainty in the method’s sensitivity dominates the overall uncertainty. Improving the signal’s uncertainty
will not improve the overall uncertainty of the analysis. To achieve an overall uncertainty of 0.8% we must improve the uncertainty
in kA to ±0.0015 ppm–1.

 Exercise 3.4.3

Verify that an uncertainty of ±0.0015 ppm–1 for kA is the correct result.

Answer
An uncertainty of 0.8% is a relative uncertainty in the concentration of 0.008; thus, letting u be the uncertainty in kA
−−−−−−−−−−−−−−−−−−
2
2
0.028 u
0.008 = √ ( ) +( )
23.41 0.186

Squaring both sides of the equation gives


2
0.028 u 2
−5
6.4 × 10 =( ) +( )
23.41 0.186

3.4.5 https://chem.libretexts.org/@go/page/406454
Solving for the uncertainty in kA gives its value as 1.47 × 10 −3
or ±0.0015 ppm–1.

Finally, we can use a propagation of uncertainty to determine which of several procedures provides the smallest uncertainty. When
we dilute a stock solution usually there are several combinations of volumetric glassware that will give the same final
concentration. For instance, we can dilute a stock solution by a factor of 10 using a 10-mL pipet and a 100-mL volumetric flask, or
using a 25-mL pipet and a 250-mL volumetric flask. We also can accomplish the same dilution in two steps using a 50-mL pipet
and 100-mL volumetric flask for the first dilution, and a 10-mL pipet and a 50-mL volumetric flask for the second dilution. The
overall uncertainty in the final concentration—and, therefore, the best option for the dilution—depends on the uncertainty of the
volumetric pipets and volumetric flasks. As shown in the following example, we can use the tolerance values for volumetric
glassware to determine the optimum dilution strategy [Lam, R. B.; Isenhour, T. L. Anal. Chem. 1980, 52, 1158–1161].

 Example 3.4.5 :

Which of the following methods for preparing a 0.0010 M solution from a 1.0 M stock solution provides the smallest overall
uncertainty?

(a) A one-step dilution that uses a 1-mL pipet and a 1000-mL volumetric flask.

(b) A two-step dilution that uses a 20-mL pipet and a 1000-mL volumetric flask for the first dilution, and a 25-mL pipet and a
500-mL volumetric flask for the second dilution.
Solution
The dilution calculations for case (a) and case (b) are
1.000 mL
case (a): 1.0 M  × = 0.0010 M
1000.0 mL

20.00 mL 25.00 mL
case (b): 1.0 M  × × = 0.0010 M
1000.0 mL 500.0mL

Using tolerance values from Table 2.1.3.1, the relative uncertainty for case (a) is
−−−−−−−−−−−−−−−−−−−−
2 2
0.006 0.3
uR = √ ( ) +( ) = 0.006
1.000 1000.0

and for case (b) the relative uncertainty is


−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
2 2 2 2
0.03 0.3 0.03 0.2
uR = √ ( ) +( ) +( ) +( ) = 0.002
20.00 1000 25.00 500.0

Since the relative uncertainty for case (b) is less than that for case (a), the two-step dilution provides the smallest overall
uncertainty. Of course we must balance the smaller uncertainty for case (b) against the increased opportunity for introducing a
determinate error when making two dilutions instead of just one dilution, as in case (a).

This page titled 3.4: Propagation of Uncertainty is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by
Kathryn Haas.
4.3: Propagation of Uncertainty by David Harvey is licensed CC BY-NC-SA 4.0.

3.4.6 https://chem.libretexts.org/@go/page/406454
3.5: Least Squares Linear Regression
Your experimental data will usually be examined in graphic form, which can be more instructive than tabulated data. In any graphic
presentation, be sure to plot your independent variable on the abscissa (x-axis) and the dependent variable on the ordinate (y-axis).
If the data are presented in graphic form, the usual procedure is to fit the best curve through the experimental points. As Brownlee3
puts it,
"When an investigator observes simultaneously two variables, x and y, usually with good reason he plots the observations on
a graph. If there is any sign of an association, he is usually seized with the impulse to fit a line, usually a straight line or,
rather infrequently, a parabola or cubic."
When we are "seized with the impulse" to fit a straight line curve the process of curve fitting is called linear regression. We assume
that the deviations of the dependent variable y are distributed normally (that is, a Gaussian distribution). This assumption is not
required for the determination of the least squares parameters, but it is necessary for construction of confidence intervals or for tests
of hypotheses about the parameters. The parameters obtained by least squares under these conditions are unbiased and provide the
best estimates of the curve fitting parameters.
Most of the experimental data that you will encounter this semester will have a simple functional form, or can be manipulated in
such a way as to assume a simple form. We shall assume here that our data follows a linear relationship. For example, consider a
gas that is thought to follow the equation of state

PV
Z ≡ = B0 + B1 P (3.5.1)
nRT

y = a + bx

Measurements of Z, called the compressibility factor, are made at each of a series of pressures in an effort to determine B andB . 0 0

The least squares approach can be used to determine the slope (b, orB ) and y-intercept (a, or B ) which define the best straight
1 0

line that can be drawn through the data. In simple linear regression analysis, it is assumed that all the error in the measurement lies
in the dependent variable (y). This is equivalent to assuming that the precision of the determination of the independent variable (x)
is considerable higher than that of y. The measured value of some quantity (yi) will differ from the value predicted by a linear
equation (y = a + bx ) by an amount ri called the residual.

ri = yi − a − bx (3.5.2)

The least squares approach identifies the best straight line as the one for which R, the sum of the squares of the residuals, has the
smallest value:
N N

2 2
R = ∑r = ∑ (yi − a − b xi ) (3.5.3)
i

i=1 i=1

The arrows in Figure 4 represent the residuals, and all of them are equally weighted.

Figure 3.5.1 : Least Squares Linear Regression: Illustration showing the evaluation of a linear regression in which we assume that
all uncertainty is the result of indeterminate errors affecting y. The points in blue, yi, are the original data and the points in red, ŷi,
are the predicted values from the regression equation, ŷ = b0 + b1x. The smaller the total residual error, the better the fit of the
straight-line to the data (source)

3.5.1 https://chem.libretexts.org/@go/page/371024
We want to minimize the sum of the squares of the residuals. That is, we want to make least the sum of squares of residuals, and,
thus, the name of the procedure. Since this is a minimization problem, we require that R be as small as possible with respect to a
and b by taking derivatives with respect to these two parameters and setting the derivatives equal to zero.

∂R
( ) = −2 ∑(y − a − b xi ) = 0 (3.5.4)
i
∂a b

∂R
( ) = −2 ∑((yi − a − b xi )(xi )) = 0
∂b
a

Solving equations 3.5.4 simultaneously, we obtain


2
(∑ x )(∑ yi ) − (∑ xi )(∑ xi yi )
i
a = (3.5.5)
2 2
N (∑ x ) − (∑ xi )
i

N (∑ xi yi ) − (∑ xi )(∑ yi )
b = (3.5.6)
2 2
N (∑ x ) − (∑ xi )
i

While we ultimately want to fit the data to y = a + bx for a series i of N points, the theoretical development is much easier if we
i i

use y = c + b(x − x̄) where obviously a and c are related by a = c − bx̄ , and where x̄  is simply the mean of the x data,
i i
¯
¯ ¯
¯ ¯
¯

N
x . This designation of variables is better because "it turns out" that in this approach the parameters b and c are
¯¯
¯
nx = ∑ i=1 i

independent, their variances are independent, and the formulas to transform to the parameters a and b are simple. Again, we want to
minimize the sum of the squares of the residuals r = y − c − b(x − x) , R.
i i i
¯¯
¯

N N

2 ¯¯ 2
R = ∑r = ∑[ y − (c + b(xi − x̄))] (3.5.7)
i i

i=1 i=1

Now we require that R be as small as possible with respect to b and c by taking derivatives with respect to these two parameters
and setting the derivatives equal to zero.
∂R
¯¯
¯
( ) = −2 ∑ (yi − c − b (xi − x)) = 0 (3.5.8)
∂c
b

∂R
¯¯
¯ ¯¯
¯
( ) = −2 ∑ ((yi − c − b (xi − x)) (xi − x)) = 0
∂b
c

These two equations are readily rewritten as


¯¯
¯
∑ yi = ∑ c + ∑ b (xi − x) (3.5.9)

2
¯¯ ¯¯ ¯¯
∑ y (xi − x̄) = ∑ c (xi − x̄) + ∑ b (xi − x̄)
i

which, since ∑ (x i
¯¯
− x̄) =0 , simplify to
∑ yi
¯
¯¯
c = = −y (3.5.10)
N
¯¯
¯
∑ yi (xi − x)
b =
2
¯¯
∑ (xi − x̄)

The knowledge of c and b then allows us to determine the parameter, a, a = c − bx̄ = ȳ − bx̄ , a particularly simple relation. ¯
¯ ¯
¯ ¯
¯

The variance and standard deviation of N data points was introduced previously. The appropriate treatment of the parameters c and
b in the least squares procedure results in
2
s
V [c] = (3.5.11)
N
2
s
V [b] =
2
¯¯
∑ (xi − x̄)

3.5.2 https://chem.libretexts.org/@go/page/371024
where

2
1 2 1 2 Rmin
¯¯
¯
s = ∑ (yi − c − b (xi − x)) = ∑ (yi − a − b xi ) = (3.5.12)
N −2 N −2 N −2

Because c and b are independent quantities, it can be shown that


2
¯¯
¯
V [a] = V [c] + x V [b] (3.5.13)

2 2
s 2 s
¯¯
¯
= +x
2
N ¯¯
∑ (xi − x̄)

The standard deviations of a and b, σ and σ , are then obtained by taking the square roots of V[a] and V[b].
a b

Another statistic that is sometimes used is the correlation coefficient, r. It is given by the expression
N ∑ xi yi − ∑ xi ∑ yi
r = −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− (3.5.14)

2 2 2 2
√ (N ∑ x − (∑ xi ) ) (N ∑ y − (∑ yi ) )
i i

When the absolute value of r is close to unity the correlation between the y and x data is good, when it is close to zero the
correlation is poor. Because the correlation coefficient is sometimes hard to interpret, it is usually better to work with the standard
deviations of the regression parameters, σ and σ .
a b

Linear regression is a powerful tool that takes the guesswork out of obtaining best fit information. However, like all mathematical
tools, it must be used with caution. While the standard deviations and the correlation coefficient give indications of the goodness of
the fit, there is no substitute for graphing the data and looking at the result. Since the errors are assumed to be random, the yi values
should be scattered about the best fit straight line in a random fashion. The residuals should likewise be either positive or negative
(and sum to zero) without any pattern. If the data show a pattern of curvature, it is possible that the results do not conform to the
linear model proposed. A linear regression is not valid in this case. You should also be alert to the effects of systematic errors,
which may change the slope or the intercept without affecting the linearity of the data.
Linear regression of the type discussed above relies on a number of critical assumptions. Among the most important is that all the
uncertainty in each of the y values is the same. If the form of the equation required to provide a linear relationship produces an
independent variable that contains a significant uncertainty, the set of conditions that produced the minimum in the residuals of the
line may not be the set that actually minimizes the uncertainty in the data. In this case, a "non-linear" least squared technique such
as simplex is recommended.

 Note

In this course, you can use the MatLab lsq.m script to perform least squares analysis. The equations given above are included
in the lsq script.

3.5: Least Squares Linear Regression is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

3.5.3 https://chem.libretexts.org/@go/page/371024
3.5.1: Rejection of Discordant Data
There may be one or more yi values that deviate markedly from the trend of the others. One can use a Q-test criterion for
considering the rejection of one outlying datum. Arrange the residuals, ri, in order of increasing magnitude and use Equation 7-4
(Skoog, et. al., Ch. 7) and Table 7.3 (Skoog, et. al., Ch. 7) as before to decide whether the point can be rejected.

3.5.1: Rejection of Discordant Data is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

3.5.1.1 https://chem.libretexts.org/@go/page/371043
3.5.2: Weighted Least Squares Analysis
When measurements of the dependent variable incur different uncertainties, a weighted least squares technique should be used. In
this approach, each residual ri is divided by a factor proportional to its uncertainty ε . Then R, the sum of the squares of the yi

residuals, becomes
N 2 N
r 1
i 2
R =∑ =∑ (yi − a − b xi ) (3.5.2.1)
ε 2 ε 2
i=1 yi i=1 yi

Using the same minimization techniques described above, the slope and intercept for weighted least squares become
2
x yi xi xi yi
i
∑ ∑ −∑ ∑
2 2 2 2
εy εy εy εy
i i i i

a = (3.5.2.2)
2
2
x xi
1 i
∑ ∑ − (∑ )
2 2 2
εy εy εy
i i i

xi y xi y
1 i i
∑ 2
∑ 2
−∑ 2
∑ 2
εy εy εy εy
i i i i

b = (3.5.2.3)
2
2
x xi
1 i
∑ 2
∑ 2
− (∑ 2
)
εy εy εy
i i i

Weighted least squares has utility in a number of applications in the physical chemistry laboratory. One common situation arises
when a linear equation requires the log of a measured quantity to be the dependent variable. Consider for example, the
measurement of vapor pressure. Most pressure sensing devices will provide a constant error in the measurement of P. If the
dependent variable to be plotted is lnP, however, its uncertainty is not constant but decreases as P increases. Thus,
∂lnP εP
εlnP = εP = (3.5.2.4)
∂P P

An appropriate weighting factor in this case would be ε yi =


1

P
.
When you use Matlab for data analysis this semester, the command lsq(x,y) will be used for least squares analysis of x, y data, and
the command wtlsq(x,y) will be used for weighted least squares analysis. In the latter case, you will be asked for a weighting
factor, e. In either case, the program will output the least squares slope, y-intercept, correlation coefficient, standard deviation of
the slope, and standard deviation of the intercept.

3.5.2: Weighted Least Squares Analysis is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

3.5.2.1 https://chem.libretexts.org/@go/page/371045
3.5.3: References
1. Shoemaker, D.P.; Garland, C.W.; Nibler, J.W. Experiments in Physical Chemistry, 6th Edition, McGraw-Hill, New York, NY,
1996, pp 27-89.
2. Bettelheim, F. A., Experimental Physical Chemistry, W.B. Saunders Co., Philadelphia, 1971, Chapter 2.
3. Brownlee, K. A., Statistical Theory and Methodology in Science and Engineering, 2nd Edition, John Wiley and Sons, Inc., New
York, New York, 1965, pp 334-338.

3.5.3: References is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

3.5.3.1 https://chem.libretexts.org/@go/page/371046
3.6: Summary - General Guidelines For The Treatment of Experimental Errors
2.2.1. Before you leave the lab, remember to do the following:
(a) Make note of the uncertainties in each of the physical quantities you measure. Follow the suggestions in step (2.2.2) below.
Consult your TA concerning systematic errors that may arise from instrumental calibration and measurements.
(b) Be aware of the consequences of systematic errors, such as starting the clock too late in one of the kinetics experiments.
2.2.2. Estimation of uncertainties in mass, volume, pressure and temperature measurements
(a) Find and record the uncertainties associated with your equipment and instrumentation. This information is often listed right on
the device, or is easily obtained from the manufactures statement of equipment specifications. In the rare event that such
information is not available, you may assume that single experimentally measured values are uncertain in the last digit by at least ±
0.5 of the smallest unit or measurable digit. For example, if a digital voltmeter reads 5.432 V, and no other information is available,
the uncertainty may be estimated to be ± 0.005 V.
(b) The uncertainty of our analytical balances is either ± 0.0001 g or ± 0.0002 g, You can often find it stated on the individual
balance and associated with the letter d (for deviation). For example, the label might say d=0.1 mg.
(c) The position of the liquid level in a burette or a graduated 1 mL Mohr pipet can be estimated to ± 0.05 mL (half of the smallest
unit measurable). Initial and final readings must be made, however, so that the probable error in the volume measured will be,
−−−−−−−−−−−−−
2 2
√ (0.05) + (0.05) = 0.07 (3.6.1)

A similar analysis can be applied to pressure determination by measurement of mercury levels in a U-tube manometer.
(d) For non-graduated volumetric pipets, the error limits can be found right on the pipet itself, or in the " Comparison of Class A
and B Volumetric Glassware” link, found in the External Links section of your Sakai site or from other resources readily available
online. You might also look in the Baxter Scientific catalog found in the lab. This catalog is a useful source for tolerances and error
limits of thermometers, balances, pipets, burets, volumetric flasks and graduated cylinders.
2.2.3. Replicate or Single Measurements
(a) If you make a series of replicate measurements of the same quantity, you should report the mean of that quantity together with
its precision (standard deviation, variance or 95% confidence limits) calculated by statistical analysis. Unless otherwise requested
in the lab manual, use the standard deviation as a measure of precision. If the accepted value of the quantity has been reported in
the literature, you might also report the percent error or absolute error as a measure of the accuracy of your measurement(s). If you
feel that one of the replicate measurements should be rejected, use the Q-test to justify your decision.
(b) If you make a single measurement of a quantity, give your best estimate of the uncertainty in that measurement, whether
random or systematic. If appropriate (see 2.2.3(a) above), you should also report the accuracy of that measurement. The
suggestions given in Step (2.2.2) may help. Consult your TA if in doubt.
2.2.4. Graphical Analysis
(a) If you use graphical analysis to determine the value of any quantity, you must make a decision whether or not to use a weighted
least squares regression analysis. In most cases, a non-weighted analysis (using the command lsq in Matlab) will suffice. If you
perform a weighted least squares analysis, you will have to decide on what weighting factor, e, to use. If the data deviate
systematically from linearity, a linear regression analysis is meaningless, and should not be performed.
(b) Outliers may be rejected only on the basis of a Q-test performed on the residuals.
(c) The standard deviation may be used as a measure of the uncertainty in the slope and intercept calculated by linear least squares
regression.
(d) In most cases, you will use the values of the slope and/or y-intercept to calculate a physical quantity. In either case, you must
perform a propagation of error treatment to determine the uncertainty in that physical quantity.
2.2.5. Propagation of Error
Whether or not you perform a graphic analysis, you will use experimentally determined variables to calculate one or more physical
quantities. In each case you must estimate the uncertainty in each measured variable and then propagate that uncertainty through to

3.6.1 https://chem.libretexts.org/@go/page/371049
the final result. In a few simple cases you may take advantage of Equations (55) - (59) (Shoemaker et. al., 1996). Usually however,
you will use Equation (54) (Shoemaker et. al., 1996), and take partial derivatives to calculate the probable propagated error.
Physical constants, such as the gas constant, R, may be assumed to be error free.
2.2.6. Reporting the Final Numerical Result(s)
(a) Each numerical result should be reported together with the associated uncertainty, which is usually the probable propagated
error. In addition, and where possible, you should report the accuracy of your result by comparison with literature values and
calculation of the absolute error.
(b) In reporting your results, ensure that the use of significant figures is appropriate to the uncertainty in those results. Use the
conventions for significant figures and for rounding off data described earlier in this text.
2.2.7. Qualitative Discussion of Uncertainty in an Experiment
In some cases, it may be difficult to quantitatively estimate the uncertainty in an experimental variable. In this instance, you may be
asked to provide a qualitative discussion of the sources of uncertainty in the experiment. This discussion should be as thorough as
possible. You should include "reasonable" estimates of the uncertainties in measured parameters together with a discussion of how
these uncertainties will affect the final result. Any estimates of uncertainty should be substantiated by hard facts!!

3.6: Summary - General Guidelines For The Treatment of Experimental Errors is shared under a not declared license and was authored, remixed,
and/or curated by LibreTexts.

3.6.2 https://chem.libretexts.org/@go/page/371049
3.6.1: Error Analysis in a Nutshell
And finally, if you read nothing else .....
2.2.1.1. Estimate the uncertainty in each experimentally determined variable. Usually this involves one or more of the following:
(a) Calculate the standard deviation of a quantity determined from the mean of a series of replicate measurements.
(b) Find and record the uncertainties associated with your equipment and instrumentation. This information is often listed right on
the device, or is easily obtained from the manufactures statement of equipment specifications. In the rare event that such
information is not available, you may assume that single experimentally measured values are uncertain in the last digit by at least ±
0.5 of the smallest unit or measurable digit. For example, if a digital voltmeter reads 5.432 V, and no other information is available,
the uncertainty may be estimated to be ± 0.005 V.
(c) Perform a simple propagation of error analysis (usually equation (54) (Shoemaker et. al., 1996)) if a quantity (e.g., volume,
pressure) is determined from a difference of two or more readings.
(d) Extraction of a Physical Quantity from Graphical Analysis: Calculation of the standard deviation of the slope (or intercept) is
followed by a propagation of error treatment of the functional form of that slope (or intercept) to give the uncertainty in the
physical quantity to be extracted. If the data deviate systematically from linearity, a linear regression analysis is meaningless and
should not be performed.
(e) If the source or estimation of uncertainty seems more obscure, -- CONSULT YOUR TA BEFORE YOU LEAVE THE LAB.
2.2.1.2. Once you have determined the uncertainties associated with each individual measurement, perform a propagation of error
treatment to estimate their combined effect on the uncertainty of the physical quantity of interest. Use equation (54) or some
combination of equations (55) - (59) (Shoemaker et. al., 1996). Physical constants such as the Gas Constant, R, may be assumed to
be error free.
2.2.1.3. Report your final result(s), the probable propagated error associated with that result, and the percent error or absolute error
if appropriate literature values are available.

3.6.1: Error Analysis in a Nutshell is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

3.6.1.1 https://chem.libretexts.org/@go/page/371061
3.7: Homework Problems
Each practice problem below requires that you correctly apply propogation of errors to report the correct value and its preceision
(including the rror). You should attempt all seven problems prior to our schedule lab meeting on this topic (submit prior to the
meeting on Sakai). This will prepare you to participate in the lab meeting by presenting problems on the board, answering
questions from your peers, and asking your own questions. Your answers to these problems may be hand-written or types.
However, any submissions that are deemed illegible, unclear, or poorly organized will not be accepted, must be re-submitted as late
reports, and will receive the usual late grade penalty of a 10% deduction per calendar day the assignment is late.

 Note
(1) You may use whatever technology (calculators, software, etc.) you like to aid in working these problems. However, it is
important to show how your calculations are set-up and, where appropriate, give a sample calculation. Simply giving the final
answer is not acceptable.
(2) In reporting all results, you should always round your reported answer to be consistent with your estimated uncertainty. For
example: a calculated result of 56.2189 ± 0.0386 would be reported as 56.20 ± 0.04 where 0.04 is the standard deviation, 95%
confidence level, or propagated uncertainty.

PROBLEM 1
A student makes five independent measurements of each of the linear dimensions of a box. The results are tabulated as follows:

Length (cm) Width (cm) Height (cm)

1.010 5.376 2.001

1.015 5.339 2.105

1.013 5.385 1.985

0.999 5.375 1.989

1.009 5.369 2.008

A. Calculate the average values of length, width, and height, and calculate the variance in each of these three quantities.
B. Calculate the standard deviation of the length, width, and height.
C. Calculate the 95% confidence limits to the measurements of length, width, and height.
D. Calculate the 95% confidence limits to the average volume obtained from the independent measurements of length, width, and
height.

PROBLEM 2
A student weighs a sample using an analytical balance that has an estimated precision of ±0.5 mg. From the data below, calculate
the weight of the sample and the range of uncertainty in the net weight.

Mass (g)

Sample + Tare 5.1647 ± 0.0005

Tare 2.3712 ± 0.0005

PROBLEM 3
An investigator measures the photon emission in the reaction

+ − ∗
Na +O → Na +O → Na +O

in which the emitted photon (at 589 nm) is counted by means of a photomultiplier. Here are some typical "run" data

3.7.1 https://chem.libretexts.org/@go/page/371063
Signal + Background Background
Run #
(counts/sec) (counts/sec)

1 1500 1000

2 1200 1250

3 1375 1100

4 1425 900

5 1280 1200

6 1490 1363

Evaluate the data using statistical methods only and calculate:


A. The average signal value.
B. The variance in the average signal.
C. The 95% confidence limits to the average signal value.
D. Using the Q test, determine whether the results from any one run can be rejected as discordant data at the 90% confidence level

PROBLEM 4
A. Fit a straight line:
y = a + bx

to the following set of data points by the linear least-squares method


x 26 30 44 50 62 68 74
y 92 85 78 81 54 51 40
and calculate the coefficients a and b using the equations 3.1.5 and 3.1.6 from Section 3.1.
B. Plot the x,y data and the least-squares line, and submit this with your report.
C. Calculate the standard deviation of the coefficients, a and b, using Equations (3.1.11), (3.1.12) and (3.1.13) from from Section
3.1.
D. Calculate the correlation coefficient, using equation (3.1.14) from Section 3.1. What is the significance of the sign of r?
E. Apply the Q-test to the y-residuals (see Equation (3.1.2) from Section 3.1), to determine whether there is a statistical basis for
rejection of any of the data points as an outlier, (a) at the 96% confidence level (b) at the 90% confidence level.
F. If a data point is rejected, repeat the regression analysis (steps A-D) using the remaining data.
G. Comment quantitatively on the effect of data removal on the regression analysis.

PROBLEM 5
In a chemistry lab, students obtain the following data and are asked to find the molecular weight (M) of methanol using PV = nRT.
Determine the uncertainty in the result using a propagation of error treatment.

Pressure (methanol) Temperature mass (methanol)

Trial 1 72.5 ± 0.5 mmHg 295.10 ± 0.05 K 0.1331 ± 0.0002 g

Trial 2 72.0 ± 0.5 mmHg 295.05 ± 0.05 K 0.1329 ± 0.0002 g

The volume (V) in both trials was 1.05678 liter, with no uncertainty given. You must decide on the uncertainty in V. Hint: See
2.2(a) in the Summary section of this. If necessary, perform appropriate conversions of the units for each quantity, including R, so
that the final value of M is reported in the correct units.

PROBLEM 6

3.7.2 https://chem.libretexts.org/@go/page/371063
The rate constant for the reaction
2DI → D2 + I2

is measured at two temperatures with the given estimated uncertainties:


k660 = (1.2 ± 0.2) × 10
−3 liter

mol sec
at 660 ± 2K
k660 = (1.8 ± 0.2) × 10
−3 liter

mol sec
at 720 ± 2K
A. Calculate the activation energy Eact for the reaction using
k2 Eact 1 1
ln ( ) = − ( − )
k1 R T2 T1

B. Assuming the above equation is a correct expression for Eact, what is the uncertainty in Eact?

PROBLEM 7
In the method of Clement and Desormes, the heat capacity ratio of an ideal gas is determined from the manometrically determined
values of P1, the gas pressure at temperature T1; P2, the gas pressure immediately after a reversible adiabatic expansion when the
temperature is T2; and, P3, the gas pressure after the gas has warmed back up to temperature T1. The equation used to calculate the
heat capacity ratio is
C
V
+1
CP R
γ = =
C
CV V

where
P
3
P2 [1−( )]
CV P
1
=
R P3 − P2

Calculate the uncertainty in the heat capacity ratio γ, given the following data:
P1 = 763.1 ± 0.2 mmHg
P2 = 757.0 ± 0.1 mmHg
P3 = 758.7 ± 0.2 mmHg

3.7: Homework Problems is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

3.7.3 https://chem.libretexts.org/@go/page/371063
CHAPTER OVERVIEW

4: Absorption Spectrum of Conjugated Dyes


You will measure the absorption spectra of a series of conjugated dyes and then use two theoretical models, the Particle-In-A-Box
Approximation and, later, the Hückel Molecular Orbital Approximation to explain the observed spectra.

The general approach to this experiment is adapted from D. P. Shoemaker, C. W. Garland, and J. W. Nibler, Experiments in
Physical Chemistry, 6th edition, McGraw Hill Co. Inc, NY, 1996, p378. The experiment has been revised to suit our laboratory
environment.
4.1: Prelaboratory Assignment
4.2: Introduction
4.3: Experimental (Part I)
4.4: Hückel Calculations Using Hulis (Part II)
4.5: References
4.6: Appendix A - Use of the Agilent 100 Series UV-Vis Spectrophotometer
4.7: Appendix B - The Hückel Approximation

4: Absorption Spectrum of Conjugated Dyes is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

1
4.1: Prelaboratory Assignment
In addition to the general pre-lab assignment, please complete the specific tasks listed below:

 Pre-Lab assignment for Absorption Spectra of Conjugated Dyes


1. Find structures and literature values for λ
max for the four dyes to be studied in this experiment:
Dye #1: 1,1’-diethyl-2,2’-cyanine iodide
Dye #2: 1,1’-diethyl-2,2’-carbocyanine chloride
Dye #3: 1,1’-diethyl-2,2’-dicarbocyanine iodide
Dye #4: 1,1’-diethyl-4-4’-carbocyanine iodide

Some helpful resources include


The Aldrich chemical web site (http://www.sigmaaldrich.com), Fisher, N.I.,
Hamer, F.M., A Comparison of the Absorption Spectra of Some Typical Symmetrical Cyanine Dyes, Proceedings of
the Royal Society of London, Series A, 154 (883), 1936, pps. 703-723, and/or find other sources for these values.

2. Define the terms, highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO).
Describe what you understand about HOMO and LUMO.

3. Read McQuarrie, D.A., and Simon, J.D., Physical Chemistry: A Molecular Approach, University Science Books, CA, 1997
(Section 3.4). Read also, Cotton, F.A., "Chemical Applications of Group Theory", Wiley, New York, 1971, Section 7.1 and
pp. 366-367

The information regarding general expectation for preparation for the lab meetings is copied below for reference:

Pre-Lab Assignment:
You should prepare for lab meetings by reading the appropriate module in this manual, and completing the specified pre-laboratory
assignments. The pre-laboratory assignments are designed to prepare you to operate efficiently and safely during the experimental
module. You will not be allowed to conduct experiments if the pre-lab assignment is not completed before your lab meeting. If
there is an extenuating circumstance that prevents your preparation, please inform your instructors as soon as possible.
For a given experiment, you should complete the assigned reading and write out the specified pre-lab assignments. The pre-lab
assignment must be submitted as a single pdf document on Sakai at least 24 hours before the scheduled lab meeting. Late pre-
laboratory work will not be accepted and will receive a score of zero. Your TA will evaluate the work and return it no later than the
beginning of your lab section. Your preparedness will also be evaluated by your answers to oral questions pertaining to the
laboratory from the TA and other lab instructors. The pre-lab and preparation accounts for 10 points of each experiment requiring a
pre-lab.
For every experimental module there is a general pre-lab assignment. It is described here:
General Pre-lab assignment:
Write a brief introduction to the experiment and the experimental plan using the following outline:
I. Experiment Introduction: In your own words, briefly describe the underlying theory and explain the principle upon which the
experiment isbased. Keep it concise; lengthy discussions or derivations of equations are unnecessary.
a. Background (importance/relevance): A brief description and justification of the importance of the topic, with references to
source material. What is the topic and why is it interesting and important? This is meant to "hook" the reader so they want
to read about the topic.
b. Theory: A concise paragraph or more describing the general theory for the method/technique used; this section should
contain a reference to the appropriate pages of the textbook, or other source material from the literature.

4.1.1 https://chem.libretexts.org/@go/page/371214
c. Goals/Purpose: Brief statement of purpose, which should indicate what is analyzed and the technique used. Limit to three to
five sentences.
II. Experimental Plan: a short summary of the specific things that you would need to know in order to do the experiment without
having access to your lab manual, with the understanding that you will always have a TA available to teach you how to use the
instrument. The best way to approach this task is to thoroughly read the procedure for the experiment. While reading, take notes
on the specifics of the work. What chemicals will you use? If you have to make solutions, what solutions do you need to make
and of what concentration? How will you prepare them? What instrument will you use, and what specific instrument settings
will you need to input (sometimes you will find this information in the Appendices that discuss the use of the instruments)?
How many runs of each sample will you do; and over what wavelength, or temperature, etc. range? You do not need to write
down how to use the instrument; we will teach you that. Bottom Line: we want you to write down only that pertinent
information about samples, sample preparation, instrument and instrument parameters, numbers of runs, and other things, that
you will need to know in order to do your work.
Pre-lab assignments should be typed. If you prefer to write out equations by hand, they should be neat and inserted as an image in
the appropriate place. If you need advice for how to insert equations or images into your documents, please ask for help at least 48
hours before the assignment is due; an in-person (or zoom) meeting is most appropriate for this issue.

Always Wear Protective Clothing:


Please dress appropriately for a laboratory environment - this is considered an important part of preparation for the laboratory
period (see also Laboratory Safety below). You must wear close-toes shoes and clothing must cover your body from your upper
arms to your knees. You should also wear protective eyewear. If lab attire is not jiving with your wardrobe, please purchase a lab
coat and bring it with you so that you can cover up in lab. No student should be in lab without wearing protective eyewear and
appropriate protective clothing.

4.1: Prelaboratory Assignment is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.
Introductory Details by W. R. Fawcett, John Berg, P. B. Kelley, Carlito B. Lebrilla, Gang-yu Liu, Delmar Larsen, Paul Hrvatin, David
Goodin, and Brooke McMahon is licensed CC BY-NC-SA 4.0.

4.1.2 https://chem.libretexts.org/@go/page/371214
4.2: Introduction
The Particle-In-A-Box approximation
Electrons in the π-electron system of a conjugated aromatic compound are not restricted to specific nuclei but are free to move
throughout the system. In a linear conjugated system the potential energy of the electrons will vary along the chain, being lowest
near the nuclei and highest between them. To a first approximation, however, it can be assumed that the potential energy of the
electrons is constant along the length of the chain and increases to infinity one bond distance beyond the end atoms of the chain.
The electrons may be visualized as being contained in a one-dimensional potential well, of length ℓ , and occupying orbitals of
different energy within the well. The wave functions representing the orbitals and the energy levels of the orbitals can be obtained
by solving the Schrödinger equation for this system.
For a one-dimensional system the Schrödinger equation takes the form
2 2
∂ ψ(x) 8π m
+ (E − V ) ψ(x) = 0 (4.2.1)
2 2
∂x h

where ψ (x)
are wavefunctions that satisfy the equation (eigenfunctions), E is the energy of the particle of mass m , V is the
potential energy, and h is Planck's constant.
Since the potential energy, V , is a constant within the well of length ℓ , where ψ must vanish for x ≤ 0 and x ≥ ℓ , the regions
(x)

where V = ∞ , it may be conveniently set equal to zero when deriving solutions of the equation. The equation then reduces to
2 2
−h ∂ ψ(x)
= E ψ(x) (4.2.2)
2
8π m ∂x2

Solutions of this equation are ψ (x) = A ⋅ sin(


nπx


) , where A is a constant and n = 1,2,3,... etc.

Verification that this wave function is indeed a solution of the Schrödinger equation may be obtained by substitution into the left-
and right-hand sides of the above equation.
2 2 2
−h −n π nπx
Left-hand side = ( ) A ⋅ sin( )
2 2
8π m ℓ ℓ

2
−h nπx
= A ⋅ sin( )
2
8mℓ ℓ

nπx
Right-hand side = EA ⋅ sin( )

2 2

The function ψ(x) is a solution of the Schrödinger equation when the energy E of the system has the values E =
n h
2
. The
8mℓ

allowed energy levels are depicted diagrammatically below.

Figure 4.2.1 : The Particle-in-a-Box Model.

4.2.1 https://chem.libretexts.org/@go/page/371215
The ground state configuration for the conjugated system of a molecule is derived by allocating two electrons to each of the energy
levels in accord with the Pauli exclusion principle. Thus, a molecule with Nπ electrons will have the N/2 lowest levels filled and all
the higher levels empty. When the molecule absorbs UV-vis light, the energy absorbed is from an electron in the HOMO (where n =
N/2) becoming excited to the LUMO where n = (N/2) + 1. The energy change, ∆E, for this transition is
ΔE = E N
−E N
( +1) ( )
2 2

2 2
N 2 N 2
[( ) + 1] h ( ) h
2 2
= −
2 2
8mℓ 8mℓ

2
(N + 1) h
=
2
8mℓ

Since ΔE = hν =
hc

λ
where c is the speed of light and λ is the wavelength, then
2
8mc ℓ
λ = (4.2.3)
h (N + 1)

Thus if ℓ is known, this relation can be used to calculate the wavelength maximum of the band representing the transition from the
ground state to the first excited state in a π-electron system.
This free electron model can be applied to carbocyanine dyes which may be represented, for example, by the resonance forms
shown below for the 1,1'- diethyl-4,4'-carbocyanine (cryptocyanine) cation.

Figure 4.2.2 : The 1,1'-Diethyl-4,4'-carbocyanine (cryptocyanine) cation. (CC-NC-BY-DUKE CHEM)


The π-electrons in the molecule are regarded as a collection of independent particles confined to a box, whose wavefunctions are
the square-well functions. Once again, these electrons are distributed among the wavefunctions and energy levels according to
Pauli principles, so that the absorption according to this model will follow equation 4.2.3. To apply the model we need to
determine N and ℓ .
In the chain between the two nitrogen atoms, each carbon atom and each nitrogen atom forms three σ-bonds. The extra valence
electron of each carbon atom and the three remaining electrons contributed by the two nitrogen atoms form a mobile cloud of π-
electrons along the length of the chain (above and below the plane of the chain). In the chain all the bonds can be considered
equivalent with bond order 1.5, (similar to the C-C bonds in benzene). If we denote the number of carbon atoms in the polymethine
chain (between the two nitrogen atoms) by p, then N, the number of π-electrons, is p + 3. The total length of the potential well is
taken to be the length of the shortest carbon chain between the nitrogen atoms plus one bond length at each end. Thus, ℓ = (p + 3)a
where a = 1.39 x 10–8 cm, the mean distance between atoms in the chain. Using these values for N and ℓ in equation 4.2.3 yields:
2 2
8mca (P + 3)
λ = (4.2.4)
h P +4

If polarizable groups are present at the ends of the chain, then the potential energy of the π-electrons in the chain does not rise so
steeply at the ends. This results in lengthening of the box length, ℓ , which can be incorporated into equation 4.2.4 through the
addition of a new "lengthening" parameter, γ, resulting in

4.2.2 https://chem.libretexts.org/@go/page/371215
2 2
8mca (P + 3 + γ)
λ = (4.2.5)
h P +4

where γ should be a constant for a series of dyes of a given type and should lie between 0 and 2. If such a series is studied
experimentally, the empirical parameter γ may be adjusted to achieve the best fit to the data.

4.2: Introduction is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

4.2.3 https://chem.libretexts.org/@go/page/371215
4.3: Experimental (Part I)

The ELN Template for this part of the experiment is here: ELN01_1_dyes2022.docx

Materials
Equipment:
UV/Vis spectrometer (you will use either an Agilent Cary 60 Spectrometer or an Agilent Cary 100 Spectrometer)
Cuvettes (you will use glass cuvettes, which are optically transparent to visible (but not to UV) light)
Reagents:
Reagent grade methyl alcohol (methanol)
Stock solutions of four polymethine dyes dissolved in methanol (solutions are approximately 10–4 M):
1,1'-diethyl-2,2'-cyanine iodide
1,1'-diethyl-2,2'-carbocyanine chloride (also known as pinacyanol chloride)
1,1'-diethyl-2,2'- dicarbocyanine iodide
1,1'-diethyl-4,4'-carbocyanine iodide (also known as cryptocyanine)

 Note

The cyanine dyes used in this experiment are toxic irritants. Solutions of these dyes in methanol must be handled wearing
gloves at all times. Perform all dilutions in a fume hood. Read the safety data sheets for these dyes
(http://www.sigmaaldrich.com) prior to the lab.

Procedure
1. Perform instrument diagnostics tests using the Validate program. (Appendix A has instructions for using the instrument and
software)
2. Fill four glass cuvettes approximately 4/5 full with methanol (about 2.5 - 3 mL). Then add about 5-15 drops of each dye to a
series of glass cuvettes. Add enough dye that you can see a light color for each solution.
3. Use the UV-vis to record the spectrum of each dye. If needed, adjust the solution concentration of one or more of the dye
samples until the peak maximum absorbance reading is approximately 0.2 to 1.0 absorbance units.
4. Annotate the resulting overlay display of all four dyes, noting which dye is associated with each trace and the associated λ max

value.

Calculations and Discussion Questions


1. Carefully examine all spectra (absorbance versus λ ) and determine λ max , the peak wavelength for each dye. Make any
necessary corrections indicated by the validation procedure (the wavelength calibration).
2. What is the uncertainty in your experimental peak wavelength determinations?
3. Compare your results with literature values for λ max .10-12,14
4. For each dye, calculate the energy (eV and kJ/mol) of the experimentally observed λ max .
5. You will use a Matlab script called dye to determine the "best" value of γ This script will prompt you for information about the
four dyes, and the experimentally measured values of λ max . It will solve the particle in a box equation for λ , separately for each
dye, over the range of γ from 0.00 to 2.00 in steps of 0.01. It will sum the absolute deviation between the calculated λ and your
experimentally measured λ max for each value of γ, for each dye. It will output the one value of γ that produces the smallest
summed deviation over the series of the four dyes. That is, the script will calculate the value of gamma which gives the best fit
between peak wavelengths calculated using the free electron model and peak wavelengths measured experimentally. To use this
script, start up a Matlab session and type the dye command:

dye

4.3.1 https://chem.libretexts.org/@go/page/371570
Then simply enter the information as prompted by the script.

Use the output value of γ to calculate a value for λ max for each of the four dyes using the equation
2 2
8mca (P + 3 + γ)
λ = (4.3.1)
h P +4

(Remember that the value of the average bond length of a carbon-carbon bond, a, is 1.39 x 10–8 cm the mass of the electron, m,
is 9.1 x 10–31 kg.)
6. Calculate the absolute error between your experimental (from the spectrum) values and the literature values. Calculate both the
absolute error and the percent deviation between your experimental (from the spectrum) and your calculated theoretical
(Particle-in-a-box Model) λ results.
max

7. Tabulate your data in the format shown below.

WAVELENGTH (nm)

λmax λmax λmax


Dye # # drops used
Literature* Experimental Particle in a Box

Absolute Error Absolute Error


Value (nm) Value (nm) % Error
(nm) (nm)

1 ±

2 ±

3 ±

4 ±

* reference for your literature values goes here

8. Use this table as the basis for discussion of your results.


a. How well do your λ values from your spectra agree with literature λ
max ? When thinking about this look to see if the
max

absolute error between your experimental result and the literature is equal to or less than the uncertainty (± value) of the
measurement.
b. How well does the Particle in the Box model reproduce the experimental results? When thinking about this look first at your
experimental result to see how varies as a function of p. Does the PIB result give the same trend? Looking at the
individual dyes, do your PIB results agree with your experimental results within their absolute errors?
c. Are the results for each dye equally satisfactory? If your answer is no, explain why you think this is so. How is the PIB
model useful and what might be done to improve the results?
9. Calculate the wavelength (in nm) and minimum energy (eV) of the molecule 1,3,5,7-octatetrene (
CH =CHCH=CHCH=CHCH=CH ) using the particle-in-a-box model. What color does the molecule absorb from white
2 2

light? What color does it then appear to the eye?

4.3: Experimental (Part I) is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

4.3.2 https://chem.libretexts.org/@go/page/371570
4.4: Hückel Calculations Using Hulis (Part II)

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.
1. Perform 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.
2. Calculate the energy difference, ∆E, between the HOMO and LUMO levels (in terms of β) for each dye.
3. Plot the HMO ∆E's (abscissa) versus the experimentally observed ∆E's (kJ/mol) (from step (4) above) (ordinate). From the
least-squares slope of the resultant line, determine a semiempirical value for the spectroscopic β. Compare your value with the
accepted range for the spectroscopic exchange integral, – 250 to – 350 kJ/mole.9 Report the uncertainty in β.
4. Use your value of β 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 → LUMO
transition for each dye, using the HMO model.
5. Calculate the percent error between your calculated HMO results your experimental values. Tabulate your data in the format
shown below.

WAVELENGTH (nm)

λmax
Dye # λmax Experimental λmax Particle-in-box λmax HMO
Literature*

Absolute Error % Error % Error


Value (nm) Value (nm) Value (nm)
(nm) (from exp.) (from exp.)

* reference for your literature values goes here

6. 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?
7. 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 β 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

 Note
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.

1. Double-click on the Hulis icon to Start the Program


2. 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.

4.4.1 https://chem.libretexts.org/@go/page/371632
Figure 4.4.1 : The Hulis Main Screen. (CC-NC-BY-DUKE CHEM)
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.

3. Continue adding carbon atoms as appropriate to form the template molecule. When you are done your structure will look like
this:

Figure 4.4.2 : Dye 1. Note that for the purposes of a Hückel calculation the ethyl groups are omitted (CC-NC-BY-DUKE
CHEM)

4.4.2 https://chem.libretexts.org/@go/page/371632
 Note

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.

4. 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.

Figure 4.4.3 : 1,1'-Diethyl-2,2'-Cyanine Iodide. The ethyl groups are omitted to simplify the calculations. (CC-NC-BY-DUKE
CHEM)
Looking at the orbital diagram in the right-hand (orange) pane of the program screen you will see the occupied and unoccupied
orbitals, as calculated by the simple Hückel approximation. You can see that the 21 pi-electrons in the molecule occupy the first
11 pi-orbitals; the HOMO to LUMO transition occurs between the 11th and 12th pi-orbital.
5. Click on the Results button, and a Results window appears with significant information about the Hückel calculation. Please
note the following things:
a. 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.

Number of Coulomb Exchange


Atom Bond
p-Electrons Integral Integral

C 1 0.00 1.00

N Pyridine 2 1.37 0.89

N Pyridinium 1 0.51 1.02

b. 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.

4.4.3 https://chem.libretexts.org/@go/page/371632
Figure 4.4.4 : Orbital Energies for Dye 1. (CC-NC-BY-DUKE CHEM)
Notice that in this case the energy of the HOMO is α + 1.2737\Beta; and, the energy of the LUMO is α − 0.4390β.

 Note

NOTE: you can change the number of digits in the display from Preferences, under the File menu.

Close the Results Box.

c. 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.

Figure 4.4.5 : The Highest Occupied Molecular Orbital in Dye 1. (CC-NC-BY-DUKE CHEM)

4.4: Hückel Calculations Using Hulis (Part II) is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

4.4.4 https://chem.libretexts.org/@go/page/371632
4.5: References
References
1. Castellan, G. W., Physical Chemistry, Addison-Wesley, Reading, Mass., 1972, Chapters 19-22.
2. Barrow, G. Physical Chemistry, 3rd ed., McGraw-Hill, New York, NY, 1973, Chapters 3 and 12.
3. Shoemaker, D. P., Garland, C. W., and Nibler, J. W. Experiments in Physical Chemistry, 6th ed., McGraw-Hill, NY, 1996,
Chapter XIV, Experiment 34.
4. Farrell, T.T., J. Chem. Ed., 1985, 62, 351.
5. Olsson, L.F., J. Chem. Ed., 1986, 63, 756.
6. Kuhn, H., J. Chem. Phys., 1949, 17, 1198.
7. Herzfeld, K.F.; Sklar, A. L. Rev. Mod. Phys., 1942, 14, 294.
8. Silbey, R. J., Alberty, R.A., and Bawendi, M.G., Physical Chemistry, 4th ed., Wiley, NY, 2005, Chapter 9, Section 9.6; Chapter
11, Section 11.7.
9. Cotton, F.A., "Chemical Applications of Group Theory", Wiley, New York, 1971, Section 7.1 and pp. 366-367.
10. Suzuki, H., "Electronic Absorption Spectra and Geometry of Organic Molecules", Academic Press, New York, 1967, Chapter
17, p. 374.
11. Aldrich Catalog, Aldrich Chemical Company, Inc., 1992., found in FFSC 1113 or online.
12. Venkataraman, K., "The Chemistry of Synthetic Dyes", Academic Press, Inc., NY, 1952, Volume II, Chapter XXXVIII, p. 1155.
13. Streitwieser, A., "Molecular Orbital Theory for Organic Chemists", John Wiley & Sons, Inc., NY, 1967, Chapter 2.
14. Okawara, M., Kitao, T., Hirashima, T., and Matsuoka, M., Physical Sciences Data 35: Organic Colorants, Kodansha, Tokyo;
Elsevier, Amsterdam-Oxford- New York-Tokyo, 1988, pp. 318-321.
15. Taubmann, G., J. Chem. Ed., 1992, 69, 96.
16. Atkins, P.W. Physical Chemistry, 5th ed., W. H. Freeman, NY, 1994, Sections 12.1, 14.9, or, 6th ed., W. H. Freeman, NY, 1997,
Sections 12.1, 14.9.
17. Engel, T., Drobney, G and Reid, P., Physical Chemistry for the Life Sciences, Pearson Prentice Hall, NJ, 2008 (Section 14.3).

4.5: References is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

4.5.1 https://chem.libretexts.org/@go/page/371574
4.6: Appendix A - Use of the Agilent 100 Series UV-Vis Spectrophotometer
Here is the Basic procedure for using the Agilent 100 Series UV-Visible spectrophotometer in the Absorption Spectrum of
Conjugated Dyes experiment.
There are two sample holders in the instrument. The Sample is placed in the front holder; and, the Reference is placed in the back
holder.
With both sample and reference holders empty and the cover closed, turn on the instrument with the switch at the front, lower left.
Wait two minutes before proceeding further.
3.5.1. Checking the Wavelength Accuracy of the UV-Vis Spectrometer

3.5.1.1 Load the Validation Software by clicking the icon from the Windows main screen. After a few moments the software
screen will appear and will begin to communicate with the instrument. You can follow the progress by monitoring the display at the
lower left of the screen. When the traffic light in the upper center turns green, the system is ready. You will see absorbance in the
upper left and wavelength in the upper right.

Figure 4.6.1 : The Validate Main Screen. (CC-NC-BY-DUKE CHEM)

3.5.1.2. Click the button and deselect, using the button, all tests except Wavelength Accuracy - Ho2O3. The screen
should look as shown here.

Figure 4.6.2 : Selecting the Holmium Oxide Validation. (CC-NC-BY-DUKE CHEM)

4.6.1 https://chem.libretexts.org/@go/page/371578
3.5.1.3. Click the button and, when the setup screen appears, place your names in the Name line; and, the name of the
test in the comments box. In the Options section select Auto Print, Results and Graph. Data Form and Company Logo should not
be checked. % Page Height set at 50 is fine.

3.5.1.4. Click to return to the main screen. Place the holmium oxide filter in the sample holder of the instrument, leave

the reference holder empty, and close the cover. Click to begin the holmium oxide wavelength accuracy test. As the test
proceeds the spectrum of holmium oxide will be displayed. At the conclusion of the test, a report will be generated on the screen.
Include these data in your laboratory report. Check to be sure that the difference between the expected and observed peak positions
is one the order of ±0.2-0.4 nm. If it is greater than that, consult your TA or your Instructor before proceeding further.

Figure 4.6.3 : Results of the Holmium Oxide Wavelength Accuracy Test. (CC-NC-BY-DUKE CHEM)
Exit the Validation software and return to the Windows Desktop.
3.5.2. Absorption Spectrum of Conjugted Dyes

3.5.2.1. Load the Scan Software by clicking the icon on the Windows main screen. In a few moments the scan software main
screen will appear and will begin to communicate with the instrument. You can follow the progress by monitoring the display at the
lower left of the screen. When the traffic light at the top center of the screen turns green, the system is ready. You will see
absorbance at the upper left and wavelength at the upper right of the screen.

Figure 4.6.4 : The Scan Module Main Screen. (CC-NC-BY-DUKE CHEM)

4.6.2 https://chem.libretexts.org/@go/page/371578
3.5.2.2 Click the button and wait for the setup screen to appear. Select, first, the Cary tab to begin setting parameters for
your sample runs. Under Xmode, the Mode should be set to nanometers and your scan range from 750 to 400. Under Ymode, the
Mode should be set to absorbance and the range from –0.05 to 1.0. The Ave time should be 0.1 seconds and the Data interval
should be set to 1 nm. This will yield a Scan Rate of 600 nm/min.

Figure 4.6.5 : The Scan Setup Screen. (CC-NC-BY-DUKE CHEM)


3.5.2.3. Select the Options tab on the setup screen and set the Spectral Band Width (SBW) to 2 nm. Select Double Beam Mode.
Select UV-Vis as the source and be sure the Source changeover is set to 350 nm. Under Display options you can choose Overlay
Data.

Figure 4.6.6 : The Options Tab of the Scan Setup Screen. (CC-NC-BY-DUKE CHEM)
3.5.2.4 Finally, select the Reports tab. Put your names in the Name box and the title of the experiment in the Comment box. In the
Options box, check Data Form, Graph and Focused trace. Do not check Auto print, or any of the other options. % Page Height set
at 50 is fine. Under Peak Table, put a check in All peaks. Be sure the Include XY pairs table is not checked. Click the Peak
Information button, and set the Threshold to 0.200. Put Peaks in the Peak Type box. Click to close the Peak information
box, and again to return to the main screen.

4.6.3 https://chem.libretexts.org/@go/page/371578
Figure 4.6.7 : The Reports Tab of the Scan Setup Screen. (CC-NC-BY-DUKE CHEM)
3.5.2.5. Place your first sample in the sample compartment and your reference in the reference compartment and close the cover.

Click to begin the run. You will be prompted for a filename. Since you will be running four samples, give the file a
general name, for example Dyes with the date (many sample spectra can be stored in a single file). You will then be prompted for a
sample name. Place the name of the sample (Dye 1, 2, 3, or 4) in this box. As each run proceeds you will see the spectrum
generated in real time on the screen. Make sure the maximum absorbance of each sample is between 0.4 and 0.8; adjust the dye
concentration to get the proper absorbance. At the conclusion of each run, place the next sample in the sample compartment and
initiate the new scan.

 Note

If you didn't select Overlay Data, or if you want to add or remove traces from the overlaid runs on the screen, click the
button at the upper left of the main screen. Select files to add or to remove by checking or removing cheks from available
traces shown in the list. Then click OK. Your changes will be made to the main screen.

4.6.4 https://chem.libretexts.org/@go/page/371578
Figure 4.6.8 : The Absoprtion Spectra of the Dyes. (CC-NC-BY-DUKE CHEM)

4.6: Appendix A - Use of the Agilent 100 Series UV-Vis Spectrophotometer is shared under a not declared license and was authored, remixed,
and/or curated by LibreTexts.

4.6.5 https://chem.libretexts.org/@go/page/371578
4.7: Appendix B - The Hückel Approximation
Introduction:
The molecular orbital (MO) energy level diagrams of π-electron molecules, especially where there are conjugated double bonds
(that is, an alternation of single and double bonds along a chain of carbon atoms) can be constructed using a set of approximations
suggested by Erich Hückel. The HMO method is not traditionally approached as an ab-initio quantum calculation, but is usually
used as an empirical method for doing quantum calculations. This means that some (hopefully small) number of adjustable
parameters appear in the theory, and these are determined by comparison with experiment (such as the electronic spectra of dye
molecules). If the theory is successful, it can be used to predict the spectra and properties of other conjugated systems once the
values of the adjustable parameters have been set for some set of reference compounds.
The Hückel model assumes that a π-molecular orbital, Ψ , delocalized over n atoms can be written as a linear combination of n
atomic orbitals (LCAO) ϕ i

Ψ = ∑ ci ϕi (4.7.1)

where the ϕ 's are the atomic pz orbitals on the atoms i (z is perpendicular to the molecular plane). The corresponding one-electron
i

π density at atom i is given by

2
ρiπ = c (4.7.2)
i

The energy E and the coefficients ci for the ground state are determined by use of the variation principle, which says that one
should choose the coefficients such that

∫ ΨH Ψdτ
E = = minimum (4.7.3)
2
Ψ dτ

Here the Hamiltonian, H, is an effective one-electron energy operator whose explicit form need not be specified in the empirical
HMO approach. The variation principle ensures that the lowest-energy state is as close as possible to the true energy, and the
coefficients are obtained from the minimization relations
∂E
= 0; i =1 ... n (4.7.4)
∂ci

This results in n equations of the form


c1 (H11 − S11 E) + c2 (H12 − S12 E) + . . . cn (H1n − S1n E) = 0

. .

. . (4.7.5)

. .

c1 (Hn1 − Sn1 E) + c2 (Hn2 − Sn2 E) + . . . cn (Hnn − Snn E) = 0

In these equations, H = ∫ ϕ H ϕ dτ ≡ α , called the Coulomb integral, is the energy of an electron in a 2p orbital on atom i,
ii i i i

while H = ∫ ϕ H ϕ dτ ≡ β , the resonance or bond integral, represents the interaction energy of two atomic orbitals on atoms i
ij i j ij

and j. Both α and β are negative energy quantities. S = ∫ ϕ ϕ dτ is the overlap integral, which, in the simplest approximation,
i ij ij i j

is taken to be 1 if i = j, and 0 otherwise. For a π system involving only carbon atoms, α ≡ α and Equations (Appendix B.5) take
i

the form
c1 (α − E) + c2 β12 + . . . cn β1n = 0

. .

. . (4.7.6)

. .

c1 βn1 + c2 βn2 + . . . cn (α − E) = 0

These equations have a nontrivial solution only if the corresponding secular determinant vanishes:

4.7.1 https://chem.libretexts.org/@go/page/371579
∣ α −E β12 . . . β1n ∣
∣ ∣
. .
∣ ∣
∣ . . ∣ =0 (4.7.7)
∣ ∣
. .
∣ ∣
∣ βn1 βn2 . . . α −E ∣

The HMO method makes several simplifying assumptions. The first of these are as follows:
1. The π-orbitals are considered separately from the σ-orbitals. The σ-orbitals only enter indirectly in so far as they determine the
geometry of the molecule.
2. Treat all the carbon atoms as identical. This means that all the Coulomb integrals (α ) are equal.
The next stage is to express the π-orbitals as LCAOs of the C2p-orbitals, ϕ . In ethene (CH2=CH2) we would write
Ψ = cA (ϕA ) + cB (ϕB ) (4.7.8)

and in butadiene (CH2=CH–CH=CH2)


Ψ = cA (ϕA ) + cB (ϕB ) + cC (ϕC ) + cD (ϕD ) (4.7.9)

The optimum coefficients and energies are found by the variation principle. This involves solving the secular determinant, equation
(Appendix B.7), which for ethene would be
∣ α −E β − ES ∣
∣ ∣ =0 (4.7.10)
∣ β − ES α −E ∣

and its roots are E = α ± β . The state with energy α + β corresponds to the bonding combination (β is negative) and α −β

corresponds to the antibonding combination, Figure Appendix B.1.

Figure 4.7.1 : The Hückel Molecular Orbital Energy Levels of Ethene. Two electrons occupy the lower π -orbital. (CC-NC-
BY-DUKE CHEM)
Each carbon atom supplies one electron to the π-system and the bonding orbital is occupied by an electron pair. The π-electronic
energy of ethene is therefore 2α + 2β . The excited state of the molecule, when an electron is excited into the π -orbital, lies about

2β above the ground state.

In the case of butadiene, the approximations result in the determinant


∣ α −E β 0 0 ∣
∣ ∣
∣ β α −E β 0 ∣
=0 (4.7.11)
∣ ∣
0 β α −E β
∣ ∣
∣ 0 0 β α −E ∣

This expands into the quadratic equation


(α − E)
4 2
x − 3x + 1 = 0; x = (4.7.12)
β

with roots x2 = 2.62, 0.38. Therefore, the energies of the four LCAO-MOs are
E = α ± 1.62β and E = α ± 0.62β (4.7.13)

4.7.2 https://chem.libretexts.org/@go/page/371579
as shown in Figure Appendix B.2. There are four electrons to accommodate, and so the ground state configuration is 1π 2

2
.

Figure 4.7.2 : The Hückel molecular orbital energy levels of butadiene and the top view of the corresponding π-orbitals. The four p-
electrons (one supplied by each carbon atom) occupy the two lower ( \pi \)-orbitals. Note that the orbitals are delocalized. (CC-NC-
BY-DUKE CHEM)
An important point emerges when we calculate the total \9 \pi \)-electron binding energy in butadiene and compare it with what we
find in ethene. In ethene the total energy is 2 (α + β) = (2α + 2β) ; in butadiene it is
2 (α + 1.62β) + 2 (α + 0.62β) = (4α + 4.48β) . Therefore, the energy of the butadiene molecule lies lower by
–1
(4α + 4.48β) − 2 (2α + 2β) = 0.48β (roughly –36kJ mol ) than the sum of two π-bonds. This extra stabilization of a conjugated

system is called the delocalization energy.


To obtain the wavefunction for a given energy state, Ej, one substitutes E = Ej into Equation (Appendix B.6) and solves for the jth
set of coefficients cj1, cj2, ..., cjn. This process actually gives only ratios of coefficients; to determine numerical values, we add the
normalization condition

2 2
∫ Ψ dτ = ∑ c =1 (4.7.14)
j ij

The four Hückel molecular orbitals for 1,3-butadiene are

Ψ1 = 0.372 ϕ1 + 0.602 ϕ2 + 0.602 ϕ3 + 0.372 ϕ4 (4.7.15)

Ψ2 = 0.602 ϕ1 + 0.372 ϕ2 − 0.372 ϕ3 − 0.602 ϕ4

Ψ3 = 0.602 ϕ1 − 0.372 ϕ2 − 0.372 ϕ3 + 0.602 ϕ4

Ψ4 = 0.372 ϕ1 − 0.602 ϕ2 + 0.602 ϕ3 − 0.372 ϕ4

In addition to π-electron energies and wavefunctions, the HMO approach also allows one to calculate several other quantities of
chemical interest, such as electron densities, bond orders and free valence. We will not be concerned with those here but you can
read about them in reference [13].
Heteroatoms
Heteroatoms (X) may be incorporated in the HMO method by appropriate changes in the empirical α and β parameters associated
with each atom and bond. A common example would be the inclusion of the π-electron contribution from the nitrogen atom in
pyridine. These changes are incorporated in units of the standard α and β , usually of benzene, by use of the definitions
0 0

αX = α0 + hX β0 (4.7.16)

βC X = kC X β0 (4.7.17)

Calculations involving such changes can be made by solving the secular equations. Unfortunately there remains considerable
controversy concerning the specific values to be used for the h and k-parameters of various heteroatoms. Part of the problem is that
there are several variations of MO theory and each will yield somewhat different parameter values with a given application.
We should distinguish, for example, the nitrogen in pyridine from a pyridinium nitrogen. In the former case, the heteroatom
contributes two electrons to the π-system; in the latter case it contributes one. In the second part of the analysis of your
experimental data, you will use a "canned" HMO computer program to determine the energy levels and molecular orbitals of the π-
states of the dye molecules (as a function of β). By comparison of these results with your experimental data, you will determine an

4.7.3 https://chem.libretexts.org/@go/page/371579
empirical value for β. Note that this program incorporates values for the heteroatom parameters h and k. You will need only to
specify the location and type (e.g., pyridine-type nitrogen) of each heteroatom in the molecule to be studied.

4.7: Appendix B - The Hückel Approximation is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

4.7.4 https://chem.libretexts.org/@go/page/371579
CHAPTER OVERVIEW

5: Molecular Spectroscopy of Iodine


Molecular Spectroscopy of Iodine
This is a two-part experiment consisting of the following parts:
Part I. Absorption Spectrum of Iodine Vapor
Part II. Emission Spectrum of Iodine Vapor
5.1: Prelaboratory Assignment
5.2: Introduction
5.3: Potential Energy Curves
5.4: Part I - Absorption Spectrum of Iodine Vapor - Experimental
5.5: Part I - Absorption Spectrum of Iodine Vapor - Data Analysis
5.6: Part I - Absorption Spectrum of Iodine Vapor - Questions
5.7: Part II - Emission Spectrum of Iodine Vapor - Apparatus
5.8: Part II - Emission Spectrum of Iodine Vapor - Experimental
5.9: Part II - Emission Spectrum of Iodine Vapor - Data Analysis
5.10: Part II - Emission Spectrum of Iodine Vapor - Questions
5.11: Part II - Emission Spectrum of Iodine Vapor - References

5: Molecular Spectroscopy of Iodine is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

1
5.1: Prelaboratory Assignment
 Note

During the lab period there will be a data collection period of about an hour to an hour and a half where you will be permitted
to leave lab and do independent work - please plan accordingly.

As usual, submit a short description of the background theory of the experiment, and the short experimental plan as a single
document (pdf). In addition, please do the following:
1. Read the following:
McQuarrie and Simon, Chapter 13, sections 13.1, 13.6–13.7; Chapter 15, sections 15.1–15.2, 15.5–15.6
Alberty, R.A. and Silbey, R.J. "Physical Chemistry," 2nd ed. (and 4th ed.), Wiley, NY, 1996, Chapter 14, Sections 14.1–14.3,
14.8–14.10
2. Read the entire contents of the manual for this experiment, including read Appendix A.
3. Review the Matlab Introduction activity on the Emission Spectrum of Iodine Vapor that you did on the first day of class; bring
the results of that activity with you so you can use it to analyze your data.

5.1: Prelaboratory Assignment is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

5.1.1 https://chem.libretexts.org/@go/page/373229
5.2: Introduction
What does iodine look like? Visual inspection leads to one answer: it's purple. But what does it look like at the atomic level? We
know from the Lewis dot structure of I2 that it consists of two iodine atoms with a single bond between them. But the simple rigid
dumbbell model that comes to mind does not describe the dynamics of the iodine molecule. The atoms can vibrate with respect to
each other. Furthermore, the electrons can change their energy state, thus changing the bonding in the molecule.
By using spectroscopic techniques more sensitive and quantitative than simply looking at the vapor, a more detailed description of
iodine can be obtained. The vibronic structure in a UV-Visible absorption spectrum reveals the vibrational energy levels in the
excited electronic states of iodine. From these levels a detailed description of the inter-atomic potential between the iodine atoms in
these excited states can be obtained. And, the vibronic structure in the emission spectrum from the excited states of iodine reveals
the ground state vibrational levels, from which the ground state potential energy curve can be constructed. Finally, comparison of
the ground and excited state potentials determined in this way reveals how the bond length changes and how the vibrational
potential softens (or stiffens) on electronic excitation.

5.2: Introduction is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

5.2.1 https://chem.libretexts.org/@go/page/373230
5.3: Potential Energy Curves
Harmonic Oscillator
Within the Born-Oppenheimer approximation, each electronic state in a diatomic molecule has a characteristic nuclear potential
surface V(R) describing the interatomic potential energy as a function of interatomic distance R. For a bound diatomic molecule
with low vibrational energy this interatomic potential V(R) is similar to that of a harmonic oscillator
1 2
V (R) = k (R − Re ) (5.3.1)
2

where R-Re is the deviation of the interatomic separation from the equilibrium value Re, and k is a force constant describing the
strength of the "harmonic spring" that binds the diatomic molecule together.

Figure 5.3.1 : The potential energy V(R) and vibrational energy levels for a harmonic oscillator. (CC-NC-BY_DUKE CHEM)
The potential energy curve V(R) for a harmonic oscillator is a parabola, and the energy levels are equally spaced, as shown in
Figure 5.3.1. The vibrational energy E increases linearly with the vibrational quantum number υ according to the equation
υ

1
Eυ = hν (υ + ) (5.3.2)
2

where h is Planck's constant and ν is the vibrational frequency. The vibrational frequency ν is related to the force constant k and
the reduced mass μ through the relation
−−−−−
1 k
ν = √( ) (5.3.3)
2π μ

m1 m2
where μ =
m1 + m2
for a diatomic molecule where m and m are the atomic masses of the atoms.
1 2

However, the chemical bonds in real diatomic molecules are not harmonic oscillators.

Anharmonic Oscillator
As a chemical bond stretches it gets weaker, ultimately resulting in dissociation as R approaches ∞. Also, as the bond shortens,
interatomic repulsions raise the potential more rapidly than in the harmonic potential. Thus, the potential V(R) for a real diatomic
molecule is anharmonic, and this results in a decrease in the vibrational energy spacings with increasing vibrational energy.

5.3.1 https://chem.libretexts.org/@go/page/373234
Figure 5.3.2 : A typical Morse potential for a diatomic molecule. (CC-NC-BY-DUKE CHEM)
In fact, the difference in energy between neighboring vibrational energy levels approaches zero as the vibrational energy
approaches the dissociation energy. Figure 5.3.2 shows schematically a typical anharmonic potential energy surface along with its
vibrational energy levels.

Vibrational energies of an anharmonic oscillator


The vibrational energy levels of an anharmonic oscillator can be described by a power series expansion in υ + 1

2
1 1
Gυ = ωe (υ + ) − ωe χe (υ + ) + ... (5.3.4)
2 2

where G is the energy of quantum level υ in cm–1, ω is the "harmonic frequency" of the oscillator in cm–1, and χ is called the
υ e e

anharmonicity parameter. For the iodine states of interest in this experiment we can truncate the expansion at the second term, as
indicated in equation 5.3.4. The difference between successive vibrational energies is then given by

ΔG(υ) = ∣
∣G(υ+1) − G(υ) ∣
∣ = ωe − 2 ωe χe (υ + 1) (5.3.5)

By measuring G for a given state, and plotting the change in vibrational energy spacing ΔG with respect to υ + 1 , both ω and
υ (υ) e

ω χ can be found. According to equation 5.3.5, this plot should be a straight line with a slope of −2ω χ and an intercept equal
e e e e

to ω . (Such a plot is called a Birge-Sponer plot. Hertha Sponer was a well-known and highly respected professor of physics at
e

Duke University, one of only a few women physicists at major research universities during the first half of the twentieth century.)

The Morse Potential


The Morse potential is a simple empirical potential with only two adjustable parameters that incorporates all of these features. In
particular it gives rise to the two-term expression for G in equation 5.3.4, with the vibrational energies converging on the
υ

dissociation energy De (See Figure 5.3.2). The Morse Potential has the functional form
2
−a(R−Re )
V (R) = De [1 − e ] (5.3.6)

where the Dissociation energy De and Morse parameter, a, are related to the parameters ω and ω e e χe by
−−−−−−−− −
2
8 π cμωe χe
a = √ (5.3.7)
h

ωe
De =
4χe

Thus, for this potential, the constants, a, and De can be determined from a measurement of ω and ω χ . While other empirical
e e e

potential functions provide a better representation of the true shape of the potential for many electronic states of diatomic
molecules, for the electronic states of iodine considered here the Morse potential works quite well.

5.3.2 https://chem.libretexts.org/@go/page/373234
Ground (X) and excited (B) states of iodine
In general, different electronic states have different potentials and different average interatomic separations. For some excited
electronic states, the potential may not even have a minimum, or "bound" region. Excitation to one of these "repulsive" potentials
would result directly in the dissociation of the molecule, and consequently the absorption spectrum would have no structure. In this
experiment, you will focus on electronic transitions between the vibrational levels of two bound electronic states, the ground X Σ 1
g

state and the excited B Π state. Other electronic states of the iodine molecule will not be discussed here. The potential surfaces
3
u

for the X and B states can be described by Morse functions. As the molecule dissociates along the ground state X surface, it will
separate into two iodine atoms (I + I) in their ground P electronic states. However, in the electronically excited B state,
2
3

dissociation leads to one iodine atom in its ground electronic 2


P 3 state, and the other in an excited 2
P 1 electronic state (I + I*)
2 2

(see Figure 5.3.3 below). The energy for the electronic excitation of an iodine atom E(I*) is known quite accurately from atomic
spectroscopy, the value being 7603 cm–1. This energy is just the separation in energy between the iodine molecule X and B state
potential curves in the limit where R approaches ∞ (See Figure 5.3.3).
Electronic transitions caused by absorption or emission of light will occur between the various ground state vibrational levels (υ ) ′′

and the excited state vibrational levels (υ ). Such transitions involving changes in both vibrational and electronic states are often

called vibronic transitions. By convention the various parameters associated with the excited state are designated by and the ′

parameters for the ground state are designated by , as shown in Figure 5.3.3. The frequency of a vibronic transition (in cm–1)
′′

between a ground vibrational state level (υ ) and an excited state vibrational level (υ ) will be given by
′′ ′

ν(υ′′ ,υ′ ) = νel + G(υ′ ) − G(υ′′ ) (5.3.8)

where ν is the difference in electronic energies of the two potentials at their respective minima, R and R (see Figure 5.3.3).
el

e
′′
e

Note that the energy ν , the dissociation energies D and R , and the atomic excitation energy of the iodine atom E(I*) are related
el

e
′′
e

through the equation


′′ ′ ∗
De = νel + De − E (I ) (5.3.9)

This can be seen in Figure 5.3.3.

Figure 5.3.3 : Schematic representation of the ground (X) and excited (B) states of iodine. Various electronic transitions are
indicated by arrows. (CC-NC-BY-DUKE CHEM)

5.3.3 https://chem.libretexts.org/@go/page/373234
In an absorption experiment, one usually identifies a series of bands (a "progression") that originate from a single ground state
vibrational level υ , say υ = 0 , and end in a series of different excited state vibrational levels. These bands are then assigned to
′′ ′′

the appropriate υ quantum numbers (not a trivial task!). If one then calculates the spacing between these bands, it can be seen from

equation 5.3.8 that this will yield a linear equation just like equation 5.3.5, but where the independent variable is (υ + 1) . From ′

the slope and intercept of this plot the values of ω and χ can be determined, and from these the values of the Morse parameters

e

e

a and D for the excited state can be calculated using equations 5.3.7.
′ ′
e

In an emission experiment, one usually pumps a single excited state υ , and then observes a band progression of transitions back to

a number of different ground states υ . These bands can then be assigned, and the separation between them plotted in the form of
′′

equation 5.3.5, with (υ + 1) as the independent variable. This allows the ground state parameters ω , χ , a , and D to be
′ ′′
e
′′
e
′′ ′′
e

determined.
Consider the molecular potential-energy curves shown in Figure 5.3.3. They relate to different electronic states of the same
molecule. Some of the vibrational energy levels of each electronic state are marked, together with some of the vibrational
wavefunctions. Before absorption of light occurs, the molecule is most likely to be in the ground vibrational state of the lower
electronic state, υ = 0 , and so we shall concentrate on that initial level. The form of the lowest vibrational wavefunction shows
′′

that the most probable location of the nuclei is near their equilibrium separation R . However, at room temperature the υ = 1 and
′′
e
′′

υ = 2 levels of iodine also have significant thermal populations and so electronic transitions originating from the υ = 1 and
′′ ′′

υ = 2 levels can also be seen; these are known as hot bands.


′′

When an electronic transition occurs, the molecule is excited to the state represented by the upper curve in Figure 5.3.3. According
to the Franck-Condon principle, the vibronic transitions occur so fast (on the order of femtoseconds or less) that the nuclear
framework remains constant during the actual transition. We may, therefore, imagine the nuclear wavefunction of the molecule as
rising up the vertical lines shown in Figure 5.3.3. This is the origin of the term vertical transition, which is used to denote an
electronic transition that occurs without change of nuclear geometry.
The vertical transition cuts through several vibrational levels of the upper electronic state. The level marked υ = 25 is the one in

which the nuclei are most probably at the separation R because the wavefunction has maximum amplitude there. Therefore this is
′′
e

the most probable level for the termination of the excitation. It is not, however the only vibrational level at which the transition may
end, because several of its neighbors also have a high probability of having nuclei at the separation R . Therefore transitions take
′′
e

place to all vibrational levels in this region, but with greatest probability to the level with the wavefunction that overlaps the initial
vibrational state wavefunction most favorably. Transitions from the υ = 0 state of iodine would tend to "miss" upper-level
′′

vibrational states with low values of υ . For this reason, these transitions are weak or not observed.

As stated above, the excited state vibrational wavefunctions which have the greatest Franck-Condon overlap with the υ = 0 ′′

ground state level, are those with the greatest amplitude near R . For these levels the maximum amplitude in the vibrational
′′
e

wavefunction occurs at the inner "classical turning point" or the edge of the potential well. Such wavefunctions lie at high values of
υ as shown in Figure 5.3.3. Thus we see that the vibrational structure of the spectrum depends on the relative displacement of the

two molecular potential energy curves, and a long progression of vibrations (a lot of vibrational structure in the spectrum) is
stimulated if the two electronic states are significantly displaced.
This Franck-Condon principle can be used to estimate the shift in the potential minimum on excitation, i.e. R − R . The strongest
′′
e

e

absorption transition (from υ = 0 ) will be to some particular excited state level υ . The separation, R , at the inner turning point
′′ ∗′ ∗′

of this exited state level can be identified. Assuming the shapes of both the ground and excited state potential curves are known, the
two curves can be overlaid so that the minimum in the ground state R overlaps vertically with the excited state turning point value
′′
e

R . The displacement, R − R , can then be read from the plots, or calculated analytically (see references 1 and 4).
∗′ ′ ′′
e e

5.3: Potential Energy Curves is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

5.3.4 https://chem.libretexts.org/@go/page/373234
5.4: Part I - Absorption Spectrum of Iodine Vapor - Experimental

The ELN template for this entire experiment is available here: ELN01_1_dyes2022.docx

Part I - Collect the Absorption Spectrum of Iodine Vapor


Here is the Basic procedure for using the Agilent 100 Series UV-Visible spectrophotometer in the Molecular Spectroscopy of
Iodine experiment.
There are two sample holders in the instrument. The Sample is placed in the front holder; and, the Reference is placed in the back
holder.
With both sample and reference holders empty and the cover closed, turn on the instrument with the switch at the front, lower left.
Wait two minutes before proceeding further.

Validate the Wavelength Accuracy of the UV-Vis Spectrometer


The first thing to do when you use an instrument is to be sure it is working well. In this case, we want to know that our wavelength
measurements are accurate. Follow the steps below to verify that the UV-vis spectrometer is measuring wavelength as accurately as
we expect.

1. Load the Validation Software by clicking the icon from the Windows main screen. After a few moments the software
screen will appear and will begin to communicate with the instrument. You can follow the progress by monitoring the display at
the lower left of the screen. When the traffic light in the upper center turns green, the system is ready. You will see absorbance
in the upper left and wavelength in the upper right.

Figure 5.4.1 : The Validate Main Screen. (CC-NC-BY-DUKE CHEM)

2. Click the button and deselect, using the button, all tests except Wavelength Accuracy - Ho2O3. The screen
should look as shown here.

5.4.1 https://chem.libretexts.org/@go/page/373393
Figure 5.4.2 : Selecting the Holmium Oxide Validation. (CC-NC-BY-DUKE CHEM)

3. Click the button and, when the setup screen appears, place your names in the Name line; and, the name of the test in
the comments box. In the Options section deselect Auto Print, Results and Graph. Data Form and Company Logo should not be
checked. % Page Height set at 50 is fine.
4. Click to return to the main screen. Place the holmium oxide filter in the sample holder of the instrument, leave the

reference holder empty, and close the cover. Click to begin the holmium oxide wavelength accuracy test. As the test
proceeds the spectrum of holmium oxide will be displayed. At the conclusion of the test, a report will be generated on the
screen. Include these data in your laboratory report. Check to be sure that the difference between the expected and observed
peak positions is one the order of ±0.2-0.4 nm. If it is greater than that, consult your TA or your Instructor before proceeding
further.

Figure 5.4.3 : Results of the Holmium Oxide Wavelength Accuracy Test. (CC-NC-BY-DUKE CHEM)
Exit the Validation software and return to the Windows Desktop.

Collect the Absorption Spectrum of Iodine Vapor

1. Load the Scan Software by clicking the icon on the Windows main screen. In a few moments the scan software main screen
will appear and will begin to communicate with the instrument. You can follow the progress by monitoring the display at the
lower left of the screen. When the traffic light at the top center of the screen turns green, the system is ready. You will see
absorbance at the upper left and wavelength at the upper right of the screen.

5.4.2 https://chem.libretexts.org/@go/page/373393
Figure 5.4.4 : The Scan Module Main Screen. (CC-NC-BY-DUKE CHEM)

2. Click the button and wait for the setup screen to appear. Select, first, the Cary tab to begin setting parameters for
your sample runs. Under Xmode, the Mode should be set to nanometers and your scan range from 620 to 490. Under Ymode,
the Mode should be set to absorbance and the range from 0.02 to 0.25. The Ave time should be 0.1 seconds and the Data
interval should be set to 0.1 nm. This will yield a Scan Rate of 60 nm/min.

Figure 5.4.5 : The Scan Setup Screen. (CC-NC-BY-DUKE CHEM)


3. Select the Options tab on the setup screen and set the Spectral Band Width (SBW) to 0.2 nm. Select Double Beam Mode. Select
UV-Vis as the source and be sure the Source changeover is set to 350 nm. Under Display options you can choose Individual
data.

5.4.3 https://chem.libretexts.org/@go/page/373393
Figure 5.4.6 : The Options Tab of the Scan Setup Screen. (CC-NC-BY-DUKE CHEM)
4. Select the Reports tab. Put your names in the Name box and the title of the experiment in the Comment box. Under Options,
check Data Form, Graph, Focused trace. and % Page Height 50. Be sure the Include XY pairs table is not checked. Click Peak
Information. Under Peak Table, put a check in All peaks, set the Threshold to 0.0200, and put Peaks in the Peak Type box.
Click to return to the Setup screen. Choose Select for ASCII (csv). Click to return to the main screen.

Figure 5.4.7 : The Reports Tab of the Scan Setup Screen. (CC-NC-BY-DUKE CHEM)
5. Remove the 1 cm cuvette holder from the Sample (front) beam of the instrument and replace it with the 10 cm sample cell

holder. Place your iodine vapor sample in the Sample holder and close the cover. Click to begin your run. You will be
prompted for a filename. After naming the file you will be prompted for a sample name. Place the name of the sample (Iodine)
in this box. As the run proceeds you will see the spectrum generated in real time on the screen.

5.4.4 https://chem.libretexts.org/@go/page/373393
Figure 5.4.8 : The Absorption Spectrum of Iodine Vapor. (CC-NC-BY-DUKE CHEM)
6. At the conclusion of the run, find the υ = 0 to υ = 25 transition and label it on your spectrum (by using the available
′′ ′

annotator - ask your TA). Also locate the wavelength showing the largest absorption and label it on your spectrum.

 Note

Since the baseline is not straight, the largest peak is found by consideration of both the base and the peak height.

7. Remove the sample and sample holder and replace the 1 cm cuvet sample holder in the instrument, using the screw to secure it
properly. Be sure you have screenshots of the needed results, and the Excel spreadsheet generated at the end of the run. Exit the
software and turn off the instrument.

5.4: Part I - Absorption Spectrum of Iodine Vapor - Experimental is shared under a not declared license and was authored, remixed, and/or curated
by LibreTexts.

5.4.5 https://chem.libretexts.org/@go/page/373393
5.5: Part I - Absorption Spectrum of Iodine Vapor - Data Analysis
1. Given that the υ = 0 to υ = 25 transition occurs at 545.5 nm, assign the peaks in the spectrum according to their υ values.
′′ ′ ′

Do not forget to make corrections to the wavelength based on the results of the holmium oxide calibration. Be careful not to
confuse peaks originating from υ = 0 with those of hot bands originating from υ > 0 . Choose 20 peaks, and convert their
′′ ′′

wavelengths to wavenumbers.
2. Using MATLAB plot ΔG (wavenumbers) vs (υ + 1) (see Equation 5.3.5). Perform a linear least squares regression and

(υ)

overlay the fitted line on the plot of the data.


3. From the intercept calculate ω using equation 5.3.5.

e

4. From the slope calculate χ using equation 5.3.5.



e

5. From equation 5.3.7 calculate D . ′


e

6. Estimate your error in measuring the peak positions. Using the standard deviation of the slope and intercept, and equation 5.3.5,
calculate the uncertainties in ω , χ and D .

e

e

e

7. Calculate the Morse parameter, a from equation 5.3.7.


8. Compare your results with literature values (Ref. 8, p. 540, and Ref. 4).

5.5: Part I - Absorption Spectrum of Iodine Vapor - Data Analysis is shared under a not declared license and was authored, remixed, and/or
curated by LibreTexts.

5.5.1 https://chem.libretexts.org/@go/page/373397
5.6: Part I - Absorption Spectrum of Iodine Vapor - Questions
1. Determine the value of ν 00 (the energy of the transition from υ ′′
=0 to υ ′
=0 ) by appropriate substitution into Equation 5.3.4.

 Note
To do this you will need to have already done the data analysis for Part II and will also need Equation 5.3.8.

2. Why do the vibrational spacings get narrower as V(R) increases?


3. Why is the Morse potential asymmetric?
4. Why are some of the transitions called "Hot Bands"?
5. How would you calculate χ and D using the hot bands?. (See reference 6 for help.)
′′
e
′′
e

6. What are the units of the Morse parameter, a, in Equation 5.3.7?

5.6: Part I - Absorption Spectrum of Iodine Vapor - Questions is shared under a not declared license and was authored, remixed, and/or curated by
LibreTexts.

5.6.1 https://chem.libretexts.org/@go/page/373399
5.7: Part II - Emission Spectrum of Iodine Vapor - Apparatus
In this part of the experiment you will use a HeNe laser at 543.5 nm to excite iodine to the B state, and a monochromator and
photomultiplier to measure the spectrum of the emission back to the ground X state. It is worth noting that this technique of "laser-
induced fluorescence" is an important spectroscopic tool, not only for physical chemists, but also for analytical chemists. The
sensitivity is extremely high, and this has led to a wide variety of applications. These include the measurement of trace amounts of
osmium and iridium in rock samples as part of an effort to test a theory for the sudden destruction of the dinosaurs by a catastrophic
collision of the earth with an asteroid!
Since the laser used in this experiment can be harmful, and the apparatus is relatively expensive, a number of safety precautions
must be taken.

 Safety and Other Precautions


1. Don't look into the laser! This is a class IIIa/3R laser. If the beam is directed into your eye it will damage your eye before
you can blink. However you will still want to avoid looking directly into the beam, or having the beam be directed at you.
Any glass that the laser hits will reflect about 5% of the beam intensity. This specular reflection may accidentally be
directed toward your eyes. Therefore, do not move items into or out of the laser path with the laser on. Wear laser
goggles whenever the laser is on.
2. Don't vent the cell. There is no danger to you if you vent it, but you won't be able to do the experiment until the cell is
evacuated again.
3. Don't kill the PMT. The Photomultiplier Tube (PMT) is an extremely sensitive light detector. If you can see the light, that's
enough light to overload the PMT if there is voltage on it. Read and understand the experimental procedure before applying
any voltage to the PMT. Do not turn on the voltage to the PMT if it is exposed to significant amounts of light. Also, turn up
the voltage slowly as you look for signal. As long as the observed signal stays below 0.2 V the PMT will not be damaged.
Again, no physical damage to you will result from overloading the PMT, only financial.

Figure 5.7.1 below shows the setup used in this experiment.

Figure 5.7.1 : Instrumentation used for Laser Fluorescence Experiments. (CC-NC-BY-DUKE CHEM)
Laser: A Research Electro-Optics, Inc. HeNe laser (model 30968), which emits visible light (543.5 nm), is used to pump primarily
transitions from υ = 0 .
′′

Monochromator: A Triax 320 monochromator is used to disperse the emission signal from the sample. The entrance and exit slits
of the monochromator control the signal intensity and resolution. It is best to use very narrow slits when looking for signal of
unknown intensity so you don't overload the photomultiplier tube (PMT) detector with too much light. The PMT is mounted on the
exit slit, which should be set to the same width as the entrance slit. There is a slide shutter on the exit port. The manual control unit
(CD-2a) can be used for tuning the monochromator when you are setting up the experiment. The monochromator is computer
controlled during data collection.
PMT Power Supply: The PMT power supply allows you to introduce a voltage to the PMT. Typically you will use a voltage off
1100 volts for this specific PMT. Depending on the available time, your TA may show you the effect of PMT voltage on the

5.7.1 https://chem.libretexts.org/@go/page/373401
observed signal.
Lock-In Amplifier: In association with Stanford Research Systems (SRS) chopper an SRS510 Lock-in Amplifier is used to
improve the signal to noise ratio of the emission spectrum. The interested student should read the Bentham Instruments
introduction to the principles of a lock-in amplifier in the Course Documents section of your Sakai site, or at
http://www.bentham.co.uk/pdf/F225.pdf.
Spectra Acq2: This is an analog to digital (A/D) converter which takes the analog signal (current from the PMT) and converts it to
a digital signal that the SynerJY computer software can read and process.

5.7: Part II - Emission Spectrum of Iodine Vapor - Apparatus is shared under a not declared license and was authored, remixed, and/or curated by
LibreTexts.

5.7.2 https://chem.libretexts.org/@go/page/373401
5.8: Part II - Emission Spectrum of Iodine Vapor - Experimental
 Note

The Helium-Neon laser frequency of 543.5 nm corresponds to one of the larger absorption peaks from the absorption spectrum
obtained in Part I of the experiment. Notice that absorptions close to 543 nm in the absorption spectrum (Part I) are among the
most intense absorptions.

Procedure
1. Use the high vacuum line, with help from your TA, to evacuate the iodine vapor cell that will be used in the emission
experiment.

Figure 5.8.1 : A sample cell containing solid iodine is placed in a liquid nitrogen dewar and evacuated using the main manifold
of a high vacuum line. (CC-BY-SA; Kathryn Haas, Duke Chem)
2. Once the cell is well evacuated, place it on the sample holder on the laser table. The cell should be in the path of the laser and
there should be nothing obstructing the light path.

Figure 5.8.2 : A sample cell is positioned on the sample cell holder in the path of the laser. The valve on the sample chamber is
positioned toward the laser and away from the detector, as shown. (CC-BY-SA; Kathryn Haas, Duke Chem)
3. Turn on the Chopper using the power switch on the back of the power supply. It should already be set to approximately 500 Hz.
If necessary, adjust settings using the frequency adjust controls to bring it to approximately 500 Hz.

5.8.1 https://chem.libretexts.org/@go/page/373403
Figure 5.8.3 : The chopper and its power supply are labeled. Turn on the chopper using the switch on the back side of the power
supply. (CC-BY-SA; Kathryn Haas, Duke Chem)
4. Turn on the power for Amplifier and Photomultiplier Tube (PMT).
a. Turn on the Pre-Amp power supply using the toggle switch on the box.
b. Turn on the Lock-in Amplifier using the black switch on the front right of the box.
c. Turn on the PMT Power Supply, but do NOT apply voltage yet.

Figure 5.8.4 : a) The Pre-Amp power supply is the black box in the center of the picture; it is turned off using the silver
toggle switch. b) the Lock-in Amplifier is the box on the left of the photo; it is powered on using the black toggle switch at
the bottom right of the control panel. c) The PMT power supply is the box on the right; it is powered on using the press
button at the bottom right. Do NOT touch the black switch to apply voltage to the PMT until the lights in the room are off.
(CC-BY-SA; Kathryn Haas, Duke Chem)
5. Turn off the lights in the room.
6. With the lights off, turn on the laser (key on the laser power supply). Notice that while the laser beam itself is green (543.5 nm),
the emitted light, seen as a beam along the length of the sample cell, is in the yellow orange, that is, at lower energy (higher
wavelength) than the exiting laser. Of course what you see in the cell is the sum of all wavelengths of emitted light. To resolve
that into a spectrum, you will send that light into the monochromator and scan an appropriate wavelength range.

7. Click the SynerJY icon on the computer screen to load the SynerJY software.
8. From the Collect menu choose Experiment Setup.

Figure 5.8.5 : Choosing Experiment Setup from the Main SynerJY Window. (CC-NC-BY-DUKE CHEM)
9. You will briefly see two windows come and go as the software connects with the monochromator. The first to appear is a
Hardware Configuration window, and Triax320PMT should be highlighted in it.

Figure 5.8.6 : Hardware Configuration. (CC-NC-BY-DUKE CHEM)


10. This window will disappear and be replaced by an Initialization Window. The software will initialize the monochromator and
PMT and then that window will also disappear.

5.8.2 https://chem.libretexts.org/@go/page/373403
Figure 5.8.7 : The Initialization Window. (CC-NC-BY-DUKE CHEM)
11. The software will then load the Experiment Setup screen.

Figure 5.8.8 : The Setup Screen. (CC-NC-BY-DUKE CHEM)


12. The first thing we must do is optimize the emission signal. Type a name for your Experiment File (use your initials and date, or
something easily recognizable). Also give a name for the Data File. Choose a Start Wavelength of 535 nm and an End
Wavelength of 550 nm. Set the Increment to 0.5 nm. Make sure Type is set to Wavelength, and Triax320 is selected under
Scanning. Click the button at the left of the screen and set the Front_Entrance and Front_Exit Slits each to 0.04 mm, or to
a value suggested by your TA.

Click on the button at the left of the screen. Make sure there is a check in the PMT Active box and set the Integration
Time to 3000 ms (i.e. 3 seconds). Leave the Gain at x1.
13. Apply 900 volts to the PMT from the power supply. Click on the button. Make sure the Slits and Integration Time are set

properly. Change the value in the Position box to 543.5 nm. Click on the button. Wait a few seconds and look at the digital
readout on the lock-in. You should see less than a tenth of a millivolt signal at this point. Change the PMT voltage to 1000 V
and watch the signal at the lock-in. Make sure it does not go off scale; decrease the sensitivity if needed. Bring the PMT voltage
up to 1100 V, adjusting the lock-in sensitivity to obtain a good signal (several tenths of a mV) without going off scale. Click the

button to return to the main Experiment screen. Click on the button to start this short scan. You will see the spectrum
appear in the Display window as the scan progresses. You should see a peak at the pump wavelength. When the scan is
complete you may be asked about saving files. To be safe, save all files when asked. Your scan will then appear on the screen.
Select Pick Peaks from the Analysis menu. To begin the peak analysis leave the default settings in the boxes and click on Find
Peaks at the bottom of the Pick Peaks box. If you do not like the results, alter the Height and Width values until you get a good
result. In the example shown at the right the major peak is found at 543.5 nm, demonstrating that our system is working
properly.

Figure 5.8.5 : The Results of the Test Scan. (CC-NC-BY-DUKE CHEM)


14. Now you are ready to do any final optimizing of the signal (check with your TA) and start the run.
15. Change the Start wavelength to 510 nm and End wavelength to 800 nm (or some value suggested by your TA).

Figure 5.8.6 : Setting Up The Sample Run. (CC-NC-BY-DUKE CHEM)

Make sure the Increment is set to 0.2 nm. Click the button to make sure the Integration Time is still set at 3000 ms.
16. Click the Preview Button to set the monochromator to the start wavelength. Click the Run button once and wait a few seconds
(this is equivalent to setting the monochromator to the start wavelength and allowing the detector signal to relax). Then click the
Cancel Button to return to Setup screen.

17. Click the button from the setup screen to start the full scan. Watch the data collection for the first 5-10 minutes to be sure
everything is running correctly. The total scan will take a while to complete. If you leave the room, close the door behind you to
be sure that stray light does not get into the experiment.

5.8.3 https://chem.libretexts.org/@go/page/373403
18. At the end of the scan you may want to rescale the horizontal axis to see the spectrum more clearly. Repeat the Pick Peaks
procedure you did earlier on the short scan. Likely you will need to adjust the parameters in the Pick Peaks dialog box to get
most or all of the peaks.

Figure 5.8.7 : The FInal Spectrum. (CC-NC-BY-DUKE CHEM)


19. When you get a good listing of the peaks look at the bottom of the screen where you will see a list of the files associated with
these data. Click on the Peaks1 worksheet and you will get a table of the peaks found in the spectrum.

Figure 5.8.8 : The Peaks Table. (CC-NC-BY-DUKE CHEM)


You can copy the appropriate data from there to an Excel spreedsheet. You should the save the spreadsheet and any screenshots
that you need.
20. Turn off the laser. Then turn off the power supply to the PMT. Then turn off the power to the Pre-amp Power Supply and to the
Lock-in Amplifier. Turn off the Chopper. Close the SynerJY software.

5.8: Part II - Emission Spectrum of Iodine Vapor - Experimental is shared under a not declared license and was authored, remixed, and/or curated
by LibreTexts.

5.8.4 https://chem.libretexts.org/@go/page/373403
5.9: Part II - Emission Spectrum of Iodine Vapor - Data Analysis
1. Assign at least 19 peaks in the emission spectrum based on the known pump transition. Remember, this is an emission
spectrum. Longer wavelengths (lower energy) correspond to transitions to higher ground state vibrations.
a. Convert the wavelengths of these peaks to wavenumbers.
b. Create a Birge-Sponer plot (see Equation 5.3.5) by plotting ΔG vs (υ ′′
υ
′′
+ 1) .
2. From the intercept of the Birge-Sponer plot, calculate (using Equation 5.3.5).
′′
ωe

3. From the slope of the Birge-Sponer plot, calculate χ (using Equation 5.3.5).
′′
e

4. From the Morse Potential (using Equation 5.3.7), calculate D . ′′


e

5. Estimate the uncertainty in the peak positions using the standard deviations from the least squares analysis and calculate the
uncertainties in ω , χ and D .
′′
e
′′
e
′′
e

6. For comparison, calculate the value of D using Equation 5.3.9 as follows: From the wavelengths of at least three absorption
′′
e

lines, calculate the corresponding transition energies ν (where υ = 0 ) in units of cm–1. For these same transitions,
( υ′′ , υ′ )
′′

calculate (G − G

(υ ) ) using Equation 5.3.4 and your values of ω , ω , χ and χ . Now use Equation 5.3.8 to calculate the
(0 )
′′
′′
e

e
′′
e

e

average value of ν for the three transitions. Finally, use Equation 5.3.9 and your value of D from Part I to calculate D . The
el

e
′′
e

value of E(I*) is known to be 7603 cm–1 from atomic spectroscopy.1 Which value of D has the least uncertainty (neglect the ′′
e

error in E(I*))?
7. Calculate the Morse parameter (a") from Equation 5.3.7.
8. To plot the Morse potential wells, you will use the Matlab routine morse. Simply type morse <ENTER> and follow the screen
prompts. Separate ground and excited state Morse potential curves will be plotted using your experimental data, and sent to the
printer. For now, both R and R are set at zero. Vibrational quantum levels will be shown for both ground and excited states.
′′
e

e

9. The maximum intensity in the absorption spectrum occurs as the excited state vibrational wave functions overlap vertically with
the ground state wave functions. A Matlab script, overlap, is used to overlay the two plots and shift the ground and excited state
potential wells such that the center of υ = 0 overlaps with the left edge of the υ = 25 . At the Matlab prompt type overlap
′′ ′

<ENTER>. The script will animate the overlap of the ground and excited states and the vertical transition. It will also estimate
the displacement (R − R ) (Ångstroms) between the minima of these two potential wells along the R-axis. Save a copy of this

e
′′
e

graph (or take a screenshot) for your ELN, and record the displacement.

 Note
The ground and excited state potential wells are overlaid so that the displacement R − R may be measured. However, ′
e
′′
e

this gives a misleading impression of the relative positions (energies) of the two potential wells. Clearly the excited state
potential well is at much higher energy than the ground state potential well (see Figure 5.3.3 and note the discontinuity in
the V(R) axis). A realistic picture is shown in reference [10]. In the Introduction to Matlab Excercises, you plotted the
ground and excited state potential wells for iodine using literature values for De, a, and Re. You will now repeat this
assignment using your results.

10. The functional form of the Morse potential may be represented by the equation
2
−a(R−Re )
V (R) = Te + De [1 − e ] (5.9.1)

where Te is the electronic term energy (the energy at the bottom of the potential well). T and T represent the electronic term e
′′
e

energies of the ground and excited states respectively. If we set T = 0 , then T = ν . Given that R = 2.66 Å (obtained
e
′′
e

el
′′
e

from IR spectroscopy), calculate R using your value of R − R .



e

e
′′
e

11. You will now use your values of T (zero), T (= ν ), D , D , a", a', R and R to replot the Morse potential wells for the
e
′′
e

el
′′
e

e
′′
e

e

ground and excited states of iodine on the same graph. Return to your Matlab command window, and as instructed, press any
key to continue. Enter your calculated T . The two Morse potentials will be plotted on one graph and the vibrational quantum
e

levels will be shown for both ground and excited states. Save a copy of this graph (or take a screenshot) for your ELN.
12. Compare all of your results with literature values (See, for example, reference. 1, p. 430-432).

5.9: Part II - Emission Spectrum of Iodine Vapor - Data Analysis is shared under a not declared license and was authored, remixed, and/or curated
by LibreTexts.

5.9.1 https://chem.libretexts.org/@go/page/373515
5.10: Part II - Emission Spectrum of Iodine Vapor - Questions
1. Why are the vibrational spacings different in the ground and excited states?
2. Why is the shape of the Morse potential different in the ground and excited states?
3. How is it possible to observe emission peaks at higher energies than the pump energy?

5.10: Part II - Emission Spectrum of Iodine Vapor - Questions is shared under a not declared license and was authored, remixed, and/or curated by
LibreTexts.

5.10.1 https://chem.libretexts.org/@go/page/373517
5.11: Part II - Emission Spectrum of Iodine Vapor - References
1. Shoemaker, D.P., Garland, C.W., Nibler, J.W., Experiments in Physical Chemistry, 6th ed., McGraw-Hill, NY, 1996, Chapter
XIV, Experiment 40.
2. Silbey, R.J, Alberty, R.A. and. Bewendi, M.G., Physical Chemistry, 4th ed., Wiley, NY, 2005, Chapter 14, Sections 14.1–14.4,
14.8, 14.9.
3. Steinfeld, J.I., J. Chem. Educ., 1965, 42, 85.
4. McNaught, I.J., J . Chem. Educ., 1980, 57, 101.
5. Tellinghuisen, J., J. Chem. Educ., 1981, 58, 438.
6. Snadden, R.B., J. Chem. Educ., 1987, 64, 919.
7. Capelle, Gene A. and Broida, H.P., J. Chem. Phys., 1973, 58, 4212.
8. Herzberg, B., Molecular Spectra and Molecular Structure I. Spectra of Diatomic Molecules, 2nd ed., Van Nostrand, Princeton,
N.J., 1950, Chapters 2-4, 8.
9. Wright, J.C. and Wirth, M.J., Anal. Chem., 1980, 52(9), 1087A.
10. Lievin, J., J. Chem. Ed., 1982, 59(9), 777.
11. Baxter, G.P. and Grose, M.R., J. Amer. Chem. Soc., 1915, 37, 1061.
12. Atkins, P.W., Physical Chemistry, 5th ed., W.H. Freeman, NY, 1994, Sections 17.1-17.5, or, 6th ed. W.H. Freeman, NY, 1997,
Sections 17.1-17.5.

5.11: Part II - Emission Spectrum of Iodine Vapor - References is shared under a not declared license and was authored, remixed, and/or curated
by LibreTexts.

5.11.1 https://chem.libretexts.org/@go/page/373518
CHAPTER OVERVIEW

6: Rotation-Vibration Spectrum of HCl AND DCl; Vibrational Spectrum of SO2


Part I. The general approach to this experiment is adapted from Shoemaker, D.P., Garland, C.W. and Nibler, J.W., Experiments in
Physical Chemistry, 6th edition, McGraw Hill Co. Inc, NY, 1996. In preparation for this experiment please read Experiment 37 in
Shoemaker, Garland and Nibler beginning on page 397. You will want to read and understand the theory section including the
subsections involving selection rules and isotope effects. In addition to this information the following section on transition
intensities is offered.
Part II. The general approach to this experiment is adapted from D. P. Shoemaker, C. W. Garland, and J. W. Nibler, Experiments in
Physical Chemistry, 6th edition, McGraw Hill Co. Inc, NY, 1996. In preparation for this experiment please read Experiment 35 in
Shoemaker, Garland and Nibler beginning on page 383. You will want to read and understand the theory section including the
subsections on the valence-force model and the vibrational partition function.
6.1: Prelaboratory Exercise
6.2: Transition Intensities
6.3: Part I - Experimental
6.4: Part I - Data Analysis
6.5: Part I - References
6.6: Part II - Experimental
6.7: Part II - Data Analysis
6.8: Part II - References

6: Rotation-Vibration Spectrum of HCl AND DCl; Vibrational Spectrum of SO2 is shared under a not declared license and was authored,
remixed, and/or curated by LibreTexts.

1
6.1: Prelaboratory Exercise
 Important note

Students in CHEM301L will only complete Part I of this module (HCl/DCl) and its related assignments.
Students in CHEM310L will complete all of this module, including both Part I (HCl/DCl) and part II (SO ). 2

Reading assignment
This experiment involves many calculations; advanced knowledge of what you will have to do will be very helpful. In preparation
for this experiment please read Experiment 37 in Shoemaker, Garland and Nibler beginning on page 397. You will want to read and
understand the theory section including the subsections involving selection rules and isotope effects. In addition to this information
read this module (not that there are two parts, and what you do depends on your course number). Additional references below are
strongly recommended.
Part l. Instrumentation, Journal of Chemical Education, 1986, 63 (1), A5.
Shoemaker, D.P., Garland, C.W., and Nibler, J.W. Experiments in Physical Chemistry, 6th ed., McGraw-Hill, New York, 1996,
Chapter XIV, Experiments 35 and 37.
McQuarrie, D.A. and Simon, J.D. "Physical Chemistry: A Molecular Approach", University Science Books, CA, 1997, Sections
13.2-13.4, 18.4-18.5); Atkins, P., Physical Chemistry, 5th ed., sections 16.8 – 16.11; and, Silbey, R.J., Alberty, R. A.
Bawendi, M.G., Physical Chemistry, 4th ed., Sections 13.4, 13.6 and 13.7, for a discussion of rotation-vibration spectroscopy of
diatomic molecules.

Written assignment
In addition to the usual prelab work expected for all experiments in the course (Introduction and Experimental Plan), please
complete the following and submit the work prior to the lab meeting.
~ ∗ ∗
νo B
For Part 1 (CHEM310L and 301L), show which equations you would use to find: B ; α ; ν~ ; k;
e e o ~ ; and e

Be
νo

kδ ~
For Part 2 (CHEM310L only), show which equations you would use to determine: k ; 1 2
;C v (vib)

Helpful videos
These videos provide a nice introduction and demonstration of the theory and techniques you will use.

Introduction to rotovibrational spectroscopy

Rovibrational Spectroscopy

Filling a gas cell

6.1.1 https://chem.libretexts.org/@go/page/372476
HCl/DCl Lab

Pre-Lab Assignment:
You should prepare for lab meetings by reading the appropriate module in this manual, and completing the specified pre-laboratory
assignments. The pre-laboratory assignments are designed to prepare you to operate efficiently and safely during the experimental
module. You will not be allowed to conduct experiments if the pre-lab assignment is not completed before your lab meeting. If
there is an extenuating circumstance that prevents your preparation, please inform your instructors as soon as possible.
For a given experiment, you should complete the assigned reading and write out the specified pre-lab assignments. The pre-lab
assignment must be submitted as a single pdf document on Sakai at least 24 hours before the scheduled lab meeting. Late pre-
laboratory work will not be accepted and will receive a score of zero. Your TA will evaluate the work and return it no later than the
beginning of your lab section. Your preparedness will also be evaluated by your answers to oral questions pertaining to the
laboratory from the TA and other lab instructors. The pre-lab and preparation accounts for 10 points of each experiment requiring a
pre-lab.
For every experimental module there is a general pre-lab assignment. It is described here:
General Pre-lab assignment:
Write a brief introduction to the experiment and the experimental plan using the following outline:
I. Experiment Introduction: In your own words, briefly describe the underlying theory and explain the principle upon which the
experiment isbased. Keep it concise; lengthy discussions or derivations of equations are unnecessary.
a. Background (importance/relevance): A brief description and justification of the importance of the topic, with references to
source material. What is the topic and why is it interesting and important? This is meant to "hook" the reader so they want
to read about the topic.
b. Theory: A concise paragraph or more describing the general theory for the method/technique used; this section should
contain a reference to the appropriate pages of the textbook, or other source material from the literature.
c. Goals/Purpose: Brief statement of purpose, which should indicate what is analyzed and the technique used. Limit to three to
five sentences.
II. Experimental Plan: a short summary of the specific things that you would need to know in order to do the experiment without
having access to your lab manual, with the understanding that you will always have a TA available to teach you how to use the
instrument. The best way to approach this task is to thoroughly read the procedure for the experiment. While reading, take notes
on the specifics of the work. What chemicals will you use? If you have to make solutions, what solutions do you need to make
and of what concentration? How will you prepare them? What instrument will you use, and what specific instrument settings
will you need to input (sometimes you will find this information in the Appendices that discuss the use of the instruments)?
How many runs of each sample will you do; and over what wavelength, or temperature, etc. range? You do not need to write
down how to use the instrument; we will teach you that. Bottom Line: we want you to write down only that pertinent
information about samples, sample preparation, instrument and instrument parameters, numbers of runs, and other things, that
you will need to know in order to do your work.
Pre-lab assignments should be typed. If you prefer to write out equations by hand, they should be neat and inserted as an image in
the appropriate place. If you need advice for how to insert equations or images into your documents, please ask for help at least 48

6.1.2 https://chem.libretexts.org/@go/page/372476
hours before the assignment is due; an in-person (or zoom) meeting is most appropriate for this issue.

Always Wear Protective Clothing:


Please dress appropriately for a laboratory environment - this is considered an important part of preparation for the laboratory
period (see also Laboratory Safety below). You must wear close-toes shoes and clothing must cover your body from your upper
arms to your knees. You should also wear protective eyewear. If lab attire is not jiving with your wardrobe, please purchase a lab
coat and bring it with you so that you can cover up in lab. No student should be in lab without wearing protective eyewear and
appropriate protective clothing.

6.1: Prelaboratory Exercise is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.
Introductory Details by W. R. Fawcett, John Berg, P. B. Kelley, Carlito B. Lebrilla, Gang-yu Liu, Delmar Larsen, Paul Hrvatin, David
Goodin, and Brooke McMahon is licensed CC BY-NC-SA 4.0.

6.1.3 https://chem.libretexts.org/@go/page/372476
6.2: Transition Intensities
Transition Intensities
The intensity of a spectral line depends not only on the transition probability and the frequency ν but also on the number of
molecules in the initial state. Thus for a theoretical prediction of intensities, a knowledge of the numbers of molecules in the
various initial states is necessary, in addition to a knowledge of the transition probabilities. Since most infrared (IR) spectra are
observed under conditions of thermal equilibrium, we need only consider the distribution of the molecules over the different
quantum states in thermal equilibrium.
According to the Boltzmann equation, the number of molecules dNE that have a classical vibrational energy between E and
E

E + dE is proportional to e dE , where k
k
B
T
is Boltzmann's constant and T is the absolute temperature. In quantum theory,
B

only discrete values, E , are possible for the vibrational energy. The number of molecules in each of the vibrational states is then
(ν )
E
(ν)

proportional to the Boltzmann factor e . The zero-point energy can be left out, since to add this to the exponent would only
k T
b

result in a multiplication factor that would be the same for all the vibrational levels (including the zero level).
The thermal distribution of the rotational levels (unlike that of the vibrational levels) is not simply given by the Boltzmann factor.
We must account for the fact that, according to quantum theory, each state of an atomic system with total angular momentum J
consists of 2J+1 energy levels which are degenerate in the absence of an external field; that is, the state has a (2J+1)-fold
degeneracy. The number of molecules NJ in the rotational level J of the lowest vibrational state at the temperature T is thus
proportional to
F hc
(J )

k T
NJ ∝ (2J + 1) e B (6.2.1)

E(J)
where F(J) is the rotational term value (equal to hc ). For sufficiently large T or small rotational constant B, the population of the
Jth rotational level, for N molecules, is given by the following equation [see p. 125 in reference (3)]
B (J +1)hc
J
hcB −
NJ = N (2J + 1) e kT (6.2.2)
kT

Since the factor 2J+1 increases linearly with J, the number of molecules in the different rotational levels goes through a maximum
before decreasing with the rotational quantum number. It can be shown3 that this maximum lies at
−−−− − −−
kT 1 T 1
Jmax = √ − = 0.5896√ − (6.2.3)
2Bhc 2 B 2

This expression clearly shows that Jmax increases with decreasing B and increasing T. It should be noted that the number of
molecules in the lowest rotational level, J = 0, is not zero. With increasing temperature, the band extends over a wider frequency
range and the intensity maxima of the two rotational branches move outward and at the same time become flatter.

6.2: Transition Intensities is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

6.2.1 https://chem.libretexts.org/@go/page/372591
6.3: Part I - Experimental
We have a Nicolet iS50 FT-IR Fourier Transform Infrared (FTIR) Spectrometer. Below are instructions for its use.

 Note
In preparation for this experiment, you should also read your textbook (McQuarrie, D.A. and Simon, J.D. "Physical Chemistry:
A Molecular Approach", University Science Books, CA, 1997, Sections 13.2-13.4, 18.4-18.5); Atkins, P., Physical Chemistry,
5th ed., sections 16.8 – 16.11; and, Silbey, R.J., Alberty, R. A. and Bawendi, M.G., Physical Chemistry, 4th ed., Sections 13.4,
13.6 and 13.7, for a discussion of rotation-vibration spectroscopy of diatomic molecules. Each pair of students may share some
of the data analysis as described in Section 4.3.3.

Preparation of Gas Sample Cell


A mixture of HCl and DCl will be introduced into an IR gas cell. Note that because natural chlorine is a mixture of two isotopes
Cl and Cl (in a ratio of 3 to 1), this sample is a mixture of four different gases, each with a distinct spectrum.
35 37

Hydrogen chloride and deuterium chloride are colorless, pungent, and corrosive gases that fume in air. Contact will cause
burns to the skin, severe burns to the eyes and burns to the respiratory system if inhaled. A vacuum line will be used to
transfer the gases from their high-pressure containers to an evacuated IR cell.
When not in use the IR gas cell must be stored in the desiccator. The cell is constructed of a glass body with KBr windows at
each end. Ensure that the salt windows are kept free of moisture at all times. Please do not touch these windows.

 Vacuum-Line Start-Up
Most of the following will be completed by the TA before the lab period.
1. Attach an empty N2 trap to the main manifold. Hold this trap in place until a vacuum has formed.
2. Turn on the floor vacuum pump to start generating a vacuum.
3. Turn on all valves in the manifold that are marked with purple tape.
4. Turn on the diffusion pump.
Ensure that the floor pump is on BEFORE the diffusion pump. Otherwise, the diffusion pump will generate heat and the
oil within the pump will oxidize, causing long-term damage to the diffusion pump.
5. Turn on Digivac monitor. Use this to actively monitor the pressure at the end of the main manifold.
6. Check the analog pressure gauge on the main manifold to ensure the pressure is stable.
Throughout the experiment, the pressure should remain below 30 mtorr if the system is closed and stable.
7. Fill a 4L dewar with liquid N2 from the shared Chemistry department supply.
8. Place a second dewar underneath the N2 trap and fill with liquid N2 from the 4L dewar such that the N2 trap is submerged
in liquid N2
9. Check the level of N2 as the experiment progresses. Refill if necessary to ensure the trap remains cold.
10. Check to ensure the pressure at the end of the main manifold is below 30 mtorr before proceeding.

Transferring HCl and DCl to gas cell


1. Remove gas cell from the desiccator and attach to the middlemost section of the working manifold shown below in Figure
6.3.1.

6.3.1 https://chem.libretexts.org/@go/page/372593
Figure 6.3.1 : Vacuum line set-up: working manifold and main manifold
2. Evacuate the back of the gas cell and introduce a vacuum to the interior of the cell.
3. Remove cell from manifold, bring it to the FTIR instrument, and take a background spectrum (see "Obtaining an FTIR
Spectrum using the Nicolet iS50" below).
4. Replace the gas cell onto the working manifold and open it to vacuum.
5. Connect the pressurized HCl cannister to the working manifold via a vacuum-rated line. This line should contain its own needle
valve that may be closed to separate the cannister from the rest of the working manifold. Introduce a vacuum to the line and to
the pocket between this needle valve and the main valve on the cannister. The HCl cannister should remain closed at this time.
Close the needle valve once this pocket has been evacuated.
6. Close the working manifold off from the main manifold.
7. Introduce a small amount of HCl gas to the pocket between the cannister and the needle valve by turning the cannister's main
valve 180 degrees open briefly, and returning 180 degrees to close the cannister once more.
8. After ensuring the HCl cannister is closed, open the needle valve to release HCl into the working manifold. Use the mercury
pressure gauge on the working manifold to monitor this process.
If the pressure rises such that it threatens to expel mercury from the pressure gauge, open the working manifold to the main
manifold to release pressure.
9. Slowly release HCl from the working manifold using the valve to the main manifold until the pressure in the gas cell reaches 10
mmHg.
10. Close the gas cell off from the working manifold and open the working manifold to the main line to release any excess HCl
from the rest of the system. (Remember to include the pocket between the needle valve and the HCl cannister!)
11. Repeat steps 5-10 using DCl gas to introduce an additional 10 mmHg DCl into the gas cell such that the final contents of the
cell are 10 mmHg HCl and 10 mmHg DCl. (What pressure DCl needs to be introduced to the working manifold for this to
happen? Hint: it's not 10 mmHg.)
12. Record HCl/DCl spectrum using "Obtaining an FTIR Spectrum using the Nicolet iS50" below.
13. Once all data has been collected, return the gas cell to the working manifold and safely evacuate all remaining HCl/DCl.
14. Close and remove the evacuated gas cell and return it to the desiccator.

 Vacuum Line Shutdown


1. Turn off the diffusion pump. Wait 5 mins for the pump to cool.
The pump should be cool before any oxygen is introduced to the system to prevent oil in the pump from oxidizing.
2. Close all valves with purple tape.
3. Turn off Digivac monitor.
4. Turn off the floor pump.
5. Pick an external valve and open it to release the vacuum and open both the working manifold and the main manifold to the
air.
Do this quickly after step 4. The system should not be left under vacuum without an active vacuum pump.
6. Remove N2 trap and transfer it to the fume hood to safely evaporate excess HCl/DCl gas.

6.3.2 https://chem.libretexts.org/@go/page/372593
1. Do this quickly after step 5. Leaving the N2 trap cold after the vacuum has been released may cause liquid O2 to
condense, which is highly volatile.

Obtaining an FTIR Spectrum using the Nicolet iS50


1. Check the sample holder in the spectrometer to make sure there is nothing blocking the path of the infrared beam. Close the
cover of the instrument.

2. The Omnic software may already be loaded. If not, load it by clicking the icon on the Windows desktop. As the software
loads you will a blue light appear on the right side of the instrument and the main window with the menu items and button bar.

Figure 6.3.1 : The Omnic Button Bar. (CC-NC-BY-DUKE CHEM)

 Note on Collecting Data

This is a single beam instrument. You must first record a background interferogram of the evacuated gas cell. A fast Fourier
transform (FFT) will then be performed to give background transmittance spectrum. You will then introduce your gas
sample into the SAME gas cell and record a new interferogram. FFT analysis of this interferogram will give a single beam
transmittance spectrum of your sample. This spectrum is then ratioed against the single beam background spectrum to give
the transmittance spectrum of your sample. Finally the absorbance spectrum of your sample is calculated and displayed.

3. Open the sample compartment of the FTIR spectrometer and mount the evacuated gas cell on the gas cell holder. Close the
cover and allow the sample compartment to be purged with dry, CO2-free air for at 5 minutes, noting the time so that the purge
time is reasonably accurate.
4. Check the Experiment Setup from the Collect menu. It should be set to run 8 scans at 0.125 Cm-1 resolution. You should be
working in Absorbance format and a background scan should be “good” for 200 minutes.

Figure 6.3.2 : The Collect Pane of the Experiment Setup Window. (CC-NC-BY-DUKE CHEM)
Click on the Bench tab and wait for the interferogram to appear. The Max signal should be around 4.5. Click OK.

6.3.3 https://chem.libretexts.org/@go/page/372593
Figure 6.3.3 : The Bench Pane of the Experiment Setup Window. (CC-NC-BY-DUKE CHEM)
5. Select Collect Background from the Collect menu (or with the Col Bkg button). Click OK when you see the Confirmation
box. It will take about 4 minutes for the background scan to complete. Part way through, the background scan will be displayed
on the main software screen. Once complete, add the spectrum to the window when prompted. Note the presence of water vapor
and CO2 in the spectrum. This arises from the atmosphere surrounding the cell (Remember: you gas cell is evacuated at this
point). Keep a copy of the background both as a file (File → Save as), and take a screen shot for your ELN.

Figure 6.3.4 : A Typical Background Spectrum. (CC-NC-BY-DUKE CHEM)


6. Remove the evacuated gas cell from the sample compartment. Under the guidance of your TA, introduce first 1 cm Hg pressure
of DCl gas followed by about 0.8 cm Hg pressure of HCl gas into the cell. You must be careful to avoid the introduction of air
into the cell. Replace the gas cell on the sample holder, and close the door of the sample compartment. Allow the sample
compartment to purge for the same time as in the background spectrum.
7. Select Collect Sample from the Collect menu (or use the Col Smp button). You can enter a title or leave the default date/time
when prompted. Click OK, and OK again in the Confirmation box to start the run.
8. When the scan is complete, add the scan to the window and you will see both sample and background displayed. Save the data
file (File → Save as).

6.3.4 https://chem.libretexts.org/@go/page/372593
Figure 6.3.5 : Sample and Background Together. (CC-NC-BY-DUKE CHEM)
Click on the background scan and select Hide Spectra from the View menu to remove the background spectrum from view.

Figure 6.3.6 : The Sample Spectrum. (CC-NC-BY-DUKE CHEM)


9. You will see two (or three) sets of peaks in the region 3200-1800 cm–1. Which of these is the spectrum of HCl and which the
spectrum of DCl? Explain your reasoning. In your report, explain the origin of the third set of peaks (if observed). Click and
drag a box around the H35Cl/H37Cl absorption and double-click in the box to expand the scale so that the spectrum of fills the
screen.

10. Select Find Peaks from the Analyze menu or click the button. Move the horizontal threshold bar down so that it captures
as many of the peaks as possible without getting into the noise at the baseline (see Figure 4.3.7 below). You can also adjust the
sensitivity slider at the left of the screen to fine tune the peak selection. When you have a good selection of peaks, click the
Clipboard button to put the peaks table into the Windows Clipboard. Open Excel and paste the clipboard into a spreadsheet.
Save the spreadsheet along with a screenshot of the peak positions overlaid onto the speactum as seen in Figure 4.3.7.

6.3.5 https://chem.libretexts.org/@go/page/372593
Figure 6.3.7 : the Find Peaks Function. (CC-NC-BY-DUKE CHEM)
11. Click Replace at the top right of the screen to bring you back to the active window. Click and drag the white selection window
over to the DCl spectrum.

Figure 6.3.8 : The Selection Window. (CC-NC-BY-DUKE CHEM)


Re-scale the DCl spectrum in the window and repeat the peak-finding steps being sure to add the resultant peaks table to your
Excel spreadsheet.
12. Save screen shots of the spectra and the report for your ELN.
13. If you wish to record the spectrum of a new sample, say of SO2 in Part 2 of this experiment, repeat steps (4.3.2.6)-(4.3.2.12).
14. When you are done with the instrument and all of your work is saved, close the software. You will notice the blue light on the
instrument turn off as the instrument goes into standby mode.
15. Remove your sample from the sample compartment and close the cover.
16. Re-evacuate the gas cell to remove your sample and store the evacuated cell (with stopcocks closed) in the desiccator.

6.3: Part I - Experimental is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

6.3.6 https://chem.libretexts.org/@go/page/372593
6.4: Part I - Data Analysis
Data Analysis for HCl/DCl FTIR spectra

 Note
This is one of the more laborious data analyses you will do this semester. Let's minimize the number of repetitive tasks to save
you some time. The following shortcuts are recommended:
Collaborate with a partner on steps 1-3 (allowed in this specific case only):You may collaborate with your lab partners
in steps 1-3 below. One student should work up the data (steps 1-3 ONLY) for H35Cl and H37Cl, and the second student
should do the same analysis for D35Cl and D37Cl. You should then share your graphs of Δν̃ versus m and the values
(m)

determined for the slopes and y-intercepts of these plots. Each student should do one complete analysis (either HCl or DCl)
alone, and the other analysis can be the work of a partner. However, you are responsible for all work that is in your ELN, so
be sure your partner's work is correct. *Steps 4 and onward should be completed individually, as usual.
Calculate errors for one isotope only: The propagation of errors that you calculate for one isotope should similar (perhaps
not identical, but close enough) as for another isotope. Save time by propagating errors carefully for (B , α , ν̃ , Ie, re, and
e e o

k for one isotope, then use those errors as an approximation for the errors for analogous values for other isotopes.
Use what you know! Do you remember that MatLab Introduction module way back in the beginning of this course? It was
a long time ago, but it was a template for data analysis in this module! Go back and use what you already know, and what
you already created as a template, for this data analysis.

Collaborate with a partner


1. Label each peak in your spectra as "P" or "R" and its J-value as: P(1), P(2),..., R(O), R(1),... etc., as shown in Figure 2 on page
399 of Shoemaker, Garland and Nibler. (Note: Spectra recorded using your iS50 FTIR spectrometer are plotted with the
wavenumber (abscissa) decreasing from left to right.)
2. For each of the four isotopic combinations (H35Cl, H37Cl, D35Cl and D37Cl), make a table of the m values and the
corresponding peak frequencies ν̃ . For each isotopic combination, calculate the separation between adjacent peaks.
(m)

Δν̃ (m) = ν̃ (m+1) − ν̃ (m) (6.4.1)

 Note
cannot be calculated for m = –1 since there is no peak for m = 0. Therefore m = –1 must be deleted from the vector
Δν̃ (m)

of m values.

3. Plot the differences between adjacent absorption frequencies (i.e., Δν̃ (m) versus m. Then perform a linear least squares fit of
the data with the understanding that

Δν̃ (m) = ν̃ (m+1) − ν̃ (m) = (2Be − 3αe ) − 2αe m (6.4.2)

Use the standard deviations of the slope and intercept in your error analysis (Q14).

Work on your own


4. Compute B and α for each isotope from the intercept (2B − 3α ) and the slope −2α of each best fit line (Equation
e e e e e

6.4.2).

5. Calculate the frequency of the forbidden transition ν̃ using your calculated values of B and α . To do this, use equation 9
o e e

in expt 37 of Shoemaker, Garland and Nibler. The equation is ν̃ = ν̃ + (2B − 2α )m − α m . Calculate ν̃ from several
m o e e e
2
o

of the lines with low m values. Calculate the average value, and standard deviation, of ν̃ . o

 Note
You will get the best results by using only the lines closest to the center (e.g., m = ±1, ±2).

6. Create a table of ν̃ , B , and α for the four isotopic combinations The units should be cm–1.
o e e

6.4.1 https://chem.libretexts.org/@go/page/372937
7. Calculate the moment of inertia (I ), and the internuclear distance (a.k.a. equilibrium bond length, r ) for all four
e e

isotopic combinations.
8. Calculate the harmonic oscillator force constant, k, for each isotopic combination from your calculated values of ν̃ . In o

doing so you will need to assume that ν̃ and the harmonic oscillator frequency ν̃ are the same, thus ignoring the effects of
o e

anharmonicity (see the formulas and discussion on pages 400 and 401 of Shoemaker, Garland and Nibbler).
9. Compare your values of re and k with those reported in the literature. Why are these values independent of isotopic
substitution?
∗ ∗
ν̃ B
10. Determine the isotope effects: Compute the ratios o
and Be
e
for 35Cl and 37Cl in the case of HCl, and again in the case of
ν̃ o

DCl. Compare these with the predicted values given by Eqs. (11) and (13) on page 400 of Shoemaker, Garland and Nibbler.
11. From the results in 10, above, explain why is the 35Cl /37Cl isotope effect so much larger in DCl than it is in HCl? That is,
when looking at your FTIR spectrum, the 35Cl/37Cl line separation is significantly greater in the DCl case then when compared
to the HCl case. A descriptive answer is a good start, but your experimental observations should also be supported with a
calculation that justifies the experimental observation.
12. Apply the Boltzmann distribution to explain observed band intensities in terms of the populations of the ground-state levels (see
section 4.2). Again, a descriptive answer is a good start, but it should be supported with calculated results. You might create a
spreadsheet to do these calculations. Does a calculation based on the Boltzmann population of rotational states match what you
see experimentally in terms of which rotational lines show the greatest intensity?
13. From the mean line intensities for at least 10 features (for example five or more lines from the HCl spectrum and five or more
from the DCl spectrum for a total of at least 10), determine the relative abundance of the two Cl isotopes. Estimate the weighted
average atomic mass of Cl (that you might see on a periodic table for example) assuming that 35Cl and 37Cl have masses of
34.968 and 36.956 amu, respectively. Compare your result with the reported atomic mass of Cl.
14. For one isotope only, calculate the estimated uncertainty in each of the parameters determined for this experiment (B , α , ν̃ , e e o

Ie, re, and k). (Another table would be helpful)


15. Why is it important to purge the sample compartment of the FTIR spectrometer for both the background and sample runs? Why
is it important to purge both background and sample for the same amount of time?

6.4: Part I - Data Analysis is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

6.4.2 https://chem.libretexts.org/@go/page/372937
6.5: Part I - References
1. Moore, W.J., Physical Chemistry, 4th ed., Prentice-Hall, Englewood Cliffs, New Jersey, 1972, pp. 767-770; or Pauling, L.;
Wilson, E.B. Introduction to Quantum Mechanics, McGraw-Hill, New York, 1935.
2. Pitzer, K.S. Quantum Chemistry, Prentice-Hall, Englewood Cliffs, New Jersey, 1953.
3. Herzberg, G. Molecular Spectra and Molecular Structure I. Spectra of Diatomic Molecules, 2nd ed., Van Nostrand, Princeton,
New Jersey, 1950, Chapter III.
4. Hill, T.L. Introduction to Statistical Thermodynamics, Addison-Wesley, Reading, MA 1960, Chapter X.
5. Prais, M.G. J. Chem. Ed., 1986, 63, 747.
6. Shoemaker, D.P., Garland, C.W., Nibler, J.W. Experiments in Physical Chemistry, 6th ed., McGraw-Hill, New York, 1996,
Chapter XIV, Experiments 35 and 37.
7. Richards, L.W. J. Chem. Ed., 1966, 43, 552. See also, McQuarrie, D.A. and Simon, J.D., Physical Chemistry: A Molecular
Approach, University Science Books, CA, 1997, p. 169.
8. Sime, R.J., Physical Chemistry: Methods, Techniques and Experiments, Saunders, Philadelphia, PA, 1990, pp. 676-687.
9. Silbey, R.J., Alberty, R.A. and Bawendi, M.G., Physical Chemistry, 4th ed., Wiley, NY, 2005, Chapter 13, Section 13.4, 13.6,
13.7.
10. Atkins, P.W. Physical Chemistry, 5th ed., W. H. Freeman, NY, 1994, Sections 16.8 – 16.11., or 6th ed., W. H. Freeman, NY,
1997, Sections 16.9 – 16.12.

6.5: Part I - References is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

6.5.1 https://chem.libretexts.org/@go/page/372968
6.6: Part II - Experimental
 Reminder

The general approach to Part II of this experiment is adapted from D. P. Shoemaker, C. W. Garland, and J. W. Nibler,
Experiments in Physical Chemistry, 6th edition, McGraw Hill Co. Inc, NY, 1996. In preparation for this experiment please read
Experiment 35 in Shoemaker, Garland and Nibler beginning on page 383. You will want to read and understand the theory
section including the subsections on the valence-force model and the vibrational partition function.

4.6.1. A sample of SO2 gas will be introduced into an IR gas cell.


4.6.2. A vacuum line will be used to transfer the toxic SO2 gas from its high-pressure containers to an evacuated IR cell.

 Note
You will use a pressure of SO2 less than or equal to 1 cm Hg. If the pressure of SO2 in the gas cell is much higher, you will not
be able to observe the fine structure on the most intense fundamental frequency. If the gas pressure is too low, you may not be
able to observe the weak combination bands. Depending on your skill with transferring known pressures of gases on the
vacuum line, you may find it advantageous to record the spectrum using two different gas pressures so that all the features of
the spectrum are well resolved (consult your TA).

4.6.3. Record an FTIR absorbance spectrum of about 0.5-1.0 cm Hg Pressure of SO2 at 0.125 cm–1 resolution, and under the same
conditions used for the HCl/DCl sample, and ratio against the background spectrum of the evacuated cell.
4.6.4. Plot the absorbance spectrum over the range 4000-400 cm–1 (or a range suggested by your TA). Determine the frequency of
each band seen in the spectrum. Search for combination and overtone bands.
4.6.5. Use the Omnic Annotation Tool (look to the bottom of the spectrum window for the button) to label the center position of
each fundamental band and of any overtone or combination bands you observe. You use this tool by clicking on the button,
position the cursor (now with an associated T) at the point you want to measure. Then click and drag a line up (maybe up and
over). When you click again the wavelength at that position on the spectrum will be added as an annotation to the spectrum (see
picture). If you want to delete an annotation, drag a box around the annotation (the actual wavelength added to the screen) and
select Delete Annotation from the Edit menu.

Figure 6.6.1 : The Symmetric Stretch of SO2(g). (CC-NC-BY-DUKE CHEM)

6.6.1 https://chem.libretexts.org/@go/page/372965
 Note

Notice that in determining the frequency of vibration in you are estimating the center of the feature by eye and hand, rather
than having the Omnic software determine the peak position. To get an estimate of the uncertainty of your determination, do
between three and five repeated trials of estimating the frequency and find the mean and standard deviation.
Notice also that because we are using both an instrument and a gas cell with KBr windows, we will not be able to see the
bending mode clearly or completely as the KBr starts to absorb strongly in that region of the spectrum. You should, however,
see enough of the bending mode to get a crude estimate of its frequency. Get help from your TA, and comment in your ELN
why this frequency may show greater deviation from expected values than the others.

4.6.6. When you are done, please remember to evacuate the gas cell and store the evacuated cell in the desiccator.

6.6: Part II - Experimental is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

6.6.2 https://chem.libretexts.org/@go/page/372965
6.7: Part II - Data Analysis
4.7.1. Assign the fundamental bands of SO2 as: symmetric stretch, ν̃ 1 ; bend, ν̃ 2 ; asymmetric stretch, ν̃ 3 ; and, report their
frequencies. Compare your value with those reported in the literature.
4.7.2. Then assign any other bands, comparing the observed frequencies with those calculated from combinations of the
fundamental values ν̃ , ν̃ , ν̃ .
1 2 3


4.7.3. Calculate k1 and 2
from Equations (1) and (3) (on page 385 of Shoemaker, Garland and Nibler) using your values of ν̃ 1 ,

, and ν̃ . Check the valence-force model by calculating both sides of Eq. (2) (also on p. 385 of Shoemaker) and comparing them.
ν̃ 2 3

Again, compare your results with literature values.

 Note
When finding literature values for the force constants of SO2, you will likely notice that different notations are used. It may be

difficult to find the k1 and 2
notation used in Shoemaker, Garland and Nibler (and as such in this manual). For example,

f dynes
Kivelson8 reports these same vales as fd ; and, a

2
and reports the values in units of cm
rather than N

m
. Be prepared to
d

convert.

~
4.7.4 Using your values of ν̃ 1 , ν̃ 2 , and ν̃ 1 , calculate C v(vib) at 298 K and at 500 K from Eq. (8) on page 386 of Shoemaker,
Garland and Nibler.
~
4.7.5 At these same temperatures, calculate C from the expression v

~ ~
C v = 3R + C v(vib) (6.7.1)

~ ~
and compare your results with the following C values (obtained from directly measured values of C and the expression
v p

~ ~
Cv = Cp − R (6.7.2)

~
Cv = 30.5 J K–1 mol–1 at 298 K and 37.7 J K–1 mol–1 at 500 K.
4.7.6 Calculate the uncertainty of each result you report in this experiment.

6.7: Part II - Data Analysis is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

6.7.1 https://chem.libretexts.org/@go/page/372971
6.8: Part II - References
1. Herzberg, G. Molecular Spectra and Molecular Structure II. Infrared and Raman Spectra of Polyatomic Molecules, Van
Nostrand, Princeton, New Jersey, 1945, pps. 168-172, 251-269, 285.
2. Herzberg, G. Molecular Spectra and Molecular Structure II. Infrared and Raman Spectra of Polyatomic Molecules, Van
Nostrand, Princeton, New Jersey, 1966, p. 605.
3. Atkins, P.W. Physical Chemistry, 5th ed., Freeman, New York, 1994, Chapter 14; McQuarrie, D.A. Statistical Thermodynamics,
Harper & Row, New York, 1976.
4. Lewis. G.N., Randall, M. (revised by Pitzer, K.S. and Brewer, L), Thermodynamics, 2nd ed., McGraw-Hill, New York, 1961, p.
419ff.
5. Silbey, R.J., Alberty, R.A., and Bawendi, M.G., Physical Chemistry, 4th ed., Wiley, NY, 2005, Chapter 13, Sections 13.6, 13.8,
and 13.10.
6. Colthup, N.B., Daly, L.H., and Wiberley, S.E. Introduction to Infrared and Raman Spectroscopy, 3rd ed., Academic Press,
Boston, 1990.
7. Shoemaker, D.P.; Garland, C.W.; Nibler, J.W. Experiments in Physical Chemistry, 6th ed., McGraw-Hill, New York, 1996,
Chapter XIV, Experiments 35 and 37.
8. Kivelson, D., The Determination of the Potential Constants of S02 from Centrifugal Distortion Effects, J. Chem. Phys., 22 (5),
904, 1954.

6.8: Part II - References is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

6.8.1 https://chem.libretexts.org/@go/page/373182
CHAPTER OVERVIEW

7: Molecular Electronic Structure Calculations


7.1: Pre-Lab Assignment
7.2: Introduction to Electronic Structure Calculations
7.3: Accessing Spartan
7.4: Exercise 1 - Estimation of the IR Spectrum of HCl and DCl
7.5: Exercise 2 - Molecular Orbitals of Formaldehyde
7.6: Exercise 3 - Modeling Sulfur Dioxide - Comparing the Results of Different Basis Sets and Calculation Models
7.7: Exercise 4 - Modeling the Absorption Spectrum of Cyanine Dyes
7.8: Exercise 5 - Modeling an Intramolecular Rearrangement

7: Molecular Electronic Structure Calculations is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

1
7.1: Pre-Lab Assignment

Pre-Lab Assignment
This module is an introduction to computational chemistry and molecular electronic structure calculations using Spartan software.
These computational exercises circle back around to some of the "wet lab" modules you completed earlier in the semester.
To help you prepare for this final experiment we ask that you do the following:
1. Read the introduction to the lab, and the chapters in your textbook on quantum mechanics and electronic structure theory for
molecules.
2. Briefly read the Guidelines for Using Spartan® Software and look over some of the helpful materials found on your Sakai site.

Additional pre-lab assignment for Chem310L Independent Project


For CHEM 310L students, the Spartan exercises will be followed by an independent project in which you can choose to do a
computational study.
We ask that you begin to think about your independent project, and we want you to do something you are really interested in! In
310L, you can choose a topic for your independent project, and choose whether you would like to dive into more computational
work, or explore other applications of physical chemistry in the peer-reviewed literature. You should begin searching the literature
(eg using Web of Science or looking at the Table of Contents for leading Physical Chemistry journals) for physical chemistry topics
that interest you. Alternatively, you can think about simple chemical phenomena or reactions that you are interested in and that you
can model using Spartan. Come to the laboratory session prepared with a brief written proposal (1 paragraph) on at least two
ideas for what you might do for your independent project and be prepared to discuss it with your instructors.

Introduction
Computers and computational methods have become major tools of chemical research. This comes from the availability of
powerful workstations and software packages at affordable prices, the development and implementation of accurate and efficient
computational theory, and the insight and user-friendly interface provided by scientific visualization.
In this laboratory exercise, you will use the Spartan® software to explore the calculations of properties of isolated molecules based
on a whole range of modeling methods, from empirical molecular mechanics to accurate ab initio quantum mechanical methods.
You will carefully study the accuracy of each method by comparison with experimental results and with other computational
methods. Understanding the limitations of modeling methods is as important as knowing their availability and accuracy. We will
learn how such techniques can be used to help solve chemical problems.
There are two parts to the Spartan Exercise:
1. The first part of the exercise is designed to relate to the previous "wet lab" modules of this course.
2. The second part focuses on using computational methods to gain knowledge that is difficult to get from experimental methods;
that is to study frontier molecular orbitals and transition state structure in chemical reactions.

 Note

This lab is split into 5 separate exercises designed to illustrate the various capabilities of Spartan. The requirements for your
lab notebook are stated within each exercise instead of at the end of the chapter. It is recommended that you follow along with
the ELN template while completing this lab to ensure that all data is properly recorded.

Independent Project
This is for Chemistry 310L students only (301L students don't need to do this). After completing the Spartan excecises, you will
porpose and complete an independent project - for this project, you are welcome to model your own favorite molecules or
reactions. You will need to design this part of the project before you come to the laboratory and learn to ask the right questions with
the proper tools. This is an opportunity to explore something that is interesting to you.

7.1: Pre-Lab Assignment is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

7.1.1 https://chem.libretexts.org/@go/page/371809
7.2: Introduction to Electronic Structure Calculations
Calculation of the electronic structure of a molecule (i.e. the molecular orbitals) can be accomplished by the following steps:
1. Preparation of the input data: Within the Born-Oppenheimer approximation, nuclei are held fixed while solving for the
electronic motion. Before you begin your calculations, you must first define the fixed nuclear positions at which the electronic
structure calculation will be performed -- i.e. you must define a molecular geometry. In Spartan®, this is done graphically in
much the same way in which you built molecules when you took organic chemistry courses.

Many of the calculations you undertake here will use a variational approach in which the molecular wavefunction is written as a
single Slater determinant of molecular orbitals (MO), where each MO, in turn, is written as a Linear Combination of Atomic
Orbitals (LCAOMO). To begin this variational approach, you must obtain a starting "trial wavefunction." In the LCAOMO
method, this requires definition of the atomic functions that will be used to construct the MO's. Recall the variational approach
to the electronic structure of H2 -- to begin we used a trial LCAOMO wavefunction Ψ written as a linear combination of a 1s
atomic orbital on HA (Φ1sA) and a 1s atomic orbital on HB ((Φ1sB); i.e., we defined the atomic functions that were used to
construct the MO's. In Spartan® you will find a reasonable collection of such sets of atomic functions built into the software,
and consequently, selection of an appropriate set of atomic functions is usually achieved by simply choosing from a selection of
available methods and basis sets.
2. Running the calculation: After building the molecule of interest and selecting the type of calculation to be run, and other
relevant parameters, one submits this job for calculation.
3. Analysis of the results: This may involve such tasks as viewing electron densities, exploring the bonding or antibonding
character of molecular orbitals, animating vibrational modes associated with nuclear motion, viewing energy-minimized
structures and determining bond lengths and bond angles, or calculating molecular properties. Most of this analysis is done
using the graphical interface available in Spartan® thereby minimizing the need to examine multiple pages of numerical data!
Again, Spartan® simplifies the analysis by providing numerous visualizations and menu-driven data displays.

7.2: Introduction to Electronic Structure Calculations is shared under a not declared license and was authored, remixed, and/or curated by
LibreTexts.

7.2.1 https://chem.libretexts.org/@go/page/371810
7.3: Accessing Spartan
A network version of Spartan 20 is available from the chemistry server, bohr. Access to Spartan through FastX is the same as
access to Matlab, except that you will want to type...

spartan

at the prompt after logging in. Every computer in the physical chemistry lab is already set up with FastX. In addition, you will
likely want to have access from your own laptop. Complete instructions for obtaining and configuring FastX, and the Duke VPN (if
needed) can be found here.
1. Loading – Spartan® 20

1.1. Click the FastX shortcut icon, on the main screen of any computer in the Advanced Labs. Click on the
bohr.chem.duke.edu session. Enter your NetID and password and click Continue. Click on the + at the upper right of the new
screen that appears, and click on Xterm. An xterm window should appear. Find your TA or your Instructor if there are problems
connecting; it happens sometimes.
2. Type

spartan

This will start the Spartan program and the main screen will appear.
Notice the menus and buttons at the top of the screen.

Figure 7.3.1 : The Spartan Menu and Button Bar. (CC-NC-BY-DUKE CHEM)
You will use these menus and button to operate Spartan®. Note also the Help menu to the right of the menu bar. Among the many
useful things included in the Help menu is a summary of mouse commands (within General Operating Features). You will use this
summary frequently as you do your work.

Figure 7.3.2 : Spartan Help. (CC-NC-BY-DUKE CHEM)


Once in General Operating Features, scroll down to find the table of mouse functions.

7.3.1 https://chem.libretexts.org/@go/page/371811
Figure 7.3.3 : Using the Mouse. (CC-NC-BY-DUKE CHEM)

 Note

Spartan works best with a 3-button mouse. However it is easy enough to translate the mouse instructions to work equally well
with a trackpad.

7.3: Accessing Spartan is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

7.3.2 https://chem.libretexts.org/@go/page/371811
7.4: Exercise 1 - Estimation of the IR Spectrum of HCl and DCl
Exercise 1 - Estimation of the IR Spectrum of HCl and DCl
The goal of this first exercise is to compare your experimental FTIR data from the HCl and DCl diatomic molecules with
computational predictions.

1. Select New Build from the File from the menu, or click the button.
2. The right side of the screen is the model kit.

Figure 7.4.1 : The Model Kit. (CC-NC-BY-DUKE CHEM)


Notice that you are in the Organic area of the model kit. Here you will see several common atom types and their hybridization.
You might take a moment now to familiarize yourself with the other screens of the model kit, Inorganic, Substituent, Peptide and
Nucleotide.
3. From within the Organic area of the model kit, click on the –Cl atom and then double-click somewhere in the working area of
the main Spartan screen. You will see a colored line, with the Cl atom at the end. Select Save as from the File menu, or click the

button. Keep the suggested filename, or give the file a name of your choice. A hydrogen atom is automatically added to

your molecule. To see the hydrogen, click the button to leave the Build Mode and enter the View Mode.

 Note

Clicking the button returns you to the Edit Build Mode.

The exact view of the molecule will depend upon what “model” system you have selected. Use the Model menu to change how the
molecule looks.
4. Select Calculations from the Setup menu. Choose Equilibrium Geometry at Ground State in Gas, with Hartree-Fock, 3-
21G. Total charge should be Neutral, and Unpaired Electrons should be 0. Select compute: IR (with Current Model). Click on
Submit. If prompted, save the file (use the same file name). If everything is correct, you will see the “Job “filename” has been
submitted” box. Click OK.
5. After a moment a “Job “filename” completed” box appears on the screen. Click OK. Inspect the output of this job using
Display, Output. The Summary tab will give you the major results of the calculation. If you expand the IR section of the

7.4.1 https://chem.libretexts.org/@go/page/371812
Summary you will see the calculated estimate of the frequency of vibration. Record this value in your lab notebook. In the
Output tab, you will see information about the chosen basis set, the number of steps for geometry optimization, and other
information related to the calculation. Close the Output dialog box. Then select Spectra under the Display menu and click on

the , and select IR Calculated. The calculated spectrum appears on the screen. Click on the and you will again see the
harmonic oscillator frequency of HCl molecule.
6. To animate the vibrational motion, put a check in the box next to the frequency (or click on the peak in the spectrum). If you
like, click on the again and add the Experimental spectrum to the window. How well does the Hartree-Fock calculation
reproduce the experimental IR spectrum? To stop the animation, click again on its frequency, or remove the check from the
checkbox next to the vibration. Close the spectrum window.

 Note

If the animated vibration is going too fast, the speed of the animation can be changed in the following way. Click on the
icon at the left of the spectrum display window. You will two editable buttons appear.

To slow down the animation, change the number of steps to a larger value (for example, 30).
Changing the Amplitude will alter the degree of displacement of the nuclei through the course of the vibration.

7. Back in the main window, under the Geometry menu choose Measure Distance, and find the H-Cl bond length by clicking on
the two atoms (or by clicking on the bond between the two atoms). Look to the lower right corner of the screen to see the bond
distance. Record the bond length in your lab manual and compare the result with your experimental and literature data. Click
the button to close your HCl file.
8. Make a new HCl molecule. Click the to enter the View Mode. Click once on the hydrogen atom to select it. Select
Properties from the Display menu. An Atom Properties Box appears. Change the Mass Number from 1 to 2 (Deuterium).
Close the Atom Properties box. Return to the Setup Calculations box and make configure the calculation to be the same HF
calculation you did on HCl (do not forget to check the IR option). Click on Submit. Save the file under a new name and let the
calculation run.
9. Find the harmonic oscillator frequency and the bond length for DCl. Compare your results with both your HCl results and your
experimental (FTIR) DCl data. Select File, Close to close your molecule when you are done.

Exp. ν Exp. Bond


Data for HCl and Calc. ν
(cm–1) Cal. Bond Length (Å) Length (Å)
DCl (cm–1)
(FTIR) (FTIR)

HCl

DCl

10. Summary: Compare your results to those you obtained from the FTIR experiment done earlier in this course. What can you
conclude about the accuracy of computing bond lengths and vibrational frequencies using this Hartree-Fock calculation? Why
is the bond length identical for both HCl and DCl?
If you are interested, and since the calculations are relatively quick, you might try the same calculation with a different basis set
(6.31-G* for example), and/or a completely different model (semi-empirical or density functional). Does the accuracy of the results
change?

7.4.2 https://chem.libretexts.org/@go/page/371812
 Note

This section of your lab notebook should include:


A table summarizing the results of the HCl/DCl calculations
Comparison of computational results to your experimental FTIR data and literature data
What can you conclude about the accuracy of computing bond lengths and vibrational frequencies using this Hartree-Fock
calculation?
Why is the bond length identical for both HCl and DCl?
(See your ELN template for details)

7.4: Exercise 1 - Estimation of the IR Spectrum of HCl and DCl is shared under a not declared license and was authored, remixed, and/or curated
by LibreTexts.

7.4.3 https://chem.libretexts.org/@go/page/371812
7.5: Exercise 2 - Molecular Orbitals of Formaldehyde
Exercise 2 - Molecular Orbitals of Formaldehyde
This experiment explores the capability of Spartan® to calculate Molecular Orbitals and construct corresponding surfaces. The
formaldehyde molecule is used as a model.

1. Click to start a new file.


2. Build H2CO. Select Carbonyl from the Groups menu and double-click in the working area. You may want to rotate the
molecule on the screen for a better view.

Figure 7.5.1 : Building Formaldehyde. (CC-NC-BY-DUKE CHEM)

3. Click the minimize button (lower right, below Clipboard).


4. Select Setup, Calculations from the menu. Choose Equilibrium Geometry at Ground State in Gas, with Hartree-Fock, 6-
31G*. Total charge should be Neutral, and Unpaired Electrons should be 0. Click OK.
5. Choose Submit from the Setup menu. Click on the OK button. Save the file with the suggested filename, or choose one of
your own.
6. Once the calculations have completed, determine the calculated equilibrium values for the C–O and C–H bond lengths, and for
the HCH and HCO bond angles, by using Geometry, Measure Distance/Measure Angle. Record these values in your lab
notebook. Compare your calculated values with literature values.

 Note
The Computational Chemistry Comparison and Benchmark Database (CCCBDB) at NIST is one suggested place to look for
this type of information. Choose an appropriate search string in your favorite web browser. For example: formaldehyde bond
lengths and angles NIST. Look for the hit directing you to the CCCBDB site.

H2CO C=O length HCO angle HCH angle

Calc Lit Calc Lit Calc Lit

7.5.1 https://chem.libretexts.org/@go/page/371863
Reference

7. Inspect the Output of this job using the button. Find the Molecular Orbital Energies section under the Summary tab.
Calculate the energy difference between the HOMO and LUMO and record this in your lab notebook.

8. Click the button to see an abbreviated MO diagram on the left side of the screen. From here you can click on an orbital
to view it on the molecule. In this way find and visualize the π-bonding MO in H2CO. You can change the appearance of the
orbital display by changing the Style from Solid to Mesh to Transparent, etc. (at the lower right corner of the screen).
9. In the same way, visualize the LUMO surface. How might you describe this molecular orbital?
10. Now visualize the HOMO. How might you describe this one? Take a look at some of the other orbitals and record any other
interesting observations (with images) in your lab notebook.

 Note

Not every MO is shown in the diagram generated by clicking the button. You can calculate the surface for any desired
MO by clicking the button, and More Surfaces. Then, from the Surface dropdown menu select HOMO{-},N or
LUMO{+},N, where N is how many orbitals above the LUMO or below the HOMO needed to get to the orbital you wish to
see. Click on Apply, and then OK. Put a check mark in the box next to the selected orbital and then return to the Spartan
window to see the result.

 Note

This section of your lab notebook should include:


A table recording your calculated bond lengths and angles for formaldehyde and a comparison to literature values
Calculation of the energy difference between the HOMO and the LUMO
Images and descriptions of the HOMO, LUMO, and any other interesting orbitals of your choice.
(See your ELN template for details)

7.5: Exercise 2 - Molecular Orbitals of Formaldehyde is shared under a not declared license and was authored, remixed, and/or curated by
LibreTexts.

7.5.2 https://chem.libretexts.org/@go/page/371863
7.6: Exercise 3 - Modeling Sulfur Dioxide - Comparing the Results of Different Basis
Sets and Calculation Models
Exercise 3 - Modeling Sulfur Dioxide: Comparing the Results of Different Basis Sets and Calculation Models
This experiment is designed to explore the performance of various basis sets, models of Ab initio calculations and Semi Empirical
calculations, and Density Functional calculations. The SO2 molecule is used as a model.
1. Start a New Build.
2. To introduce you to another approach to building a molecule, select the Inorganic tab of the builder.
3. Build SO2 as follows: Change the element to sulfur.

Figure 7.6.1 : Selecting the Element Sulfur From the Inorganic Area of the Builder. (CC-NC-BY-DUKE CHEM)

Select atom hybridization. Double-click in the working area. Rotate the molecule to make all the bonds visible. Select O from

the periodic table and Select atom hybridization. Click at the ends of the bonds to S to make .

Select double bond from the Bond Order buttons and double-click on both S–O bonds. Be sure they have changed to double
bonds.

4. Click the button to do a quick molecular mechanics minimization of the energy.


5. Select Setup, Calculations from the menu. Choose calculate: Equilibrium Geometry with Hartree-Fock, 6-31G*. Total
charge should be Neutral, and Unpaired Electrons should be 0. Select compute: IR. Click submit. Save the file when prompted
and run the calculation.
6. From the Summary section of the Output (Display menu), record the energy of the molecule in your lab notebook. Record the
frequencies of the SO2 vibrations from the Raman/IR Table section of the Output Summary in your lab notebook.
7. Go to the Display Spectra area and view the calculated IR spectrum as you did for HCl and DCl. Animate the vibrations by
highlighting the different frequencies. Assign frequencies as ‘symmetric stretch’, ‘bend’, or ‘asymmetric stretch’.
8. Close the spectrum window.
9. Measure and record the O–S–O angle and the S–O bond length.

7.6.1 https://chem.libretexts.org/@go/page/371866
10. Return to the calculation setup window and repeat the HF calculation with a 3-21G basis set.
11. Submit the calculation. Record and compare your results for frequencies, bond length, and the angle with the previous
calculation.
12. Repeat again using the STO3G basis set, the Semi-Empirical theories PM3, RM1, and AM1, and the Density Functional theory
ωB97X-D with 6-31G* basis set. It might be a good idea to Save as each time to save each calculation in a separate file. Or if

not, at least save a new file for each calculation model (HF, Semi-Empirical, Density Functional.

 Note

If you do not save each calculation under a new filename, the new calculation overwrites the old in the current file.

13. Fill out the table of Theories and Basis sets for each combination.
14. Compare your results with your experimental (FTIR) data for SO2 (CHEM 310 only), and literature data for frequencies, bond
length and bond angle. Which combination of theory and basis set gives the most accurate prediction? Close your molecule.

 Note

The NIST CCCBDB database is a good source for literature values for the SO2 molecule.

S–O bond Bond Angle,


Table for SO2 Theory Basis set 1
nu nu2 nu 3
length(Å) (˚)

Ab initio HF STO-3G

HF 3-21G

HF 6-31G*

Semi-empirical PM3 N/A

RM1 N/A

AM1 N/A

ωB97X-D, 6-
Density Functional
31G*

Your exptl. data

Literature data*

* Reference

 Note

This section of your lab notebook should include:


A table that records, for each basis set:
Calculated energy of SO2
Frequencies of SO2 oscillations
Assign each frequency ‘symmetric stretch’, ‘bend’, or ‘asymmetric stretch’. You only have to do this for one basis
set.
S-O bond length
O-S-O bond angle
Comparison of calculated results to experimental (CHEM 310 only) and literature values
Which combination of theory and basis set gives the most accurate prediction?
(See your ELN template for details)

7.6.2 https://chem.libretexts.org/@go/page/371866
7.6: Exercise 3 - Modeling Sulfur Dioxide - Comparing the Results of Different Basis Sets and Calculation Models is shared under a not declared
license and was authored, remixed, and/or curated by LibreTexts.

7.6.3 https://chem.libretexts.org/@go/page/371866
7.7: Exercise 4 - Modeling the Absorption Spectrum of Cyanine Dyes
Exercise 4 - Modeling the Absorption Spectrum of Cyanine Dyes
This experiment is designed to build and study the series of conjugated dyes used in the UV-Vis experiment earlier in the course.
1. The first dye you will build and work with is the 1,1’-diethyl-4,4’-carbocyanine cation. (See Figure 3.2.2 in Experiment 1).
Start a new file. Select the Inorganic tab of the Builder. Under Rings select Naphthalene. Click in the workspace to place a

naphthalene ring on the screen. Now select the sp2 hybridization option and check to see that the element in the element
2
selection button is Carbon. Move to the working area and add three sp carbons to the para-position of the naphthylene ring. At
this point your structure should look something like this.

Figure 7.7.1 : Adding sp Carbons to the Naphthalene Ring. (CC-NC-BY-DUKE CHEM)


2

Return to the Naphthalene selection under Rings again. Notice the small yellow circle on the structure in the Rings window

. This indicates the “active” valence of the naphthyl group. Make the para-position of naphthyl active by clicking on the
valence in the small structure window. Click on the free valence of the 3rd carbon atom in the working area, to add the second ring
to your structure. At this point it might look something like this.

7.7.1 https://chem.libretexts.org/@go/page/371872
Figure 7.7.2 : Adding the Second Naphthalene Ring. (CC-NC-BY-DUKE CHEM)
Now click on the delocalized bond . Double click on each C–C bond connecting the three carbons to each other and to the two
rings. Look closely to be sure that all C–C bonds are delocalized.

Figure 7.7.3 : Changing the Bond Order of the Connecting Carbons. (CC-NC-BY-DUKE CHEM)

Click on the button to minimize the structure. Chances are good that not all of the connecting carbons between the rings are
in the trans- orientation, and that the overall structure is not yet planar. You can rotate around specific bonds to alter the
conformation of the molecule. To do this, click on a bond you wish to rotate around. A red arrow appears, wrapping around the
bond.

7.7.2 https://chem.libretexts.org/@go/page/371872
Figure 7.7.4 : The Rotate Fragment Arrow. (CC-NC-BY-DUKE CHEM)
Hold down both the control and shift keys and click and drag the left mouse button to rotate a fragment of the molecule around the
selected bond. In this way put all of the connecting carbons into the trans-orientation, and position the rings so that they point up

and outward away from each other. Click again on the button to minimize the structure. It should now be planar and look like
this.

Figure 7.7.5 : The Nearly Completed Structure. (CC-NC-BY-DUKE CHEM)


Return to the element section button and choose Nitrogen. Double-click on the carbon atoms of the naphthalene rings that are para-
to the connecting three-carbon chain. The color of the atom should change to blue as the element changes from carbon to nitrogen.
Change the element selection button back to Carbon, and change the hybridization to sp3 . Add ethyl groups to the free

valence of the nitrogen atoms. Minimize the structure using . Check your structure. Is it planar? The final molecule should
look something like this (rotated some to show the planar nature).

7.7.3 https://chem.libretexts.org/@go/page/371872
Figure 7.7.6 : The Completed Structure. (CC-NC-BY-DUKE CHEM)
Notice in this picture that the ethyl groups are positioned axial to the plane of the rings. Save this molecule as dye4.
2. At this point you might see that it take more work to build a larger molecule. It might be prudent, therefore to keep the dye4 file
safe, just in case something happens and you need to recall it. To do so, select Save as from the File menu and save the file as
dye4_MNDO. Select Setup, Calculations from the menu. Choose calculate: Equilibrium Geometry at Ground State in Gas
with Semi-Empirical, MNDO. Set Total charge : Cation, Unpaired Electrons: 0. Click Submit.
3. Once the calculation is finished, inspect the molecule. Did the conformation change? Inspect the Output and record the
equilibrium geometry Energy from the Summary section in your lab notebook.
4. In preparation for the next calculation select Save As from the File menu and name the file dye4-ci. Then select Setup,
Calculations again, and change the parameters to, Calculate: Energy at Ground State in Gas with Semi-Empirical, MNDO. Set
Total charge : Cation, Unpaired Electrons: 0. In the Options box: type CI PRINTLEV=2. Click OK.
*Make sure the above is correctly entered: note Energy instead of Equilibrium Geomentry*
5. Submit the calculation.
6. Inspect the Output when the calculation is finished. From the Output tab, look for the table of Spin Allowed Singlet
Excitations (immediately below the very long list of Slater Determinants) find the wavelength (in nm) and energy (in eV) of the
first transition (it should have a large oscillator strength, OS). Record these values in your lab notebook.
7. Compare your computational results with your experimental and theoretical data generated in the previous "wet lab" dye
experiment. Repeat steps 2-6 on dye4 using AM1 and PM3 models.

7.7.4 https://chem.libretexts.org/@go/page/371872
E
,
k
c
a
l
/
m
o
l
(
f
r
o
m
t
h
e
E
q
u
i
l
D
i
y
b
Energy of transition (eV) λmax (nm)
e
r
4
i
u
m
G
e
o
m
e
t
r
y
c
a
l
c
u
l
a
t
i
o
n
)

M
N
D
O

A
M
1

7.7.5 https://chem.libretexts.org/@go/page/371872
P
M
3

For Comparison: λ max from Expt 1 =

 Things to Think About


A Change in Conformation Leads to a Change in Calculated λ : You might have noticed during Exercise 3 with
max

SO2 that with each change in calculation theory (HF, Semi-Empirical, Density Functional) and with each basis set within a
given theory, the calculated bond lengths, bond angles, and frequencies of vibration changed. In the same way in this
current exercise with the dyes, if you look closely you will see that changing the approximations within the Semi-Empirical
calculation (MNDO, PM3, AM1) each lead to slightly different conformations of the molecule. As the result, the predicted
value of both the Energy and λ changes.
max

Seeing Trends in the Calculated Results, and that λ Changes with Conformation is Likely More Important than
max

Seeing if the Calculation Gave the "Right" λ : It would be nice if your calculation could give a calculated λ
max in max

complete agreement with what you saw for the dye when you did the UV-Vis spectrum in Experiment 1. But thinking about
it, there are good reasons why it would not be expected to. For one, your calculation is for a single molecule in the gas
phase. In Experiment 1 you measured the spectrum of a large number of molecules in solution in methanol. Indeed, one
would expect the results to be different. Does that mean that the calculation is of little value?
Thinking about This: Explain why the actual UV-Vis spectrum of the dye in Experiment 1 had such broad peaks. (Perhaps
this understanding, coming from the results of these calculations, is more important than seeing if the calculation gave the
"right" λmax .)

8. Look over the results from Dye 4. Which semi-empirical calculation (MNDO, AM1 or PM3) gave the best result when
comparing λ to your Experiment 1 result? Build the other three dyes. Perform that same semi-empirical calculation for each.
max

Determine the energy of the equilibrium geometry, and the energy (in eV) and λ of the spin allowed singlet excitation. Record
max

the results.

Semi-empirical Calculation did you used

Equil. Geom. E, Energy of transition Calc. λmax Exptl. λ max


Name
(kcal/mol) (eV) (nm) from Expt. 1

Dye 1 1,1'-diethyl-2,2'-cyanine

Dye 2 1,1'-diethyl-2,2'-carbocyanine

Dye 3 1,1'-diethyl-2,2'-dicarbocyanine

Dye 4 1,1'-diethyl-4,4'-carbocyanine

9. What can you conclude at the end of this exercise? Do you see the expected trend in λ as a function of the conjugation length
max

of the connecting carbons between the cyanine rings? Does this more sophisticated calculation give more accurate results in terms
of predicting λ max for the series of dyes as compared to the theoretical, model-based values you calculated in the first dye
experiment? Looking only at Dye 4, which you might recall from Experiment 1 did not show the expected increase in λ even max

though the conjugation length was greater than Dye 3, did the semi-empirical calculation do a better job of accurately predicting
the change in λ maxas you move from Dye 3 to Dye 4?

 Note

This section of your lab notebook should include:


A table recording, for each dye:
Equilibrium geometry energy
Energy of the first spin-allowed singlet transition
Wavelength of the first spin-allowed singlet transition
Comparison of computational λ max to experimental λ
max and theoretical λ max from the previous dye experiment

7.7.6 https://chem.libretexts.org/@go/page/371872
Which semi-empirical calculation gave the most accurate result?
Do you see the expected trend in λ maxas a function of the conjugation length of the connecting carbons between the
cyanine rings?
Does this more sophisticated calculation give more accurate results in terms of predicting λ for the series of dyes as
max

compared to the theoretical, model-based values you calculated in the first dye experiment?
Looking only at Dye 4, which you might recall from Experiment 1 did not show the expected increase in λ even though
max

the conjugation length was greater than Dye 3, did the semi-empirical calculation do a better job of accurately predicting
the change in λ max as you move from Dye 3 to Dye 4?
(See your ELN template for details)

7.7: Exercise 4 - Modeling the Absorption Spectrum of Cyanine Dyes is shared under a not declared license and was authored, remixed, and/or
curated by LibreTexts.

7.7.7 https://chem.libretexts.org/@go/page/371872
7.8: Exercise 5 - Modeling an Intramolecular Rearrangement
This exercise is designed to calculate the transition state of the thermal conversion of vinyl alcohol to acetaldehyde using a semi-
empirical calculation method. Along the way it will show you how you can model a chemical reaction (in this case an
intramolecular rearrangement), and using spreadsheet functionality and graphing to create a reaction profile.

1. Start a new file. Build vinyl alcohol. In the Organic tab, use to construct the carbon skeleton. Select the oxygen and
click on one of the free valences on the vinyl skeleton. If needed rotate the CO bond into the cis-configuration (i.e., bringing the

alcohol hydrogen into position over the C=C double bond). Minimize the structure using .

Figure 7.8.1 : Vinyl Alcohol. (CC-NC-BY-DUKE CHEM)


2. Setup the following calculation: Equilibrium Geometry from Ground State in Water, with Semi-Empirical, AM1. Total
charge: Neutral, Unpaired Electrons : 0. Click Submit. Save the file as vinyl alcohol, and let the calculation run
3. When it has finished, go the Summary tab of the Output and record the Energy in kJ/mol in your lab notebook.
4. Select Save As from the File menu and create a copy of this molecule under the name, vinyl alcohol_ts. Select Guess
Transition State from the Build menu.
5. Click on the C=C bond and then the C-O bond. Next, click on the O-H bond and, while holding down the SHIFT key, click on
the O-H hydrogen to be transferred and then on the carbon atom to receive this hydrogen. Another electron-pushing arrow
should appear along with a faint blue dotted line between the O-H hydrogen and the carbon. It should look like this.

Figure 7.8.2 : Building the Transition State Query. (CC-NC-BY-DUKE CHEM)

6. Click on the button at the lower right of the screen. In a few seconds your initial structure will be replaced by a calculated
transition state structure. There should be a name (1_3-H-Shift_ethenol_AM1.ts (Exact)) at the lower right corner of the screen,
and the structure should look like this.

7.8.1 https://chem.libretexts.org/@go/page/372409
Figure 7.8.3 : The Guessed Transition State. (CC-NC-BY-DUKE CHEM)
7. Setup the following calculation: Calculate Transition State Geometry with Semi-empirical with AM1 in Water. Charge
should be Neutral; Unpaired Electrons, 0; compute: IR. Submit the job.
8. When it has finished, examine the geometry of the transition state using the Geometry functions (bond lengths and angles). Go
to Display, Output and record the energy of the transition state structure.

9. Select Spectra under the Display menu. Use the and to display the calculated IR and frequency table. Animate the

vibrational mode corresponding to the “imaginary” (i) frequency. You might click the button and change the Steps to 30
to slow down the animation.

 The Imaginary (i) Vibration


The idea is this: By now I am sure you are familiar with a potential energy well showing how the energy of a bond varies as a
function of internuclear distance. The minimum of this well represents the equilibrium bond length. Additionally, you know
that the tie-lines across this potential well represent vibrational modes of the bond. Now a similar tie line across the potential
energy curve of a reaction profile, near where the maximum in the curve represents the transition state of the molecule. We can
imagine, and calculate a value for, such a vibrational mode associated with this idea. This calculated "imaginary" vibration is
typically denoted with an i, and models the dynamics of the molecule as it traverses over the top of the activation energy
barrier through the transition state structure.

Notice in the animation that the C=C double bond lengthens as the hydrogen moves toward the carbon, and that the carbon
accepting this hydrogen shifts from sp2 toward sp3 geometry. Notice also that the C-O bond shortens with the hydride shift as it
starts to become a C=O bond. And also that the forming carbonyl carbon shows the shift from sp3 to sp2 geometry. That is, this one
calculated imaginary vibration shows the structural changes in the transition state and across the top of the activation energy barrier
in a reaction profile. Stop the animation and close the Spectrum window.
10. Build acetaldehyde as a New file. Set up the same Equilibrium Geometry calculation as you did for vinyl alcohol (Semi-
empirical, AM1, Water). Save the file and Submit the calculation. Record the energy as you did for vinyl alcohol and the
transition state. Close your molecule.

Figure 7.8.4 : Acetaldehyde. (CC-NC-BY-DUKE CHEM)


Creating the Reaction Profile:

7.8.2 https://chem.libretexts.org/@go/page/372409
11. Open your vinyl alcohol molecule made and optimized in steps 4.8.1-4.8.2. From the File menu, select Save-as and save the
file as Reaction Profile. Select Spreadsheet from the Display menu. A spreadsheet appears on the screen with one entry. This
entry is your vinyl alcohol molecule. Double-click on the M001 entry and change the name to Vinyl Alcohol.

Figure 7.8.5 : Adding Molecules to a Spreasheet. (CC-NC-BY-DUKE CHEM)


12. Now open the transition state molecule created and optimized in steps 4.8.4-4.8.7 above. Notice that this molecule opens in its
own window (you can see the tabs at the lower left of the Spartan screen). Use Select All from the File menu to highlight the
molecule (or right click on the molecule and select all). Use Copy from the File menu (or right mouse click – Copy) to copy the
molecule to the Spartan clipboard. Return to the vinyl alcohol/spreadsheet screen by clicking on its tab. Click once on the cell
immediately under the vinyl alcohol entry. Right-click on that same cell and paste the transition state into the spreadsheet.
Change the name to Transition State.
13. Open and paste your acetaldehyde molecule into the third spreadsheet cell. At this point your spreadsheet will look like this.

Figure 7.8.6 : The Three Molecules in the Spreadsheet. (CC-NC-BY-DUKE CHEM)


14. Open the Setup Calculations box for each of the three molecules in the spreadsheet and confirm that the desired calculation is
correct for each. If it is not, correct them. It should be Equilibrium Geometry in water, with semi-empirical, AM1 for vinyl
alcohol and acetaldehyde; and, Transition State Geometry in water with AM1 for transition state molecule. Click OK (not
Submit) to close the Setup Calculations box each time.
15. Click the Add button in the spreadsheet window. Select Energy (kJ/mol) from the list of choices. Click OK. You will notice a
new Energy column appear in the spreadsheet, with Pending next to each molecule.
16. Return to the Setup Calculations box. Click Submit. The calculations will run and when complete the spreadsheet will show
the aqueous energies for all three entries.

7.8.3 https://chem.libretexts.org/@go/page/372409
Figure 7.8.7 : The Completed Calculations. (CC-NC-BY-DUKE CHEM)

17. Select Plots from the Display menu. Click on the button and highlight Energy in the window that appears. Click Create.

A plot appears in a window to the right of the molecule screen. Click and select Smooth from the Curve options. Add a

Title, and click Done. Click , and the plot, showing the reaction profile (Energy vs molecule) will move to the main screen.
You can use the Animation bar, , to step through the profile, showing the structure at each point
in the process. Record images of the reaction profile in your lab notebook. Use the energy values from the spreadsheet to
calculate the activation energy and ΔH for the reaction.

Figure 7.8.8 : The Reaction Profile. (CC-NC-BY-DUKE CHEM)

 Note

This section of your lab notebook should include:


The equilibrium geometries energies of the initial, transition, and final states
Images of all reaction profile plots
Calculation of the activation energy and ΔH of the tautomerization reaction

7.8.4 https://chem.libretexts.org/@go/page/372409
(See your ELN template for details)

7.8: Exercise 5 - Modeling an Intramolecular Rearrangement is shared under a not declared license and was authored, remixed, and/or curated by
LibreTexts.

7.8.5 https://chem.libretexts.org/@go/page/372409
CHAPTER OVERVIEW

8: Independent Project
8.1: Introduction
8.2: Literature
8.3: Computational

8: Independent Project is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

1
8.1: Introduction
This project is only assigned to students in CHEM310L (not CHEM301L). Its time to see how physical chemistry can be
applied to a topic that really interests you.
In 310L, you can choose a topic for your independent project, and choose whether you would like to dive into more computational
work, or explore other applications of physical chemistry in the peer-reviewed literature. You should begin searching the literature
(eg using Web of Science or looking at the Table of Contents for leading Physical Chemistry journals) for physical chemistry topics
that interest you. Alternatively, you can think about simple chemical phenomena or reactions that you are interested in and that you
can model using Spartan. You should prepar a brief written proposal (1 paragraph) on at least two ideas for what you might do for
your independent project and be prepared to discuss it with your instructors.
You will choose a topic, study it using computation or literature, and present your work to your classmates. Please prepare two
proposals for this project and discuss them with your instructors.

Objectives of the Independent Project:


1. To apply what you're learned in class and lab to learn more about computation and/or a cutting edge research application!
2. To teach your classmates about a cool new topic in physical chemistry!
3. To practice communicating technical scientific topics through an oral presentation.

Start here:
Start by seeing what's out there. Find the table of contents for a journal that focuses on physical chemistry or find the most popular
articles for the last 12 months from one of the ACS journals. Links for some of the Journals of Physical Chemistry from ACS are
linked below.
Journal of Physical Chemistry A
Current Issue (https://pubs.acs.org/toc/jpcafh/current)
Top articels (https://pubs.acs.org/action/showMost...nalCode=jpcafh)
Journal of Physical Chemistry B
Current Issue (https://pubs.acs.org/toc/jpcbfk/current)
Top articles (https://pubs.acs.org/action/showMost...nalCode=jpcbfk)
Journal of Physical Chemistry C
Current Issue (https://pubs.acs.org/toc/jpccck/current)
Top articles (https://pubs.acs.org/action/showMost...nalCode=jpccck)
The Journal of Physical Chemistry Letters
Current Issue (https://pubs.acs.org/toc/jpclcd/current)
Top articels (https://pubs.acs.org/action/showMost...nalCode=jpclcd)
With one of the lists, do the following.
1. Spend 10 minutes to scan the titles and abstracts of the articles in the latest issue(s).
a. What subdisciplines or applications of physical chemistry do you notice?
b. Are there any applications that surprise you?
c. Is there anything that piques your interest?
2. Look at the synopsis figure that accompanies each title.
1. Do the figures change your assignment of subdiscipline?
2. Do the figures change your interest in reading the article?
3. Choose one paper that you would be interesting in learning more about and examine the abstract more closely. Decide if you
want to actually read this paper. If not, move on to another, check out the Journal's ASAP articles, or go to a different journal to
find something that is interesting to you. Alternatively, you can check out the list of top articles read for a specific year (eg: Top
articles for 2021 are here: https://axial.acs.org/2022/01/20/mos...articles-2021/)

8.1.1 https://chem.libretexts.org/@go/page/417545
8.1: Introduction is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

8.1.2 https://chem.libretexts.org/@go/page/417545
8.2: Literature
8.2: Literature is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

8.2.1 https://chem.libretexts.org/@go/page/417546
8.3: Computational
8.3: Computational is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

8.3.1 https://chem.libretexts.org/@go/page/417547
Index

1 https://chem.libretexts.org/@go/page/369889
Glossary
Sample Word 1 | Sample Definition 1

1 https://chem.libretexts.org/@go/page/369890
Detailed Licensing
Overview
Title: CHEM310L - Physical Chemistry I Lab Manual
Webpages: 105
Applicable Restrictions: Noncommercial
All licenses found:
Undeclared: 97.1% (102 pages)
CC BY-NC-SA 4.0: 1.9% (2 pages)
CC BY-SA 4.0: 1% (1 page)

By Page
CHEM310L - Physical Chemistry I Lab Manual - 2.4.6: Output Format - Undeclared
Undeclared 2.4.7: The Help Facility - Undeclared
Front Matter - Undeclared 2.4.8: Array Operations - Undeclared
TitlePage - Undeclared 2.4.8.1: Array Addition and Subtraction -
InfoPage - Undeclared Undeclared
Table of Contents - Undeclared 2.4.8.2: Array Multiplication and Division -
Licensing - Undeclared Undeclared
2.4.8.3: Array Powers - Undeclared
1: Orientation to this course - Undeclared
2.4.8.4: Math Functions - Undeclared
1.1: Pre-lab orientation assignment - Undeclared
2.4.9: Graphics - Undeclared
1.2: Introductory Details - Undeclared
1.3: Rubrics - Undeclared 2.4.9.1: Basic Plots - Undeclared
1.4: Statements - CC BY-SA 4.0 2.4.9.2: Line Types - Undeclared
2: Introduction to Matlab - Undeclared 2.4.9.3: Multiple Plots - Undeclared
2.4.9.4: Manual Axis Scaling - Undeclared
2.1: Pre-lab Assignment - Undeclared
2.4.9.5: Overlaying Plots - Undeclared
2.2: Matlab Practice (In Class Activity) - Undeclared
2.4.10: Linear Least Squares Analysis -
2.2.1: Practice using a custom script - Undeclared
Undeclared
2.2.2: FTIR spectrum of HCl - Undeclared
2.4.11: Last Line Editing and Recall - Undeclared
2.2.3: Emission Spectrum of Iodine Vapor -
2.4.12: Saving Data - Undeclared
Undeclared
2.4.13: Deleting Matrices - Undeclared
2.3: How to Access Matlab - Undeclared 2.4.14: Quitting Matlab - Undeclared
2.3.1: Invoking Matlab From the Lab Computers - 2.4.15: Transferring Files Between Your
Undeclared Computer and Your Duke NetID Account -
2.3.2: Install FastX and access software with your Undeclared
own computer - Undeclared 3: The Treatment of Experimental Error - Undeclared
2.3.3: Access from off-campus (The Duke VPN) -
3.1: Pre-lab Assignment - Undeclared
Undeclared
3.2: Learning Objectives, Important Info -
2.3.4: Install Matlab on your own computer -
Undeclared
Undeclared
3.3: Characterizing Experimental Errors - CC BY-NC-
2.4: Using Matlab (Supplemental) - Undeclared SA 4.0
2.4.1: Entering Matricies - Undeclared 3.4: Propagation of Uncertainty - CC BY-NC-SA 4.0
2.4.2: Matrix Computation - Undeclared 3.5: Least Squares Linear Regression - Undeclared
2.4.3: Vectors - Undeclared 3.5.1: Rejection of Discordant Data - Undeclared
2.4.4: Subscripting - Undeclared 3.5.2: Weighted Least Squares Analysis -
2.4.5: Finding, Deleting, Adding or Replacing Undeclared
Data - Undeclared 3.5.3: References - Undeclared

1 https://chem.libretexts.org/@go/page/417348
3.6: Summary - General Guidelines For The 6: Rotation-Vibration Spectrum of HCl AND DCl;
Treatment of Experimental Errors - Undeclared Vibrational Spectrum of SO2 - Undeclared
3.6.1: Error Analysis in a Nutshell - Undeclared 6.1: Prelaboratory Exercise - Undeclared
3.7: Homework Problems - Undeclared 6.2: Transition Intensities - Undeclared
6.3: Part I - Experimental - Undeclared
4: Absorption Spectrum of Conjugated Dyes -
6.4: Part I - Data Analysis - Undeclared
Undeclared
6.5: Part I - References - Undeclared
4.1: Prelaboratory Assignment - Undeclared
6.6: Part II - Experimental - Undeclared
4.2: Introduction - Undeclared
6.7: Part II - Data Analysis - Undeclared
4.3: Experimental (Part I) - Undeclared
6.8: Part II - References - Undeclared
4.4: Hückel Calculations Using Hulis (Part II) -
7: Molecular Electronic Structure Calculations -
Undeclared
Undeclared
4.5: References - Undeclared
4.6: Appendix A - Use of the Agilent 100 Series UV- 7.1: Pre-Lab Assignment - Undeclared
Vis Spectrophotometer - Undeclared 7.2: Introduction to Electronic Structure Calculations
4.7: Appendix B - The Hückel Approximation - - Undeclared
Undeclared 7.3: Accessing Spartan - Undeclared
5: Molecular Spectroscopy of Iodine - Undeclared 7.4: Exercise 1 - Estimation of the IR Spectrum of
HCl and DCl - Undeclared
5.1: Prelaboratory Assignment - Undeclared
7.5: Exercise 2 - Molecular Orbitals of Formaldehyde
5.2: Introduction - Undeclared
- Undeclared
5.3: Potential Energy Curves - Undeclared
7.6: Exercise 3 - Modeling Sulfur Dioxide -
5.4: Part I - Absorption Spectrum of Iodine Vapor -
Comparing the Results of Different Basis Sets and
Experimental - Undeclared
Calculation Models - Undeclared
5.5: Part I - Absorption Spectrum of Iodine Vapor -
7.7: Exercise 4 - Modeling the Absorption Spectrum
Data Analysis - Undeclared
of Cyanine Dyes - Undeclared
5.6: Part I - Absorption Spectrum of Iodine Vapor -
7.8: Exercise 5 - Modeling an Intramolecular
Questions - Undeclared
Rearrangement - Undeclared
5.7: Part II - Emission Spectrum of Iodine Vapor -
Apparatus - Undeclared 8: Independent Project - Undeclared
5.8: Part II - Emission Spectrum of Iodine Vapor - 8.1: Introduction - Undeclared
Experimental - Undeclared 8.2: Literature - Undeclared
5.9: Part II - Emission Spectrum of Iodine Vapor - 8.3: Computational - Undeclared
Data Analysis - Undeclared Back Matter - Undeclared
5.10: Part II - Emission Spectrum of Iodine Vapor - Index - Undeclared
Questions - Undeclared Glossary - Undeclared
5.11: Part II - Emission Spectrum of Iodine Vapor - Detailed Licensing - Undeclared
References - Undeclared

2 https://chem.libretexts.org/@go/page/417348

You might also like