DOI:
10.1039/C6AN02674E
(Paper)
Analyst, 2017,
142, 14921498
Simultaneous viscosity and density measurement of small volumes of liquids using a vibrating microcantilever†
Received
17th December 2016
, Accepted 18th March 2017
First published on 20th March 2017
Many industrial and technological applications require precise determination of the viscosity and density of liquids. Such measurements can be time consuming and often require sampling substantial amounts of the liquid. These problems can partly be overcome with the use of microcantilevers but most existing methods depend on the specific geometry and properties of the cantilever, which renders simple, accurate measurement difficult. Here we present a new approach able to simultaneously quantify both the density and the viscosity of microliters of liquids. The method, based solely on the measurement of two characteristic frequencies of an immersed microcantilever, is completely independent of the choice of a cantilever. We derive analytical expressions for the liquid's density and viscosity and validate our approach with several simple liquids and different cantilevers. Application of our model to nonNewtonian fluids shows that the calculated viscosities are remarkably robust when compared to measurements obtained from a standard rheometer. However, the results become increasingly dependent on the cantilever geometry as the frequencydependent nature of the liquid's viscosity becomes more significant.
Introduction
Accurate and rapid determination of the density and viscosity of liquids is central to countless industrial, technological and scientific processes. Applications range from oil and lubricant characterization in the petroleum industry^{1} to chemical engineering,^{2} quality control in food science^{3} and biomedical research, in particular for the detection and diagnosis of diseases from bodily fluids.^{4–6} One of the challenges faced by conventional measurement methods is the need for large volumes of liquid. Standard rheometers can provide accurate viscosity measurements over an extensive range of temperatures and pressures^{7} but they require relatively large samples, typically several millilitres or more. Measurement methods based on acoustic waves,^{8} tuning forks^{9} or microfluidics^{10} have made it possible to probe smaller liquid volumes, but the liquid's density and viscosity cannot be measured simultaneously; one quantity is needed in order to deduce the other from the experimental data. To overcome this limitation, sensors based on microcantilevers have been proposed.^{11–21} These sensors typically require only small volumes (tens of microlitres) of fluid^{22} and are able to determine the viscosity and density simultaneously,^{11,12,15,18–20} making them particularly attractive for labonchiptype diagnostic devices. Measurements effectively quantify changes in the dynamic response of the microcantilever upon immersion into the liquid examined. Fitting the experimental results with theoretical models yields the rheological parameters of the liquid, but the accuracy of the results depends crucially on the quality of the theoretical model, and the ability to implement it fast and robustly. Developing a suitable model is hence far from trivial, because it requires taking into account the coupling between a vibrating cantilever of a given geometry and the surrounding liquid. This is usually characterized by the socalled hydrodynamic function of the cantilever, which in turn depends on the rheological properties of the liquid.
Early developments used significant simplifications such as a spherical model for cantilevers^{23} or an inviscid fluid^{24–26} resulting in large errors or limited applicability. Part of the difficulty comes from the need to take into account the exact geometry of the cantilever and its mechanical properties to precisely determine its hydrodynamic function. By measuring a cantilever's thermal spectrum (that is, its frequency response to the thermal excitation in the liquid) and fitting it to a simple harmonic oscillator model,^{27,28} it is in principle possible to determine the unknown geometrical factors,^{29} but this method fails in highly viscous environments. The first semianalytical model explicitly taking into account the geometry and properties of the cantilever to describe its behaviour in liquid was developed by Sader.^{30} The socalled Sader model provides acceptable agreement with experimental measurements,^{31} but at the cost of computationally intensive calculations and a detailed knowledge of the cantilever. Extending the Sader approach to higher resonances of the microcantilever^{32,33} makes it possible to overcome these difficulties and derive analytical expressions for the cantilever's hydrodynamic function solely based on the frequency of the different resonances.^{34} However, to experimentally reconstruct the hydrodynamic function of a fluid, at least the first three resonance frequencies of the immersed cantilever are needed, something often challenging to measure experimentally in highly viscous liquids.^{35} Once the hydrodynamic function is known, it can be inverted to derive the viscosity and density of the liquid. This step requires some approximations and different methods have been proposed^{15–20} with a typical reported accuracy of 20% when both rheological parameters are determined simultaneously.^{15}
In this paper, we derive novel analytical expressions for the hydrodynamic function of the immersed cantilever. This allows us to propose analytical expressions for the viscosity and the density of the surrounding liquid using only two resonance frequencies of the microcantilever. Significantly, the expressions are fully independent of the type or geometry of the cantilever used, greatly simplifying the derivation of liquid's properties. Our derivation, based on the Euler–Bernoulli beam theory for an immersed cantilever, extends the expression of the hydrodynamic function while considering water as a reference. We test experimentally our analytical expressions over a range of liquids and with several different cantilevers, demonstrating differences smaller than 10% between predictions and experimental results for both density and viscosity. We also investigate the limitations of these models by probing idealised nonNewtonian mixtures of water and poly(ethylene) oxide (PEO). PEO is an uncrosslinked polymer and as such represents a model bodily fluid that displays viscoelastic properties at high molecular weight due to the overlapping of neighbouring polymer coils.^{36–40}
For the purpose of this study, we used an atomic force microscope (AFM) to conduct our measurements. Our results are general and can be implemented in any type of microcantileverbased device. We also expect our findings to contribute towards significant improvements in the booming field of AFM in liquid, in particular for the analysis of surfacecoupled effects on the cantilever vibrations and for the investigation of liquid flow near liquid–solid interfaces.^{41–58}
Materials and methods
In order to verify the accuracy of the proposed expressions and analyse the effect of cantilever parameters and liquid properties on the fluid dissipation mechanisms, thermal spectra were recorded in six wellcharacterised liquids with different viscosities and densities, namely isopropanol, acetone, butanol, decane, bromoform and hexanol. All the liquids were purchased from SigmaAldrich (Dorset, UK) with a purity > 99% and used without further purification. The reference liquid was ultrapure water (18.2 MΩ, MerckMilipore, Dorset, UK). Our model nonNewtonian fluid consisted of varying concentrations of 300000 g mol^{−1} poly(ethylene) oxide (PEO) (SigmaAldrich, Dorset, UK) in water. The PEO was immersed in ultrapure water to a concentration of 3 wt% and dissolved using a magnetic stirrer at 700 RPM for 24 hours until a uniform milky solution was obtained. The solution was then decanted into several 2 mL Eppendorf tubes and centrifuged at 2500 RPM for 10 minutes to separate the PEO solution from the insoluble butylated hydroxytoluene (BHT) that the initial powder contained as an inhibitor. The thenclear fluid was removed with a pipette and bathsonicated for 10 minutes to remove any dissolved air. The PEO solution was then diluted to the required concentration with ultrapure water and the resulting mixture was bathsonicated for 10 minutes to ensure uniformity.
The measurements on pure liquids were conducted on an MFP3D Infinity AFM, and the PEO mixtures were measured on a Cypher ES AFM with a temperaturecontrolled sample stage (both AFMs from Asylum Research, Santa Barbara, CA). For each liquid, we took thermal spectra with four different cantilevers (OMCLRC800PSA, Olympus, Japan). The cantilevers are made of siliconnitride and have different lengths and widths. The nominal geometrical and physical characteristics of the different cantilevers (hereafter referred to as C1–C4) are summarized in Table 1. For each measurement, the cantilever was naturally excited by the Brownian motion of the fluid surrounding it, and its vertical motion was detected using the AFM laser.
Table 1 Summary of the physical characteristics of different cantilevers (C1–C4) used for this study. The cantilevers have a 3 μmhigh tip mounted at one extremity
Cantilever reference 
Width (μm) 
Length (μm) 
Thickness (μm) 
Spring const. (N m^{−1}) 
C1 
40 
100 
0.8 
0.76 
C2 
20 
100 
0.8 
0.39 
C3 
40 
200 
0.8 
0.10 
C4 
20 
200 
0.8 
0.05 
Results and discussion
Theory
The derivation of the analytical expressions for viscosity and density is presented hereafter. For the sake of clarity, only the main steps and the results are presented. The detailed stepbystep derivation of the different results is given in the ESI.†
The immersed microcantilever resonator is assumed to follow the Euler–Bernoulli formalism:

 (1) 
where
E is the cantilever's Youngmodulus,
I is the rotary inertia of cantilever,
ρ_{c} is the cantilever density,
L,
b and
h are the length, width and thickness of the cantilever, respectively,
W(
x,
t) is the timedependent displacement of the cantilever,
F_{exc} is the excitation force and
F_{h} is the hydrodynamic force which can be described by a separate added mass and damping. Considering the added mass and damping per length of the cantilever
^{42,53,54} and assuming a hydrodynamic function characterized by two real (
a_{1},
a_{2}) and two imaginary (
b_{1},
b_{2}) regression coefficients,
^{34} we can relate the resonance frequencies of the cantilever in air,
ω_{an}, and in liquid,
ω_{fn},
^{35,42,54} for any given mode
n (see the ESI
† for details):

 (2) 
where
ρ_{f} and
η are the density and the viscosity of the fluid respectively. Using
eqn (2), the coefficients of the real part of the hydrodynamic function can be calculated if the viscosity and density are known for a reference liquid (usually water). It requires measuring two resonance frequencies of the cantilever in air and in the liquid of interest (see the ESI
†). Experimentally, the measurement of the lower resonance frequencies of the thermal spectra is easier, especially for stiffer cantilevers whose resonance frequencies can be relatively high. We therefore propose using the first two resonance frequencies, obtained directly from the cantilever's thermal spectrum (
Fig. 1).

 Fig. 1 (a) Schematic of a cantilever with relevant dimensions (upper) and illustration of its first and second vibrational modes (lower left and lower right respectively). (b) Example of thermal spectra obtained from a cantilever immersed in ultrapure water (red) and in butanol (blue). Both spectra are displayed on the same log–log plot. The cantilever's first three eigenmodes (ω_{1}, ω_{2}, ω_{3}) can be identified in water. In butanol, a shifting of the resonance frequencies to lower frequencies, together with the broadening of the resonance peaks is visible. The third eigenmode, ω_{3}, can no longer be identified in the baseline noise, highlighting the importance of our model only requiring the first two modes.  
Using the resonance frequencies obtained experimentally for the reference liquid, here water, analytical expressions for the viscosity and density of an unknown fluid can then be derived:

 (3) 

 (4) 
where the indices w, a, 1 and 2 correspond to water, air, the first resonance frequency and the second resonance frequency, respectively.
We note that any direct dependence on the geometrical parameters of the cantilever is cancelled out; the dependence is implicit in the hydrodynamic regression coefficients. Expressions (3) and (4) provide the core results for this paper but the derivation also provides analytical expressions for the quality factor (Q) of the different modes, as well as the added mass and damping to the cantilever at each frequency. Unlike the density and viscosity of the liquid, these quantities are expected to depend on the properties of the cantilever used for the measurement and hence provide a good opportunity to test the quality of the model.
Experiments
In order to calculate the density and viscosity of a fluid from eqn (3) and (4), the two first resonance frequencies of the cantilever while immersed in the unknown fluid had to be measured (Fig. 1(b)). We therefore recorded the thermal spectra of the cantilever in the six liquids of interest, and in ultrapure water for reference. For each liquid, the measurement was repeated with four different cantilevers that exhibited different lengths, widths, and stiffnesses. The cantilevers are hereafter referred to as C1–C4 and their respective characteristics are presented in Table 1 (Materials and methods). This allowed us to examine the impact of the cantilevers’ properties on the calculated ρ and η. Fig. 2 shows the resonance frequencies of each cantilever in different fluids, plotted against the accepted^{25,54,59,60} viscosity (a) and density (b) of each fluid. From Fig. 2, it is clear that there is no obvious relationship between the resonance frequencies and either the viscosity or density of the surrounding liquid. Although increasing density and viscosity tends to reduce the resonance frequencies, the relationship is nonmonotonic, and often different for the two resonances measured. The influence of the cantilever geometry is also evident, with shorter lengths exhibiting higher resonance frequencies. The effect of the cantilever's width is however less pronounced.

 Fig. 2 Measured resonance frequencies of the four different cantilevers in liquids of varying density and viscosity. The first and second resonance frequencies are plotted against the liquids’ viscosities (a) and densities (b). Overall, no obvious correlation is visible, except for a decrease in resonance frequencies with increasing density and viscosity. This decrease is however nonmonotonic and its extent depends on the geometry of the cantilever used for the measurement.  
The accepted density and viscosity values for water at 25 °C are ρ_{w} = 997 kg m^{−3} and η_{w} = 8.94 × 10^{−4} kg m^{−1} s^{−1}.^{59} For each cantilever, it is possible to calculate the hydrodynamic coefficients in water and use the accepted value as a calibration (see the ESI†). This can in turn be used to evaluate the added mass and added damping for the cantilevers in the different liquids investigated (ESI Fig. S1 and S2† respectively).
The accuracy of the derived expressions with water as a reference can be directly evaluated by comparing the measured resonance frequencies and quality factors of the cantilevers with those calculated (eqn (2) and S10†) using the accepted density and viscosity values (Fig. 3). We hereafter refer to the values calculated from accepted densities and viscosities as “literaturecalculated”. Overall, the results show good agreement between the measured and literaturecalculated values; a difference in frequencies of less than 3% except for hexanol, which shows the largest deviation (still <10%). Similarly, quality factors all show deviations smaller than 10% between the derived and measured values.

 Fig. 3 Comparison between the literaturecalculated and measured frequencies ((a), (b)) and quality factors ((c), (d)) for different cantilevers in the test liquids. Lines with a gradient of unity are given as a guide to the eye to facilitate the evaluation. All points measured deviated from the calculated values by less than 3%, apart from hexanol, which has an error of less than 10%.  
Finally, the expressions given by eqn (3) and (4) are used to derive the viscosities and densities of the different test liquids. The values obtained are directly compared with the accepted values for each liquid at 25 °C (ref. 25, 54, 59 and 60) in Fig. 4. As for Fig. 3, the error in the derived values is always smaller than 10%, which represents a significant improvement over previous approaches. Larger errors are incurred for liquids with higher viscosity, as illustrated by the data point for hexanol. This is somewhat to be expected, since the approximation for the hydrodynamic function^{34} is optimized for lower viscosities.

 Fig. 4 Comparison between the accepted and calculated viscosities (a) and densities (b) of the probed fluids. The values of η and ρ derived from the measurements compare well with the accepted values, as evidenced by their collapse onto the line of unity gradient. The inset in (b) highlights the data points at lower densities.  
In practice, it is possible to improve the accuracy of the calculated viscosity and density by including the third resonance frequency in the hydrodynamic function. The density and viscosity of hexanol calculated from the second and third cantilever eigenmodes are shown in Table 2. Comparison with derivations using respectively the 1^{st}/2^{nd} (ρ_{12}, η_{12}) and 2^{nd}/3^{rd} (ρ_{23}, η_{23}) resonance frequencies shows that the latter works better for more viscous fluids.
Table 2 Percentage errors between the calculated and accepted values of density and viscosity for the more viscous liquids using the 1^{st}/2^{nd} and 2^{nd}/3^{rd} resonance frequencies. In all cases the long/thin cantilever was used. The error is reduced by considering the 2^{nd}/3^{rd} resonance frequencies except for the density of butanol where the third resonance is difficult to identify (see Fig 1b)

Error in ρ_{12} (%) 
Error in ρ_{23} (%) 
Error in η_{12} (%) 
Error in η_{23} (%) 
Isopropanol 
3.8 
1 
10 
8 
Butanol 
0.62 
5 
8.4 
5.8 
Hexanol 
7.3 
4.5 
28.2 
5 
The results shown in Table 2 indicate that using the frequencies from the second and third eigenmodes of the cantilever tends to reduce the error. This is due to the fact that higher resonance frequencies are less sensitive to thermal noise when compared to lower modes. A thorough error analysis (ESI Fig S4†) validates this observation. However, using the third eigenmode requires the ability to measure it. Practically, this can be challenging in highdensity fluids, especially for soft cantilevers. The current expression therefore provides a good compromise between accuracy, simplicity and practicality in most technological applications, and can readily be adapted to more viscous liquids if needed.
There is, however, an important point that has not been considered so far. Our model assumes that the viscosity is a scalar quantity and not a function of the probing frequency. In other words, we make the implicit assumption that the liquids probed are Newtonian. This assumption is mostly justified for the test liquids used to validate our model, but this may not hold, for example, for bodily fluids^{4–6} or lubricants.^{1} A deviation from Newtonian behaviour will induce some error in our predictions since the liquid is probed simultaneously at different frequencies, with the second frequency typically 5–6 times higher than the first. This could partially explain the poorer results obtained in the more viscous hexanol. In order to tackle this issue up front, we tested the model in ultrapure water solutions containing increasing concentrations of poly(ethylene) oxide (PEO), a simple uncrosslinked polymer that has been shown to exhibit nonNewtonian properties in aqueous solutions.^{37,39} Specifically, these solutions are shear thinning across a broad range of molecular weights and concentrations,^{40} similarly to most bodily fluids. Practically, this means the cantilevers of different lengths will experience different rheological environments, with the shorter cantilever effectively subjected to a lower viscosity than its longer counterpart, as it oscillates at higher frequencies.
Fig. 5 shows the density and viscosity of various dilutions of PEO in ultrapure water, as calculated using our model, with two cantilevers of different lengths (C1 and C3; subscripts “short” and “long” respectively). As the concentration of PEO increases, the viscosity and density derived from both cantilevers behave similarly; increasing and decreasing, respectively.

 Fig. 5 The calculated density, ρ, and viscosity, η, of different concentrations of PEO in ultrapure water as measured by two different cantilevers. Both η_{Short} and η_{Long} increase with increasing PEO concentration, but the effective viscosity measured by the shorter lever is always higher. The discrepancy between the two cantilevers increases with the PEO concentration due to a frequency dependence of the actual viscosity that is not accounted for by our model. However, η_{Short} agrees with standard rheometer measurements even for the highest PEO concentration measured.^{61} This demonstrates the validity of microcantilever measurements for nonNewtonian liquids, provided there is a suitable choice of a cantilever. The calculated density decreases as the concentration of PEO increases, and the discrepancy between the two cantilevers, although nonmonotonic, increases in a similar manner to the viscosity. Comparison with the directly measured densities of the same weightpercent of PEO (see Fig. S5†) does not show a similar dramatic reduction of both ρ_{Short} and ρ_{Long}, implying that the model's calculated densities are less robust than its viscosities.  
For relatively low PEO concentrations (<1.0 wt%), ρ and η as measured by each cantilever are similar, but at greater concentrations, the discrepancy increases dramatically. This indicates a strong dependence of the calculated values on the cantilever geometry and therefore the resonance frequency, as expected for nonNewtonian fluids. The fact that the observed discrepancy increases with PEO concentration is to be expected given that cantilevers are of different lengths (see Table 1) and therefore resonate at quite different frequencies – for example in ultrapure water, the resonance frequency of the first mode of the short cantilever is more than four times that of the long one. Our model depends on the ratio between the cantilever's eigenmodes and hence is particular sensitive to viscosity and density variations in nonNewtonian liquids.
For the shorter cantilever, the derived viscosity agrees very well with the standard rheometer measurements of PEO in fluid at concentrations of 2 and 3 wt%^{61} while for the longer cantilever the calculated viscosity values are significantly lower. The viscosity as measured by the shorter cantilever is always greater than that obtained from the longer cantilever, for solutions containing PEO. This reflects the fact that for longer cantilevers, the second and third modes of vibration were used as part of our model, due to the first mode being not measureable at high PEO concentration. This is in line with PEO's shearthinning behaviour^{40} since the frequencies used for the longer cantilever are in fact higher than those used for the short cantilever. We therefore expect to find η_{long} < η_{short} as observed.
The measured density decreases monotonically with PEO concentration for both cantilevers, but the relationship between the two measurements is less straightforward than for viscosity. At PEO concentrations of nearly 0.5 wt% and 1.1 wt% the density measured by both cantilevers coincides. These concentrations are greater than the socalled overlap concentration i.e. the concentration above which the polymer coils are dense enough to form transient meshes.^{37} This suggests that the agreement may be due to nonlinearities in our model or possibly errors in determining the resonant frequency (see Fig. S4†), rather than reflecting the intrinsic properties of the polymer solution. At higher concentrations, the apparent reduction in density by a factor of over 3.5 (ρ_{Short}) or over 2 (ρ_{Long}) cannot be correct given the inclusion of only a few weight percent of the polymer. Indeed, independent density measurements, found no such dramatic change (see Fig. S5†). This suggests our model's calculated viscosity to be more reliable than the derived density, the latter becoming unphysical when probing nonNewtonian fluids.
Overall, we find that the discrepancy in our measurements in nonNewtonian PEO is comparable to that of previous methods in pure liquids. Furthermore, the use of several cantilevers provides an effective method for quantifying any deviation from the liquid's Newtonian behaviour, opening new possibilities of diagnostic devices.
Conclusions
In this paper, we present a method to quantitatively determine the viscosity and density of different liquids from the thermal vibration of an immersed microcantilever. We derive analytical expressions based on the Euler–Bernoulli beam theory to quantify the surfacecoupled hydrodynamic effect and deduce the viscosity and density of unknown fluids. Our method only requires measurement of the first two resonance frequencies of the immersed cantilever after calibration in air and in a reference liquid, here water. Significantly, the method implicitly takes into account the cantilever geometry in the calibration process, hence providing analytical expressions for the viscosity and density that are completely independent of the cantilever's characteristics. Experimental validation of the method over an extensive range of liquids yields errors of less than 10% with accepted values. The validity of the model in fluids with frequencydependent viscosities, η(ω), was also investigated using PEO in different concentrations as a model nonNewtonian shearthinning fluid. As expected, the method becomes progressively dependent on the cantilever geometry as the concentration of PEO increases. This is due to the fluid's viscosity becoming more dependent on the frequency as the density of the polymer chains increases. However, this dependence on the cantilever properties can be exploited to quantify the relative error of the measurement. Here, this error was in most cases less than 10% even for the liquid with behaviour comparable to bodily fluids.
We expect our results to contribute primarily to the development of labon chip devices and in nanofluidics. The method could also be used in the field of AFM in liquid, in particular in the analysis of surfacecoupled effects on the cantilever vibrations and for the investigation of liquid flow near liquid–solid interfaces.
Acknowledgements
We acknowledge funding from the Biotechnology and Biological Sciences Research Council (grant BB/M024830/1). We are also grateful to Dr Qing He for providing access to her equipment and Mr Ethan Miller for help with some of the measurements. We are indebted to Dr Richard Thompson and Ms AnnaMarie Stobo for their contributions regarding the viscosity of PEO solutions using standard rheology techniques.
References
 O. Amund and A. Adebiyi, Tribol. Int., 1991, 24, 235–237 CrossRef CAS .
 Y. Xie, H. Dong, S. Zhang, X. Lu and X. Ji, J. Chem. Eng. Data, 2014, 59, 3344–3352 CrossRef CAS .
 L. Phillips, M. Mcgiff, D. Barbano and H. Lawless, J. Dairy Sci., 1995, 78, 1258–1266 CrossRef CAS .
 M. J. Simmonds, H. J. Meiselman and O. K. Baskurt, J. Geriatr. Cardiol., 2013, 10, 291–301 Search PubMed .
 G. Young, B. Sorensen, Y. Dargaud, C. Negrier, K. BrummelZiedins and N. S. Key, Blood, 2013, 121, 1944–1950 CrossRef CAS PubMed .
 P. Ruef, J. Gehm, L. Gehm, J. C. Felbinger, J. Poschl and N. Kuss, Gen. Physiol. Biophys., 2014, 33, 285–293 CrossRef PubMed .

C. Tropea, A. L. Yarin and J. F. Foss, Springer handbook of experimental fluid mechanics, Springer, Berlin, 2007 Search PubMed .
 A. J. Ricco and S. J. Martin, Appl. Phys. Lett., 1987, 50, 1474–1476 CrossRef .
 J. Zhang, C. Dai, X. Su and S. J. O′Shea, Sens. Actuators, B, 2002, 84, 123–128 CrossRef CAS .
 C. J. Pipe and G. H. McKinley, Mech. Res. Commun., 2009, 36, 110–120 CrossRef .
 F. Audonnet and A. A. H. Padua, Fluid Phase Equilib., 2001, 181, 147–161 CrossRef CAS .
 J. L. C. G. da Mata, J. M. N. A. Fareleira, C. M. B. P. Oliveira, F. J. P. Caetano and W. A. Wakeham, High Temp.  High Pressures, 2001, 33, 669–676 CrossRef .
 J. Lee and A. Tripathi, Anal. Chem., 2005, 77, 7137–7147 CrossRef CAS PubMed .
 N. McLoughlin, S. L. Lee and G. Hahner, Lab Chip, 2007, 7, 1057–1061 RSC .
 M. Youssry, N. Belmiloud, B. Caillard, C. Ayela, C. Pellet and I. Dufour, Sens. Actuators, A, 2011, 172, 40–46 CrossRef CAS .
 M. Papi, G. Arcovito, M. D. Spiritio, M. Vassalli and B. Tiribilli, Appl. Phys. Lett., 2006, 88, 194102 CrossRef .
 P. B. Abel, S. J. Eppell, a. M. Walker and F. R. Zypman, Measurement, 2015, 61, 67–74 CrossRef .
 W. Y. Shih, X. Li, H. Gu, W. H. Shih and I. A. Aksay, J. Appl. Phys., 2001, 89, 1497–1505 CrossRef CAS .
 M. F. Khan, S. Schmid, P. E. Larsen, Z. J. Davis, W. Yan, E. H. Stenby and A. Boisen, Sens. Actuators, B, 2013, 185, 456–461 CrossRef CAS .
 B. A. Bircher, L. Duempelmann, K. Renggli, H. P. Lang, C. Gerber, N. Bruns and T. Braun, Anal. Chem., 2013, 85, 8676–8683 CrossRef CAS PubMed .
 M. K. Ghatkesara, E. Rakhmatullinab, H. P. Langa, C. Gerbera, M. Hegnera and T. Brauna, Sens. Actuators, B, 2008, 135, 133–138 CrossRef .
 S. Kim, K. D. Kihm and T. Thundat, Exp. Fluids, 2010, 48, 721–736 CrossRef CAS .
 G. Y. Chen, R. J. Warmack, T. Thundat, D. P. Allison and A. Huang, Rev. Sci. Instrum., 1994, 65, 2532 CrossRef .

W. H. Chu, Vibration of fully submerged cantilever plates I water, Technical report no. 2, DTMB, Southwest Research Institute, San Antonio, Texas, 1963 .
 F. J. Elmer and M. Derier, J. Appl. Phys., 1997, 81, 7709–7715 CrossRef CAS .
 S. Inaba and K. Hane, J. Vac. Sci. Technol., A, 1991, 9, 2138–2139 Search PubMed .
 P. I. Oden, G. Y. Chen, R. A. Steele, R. J. Warmack and T. Thundat, Appl. Phys. Lett., 1996, 68, 3814–3816 CrossRef CAS .
 N. Ahmed, D. F. Nino and V. T. Moy, Rev. Sci. Instrum., 2001, 72, 2731–2734 CrossRef CAS .
 W. Y. Shih, X. Li, H. Gu, W. H. Shih and I. A. Aksay, J. Appl. Phys., 2001, 89, 1497–1505 CrossRef CAS .
 J. E. Sader, J. Appl. Phys., 1998, 84, 64–76 CrossRef CAS .
 J. W. M. Chon, P. Mulvaney and J. E. Sader, J. Appl. Phys., 2000, 87, 3978–3988 CrossRef CAS .
 C. A. Van Eysden and J. E. Sader, J. Appl. Phys., 2006, 100, 114916 CrossRef .
 C. A. Van Eysden and J. E. Sader, J. Appl. Phys., 2007, 101, 044908 CrossRef .
 A. Maali, C. Hurth, R. Boisgard, C. Jai, T. C. Bouhacina and J. P. Aime, J. Appl. Phys., 2005, 97, 074907 CrossRef .
 R. C. Tung, J. P. Killgore and D. C. Hurley, J. Appl. Phys., 2014, 115, 224904 CrossRef .
 B. R. Dasgupta, S. Y. Tee, J. C. Crocker, B. J. Frisken and D. A. Weitz, Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top., 2002, 65, 051505 CrossRef PubMed .
 E. C. Cooper, P. Johnson and A. M. Donald, Polymer, 1991, 32, 2815–2822 CrossRef CAS .
 D. M. Woodley, C. Dam, H. Lam, M. LaCave, K. Devanand and J. C. Selser, Macromolecules, 1992, 25, 5283–5286 CrossRef CAS .
 J. H. van Zante, S. Amin and A. A. Abdala, Macromulecules, 2004, 37, 3874–3880 CrossRef .
 K. W. Ebagninina, A. Benchabaneb and K. Bekkour, J. Colloid Interface Sci., 2009, 36, 360–367 CrossRef PubMed .
 A. F. Payam, Ultramicroscopy, 2013, 135, 84–88 CrossRef PubMed .
 A. F. Payam and M. Fathipour, Micron, 2015, 70, 50–54 CrossRef PubMed .
 H. Lee and W. J. Chang, Micron, 2016, 80, 1–5 CrossRef CAS PubMed .
 R. C. Tung, J. P. Killgore and D. C. Hurley, Rev. Sci. Instrum., 2013, 84, 073703 CrossRef PubMed .
 A. Raman, J. Melcher and R. C. Tung, Appl. Phys. Lett., 2013, 103, 263702 CrossRef .
 K. Voïtchovsky, Nanotechnol., 2015, 26, 100501 CrossRef PubMed .
 K. Voïtchovsky, J. J. Kuna, S. A. Contera, E. Tosatti and F. Stellacci, Nat. Nanotechnol., 2010, 5, 401–405 CrossRef PubMed .
 K. Voïtchovsky and M. Ricci, Proc. SPIEInt. Soc. Opt. Eng., 2012, 8232, 82320O CrossRef .
 M. Ricci, P. Spijker and K. Voïtchovsky, Nat. Commun., 2014, 5, 4400 CrossRef CAS PubMed .
 D. OrtizYoung, C. H. Chih, S. Kim, K. Voïtchovsky and E. Riedo, Nat. Commun., 2013, 4, 2482 CrossRef PubMed .
 J. Melcher, C. Carrasco, X. Xu, J. L. Carrascosa, J. G. Herrero, P. J. de Pablo and A. Raman, Proc. Natl. Acad. Sci. U. S. A., 2009, 106, 13655–13660 CrossRef CAS PubMed .
 A. F. Payam, J. R. Ramos and R. Garcia, ACS Nano, 2012, 6, 4663–4670 CrossRef CAS PubMed .
 D. Kiracofe and A. Raman, J. Appl. Phys., 2010, 107, 033506 CrossRef .
 S. Basak, A. Raman and S. V. Garimella, J. Appl. Phys., 2006, 99, 114906 CrossRef .
 M. Lee, B. Kim, Q. H. Kim, J. G. Hwang, S. An and W. Jhe, Phys. Chem. Chem. Phys., 2016, 18, 27684–27690 RSC .
 R. Garcia and R. Proksch, Eur. Polym. J., 2013, 49, 1897–1906 CrossRef CAS .
 M. H. Korayem and N. Ebrahimi, J. Appl. Phys., 2011, 109, 084301 CrossRef .
 M. H. Korayem and R. Ghaderi, Sci. Iran., 2013, 20, 1 Search PubMed .
 J. Kestin, M. Sokolov and W. A. Wakeham, J. Phys. Chem. Ref. Data, 1978, 7, 941–948 CrossRef CAS .
 B. Garcia, R. Alcalde, S. Aparicio and J. M. Leal, Phys. Chem. Chem. Phys., 2002, 4, 1170–1177 RSC .

R. Thompson, Private communication, 2017.
Footnote 
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6an02674e 

This journal is © The Royal Society of Chemistry 2017 