Phase-sensitive narrowband heterodyne holography

Phase-sensitive narrowband heterodyne holography

Francois Bruno Institut Langevin. Fondation Pierre-Gilles de Gennes. Centre National de la Recherche Scientifique (CNRS) UMR 7587, Institut National de la Santé et de la Recherche Médicale (INSERM) U 979, Université Pierre et Marie Curie (UPMC), Université Paris 7. Ecole Supérieure de Physique et de Chimie Industrielles (ESPCI) - 1 rue Jussieu. 75005 Paris. France    Jean-Baptiste Laudereau Institut Langevin. Fondation Pierre-Gilles de Gennes. Centre National de la Recherche Scientifique (CNRS) UMR 7587, Institut National de la Santé et de la Recherche Médicale (INSERM) U 979, Université Pierre et Marie Curie (UPMC), Université Paris 7. Ecole Supérieure de Physique et de Chimie Industrielles (ESPCI) - 1 rue Jussieu. 75005 Paris. France    Max Lesaffre Institut Langevin. Fondation Pierre-Gilles de Gennes. Centre National de la Recherche Scientifique (CNRS) UMR 7587, Institut National de la Santé et de la Recherche Médicale (INSERM) U 979, Université Pierre et Marie Curie (UPMC), Université Paris 7. Ecole Supérieure de Physique et de Chimie Industrielles (ESPCI) - 1 rue Jussieu. 75005 Paris. France    Nicolas Verrier Institut Langevin. Fondation Pierre-Gilles de Gennes. Centre National de la Recherche Scientifique (CNRS) UMR 7587, Institut National de la Santé et de la Recherche Médicale (INSERM) U 979, Université Pierre et Marie Curie (UPMC), Université Paris 7. Ecole Supérieure de Physique et de Chimie Industrielles (ESPCI) - 1 rue Jussieu. 75005 Paris. France    Michael Atlan Institut Langevin. Fondation Pierre-Gilles de Gennes. Centre National de la Recherche Scientifique (CNRS) UMR 7587, Institut National de la Santé et de la Recherche Médicale (INSERM) U 979, Université Pierre et Marie Curie (UPMC), Université Paris 7. Ecole Supérieure de Physique et de Chimie Industrielles (ESPCI) - 1 rue Jussieu. 75005 Paris. France
July 10, 2019

We report on amplitude and phase imaging of out-of-plane sinusoidal surface vibration at nanometer scales with a heterodyne holographic interferometer. The originality of the proposed method is to make use of a multiplexed local oscillator to address several optical sidebands into the temporal bandwidth of a sensor array. This process is called coherent frequency-division multiplexing. It enables simultaneous recording and pixel-to-pixel division of sideband holograms, which permits quantitative wide-field mapping of optical phase modulation depths. Additionally, a linear frequency chirp ensures the retrieval of the local mechanical phase shift of the vibration with respect to the excitation signal. The proposed approach is validated by quantitative motion characterization of the lamellophone of a musical box, behaving as a group of harmonic oscillators, under weak sinusoidal excitation. Images of the vibration amplitude versus excitation frequency show the resonance of the nanometric flexural response of one individual cantilever, at which a phase hop is measured.

Imaging nanometric optical path length modulations is of interest in non-destructive testing of micro electro-mechanical systems Bosseboeuf2003 (); Rembe2006 (); Kokkonen2008 (); BramhavarPouet2009 (). Among the approaches used for single-point vibration spectra measurements, laser Doppler schemes are commonly found; they enable non-contact out-of-plane vibration monitoring WhitmanKorpel1969 (); DeLaRue1972 (); Stegeman1976 (). In these measurement schemes, the resolution of the minimum measurable vibration amplitude is limited by the detection sensitivity WhitmanKorpel1969 (); Monchalin1985 (); WagnerSpicer1987 (); RoyerDieulesaint1989 (); Rembe2006 (); Kokkonen2008 (). These methods exhibit high reliability for wideband RoyerDieulesaint1986 (); RoyerDieulesaint1989 (); JiaBoumiz1993 (); RoyerKmetik1992 () and narrowband DeLaRue1972 (); Monchalin1985 (); BramhavarPouet2009 () measurements. Narrowband measurements facilitate high frequency monitoring without the need for specialized electronics Kokkonen2008 (); BramhavarPouet2009 (), while wideband measurements allow for transient signal analysis RoyerDieulesaint1986 (). Nevertheless, imaging is a tedious and time-consuming process because laser probe scanning is required. To circumvent scanning issues, wide-field heterodyne optical detection schemes with dedicated sensor arrays were introduced recently Kimachi2010 (); PatelAchamfuoYeboah2011 (). Another approach involves time-averaged holographic interferometry, which permits high sensitivity measurements Powell1965 (); Aleksoff1969 (); Stetson1970 (); Levitt1976 (); UedaMiida1976 (); PicartLeval2003 (). This method enables narrowband measurements of the optical pathlength modulation depth when the exposure time of the recorded interferogram is much longer than the vibration period Aleksoff1969 (); PsotaLedl2012 (); VerrierAtlan2013 (). With this approach, the optical phase retardation of the recorded scattered field beating against a reference optical field can also be exploited for time-resolved (wideband) determination of the object motion Pedrini2006 (); PerezLopez2006 (). Nevertheless, typical vibration frequencies of interest in non-destructive monitoring are still much higher than camera frame rates, and hence the time-averaging condition holds. In this regime, the phase shift with respect to the excitation signal is usually retrieved from phase-stepped stroboscopic schemes Lokberg1976 (); Petitgrand2001 (); Bosseboeuf2003 (); LevalPicart2005 (); VerrierGrossAtlan2013 (). Quantitative imaging of optical path length modulation and phase retardation in stationary regimes, at the time scale of the frame exposure, may thus find applications.

Figure 1: Experimental set-up.

In this letter, we report on a heterodyne holographic arrangement used for amplitude and phase imaging of optical path length modulation in time-averaged recording conditions, without strobe light. The proposed arrangement consists of a Mach-Zehnder holographic interferometer in an off-axis and frequency-shifting configuration AtlanGross2007JOSAA (); VerrierAtlan2013 (), with a multiplexed local oscillator (Fig. 1). Frequency-tunable narrowband detection is achieved in time-averaged heterodyne detection conditions. Time-averaged holograms of a phase-modulated optical field are recorded with an exposure time much longer than the modulation period. Optical heterodyning is performed with a frequency-shifted optical local oscillator, a frequency-conversion process aimed at shifting a given radiofrequency (RF) optical sideband in the sensor bandwidth. Optical local oscillator multiplexing is realized by the coherent addition of RF signals driving an acousto-optic modulator in the reference beam channel. Coherent frequency-division multiplexing is a technique by which the total available bandwidth is divided into non-overlapping frequency sub-bands carrying separate signals. Addressing simultaneously distinct domains of the spectral bandwidth of holograms was performed with success for angular PaturzoMemmolo2009 (), wavelength Kiire2012 (), space-division Tahara2012 (), and sideband VerrierAtlan2013 () multiplexing schemes. Coherent frequency-division multiplexing in heterodyne holography offers the opportunity to measure complex-valued maps of the RF spectral components of the optical radiation field with accuracy and sensitivity, as long as they can be isolated in sub-bands of the sensor’s temporal bandwidth. In the reported experiment, high-quality factor flexural responses at audio frequencies of a lamellophone (the metallic comb of cantilevers from a musical box) excited in a sympathetic resonance regime are sought. Sympathetic vibration is a harmonic phenomenon wherein a vibratory body responds to external excitation in the vicinity of a resonance frequency.

The optical setup used to perform non-destructive testing of the cantilevers’ out-of-plane flexural responses at nanometric scales is sketched in Fig. 1. It is designed to monitor quantitatively steady-state vibration amplitudes and to detect mechanical phase jumps with respect to the excitation signal, by interferometric detection of an optical field scattered by the lamellophone, referred to as the object field, beating against a reference optical field. The out-of-plane motion of the cantilevers provokes an optical path length modulation of the backscattered object field. Hence the temporal part of the object field undergoes a sinusoidal phase modulation . The equation of the out-of-plane motion of a cantilever, modeled as a damped harmonic oscillator of spring constant , mass , damping coefficient , driven sinusoidally by the force provided by a piezo-electric actuator is


We seek a steady-state solution at the excitation frequency with an induced phase change of with respect to the excitation , as a result of viscous damping


Its amplitude, , is proportional to the driving force


Two intermediate variables were introduced : the angular resonance frequency of the undamped oscillator, and the dimensionless damping ratio . They satisfy the relationships and . The phase shift of the sinusoidal response with respect to the driving force is


The resonance frequencies of each cantilever are almost equal to the resonant frequencies of the undamped system. The monochromatic optical field of amplitude , oscillating at the angular frequency , backscattered by the object in vibration can be written


The phase-modulated field can be decomposed on a basis of Bessel functions of the first kind , via the Jacobi-Anger identity


where the quantity


is the complex weight of the optical sideband of order , and is the modulation depth of the optical phase of the radiation backscattered by the cantilever. The local amplitude of the out-of-plane vibration and phase retardation with respect to the sinusoidal driving force can be derived from the complex weight of two optical sidebands with Eq. 7 and the first-order Taylor developments of and , near , for

Figure 2: Magnitude of a statically-scattered light hologram (a). Magnitude of sideband holograms (b) and (c). Amplitude map of the flexural response of the first cantilever for the first resonance at (d). Phase images calculated from 128 raw interferograms around 540.9 Hz (e) and 541.8 Hz (f), in the neighborhood of the first resonance. Movies of the amplitude and phase maps of the out-of-plane vibration are reported in media 1 and 2.

The holographic interferometer is used for the detection of RF phase modulation sidebands of an object field in reflective geometry, beating against a LO field , as reported in VerrierAtlan2013 (). Its purpose is to enable a robust quantitative measurement of the nanometric modulation amplitudes of cantilever resonances, as well as the associated phase shift with respect to the sinusoidal excitation signal. To do so, the adopted strategy is to separate the sensor bandwidth into three non-overlapping frequency sub-bands, centered at , , and , each of which is used to carry a separate optical modulation sideband. The LO signal is composed of the sum of three phase-locked RF signals, to yield an optical LO field of the form


where is the frequency shift and is the complex weight of the LO component of rank . The positive parameters are the normalized weights of each LO component; in this experiment . The LO components are tuned around the optical modulation bands of order , , and . Precisely, the LO component of rank beats against the optical modulation band of rank at the apparent frequency , which is set within the camera bandwidth. The sensor array used for detection is an EMCCD camera, (Andor IXON 885+, frame rate ). The main optical radiation field is provided by a doubled Nd:YAG laser (Oxxius SLIM 532, power 100 mW, wavelength nm, optical pulsation ). Image rendering from the raw interferograms impinging on the sensor array involves a diffraction calculation performed with a numerical Fresnel transform. Its implementation in off-axis holographic conditions Cuche2000 () yields complex-valued time-averaged Powell1965 (); Aleksoff1969 (); Stetson1970 (); Levitt1976 () holograms from which the off-axis region


is the support of the object-reference fields cross-term. The demodulation process (Eq. 12) from of consecutive recordings of finite exposure times acts as a band-pass filter of width in the temporal frequency domain. Only the contributions within the sensor temporal bandwidth, between the Nyquist frequencies , have to be taken into account in , which can hence be expressed as

Figure 3: Magnitude of the FFT-spectrum (Eq. 12 and Eq. 13), where three modulation sidebands are addressed.

where the quantities satisfy . The FFT-spectrum of consecutive off-axis holograms


exhibits three peaks at the frequencies , , and (Fig. 3). Their complex values yield holograms , of rank , whose magnitude is proportional the modulation depth at the frequency


where is a complex constant. Thus, is a measure of the complex weight . For small-amplitude vibrations, a quantitative map of the local out-of-plane vibration amplitude can be assessed from the ratio of the first-order sideband and the non-shifted light component


In the same manner, a map of the local mechanical phase retardation with respect to the excitation signal can be assessed from , and


where . Relationships 14 and 15 are derived from Eq. 8, and Eq. 13. One can also make use of the third recorded band and the property to derive .

To perform the measurement, an optical LO consisting of the combination of three frequency components shifted by , , and (Eq. 9) was realized. Those shifts were chosen to yield non-symmetrical (non-opposed) algebraic modulation frequencies , , and in the detector bandwidth in order to avoid spurious aliasing or cross-talk effects. To do so, we made use of phase-locked RF signals generated by frequency-synthesizers at the carrier frequency , , , and , where the shifts used to generate the three RF bands were , , (see Tab. 1). As sketched in Fig. 1, those signals were mixed and combined to drive a set of acousto-optic modulators (AA Opto Electronic) designed to operate around , from which opposite diffraction orders were selected. The resulting modulation frequencies of the holograms (Eq. 13) were , , . The three sideband components , , and , encoded at 5 Hz, -2.5 Hz, and -7.5 Hz respectively, are visible on the magnitude of the discrete Fourier spectrum (Eq. 13) reported in Fig. 3. Spectral maps of the vibration amplitude were calculated from a set of sequential measurements at each frequency. We took 3000 sequences of raw interferograms sampled at a frame rate of , for excitation frequencies ranging from 0 Hz to 3 kHz, in 1 Hz steps (supply voltage of the piezo electric actuator : 10 mV). The first resonance frequencies of the lamellophone’s cantilevers are observed between 500 Hz and 2500 Hz (Media 1). Plots of the vibration amplitude of the 1st, 5th, and 17th cantilevers versus the excitation frequency are reported in Fig. 4. In particular, individual resonances of the 1st, 5th, and 17th cantilevers were identified around 541 Hz (Fig. 4(a), inset), 1006 Hz (Fig. 4(b), inset), and 2211 Hz (Fig. 4(c), inset).

angular frequency value
Table 1: Frequency shifts used in the experiments.
Figure 4: Vibration amplitude versus excitation frequency , averaged over the 1st (a), 5th (b), and 17th cantilevers. Insets: vibration amplitude maps at 541 Hz (a), 1006 Hz (b), and 2211 Hz (c).

A finer screening around the resonance of the first cantilever was performed, from 248 sequences of raw interferograms, for excitation frequencies ranging from 536 Hz to 546 Hz. The values of , spatially averaged over the extent of the resonant cantilever, arrowed in Fig. 2(b), are marked as points in the plot reported in Fig. 5(a). But is randomized from the measurement of one frequency point to the next. Hence the phase retardation cannot be retrieved. To prevent reference phase randomization between measurements, a chirp of the frequency was performed. It consisted of a linear drift with time of the excitation frequency from to in


Two out of the three components of the optical LO (Eq. 9) were chirped as well via throughout the acquisition. During the time lapse , a sequence of raw interferograms was steadily recorded. In that manner, remained constant. Amplitude and phase maps at each frequency were calculated from the demodulation of consecutive raw frames, reported in media 2. The complex values of these maps were averaged over the spatial extent of the first cantilever and reported in Fig. 5 (blue lines). In agreement with Eq. 4, for non-zero viscosity, a phase shift of is observed at the resonance frequency . The theoretical resonance lines in Fig. 5 were determined for a resonance frequency of the undamped oscillator and a dimensionless damping ratio , for which Eq. 3 and Eq. 4 best fit the experimentally-measured vibration amplitude.

Figure 5: Vibration amplitude (a) and phase (b) averaged over the first cantilever, arrowed in Fig. 2(b), versus excitation frequency . Insets : retrieved vibration amplitude and phase maps in the neighborhood of the resonance, reported in Fig. 2(b-f). The points were obtained from a sequential measurement, the lines were obtained with the linear chirp (Eq. 16). The theoretical resonance lines in gray were determined from Eq. 3 and Eq. 4.

In conclusion, we performed optical path length modulation imaging by time-averaged heterodyne holography in off-axis and frequency-shifting conditions. To do so, modulation sidebands were acquired simultaneously through coherent frequency-division multiplexing. Additionally, a linear chirp of the excitation and the detection frequencies allowed to map the local mechanical phase retardation of a resonant cantilever with respect to the excitation signal without stroboscopy. The proposed method enabled robust and quantitative narrowband imaging of out-of plane vibration amplitude and phase.

We gratefully acknowledge support from Fondation Pierre-Gilles de Gennes (FPGG014), Agence Nationale de la Recherche (ANR-09-JCJC-0113, ANR-11-EMMA-046), région Île-de-France (C’Nano, AIMA), and the ”investments for the future” program (LabEx WIFI: ANR-10-LABX-24, ANR-10-IDEX-0001-02 PSL*).


  • (1) A. Bosseboeuf and S. Petitgrand. Characterization of the static and dynamic behaviour of m (o) ems by optical techniques: status and trends. Journal of micromechanics and microengineering, 13:S23, 2003.
  • (2) C. Rembe and A. Drabenstedt. Laser-scanning confocal vibrometer microscope: Theory and experiments. Rev. Scient. Instr., 77(8):083702–083702, 2006.
  • (3) K. Kokkonen and M. Kaivola. Scanning heterodyne laser interferometer for phase-sensitive absolute-amplitude measurements of surface vibrations. Appl. Phys. Lett., 92(6):063502–063502, 2008.
  • (4) S. Bramhavar, B. Pouet, and T. W. Murray. Superheterodyne detection of laser generated acoustic waves. Appl. Phys. Lett., 94(11):114102–114102, 2009.
  • (5) RL Whitman and A Korpel. Probing of acoustic surface perturbations by coherent light. Appl. Opt., 8(8):1567–1576, 1969.
  • (6) R.M. De La Rue, R.F. Humphryes, I.M. Mason, and E.A. Ash. Acoustic surface wave amplitude and phase measurements using laser probes. Proc. IEE., 119(2):117–126, 1972.
  • (7) G. Stegeman. Optical probing of surface waves and surface wave devices. IEEE Transactions on Sonics Ultrasonics, 23:33–63, 1976.
  • (8) Jean-Pierre Monchalin. Heterodyne interferometric laser probe to measure continuous ultrasonic displacements. Rev. Scient. Instr., 56(4):543–546, 1985.
  • (9) J.W. Wagner and J.B. Spicer. Theoretical noise-limited sensitivity of classical interferometry. J. Opt. Soc. Am. B, 4(8):1316–1326, 1987.
  • (10) D Royer and E Dieulesaint. Mesures optiques de déplacements d’amplitude 10^-4 à 10^ 2 angström. application aux ondes élastiques. Rev. Phys. Appl., 24(8):833–846, 1989.
  • (11) D Royer and E Dieulesaint. Optical probing of the mechanical impulse response of a transducer. Appl. Phys. Lett., 49(17):1056–1058, 1986.
  • (12) X Jia, A Boumiz, and G Quentin. Laser interferometric detection of ultrasonic waves propagating inside a transparent solid. Appl. Phys. Lett., 63(16):2192–2194, 1993.
  • (13) D. Royer and V. Kmetik. Measurement of piezoelectric constants using an optical heterodyne interferometer. Electronics Letters, 28(19):1828–1830, 1992.
  • (14) Akira Kimachi. Real-time heterodyne speckle pattern interferometry using the correlation image sensor. Applied optics, 49(35):6808–6815, 2010.
  • (15) Rikesh Patel, Samuel Achamfuo-Yeboah, Roger Light, and Matt Clark. Widefield heterodyne interferometry using a custom cmos modulated light camera. Optics Express, 19(24):24546–24556, 2011.
  • (16) R. L. Powell and K. A. Stetson. Interferometric vibration analysis by wavefront reconstruction. J. Opt. Soc. Am., 55:1593, 1965.
  • (17) C. C. Aleksoff. Time average holography extended. Appl. Phys. Lett, 14:23, 1969.
  • (18) Karl A. Stetson. Effects of beam modulation on fringe loci and localization in time-average hologram interferometry. J. Opt. Soc. Am., 60(10):1378, 1970.
  • (19) JA Levitt and KA Stetson. Mechanical vibrations: mapping their phase with hologram interferometry. Applied Optics, 15(1):195–199, 1976.
  • (20) M. Ueda, S. Miida, and T. Sato. Signal-to-noise ratio and smallest detectable vibration amplitude in frequency-translated holography: an analysis. Appl. Opt., 15(11):2690–2694, 1976.
  • (21) P. Picart, J. Leval, D. Mounier, and S. Gougeon. Time-averaged digital holography. Opt. Lett., 28:1900–1902, 2003.
  • (22) P. Psota, V. Ledl, R. Dolecek, J. Erhart, and V. Kopecky. Measurement of piezoelectric transformer vibrations by digital holography. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 59(9):1962–1968, 2012.
  • (23) Nicolas Verrier and Michael Atlan. Absolute measurement of small-amplitude vibrations by time-averaged heterodyne holography with a dual local oscillator. Opt. Lett., 38(5):739–741, 2013.
  • (24) Giancarlo Pedrini, Wolfgang Osten, and Mikhail E. Gusev. High-speed digital holographic interferometry for vibration measurement. Appl. Opt., 45(15):3456–3462, May 2006.
  • (25) C. Perez-Lopez, M.H. De la Torre-Ibarra, and F. Mendoza Santoyo. Very high speed cw digital holographic interferometry. Optics Express, 14(21):9709–9715, 2006.
  • (26) Ole J. Løkberg and Kåre Høgmoen. Vibration phase mapping using electronic speckle pattern interferometry. Appl. Opt., 15(11):2701–2704, Nov 1976.
  • (27) S. Petitgrand, R. Yahiaoui, K. Danaie, A. Bosseboeuf, and JP Gilles. 3d measurement of micromechanical devices vibration mode shapes with a stroboscopic interferometric microscope. Optics and lasers in engineering, 36(2):77–101, 2001.
  • (28) Julien Leval, Pascal Picart, Jean Pierre Boileau, and Jean Claude Pascal. Full-field vibrometry with digital fresnel holography. Appl. Opt., 44(27):5763–5772, Sep 2005.
  • (29) Nicolas Verrier, Michel Gross, and Michael Atlan. Phase-resolved heterodyne holographic vibrometry with a strobe local oscillator. Optics letters, 38(3):377–379, 2013.
  • (30) M. Atlan and M. Gross. Spatiotemporal heterodyne detection. Journal of the Optical Society of America A, 24(9):2701–2709, 2007.
  • (31) M. Paturzo, P. Memmolo, A. Tulino, A. Finizio, and P. Ferraro. Investigation of angular multiplexing and de-multiplexing of digital holograms recorded in microscope configuration. Opt. Express, 17(11):8709–8718, 2009.
  • (32) Tomohiro Kiire, Daisuke Barada, Jun ichiro Sugisaka, Yoshio Hayasaki, and Toyohiko Yatagai. Color digital holography using a single monochromatic imaging sensor. Opt. Lett., 37(15):3153–3155, Aug 2012.
  • (33) Tatsuki Tahara, Akifumi Maeda, Yasuhiro Awatsuji, Takashi Kakue, Peng Xia, Kenzo Nishio, Shogo Ura, Toshihiro Kubota, and Osamu Matoba. Single-shot dual-illumination phase unwrapping using a single wavelength. Opt. Lett., 37(19):4002–4004, Oct 2012.
  • (34) Etienne Cuche, Pierre Marquet, and Christian Depeursinge. Spatial filtering for zero-order and twin-image elimination in digital off-axis holography. Applied Optics, 39(23):4070, 2000.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description