Physics:Sauerbrey equation

The Sauerbrey equation was developed by the German Günter Sauerbrey in 1959, while working on his doctoral thesis at the Technical University of Berlin, Germany. It is a method for correlating changes in the oscillation frequency of a piezoelectric crystal with the mass deposited on it. He simultaneously developed a method for measuring the characteristic frequency and its changes by using the crystal as the frequency determining component of an oscillator circuit. His method continues to be used as the primary tool in quartz crystal microbalance (QCM) experiments for conversion of frequency to mass and is valid in nearly all applications.

The equation is derived by treating the deposited mass as though it were an extension of the thickness of the underlying quartz. Because of this, the mass to frequency correlation (as determined by Sauerbrey’s equation) is largely independent of electrode geometry. This has the benefit of allowing mass determination without calibration, making the set-up desirable from a cost and time investment standpoint.

The Sauerbrey equation is defined as:

[math]\displaystyle{ \Delta f = -\frac{2f_0^2}{A \sqrt{ \rho_q \mu_q } }\Delta m }[/math]

where:

[math]\displaystyle{ f_0 }[/math] – Resonant frequency of the fundamental mode (Hz)

[math]\displaystyle{ \Delta f }[/math] – normalized frequency change (Hz)

[math]\displaystyle{ \Delta m }[/math] – Mass change (g)

[math]\displaystyle{ A }[/math] – Piezoelectrically active crystal area (Area between electrodes, cm²)

[math]\displaystyle{ \rho_q }[/math] – Density of quartz ([math]\displaystyle{ \rho_q }[/math] = 2.648 g/cm³)

[math]\displaystyle{ \mu_q }[/math] – Shear modulus of quartz for AT-cut crystal ([math]\displaystyle{ \mu_q }[/math] = 2.947x10¹¹ g·cm⁻¹·s⁻²)

The normalized frequency [math]\displaystyle{ \Delta f }[/math] is the nominal frequency shift of that mode divided by its mode number (most software outputs normalized frequency shift by default). Because the film is treated as an extension of thickness, Sauerbrey’s equation only applies to systems in which the following three conditions are met: the deposited mass must be rigid, the deposited mass must be distributed evenly and the frequency change [math]\displaystyle{ \Delta f /f }[/math] < 0.05.

If the change in frequency is greater than 5%, that is, [math]\displaystyle{ \Delta f /f }[/math] > 0.05, the Z-match method must be used to determine the change in mass. The formula for the Z-match method is:

[math]\displaystyle{ \frac{\Delta m}{A}\ = \frac{N_q \rho_q}{\pi Z f_L}\tan^{-1} \left [ Z\tan \left ( \pi \frac{f_U-f_L}{f_U} \right ) \right ] }[/math]

Equation 2 – Z-match method

[math]\displaystyle{ f_L }[/math] – Frequency of loaded crystal (Hz)

[math]\displaystyle{ f_U }[/math] – Frequency of unloaded crystal, i.e. Resonant frequency (Hz)

[math]\displaystyle{ N_q }[/math] – Frequency constant for AT-cut quartz crystal (1.668x10¹³Hz·Å)

[math]\displaystyle{ \Delta m }[/math] – Mass change (g)

[math]\displaystyle{ A }[/math] – Piezoelectrically active crystal area (Area between electrodes, cm²)

[math]\displaystyle{ \rho_q }[/math] – Density of quartz ([math]\displaystyle{ \rho_q }[/math] = 2.648 g/cm³)

[math]\displaystyle{ Z }[/math] – Z-Factor of film material [math]\displaystyle{ = \sqrt{ \left ( \frac{\rho_q\mu_q}{\rho_f\mu_f}\ \right ) } }[/math]

[math]\displaystyle{ \rho_f }[/math] – Density of the film (Varies: units are g/cm³)

[math]\displaystyle{ \mu_q }[/math] – Shear modulus of quartz ([math]\displaystyle{ \mu_q }[/math] = 2.947x10¹¹ g·cm⁻¹·s⁻²)

[math]\displaystyle{ \mu_f }[/math] – Shear modulus of film (Varies: units are g·cm⁻¹·s⁻²)

Limitations

The Sauerbrey equation was developed for oscillation in air and only applies to rigid masses attached to the crystal. It has been shown that quartz crystal microbalance measurements can be performed in liquid, in which case a viscosity related decrease in the resonant frequency will be observed:

[math]\displaystyle{ \Delta f = { -\ f_0^{3/2} ( \eta_l \rho_l / \pi \rho_q \mu_q n )^{1/2} } }[/math]

where [math]\displaystyle{ \rho_l }[/math] is the density of the liquid, [math]\displaystyle{ \eta_l }[/math] is the viscosity of the liquid, and [math]\displaystyle{ n }[/math] is the mode number.