How do I calculate the force constant, zero-point energy, and the energy level spacings for ""^(12) "C"""^(16)"O"12C16O if tildeomega_e = "2170 cm"^-1˜ωe=2170 cm1?

1 Answer
Dec 31, 2016

I got, for ""^(12) "C"""^(16)"O"12C16O:

  • k = "1856.92 N/m"k=1856.92 N/m
  • E_0 = 1/2hnu_0 = 2.129xx10^(-20) "J"E0=12hν0=2.129×1020J
  • DeltaE = hnu_0 = 4.258xx10^(-20) "J"

We're treating the system as an anharmonic oscillator, which is basically two hard spheres connected by a spring, with a force constant k, a vibrational zero-point energy E_0, and energy level spacings DeltaE:

![https://upload.wikimedia.org/](useruploads.socratic.org)

Of the above terms, zero-point energy and energy level spacings need further explanation:

  • E_0 is the lowest possible energy of the molecule in its ground state, i.e. the energy E_(upsilon) as a function of the vibrational quantum number upsilon, at the energy level upsilon = 0.
  • The energy levels for this simple harmonic oscillator model go as bb(E_(upsilon) = hnu_0(upsilon + 1/2)), and the energy level spacings are assumed even, because there is no dampening of this ideal spring over time.

The fundamental vibrational frequency is related to all of this, and it is the first harmonic frequency (rather than the second, third, etc. harmonics, which are the overtones) at which the molecule vibrates.

It is the frequency nu_0 that you need to use. However, this is not to be confused with tildenu_0, which is in "cm"^(-1), whereas nu_0 is in "s"^(-1).

Furthermore, you have actually given tildenu_e, the equilibrium frequency. NIST gives

  • tildeomega_e -= tildenu_e = "2169.81358 cm"^(-1)
  • tildeomega_echi_e -= tildenu_e chi_e = "13.28831 cm"^(-1)

And we should recall that for an anharmonic oscillator,

tildenu_0 = tildenu_e - 2tildenu_echi_e.

First, convert over to tildenu_0:

tildenu_0 = "2170 cm"^(-1) - 2 cdot "13.28831 cm"^(-1)

= "2143.4 cm"^(-1)

Then, you therefore need to convert from "cm"^(-1) to "s"^(-1) to get the units to work out for E_(upsilon):

nu_0 = tildenu_0c,

where c = 2.998xx10^10 "cm/s" is the speed of light in the proper units. If you notice, "cm"^(-1) cdot "cm/s" = "s"^(-1), as required for the units of nu_0.

So, we obtain the fundamental vibrational frequency in the correct units so far as:

color(green)(nu_0) = ("2143.4 cm"^(-1))(2.998xx10^10 "cm/s")

= color(green)(6.426xx10^13 "s"^(-1))

So, we can now calculate the zero-point energy. E_0 is the energy defined at upsilon = 0, so:

color(blue)(E_0) = hnu_0(0 + 1/2) = 1/2hnu_0

= 1/2(6.626xx10^(-34) "J"cdot"s")(6.426xx10^13 "s"^(-1))

= color(blue)(2.129xx10^(-20) "J")

compared to an experimental value of 2.149 xx 10^(-20) "J" (the reference shows it in "cm"^(-1) on pg. 7).

Then, the energy level spacings can be calculated without much effort if we use the definition of DeltaE:

color(blue)(DeltaE) = E_1 - E_0 = [ . . . ] = E_(upsilon+1) - E_(upsilon)

= hnu_0(upsilon + 1 + 1/2) - hnu_0(upsilon + 1/2)

= hnu_0(upsilon + 3/2) - hnu_0(upsilon + 1/2)

= hnu_0[(upsilon + 3/2) - (upsilon + 1/2)]

= hnu_0(cancel(upsilon) + 3/2 - cancel(upsilon) - 1/2)

= hnu_0

= (6.626xx10^(-34) "J"cdot"s")(6.426xx10^13 "s"^(-1))

= color(blue)(4.258xx10^(-20) "J")

And now the last thing we need is the force constant k! Recall that the angular frequency from physics is defined as:

omega = sqrt(k/m),

where omega is the angular frequency and m is the mass.

For this system however, we instead utilize the reduced mass, mu = (m_1m_2)/(m_1 + m_2), and treat the system as one "effective mass". Also, recall that omega/(2pi) = nu_0 = 1/T, where T is the period in "s".

Combine that with the equation nu_0 = stackrel(~)nu_0c to find that:

bb(stackrel(~)nu_0 = 1/(2pic)sqrt(k/mu)),

or

bb(nu_0 = 1/(2pi)sqrt(k/mu)).

First, we should find the reduced mass using the ISOTOPIC molar masses of each element:

""_(8)^(16) "O": "15.994915 g/mol"
""_(6)^(12) "C": "12 g/mol"

Note that in order for the units to work out, these must be in bb("kg/mol")!

mu = (m_(C)m_(O))/(m_C + m_O)

= (("0.012000 kg/mol")("0.015994915 g/mol"))/("0.012000 + 0.015994915 kg/mol")

= "0.006856 kg/mol"

Since we already have the value of nu_0, we can now find the force/spring constant:

2pinu_0 = sqrt(k/mu)

=> color(blue)(k) = mu(2pinu_0)^2

= "0.006856 kg"/cancel"mol"xx(cancel"1 mol")/(6.0221413xx10^(23) "things")[(2pi)(6.426xx10^13 "s"^(-1))]^2

= color(blue)("1856.92 kg/s"^2 = "1856.92 N/m")

whereas the experimental force constant in this reference is "1901.82 N/m" (derived from tildeomega_e = 2169.76689(8) "cm"^(-1), which is just another notation for tildenu_e), while in this reference it is "1860 N/m" for tildenu_0 ~~ "2141.42 cm"^(-1).