We obtain quantitative measurements of the oscillation amplitude of vibrating objects by using sideband digital holography. The frequency sidebands on the light scattered by the object, shifted by n times the vibration frequency, are selectively detected by heterodyne holography, and images of the object are calculated for different orders n. Orders up to n = 120 have been observed, allowing the measurement of amplitudes of oscillation that are significantly larger than the optical wavelength. Using the positions of the zeros of intensity for each value of n, we reconstruct the shape of vibration the object.
© 2009 Optical Society of America
Holographic imaging is based on interferences between a signal optical field and a reference beam; it provides an accurate method to image the vibration of various objects. Powell and Stetson  have shown that the time-averaged hologram of an harmonically vibrating object involves the Bessel function J 0(z), where z is the phase modulation amplitude. In the backscat-tering geometry, z = 4πA/λ where A is the mechanical amplitude of vibration and λ the optical wavelength. Image reconstruction then provides a direct mapping of the amplitude A at various points of the object, in the form of dark fringes appearing around points where J 0(z) is zero.
Picard et al.  have simplified the processing of the data by performing Time Averaged Digital Holography (TADH) with a CCD camera, leading to numerous recent developments [3, 4, 5, 6, 7]. Nevertheless, the quantitative measurement of vibration amplitudes remains limited to situations where fringes can be counted, which implies that the amplitudes of vibration must be relatively small; otherwise, it becomes difficult to count many narrow fringes, and even impossible when they are smaller than the optical resolution.
In this letter, we show how Sideband Digital Holography (SDH) can overcome the problem. Sideband holography was already demonstrated in 1971 by C.C. Aleksoff, in a pioneering experiment  where the reference beam was frequency-shifted by an ultrasonic diffraction cell; selecting a given frequency resulted in a selection of a given sideband of the signal backscattered by the object. The experiment was done with holographic plates. Here, we use digital heterodyne holography ; both reference and illumination laser beams are frequency shifted, and the selection of one sideband is obtained by a proper detuning of one of the beams followed by heterodyne detection on a CCD camera; this technique provides much more flexibility and leads to accurate quantitative measurements. As a demonstration, inspired by the work described in [10, 11] we perform an experiment where the vibrating object is the reed of a clarinet; in addition to the intrinsic interest of such an object, with its possible musical implications, a reed provides a test system that is particularly well adapted to our purposes, with typical vibration amplitudes of the order of 0.1 mm. We will see that, even with smaller amplitudes, there is no difficulty in obtaining holographic images corresponding to the n th sideband with n up to 100 or more.
The experimental setup is schematically shown in Fig. 1. A laser beam, with wavelength λ ≃ 650 nm (angular frequency ωL) is split into a local oscillator beam (ELO) and an illumination beam (EI); their angular frequencies ωLO and ωI are tuned by using two acousto-optic modulators (Bragg cells with a selection of the first order diffraction beam) AOM1 and AOM2: ωLO = ωL + ω AOM1 and ωI = ωL + ω AOM2 where ω AOM1,2/2π ≃ 80 MHz. The clarinet reed is attached to a clarinet mouthpiece and its vibration is driven by a sound wave propagating inside the mouthpiece, as in playing conditions, but in our experiment the sound wave is created by a loudspeaker excited at frequency ωA and has a lower intensity than inside a clarinet; no attempt has been made to reproduce the mechanical effect of the lip of the player. The excitation frequency is adjusted to be resonant with the first flexion mode (2143 Hz) of the reed. The phase of the field E backscattered by any element of the reed is then phase modulated at frequency ωA, so that E can be expanded into carrier (E 0) and sideband components (En) as:
where ℰ is the complex field of the object in the absence of modulation, z = 4πA/λ is the phase modulation amplitude, A the amplitude of vibration of the particular element, and Jn is the n-th order Bessel function of the first kind (with J -n(z) = (-1)n Jn(z) for integer n and real z).
To selectively detect the n th sideband we use 4 phase shifting holography  and adjust the ω AOM1 to fulfil the condition:
where ω CCD is the CCD camera frame frequency and n is a (positive or negative) integer; the phase of the modulated signal is then shifted by 90° from one CCD image to the next . The sideband complex hologram signal Hn provided by each pixel of the camera, proportional to the local sideband complex field, is obtained by 4-phase demodulation:
where I 0, ⋯I 3 are 4 consecutive intensity images digitally recorded by the CCD camera (Pike: Allied vision Technology Inc, 12 bits, ω CCD/2π = 10 Hz, exposure time 100 ms, 1340 × 1024 pixels). The LO beam is slightly angularly shifted (~ 1°) with respect to the beam originating from the reed; in this way, the zero order and twin-image aliases  can be suppressed and an ultimate detection sensitivity can be reached [15, 16].
From the complex holograms, the images of the reed are reconstructed by the standard convolution method . To simplify the Fourier transform calculation, the 1340 × 1024 pixels data are truncated to 1024 x 1024. Because of the angular shift, the image of the reed lies in the upper side of the reconstruction grid. In Fig. 2, for pixels ranging from 200 to 1023 along the x axis and 0 to 511 on the y axis, we show the reconstructed images for different sideband orders n; only the square modulus of the complex signal is retained (the phase coming from ℰ is ignored), so that we obtain intensity images, proportional to ∣En∣2 for each different sideband.
Figure 2(a) is obtained at the unshifted carrier frequency (n = 0). The left side of the reed is attached to the mouthpiece, and the amplitude of vibration is larger at the tip of the reed on the right side; in this region the fringes become closer and closer and difficult to count. The mouthpiece is also visible, but without fringes since it does not vibrate. Similar images of clarinet reeds have been obtained in [10, 11], with more conventional techniques and lower image quality. Figures 2(b)–2(d) show images obtained for three sideband frequencies. As expected, the non-vibrating mouthpiece is no longer visible. Figure 2(b) shows the n = 1 sideband image, with J 1 fringes that are slightly shifted with respect to those of J 0. Figure 2(c) shows the order n = 10 and (d) the order n = 100 (e). The left side region of the image remains dark because, in that region, the vibration amplitude is not sufficient to generate these sidebands, Jn(z) being evanescent for z < n.
We have performed cuts of the reconstructed images signal ∣En(x,y∣2 obtained for different sideband orders n along the horizontal line y = 174 (this value has been chosen because it corresponds to a region where the fringes are orthogonal to the y axis). To reduce the effect of speckle fluctuations, the intensity signal ∣En(x,y = 174) ∣2 has been averaged over 20 pixels in the y direction (y = 164… 184). Moreover, since the illumination of the reed is not uniform, we have normalized the sideband signal by the reconstructed image intensity ∣E∣2 obtained at the carrier frequency (n = 0) without mechanical excitation of the reed (loud-speaker off). One example of the normalized signals ∣En∣2/∣E∣2 along the cut is shown for n = 5, in the lower part of Fig. 3, with a downwards intensity axis. In this case, and for many other values of n, we have checked that the normalized curves vary as the square of Bessel functions, with a first fringe that moves to higher values when n increases (it occurs when z ≃ n and the amplitude of vibrations increases with x). For instance, for n = 0, which corresponds to standard TADH, the variations of ∣J 0(z)∣2 show a large number of nodes and anti-nodes, with decreasing visibility for increasing x; above x = 700, the fringes are no longer visible. For n = 20, the first fringe, which corresponds to z ≃ n = 20, is located near x ≃ 560. Going to n = 60 pushes the first fringe to x ≃ 755, where the counting of fringes with TADH would no longer be possible. In a more general way, the first fringe on the nth sideband image provides a convenient marker for the region z ≃ n.
For a more precise treatment of the data, we use the position of the antinodes, which is insensitive to inhomogeneities of illumination; no normalization is then required. To build the central part of Fig. 3, a similar linear cut at y= 174 is processed, for the values n = 0,1,5,10,… 100. The x locations of the successive minima of the signal ∣En(x)∣2 are plotted against the zeros of the corresponding ∣Jn(z)∣2. To illustrate the way the points are obtained, we show in the left part of the Fig. 3 the expected signal ∣J 5(z)∣2, to be compared to the signal in the lower part. For all observed sideband order n ranging from 0 to 100, all the points displayed on Fig. 3 fall on a well defined curve, which provides the vibration amplitude z as function of the location x. The excellent consistency of the overlapping data sets demonstrates the validity of expansions Eq. 1 for any order n. Therefore the figure provides a direct and accurate visualization of the map of the maximal elongations of the reed, from the left part clamped on the mouthpiece to the tip on the right; here, the maximum amplitude is z ≃ 120 radians, corresponding to A ≃ 6.2 μm.
To conclude, recording separately different sideband Fourier components of the light back-scattered by a vibrating object gives access to a large amount of information; this can be used to remove ambiguities and inaccuracies in measurements of large modulation amplitude. Using large sideband order n actually makes the calibration of the vibration amplitude z straightforward: the first fringe of a high order image can be used as a marker. This marker is easy to locate since it corresponds to a fast transition from zero signal, and moreover is brighter and broader. Performing sideband holography with many orders n gives a direct and accurate access to the shape of the vibration amplitude and its position dependence. Figure 3 shows only a representation of a slice of the data along one given axis, but more information is available if one changes the value of y or chooses another direction, for instance parallel to y to accurately map the transverse variations of the amplitude A. Sideband digital holography therefore opens up a variety of new possibilities in holography.
The authors acknowledge the support of French National Research Agency (ANR-05-NANO-031) and the “Centre de compétence NanoSciences Île de France” (C’nano IdF).
References and links
1. R. Powell and K. Stetson, “Interferometric vibration analysis by wavefront reconstruction,” J. Opt. Soc. Am 55, 1593–1598 (1965). [CrossRef]
3. F. Zhang, J. Valera, I. Yamaguchi, M. Yokota, and G. Mills, “Vibration Analysis by Phase Shifting Digital Holography,” Opt. Rev. 11, 297–299 (2004). [CrossRef]
5. V. Singh, J. Miao, Z. Wang, G. Hegde, and A. Asundi, “Dynamic characterization of MEMS diaphragm using time averaged in-line digital holography,” Opt. Commun. 280, 285–290 (2007). [CrossRef]
6. N. Demoli and D. Vukicevic, “Detection of hidden stationary deformations of vibrating surfaces by use of time-averaged digital holographic interferometry,” Opt. Lett. 29, 2423–2425 (2004). [CrossRef] [PubMed]
8. C. C. Aleksoff, “Temporally modulated holography,” Appl. Opt. 10, 1329–1341 (1971). [PubMed]
9. F. LeClerc, L. Collot, and M. Gross, “Numerical Heterodyne Holography Using 2D Photo-Detector Arrays,” Opt. Lett. 25, 716–718 (2000).
14. E. Cuche, P. Marquet, and C. Depeursinge, “spatial filtering for zero-order and twin-image elimination in digital off-axis holography,” Appl. Opt. 39, 4070–4075 (2000). [CrossRef]