Thermal Behavior of CO2

Carbon dioxide is a triatomic molecule, but because it is linear, it behaves in many respects like a diatomic molecule.  In fact, the similarity between CO2 and other diatomic molecules is so striking that the molecule can be described as a Hooke’s Law harmonic oscillator, a model that generally works well only for diatomic molecules.

It is well-known that the heat capacity ratio,CP /CV, for most small molecule gases at room temperature or below, calculated using classical equipartition theory, does not agree well with the ratios determined experimentally.  For example, the calculated heat capacity ratio of N2 gas at 25˚C is 1.28, whereas the experimentally determined value is 1.39.  Contrast this to a monatomic gas like argon.  The calculated heat capacity ratio of Ar at 25˚C is 1.67, while the experimentally determined value is 1.62, very good agreement.  The general explanation for the difference between calculated and experimental heat capacity ratios in molecular systems is justified by the difference between the classical and quantum mechanical descriptions of the systems.  This is supported by the fact that as the temperature of N2 gas, for example, is increased, the capacity ratio begins to approach the calculated value of 1.28.

To understand this, consider the following:  according to classical equipartition theory, at any temperature molecules absorb thermal energy through molecular collisions and distribute this energy equally into their various molecular degrees of freedom.  In order for this to work, molecules must be able to absorb energy continuously.  The important thermal molecular degrees of freedom for gases are translational degrees, the motion of the molecule from one point in space to another as specified by its center of mass, rotations of atoms about various axes in the molecule, as specified by moments of inertia, and vibrations of atoms in molecules along certain allowed normal coordinate axes.  Additionally, one must consider that translational motion and rotational motion only involve kinetic energy, whereas vibrational motion involves both kinetic energy and potential energy.  Classically, all these modes of motion absorb energy continuously, and, according to the equipartition theory, are assigned an amount of energy equal to 1/2 kT for each degree, where k is the Boltzmann constant (1.381 x 10-23 J/K) and is the absolute temperature.  As a reminder, remember that monatomic gases possess only translational degrees, and solids possess (for the most part) only vibrational degrees.

To assign the position of the center of mass of a molecule, we need three coordinates.  Consequently, there are three translational degrees of freedom.  Assigning a value of 1/2 kT for each gives us a total of 3/2 kT energy for the translational degrees of freedom.  If a molecule is linear, there is no moment of inertia about the bond axis (assumed to be the x-axis).  Therefore, rotational degrees only involve the y– and z-axes, giving 2/2 kT energy for the rotational degrees of freedom (linear molecule).  Of course, for a non-linear molecule there are three rotational degrees of freedom, giving 3/2 kT energy for the rotational degrees of freedom (non-linear molecule).  We generally determine the vibrational degrees by difference.  Originally, equipartition theory was described in terms of a single point mass, giving a total of 3 degrees of freedom.  To extend this to a molecule, thought of as a collection of point masses, one needs to assign an x-, y, and z-coordinate for each atom in the molecule.  Consequently, a molecule has a total of 3N degrees of freedom, where N represents the total number of atoms in the molecule.

Let us apply this to N2 gas to see how all this works.  Nitrogen gas is a linear molecule and should have 3N  = 3(2) = 6 total degrees of freedom.  But 3 of those degrees are translational degrees and 2 of those are rotational degrees, so 6 – 5 =  1 must be a vibrational mode.  Generalizing these rules:

Translation:  3 degrees of freedom

Rotation:  (linear molecule) 2 degrees of freedom; (non-linear molecule) 3 degrees of freedom

Vibration:  3N – 5  (linear);  3– 6  (non-linear)

We now can calculate the heat capacity ratio of N2:

Translational Energy:    3/2 kT

Rotational Energy:  2/2 kT

Vibrational Energy:  (1 mode)  1/2 kT for kinetic energy; 1/2 kT for potential energy

Total Energy:  7/2 kT

To find the heat capacity at constant volume:

CodeCogsEqn-6

If we assume that N2 is an ideal gas, then  CP  =  CV  +  k  =  9/2 k.  The heat capacity ratio then is

γ  =  CP /CV  =  9/7  =  1.28

Application of Quantum Theory

We know that for all practical purposes translational energy states for all gases are so closely spaced that we can consider them as being continuous.  The same can be said for the rotational energy of most gases, particularly at room temperature.  However, this is not the case for vibrational energy states.  For many diatomic gases, such as nitrogen, the first excited vibrational energy state is separated from the ground state by a value that is too large for the thermal energy around room temperature, approximately 2500 J/mol, to excite nitrogen molecules to higher vibrational states.  Consequently, in the case of nitrogen, thermal energy is partitioned only between translational and rotational degrees.

Translational energy:  3/2 kT

Rotational energy:  2/2 kT

Total energy:  5/2 kT

Therefore:       CV  =  5/2 kCP  =  7/2 k

Ratio  =  7/5  =  1.40, which is very close to the experimentally determined value.

Application of Quantum Statistical Mechanics to Carbon Dioxide

As previously described, carbon dioxide can be treated as a Hooke’s Law diatomic molecule, because it is linear. We find that CO2 has four vibrational modes:  a symmetric stretching mode at 1388 cm-1  , an asymmetric stretching mode at 2349 cm-1  , and a degenerate pair of bending modes at 667 cm-1  .  We can apply the Boltzmann equation to determine the fraction of CO2 molecules in the first excited state for each vibrational mode.  Keep in mind that we are dealing with thermal excitation here, not light excitation.  It does not matter in this analysis whether the modes are infrared active or not.  Recognizing that ∆E  =  hc/λ, where h is Planck’s constant, c is the speed of light, and 1/λ is the wave number (in reciprocal meters), we have at 15˚C (288K)

Symmetric Stretch:  fraction in first excited state  =   CodeCogsEqn-14

  =    CodeCogsEqn-15

=  CodeCogsEqn-16   =   9.94 x 10-4

Asymmetric Stretch:  fraction  =  8.22 x 10-6

Each Bending Mode:     fraction  =  3.60 x 10-2

This is a surprising result.  While the asymmetric stretching mode is similar to most diatomic molecules, indicating that thermal energy would not promote CO2 molecules to higher energy states in this mode, this is not the case for the other two modes.  In fact, for the bending modes (3.6 x 2) about 7% of the molecules are already in the first excited state due to thermal collisions, irrespective of the presence of infrared radiation.  To check these numbers for accuracy, let us calculate the heat capacity ratio of CO2 and compare it to the experimental value.  The contributions for translation and rotation are the same as they are in the classical case.  The contributions for vibration are not that straight forward.  Certainly, we can eliminate the contribution from asymmetric stretch.  The contributions from the other modes must be calculated using statistical mechanics.  The appropriate equation to do this is

CodeCogsEqn-13

where R is the gas constant 8.314 J/mol·K.  This will give the heat capacity for a vibrational mode in joules per mole rather than in joules per molecule.  The following contributions are:

Symmetric stretch:    CV  =  0.397 J/mol·K

Bending modes:    CV   =  3.57  J/mol·K

The contributions from translation and rotation (per mole) are 3/2 R + 2/2 R  =  20.79 J/mol·K.  We now add all contributions, keeping in mind that there are two bending modes:

CV(total)  =  20.79  +  0.397  +  3.57  +  3.57  =  28.32 J/mol·K

Assuming CO2 is an idea gas (not a great assumption):    C =  36.64 J/mol·K

Ratio  =  1.294 at 15˚C

The actual published heat capacity ratio of CO2 at 15˚C is 1.304, which is quite close, considering we assumed that the gas was ideal.

Discussion

The band of radiation that is infrared active and important in the greenhouse gas scenario is the bending mode at 667 cm-1 (15 microns).  Symmetric stretch, while able to store a small amount of energy (about 229 J/mol) is not IR active.  The fact that about 7% of the CO2 molecules are already in the first excited state at 15˚C due to molecular collisions will affect their ability to absorb IR radiation.  In fact, even in the absence of a source of IR radiation at 15 microns, CO2 will be radiating a small amount of 15 micron radiation.

If we look at the variation of this effect with respect to temperature, we find that from approximately 250K to 300K the fraction of molecules of CO2 in the first excited vibrational bending mode states changes from about 4.4% to 8.0%, which is relatively small, and most likely would not have a large impact on the molecule’s ability to absorb 15 micron IR radiation.  Nevertheless, this does not justify ignoring it, and it should be included in any calculations dealing with the greenhouse gas behavior of CO2.

Advertisements

1 Comment

Filed under Global Warming, James R. Barrante

One response to “Thermal Behavior of CO2

  1. I love what you guys are usually up too. This sort of clever
    work and reporting! Keep up the terrific works guys I’ve added you guys to my own blogroll.

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s