Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Self-mixing interferometer based on sinusoidal phase modulating technique

Open Access Open Access

Abstract

A new self-mixing interferometry based on sinusoidal phase modulating technique is presented. Self-mixing interference occurs in the laser cavity by reflecting the light from a mirror-like target in front of the laser. Sinusoidal phase modulation of the beam is obtained by an electro-optic modulator (EOM) in the external cavity. The phase of the interference signal is calculated by Fourier analysis method. The interferometer is applied to measure the displacement of a high-precision commercial PZT with an accuracy of <10nm. The measurement range of the system mainly depends on the maximum operating frequency of EOM and the maximum sampling rate of A/D converter.

©2005 Optical Society of America

1. Introduction

Self-mixing interferometer (SMI) demonstrates unique features since it was first reported in the 1960’s to measure changes in the optical length (OPL) and the behavior of lasers [1]. These features include a simple, single-axis optical arrangement that requires minimal optical components and high sensitivity at low light level. SMI has been widely used for displacement, velocity measurement and ranging [24].

SMI was used to measure displacement with an accuracy of λ/2 by counting interference signal peaks. In order to increase the measurement accuracy beyond λ/2, some methods for analysis of SMI signal have been reported. J. Kato et al measured displacement by calculating the ratio of a phase shift and a period based on pseudo-heterodyne method [5]. The phase-locked technique was introduced in the SMI using a feedback loop [6]. The double external cavity method combined with FFT technique was reported to improve the displacement measurement accuracy [7]. Ovryn et al. used the phase shifting technique in the external cavity to analyze the SMI signal [8].

In this article a new self-mixing interferometry based on sinusoidal phase modulating technique is proposed. An EOM situated in the beam path between the front mirror of the laser and the external target is used to introduce phase modulation in the OPL. Considering that the characteristic of a He-Ne laser beam is quite better than that of a semiconductor laser and the polarization characteristic which is appropriate for EOM, we use a linearly polarized He-Ne laser as the coherent source. Fourier analysis method is proposed to extract phase from the modulated interference signal. Theoretical analysis and simulation results are presented. The optimum range of the modulation depth and the modulation frequency are discussed. Experimentally, the movement of a high-precision commercial PZT is reconstructed and a displacement measurement accuracy of <10nm can be obtained.

2. Experimental setup

 figure: Fig. 1.

Fig. 1. Experimental setup.

Download Full Size | PDF

The experimental setup is shown in Fig. 1. The light source is a 1mW linearly polarized He-Ne (632.8nm, TEM00) laser. The laser cavity consists of a plane mirror (front mirror) and a concave mirror (back mirror) which has a radius of 1m. Their reflectivities are 0.98 and 0.99 respectively. The reflectivity of the external target is 0.04. The length of the laser cavity and the external cavity are 180mm and 200mm, respectively. The beam emitted from the front mirror of the laser is sent through a various neutral density (ND) filter, an EOM (New Focus, 4002), a linear polarizer with a fixed orientation parallel to the polarization of the laser and then reflected back to the laser cavity by the external target. The interference signal is monitored by a photodiode (PD) situated behind the back mirror of the laser. The angle between the polarization direction of the laser and the electro-optically active axis of EOM is 00. Then the EOM can provide pure phase modulation with extremely low amplitude modulation. The Brewster window in the laser cavity and the linear polarizer in the external cavity make the interferometer more immune to a little change in the polarization of the laser beam caused by EOM. The external target is fixed on a high-precision commercial PZT (PI, P-841.10) which can obtain displacement accuracy of 0.15nm based on the close-loop design. The voltage signal from the PD is amplified and then digitized with a 200-KHZ, 12-bit analog-to-digital board on a PC bus (National Instrument, NI 6024E); this board is also used to send voltage to the PZT driver.

3. Theoretical analysis and simulation results

The basic theory of a He-Ne laser SMI can be explained by a three-mirror cavity model. The light emitted from the laser is reflected from the external target, a portion of the laser output back into the laser cavity and mix with the original light in the laser cavity forming the self-mixing effect. We assume that the laser’s frequency is unperturbed by the optical feedback. This assumption is appropriate for gas lasers with symmetric gain and dispersion curves. The output intensity of a He-Ne laser SMI can be expressed as [89]

I=I0[1+mcos(ϕ)j=0(η)jcos(jϕ)]

Where ϕ=4πL/λ, m=/l, I 0 is the power of the laser output without feedback, η is the parameter denoting coupling effect from the external cavity back into the laser cavity, l is the length of the laser cavity, L is the length of the external cavity, λ is the wavelength of the laser, K is a constant relating to the laser’s operating parameters. Also included in Eq. (1) is the effect of multiple reflections that can occur between the external target and the laser’s output mirror. Clearly, m/η has no relation with the optical feedback.

Assuming that SMI operates in a low-feedback regime where the effect of multiple reflections in the external cavity is small (i.e., η is small), Eq. (1) can be expressed as

I=I0[1+mcos(ϕ)]

This assumption forms the basis of our subsequent research. Error caused by multiple reflections will be considered later. Using this simple configuration, the displacement of the external target can be measured by counting interference signal peaks. To achieve higher accuracy, an experimentally controlled additive phase ψ(t)=asin(2πfmt+β) is introduced by EOM situated in the external cavity, where a is the modulation depth, fm is the modulation frequency, β is the initial phase of the modulation. Considering that the beam pass through the EOM twice in the external cavity, the modulated interference signal can be written as

I(t)=I0{1+mcos[ϕ(t)+2asin(2πfmt+β)]}

Expanding Eq. (3) in a Fourier series, harmonics at frequency fm and 2fm have following expressions:

I(fm,t)=2mI0sinϕ(t)J1(2a)sin(2πfmt+β)
=A1(t)sin(2πfmt+β)
I(2fm,t)=2mI0cosϕ(t)J2(2a)cos(4πfmt+2β)
=A2(t)cos(4πfmt+2β)

Where Jn(2a) is the nth-order Bessel function. Then we can calculate the phase ϕ(t) from A 1(t) and A 2(t)using the subsequent relation:

ϕ(t)=arctan[A1(t)A2(t)·J2(2a)J1(2a)]

The Fourier analysis method is proposed to determine the values of A 1(t) and A 2(t). It mainly takes following steps: (1) taking the Fourier transform of the interference signal; (2) selecting the first harmonic I(fm) and second harmonic components I(2fm) of the Fourier spectra with two windows: fm/2<f<3fm/2 and 3fm/2<f<5fm/2 ; (3) taking the inverse Fourier transform of the filtered component, represented by Ifm(t) and I2fm(t) respectively; and (4) computing A 1(t) and A 2(t) using the subsequent relation:

A1(t)=imag{Ifm(t)ej(2πfmt+β)}
A2(t)=real{I2fm(t)ej(4πfmt+2β)}

As for the modulation depth, a question which must be taken into account is the null points of the Bessel function. When J 1(2a) or J 2(2a) is zero, we can not use Eq. (6) to calculate the phase ϕ(t). As we all know, J 1(2a) and J 2(2a) are oscillating functions of parameter a, with a series of zero-crossings whose values are listed in Table 1. These values of a should be avoided when choosing the modulation depth of EOM.

Tables Icon

Table 1. Values of a which null J1(2a) or J2(2a) (represented by X)

As the modulation frequency is concerned, the occurrence of measurement error due to band overlapping in the frequency domain should be considered. A 1(t) and A 2(t) have instantaneous frequency F=(t)/(2πdt) which is proportional to the speed of the external target. Therefore, the first harmonic and the second harmonic centered at fm and 2fm present a spectral width proportional to the maximum velocity of the external target. To avoid overlapping problems, ϕ(t) can only have spectral components F limit to fm/2 :

Ffm2

For example, if the external target has a motion d=d 0 sin(2πft), then ϕ(t)=4π[L+d]/λ. The instantaneous frequency F of A 1(t) and A 2(t)is

F=dϕ(t)(2πdt)=4πfd0cos(2πft)λ

In order to satisfy the relationship in Eq. (9) all the time, the following relation should be satisfied:

Fmax=4πfd0λfm2fm8πfd0λ

From Eq. (11) we can see that the minimum modulation frequency needed to demodulate the phase is directly proportional to the product of the amplitude and the frequency of the external target. So the measurement range of the system has no theoretical limit if the EOM has an enough high modulation frequency.

The phase ϕ obtained using the demodulation technique proposed is wrapped within the region of -π/2 and π/2. After a phase unwrapping process and following the relationship between the phase ϕ and the length of the external cavity, the movement of the external target can be reconstructed.

Assuming that the external target move at a sinusoidal form with amplitude 1µ m and frequency 20HZ, the phase retrieved process is simulated in Fig. 2. The modulation depth of EOM is 1.1rad. The modulation frequency is 2KHZ. The output intensity of the laser is 1mW in the absence of feedback and the fringe visibility m is 0.08. The phase change due to the motion of the external target is demonstrated in Fig. 2(a). The modulated interference signal is shown in Fig. 2(b). The Fourier spectra of the SMI signal is shown in Fig. 2(c). Following the phase calculation steps, the phase is demodulated as in Fig. 2(d). After an unwrapping process, the continuous phase is obtained as in Fig. 2(e). The phase retrieved error due to calculation is shown in Fig. 2(f). The max error is 2.5×10-13 rad in Fig. 2(f) and the standard deviation is 6.3×10-14 rad.

 figure: Fig. 2.

Fig. 2. (a) Simulated phase change due to the motion of the external target. (b) Simulated modulated interference signal. (c) Fourier spectra of the interference signal. (d) Calculated phase (wrapping). (e) Retrieved phase (unwrapped). (f) Phase retrieved error.

Download Full Size | PDF

When the effect of multiple reflections can not be neglected completely, the additional terms in Eq. (1) (i.e., j=1, 2 …) will cause phase retrieved error. Since the interferometer is restricted in a low feedback regime during measurement, we only consider the phase retrieved error caused by the second feedback in the external cavity (i.e., error caused by additional term j=1). Figure 3 shows the phase retrieved error due to the second feedback in the external cavity corresponding to Fig. 2 in the condition of m/η=0.5. The maximum error is 0.05rad and the standard deviation is 0.03rad. Figure 4 shows the maximum error and the standard deviation caused by j=1 with respect to m in the condition of m/η=0.5. Other simulation parameters are the same as in Fig. 2.

 figure: Fig. 3.

Fig. 3. Phase retrieved error due to the second feedback.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. Error caused by second feedback with respect to m.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. Standard deviation caused by second feedback versus the modulation depth a.

Download Full Size | PDF

For two reflections in the external cavity, the modulated interference signal is given by

I(t)=I0{1+mcos[ϕ(t)+2asin(2πfmt+β)]mηcos2[ϕ(t)+2asin(2πfmt+β)]}

Then the demodulated phase ϕ(t) j=0,1 is related to the actual phase ϕ(t) by

tan(ϕ(t)j=0,1)=A1(t)A2(t)·J2(2a)J1(2a)=1ηcosϕ(t)J1(4a)J1(2a)1ηcosϕ(t)J2(4a)J2(2a)+η2cosϕ(t)J2(4a)J2(2a)tanϕ(t)

From Eq. (13) we can see that the modulation depth of the EOM influence the phase retrieved error due to second feedback in the external cavity. So the appropriate modulation depth a can be selected to decrease the influence of the multiple reflections. Here the modulation depth is limited in the region of a<π considering the operating characteristic of the EOM we used. Figure 5 shows the standard deviation caused by the second feedback in the external cavity versus modulation depth a in the condition of m/η=0.5. Other simulation parameters are the same as in Fig. 2. It is obvious that the appropriate range of the modulation depth is a∈[1.1rad,1.8rad]∪[2.8rad, 3.14rad]. For simulation of other parameters (not reported here for brevity), the same range can be obtained.

4. Experimental results

Experiments have been done to confirm the validity of the sinusoidal phase modulating self-mixing interferometer and the phase demodulation method. In our experiment the first rising edge of TTL level from the signal generator which drives the EOM driver is introduced to the A/D card as the trigger of data acquisition. Then a sinusoidal phase modulation with the initial phase π/2 can be obtained.

 figure: Fig. 6.

Fig. 6. (a) The obtained SMI signal. (b) Fourier spectra of the SMI signal. (c) Extracted phase (wrapping). (d) Displacement reconstruction result.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. Displacement reconstruction result of PZT with frequency 40HZ and amplitude of 4µ m (peak to peak).

Download Full Size | PDF

First, the PZT is controlled to move at a sinusoidal form with frequency10HZ and amplitude 2µm (peak to peak). The modulation frequency is1KHZ and the modulation depth is 1.2rad. The sampling rate of the A/D card is 50KHZ. Figure 6(a) is the interference signal obtained by A/D card. Figure 6(b) is the Fourier spectra of the detected signal in which the first harmonic component centered at 1KHZ and the second harmonic component centered at 2KHZ. Following the phase extraction steps, the phase is retrieved as in Fig. 6(c). After unwrapping the phase and using the relationship between the phase ϕ and the length of the external cavity, the movement of PZT is reconstructed as in Fig. 6(d). In displacement measurement, only the relative value of the two positions is important. So, we examine the peak-to-peak amplitude in Fig. 6(d) and get the displacement measurement accuracy. The peak-to-peak amplitude of the reconstructed waveform in Fig. 6(d) is 2.002µm. Second the PZT is controlled to move at a sinusoidal form with frequency 40HZ and amplitude 4µm (peak to peak). The modulation frequency is 6KHZ and the modulation depth is1.2rad. The sampling rate of A/D card is 100KHZ. The reconstruction result is in Fig. 7 and the max error for displacement measurement is 7nm.

5. Conclusion

The numerical simulation and experimental results show that the proposed new method can effectively demodulate the phase of the interference signal with high accuracy. The combination of modulation and demodulation decreases the sensitivity of the instrument to fluctuations of the laser power and the noise induced by environment. Since the minimum modulation frequency needed to demodulate the phase is directly proportional to the product of the amplitude and the frequency of the external target, the measurement range of the setup mainly depends on the maximum operating frequency of EOM and the maximum sampling rate of A/D converter. Considering that EOM can operate at a frequency up to 100MHZ, the system presented in this paper has a lager measurement range regarding the frequency and the amplitude of the motion of the measured object than the traditional SMI. In our experiments, we only reconstructed the movement of PZT with amplitude (peak to peak) up to 8µm and frequency up to 60HZ due to the limitation of the PZT used and a displacement measurement precision of <10 nm can be obtained. So the sinusoidal phase modulating SMI offers both the resolution and measurement range increase in the self-mixing displacement sensor. Besides the application of displacement measurement presented in this paper, the sinusoidal phase modulating SMI can also be used to microvibration measurement and surface profile measurement. The related works will be reported in the future.

Acknowledgments

This work was supported by the National Natural Science Foundation of China, No. 50375074.

References and links

1. Th. H. Peek, P. T. Bolwjin, and C. Th. Alkemade, “Axial mode number of gas lasers from moving-mirror experiments,” Am. J. Phys. 35, 820–831 (1967). [CrossRef]  

2. T. Yoshino, M. Nara, and S. Mnatzakanian, “Laser diode feedback interferometer for stabilization and displacement measurements,” Appl Opt. 26, 872–877 (1987). [CrossRef]  

3. S. Shinohara, A. Mochizuki, and H. Yoshida et al. “Laser Doppler velocimeter using the self-mixing effect of a semiconductor laser diode,” Appl Opt. 25, 1417–1419 (1986). [CrossRef]   [PubMed]  

4. S. Shinohara, H. Yoshida, H. Ikeda, K. Nishide, and M. Sumi, “Compact and high-precision range finder with wide dynamic range and its application,” IEEE Trans Instrum. Meas. 41, 40–44 (1992). [CrossRef]  

5. J Kato, N Kikuchi, I Yamaguchi, and S. Ozono “Optical feedback displacement sensor using a laser diode and its performance improvement,” Meas. Sci. Technol. 6, 45–52 (1995). [CrossRef]  

6. T. Suzuki, S. Hirabayashi, and O. Sasaki et al “Self-mixing type of phase-locked laser diode interferometry,” Opt. Eng. 38, 543–548 (1999). [CrossRef]  

7. M. Wang and Guangming Lai, “Displacement measurement based on Fourier transform method with external cavity modulation,” Rev. Sci. Instrum. 72, 3440–3445 (2001). [CrossRef]  

8. B. Ovryn and J. H. Andrews, “Phase-shifted laser feedback interferometry,” Opt. Lett. 23, 1078–1080 (1998). [CrossRef]  

9. G. Liu and S. L. Zhang et al, “Theoretical and experimental study of intensity branch phenomena in self-mixing interference in a He-Ne laser,” Opt. Commun. 221, 387–393 (2003). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. Experimental setup.
Fig. 2.
Fig. 2. (a) Simulated phase change due to the motion of the external target. (b) Simulated modulated interference signal. (c) Fourier spectra of the interference signal. (d) Calculated phase (wrapping). (e) Retrieved phase (unwrapped). (f) Phase retrieved error.
Fig. 3.
Fig. 3. Phase retrieved error due to the second feedback.
Fig. 4.
Fig. 4. Error caused by second feedback with respect to m.
Fig. 5.
Fig. 5. Standard deviation caused by second feedback versus the modulation depth a.
Fig. 6.
Fig. 6. (a) The obtained SMI signal. (b) Fourier spectra of the SMI signal. (c) Extracted phase (wrapping). (d) Displacement reconstruction result.
Fig. 7.
Fig. 7. Displacement reconstruction result of PZT with frequency 40HZ and amplitude of 4µ m (peak to peak).

Tables (1)

Tables Icon

Table 1. Values of a which null J1(2a) or J2(2a) (represented by X)

Equations (15)

Equations on this page are rendered with MathJax. Learn more.

I = I 0 [ 1 + m cos ( ϕ ) j = 0 ( η ) j cos ( j ϕ ) ]
I = I 0 [ 1 + m cos ( ϕ ) ]
I ( t ) = I 0 { 1 + m cos [ ϕ ( t ) + 2 a sin ( 2 π f m t + β ) ] }
I ( f m , t ) = 2 m I 0 sin ϕ ( t ) J 1 ( 2 a ) sin ( 2 π f m t + β )
= A 1 ( t ) sin ( 2 π f m t + β )
I ( 2 f m , t ) = 2 m I 0 cos ϕ ( t ) J 2 ( 2 a ) cos ( 4 π f m t + 2 β )
= A 2 ( t ) cos ( 4 π f m t + 2 β )
ϕ ( t ) = arc tan [ A 1 ( t ) A 2 ( t ) · J 2 ( 2 a ) J 1 ( 2 a ) ]
A 1 ( t ) = imag { I f m ( t ) e j ( 2 π f m t + β ) }
A 2 ( t ) = real { I 2 f m ( t ) e j ( 4 π f m t + 2 β ) }
F f m 2
F = d ϕ ( t ) ( 2 π d t ) = 4 π f d 0 cos ( 2 π f t ) λ
F max = 4 π f d 0 λ f m 2 f m 8 π f d 0 λ
I ( t ) = I 0 { 1 + m cos [ ϕ ( t ) + 2 a sin ( 2 π f m t + β ) ] m η cos 2 [ ϕ ( t ) + 2 a sin ( 2 π f m t + β ) ] }
tan ( ϕ ( t ) j = 0 , 1 ) = A 1 ( t ) A 2 ( t ) · J 2 ( 2 a ) J 1 ( 2 a ) = 1 η cos ϕ ( t ) J 1 ( 4 a ) J 1 ( 2 a ) 1 η cos ϕ ( t ) J 2 ( 4 a ) J 2 ( 2 a ) + η 2 cos ϕ ( t ) J 2 ( 4 a ) J 2 ( 2 a ) tan ϕ ( t )
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.