## Abstract

An optical feedback laser diode (OFLD) operating in period-one oscillation (POO) with a moving external target is investigated by exploring its potential sensing capability. First, the modeling of an OFLD-POO sensing system is presented. An analytical expression is derived for OFLD-POO sensing signal, from which a new displacement measurement method is developed. The proposed sensing model is verified by the well-known Lang-Kobayashi equations used to describe the dynamic behavior of a laser with optical feedback. Then, an experimental OFLD-POO system is built in order to demonstrate an application example for displacement sensing. The measurement results show that the OFLD-POO sensing system can achieve displacement measurement with large measurement range, high sensitivity, and resolution.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

Optical feedback interferometry (OFI) or Self-mixing interferometry (SMI) is a promising non-contact sensing technology, which uses the self-mixing effect that occurs when a fraction of the light back-reflected or back-scattered by an external target reenters the laser cavity [1,2]. In this case, the steady-state intensity of the lasing light is modulated by the external optical feedback. A typical OFI system consists of a laser diode (LD), a photodiode (PD) packaged at the rear of LD, a lens and an external target [3–6]. We call the laser diode in this case as optical feedback laser diode (OFLD) and the modulated intensity as OFLD signal. As a preferred, minimum part-count scheme useful for engineering implementation, various OFLD-based applications have been developed in the industrial and laboratory environment, such as measurement of displacement [7,8], vibration [9], velocity [10], alpha factor [11,12], Young’s modulus [13], strain [14], acoustic emission [15], acoustic field imaging [16], etc.

Most of the OFLD-based applications are based on the stable OFI mathematic model, which is derived from the steady-state solution of the Lang and Kobayashi (L-K) equations [17], or from the classical three-mirror model [2], by assuming that the OFLD system operates in stable mode, i.e., both the electric field and carrier density in an LD with a stationary external cavity can reach a constant state after a transient period. However, many works [18–20] have demonstrated that an OFLD system with a stationary external target may be unstable under some certain operation conditions, such as in periodic, quasi-periodic or chaotic oscillation state. We studied OFLD signals with a moving target in different conditions and our recent works [21,22] shows that such sensing system may produce different types of OFLD signals which are classified into three regions defined in [21], i.e. stable, semi-stable and unstable region under different feedback level factors C. C is defined as: $C\text{=}\kappa \tau \sqrt{1\text{+}{\alpha}^{2}}/{\tau}_{in}$, where $\kappa $is the optical feedback strength, $\tau $and ${\tau}_{in}$ is the laser round-trip time in the external and internal cavity respectively, and $\alpha $is the linewidth enhancement factor [2]. In stable region, for a stationary target, the laser intensity will be constant after the transient state, whereas in semi-stable and unstable region, relaxation oscillation will be undamped and the laser may enter period-one, quasiperiodic or chaotic oscillation under different feedback strength. The work [21,22] discovers that an OFLD can easily and sometimes inevitably enter the period-one oscillation, e.g. for the case when the feedback level factor C is around 2.5, which is quite commonly reported for the moderate feedback regime in the literature. In this case, the OFLD signals (laser intensity) generated by an OFLD system cannot be described by the existing OFI model. The laser intensity is modulated in a quite complicated form and contains very high-frequency components. In order to differentiate the conventional OFLD signals, we name the laser intensity signals generated by the OFLD system operating in period-one oscillation (POO) as OFLD-POO signals. The work in [23] also depicts some high-frequency components in the sawtooth-like OFLD signals, but the LD used in [23] still operates in conventional self-mixing conditions rather than in POO state. Usually, the PD packaged inside an LD is used to detect the laser intensity, but it has limited bandwidth less than 1GHz. However, the frequency components contained in an OFLD-POO signal can even higher than several GHz [21,22]. Hence, using existing OFI configuration for detection of OFLD-POO signal, many frequency components are missing. The work in [22] has experimentally revealed the features of such signals and indicated the possible sensing application but it did not make analytic expression on the signal waveform. As an OFLD-POO signal contains rich information associated with the laser and its external target, it is of great interest and significance to investigate the signal features by modelling the signal waveform, from which we can develop a suitable signal processing algorithm to retrieve the useful information contained in the signals.

In this work, an OFLD operating in POO with a moving external target is investigated by exploring its potential sensing capability. Starting from the L-K equations, we firstly derive an approximate analytical model to describe the OFLD-POO signals. The numerical simulations for original L-K equations are then conducted to verify the derived analytical expression. Finally, a displacement sensing method by simultaneously using the time-domain OFLD-POO waveform and relaxation oscillation (RO) frequency contained in this signal is proposed. From both simulation and experiment, the feasibility of proposed displacement sensing method has been verified.

## 2. Modelling and analysis

#### 2.1 Theoretical modelling

The behavior of an LD with optical feedback can be described by the well-known L-K equations [17], as shown in Eqs. (1)-(3). Here, $E(t)$ is the amplitude of the electric filed, $\varphi (t)$ is the phase of the electric field, and $N(t)$ is the carrier density. ${E}^{2}(t)$ is usually considered as the laser intensity or power, denoted by $P(t)={E}^{2}(t)$.

Starting from the L-K equations, we aim to get an approximate analytical expression for OFLD-POO signal (i.e.,${E}^{2}(t)$). Based on the numerical simulation of L-K equations and the work in [18,20], when an OFLD operates in POO, $E(t)$, $N(t)$ and $\varphi (t)$ can be written in the form of below:

Regarding$\overline{E}$and$\overline{N}$, they can be determined by solving Eqs. (1)-(3) through setting $dE(t)/t=0,dN(t)/t=0,d\varphi (t)/dt={\omega}_{s}-{\omega}_{0}$ [18], as shown below:

Substituting $E(t)$ and$N(t)$ in Eq. (3) by Eq. (4) and Eq. (5), and neglecting the terms with$\Delta E\cdot \Delta N$, $\Delta {E}^{2}$ and $\Delta {N}^{2}$(as $\Delta E/\overline{E}<<1$ and $\Delta N/\overline{N}<<1$) [18], we have:

^{th}-order and 1

^{st}-order Bessel function respectively. Note that higher harmonic components of ${\omega}_{r}$ are neglected [18]. Then, we have:

As Eq. (18) is an identical equation as well, letting the corresponding terms in both sides are equal, we have:

Therefore, Eq. (34) is the laser intensity when the OFLD system is in period-one oscillation. When the external cavity is moving, it will cause a varying feedback phase ${\varphi}_{F}$ (${\varphi}_{F}\text{=}{\omega}_{s}\tau =\text{2}{\omega}_{s}L/c$, where c is the vacuum light speed) [2,11–13] as well as ${\omega}_{r}$ [24]. We re-write Eq. (34) as below to include displacement sensing:

Equation (36) is the conventional OFLD signal reported in the literature [2,11,12]. This signal is in the form of sinusoidal (at weak feedback case) or sawtooth-like (moderate or above feedback level) waveform when the external cavity length $L$changes. Each periodical variation (called a fringe) corresponds to a displacement with half laser wavelength (${\lambda}_{0}/2$) [2]. For a varying external target, the external cavity length can be expressed as $L\text{=}{L}_{0}\text{+}\Delta L$, where ${L}_{0}$ is the initial external cavity length, and $\Delta L$ is the displacement.

We are more interested in the last term appeared in Eq. (35). Since this term also contains the information of displacement $\Delta L$, we can make use of it to achieve displacement sensing. The RO frequency of an OFLD has been investigated and an approximate analytical expression for RO frequency has been derived for an OFLD is in steady state [20]. The analysis on the RO frequency and its application for measuring linewidth enhancement factor of the LD has been discussed in [26]. When the OFLD operates in POO, the expression is invalid and the RO frequency is different from that of the OFLD in steady state. In this case, the frequency of POO can be treated as the RO frequency, i.e. ${\omega}_{r}$ contained in Eqs. (35) and (37), as in our previous work [24], which shows that ${\omega}_{r}$changes with displacement $\Delta L$ in the form of a saw-tooth waveform with a period of ${\lambda}_{0}/2$ when the OFLD system is in period-one oscillation. Thus, we have ${\omega}_{r}(\Delta L)\text{=}{\omega}_{r}(\Delta L\text{+}{\lambda}_{0}/2)$. Denoting the maximum and minimum value of ${\omega}_{r}$within a period by ${\omega}_{r\mathrm{max}}$ and ${\omega}_{r\mathrm{min}}$respectively, the expression for ${\omega}_{r}$within one period ($0<\Delta L<{\lambda}_{0}/2$) can be expressed as:

It can be seen that $\Delta {P}_{OFLD-POO-Ep}(L)$ also periodically change with displacement $\Delta L$ by${\lambda}_{0}/2$ . Figure 1(a) shows an example of the relationship between the angular RO frequency and displacement of external target when ${L}_{0}\text{=}24cm$, $I=1.3{I}_{th}$ and $C=4.5$, where ${I}_{th}$is the threshold injection current. Figure 1(b) shows the waveform of $\Delta {P}_{OFLD-POO-Ep}(L)$ for a varying external cavity with the same operation conditions.

In the following, we only consider the last term in Eq. (35), which is called OFLD-POO signal denoted by$\Delta {P}_{OFLD-POO}(L,t)$, and re-write it as below:

Now, we can describe that the amplitude and frequency of an OFLD-POO signal is modulated by the displacement to be measured. Both the envelope ($\Delta {P}_{OFLD-POO-Ep}(L)$) and frequency of the signal contain displacement information. So we can explore a new sensing approach by using this signal. We denote the maximum and minimum values of $\Delta {P}_{OFLD-POO-Ep}$ as $\Delta {P}_{OFLD-POO-Ep-\mathrm{max}}$ and$\Delta {P}_{OFLD-POO-Ep-\mathrm{min}}$, respectively. Based on the simulation and the work in [24], ${\omega}_{r\text{-}offset}$ and ${\omega}_{r\mathrm{max}}-{\omega}_{r\mathrm{min}}$ is much smaller than ${\omega}_{r0}$ and we have${\omega}_{r}\approx {\omega}_{r0}$ . In the case of $\kappa <0.01$, we get:

Based on the simulation and the work in [24], ${\omega}_{r\text{-}offset}$and ${\omega}_{r\mathrm{max}}-{\omega}_{r\mathrm{min}}$ are usually in the order of ${10}^{8}rad/s$. e.g. in the case in Fig. 1, ${\omega}_{r\text{-}offset}\text{+}{\omega}_{r\mathrm{max}}-{\omega}_{r\mathrm{min}}\approx 6\times {10}^{8}rad/s$. Substituting the typical values in Table1 into Eq. (44), i.e.${\tau}_{in}=8.0\times {10}^{-12}s$, ${\tau}_{p}=2.0\times {10}^{-12}s$, $\alpha \text{=}3$${\tau}_{s}=2.0\times {10}^{-9}s$ and $\kappa \text{=}0.007$ in the case in Fig. 1, we get $\Delta {P}_{OFLD-POO-Ep-pp}/\Delta {P}_{OFLD-pp}\approx 154$ . This reveals that the envelope is much more sensitive to the displacement compared to using conventional OFLD signal, which motivates the use of the envelop for displacement detection.

Since the frequency of the OFLD-POO signal is also modulated by the displacement, we can detect a micro-displacement within ${\lambda}_{0}/2$ through frequency variation. This enables us to achieve very high sensing resolution for displacement. Regarding this point, our work in [24] has demonstrated that the sensing resolution can be as high as ${\lambda}_{0}/1280$ when the spectrum analyzer has a resolution bandwidth of 62.5 KHz. If we use envelop for large range displacement by fringe-counting with a resolution of ${\lambda}_{0}/2$ and cooperate it with frequency sensing, a new displacement sensor can be designed with large measurement range, high sensitivity and resolution. The restriction on the speed of the target movement (denoted by ν) depends on the value of the RO frequency. Each OFLD-POO fringe should at least cover a few RO periods. Thus we can evaluate that the speed meet: $\nu <{\lambda}_{0}{f}_{RO}/2K$, where $K$ is number of the RO periods within one OFLD-POO fringe. Based on our simulation, to get a clear envelop of OFLD-POO signal, we need to have $K>20$. Theoretically, the dynamic measurement range is the maximum measurable displacement that guarantees OFLD always operating in POO during the displacement measurement. This limits the displacement measurement range. The work in [27] has investigated how POO region change with external target position in a certain injection current, from which we can estimate the dynamic displacement range is around several centimeters.

#### 2.2 Verification for the proposed sensing model

Figure 2 shows the simulation results by numerically solving the original L-K Eqs. (1)-(3) and the approximate analytical expression of Eq. (40) when the external cavity length has a linear displacement of $3{\lambda}_{0}$, the initial external cavity length ${L}_{0}\text{=}24cm$, and the injection current $I=1.3{I}_{th}$. Note that ${\omega}_{r}\tau \approx 23\pi $ in this case, which satisfies the requirement for Eq. (40). Figure 2(a) is the displacement; Figs. 2(b) and 2(c) show the numerically simulation results of L-K equations with $C=2.8$ and $C=4.5$; Figs. 2(d) and 2(e) are the simulation results of the Eq. (40) with the same operation conditions as those in Figs. 2(b) and 2(c). For the simulation of Eq. (40), we firstly obtain the relationship between ${\omega}_{r}$and $\Delta L$by numerically solving the L-K equations as shown in Fig. 1(a). From Fig. 2, it can be concluded that the simulation results from the L-K equations and the approximate analytical expression match closely, demonstrating that the approximate expression is valid when the OFLD is in period-one oscillation. Note that, in Figs. 2(b) and 2(d), there are regions between the fringes where the OFLD is stable. Therefore, it is essential to ensure that the OFLD system is always in period-one oscillation with the change of external cavity length when using RO frequency for displacement sensing, as shown in Figs. 2(c) and 2(e). To guarantee the robustness of the so-called OFLD-POO regime, we need to choose suitable injection current the initial external cavity length. The method of how to choose them is detailed in our previous work [24].

## 3. Displacement sensing

As discussed in Section 2, the displacement of the external target can be retrieved by using both time-domain waveform (OFLD-POO signal) and the RO frequency contained in the signal. Let us express this displacement to be measured as below:

where $\Delta {L}_{t}$ is measured by fringe counting with ${\lambda}_{0}/2$ resolution by using the envelop in OFLD-POO signal, and $\Delta {L}_{f}$is measured by the RO frequency for the fractional displacement within ${\lambda}_{0}/2$. As an example, Fig. 3 shows an OFLD-POO signal obtained by Eq. (40) as in Fig. 2(e) and the corresponding displacement information is contained in it.It can be seen that an OFLD-POO signal contains a few complete fringes and two incomplete fringes (or called fractional fringes) at the beginning and end part of the waveform. The measurement of the fractional fringe enable a displacement measurement with resolution higher than ${\lambda}_{0}/2$. As described in Fig. 1, there is a certain linear relationship between ${\omega}_{r}$ and the displacement within each fringe. As in Eq. (38), the maximum and minimum RO frequency within ${\lambda}_{0}/2$ are ${\omega}_{r\mathrm{max}}$ and ${\omega}_{r\mathrm{min}}$.We denote the RO angular frequency at the start and end point in a fractional fringe by ${\omega}_{r1}$ and ${\omega}_{r2}$respectively. Then the displacement corresponding to fractional fringes can be obtained as:

A detailed procedure of retrieving the displacement from an OFLD-POO signal is presented as below. Assuming that the external target has a linear movement with $\Delta L\text{=}3{\lambda}_{0}$ as shown in Fig. 4(a), the corresponding OFLD-POO signal is shown in Fig. 4(b) which can be obtained by Eq. (40). The upper envelope of the OFLD-POO signal can be retrieved by upper-peak detection shown in Fig. 4(c). Then, applying differentiation on Fig. 4 (c), we can get pulse train shown in Fig. 4(d). By counting the number of pulses, we can obtain $\Delta {L}_{t}\text{=(}N-1)\cdot {\lambda}_{0}/2$ (here *N* is the pulse train number). It can be seen from Fig. 4(d) that there are 6 pulses. Thus $\Delta {L}_{t}=2.5{\lambda}_{0}$. Then we need to obtain the fractional displacement determined by the first fractional fringe and last fractional fringe denoted by $\Delta {L}_{f1}$and $\Delta {L}_{f2}$ as shown in Fig. 4(b). In order to get $\Delta {L}_{f1}$and $\Delta {L}_{f2}$, the RO angular frequency at the start and end of the displacement should be obtained firstly. Before the displacement starts and after the displacement ends, for the stationary target, we get the laser intensity and then the RO angular frequency by applying FFT on the intensity as done in Fig. 1(a). Using this method, we get that the RO angular frequencies at these two time points are the same, which are$1.399\times {10}^{10}rad/s$. According to Fig. 1(a), ${\omega}_{r\mathrm{max}}\text{=}1.436\times {10}^{10}rad/s$ and ${\omega}_{r\mathrm{min}}\text{=}1.394\times {10}^{10}rad/s$. Hence, based on Eq. (42), for $\Delta {L}_{f1}$, ${\omega}_{r2}\text{=}1.436\times {10}^{10}rad/s$, and ${\omega}_{r1}\text{=}1.399\times {10}^{10}rad/s$, whereas for $\Delta {L}_{f2}$, ${\omega}_{r2}\text{=}1.399\times {10}^{10}rad/s$, and ${\omega}_{r1}\text{=}1.394\times {10}^{10}rad/s$ . Accordingly, we can get $\Delta {L}_{f1}\text{=}0.441{\lambda}_{0}$, and$\Delta {L}_{f2}\text{=}0.059{\lambda}_{0}$. Thus, $\Delta L\text{=}\Delta {L}_{t}+\Delta {L}_{f1}\text{+}\Delta {L}_{f2}=3{\lambda}_{0}$, which is the same as the pre-set value.

In the following, an experimental system is built and an application example of displacement sensing by using the proposed method is presented. Figure 5 shows the experimental system. The laser source is a single mode laser diode (Hitachi, HL8325G) with a wavelength of 830 nm, maximum output power of 40 mW, and threshold current of 35 mA. The laser diode is driven and temperature-stabilized by an LD controller (Thorlabs, ITC4001). The temperature of the LD is stabilized at 21 ± 0.01${}^{o}C$. A mirror is used as the external cavity which is affixed on the surface of a piezoelectric transducer (PZT) (PI, P-841.20) in order to provide sufficient optical feedback to the LD. The PZT actuator with an integrated displacement sensor is used to generate the displacement, and it has a displacement resolution of 9 nm with a traveling range up to 30 μm, which is driven by a PZT controller (PI, E-625). A high-speed digital oscilloscope (Tektronix DSA70804) with sampling rate 25 GS/s is used to observe OFLD-POO signals. An attenuator is used to adjust the optical feedback level of the OFLD system. As the OFLD-POO signals contain high-frequency components, high speed detection devices are needed. A beam splitter (BS) with splitting ratio of 50:50 is used to direct a part of the self-mixing light into the external fast PD, with a bandwidth of 9.5GHz, through the fiber port coupler. The detected OFLD/OFLD-POO signals are then collected and recorded by the digital oscilloscope.

To make the system operate in period-one oscillation, the initial external cavity length and injection current should be properly chosen, which is referred to the work in [24]. We choose $I=50mA$,and ${L}_{0}=18cm$ in this work. Before staring the measurement, we firstly need to get the relationship between RO frequency (${\omega}_{r}$) and external target displacement ($\Delta L$) when the external cavity length changing by ${\lambda}_{0}/2$.

we apply 30 nm displacement on the PZT each time, then use the Tektronix DSA 70804 oscilloscope to record the laser intensity waveform and thus obtain the corresponding peak frequency (denoted by ${f}_{RO}={\omega}_{r}/2\pi $) by using the FFT function provided by the oscilloscope [24]. Figure 6 shows an example of the laser intensity waveform and the corresponding spectrum. Thus, in each half wavelength, we will get about 13 points. Figure 7 illustrates the experimental results for the relationship between ${f}_{RO}$and displacement $\Delta L$.

Then we apply a dynamic linear displacement on the PZT by applying a controlling voltage through the PZT driver. Figure 8 (a) is the PZT controlling signal applied on the PZT, Fig. 8(b) is the corresponding OFLD-POO signal. The displacement generated by the PZT (PI, P-841.20) can be read from its internal integrated sensor which indicates that the displacement is $\Delta L\text{=}2304nm$.

The experimental OFLD-POO signal in Fig. 8(b) is processed according to the procedure described in Fig. 4 to get the displacement. 5 complete fringes can be found which corresponds to $\Delta {L}_{t}=5{\lambda}_{0}/2$. In order to obtain the displacement of the first and last fractional fringes, we need to obtain the RO frequency at the start and end points. By using the method in Fig. 6, the RO frequency at these two points are measured as $2.634GHz$ and $2.595GHz$ respectively. From Fig. 7, it can be seen the variation range of ${f}_{RO}$ in each period is from $2.579GHz$ to $2.645GHz$. Thus, we can get $\Delta {L}_{f1}\text{=}0.083{\lambda}_{0}$$$and $\Delta {L}_{f2}\text{=}0.121{\lambda}_{0}$ based on Eq. (46). So $\Delta L\text{=}\Delta {L}_{t}+\Delta {L}_{f1}+\Delta {L}_{f2}=2.704{\lambda}_{0}$ (2244 nm for${\lambda}_{0}\text{=}830nm$), which matches the value from the integrated sensor in the PZT. We then repeat the experiments under the same laser operation conditions for different displacements. The results are shown in Table 2, from which, it can be found that the results from these two methods are consistent and the average difference is about 60nm. The difference between them may come from the accuracy of the PZT integrated sensor, laser wavelength shift induced by external optical feedback, and the accuracy of the measured RO frequencies.

The resolution of this method is determined by the measured RO frequency and its variation range within ${\lambda}_{0}/2$. In this experiment, the frequency resolution of the oscilloscope in the experiments is 62.5 KHz and the variation of ${f}_{RO}$ in each period is measured as 66 MHz. Hence, in theory, if a perfect linearity can be guaranteed, the proposed experimental system can achieve a resolution of ${\lambda}_{0}/2112$ for the displacement measurement, which is 0.39 nm for λ_{0} = 830 nm. Therefore, the proposed method can be used for displacement measurement with high resolution.

## 4. Conclusion

Precision displacement measurement and positioning have become increasingly important in many areas, e.g. precision manufacturing, semiconductor industry, biotechnology, etc. With the progress of these fields, there is great demand for development of displacement sensing technique with high sensitivity, resolution and relatively simple configuration. In this work, a sensing system with an OFLD operating in POO is implemented for measuring displacement of a target. An analytical expression is presented for describing the sensing signals generated by the OFLD-POO system. It can be seen that both the amplitude and frequency of an OFLD-POO sensing signal are modulated by the displacement to be measured. Hence, we developed a new algorithm for measuring the displacement by simultaneously using the time-domain OFLD-POO signal waveform and RO frequency, which enables us to achieve displacement sensing with large measurement range, high sensitivity and resolution.

## References

**1. **S. Donati, G. Giuliani, and S. Merlo, “Laser diode feedback interferometer for measurement of displacements without ambiguity,” IEEE J. Quantum Electron. **31**(1), 113–119 (1995). [CrossRef]

**2. **T. Taimre, M. Nikolić, K. Bertling, Y. L. Lim, T. Bosch, and A. D. Rakić, “Laser feedback interferometry: a tutorial on the self-mixing effect for coherent sensing,” Adv. Opt. Photonics **7**(3), 570–631 (2015). [CrossRef]

**3. **D. Guo, M. Wang, and S. Tan, “Self-mixing interferometer based on sinusoidal phase modulating technique,” Opt. Express **13**(5), 1537–1543 (2005). [CrossRef] [PubMed]

**4. **M. Wang and G. Lai, “Self-mixing microscopic interferometer for the measurement of microprofile,” Opt. Commun. **238**(4–6), 237–244 (2004). [CrossRef]

**5. **S. Zhang, S. Zhang, L. Sun, and Y. Tan, “Fiber self-mixing interferometer with orthogonally polarized light compensation,” Opt. Express **24**(23), 26558–26564 (2016). [CrossRef] [PubMed]

**6. **M. Norgia, A. Pesatori, and L. Rovati, “Self-mixing laser Doppler spectra of extracorporeal blood flow: a theoretical and experimental study,” IEEE J. Sensors **12**(3), 552–557 (2012). [CrossRef]

**7. **U. Zabit, F. Bony, T. Bosch, and A. D. Rakic, “A self-mixing displacement sensor with fringe-loss compensation for harmonic vibrations,” IEEE Photonics Technol. Lett. **22**(6), 410–412 (2010). [CrossRef]

**8. **U. Zabit, O. D. Bernal, and T. Bosch, “Self-mixing laser sensor for large displacements: Signal recovery in the presence of speckle,” IEEE J. Sensors **13**(2), 824–831 (2013). [CrossRef]

**9. **L. Scalise and N. Paone, “Laser Doppler vibrometry based on self-mixing effect,” Opt. Lasers Eng. **38**(3–4), 173–184 (2002). [CrossRef]

**10. **M. Nikolić, E. Hicks, Y. L. Lim, K. Bertling, and A. D. Rakić, “Self-mixing laser Doppler flow sensor: an optofluidic implementation,” Appl. Opt. **52**(33), 8128–8133 (2013). [CrossRef] [PubMed]

**11. **Y. Yu, G. Giuliani, and S. Donati, “Measurement of the linewidth enhancement factor of semiconductor lasers based on the optical feedback self-mixing effect,” IEEE Photonics Technol. Lett. **16**(4), 990–992 (2004). [CrossRef]

**12. **Y. Yu, J. Xi, and J. F. Chicharo, “Measuring the feedback parameter of a semiconductor laser with external optical feedback,” Opt. Express **19**(10), 9582–9593 (2011). [CrossRef] [PubMed]

**13. **K. Lin, Y. Yu, J. Xi, H. Li, Q. Guo, J. Tong, and L. Su, “A Fiber-Coupled Self-Mixing Laser Diode for the Measurement of Young’s Modulus,” Sensors (Basel) **16**(6), 928 (2016). [CrossRef] [PubMed]

**14. **M. Dabbicco, A. Intermite, and G. Scamarcio, “Laser-self-mixing fiber sensor for integral strain measurement,” J. Lightwave Technol. **29**(3), 335–340 (2011). [CrossRef]

**15. **B. Liu, Y. Ruan, Y. Yu, J. Xi, Q. Guo, J. Tong, and G. Rajan, “Laser self-mixing fiber Bragg grating sensor for acoustic emission measurement,” Sensors (Basel) **18**(6), 1956 (2018). [CrossRef] [PubMed]

**16. **K. Bertling, J. Perchoux, T. Taimre, R. Malkin, D. Robert, A. D. Rakić, and T. Bosch, “Imaging of acoustic fields using optical feedback interferometry,” Opt. Express **22**(24), 30346–30356 (2014). [CrossRef] [PubMed]

**17. **R. Lang and K. Kobayashi, “External optical feedback effects on semiconductor injection laser properties,” IEEE J. Quantum Electron. **16**(3), 347–355 (1980). [CrossRef]

**18. **A. Ritter and H. Haug, “Theory of laser diodes with weak optical feedback. II. Limit-cycle behavior, quasi-periodicity, frequency locking, and route to chaos,” J. Opt. Soc. Am. B **10**(1), 145–154 (1993). [CrossRef]

**19. **B. Tromborg, J. H. Osmundsen, and H. Olesen, “Stability analysis for a semiconductor laser in an external cavity,” IEEE J. Quantum Electron. **20**(9), 1023–1032 (1984). [CrossRef]

**20. **A. Uchida, *Optical communication with chaotic lasers: applications of nonlinear dynamics and* synchronization (John Wiley & Sons, 2012).

**21. **Y. Fan, Y. Yu, J. Xi, and Q. Guo, “Dynamic stability analysis for a self-mixing interferometry system,” Opt. Express **22**(23), 29260–29269 (2014). [CrossRef] [PubMed]

**22. **B. Liu, Y. Yu, J. Xi, Y. Fan, Q. Guo, J. Tong, and R. A. Lewis, “Features of a self-mixing laser diode operating near relaxation oscillation,” Sensors (Basel) **16**(9), 1546 (2016). [CrossRef] [PubMed]

**23. **R. Teysseyre, F. Bony, J. Perchoux, and T. Bosch, “Laser dynamics in sawtooth-like self-mixing signals,” Opt. Lett. **37**(18), 3771–3773 (2012). [CrossRef] [PubMed]

**24. **B. Liu, Y. Yu, J. Xi, Q. Guo, J. Tong, and R. A. Lewis, “Displacement sensing using the relaxation oscillation frequency of a laser diode with optical feedback,” Appl. Opt. **56**(24), 6962–6966 (2017). [CrossRef] [PubMed]

**25. **D. Lenstra, “Relaxation oscillation dynamics in semiconductor diode lasers with optical feedback,” IEEE Photonics Technol. Lett. **25**(6), 591–593 (2013). [CrossRef]

**26. **Y. Ruan, B. Liu, Y. Yu, J. Xi, Q. Guo, and J. Tong, “Measuring Linewidth Enhancement Factor by Relaxation Oscillation Frequency in a Laser with Optical Feedback,” Sensors (Basel) **18**(11), 4004 (2018). [CrossRef] [PubMed]

**27. **J. Ohtsubo, Semiconductor *lasers: stability, instability and chaos* (Springer, 2012), Vol. 111.