Statistical mechanics

Statistical mechanics is the application of statistics, which includes mathematical tools for dealing with large populations, to the field of Mechanics, which is concerned with the motion of particles or objects when subjected to a force. It provides a framework for relating the microscopic properties of individual atoms and molecules to the macroscopic or bulk properties of materials that can be observed in every day life, therefore explaining thermodynamics as a natural result of statistics and mechanics (classical and quantum). In particular, it can be used to calculate the thermodynamic properties of bulk materials from the spectroscopic data of individual molecules.

At the heart of statistical mechanics is the partition function (see Derivation of the partition function):

<math>Q = \sum_i \exp\left(\frac{-E_i}{kT}\right)</math>

where <math>k</math> is Boltzmann's constant, <math>T</math> is the temperature and <math>E_i</math> reflects each possible energetic state of the system. The partition function provides a measure of the total number of energetic states available to the system at a given temperature. Similarly,

<math>\exp\left(\frac{-E_i}{k T}\right)</math>

provides a measure of the number of energetic states of a particular energy that are likely to be occupied at a given temperautre.

Dividing the second equation by the first equation gives the probability of finding the system in a particular energetic state, <math>i</math>:

<math>p_i = \frac{\exp(-E_i / kT)}{Q}</math>

This probability can be used to find the average value, which corresponds to the macroscopic value, of any property, <math>J</math>, that depends on the energetic state of the system by using the formula:

<math> \langle J \rangle = \sum_i p_i J_i = \sum_i J_i \frac{\exp(-E_i / kT)}{Q}</math>

where <math><J></math> is the average value of property <math>J</math>. This equation can be applied to the internal energy, <math>U</math>, and pressure, <math>P</math>:

<math>U = \sum_i E_i \frac{\exp(-E_i / kT)}{Q}</math>

<math>P = \sum_i P_i \frac{\exp(-E_i / kT)}{Q}</math>

Subsequently, these equations can be combined with known thermodynamic relationships between <math>U</math> and <math>P</math> to arrive at an expression for <math>P</math> in terms of only temperature, volume and the partition function. Similar relationships in terms of the partition function can be derived for other thermodynamic properties as shown in the following table.

Helmholtz free energy:	<math>A = -kT \ln Q</math>
Internal Energy:	<math>U = kT^2 \left( \frac{d\ln Q}{dT} \right)_{N,V}</math>
Pressure:	<math>P = kT \left( \frac{d \ln Q}{dV} \right)_{N,T}</math>
Entropy:	<math>S = k \ln Q + U / T</math>
Gibbs free energy:	<math>G = -kT \ln Q + kTV \left( \frac{d \ln Q}{dV}\right)_{N,T}</math>
Enthalpy:	<math>H = U + PV</math>
Constant Volume Heat Capacity:	<math>C_V = \left( \frac{dU}{dT} \right)_{N,V}</math>
Constant Pressure Heat Capacity:	<math>C_P = \left( \frac{dH}{dT} \right)_{N,P}</math>
Chemical Potential:	<math>\mu_i = -kT \left( \frac{d \ln Q}{dN_i} \right)_{T,V,N}</math>

It is often useful to consider the energy of a given molecule to be distributed among a number of modes. For example, translational energy refers to that portion of energy associated with the motion of the center of mass of the molecule. Configurational energy refers to that portion of energy associated with the various attractive and repulsive forces between molecules in a system. The other modes are all considered to be internal to each molecule. They include rotational, vibrational, electronic and nuclear modes. If we assume that each mode is independent the total energy can be expressed as the sum of each of the components:

<math>E = E_t + E_c + E_n + E_e + E_r + E_v</math>

Where the subscripts t, c, n, e, r, and v correspond to translational, configurational, nuclear, electronic, rotational and vibrational modes, respectively. The relationship in this equation can be substituted into the very first equation to give:

<math>Q = \sum_i \exp\left(-\frac{E_{ti} + E_{ci} + E_{ni} + E_{ei} + E_{ri} + E_{vi}}{kT}\right)</math>

<math>= \sum_i

\exp\left(\frac{-E_{ti}}{kT}\right) \exp\left(\frac{-E_{ci}}{kT}\right) \exp\left(\frac{-E_{ni}}{kT}\right) \exp\left(\frac{-E_{ei}}{kT}\right) \exp\left(\frac{-E_{ri}}{kT}\right) \exp\left(\frac{-E_{vi}}{kT}\right)</math>

<math>= Q_t Q_c Q_n Q_e Q_r Q_v</math>

Thus a partition function can be defined for each mode. Simple expressions have been derived relating each of the various modes to various measurable molecular properties, such as the characteristic rotational or vibrational frequencies.

Expressions for the various molecular partition functions are shown in the following table.

Nuclear	<math>Q_n = 1 \qquad (T < 10^8 K)</math>
Electronic	<math>Q_e = W_0 \exp(kT D_e + W_1 \exp(-\theta_{e1}/T) + \cdots)</math>
vibrational	<math>Q_v = \prod_j \frac{\exp(-\theta_{vj} / 2T)}{1 - \exp(-\theta_{vj} / T)}</math>
rotational (linear)	<math>Q_r = \frac{T}{\sigma} \theta_r</math>
rotational (non-linear)	<math>Q_r = \sqrt{\frac{\pi}{\sigma} \frac{T^3}{\theta_A \theta_B \theta_C)^{1/2}}}</math>
Translational	<math>Q_t = \frac{(2 \pi mkT)^{3/2}}{h^3}</math>
Configurational (ideal gas)	<math>Q_c = V</math>

These equations can be combined with those in the first table to determine the contribution of a particular energy mode to a thermodynamic property. For example the "rotational pressure" could be determined in this manner. The total pressure could be found by summing the pressure contributions from all of the individual modes, ie:

<math>P = P_t + P_c + P_n + P_e + P_r + P_v</math>