Einstein-DEbye Function and Tarasov Model

The temperature dependence of the solid heat capacity of polymers can be explained with a microscopic model of thermal motion. The equation that links the vibrational spectrum to the heat capacity is the Einstein function for independent oscillators. It describes the atomic group vibrations in the polymer.

E(θ/T) = [(θ/T)2exp(θ/T)] / [exp(θ/T) - 1]2

CE = NkNE E(θ/T)

and

NE = 3NAtoms - N

where NAtoms is the number of atoms in the characteristic mer unit and N is the number of skeletal modes of vibration in the same structural unit. In the case of polyethylene, the structural unit is CH2 so that 3NAtoms = 9 and N = 2.

The distribution of the skeletal vibrations are coupled in a way that they stretch to zero frequency, that is to the acoustical vibrations. In the lower part of the frequency spectrum, the vibrations will couple intermolecularly, so that the wavelengths of the vibrations become larger than the molecular anisotropy of the chain structure. As a result, the detailed molecular arrangement is of little consequence to these frequencies. Therefore, a three-dimensional Debye function applies

D3(θ/T3) = 3 (θ/T3)30θ/T3 {[(θ/T)4exp(θ/T)] / [exp(θ/T) - 1]2} d(θ/T)

and

Cv / 3Nk = D3 (θ/T3)

where θ3 = h ν3/ k is the reference temperature of skeletal mode vibrations normal to the chain axis. The function above describes only the weak oscillations normal the the chain axis that takes place in the three dimensional force field created by the intermolecular van der Waals forces that act mainly normal to the chain axis. The vibrations along the chain axis, i.e. the one-dimensional oscillation of the strong one-dimensional force field created by the chemically bonded polymer groups is described by a Debye function in one dimension:

D1(θ/T1) = (θ/T1)0θ/T1 {[(θ/T)2exp(θ/T)] / [exp(θ/T) - 1]2} d(θ/T)

and

Cv / 3Nk = D1(θ/T1)

where θ/T = h ν1 / k is the reference temperature of skeletal mode vibrations in the chain axis. These vibration are responsible for the nearly linear increase of the heat capacity between 0 and 200 K.

To approximate all skeletal vibrations of linear polymers, one should start out at low frequency with a three-dimensional Debye function and then switch to a one-dimensional Debye function. Such an approach was derived by Tarasov (1950-1965). The skeletal vibration frequencies are, thus, separated into two groups; the intermolecular group between zero and ν3, and the intramolecular group between ν3 and ν1:

Cv(Tarasov) = NR/3·{D11/T) - (θ31)[D13/T) - D33/T)] }

Knowing the group vibration frequencies (the frequenzy spectrum), the two θ temperatures and the number of skeletal vibrators, N, it is possible to calculate Cv by applying the Tarasov and Einstein equation. Specific heat capacities are usually measured at constant pressure. The thermodynamic relation between Cv and Cp is:

Cp = Cv + T V α2 B 

where B is the bulk modulus and α the volumetric coefficient of thermal expansion. The calculation is quite cumbersome and usually done with a computer program. To understand how the method works, go to Heat Capacity of Polystyrene

References
  1. B. Wunderlich, Pure and Appl. Chem., Vol. 67, No. 6, pp. 1019-1026 (1995)
  2. M. Pyda, and B. Wunderlich, Macomolecules, Vol. 32, 2044 - 2050 (1999)
  3. FD. Porter, Group Interaction Modelling of Polymer Properties, Marcel Dekker 1995
  4. M. Pyda and B. Wunderlich, Macromolecules, 32, 2044-2050 (1999)
  • Summary

    Einstein and Debye Function

    The equation that links the vibrational spectrum to the heat capacity is the Einstein and Debye function.

  • The frequenzy spectrum of a polymer can be separated into group vibrations and skeletal vibrations.

  • The skeletal mode contribution to the molar heat capacity can be calculated with the Tarasov equation.

  • Each group vibration in the polymer chain is assumed to be an independent Einstein oscillator.

  • Group vibrations that spread over a wider frequency range are often approximated by a box-distribution.