Analysis of the IR spectrum of carbon monoxide
Comments to Tandy Grubbs - wgrubbs@stetson.edu
Goal:  The infrared (IR) spectrum of carbon monoxide is analyzed in order to determine the rotational constant (B) of this diatomic.

Prerequisites:  A basic knowledge of IR spectroscopy and the quantum mechanical harmonic oscillator and rigid rotor models.

Resources you will need:  This exercise should be carried out within a quantitative analysis software environment that is capable data manipulation and can generate a best-fit line for an x-y data set.  You will also be graphing the data along with the fitted function.


Background:

The infrared (IR) spectrum of a simple heteronuclear diatomic reveals a series of regularly spaced peaks.  A typical spectrum is illustrated in Figure 1.

Figure 1


This exercise will involve an analysis of the IR spectrum of carbon monoxide, which looks much like the spectrum in Figure 1.  To understand the information contained in the CO spectrum, we must take into account the vibrational and rotational energy levels of the molecule.
 
To a first approximation, the vibration of the carbon-oxygen bond is well described by the harmonic oscillator model.  The energy levels of a harmonic oscillator are given by

                                                        (1)

where υ is the vibrational quantum number (0, 1, 2 ….) and ω is the harmonic frequency of oscillator.  This frequency can in turn be calculated from

                                                                        (2)

where k is the spring constant of the bond and μ is the reduced mass of the diatomic,

                                                               (3)

The rotational energy levels of CO can be described using the rigid rotor model, which has energies

                                                            (4)

where J is the rotational quantum number (0, 1, 2, 3 …) and B is the rotational constant,

                                                          (5)

In equation (5), I represents the moment of inertia and R the equilibrium bond distance.
 
The combined rotational-vibrational energy (ro-vibrational energy) can be obtained by adding equations (1) and (4) which yields

                                   (6)

When a diatomic absorbs a photon of IR light, there can be a simultaneous change in vibrational and rotational quantum number.  However, not all transitions are permitted.  The selection rules that apply to this system tell us that only transitions whereby Δυ = ±1 and ΔJ = ±1 are allowed.  At room temperature, nearly all of the CO molecules reside in their ground vibrational state (υ = 0), so the vibrational transition that predominates is the υ = 0 to 1 transition.  But a wide range of rotational states are thermally occupied within the ground vibrational state, so several ΔJ = ±1 rotational transitions can coincide with the υ = 0 to 1 jump. 
 
Figure 2 shows a number of ro-vibrational energy levels for CO.  By convention, ground state parameters are identified with a double prime, “, and excited state parameters with a single prime,’.  Of the transitions that can occur between the υ” and υ’ level, some will involve a decrease in rotational quantum number (ΔJ = -1) and the others will involve an increase (ΔJ = +1).  The peaks in the IR spectrum that arise from ΔJ = -1 transitions are called P-branch transitions and the ΔJ = +1 transitions are called R-branch transitions (three possible P-branch and three possible R-branch transitions are shown in Figure 2).

Figure 2


Referring back to Figure (1), the P-branch transitions give rise to a series of peaks that appear at lower wave numbers.  The R-branch series appears at higher wave numbers.  There appears to be a peak missing between these two series; no peak appears here because ΔJ = 0 transitions are not allowed.
 
In this assignment, you will be given the peak positions of several υ = 0 to 1 transitions in the carbon monoxide IR spectrum with a goal of determining the rotational constant of this diatomic.  To accomplish this, you will first identify pairs of transition peaks that terminate in the same J’ level (one peak each from the P and R branches).  Once we have identified the appropriate pairs, we will subtract their wave number values to obtain an energy difference.  As an example, consider the R-1 and P-2 transitions in Figure 2, which both terminate in the J’=1 level.  By subtracting the wave number values for this pair of transitions we obtain the energy difference between the J”=0 and J”=2 levels (in the ground vibrational state).  This can be represented symbolically by writing

                         (7)

In a similar fashion, the subtraction of successive pairs of peak positions will yield the energy difference between the J”=1 and J”=3 levels, and the J”=2 and J”=4 levels, which can be represented by writing





or more generally as

                                            (8)

Using equation (6), we can derive a theoretical expression for this energy difference as follows:



which simplifies to


                                                    (9)

According to equation (9), a plot of successive energy differences versus (2J"+1) is linear and the slope equals 2B.


Experimental Data:

Click on the following link and save the numerical information to an appropriate folder.  The data corresponds to the fundamental gas phase IR spectrum of 12C16O; the absorption peak wave number positions are listed from smallest to largest in one long column (the P-branch peaks are listed first, the word ‘center’ indicates the relative placement of the missing peak, and then the R-branch peaks are listed).  All wave number values have been corrected for the index of refraction in air (yielding peak positions that would be observed in a vacuum).  The data was obtained from Rao, K.N. and Mathews, C.W., editors, Molecular Spectroscopy: Modern Research, Academic Press, New York, 1972.

                                      Wave number positions for the υ = 0 to 1 infrared transitions in 12C16O
 

Exercise:

1.   Import the data into an appropriate quantitative analysis software package.  For convenience, separate the data into two new columns (or two arrays); one array will contain R-branch peak positions starting with the peak closest to the center of the spectrum (start the array with the peak position originating from J” = 0) and the other array will contain P-branch peak positions starting with the peak closest to the center of the spectrum (start the array with the peak position originating from J” = 1).
 
2.   Which two peaks (one from the P-branch and one from the R-branch) correspond to transitions that terminate in the J’=1 level?  What two peaks terminate in the J’=2 level?  The J’=3 level?
 
3.   Once you recognize the pattern of successive peak pairs that terminate in a particular J’ level, subtract the successive wave number positions so that you generate two new columns of data corresponding to ΔEJ and (2J"+1).  Plot this data along with a best fit line, and determine the rotational constant for this diatomic.  Look up an appropriate literature value for B and calculate the percent difference.
 
4.   Using your result, calculate the moment of inertia (I) and the equilibrium bond length (R) of CO.


5.   MORE CHALLENGING:  The analysis described here yields the rotational constant for the ground vibrational state.  When CO transitions to the υ = 1 state, the equilibrium bond length becomes slightly longer and the rotational constant becomes correspondingly smaller, as dictated by equation (5).  In other words, the rotational constant varies with vibrational quantum number (these rotational constants are distinguished from one another by including the vibrational quantum number as a subscript; B0, B1, B2 …).  How would you alter the analysis described here to obtain B1 (HINT: you will need to find pairs of transitions that originate from the same rotational level in the ground vibrational state)?  Carry out this analysis and report the new rotational constant.  Based on your result, how much does the equilibrium bond length differ between the υ = 0 and υ = 1 levels?
 


Suggestions for improving this web site are welcome.  You are also encouraged to submit your own data-driven exercise to this web archive. All inquiries should be directed to the curator: Tandy Grubbs, Department of Chemistry, Unit 8271, Stetson University, DeLand, FL 32720.

stetsonlogo1