Vibrational partition function - Calculate from several frequencies

In summary, the conversation discusses the calculation of the vibrational partition function from a list of wavenumbers for molecular systems. The participants discuss the use of different equations to calculate the partition function and the importance of the zero of energy scale. It is concluded that the choice of the zero of energy is arbitrary and does not affect the probability distribution between modes. However, when comparing partition functions, the choice of zero of energy should be consistent to avoid confusion.
  • #1
Giogio
5
0
Hello everybody,

I registered here hoping to finally find a fundated answer about what I by myself seem not be able to figure out.

Question in short:
We have calculated a list of wavenumbers for some molecular systems. How do you get the vibrational partition function from that?

My problem in detail:
Let's take this list of wave numbers as an example:

[cm^-1]
1602,13
3710,78
3829,60
(this is for one water molecule)

Now my colleague explained to me that it's correct to calculate the vibrational partition function as
[itex]q_{vib} = \frac{e^{{-h \nu}/{2kT}}}{1-e^{{-h \nu}/{kT}}}[/itex]
for each frequency and then multiply them all.
Which would give

[itex]q_{vib,1} = 0,020806[/itex]
[itex]q_{vib,2} = 0,000127[/itex]
[itex]q_{vib,3} = 0,000096[/itex]
[itex]q_{vib,tot} = 2,53E-10[/itex]

which seems really weird to me because the partition function becomes lower and lower with every vibrational mode. Can a partition function even be lower than 1?

One of the few things I could find about this is this one, page 5, eq. 3.21 to 3.22: http://www.chem.iitb.ac.in/~bltembe/pdfs/ch_3.pdf
But isn't the zero of energy scale kind of arbitrary?
Anyway, using
[itex]q_{vib} = \frac{1}{1-e^{{-h \nu}/{kT}}}[/itex]
gives

[itex]q_{vib,1} = 1,000433438[/itex]
[itex]q_{vib,2} = 1,000000016[/itex]
[itex]q_{vib,3} = 1,000000009[/itex]
[itex]q_{vib,tot} = 1,000433463[/itex]

which seems way more reasonable to me. (Right? :blushing:)

But it leaves open the question of the zero of energy scale for me. My colleague insists that it's arbitrary where you put it, so both solutions should be correct. But since it gives a factor into each of the single partition functions, how can it be equivalent? Also, would [itex]{h \nu}/{2kT}[/itex] be the correct zero?
Is it even correct to just plainly multiply the single vibrational partition functions?

Thanks for any help. Sorry for the long question, I'm just really confused about this.
Ciao!

Giogio
 
Physics news on Phys.org
  • #2
Both the point of zero energy and the normalization of the partition function are arbitrary, they have no physical meaning.
I have never seen this way to calculate partition functions from spectroscopic lines like that... I would calculate the energy levels as intermediate result first, I think.
 
  • #3
Hey mfb,

Thank you!
How would you do that?
In my example, the energies would be

[itex]\epsilon_1 = h\nu_1 = 3,18E-20 J[/itex]
[itex]\epsilon_2 = h\nu_2 = 7,38E-20 J[/itex]
[itex]\epsilon_3 = h\nu_2 = 7,61E-20 J[/itex]

and the only equation I'm finding for those would be [itex]q_{vib} = \prod_i \frac{e^{{-\epsilon_i}/{2kT}}}{1-e^{{-\epsilon}/{kT}}}[/itex] , which is the same as above.

What exactly do you mean?

Ciao,
Giogio

Edit:
[itex]q_{vib} = \sum_i q_{vib,i}[/itex] ?! Is it that simple? It... would sort of make sense, I'm only slowly beginning to understand what the partition function actually is. But then what's whith that other equation, the product one above?
 
Last edited:
  • #4
Those ϵ are energy differences - but differences between what?
Based on the formulas you use, I think they correspond to the distances between adjacent states, for three different types of vibration. In that case, I would expect a product, where each factor is a sum over all excitations of that specific type of vibration (that sum leads to those e^(...)/(1-e^(...)-terms).
Therefore, I think you can use your formula with the product.

If you use a specific temperature to get numerical values, please write that down.
 
  • #5
Oh, sorry, forgot about the temperature, it was 298K.
Energy differences is a good point, I'll have to think about this some more.
 
  • #6
Hi,

Okay, thanks a lot for your help so far, mfb!

I've now retraced and understood that it is of no importance for the probability distribution between the modes if you use
[itex]q_{vib} = \prod_i \frac{e^{{-\epsilon_i}/{2kT}}}{1-e^{{-\epsilon}/{kT}}}[/itex]

or

[itex]q_{vib} = \prod_i \frac{1}{1-e^{{-\epsilon}/{kT}}}[/itex]

to calculate the partition function, it's just that in the latter case, where the lowest energy state for each mode is set to zero, partition functions close to unity indicate most particles being in the vibrational ground state.

Okay, now I just have one last problem/question:

If you calculate the total partition function:
[itex]q_{tot} = q_{rot}q_{vib}q_{tra}q_{ele}[/itex],

Does it then matter where you set the zero of energy for the vibration?

For instance, for the example from above:

[cm^-1]
1602,13
3710,78
3829,60

I get

[itex]q_{vib} = \prod_i \frac{e^{{-\epsilon_i}/{2kT}}}{1-e^{{-\epsilon}/{kT}}}= 2,60E-10
[/itex]

or

[itex]q_{vib} = \prod_i \frac{1}{1-e^{{-\epsilon}/{kT}}} = 1,000433463
[/itex]

whereas for a system of two water molecules with the wavenumbers

123,05
141,42
144,3
174,76
385,44
633,67
1601,73
1617,01
3559,79
3704,74
3793,53
3818,69
(more, and smaller ones too)

the results are

[itex]q_{vib} = \prod_i \frac{e^{{-\epsilon_i}/{2kT}}}{1-e^{{-\epsilon}/{kT}}} = 4,37E-20[/itex]

or

[itex]q_{vib} = \prod_i \frac{1}{1-e^{{-\epsilon}/{kT}}} = 19,55[/itex]

If I want to compare these two systems with respect to the partition function, I can't imagine it's arbitrary which zero of energy I use, that is, if the vibrational contribution to the total partition function becomes smaller or larger for the bigger system. Is it?

It would be great if you could give me a tip about that as well! :D
 
  • #7
I would always use 1 for the ground state, but again that is an arbitrary choice. A different choice multiplies q by a constant factor, and drops out if you calculate any physical observable (like the probability to find a molecule in a state).
 
  • #8
Ahhh, of course, when comparing the partition functions, they get exponents so that the factor for the zero of energy falls out again. Now it's clear. Thank you very much!
 

1. What is a vibrational partition function?

A vibrational partition function is a mathematical tool used in statistical mechanics to describe the distribution of energy among the different vibrational modes of a molecule. It takes into account the different vibrational frequencies and their associated energies.

2. How do you calculate the vibrational partition function?

The vibrational partition function can be calculated by summing over all the possible energy states of the molecule. This involves multiplying the Boltzmann factor (e^(-E/kT)) with the degeneracy factor (g), which takes into account the number of ways a particular energy state can be achieved. The sum should be taken over all the vibrational frequencies of the molecule.

3. Can the vibrational partition function be calculated from several frequencies?

Yes, the vibrational partition function can be calculated from several frequencies. This is because it takes into account all the vibrational modes of the molecule, which can be described by multiple frequencies.

4. Why is the vibrational partition function important?

The vibrational partition function is important because it allows us to calculate the thermodynamic properties of a molecule, such as its heat capacity and entropy, which are essential for understanding the behavior and stability of a system.

5. How does temperature affect the vibrational partition function?

The vibrational partition function is directly proportional to temperature. As the temperature increases, the vibrational energy levels become more populated, resulting in a higher value for the vibrational partition function. This relationship is described by the Boltzmann distribution.

Similar threads

  • Thermodynamics
Replies
1
Views
666
  • Atomic and Condensed Matter
Replies
2
Views
7K
  • Advanced Physics Homework Help
Replies
1
Views
930
  • Advanced Physics Homework Help
Replies
2
Views
1K
  • Advanced Physics Homework Help
Replies
3
Views
1K
  • Advanced Physics Homework Help
Replies
12
Views
2K
  • Other Physics Topics
Replies
6
Views
1K
  • Advanced Physics Homework Help
Replies
1
Views
1K
  • Advanced Physics Homework Help
Replies
1
Views
1K
  • Advanced Physics Homework Help
Replies
5
Views
2K
Back
Top