Abstract
Selfsensing techniques for atomic force microscope (AFM) cantilevers have several advantageous characteristics compared to the optical beam deflection method. The possibility of down scaling, parallelization of cantilever arrays and the absence of optical interference associated imaging artifacts have led to an increased research interest in these methods. However, for multifrequency AFM, the optimization of the transducer layout on the cantilever for higher order modes has not been addressed. To fully utilize an integrated piezoelectric transducer, this work alters the layout of the piezoelectric layer to maximize both the deflection of the cantilever and measured piezoelectric charge response for a given mode with respect to the spatial distribution of the strain. On a prototype cantilever design, significant increases in actuator and sensor sensitivities were achieved for the first four modes without any substantial increase in sensor noise. The transduction mechanism is specifically targeted at multifrequency AFM and has the potential to provide higher resolution imaging on higher order modes.
Introduction
The invention of the atomic force microscope (AFM) [1] provided for the observation of the nanoscale like no other tool before it [2]. The technologies derived from research into the AFM have led to developments in nanomachining [3], nanometrology [4], material science [5], semiconductor manufacturing [6,7] and highdensity data storage systems [810].
The AFM uses a sharp probe tip at the free end of a cantilever to interrogate and image the surface of a sample [1113]. When using the AFM in dynamic mode [14], the cantilever is excited at its fundamental modal frequency and the probe lightly taps the surface of the sample. Observed changes in the amplitude, phase or frequency shift of the cantilever’s motion correlate to properties of the sample [15]. When closing a feedback loop around these observables with the zaxis nanopositioner, the controller output is routinely used to map the surface topography of the sample. Recently, the additional excitation and detection with multiple frequencies has led to vast improvements in the nanomechanical characterization of the sample beyond it’s topography [16]. For these multifrequency AFM (MFAFM) methods, higher order modes provide enhanced imaging properties such as higher modal stiffnesses and faster response times. It was shown that these higher modes can be more sensitive to material properties such as elastic modulus and damping coefficients [1719]. Additionally, stiff cantilevers have proven to provide high resolution imaging in ambient and liquid environments using quartz resonators [20,21].
Traditional AFM cantilever instrumentation requires a piezoelectric stack actuator at the base of the cantilever for excitation [3] inevitably adding additional resonances as is visible from the so called forest of peaks [22]. These additional frequency components make cantilever resonance tuning almost impossible in liquids [23,24] and can alter the cantilever response rendering the identification and subsequent analysis of higher modes exceedingly difficult. For this reason, numerous integrated actuation methods such as magnetic [25], photothermal [26], resistive thermal [27], ultrasonic [28] or via a piezoelectric layer [29] have been devised.
In order to measure the cantilever deflection, the optical beam deflection (OBD) method [30] is still the widely used standard. However, the measurement setup for the OBD method has limitations, such that it requires frequent laser alignment, a cantilever with a reflective surface and a certain minimum dimension as dictated by the laser spot size. Further, the method suffers from imaging artifacts due to optical interferences originating from stray light reflected by the sample surface [31,32] and bandwidth limitation of the readout circuit [33]. In contrast, a strainbased deflection measurement offers several advantages including a much more compact measurement setup, potential for scalability to cantilever arrays as well as increased sensitivity for smaller cantilever dimensions [3438].
Among the existing integrated actuation and sensing methods, piezoelectric transduction seems to be the only one capable of simultaneously serving as an actuator and a sensor even with a single active layer [39,40]. A set of cantilever designs exist which integrate a piezoelectric transducer onto the cantilever as part of the microfabrication process [41]. While good imaging performance is achieved when operating these cantilevers at higher order modes, the topology of the piezoelectric layer is designed with no consideration of the modal response of the cantilever.
A number of researchers have investigated shaping the piezoelectric layer to actuate or sense a single mode of beam and plate structures whilst filtering the responses of the other modes [4246]. The design of these transducers, denoted modal sensors/actuators, encompasses the modeling and modal analysis of the structure in order to determine the modal frequencies and deflection mode shapes. Here, the overall deflection is the linear combination of mode shapes which have a fixed spatial distribution. Moreover, there exists a linear mapping from the deflection of the structure to the charge developed on the piezoelectric layer and hence the charge response is related to the spatial distribution of each mode shape. If the mode shapes are orthogonal, an analytical approach can be used to shape the piezoelectric layer. Otherwise optimization is used to shape the piezoelectric layer to minimize response to the undesired modes.
This work formulates the design of the topology of the piezoelectric layer on an AFM cantilever to maximize the actuator gain and sensor sensitivity with respect to the cantilever’s higher order modes. Compared to previous work on modal sensor/actuators [4246], the design specification of the presented work is to enhance the desired modes rather than suppress the undesirable modes. This difference leads to fundamental changes in the design strategy, resulting piezoelectric layer topology, instrumentation, actuator characteristics and sensor characteristics. Indeed, the justification for this altered approach comes from the fundamental reasons for multifrequency AFM which is based on the assumption that additional information is encoded in these higher modes. To enable the optimization of the piezoelectric response to higher order modes, plate theory with finite element analysis is used to determine the spatial distribution and polarity of the transducers response for a given mode shape. Using this result, the piezoelectric layer is split into isolated regions whose individual responses constructively combine in order to maximize the actuator gain/sensor output.
The remainder of the paper is organized as follows. Section ’Modal analysis of the piezoelectric cantilever’ outlines the modeling approach to determine the spatial charge distribution of the piezoelectric transducer as a function of its modal response. In section ’Proposed piezoelectric cantilever designs’, results of this analysis are used to determine the design of the piezoelectric actuator arrangements to maximize the transducer response for the first four modes of a cantilever with a stepped geometry. In section ’Instrumentation of the cantilever’ the working principle and modeling of the instrumentation is presented. The experimentally determined actuator and sensor transfer functions are presented in section ’Experimental Results’, which highlight the actuator gain and sensor sensitivity improvements of the proposed designs. In addition, this section presents and discusses the noise characterization of the sensor. The cantilever designs presented in this work target a single mode each. Section ’Instrumentation for multifrequency AFM’ outlines a method to design and instrument the cantilevers to target multiple modes.
Results and Discussion
Modal analysis of the piezoelectric cantilever
Figure 1 shows the silicon cantilever analyzed in this work. The dimensions in the diagram are stated in Table 1. The benefit of the stepped geometry of the cantilever is that higher modes are more closely spaced compared to rectangular cantilevers [47,48] and higher mode deflections are amplified [49] which benefits higher harmonic/higher mode applications [41].
Table 1: Cantilever dimensions shown in Figure 1.
Parameter  Value 

l_{1}  400 μm 
l_{2}  400 μm 
w_{1}  100 μm 
w_{2}  500 μm 
h  10 μm 
The cantilever is modeled using Mindlin plate theory and a finite element (FE) model is developed to perform modal analysis [50,51]. Modal analysis using the FE model provides a solution to the outofplane deflection w(x,y,t) and the rotations of the normal of the cantilever’s neutral plane around the xaxis and yaxis, θ_{x}(x,y,t) and θ_{y}(x,y,t) respectively. These quantities are shown in Figure 1. Assuming a thin piezoelectric layer, the response of the piezoelectric transducer is proportional to the strain at the surface of the cantilever. The inplane strains at the surface of the cantilever are [51]
The electrodes are uniformly distributed on both sides of the piezoelectric layer to generate electric fields only in the zaxis. Therefore, the electric displacement in the piezoelectric material is [52]
where d_{31}, d_{32} and d_{36} are the piezoelectric coefficients. Assuming the piezoelectric material is poled along the zaxis and is homogeneous, the coefficients d_{31} = d_{32} = d and d_{36} = 0 [52,53]. The charge produced is the integral of the electric displacement, that is
The domain Ω is the area of the piezoelectric layer.
Modal analysis with the FE model evaluates harmonic solutions for θ_{x} and θ_{y} of the form
The spatial functions (x,y) and (x,y) are the mode shapes of cantilever.
In the following analysis, the domain Ω from Equation 5 is restricted to the domain A_{e} of a single rectangular element from the mesh used in the FE model. The modal analysis calculates the rotations and at the four nodes of the element. The mode shapes over an element are
where N(x,y) are the shape functions [50,51]
where the dimensions of the rectangular element are 2a× 2b and the origin is placed at the center of the rectangular element.
By substituting the harmonic solution into Equation 5 the charge response of the piezoelectric transducer over the element is
where
B^{e} represents the response due to the spatial nature of the piezoelectric layer in the xyplane for a given mode. This expression is evaluated using Gaussian quadrature. Since the shape functions are quadratic, the derivatives are linear. This allows the exact integral to be evaluated with Gaussian quadrature at the midpoint of the rectangular element. Evaluating B^{e} for each element over the entire cantilever provides the charge response for a given mode.
Proposed piezoelectric cantilever designs
The aim of the proposed cantilever designs is to increase the actuator gain and sensor sensitivity of the piezoelectric transducer to flexural and torsional modes. First, using the finite element method, modal analysis is performed on the cantilever topology shown in Figure 1 to calculate the mode shapes. For mode 1 to mode 4 (M1–M4), the simulated mode shapes of the cantilever are shown in Figure 2a–d. The modal frequencies are 38.3 kHz, 119 kHz, 176 kHz, and 342 kHz. M1, M2 and M4 are flexural modes while M3 is a torsional mode. The modal analysis provides the mode shapes for the deflection and rotations at the nodes of the FE mesh. Using these values, the quantity B^{e} is calculated for each element in the mesh.
If the sign of B^{e} for two elements are the same, the response over the two elements adds constructively. If of opposite sign, the response over the two elements adds destructively. Based on the sign of B^{e}, the piezoelectric layer is split into two, denoted the positive and negative transducer. By actuating and sensing each separately and combining the responses with opposite polarities, the combination of responses is purely constructive. This improves the actuation and sensing by the transducer to its targeted mode.
This design procedure can be posed as an optimization problem. For the ith finite element in the cantilever mesh, the associated B^{e} is denoted The piezoelectric material on each finite element is observed in either a positive polarity or negative polarity. A design parameter is introduced indicating this polarity. If χ_{i} = 1 the finite element is part of the positive transducer otherwise if χ_{i} = −1 it is connected to the negative transducer. Therefore, maximizing the response of the piezoelectric layer is equivalent to the optimization problem
The solution to this optimization problem is χ_{i} = sign(). Figure 2i–l shows the split piezoelectric arrangement for the first four modes of the cantilever. In flexural modes, a significant θ_{y} is observed while θ_{x} is comparatively small. From Equation 1 and Equation 2 this results in the strain ε_{xx} dominating the charge response while the effect of ε_{yy} is insignificant over most of the cantilever area. Though in a few small areas on the cantilever the opposite occurs. In these locations, particularly at the tip of the cantilever, ε_{yy} dominates the charge response while ε_{xx} becomes insignificant. This causes the presence of the small electrodes seen in Figure 2i–l.
The fabricated cantilever designs are shown in Figure 2m–p. The cantilevers are fabricated using the PiezoMUMPs microfabrication process available from the company MEMSCAP Inc [54]. The device layer is a 10 μm thick layer of singlecrystalsilicon deposited on a (100) oriented wafer. A 0.5 μm layer of AlN and a 1 μm layer of aluminium is deposited on the device layer. A particular limitation of this process in the context of AFM is that it does not allow for the fabrication of tips preventing the demonstration of imaging using these cantilevers.
The material properties of the silicon used in the analysis are an elastic modulus of 169 GPa, density of 2500 kg m^{−3} and Poisson’s ratio of 0.29. To account for inaccuracies in these parameters, the design routine was executed with varying parameters to evaluate their affect on the piezoelectric layer topology. Since the piezoelectric material was not included in the mechanical modeling, the thickness of the silicon layer was also varied. It was found that changes in the mode shape, and thus boundary between the two transducers, was insignificant for densities of 1000–4000 kg m^{−3}, elasticities of 120–280 GPa, for Poisson’s ratio of 0.2–0.4 and for a silicon thickness of 8–15 μm. This results from the relative invariance of the mode shape with material properties. Rather, the mode shape is more strongly associated with the geometric shape of the cantilever.
The four cantilevers in Figure 2m–p are denoted C1, C2, C3 and C4. The piezoelectric layer on each cantilever is designed to optimally actuate and sense M1–M4 respectively. Metal traces forming electrical connections run down the center of the cantilever. This splits some piezoelectric layers in two and these split layers are wirebonded back together.
Instrumentation of the cantilever
Instrumentation design
The microfabrication process used to fabricate the cantilevers requires the two piezoelectric transducers share a common terminal [54]. The common terminal has to be grounded to electrically isolate them from each other. Therefore, the actuation and sensing circuits are applied to a grounded load. Two instrumentation arrangements are examined, both shown in Figure 3. The first denoted the voltage driven arrangement, is a grounded load charge sensor with the voltage across the device controlled for actuation. The second denoted the charge driven arrangement, is a grounded load charge amplifier and the voltage across the transducer provides the sensor output.
In the voltage driven arrangement [55], an opamp controls the voltage across the piezoelectric actuator. The charge which flows from the piezoelectric sensor flows into the capacitor C_{s}. A differential amplifier at the output measures the voltage across C_{s} to provide a measurement of the charge. The resistors R_{s} and R_{p} set the DC biases in the circuit. The FET input opamp (OPA656 from Texas Instruments) is used to prevent loading of the piezoelectric transducer. The component values used are C_{s} = 10 pF, R_{s} = 1 MΩ and R_{p} = 10 MΩ.
In the charge driven arrangement [56], the circuit controls the charge across the fixed capacitor C_{s}. An equal charge flows into the piezoelectric transducer as it is in series with C_{s}. Similar to the voltage driven arrangement, the resistors bias the circuit and the FET opamp prevents loading of the piezoelectric device. The component values which determine cutoff frequency and gain are chosen as C_{s} = 100 pF, R_{s} = 10 MΩ and R_{p} = 1 MΩ. To prevent oscillations in the instrumentation circuit, an integrator is used to control the charge on C_{s}. The integrator’s high gain at low frequencies allows the charge to be accurately controlled and their low gain at high frequencies prevents oscillations. The gain of the integral controller k_{i} is set to maintain V_{i} = V_{s} over a bandwidth that contains the modes of interest.
Instrumentation modeling
An applied voltage to the piezoelectric transducer excites motion in the cantilever. The mapping from voltage to displacement is modeled as a set of second order modes. The transfer function from voltage V to displacement d is [52]
where for the ith mode, ω_{i} is the natural frequency, Q_{i} is the quality factor and α_{i} is the gain.
While under motion, the strain on the piezoelectric transducer induces charge on its electrodes. This effect is modeled as an internal voltage source V_{p} in series with a capacitor C_{p} as shown in Figure 4. The transfer function from the applied voltage V to the piezoelectric voltage V_{p} is
The piezoelectric voltage allows for the electric sensing of the motion of the cantilever. Considering the model in Figure 4, the mapping from the voltage applied to the charge Q generated is [22]
There are two terms in this transfer function. The charge associated with the first term is called the feedthrough charge and the charge associated with the second term is called the motional charge. The feedthrough charge flows due to the capacitive structure of the piezoelectric transducer and the motional charge flows due to the strain and is used to observe the motion of the cantilever.
With the piezoelectric transducer incorporated into the voltage driven circuit, the transfer function of the instrumentation is
The model can be simplified by considering the dynamics in the neighborhood of the ith cantilever mode. First, the transfer function of the feedthrough component is identified by letting G_{vv} = 0. The resistance R_{p} and R_{s} are chosen such that the pole and zero in the feedthrough transfer function are much lower than the modal frequencies of the cantilever. Then by considering only the frequencies in the passband (i.e., for ω > 1/R_{s}C_{s} and ω > 1/R_{p}C_{p}, s = jω), the system G_{va} becomes
When using the piezoelectric transducer for realtime sensing such as during AFM imaging, the feedthrough component has to be estimated and removed from the sensor response to maximize the dynamic range of the sensor. This can be done by using modelbased feedforward compensators, implemented in either analog or using switched capacitor prototyping systems such as a Field Programmable Analog Arrays (FPAAs) [39,40]. Since the cantilevers proposed in this work do not feature tips for AFM imaging, the feedthrough is identified and removed offline to highlight the increase in sensor sensitivity.
The charge driven arrangement is the inverse of the voltage driven arrangement. The inversion maintains the same structure as in Equation 21, however the resulting transfer function shows flipped poles and zeros (compare Figure 7) as well as slightly differing gains, quality factors and resonance frequencies due to the internal feedback nature in Equation 20 [40]. The transfer function in the neighborhood of the cantilever’s ith mode is
Experimental results
Experimental setup
Figure 5 shows the experimental setup to characterize the performances of cantilever designs (C1 to C4). The positive and negative transducers are connected to separate instrumentation circuits which are actuated and sensed in the opposite polarity to constructively combine the two responses. For displacement measurement, a vibrometer (Polytec MSA400) is used to detect the motion of the cantilever. The actuator gains (V_{i}→d) are measured in two locations which are shown in Figure 2m. Measurements are made at location 1 for flexural modes M1, M2 and M4 and location 2 for the torsional mode M3.
The system response (V_{i}→V_{o}) is the combination of a motional and feedthrough component. To observe the motional component feedthrough cancellation is performed offline. A third order transfer function is fitted to the measured frequency response in a small band around the mode of interest. The thirdorder model accounts for a second order mechanical system with a first order feedthrough system in parallel. Identification is performed using the subspace method [57]. The identified feedthrough is subtracted from the measurements to produce the system response with feedthrough cancellation (V_{i}→V_{d}). Due to small phase shifts from the opamp dynamics and unmodeled electrical parasitics, feedthrough cancellation can only be accurately performed in a narrowband in the vicinity of the mode of interest.
To evaluate the effect of the proposed piezoelectric topologies, the magnitude responses from V_{i} to both V_{d} and d are measured for each cantilever and compared to the response of C1. C1 is used as the reference cantilever because it is considered a standard topology with the piezoelectric layer covering the entire cantilever.
Discussion of results
From the magnitude frequency response of C1 (Figure 6a), the frequency of the cantilever’s first four modes are at 44.02 kHz, 133.7 kHz, 186.8 kHz and 402.9 kHz respectively. Fabrication tolerances and the mechanical action of the piezoelectric layer account for the frequency differences compared to the FE model in section ’Proposed piezoelectric cantilever designs’.
In Figure 6a–d (voltage driven) and Figure 7a–d (charge driven), the frequency responses from the input voltage V_{i} to the displacement d are shown. The magnitude of the response at each mode is tabulated in Table 2a and c. Compared to C1 in the voltage driven arrangement, C2 provides a 20.45 dB increase in actuation gain for mode 2. C3 provides a 42.05 dB increase in actuation gain for mode 3. C4 provides an 10.04 dB increase in gain for mode 4. Compared to C1 for the charge driven arrangement C2 provides a 22.04 dB increase in actuation gain for mode 2. C3 provides a 47.01 dB increase in actuation gain for mode 3. C4 provides an 14.56 dB increase in gain for mode 4.
Table 2: Sensitivities of the reference cantilever C1 compared to the cantilevers with optimal transducer topologies.
(a) Voltage driven actuator gain (V_{i}→d).  

mode  Reference gain [μm V^{−1}]  Optimized gain [μm V^{−1}]  Improvement [dB] 
M1  (C1) 16.97  –  – 
M2  (C1) 2.968 × 10^{−1}  (C2) 3.125  20.45 
M3  (C1) 1.522 × 10^{−2}  (C3) 1.927  42.05 
M4  (C1) 2.667 × 10^{−1}  (C4) 8.475 × 10^{−1}  10.04 
(b) Voltage driven magnitude of system responses with feedthrough cancellation (V_{i}→V_{d}).  
mode  Reference gain [V V^{−1}]  Optimized gain [V V^{−1}]  Improvement [dB] 
M1  (C1) 3.265  –  – 
M2  (C1) 1.910 × 10^{−2}  (C2) 1.514  37.98 
M3  –  (C3) 7.247 × 10^{−1}  – 
M4  (C1) 1.471 × 10^{−1}  (C4) 1.061  17.16 
(c) Charge driven actuator gain (V_{i}→d).  
mode  Reference gain [μm V^{−1}]  Optimized gain [μm V^{−1}]  Improvement [dB] 
M1  (C1) 25.08  –  – 
M2  (C1) 5.034 × 10^{−1}  (C2) 6.373  22.04 
M3  (C1) 2.626 × 10^{−2}  (C3) 5.890  47.01 
M4  (C1) 4.999 × 10^{−1}  (C4) 2.673  14.56 
(d) Charge driven magnitude of system responses with feedthrough cancellation (V_{i}→V_{d}).  
mode  Reference gain [V V^{−1}]  Optimized gain [V V^{−1}]  Improvement [dB] 
M1  (C1) 9.385 × 10^{−1}  –  – 
M2  (C1) 5.687 × 10^{−3}  (C2) 7.423 × 10^{−1}  42.31 
M3  –  (C3) 6.251 × 10^{−1}  – 
M4  (C1) 4.128 × 10^{−2}  (C4) 7.910 × 10^{−1}  25.65 
In Figure 6e–h and Figure 7e–h, the frequency responses from input voltage V_{i} to output voltage V_{o} are shown. Here substantial feedthrough is observed which dominates in comparison to the motional response of the system. In a neighborhood around the modes of interest, a third order system is identified. Transfer functions were identified around C1 mode 2, C1 mode 4, C2 mode 2, C3 mode 3 and C4 mode 4. The motional response for C1 mode 3 was unobservable due to the small magnitude of resonance response. The parameters of Equation 21 and Equation 23 for each of these transfer functions are tabulated in Table 3.
Table 3: Parameters of the identified transfer functions around the modes.
(a) Voltage driven system parameters.  

C_{p}/C_{s}  Q_{i}  f_{i} [kHz]  δ_{i}α_{i}  
C1 (M1)  5.433  408.6  43.92  1.395 × 10^{−3} 
C1 (M2)  5.716  316.6  133.5  9.630 × 10^{−6} 
C1 (M4)  5.756  358.9  402.9  7.074 × 10^{−5} 
C2 (M2)  8.698  309.0  130.7  5.545 × 10^{−4} 
C3 (M3)  6.718  608.9  175.9  1.761 × 10^{−4} 
C4 (M4)  6.728  381.0  404.3  4.144 × 10^{−4} 
(b) Charge driven system parameters.  
C_{s}/C_{p}  Q_{i}^{*}  f_{i}^{*} [kHz]  δ_{i}^{*}α_{i}^{*}  
C1 (M1)  1.599  390.3  43.94  1.506 × 10^{−3} 
C1 (M2)  1.602  277.4  133.8  1.260 × 10^{−5} 
C1 (M4)  1.598  371.1  402.9  6.562 × 10^{−5} 
C2 (M2)  4.370  290.3  130.7  5.811 × 10^{−4} 
C3 (M3)  5.613  578.7  176.0  1.869 × 10^{−4} 
C4 (M4)  5.434  386.1  404.4  3.077 × 10^{−4} 
The resulting magnitude responses with the feedthrough cancellation from input voltage V_{i} to sensor output V_{d}, are shown in Figure 6i–l (voltage driven) and Figure 7i–l (charge driven). Around each mode the feedthrough is removed and the motional component in the neighborhood of the modal frequency is revealed. The magnitudes of the motional components at each modal frequency are tabulated in Table 2b,d. For the voltage driven arrangement C2 increases the system response of mode 2 by 37.98 dB. C4 increases the system response to mode 4 by 17.16 dB. For the charge driven arrangement C2 increases the system response of mode 2 by 42.31 dB. C4 increases the system response to mode 4 by 25.65 dB.
Noise discussion
Amplitude modulation AFM always requires an actively driven cantilever and subsequent demodulation using a lockin amplifier. Therefore, the subsequent noise characterization is for the demodulated amplitude signal. A 4thorder lowpass filter with cutoff frequency f_{c} = 1 kHz is employed in the lockin amplifier (Zurich Instruments, HF2LI). The voltage noise density plot of the sensor output is obtained by sampling the demodulated amplitude at f_{s} = 57.6 kHz and calculating a power spectral density estimate using Welch’s segment averaging estimator with 64 segments. The results are presented in Figure 8.
The RMS voltage noise present in the sensor output voltage is obtained by integrating the noise density (Figure 8) from 0 to f_{s}/2. The deflection noise is obtained by dividing the voltage noise by the identified sensor sensitivity. The sensor sensitivities are obtained by dividing the system response (Table 2b,d) by the actuator gain (Table 2a,c) for each mode. The voltage noise measurements are made before the feedthrough cancellation, therefore the deflection noise results are for the case where the feedthrough cancellation is noiseless. The sensor sensitivities, sensor output voltage noise and sensor output deflection noise are tabulated in Table 4.
Table 4: Sensor sensitivities and noise performance of nonoptimized and optimized higher order modes.
(a) Voltage driven arrangement.  

C1 M1  C1 M2  C1 M4  C2 M2  C3 M3  C4 M4  
Sensor sensitivity [V/μm]  0.192  0.0644  0.551  0.484  0.376  1.25 
RMS voltage noise [μV]  9.6959  5.4284  6.2378  7.0169  13.9297  8.2834 
RMS deflection noise [pm]  50.4  84.3  11.3  14.5  37.0  6.62 
(b) Charge driven arrangement.  
C1 M1  C1 M2  C1 M4  C2 M2  C3 M3  C4 M4  
Sensor sensitivity [V/μm]  0.0374  0.0113  12.11  8.586  9.422  3.379 
RMS voltage noise [μV]  5.7846  2.5007  2.6202  5.4924  5.3145  4.5816 
RMS deflection noise [pm]  154.6  221.4  31.73  47.16  50.07  15.48 
Due to the higher sensor sensitivities, the deflection noise has decreased demonstrating an improved sensor performance resulting from shaping the piezoelectric layer. This result occurs when noise sources from the sensor and instrumentation dominate the noise output rather than thermomechanical noise. Thermomechanical noise is amplified by the sensor sensitivity and since an associated increase in noise was not observed, this noise source must be insignificant compared to the noise from the instrumentation.
Instrumentation for multifrequency AFM
The cantilevers fabricated and characterized in this work each enhance the sensitivity of a single mode. To facilitate multifrequency AFM, it is preferable to have a cantilever optimized to actuate and sense multiple modes and in addition allow them to be transduced simultaneously with maximized sensitivities. Given the linear nature of the cantilever and instrumentation, superposition can be used to achieve this. For example, the cantilever C2 has the potential to optimally actuate and sense modes M1 and M2 simultaneously. In the setup shown in Figure 9, the M1 input voltage V_{i1} is driven at the M1 resonance frequency and is applied to both piezoelectric transducers with the same polarity to optimally actuate the cantilever for M1. The outputs of the two instrumentation circuits are combined with the same polarity to produce the optimized M1 output voltage V_{o1} which maximizes the sensed response for M1. The M2 input voltage V_{i2} is combined with V_{i1}, however, for the second transducer V_{i2} is combined with a negative polarity to optimally actuate M2. The optimal output voltage for M2, V_{o2}, results from taking the difference between the outputs of the two instrumentation circuits.
This principle can be extended for the design of a piezoelectric cantilever used to actuate and sense multiple modes. The design example for the first four modes is shown in Figure 10. By considering the union of cantilevers C1–C4, the resulting cantilever has 10 separate piezoelectric transducers each with its own instrumentation circuit. The M1–M4 input voltages V_{i1}–V_{i4}, each driven at their respective modal resonance frequencies, are applied to the cantilever on each transducer with a positive polarity for black electrodes (refer to Figure 10) and a negative polarity for white electrodes. The outputs of the instrumentation circuits would be combined in the same fashion for the optimized output voltages V_{o1}–V_{o4}. In comparison to using a single piezoelectric transducer with a sum of sinusoids excitation, the multielectrode design provides increased amplitudes at the expense of more complex instrumentation.
Conclusion
An AFM cantilever with a piezoelectric layer is a versatile transducer for both actuation and displacement sensing. As the response of the active layer is a function of the strain over the surface of the cantilever, careful electrode layout has to be employed. Specifically in multifrequency AFM, the cantilever is excited at higher order modes for which the spatial distribution of the strain cause portions of the piezoelectric layer to contribute destructively to the transducer’s overall response. In this work, we have outlined a design method which provides a systematic way of increasing the actuator gain and sensor sensitivity of the selfsensing piezoelectric cantilever by considering the spatial distribution of the strain for a given mode. The design consists of splitting up the piezoelectric layer into several transducers, whose individual responses are constructively combined. The resulting three terminal piezoelectric device requires additional circuitry for instrumentation compared to a typical two terminal piezoelectric cantilever. For this reason, we have proposed a grounded load charge sensor and a grounded load charge amplifier to realize the selfsensing implementation. The experimental results show that by shaping the electrodes on the piezoelectric layer, significant increases in actuator gain and sensor sensitivity are attained. Furthermore, torsional modes are strongly observed in contrast to a cantilever with an evenly distributed piezoelectric layer. Despite the additional circuitry, no significant increase in noise was observed. Future work will focus on the fabrication of cantilevers with tips with the aim of exploiting the proposed technique to potentially perform higher precision imaging in multifrequency AFM.
References

Binnig, G.; Quate, C. F.; Gerber, C. Phys. Rev. Lett. 1986, 56, 930–933. doi:10.1103/physrevlett.56.930
Return to citation in text: [1] 
Wiesendanger, R. Scanning probe microscopy and spectroscopy; Cambridge University Press: Cambridge, United Kingdom, 1994.
Return to citation in text: [1] 
Bhushan, B. Springer Handbook of Nanotechnology; Springer: Berlin, Germany, 2010.
Return to citation in text: [1] [2] 
Mazzeo, A. D.; Stein, A. J.; Trumper, D. L.; Hocken, R. J. Precis. Eng. 2009, 33, 135–149. doi:10.1016/j.precisioneng.2008.04.007
Return to citation in text: [1] 
Yamanaka, K.; Noguchi, A.; Tsuji, T.; Koike, T.; Goto, T. Surf. Interface Anal. 1999, 27, 600–606. doi:10.1002/(sici)10969918(199905/06)27:5/6<600::aidsia508>3.0.co;2w
Return to citation in text: [1] 
Oliver, R. A. Rep. Prog. Phys. 2008, 71, 076501. doi:10.1088/00344885/71/7/076501
Return to citation in text: [1] 
Yao, T.F.; Duenner, A.; Cullinan, M. Precis. Eng. 2017, 47, 147. doi:10.1016/j.precisioneng.2016.07.016
Return to citation in text: [1] 
Mamin, H. J.; Rugar, D. Appl. Phys. Lett. 1992, 61, 1003–1005. doi:10.1063/1.108460
Return to citation in text: [1] 
Minne, S. C.; Yaralioglu, G.; Manalis, S. R.; Adams, J. D.; Zesch, J.; Atalar, A.; Quate, C. F. Appl. Phys. Lett. 1998, 72, 2340–2342. doi:10.1063/1.121353
Return to citation in text: [1] 
Vettiger, P.; Cross, G.; Despont, M.; Drechsler, U.; Durig, U.; Gotsmann, B.; Haberle, W.; Lantz, M. A.; Rothuizen, H. E.; Stutz, R.; Binnig, G. K. IEEE Trans. Nanotechnol. 2002, 99, 39–55. doi:10.1109/tnano.2002.1005425
Return to citation in text: [1] 
Albrecht, T. R.; Akamine, S.; Carver, T. E.; Quate, C. F. J. Vac. Sci. Technol., A 1990, 8, 3386–3396. doi:10.1116/1.576520
Return to citation in text: [1] 
Alves, M. A. R.; Takeuti, D. F.; Braga, E. S. Microelectron. J. 2005, 36, 51–54. doi:10.1016/j.mejo.2004.10.004
Return to citation in text: [1] 
Utke, I.; Hoffmann, P.; Melngailis, J. J. Vac. Sci. Technol., B: Microelectron. Nanometer Struct.–Process., Meas., Phenom. 2008, 26, 1197. doi:10.1116/1.2955728
Return to citation in text: [1] 
García, R.; Pérez, R. Surf. Sci. Rep. 2002, 47, 197–301. doi:10.1016/S01675729(02)000778
Return to citation in text: [1] 
Bhushan, B. Scanning Probe Microscopy in Nanoscience and Nanotechnology; Springer: Berlin, Germany, 2010.
Return to citation in text: [1] 
Garcia, R.; Herruzo, E. T. Nat. Nanotechnol. 2012, 7, 217–226. doi:10.1038/nnano.2012.38
Return to citation in text: [1] 
Rodríguez, T. R.; García, R. Appl. Phys. Lett. 2004, 84, 449–451. doi:10.1063/1.1642273
Return to citation in text: [1] 
Martínez, N. F.; Patil, S.; Lozano, J. R.; García, R. Appl. Phys. Lett. 2006, 89, 153115. doi:10.1063/1.2360894
Return to citation in text: [1] 
Herruzo, E. T.; Perrino, A. P.; Garcia, R. Nat. Commun. 2014, 5, 3126. doi:10.1038/ncomms4126
Return to citation in text: [1] 
Giessibl, F. J. Mater. Today 2005, 8, 32–41. doi:10.1016/s13697021(05)008448
Return to citation in text: [1] 
Wutscher, E.; Giessibl, F. J. Rev. Sci. Instrum. 2011, 82, 093703. doi:10.1063/1.3633950
Return to citation in text: [1] 
Ruppert, M. G.; Moheimani, S. O. R. Beilstein J. Nanotechnol. 2016, 7, 284–295. doi:10.3762/bjnano.7.26
Return to citation in text: [1] [2] 
Schäffer, T. E.; Cleveland, J. P.; Ohnesorge, F.; Walters, D. A.; Hansma, P. K. J. Appl. Phys. 1996, 80, 3622–3627. doi:10.1063/1.363308
Return to citation in text: [1] 
Putman, C. A. J.; Van der Werf, K. O.; De Grooth, B. G.; Van Hulst, N. F.; Greve, J. Appl. Phys. Lett. 1994, 64, 2454–2456. doi:10.1063/1.111597
Return to citation in text: [1] 
Han, W.; Lindsay, S. M.; Jing, T. Appl. Phys. Lett. 1996, 69, 4111–4113. doi:10.1063/1.117835
Return to citation in text: [1] 
Umeda, N.; Ishizaki, S.; Uwai, H. J. Vac. Sci. Technol., B: Microelectron. Nanometer Struct.–Process., Meas., Phenom. 1991, 9, 1318–1322. doi:10.1116/1.585187
Return to citation in text: [1] 
Fantner, G. E.; Burns, D. J.; Belcher, A. M.; Rangelow, I. W.; YoucefToumi, K. J. Dyn. Syst., Meas., Control 2009, 131, 061104. doi:10.1115/1.4000378
Return to citation in text: [1] 
Yamanaka, K.; Nakano, S. Jpn. J. Appl. Phys., Part 1 1996, 35, 3787–3792. doi:10.1143/JJAP.35.3787
Return to citation in text: [1] 
Indermühle, P.F.; Schürmann, G.; Racine, G.A.; de Rooij, N. F. Sens. Actuators, A 1997, 60, 186. doi:10.1016/s09244247(96)014409
Return to citation in text: [1] 
Meyer, G.; Amer, N. M. Appl. Phys. Lett. 1988, 53, 1045–1047. doi:10.1063/1.100061
Return to citation in text: [1] 
Kassies, R.; van der Werf, K. O.; Bennink, M. L.; Otto, C. Rev. Sci. Instrum. 2004, 75, 689–693. doi:10.1063/1.1646767
Return to citation in text: [1] 
Eaton, P.; West, P. Atomic force microscopy; Oxford University Press: Oxford, United Kingdom, 2010.
Return to citation in text: [1] 
Nievergelt, A. P.; Adams, J. D.; Odermatt, P. D.; Fantner, G. E. Beilstein J. Nanotechnol. 2014, 5, 2459–2467. doi:10.3762/bjnano.5.255
Return to citation in text: [1] 
Tortonese, M.; Barrett, R. C.; Quate, C. F. Appl. Phys. Lett. 1993, 62, 834–836. doi:10.1063/1.108593
Return to citation in text: [1] 
Itoh, T.; Suga, T. Nanotechnology 1993, 4, 218–224. doi:10.1088/09574484/4/4/007
Return to citation in text: [1] 
Arlett, J. L.; Maloney, J. R.; Gudlewski, B.; Muluneh, M.; Roukes, M. L. Nano Lett. 2006, 6, 1000–1006. doi:10.1021/nl060275y
Return to citation in text: [1] 
Li, M.; Tang, H. X.; Roukes, M. L. Nat. Nanotechnol. 2007, 2, 114–120. doi:10.1038/nnano.2006.208
Return to citation in text: [1] 
Dukic, M.; Adams, J. D.; Fantner, G. E. Sci. Rep. 2015, 5, 16393. doi:10.1038/srep16393
Return to citation in text: [1] 
Ruppert, M. G.; Moheimani, S. O. R. Rev. Sci. Instrum. 2013, 84, 125006. doi:10.1063/1.4841855
Return to citation in text: [1] [2] 
Ruppert, M. G.; Moheimani, S. O. R. Novel Reciprocal SelfSensing Techniques for TappingMode Atomic Force Microscopy. In Proceedings of 19th IFAC World Congress, IFAC World Congress, Cape Town, South Africa; 2014.
Return to citation in text: [1] [2] [3] 
Sebastian, A.; Shamsudhin, N.; Rothuizen, H.; Drechsler, U.; Koelmans, W. W.; Bhaskaran, H.; Quenzer, H. J.; Wagner, B.; Despont, M. Rev. Sci. Instrum. 2012, 83, 096107. doi:10.1063/1.4755749
Return to citation in text: [1] [2] 
Lee, C.K.; Moon, F. C. ASME J. Appl. Mech. 1990, 57, 434–441. doi:10.1115/1.2892008
Return to citation in text: [1] [2] 
Jian, K.; Friswell, M. I. J. Intell. Mater. Syst. Struct. 2007, 18, 939–948. doi:10.1177/1045389x06070589
Return to citation in text: [1] [2] 
Tanaka, N.; Sanada, T. Smart Mater. Struct. 2007, 16, 36–46. doi:10.1088/09641726/16/1/004
Return to citation in text: [1] [2] 
SanchezRojas, J. L.; Hernando, J.; Donoso, A.; Bellido, J. C.; Manzaneque, T.; Ababneh, A.; Seidel, H.; Schmid, U. J. Micromech. Microeng. 2010, 20, 055027. doi:10.1088/09601317/20/5/055027
Return to citation in text: [1] [2] 
Donoso, A.; Bellido, J. C. Struct. Multidisc. Optim. 2009, 38, 347–356. doi:10.1007/s0015800802797
Return to citation in text: [1] [2] 
Sadewasser, S.; Villanueva, G.; Plaza, J. A. Appl. Phys. Lett. 2006, 89, 033106. doi:10.1063/1.2226993
Return to citation in text: [1] 
Sadewasser, S.; Villanueva, G.; Plaza, J. A. Rev. Sci. Instrum. 2006, 77, 073703. doi:10.1063/1.2219738
Return to citation in text: [1] 
SalehiKhojin, A.; Bashash, S.; Jalili, N. J. Micromech. Microeng. 2008, 18, 085008. doi:10.1088/09601317/18/8/085008
Return to citation in text: [1] 
Petyt, M. Introduction to Finite Element Vibration Analysis; Cambridge University Press: Cambridge, United Kingdom, 1990.
Return to citation in text: [1] [2] 
Liu, G. R.; Quek, S. S. Finite Element Method: A Practical Course; Elsevier Science: Amsterdam, Netherlands, 2003.
Return to citation in text: [1] [2] [3] 
Moheimani, S. O.; Fleming, A. J. Piezoelectric Transducers for Vibration Control and Damping; Springer: Berlin, Germany, 2006.
Return to citation in text: [1] [2] [3] 
Jian, K.; Friswell, M. I. Mech. Syst. Signal Process. 2006, 20, 2290–2304. doi:10.1016/j.ymssp.2005.05.010
Return to citation in text: [1] 
Cowen, A.; Hames, G.; Glukh, K.; Hardy, B. PiezoMUMPs Design Handbook, 1st ed.; MEMSCAP Inc., 2014.
Return to citation in text: [1] [2] 
Fleming, A. J.; Moheimani, S. O. R. IEEE Trans. Control Syst. Technol. 2005, 13, 98–112. doi:10.1109/tcst.2004.838547
Return to citation in text: [1] 
Fleming, A. J.; Moheimani, S. O. R. Rev. Sci. Instrum. 2005, 76, 073707. doi:10.1063/1.1938952
Return to citation in text: [1] 
McKelvey, T.; Akcay, H.; Ljung, L. IEEE Trans. Autom. Control 1996, 41, 960–979. doi:10.1109/9.508900
Return to citation in text: [1]
42.  Lee, C.K.; Moon, F. C. ASME J. Appl. Mech. 1990, 57, 434–441. doi:10.1115/1.2892008 
43.  Jian, K.; Friswell, M. I. J. Intell. Mater. Syst. Struct. 2007, 18, 939–948. doi:10.1177/1045389x06070589 
44.  Tanaka, N.; Sanada, T. Smart Mater. Struct. 2007, 16, 36–46. doi:10.1088/09641726/16/1/004 
45.  SanchezRojas, J. L.; Hernando, J.; Donoso, A.; Bellido, J. C.; Manzaneque, T.; Ababneh, A.; Seidel, H.; Schmid, U. J. Micromech. Microeng. 2010, 20, 055027. doi:10.1088/09601317/20/5/055027 
46.  Donoso, A.; Bellido, J. C. Struct. Multidisc. Optim. 2009, 38, 347–356. doi:10.1007/s0015800802797 
42.  Lee, C.K.; Moon, F. C. ASME J. Appl. Mech. 1990, 57, 434–441. doi:10.1115/1.2892008 
43.  Jian, K.; Friswell, M. I. J. Intell. Mater. Syst. Struct. 2007, 18, 939–948. doi:10.1177/1045389x06070589 
44.  Tanaka, N.; Sanada, T. Smart Mater. Struct. 2007, 16, 36–46. doi:10.1088/09641726/16/1/004 
45.  SanchezRojas, J. L.; Hernando, J.; Donoso, A.; Bellido, J. C.; Manzaneque, T.; Ababneh, A.; Seidel, H.; Schmid, U. J. Micromech. Microeng. 2010, 20, 055027. doi:10.1088/09601317/20/5/055027 
46.  Donoso, A.; Bellido, J. C. Struct. Multidisc. Optim. 2009, 38, 347–356. doi:10.1007/s0015800802797 
47.  Sadewasser, S.; Villanueva, G.; Plaza, J. A. Appl. Phys. Lett. 2006, 89, 033106. doi:10.1063/1.2226993 
48.  Sadewasser, S.; Villanueva, G.; Plaza, J. A. Rev. Sci. Instrum. 2006, 77, 073703. doi:10.1063/1.2219738 
1.  Binnig, G.; Quate, C. F.; Gerber, C. Phys. Rev. Lett. 1986, 56, 930–933. doi:10.1103/physrevlett.56.930 
5.  Yamanaka, K.; Noguchi, A.; Tsuji, T.; Koike, T.; Goto, T. Surf. Interface Anal. 1999, 27, 600–606. doi:10.1002/(sici)10969918(199905/06)27:5/6<600::aidsia508>3.0.co;2w 
22.  Ruppert, M. G.; Moheimani, S. O. R. Beilstein J. Nanotechnol. 2016, 7, 284–295. doi:10.3762/bjnano.7.26 
50.  Petyt, M. Introduction to Finite Element Vibration Analysis; Cambridge University Press: Cambridge, United Kingdom, 1990. 
51.  Liu, G. R.; Quek, S. S. Finite Element Method: A Practical Course; Elsevier Science: Amsterdam, Netherlands, 2003. 
4.  Mazzeo, A. D.; Stein, A. J.; Trumper, D. L.; Hocken, R. J. Precis. Eng. 2009, 33, 135–149. doi:10.1016/j.precisioneng.2008.04.007 
23.  Schäffer, T. E.; Cleveland, J. P.; Ohnesorge, F.; Walters, D. A.; Hansma, P. K. J. Appl. Phys. 1996, 80, 3622–3627. doi:10.1063/1.363308 
24.  Putman, C. A. J.; Van der Werf, K. O.; De Grooth, B. G.; Van Hulst, N. F.; Greve, J. Appl. Phys. Lett. 1994, 64, 2454–2456. doi:10.1063/1.111597 
54.  Cowen, A.; Hames, G.; Glukh, K.; Hardy, B. PiezoMUMPs Design Handbook, 1st ed.; MEMSCAP Inc., 2014. 
3.  Bhushan, B. Springer Handbook of Nanotechnology; Springer: Berlin, Germany, 2010. 
20.  Giessibl, F. J. Mater. Today 2005, 8, 32–41. doi:10.1016/s13697021(05)008448 
21.  Wutscher, E.; Giessibl, F. J. Rev. Sci. Instrum. 2011, 82, 093703. doi:10.1063/1.3633950 
52.  Moheimani, S. O.; Fleming, A. J. Piezoelectric Transducers for Vibration Control and Damping; Springer: Berlin, Germany, 2006. 
2.  Wiesendanger, R. Scanning probe microscopy and spectroscopy; Cambridge University Press: Cambridge, United Kingdom, 1994. 
3.  Bhushan, B. Springer Handbook of Nanotechnology; Springer: Berlin, Germany, 2010. 
52.  Moheimani, S. O.; Fleming, A. J. Piezoelectric Transducers for Vibration Control and Damping; Springer: Berlin, Germany, 2006. 
53.  Jian, K.; Friswell, M. I. Mech. Syst. Signal Process. 2006, 20, 2290–2304. doi:10.1016/j.ymssp.2005.05.010 
14.  García, R.; Pérez, R. Surf. Sci. Rep. 2002, 47, 197–301. doi:10.1016/S01675729(02)000778 
16.  Garcia, R.; Herruzo, E. T. Nat. Nanotechnol. 2012, 7, 217–226. doi:10.1038/nnano.2012.38 
50.  Petyt, M. Introduction to Finite Element Vibration Analysis; Cambridge University Press: Cambridge, United Kingdom, 1990. 
51.  Liu, G. R.; Quek, S. S. Finite Element Method: A Practical Course; Elsevier Science: Amsterdam, Netherlands, 2003. 
11.  Albrecht, T. R.; Akamine, S.; Carver, T. E.; Quate, C. F. J. Vac. Sci. Technol., A 1990, 8, 3386–3396. doi:10.1116/1.576520 
12.  Alves, M. A. R.; Takeuti, D. F.; Braga, E. S. Microelectron. J. 2005, 36, 51–54. doi:10.1016/j.mejo.2004.10.004 
13.  Utke, I.; Hoffmann, P.; Melngailis, J. J. Vac. Sci. Technol., B: Microelectron. Nanometer Struct.–Process., Meas., Phenom. 2008, 26, 1197. doi:10.1116/1.2955728 
17.  Rodríguez, T. R.; García, R. Appl. Phys. Lett. 2004, 84, 449–451. doi:10.1063/1.1642273 
18.  Martínez, N. F.; Patil, S.; Lozano, J. R.; García, R. Appl. Phys. Lett. 2006, 89, 153115. doi:10.1063/1.2360894 
19.  Herruzo, E. T.; Perrino, A. P.; Garcia, R. Nat. Commun. 2014, 5, 3126. doi:10.1038/ncomms4126 
51.  Liu, G. R.; Quek, S. S. Finite Element Method: A Practical Course; Elsevier Science: Amsterdam, Netherlands, 2003. 
8.  Mamin, H. J.; Rugar, D. Appl. Phys. Lett. 1992, 61, 1003–1005. doi:10.1063/1.108460 
9.  Minne, S. C.; Yaralioglu, G.; Manalis, S. R.; Adams, J. D.; Zesch, J.; Atalar, A.; Quate, C. F. Appl. Phys. Lett. 1998, 72, 2340–2342. doi:10.1063/1.121353 
10.  Vettiger, P.; Cross, G.; Despont, M.; Drechsler, U.; Durig, U.; Gotsmann, B.; Haberle, W.; Lantz, M. A.; Rothuizen, H. E.; Stutz, R.; Binnig, G. K. IEEE Trans. Nanotechnol. 2002, 99, 39–55. doi:10.1109/tnano.2002.1005425 
49.  SalehiKhojin, A.; Bashash, S.; Jalili, N. J. Micromech. Microeng. 2008, 18, 085008. doi:10.1088/09601317/18/8/085008 
6.  Oliver, R. A. Rep. Prog. Phys. 2008, 71, 076501. doi:10.1088/00344885/71/7/076501 
7.  Yao, T.F.; Duenner, A.; Cullinan, M. Precis. Eng. 2017, 47, 147. doi:10.1016/j.precisioneng.2016.07.016 
15.  Bhushan, B. Scanning Probe Microscopy in Nanoscience and Nanotechnology; Springer: Berlin, Germany, 2010. 
41.  Sebastian, A.; Shamsudhin, N.; Rothuizen, H.; Drechsler, U.; Koelmans, W. W.; Bhaskaran, H.; Quenzer, H. J.; Wagner, B.; Despont, M. Rev. Sci. Instrum. 2012, 83, 096107. doi:10.1063/1.4755749 
27.  Fantner, G. E.; Burns, D. J.; Belcher, A. M.; Rangelow, I. W.; YoucefToumi, K. J. Dyn. Syst., Meas., Control 2009, 131, 061104. doi:10.1115/1.4000378 
25.  Han, W.; Lindsay, S. M.; Jing, T. Appl. Phys. Lett. 1996, 69, 4111–4113. doi:10.1063/1.117835 
54.  Cowen, A.; Hames, G.; Glukh, K.; Hardy, B. PiezoMUMPs Design Handbook, 1st ed.; MEMSCAP Inc., 2014. 
26.  Umeda, N.; Ishizaki, S.; Uwai, H. J. Vac. Sci. Technol., B: Microelectron. Nanometer Struct.–Process., Meas., Phenom. 1991, 9, 1318–1322. doi:10.1116/1.585187 
55.  Fleming, A. J.; Moheimani, S. O. R. IEEE Trans. Control Syst. Technol. 2005, 13, 98–112. doi:10.1109/tcst.2004.838547 
56.  Fleming, A. J.; Moheimani, S. O. R. Rev. Sci. Instrum. 2005, 76, 073707. doi:10.1063/1.1938952 
39.  Ruppert, M. G.; Moheimani, S. O. R. Rev. Sci. Instrum. 2013, 84, 125006. doi:10.1063/1.4841855 
40.  Ruppert, M. G.; Moheimani, S. O. R. Novel Reciprocal SelfSensing Techniques for TappingMode Atomic Force Microscopy. In Proceedings of 19th IFAC World Congress, IFAC World Congress, Cape Town, South Africa; 2014. 
41.  Sebastian, A.; Shamsudhin, N.; Rothuizen, H.; Drechsler, U.; Koelmans, W. W.; Bhaskaran, H.; Quenzer, H. J.; Wagner, B.; Despont, M. Rev. Sci. Instrum. 2012, 83, 096107. doi:10.1063/1.4755749 
33.  Nievergelt, A. P.; Adams, J. D.; Odermatt, P. D.; Fantner, G. E. Beilstein J. Nanotechnol. 2014, 5, 2459–2467. doi:10.3762/bjnano.5.255 
57.  McKelvey, T.; Akcay, H.; Ljung, L. IEEE Trans. Autom. Control 1996, 41, 960–979. doi:10.1109/9.508900 
34.  Tortonese, M.; Barrett, R. C.; Quate, C. F. Appl. Phys. Lett. 1993, 62, 834–836. doi:10.1063/1.108593 
35.  Itoh, T.; Suga, T. Nanotechnology 1993, 4, 218–224. doi:10.1088/09574484/4/4/007 
36.  Arlett, J. L.; Maloney, J. R.; Gudlewski, B.; Muluneh, M.; Roukes, M. L. Nano Lett. 2006, 6, 1000–1006. doi:10.1021/nl060275y 
37.  Li, M.; Tang, H. X.; Roukes, M. L. Nat. Nanotechnol. 2007, 2, 114–120. doi:10.1038/nnano.2006.208 
38.  Dukic, M.; Adams, J. D.; Fantner, G. E. Sci. Rep. 2015, 5, 16393. doi:10.1038/srep16393 
30.  Meyer, G.; Amer, N. M. Appl. Phys. Lett. 1988, 53, 1045–1047. doi:10.1063/1.100061 
39.  Ruppert, M. G.; Moheimani, S. O. R. Rev. Sci. Instrum. 2013, 84, 125006. doi:10.1063/1.4841855 
40.  Ruppert, M. G.; Moheimani, S. O. R. Novel Reciprocal SelfSensing Techniques for TappingMode Atomic Force Microscopy. In Proceedings of 19th IFAC World Congress, IFAC World Congress, Cape Town, South Africa; 2014. 
31.  Kassies, R.; van der Werf, K. O.; Bennink, M. L.; Otto, C. Rev. Sci. Instrum. 2004, 75, 689–693. doi:10.1063/1.1646767 
32.  Eaton, P.; West, P. Atomic force microscopy; Oxford University Press: Oxford, United Kingdom, 2010. 
40.  Ruppert, M. G.; Moheimani, S. O. R. Novel Reciprocal SelfSensing Techniques for TappingMode Atomic Force Microscopy. In Proceedings of 19th IFAC World Congress, IFAC World Congress, Cape Town, South Africa; 2014. 
28.  Yamanaka, K.; Nakano, S. Jpn. J. Appl. Phys., Part 1 1996, 35, 3787–3792. doi:10.1143/JJAP.35.3787 
52.  Moheimani, S. O.; Fleming, A. J. Piezoelectric Transducers for Vibration Control and Damping; Springer: Berlin, Germany, 2006. 
29.  Indermühle, P.F.; Schürmann, G.; Racine, G.A.; de Rooij, N. F. Sens. Actuators, A 1997, 60, 186. doi:10.1016/s09244247(96)014409 
22.  Ruppert, M. G.; Moheimani, S. O. R. Beilstein J. Nanotechnol. 2016, 7, 284–295. doi:10.3762/bjnano.7.26 
© 2017 Moore et al.; licensee BeilsteinInstitut.
This is an Open Access article under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (http://www.beilsteinjournals.org/bjnano)