Basics of a Quartz Crystal Microbalance

Similar documents
Online supplementary information

THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER

MODELLING OF FRICTIONAL SOIL DAMPING IN FINITE ELEMENT ANALYSIS

Fluid-structure coupling analysis and simulation of viscosity effect. on Coriolis mass flowmeter

Convective Heat Transfer

Liquid cooling

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Lecture 6: Resonance II. Announcements

Lecture 23 Date:

Mechanics. Free rotational oscillations. LD Physics Leaflets P Measuring with a hand-held stop-clock. Oscillations Torsion pendulum

The Influence of Landau Damping on Multi Bunch Instabilities

Module: 8 Lecture: 1

Tuning of High-Power Antenna Resonances by Appropriately Reactive Sources

FUNDAMENTALS OF POWER SYSTEMS

EP225 Note No. 5 Mechanical Waves

Acoustic Second Harmonic Generation from Rough Surfaces under Shear Excitation in Liquids

Introduction to Laplace Transform Techniques in Circuit Analysis

Question 1 Equivalent Circuits

The Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor

EE C128 / ME C134 Problem Set 1 Solution (Fall 2010) Wenjie Chen and Jansen Sheng, UC Berkeley

Time [seconds]

No-load And Blocked Rotor Test On An Induction Machine

Bogoliubov Transformation in Classical Mechanics

Supporting Information

EELE 3332 Electromagnetic II Chapter 10

Final Comprehensive Exam Physical Mechanics Friday December 15, Total 100 Points Time to complete the test: 120 minutes

SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005.

BASIC INDUCTION MOTOR CONCEPTS

Lateral vibration of footbridges under crowd-loading: Continuous crowd modeling approach

Sampling and the Discrete Fourier Transform

Thermal Resistance Measurements and Thermal Transient Analysis of Power Chip Slug-Up and Slug-Down Mounted on HDI Substrate

Simulation Study on the Shock Properties of the Double-Degree-of-Freedom Cushioning Packaging System

S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS

Chapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem

Finite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization

A) At each point along the pipe, the volume of fluid passing by is given by dv dt = Av, thus, the two velocities are: v n. + ρgy 1

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays

LTV System Modelling

CHAPTER 3 LITERATURE REVIEW ON LIQUEFACTION ANALYSIS OF GROUND REINFORCEMENT SYSTEM

Estimating floor acceleration in nonlinear multi-story moment-resisting frames

THEORETICAL CONSIDERATIONS AT CYLINDRICAL DRAWING AND FLANGING OUTSIDE OF EDGE ON THE DEFORMATION STATES

At the end of this lesson, the students should be able to understand:

Q.1 to Q.30 carry one mark each

The Measurement of DC Voltage Signal Using the UTI

A Linear Model of Magnetostrictive Actuators for Active Vibration Control

Horizontal Biaxial Loading Tests on Sliding Lead Rubber Bearing System

Mucus Transport in the Human Lung Airways: Effect of Porosity Parameter and Air Velocity

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

Electromechanical Dynamics for Micro Film

Open Access Response of Quartz Crystal Microbalance Loaded with Single-drop Liquid in Gas Phase

III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES

A novel protocol for linearization of the Poisson-Boltzmann equation

ECEN620: Network Theory Broadband Circuit Design Fall 2018

Study on the effect of vent on the electroacoustic absorber

Finite Element Truss Problem

Effects of soil structure interaction on behavior of reinforced concrete structures

A Single Particle Thermal Model for Lithium Ion Batteries

Frequency dependent attenuation and dispersion in patchysaturated

Pulsed Magnet Crimping

THE RATIO OF DISPLACEMENT AMPLIFICATION FACTOR TO FORCE REDUCTION FACTOR

Numerical Simulations of Coriolis Flow Meters for Low Reynolds Number Flows

Unified Design Method for Flexure and Debonding in FRP Retrofitted RC Beams

SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. R 4 := 100 kohm

Determination of Flow Resistance Coefficients Due to Shrubs and Woody Vegetation

ME 375 EXAM #1 Tuesday February 21, 2006

A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: CORRESPONDENCE: ABSTRACT

Lecture Notes II. As the reactor is well-mixed, the outlet stream concentration and temperature are identical with those in the tank.

Comparison of Low Field Electron Transport Properties in Compounds of groups III-V Semiconductors by Solving Boltzmann Equation Using Iteration Model

EE 508 Lecture 16. Filter Transformations. Lowpass to Bandpass Lowpass to Highpass Lowpass to Band-reject

ANALYTICAL BEARING MODEL FOR ANALYSIS OF INNER LOAD DISTRIBUTION AND ESTIMATION OF OPERATIONAL LUBRICATION REGIME

Section Induction motor drives

Midterm Test Nov 10, 2010 Student Number:

Stresses near a plate vertex due to a shear force on one of the edges

UNITS FOR THERMOMECHANICS

Effects of vector attenuation on AVO of offshore reflections

Finite-element simulations of Stoneley guided-wave reflection and scattering at the tips of fluid-filled fractures

Seismic Loads Based on IBC 2015/ASCE 7-10

Nonlinear Single-Particle Dynamics in High Energy Accelerators

PHYS 110B - HW #2 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased

SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)

The machines in the exercise work as follows:

We recommend you cite the published version. The publisher s URL is:

FIG. 1: Droplet behaviour with increasing excitation power.

Macromechanical Analysis of a Lamina

EE Control Systems LECTURE 6

LECTURE 22. Collective effects in multi-particle beams: Parasitic Losses. Longitudinal impedances in accelerators (continued)

Reference:W:\Lib\MathCAD\Default\defaults.mcd

CHAPTER 14 SIGNAL GENERATORS AND WAVEFORM-SHAPING CIRCUITS

Contents. Associated Editors and Contributors...XXIII

HOMEWORK ASSIGNMENT #2

SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD

A Method for Assessing Customer Harmonic Emission Level Based on the Iterative Algorithm for Least Square Estimation *

A Buckling Problem for Graphene Sheets. J. Galagher 1, Y. Milman 2, S. Ryan 3, D. Golovaty 3, P. Wilber 3, and A. Buldum 4

B. A. Ivanov 1 and C. E. Zaspel , USA. PACS: a, Gb, Hk

EE 508 Lecture 16. Filter Transformations. Lowpass to Bandpass Lowpass to Highpass Lowpass to Band-reject

Supplementary Figures

These are practice problems for the final exam. You should attempt all of them, but turn in only the even-numbered problems!

Piezoelectric Constitutive Equations. A.1 Three-Dimensional Form of the Linear Piezoelectric Constitutive Equations. S ij = s E. ijkl T kl + d kij E k

FLEXOELECTRIC SIGNALS ON RINGS IN TRANSVERSE MOTIONS

c = ω k = 1 v = ω k = 1 ²µ

Transcription:

Baic of a Quartz Crytal Microbalance Introduction Thi document provide an introduction to the quartz crytal microbalance (QCM) which i an intrument that allow a uer to monitor mall ma change on an electrode. The reader i directed to the numerou review 1 and book chapter 2 for a more in-depth decription concerning the theory and application of the QCM. A baic undertanding of electrical component and concept i aumed. The two major point of thi document are: Explanation of the Piezoelectric Effect Equivalent Circuit Model Explanation of the Piezoelectric Effect The application of a mechanical train to certain type of material (motly crytal) reult in the generation of an electrical potential acro that material. Converely, the application of a potential to the ame material reult in a mechanical train a deformation. Removal of the potential allow the crytal to retore to it original orientation. The igniter on ga grill are a good example of everyday ue of the piezoelectric effect. Depreing the button caue the pring-loaded hammer to trike a quartz crytal thereby producing a large potential that dicharge acro a gap to a metal wire igniting the ga. Quartz i by far the mot widely utilized material for the development of intrument containing ocillator partly due to hitorical reaon (the firt crytal were harveted naturally) and partly due to it commercial availability (ynthetically grown nowaday). There are many way to cut quartz crytal and each cut ha a different vibrational mode upon application of a potential. The AT-cut ha gained the mot ue in QCM application due to it low temperature coefficient at room temperature. Thi mean that mall change in temperature only reult in mall change in frequency. It ha a vibrational mode of thickne hear deformation a hown below in Figure 1. Quartz Crytal No Applied Potential ++++++++++++++++++++++++++ - - - - - - - - - - - - - - - - - - - - - - - - - - Quartz Crytal Under Applied Potential Figure 1. Graphical Repreentation of Thickne Shear Deformation. The application of an alternating potential (a ine wave in nearly all cae) to the crytal face caue the crytal to ocillate. When the thickne of the crytal (t q ) i twice the acoutical wavelength, a tanding wave can be etablihed where the invere of the frequency of the applied potential i ½ of the period of the tanding wave. Thi frequency i called the

reonant frequency, f 0, and i given by the equation ocillation and the reitance repreent energy diipation during ocillation. μq f0 = / 2tq (1) ρq where μ q i the hear modulu (a ratio of heer tre to hear train), ρ q i the denity, and t q i the crytal thickne. The amount of energy lot during ocillation at thi frequency i at a minimum. The ratio of peak energy tored to energy lot per cycle i referred to a the quality factor, Q and i given by the equation Q f c = (2) Δ f FWHM where f c i the center frequency and Δf FWHM i the full width at half max. Thi full width at half max i alo called the bandwidth. For quartz crytal in air, Q can exceed 100,000 while in olution Q decreae to ~3000. Thi i becaue the crytal ha been damped by the olution. Thi damping increae the amount of energy lot per cycle, decreaing Q a hown in Figure 2 below. 30 Figure 3. Quartz Crytal Microbalance Equivalent Mechanical Model. The QCM mechanical model can be electrically modeled in everal different way. The eaiet model to undertand might be an RLC circuit a hown in Figure 4. Figure 4. RLC Circuit. Amplitude 20 10 0 4.4 4.6 4.8 5 5.2 5.4 5.6 Frequency Here R 1 repreent the energy diipated during ocillation, C 1 repreent the energy tored during ocillation, and L 1 repreent the inertial component related to the diplaced ma. At the reonant frequency, f, the impedance of the circuit i at a minimum and i equal in magnitude to R 1 a hown in Figure 5. Figure 2. Comparion of High Q and Low Q. Equivalent Circuit Model The mechanical model (Figure 3) of an electroacoutical ytem conit of a ma (M), a compliance (C m ), and a reitance (r f ). The compliance repreent energy tored during

Figure 5. Impedance Spectrum for a Serie RLC Circuit. Making electrical contact to a quartz crytal i mot eaily done be addition of an electrode to each face of the crytal. Thee electrode introduce an additional capacitance (C0) in parallel with the erie RLC a hown in Figure 6. Thi circuit i commonly referred to a the Butterworth van Dyke (BvD) model. Figure 7. Impedance Spectrum for the Butterworth van Dyke Equivalent Circuit Model. Mot commercial QCM that rely upon a phae lock ocillator, manually cancel out C 0, and only report the erie reonant frequency, f ince f ~= f 0 and f i obviouly dependent upon L 1. The eqcm 10M report both frequencie, f and f p, and a relative impedance pectrum. Figure 6. Butterworth van Dyke Equivalent Circuit Model. The circuit hown above now ha two reonant frequencie, f and f p, which tand for the erie reonant frequency (a in the original RLC circuit) and the parallel reonant frequency, repectively. The impedance pectrum for the BvD model i hown below with a minimum at f and a maximum at f p. Since thee two reonant frequencie are dependent upon L 1, ma change on the urface of the electrode will reult in frequency change. When a depoited film i thin and rigid, the decreae in frequency can be directly correlated to the increae in ma uing the Sauerbrey 4 equation. 2 0 2 f mn Δ f = (3) ( μ ρ ) 1/ 2 q q where f 0 i the fundamental frequency of the crytal a defined in Equation (1), m i the ma added, n i the harmonic number (e.g. n = 1 for a 5 MHz crytal driven at 5 MHz), and μ q and ρ q are a defined above. Equation (3) can be reduced to Δ f = C m (4) f where C f i the calibration contant. The calibration contant for a 5 MHz At-cut quartz crytal in air i 56.6 Hz cm 2 μg -1.

The majority of electrochemical experiment will correpond to a low-loading cae, allowing you to directly correlate change in frequency with change in ma uing the Sauerbrey equation (4). Once a film ha been depoited and the quartz crytal i immered in a liquid, the BvD model can be modified a hown below to include coupling to the liquid. Knowing the liquid vicoity and denity allow one to calculate 3 the expected frequency decreae upon immering the crytal in that liquid uing the equation Δf = 3/ 2 f a ηlρ L πμqρq 1/ 2 (5) where f a i the frequency of the crytal in air, η L i the vicoity of the liquid, ρ L i the denity of the liquid, μ q i the hear modulu of the quartz crytal, and ρ q i the denity of the quartz crytal. Figure 8. Equivalent Circuit Model for a Quartz Crytal Immered in a Liquid Three new component account for the ma loading of the film, L f, and the liquid loading (baed on η L and ρ L ), L l and R l. Both new ma loading L f and L l have the effect of reducing the frequency a indicated with the black arrow in the chart below. The original BvD model (black curve) how reonance approximately 300 khz higher than the modified BvD model (red curve). Notice that the hape of the reonance curve are identical depite the addition of the film and liquid loading. Zmod / Ω 10 7 10 6 10 5 10 4 10 2 BvD modified BvD Serie4 10 1 5 5.2 5.4 5.6 5.8 6 f / MHz Figure 9. Comparion of BvD and Modified BvD Model. For example, upon immering the crytal in pure water, a frequency decreae of approximately 800 Hz i expected. Equation (4) hold a long a the film i aumed to be thin and rigid. When the film i no longer acoutically thin or i not rigid the BvD model can be modified further a hown below. Quartz Ma Loading (film) η L ρ L η f ρ f Figure 10. One Poible Equivalent Circuit Model for a Quartz Crytal Coated with a Polymer Film and Immered in a Liquid μ f (elaticity) Two new feature have been added Film Loading (η f ρ f ) and Elaticity (μ f ). A vicoelatic polymer will influence the reonant frequencie baed on the film vicoity (η f ), denity (ρ f ) and elaticity (μ f ). If the polymer i rigid or η f ρ f doe not change during the experiment, there will be no change in peak hape and contribution from η f ρ f can be ignored. In cae where η f ρ f doe change during the experiment, the frequency, magnitude and hape of the peak will alo change a hown in Figure 11.

Zmod / Ohm 10 7 10 6 10 5 10 4 10 2 BvD Vicou Polymer Modified BvD Serie4 Another way to viualize if η f ρ f i changing during an experiment i to look at the reduced quality factor, Q R a a function of time. Q R ( f + f p ) / 2 = (6) ( f f ) p 10 1 5 5.2 5.4 5.6 5.8 6 f / MHz Figure 11. Comparion of BvD and Vicou Polymer Modified BvD Model. Change in Q R correpond to changing η f ρ f. Quantifying change in η f ρ f i challenging mathematically and will not be explored further here. The reader i directed to the literature 5 for additional information at thi time. Reference 1. Hillman, A. R. The EQCM: electrogravimetry with a light touch, J. Solid State Electrochem., 15, 2011, 1647-1660. 2. Buttry, D. A. The Quartz Cytal Microbalance a an In Situ Tool for Electrochemitry, in Electrochemical Interface. Modern Technique for In-Situ Interface Characterization, H. D. Abruña, Ed., VCH Publiher, Inc., New York, 1991, 531-66. Hillman, A. R. The Electrochemical Quartz Crytal Microbalance. In Intrumentation and Electroanalytical Chemitry; Bard, A. J., Stratmann, M., Unwin, P. R., Ed.; Encyclopedia of Electrochemitry; Wiley: New York, 2003; Vol. 3, 230-289. 3. Glaford, A. P. M. J. Vac. Sci. Technol. 1978, 15, 1836. Kanazawa, K. K.; Gordon II, J. Anal. Chem. 1985, 57, 1770. Kanazawa, K. K.; Gordon II, J. Analytica Chimica Acta 1985, 175, 99-105. 4. Sauerbrey, G. Z. Phy. 1959, 155, 206. 5. Muramatu, H.; Tamiya, E.; Karube, I. Anal. Chem. 1988, 60, 2142-2146. Hillman, A. R.; Jackon, A.; Martin, S. J. Anal. Chem. 2001, 73, 540-549. Martin; Hillman; Etchnique Baic of A Quartz Crytal Microbalance. Rev. 1.3 12/2/2014 Copyright 1990-2014 Gamry Intrument, Inc.