QCM Advanced Calculator

Quartz Crystal Microbalance — Multi-Model Frequency Analysis

Sauerbrey Kanazawa–Gordon Voigt–Voinova
← Back to openQCM.com
🔬
Overview
How to use this calculator

This calculator provides real-time analysis for Quartz Crystal Microbalance (QCM) experiments using three complementary physical models. All input fields update results instantly — no button click required. Hover over the ? icons for contextual explanations of each parameter and result.

Choose the appropriate model tab based on your experimental conditions: Sauerbrey for rigid thin films in air, Kanazawa–Gordon for Newtonian liquid loading, and Voigt–Voinova for viscoelastic films with or without liquid contact. All three models share the same AT-cut quartz constants (ρQ = 2.648 g/cm³, μQ = 2.947 × 10¹¹ g/(cm·s²), ZQ ≈ 8.828 × 10⁵ g/(cm²·s)), displayed at the bottom of each panel for reference.

⚖️
1 — Sauerbrey Model
Rigid thin film in air

Use this tab for rigid, acoustically thin films deposited in air (or vacuum). Enter the crystal's resonant frequency f₀, active electrode diameter, observed frequency shift Δf, and film density ρf.

The calculator returns:

  • Mass sensitivity Cf — Sauerbrey constant [ng/(Hz·cm²)]
  • Mass per Hz — total mass per 1 Hz shift over the active area
  • Saturation mass — maximum mass at the 1% Δf/f₀ limit
  • Film thickness — in nm, μm and Å
  • Total mass change Δm and areal mass Δm/A

Validity criterion: |Δf/f₀| < 1/1000 for reliable results; theoretical Sauerbrey limit ≈ 1/100. Beyond this, the acoustic impedance mismatch between film and quartz causes nonlinear deviations.

Switch between QCM Characterization (sensor properties at a given f₀ and electrode geometry) and Film Measurement (thickness and mass from a measured Δf).

💧
2 — Kanazawa–Gordon Model
Newtonian liquid loading

Select this tab when the crystal contacts a semi-infinite Newtonian liquid with no adsorbed film. Two operating modes are available:

  • Forward: given liquid density ρL and viscosity ηL, predict Δf, ΔD, ΔΓ and the viscous penetration depth δ.
  • Inverse: given a measured Δf, extract the ρL·ηL product. If ρL is known, extract ηL individually.

The Kanazawa signature is |Δf| = ΔΓ, equivalently ΔD·f₀/(2·|Δf|) = 1. Significant deviation from this ratio indicates non-Newtonian behaviour, viscoelastic film presence, or interfacial slip.

Penetration depth δ = √(ηL/(π·f₀·ρL)) defines the liquid layer thickness sensed by the shear wave. In water at 10 MHz, δ ≈ 178 nm.

🌊
3 — Voigt–Voinova Model
Viscoelastic film ± liquid

Use this tab for viscoelastic films — hydrogels, polymer brushes, protein adlayers, lipid bilayers, cell monolayers. Provide the film's density ρf, thickness hf, storage modulus G′, and a loss parameter (either film viscosity ηf or loss modulus G″). Choose between Air and Liquid contact medium.

The calculator predicts:

  • Δf — frequency shift (includes both mass and viscoelastic contributions)
  • ΔD — dissipation change (energy lost per oscillation cycle)
  • ΔD/|Δf| ratio — key diagnostic for regime classification
  • |G*| — magnitude of the complex shear modulus
  • tan δ — loss tangent (G″/G′), ratio of dissipated to stored energy
  • ΔfS — Sauerbrey rigid-film prediction for comparison

A regime assessment automatically classifies the film behaviour:

Rigid
Soft-Rigid
Viscoelastic
ΔD/|Δf| < 0.05 0.05 – 0.2 > 0.2

Limiting cases: G′ → ∞ recovers the Sauerbrey equation; hf → 0 with liquid contact recovers Kanazawa. Use this to verify continuity across models.

💡
Practical Tips
For QCM-D and multi-harmonic users
  • All calculations update automatically with a short debounce delay (280 ms) as you type.
  • Negative Δf means mass was added to the crystal; positive Δf means mass was removed or liquid was displaced.
  • For QCM-D users: ΔD = 2ΔΓ / f₀, reported in units of 10⁻⁶.
  • The Voigt panel includes a Sauerbrey rigid-film comparison (ΔfS), so you can quantify how much viscoelasticity shifts the result.
  • Use the loss parameter togglef ↔ G″) in the Voigt panel for convenience. The calculator auto-converts between the two using G″ = ω·ηf at the fundamental frequency.
  • For multi-harmonic analysis, run the Voigt model at each overtone frequency (e.g. 5, 15, 25, 35 MHz for n = 1, 3, 5, 7) and compare Δf/n and ΔDn to diagnose frequency-dependent viscoelastic behaviour.
  • Unit conventions: viscosity in mPa·s (= cP), density in g/cm³, frequency in Hz, moduli in Pa. The calculator uses CGS internally (standard for QCM literature) and converts for display.
📚
Key References
Foundational literature for the three models
  • Sauerbrey: G. Sauerbrey, "Verwendung von Schwingquarzen zur Wägung dünner Schichten und zur Mikrowägung", Z. Phys. 155, 206–222 (1959).
  • Kanazawa–Gordon: K.K. Kanazawa & J.G. Gordon II, "Frequency of a quartz microbalance in contact with liquid", Anal. Chem. 57, 1770–1771 (1985).
  • Voigt–Voinova: M.V. Voinova, M. Rodahl, M. Jonson, B. Kasemo, "Viscoelastic acoustic response of layered polymer films at fluid-solid interfaces: Continuum mechanics approach", Phys. Scr. 59, 391–396 (1999).
  • Comprehensive review: M. Rodahl, F. Höök, B. Kasemo et al., "Quartz crystal microbalance setup for frequency and Q-factor measurements in gaseous and liquid environments", Rev. Sci. Instrum. 66, 3924 (1995).
  • QCM-D technique: M. Rodahl & B. Kasemo, "On the measurement of thin liquid overlayers with the quartz-crystal microbalance", Sens. Actuators A 54, 448–456 (1996).
AT-cut Quartz Rigid Film (Δm) Sauerbrey Model — Rigid Film in Air / Vacuum f₀ → f₀ + Δf Δf < 0 t Air / Vacuum Air / Vacuum
Input Parameters
Valid range: 1 kHz – 1 GHz
Valid range: 0.01 – 10 cm
Negative → mass deposition  |  Positive → mass removal
Valid range: 0.01 – 30 g/cm³
Validity: Rigid, elastically uniform film in air. Keep |Δf/f₀| < 1/1000 for reliable results. Theoretical Sauerbrey limit ≈ 1/100.
Results
Mass Sensitivity Cf ?
Mass per Hz over active area ?
Mass Saturation (at Δf/f₀ = 1%) ?
Maximum Sample Thickness (at saturation) ?
Film Thickness ?
Mass Change Δm ?
Thickness — alternate units
Mass per Unit Area Δm/A ?
Quartz Constants (AT-cut)
Density ρQ2.648 g/cm³
Shear modulus μQ2.947 × 10¹¹ g/(cm·s²)
Acoustic impedance ZQ8.828 × 10⁵ g/(cm²·s)
Sauerbrey Equation
Δm = −Δf · A · √(ρQ·μQ) / (2·f₀²)
Cf = √(ρQ·μQ) / (2·f₀²)   [g/(Hz·cm²)]
t = −Δf · Cf / ρf
Limits: For a 10 MHz crystal, Cᶠ ≈ 4.42 ng/(Hz·cm²). Mass resolution is sub-nanogram per cm². Deviations from linearity (Sauerbrey breakdown) begin when the film acoustic thickness approaches that of the quartz. Monitor ΔD to verify the rigid-film assumption holds.
δ AT-cut Quartz Kanazawa–Gordon Model — Newtonian Liquid Loading Newtonian Liquid (ρL, ηL) |Δf| = ΔΓ Newtonian signature Air
Input Parameters
Water 20°C ≈ 1.002 mPa·s  |  Glycerol ≈ 1412 mPa·s
Kanazawa model: QCM in contact with a semi-infinite Newtonian liquid. Characteristic signature: |Δf| = ΔΓ, i.e. ΔD = 2|Δf|/f₀. Higher ΔD/Δf ratios indicate non-Newtonian behaviour.
Results Forward
Frequency Shift Δf ?
Dissipation Change ΔD (× 10⁻⁶) ?
Half-Bandwidth Change ΔΓ ?
ρL · ηL product ?
Viscous Penetration Depth δ ?
Quartz Constants (AT-cut)
ρQ2.648 g/cm³
μQ2.947 × 10¹¹ g/(cm·s²)
ZQ = √(ρQ·μQ)8.828 × 10⁵ g/(cm²·s)
Kanazawa – Gordon Equation
Δf = −f₀3/2 · √(ρL·ηL / (π·ρQ·μQ))
ΔΓ = |Δf|  ⟹  ΔD = 2|Δf|/f₀
δ = √(ηL / (π·f₀·ρL))
Inverse: ρL·ηL = Δf² · π·ρQ·μQ / f₀³

Validity: Semi-infinite, homogeneous Newtonian liquid, no-slip at interface, no film present. Combined ρ·η characterisation of a liquid from a single Δf measurement.
AT-cut Quartz Viscoelastic Film (G', G'') ρf, hf, ηf G' ηf Quartz: Film: Liquid: amplitude hf Voigt–Voinova Model — Viscoelastic Film ± Liquid Liquid (optional) ΔD/|Δf| → regime G' → ∞ ⟹ Sauerbrey hf → 0 ⟹ Kanazawa Air
Input Parameters

Film Properties
Polymer ~1.0–1.5  |  Protein layer ~1.3–1.4
Soft gel ~10³–10⁴ Pa  |  Stiff polymer ~10⁶–10⁸ Pa
Loss Parameter
Soft gel ~ 0.001–0.1 Pa·s

Contact Medium
Results In Air
Frequency Shift Δf ?
Dissipation Change ΔD (× 10⁻⁶) ?
Half-Bandwidth Change ΔΓ

ΔD / |Δf| ratio (× 10⁻⁶ / Hz) ?
Complex Modulus |G*| ?
Loss Tangent tan δ = G'' / G' ?
Sauerbrey rigid-film limit ΔfS ?
Enter parameters to see viscoelastic regime assessment.
Quartz Constants (AT-cut)
ρQ2.648 g/cm³
μQ2.947 × 10¹¹ g/(cm·s²)
ZQ8.828 × 10⁵ g/(cm²·s)
ReferenceVoinova et al., 1999
Voigt – Voinova Model
Δf̃ = Δf + iΔΓ = −(f₀/πZQ) · Z̃load

load = Z̃f · (Z̃L + Z̃f·tanh ξ̃f) / (Z̃f + Z̃L·tanh ξ̃f)

f = √(ρf·G̃f)    G̃f = G' + iG'' = G' + iωηf
L = √(iω·ρL·ηL)   [0 in air]
ξ̃f = ω·hf·√(ρf/G̃f)
Limits: G'→∞ ⟹ Sauerbrey  |  hf→0 ⟹ Kanazawa
ΔD = 2ΔΓ/f₀

High ΔD/|Δf| ratio indicates strong viscoelastic losses. Typical threshold: ΔD/|Δf| > 0.2×10⁻⁶/Hz signals Sauerbrey breakdown.