## 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*)^{2}exp(θ/*T*)] / [exp(θ/*T*) - 1]^{2}

*C _{E}* =

*Nk*∑N

_{E}

**E**(θ/

*T*)

and

N_{E} = 3*N _{Atoms}* -

*N*

where *N _{Atoms}* 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 CH

_{2}so that 3

*N*= 9 and

_{Atoms}**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

**D**_{3}(θ/*T*_{3}) = 3 (θ/*T*_{3})^{3}∫_{0}^{θ/T3}
{[(θ/*T*)^{4}exp(θ/*T*)] / [exp(θ/*T*) - 1]^{2}} d(θ/*T*)

and

*C _{v} / *3

*Nk*=

**D**

_{3 }(θ/

*T*

_{3})

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:

**D**_{1}(θ/*T*_{1}) = (θ/*T*_{1})∫_{0}^{θ/T1}
{[(θ/*T*)^{2}exp(θ/*T*)] / [exp(θ/*T*) - 1]^{2}} d(θ/*T*)

and

*C _{v} / *3

*Nk*=

**D**

_{1}(θ/

*T*

_{1})

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}:

*C _{v}*(Tarasov) =

*NR*/3·{

**D**

_{1}(θ

_{1}/

*T*) - (θ

_{3}/θ

_{1})[

**D**

_{1}(θ

_{3}/

*T*) -

**D**

_{3}(θ

_{3}/

*T*)] }

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

*C*and

_{v}*C*is:

_{p} *C _{p}* =

*C*+

_{v}*T V*α

*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

- B. Wunderlich,
*Pure and Appl. Chem*., Vol. 67, No. 6, pp. 1019-1026 (1995) - M. Pyda, and B. Wunderlich,
*Macomolecules*, Vol. 32, 2044 - 2050 (1999) - FD. Porter,
*Group Interaction Modelling of Polymer Properties*, Marcel Dekker, New York (1995)