Detection of Instantaneous Phase. The cantilever signal is
where \(x_{\text{rms}}\) is the cantilever root mean square amplitude, \(\omega_0\) is the cantilever frequency, and \(\phi\) is the cantilever phase. Here \(\delta x(t)\) is random noise which includes contributions from cantilever thermomechanical fluctuations as well as detector noise.
In order to detect the cantilever frequency we create a quadrature signal by taking the Hilbert transform of the cantilever signal. This procedure gives
where \(\delta y(t)\) is the Hilbert transform of \(\delta x(t)\). An expression for \(\delta y(t)\) can be written down, but it is not instructive. There is a simple relation, however, between \(y\) and \(x\) in the Fourier domain:
where \(\widehat{\delta x}(f)\) indicates the Fourier transform of \(\delta x(t)\). The function \(H\) implements the Hilbert transform in Fourier space:
Since \(H(f) H^{*}(f) = 1\) (except for the single point at \(f=0\)), it follows that \(\delta y(t)\) has essentially the same power spectrum as \(\delta x(t)\).
In our frequency-detection algorithm we measure the instantaneous phase of the cantilever using
Substituting equations (1) and (2) into equation (3),
Let us now, with the help of Mathematica, expand \(\phi(t)\) in a Taylor series to first order in both \(\delta y(t)\) and \(\delta x(t)\). The result is
We can extract the instantaneous frequency as the slope of the \(\phi(t)\) versus \(t\) line. After subtracting away the best-fit line, we are left with phase noise
given by
Phase Noise Power Spectrum. Taking the Fourier transform of \(\delta \phi(t)\), and switching frequency units
Which can be simplified to
We can eliminate \(\widehat{\delta y}\) from equation (5) by recognizing
which holds when \(-f_0 < f < f_0\); we can arrange for this condition to be met by applying a bandpass filter to the cantilever signal. Substituting equations (6) and (7) into equation (5) gives
Passing to the power spectrum requires a limiting procedure, as follows. We should consider that \(x(t)\) is only sampled for a finite amount of time \(T\), which we can indicate with a subscript: \(x(t) \rightarrow x_{T}(t)\) where
Equation (4) holds with \(\delta x \rightarrow \delta x_T\), \(\delta x \rightarrow \delta y_T\), and \(\delta \phi \rightarrow \delta \phi_T\). Time correlation functions are defined in terms of \(x_T(t)\), not \(x(t)\),
where \(\langle \cdots \rangle\) indicates a statistical average. The manipulations leading to equation (8) are still valid with the \(T\)-subscripted variables, with the result that
The next step to computing the power spectrum is to calculate
We may now pass to the power spectrum by taking the limit
with the power spectrum \(P_{\delta \phi}(f)\) analogously defined. Carrying out this limiting procedure on both sides of equation (11) yields
where \(\text{Re} \! \left( \cdots \right)\) indicates taking the real part. The last term will not survive statistical averaging over the phase \(\phi\) since
Implicit in this average is the assumption that \(\phi\) is randomly distributed, that is, there is no correlation between the phase of the cantilever and the cantilever noise. After statistical averaging over \(\phi\), the power spectrum of cantilever phase noise becomes
Frequency Shift Power Spectrum. Let us define the instantaneous frequency shift as
and compute the power spectrum of the instantaneous frequency shift. Let us define \(\delta f_{T}(t)\) as in equation (9). The time-correlation function of the frequency shift is then
with \(C_{\delta \phi}\) defined likewise. Substituting, and dropping \(\langle \cdots \rangle\) for notational convenience,
The time derivative \(\delta \dot{\phi}\) may be computing from its Fourier transform. With
we can compute the time derivative of the instantaneous phase shift as
If we substitute equation (14) into equation (13) and use
where \(\delta(t)\) is the Kroenecker delta function, then
where we have passed the limit into the integral. Because \(\delta \phi_T(t)\) is a real function,
The term in braces is thus \(P_{\delta \phi}(f)\), the power spectrum of phase fluctuations. We find
Comparing this to the usual relation between the correlation function and the power spectrum
we see that
We have used that \(P_{\delta \phi}(-f) = P_{\delta \phi}(f)\). Substituting equation (15) into equation (12) we conclude that position fluctuations lead to frequency noise having a power spectrum
Instrument Noise. Equation (16) is a general relation between the position-fluctuation power spectrum and the frequency-fluctuation power spectrum. The power spectrum of detector noise is typically flat:
Within this approximation,
This relation holds whether the power spectra are defined as one-sided or two-sided, as long as the power spectrum is computed consistently on both sides of equation. We typically work up data using a one-sided power spectrum. The more general equation (16) can be used when the detector noise spectrum is not independent of frequency.
Cantilever Thermomechanical Fluctuations. We have previously shown that the (one sided) power spectrum of cantilever position fluctuation is
where \(T\) is temperature, \(k_b\) is Boltzmann’s constant, and \(f_0\), \(\tau_0\), and \(\Gamma\) are cantilever frequency, ring-down time, and dissipation constant, respectively. For frequency offsets \(f \gg f_0 / Q\) we find that
Substituting this result into equation (16) gives
The last term equals \(Q^{-2}\), where \(Q\) is the cantilever quality factor. Using \(\Gamma = k /(2 \pi f_0 Q)\) we can rewrite the one-sided power spectrum of cantilever frequency fluctuations as
Discussion. Equations (17) and (19) agree exactly with what Loring and co-workers have derived [1]. Together, thermomechanical fluctuations and detector noise lead to cantilever frequency noise with a one-sided power spectrum of
This equation is valid for offset frequencies \(f \gg f_0/Q\) and assumes for simplicity that detector noise is frequency independent in the vicinity of the cantilever resonance frequency.
References
[1] | Yazdanian, S. M.; Marohn, J. A. & Loring, R. F. Dielectric Fluctuations in Force Microscopy: Noncontact Friction and Frequency Jitter. J. Chem. Phys., 2008, 128: 224706 [http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2674627/] [http://dx.doi.org/10.1063/1.2932254] . See equations 6.7 through 6.9. |