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

Suppression of undamped relaxation oscillation in a laser self-mixing interferometry sensing system

Open Access Open Access

Abstract

Undamped relaxation oscillation (RO) in a laser self-mixing interferometry (SMI) system may occur in some common application conditions, which may impact the stable operation of the system and degrade its sensing performance. In this work, we proposed to suppress the undamped RO by controlling the system operation parameters in a laser SMI sensing system. By numerically solving the famous Lang Kobayashi equations, the stability of a laser SMI system in a 3-parameter space of external cavity length, injection current and optical feedback factor were investigated. Based on the stability analyses, we determined the system operation conditions required for suppressing the undamped RO and derived an analytical expression for describing the relationship between the operation parameters. An experimental SMI system based on a laser diode (Sanyo, DL4140-001s) was implemented and verified the suppressing method. The experimental results showed that the SMI system in a moderate feedback regime can operate in steady state without undamped RO by setting proper operation parameters. This work provides useful guidance to design a stable SMI sensing system for practical applications.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Relaxation oscillation (RO) is one of the important characteristics of semiconductor lasers (SLs), which often occurs when the laser is turned on [1]. The basic physical mechanism of RO is the interaction between the in-cavity laser intensity and the population inversion, i.e., the cavity damping rate of the laser is higher than the carrier damping rate [2]. For an SL without external perturbations, RO is often damped after excitation. It can also be stimulated by perturbations, like optical injection or external optical feedback. With these perturbations, RO may evolve to be undamped [1]. In this case, SLs would experience abundant complex dynamics, e.g., period-one oscillation, quasi-periodical oscillation and chaos, which can be developed into different applications, e.g., chaotic communication [1], random bit generator [3], microwave photonics [4] and sensing of physical quantities [5,6]. On the other hand, undamped RO has great impact on the stable operation of SLs, which may degrade the system performance for those that require stable operation of SLs [7]. One example is the case of self-mixing interferometry (SMI), which is a promising high-precision measurement and sensing technology. The intrinsic mechanism of this technology is the phenomenon known as self-mixing effect, i.e., the laser output power and laser frequency are modulated when a proportion of the emitted laser light reflected by an external object re-injects the laser in-cavity, and mixes with the internal laser beam [8]. A typical SMI sensing system is composed of an SL with a build-in monitor photodiode (PD), a lens and an external target, which represents a minimum-part-count configuration as shown in Fig. 1. As compared with traditional two-beam interferometers, e.g., Michelson interferometer, SMI system has many unique characteristics and merits, such as the structure is simple and compact, the cost of implementation is low, and sensing signals are able to be detected anywhere the laser beam can reach [8]. To date, different SMI-based applications have been proposed, e.g., measurement of displacement [9,10], distance [1113], speed [1416], mechanical parameters of materials [17,18], linewidth enhancement factor [19,20], angle [21,22], acoustic emission [23], biomedical signals [2426], and monitoring plasma accumulation in laser ablation [27,28].

 figure: Fig. 1.

Fig. 1. Typical structure of a laser SMI sensing system

Download Full Size | PDF

In SMI systems, the SLs should operate in steady state with damped RO. The stable SMI model is the theory basis to analyze the system behavior, which is a set of stationary solutions to the famous Lang Kobayashi (L-K) equations [29], or is derived from the three-mirror model with the assumption that the SL is in steady state [8]. In the SMI model, there is an important parameter called feedback level factor (C). It determines the shape of SMI signals, and its values discriminate between different feedback regimes, i.e., weak feedback regime ($C\textrm{ } \le \textrm{ }1$), moderate/strong feedback regime (C > 1). In weak regime, a typical SMI signal is closely a sinusoidal curve. Whereas it becomes to be sawtooth-like in moderate/strong feedback regimes, and interference fringes may lose due to the existence of multi-solutions to the SMI phase equation and the induced laser mode competition and selection [8,30]. SMI-based applications are usually researched in weak or moderate feedback regime [8]. Strong feedback regime is often avoided in sensing applications because the system may be unstable or even enter chaotic state [8]. As compared with weak feedback regime, moderate feedback regime has unique characteristics and advantages in applications, e.g., SMI signals in moderate feedback regime show hysteresis and displacement can be retrieved by a one-channel interferometric signal without ambiguity [31]. Thus, many applications, e.g., measurement of displacement [31,32], velocity [33], and linewidth enhancement factor [19], are in moderate feedback regime.

However, some works have indicated that an SL with external optical feedback may suffer instability due to the existence of undamped RO with C > 1 [7,34,35]. More recently, we have found the a laser SMI system in moderate feedback regime may exhibit undamped RO and show complex SMI signals instead of conventional SMI signals by numerically solving the Lang Kobayashi equations [36]. Then we analyzed the SMI signal waveform features in details when undamped RO occurs by experiments, which shows SMI signals with undamped RO still contains some features that can be used for sensing potentially [37]. After that, we modeled the SMI signals with period-one oscillation, which is a special case when undamped RO happens [38]. We have found that SMI signals with undamped RO are in complicated forms, which requires extra signal processing to retrieve the required information [3638]. The widely-accepted stable SMI model would fail to completely describe the behaviors of a laser SMI system [38], and the sensing performance may be degraded [36]. Therefore, it is of great necessity and interest to suppress the undamped RO in a laser SMI sensing system. In [36], we have preliminarily discussed how to achieve a stable SMI signal by considering the external cavity length or the injection current separately but their relationship was not discussed for a certain feedback regime. In fact, we need to consider the feedback level factor, the external cavity length and the injection current simultaneously when building a laser SMI system.

In this work, we study to suppress the undamped RO by controlling the system operation parameters to design a stable SMI sensing system. Firstly, we investigated the stability of a laser SMI system in a 3-parameter space of external cavity length, injection current and optical feedback factor and discussed the combined influence of operation parameters on its stability and SMI signals by numerically solving the L-K equations. Based on the stability analyses, we determined the required operation parameters for suppressing the undamped RO. Then, an analytical expression for the required operation parameters was derived. Finally, an experimental system was built to verify the simulation analyses. The results are useful for designing a stable SMI sensing system with suppression of undamped RO.

2. Theory and simulation

2.1 Fundamentals of SMI

The behavior of a laser SMI system is usually described by the following mathematic model, which is the basis of SMI-based applications [8].

$${\phi _F}(\tau ) = {\phi _0}(\tau ) - \textrm{Csin}[{\phi _F}(\tau ) + \arctan \alpha ],$$
$$P(\tau ) = {P_0}[1 + m \times \cos ({\phi _F}(\tau )],$$
Here, Eq. (1) is called the phase equation. ${\phi _F}(\tau ) = {\omega _F}(\tau ) \cdot \tau$ and ${\phi _0} = {\omega _0}\tau$, where ${\omega _0}$ is the laser angular frequency without perturbation, ${\omega _F}(\tau )$ is that with optical feedback, $\tau = 2L/c$ is the roundtrip time of the laser light within the external cavity, where L is the distance between the laser emitting facet and the target surface, which is also called the external cavity length, c is light speed in the external cavity. $P(\tau )$ and ${P_0}$ are the laser emitting power or intensity with and without optical feedback respectively. m is the modulation index. The feedback level factor C in Eq. (1) is defined as $C = \eta \sqrt {1 + {\alpha ^2}} \kappa L\textrm{/}{L_{in}}{n_{in}}$, where $\kappa $ is the feedback strength which is governed by the reflectivity of the laser emitting facet and target surface, Lin is the length of laser internal cavity, nin is the refractive index of the laser internal cavity, $\alpha $ is the linewidth enhancement factor, $\eta $ is related to the possible loss on re-injection [8,39].

The SMI model can be derived from the L-K equations [29], shown as Eqs. (3)– (5), with the assumptions of the system operating in steady state with damped RO, i.e., by setting $dE(t)/dt = 0$, $d\phi (t)/dt = {\omega _F}(\tau ) - {\omega _0}$, and $dN(t)/dt = 0$[36]. Here $E(t)$, $\phi (t)$ is the amplitude and phase of the electric field respectively, and $N(t)$ is the carrier density. ${E^2}(t)$ is considered as the laser intensity or output power in the simulations as commonly treated in the literatures [35,36]. I is the injection current to the SL. The physical meanings of other parameters in L-K equations and their typical values for simulations in this work are in Table 1 [36,40].

$$\begin{aligned} \frac{{dE(t)}}{{dt}} &= \frac{1}{2}\{{{G_N}[N(t) - {N_0}][1 - \varepsilon \Gamma {E^2}(t)] - 1/{\tau_p}} \}E(t)\\ &+ \frac{\kappa }{{{\tau _{in}}}} \cdot E(t - \tau ) \cdot \cos [{\omega _0}\tau + \phi (t) - \phi (t - \tau )]\textrm{ }, \end{aligned}$$
$$\begin{aligned}\frac{{d\phi (t)}}{{dt}} &= \frac{1}{2}\alpha \{{{G_N}[N(t) - {N_0}][1 - \varepsilon \Gamma {E^2}(t)] - 1/{\tau_p}} \}\\ &- \frac{\kappa }{{{\tau _{in}}}} \cdot \frac{{E(t - \tau )}}{{E(t)}} \cdot \sin [{{\omega_0}\tau + \phi (t) - \phi (t - \tau )} ], \end{aligned}$$
$$\frac{{dN(t)}}{{dt}} = \frac{I}{{eV}} - \frac{{N(t)}}{{{\tau _s}}} - {G_N}[{N(t) - {N_0}} ]{E^2}(t).$$

Tables Icon

Table 1. Typical values of parameters for simulations

2.2 Stability analysis on a laser SMI system

We performed stability analyses on a laser SMI sensing system based on the L-K equations through Hopf bifurcation diagram, which is a useful tool to investigate the evolution of the system dynamics by sampling the local maximum and minimum points of the temporal waveforms. The bifurcation point in the diagram is the critical point that discriminates the steady state and period-one oscillation [1]. The Hopf bifurcation diagram was obtained by numerically solving the L-K equations. Since we focus on the SMI applications, the bifurcation diagram was presented in the plane of laser intensity and optical feedback factor C. Figure 2 shows an example of the Hopf bifurcation diagram of a laser SMI system when $I = 1.3{I_{th}}$, $L = 30.0\textrm{ cm}$, where Ith is the threshold injection current of an SL with ${I_{th}}\textrm{ = eV(}{N_0} + 1/{G_N}{\tau _p})/{\tau _s}$ [1] and it can be calculated as ${I_{th}} = \textrm{14 mA}$ with the values in Table 1. From Fig. 2, we can find a laser SMI system may operate in different states, i.e., steady state, period-one oscillation, quasi-periodical oscillation or chaos. For SMI-based applications, the system should be in steady state. The bifurcation point in Fig. 2 decides the maximum C that ensures the system is stable with damped RO. The RO evolves to be undamped when the feedback level factor C is larger than the value at the bifurcation point. In this case, the SMI system no longer operates at steady state and the SMI model of Eqs. (1) and (2) is inaccurate to describe the complete behavior of SMI system. Complex dynamics would appear in the SMI signals and influence the sensing performance of the SMI system.

 figure: Fig. 2.

Fig. 2. Bifurcation diagram of a laser SMI system with I = 1.3Ith, L = 30.0 cm

Download Full Size | PDF

A well-accepted method for describing the feedback regimes of a laser SMI system is based on the value of C, i.e., weak feedback regime ($C \le 1$), moderate feedback regime (1 < C < 4.6), and strong feedback regime (C > 4.6) [39,4144]. Because strong feedback regime is often avoided in sensing applications [8], we focus on the case of a laser SMI system in weak and moderate feedback regime. Additionally, because a laser SMI system in weak feedback regime is unconditionally stable [7,35], we mainly investigate a laser SMI system in moderate feedback regime. Although the work [30] has made some doubts on defining strong feedback regime with C > 4.6 from the aspect of SMI fringe loss, we still adapted the common definition of moderate feedback regime in the community, i.e., with 1 < C < 4.6, because the proposed analyses can be easily to extend to C > 4.6. We denote the C at the bifurcation point by Cbp. We can obtain that Cbp is about 2.85 in Fig. 2, where the dash line is for C = 4.6. It can be predicted that the SMI system in the case as Fig. 2 is not always stable in moderate feedback regime, i.e., it is unstable with 2.85 < C <4.6. In order to confirm it, we then numerically solved the L-K equations to obtain the laser intensity ${E^2}(t)$, i.e., the SMI signals, when a sinusoidal displacement was applied on the external target. Figure 3 presents the simulation results, where Fig. 3(a) is the displacement, Figs. 3(b) - (d) are the corresponding SMI signals when C = 2.80, C = 2.90 and C = 4.00 respectively. We can see that Fig. 3(b) is the typical SMI signal which can be described by the stable SMI model of Eqs. (1) and (2), while Figs. 3(c) and (d) are in complicated patterns containing complex dynamics. As a result, Cbp should be larger than 4.6, i.e., Cbp > 4.6, to keep a moderate feedback SMI system stable. Therefore, the stability of a laser SMI system is described by Cbp. The bifurcation point in the bifurcation diagram as in Fig. 2 should be on the right side of the dashed line to maintain the system stable.

 figure: Fig. 3.

Fig. 3. (a) Displacement of the target; (b), (c) and (d) the corresponding SMI signals with C = 2.80, C = 2.90 and C = 4.00, respectively

Download Full Size | PDF

We then investigated the influence of external cavity length L and injection current I on Cbp. In this work, we mainly studied a laser SMI system with long external cavity length, i.e., L > c/2fRO, where fRO is the RO frequency of the solitary SL [34,35]. We can calculate with fRO > 2 GHz with the values in Table 1 [1], and thus a long external cavity length corresponds to L >7.5 cm. As a result, we varied the external cavity length with a step of 0.5 cm in the range from 10.0 cm to 100.0 cm, the injection current I from 1.1Ith to 2.0Ith with a step of 0.01Ith, then got $180 \times 90 = 16200$ pairs of L and I. For each pair of L and I, we numerically solved the L-K equations to get the bifurcation point and recorded the Cbp. Figure 4 presents the simulation results, where Fig. 4(a) is the Cbp described by both L and I. It can be seen that Cbp increases with fluctuations when L increases for a fixed injection current. The similar variation trend can be found for Cbp with respect to I for a fixed external cavity length. Figures 4(b) and (c) are two examples for a fixed injection current of I = 1.5Ith and a fixed external cavity length of L = 30.0 cm respectively. It can be concluded that a longer external cavity and a larger injection contribute to a larger Cbp.

 figure: Fig. 4.

Fig. 4. The feedback level factor at bifurcation point, (a) described by external cavity length L and injection current I, (b) described by L with a fixed injection current I = 1.5Ith, (c) described by injection current with a fixed external cavity length L = 30.0 cm

Download Full Size | PDF

2.3 Operation requirements for suppression of undamped RO

Based on the above simulations, we set Cbp= 4.6 and are able to get the required critical injection current and external cavity length that ensure a moderate feedback SMI sensing system is stable. The simulation results were obtained as shown in Fig. 5. It can be seen that the relationship between the required critical injection current and external cavity length is a hyperbolic curve approximately. The SMI system remains stable when it operates with the parameters in the area above the curve. Note that, for the numerical results in Fig. 5, there is an average uncertainty of about 5% when we use the bifurcation diagram to discriminate the operation states of the SMI system due to the limit of simulation step of time. In fact, if we intend to apply SMI technique in strong feedback regime [45] and maintain the system stable with damped RO, we just need to set Cbp to be equal to the required feedback level factor C. Then the curve for the required external cavity length and injection current in strong feedback regime can be easily obtained with the same procedure.

 figure: Fig. 5.

Fig. 5. The critical injection current and external cavity length for damped RO in moderate feedback regime.

Download Full Size | PDF

We then try to find an analytical expression for the relationship between the required injection current and external cavity length. Our recent work has found an approximate expression for the critical feedback level factor for a stable SMI system as below [36]:

$${C_{critical}} = {\left( {\frac{\Omega }{{{\omega_R}}}} \right)^2}\frac{{L\sqrt {1 + {\alpha ^2}} }}{{{\tau _R}c\,{{\sin }^2}({{{\Omega L} / c}} )}}\frac{{\sqrt {{{[{2{{({{\Omega / {{\omega_R}}}} )}^2} - 1} ]}^2} + {\alpha ^2}} }}{{{{[{2{{({{\Omega / {{\omega_R}}}} )}^2} - 1} ]}^2} + {\alpha ^2}}}.$$
Here ${\omega _R}$ and $\Omega $ is the RO angular frequency of an SL without and with external optical feedback respectively, ${\tau _R}$ is the damping time of RO for an SL without perturbations. They can be expressed as [1,40]:
$${\omega _R} = \sqrt {\frac{{(I/{I_{th}} - 1) (1 + {G_N}{N_0}{\tau _p})}}{{{\tau _s}{\tau _p}}}}, $$
$${\Omega ^2} - \omega _R^2 = \frac{\Omega }{{{\tau _R}}}\cot (\frac{{\Omega \tau }}{2}), $$
$$\frac{1}{{{\tau _R}}} = \frac{1}{{{\tau _s}}} + \left( {{\tau_p} + \frac{{\varepsilon \Gamma }}{{{G_N}}}} \right)\omega _R^2.$$

The critical feedback level factor has been proposed to distinguish the stable and unstable regions of a laser SMI system in the plane of optical feedback factor and optical phase which has the similar meaning with Cbp, but it has been derived under the assumption of ${\omega ^2}_R > > \kappa /({\tau _{in}}{\tau _R})$ and ${\omega ^2}_R > > {(\kappa /{\tau _{in}})^2}$[36]. We take ${C_{bp}} \approx {C_{critical}}$ for approximation. Then, for a long external cavity (i.e., $L > \pi {c / {{\omega _R}}}$), we have $\Omega \approx {\omega _R}$[40]. We insert Eqs. (7)– (9) into Eq. (6), and set Cbp =4.6 by considering the moderate feedback regime, and then get:

$$\frac{L}{c}\left[ {\left( {\frac{{(I/{I_{th}} - 1)({N_0}{G_N}{\tau_p} + 1)}}{{{G_N}{\tau_p}{\tau_s}}}} \right)({{G_N}{\tau_p} + \varepsilon \Gamma } )+ \frac{1}{{{\tau_s}}}} \right]\textrm{ = }4.6.$$

Equations (10) gives the relationship between the required injection current I and external cavity length L for a laser SMI system in moderate feedback regime, where other parameters are constant for a certain SMI system except I and L. The result is also plotted in Fig. 5 for comparison. It can be seen that the curve of Eq. (10) is consistent with the simulation results. The relationship between the required critical injection current and external cavity length is a hyperbolic curve. In order to guarantee a laser SMI system stable in the whole moderate feedback regime, the injection current and external cavity length should be set above the hyperbolic curve.

3. Experiment

Experiments were then carried out to verify the simulation results. Figure 6 is the structure diagram of the experimental setup. The SL is a 790 nm single mode laser diode (Sanyo, DL4140-001s). Its threshold value of injection current is 30 mA. The SL is driven by an SL controller (Thorlabs, ITC4001) with integrated temperature controller. The laser beam transmits through a collimating lens and a 50:50 beam splitter (BS). One proportion of the beam transmits via an optical attenuator (Thorlabs, NDC-50C-2M-B) and is then reflected by a mirror affixed on a piezoelectric transducer (PZT) (Thorlabs, PAS009). The PZT is driven by a PZT controller (Thorlabs, MDT694) and is used to generate a small displacement. The PZT with the mirror assembled on a linear translation stage acts as the external target. Using the linear stage, different external cavity lengths can be obtained. The optical attenuator is used to adjust the feedback strength. By using the attenuator, mirror and the linear stage, different feedback level factors can be achieved. The other part of the light through the BS transmits to a fast external photodetector (Thorlabs, PDA8GS) via an optical isolator (Thorlabs, IO-5-780) and it is detected as the SMI signals. Finally, we use an oscilloscope with an 8 GHz analog bandwidth (Tektronix DSA 70804) to capture and observe the SMI signals.

 figure: Fig. 6.

Fig. 6. Experimental setup

Download Full Size | PDF

It is noteworthy that the structure of the experimental setup is slightly different from the typical SMI configuration in Fig. 1. We used an external photodetector with a frequency bandwidth of 9.5 GHz rather than the photodiode packaged in the SL to detect the SMI signals. The reasons are twofold. One is that the RO frequency of an SL is usually in GHz scale. The response speed of the internal packaged photodiode and the common trans-impedance amplifier is not fast enough to detect the complete SMI signals with undamped RO, e.g,. the overall bandwidth of the integrated photodiode with the detection circuit in our experimental system is about 10 MHz. In fact, if we use the integrated photodiode to detect the SMI signals with undamped RO, the detected SMI signals could be like the normal ones but with distortions for different detection bandwidths. In order to keep the detected SMI signals consistent with the simulation results in Fig. 3 and recognize undamped RO, we need a fast photodetector with a wide frequency response bandwidth. The other is that SMI signals can be detected anywhere along the optical path, which means the structure in Fig. 6 can be still treated as a laser SMI configuration [8].

The SL operation temperature was stabilized at $25 \pm 0.01^\circ \textrm{C}$ during the experiments by using the SL controller. According to the simulations, Cbp increases with the injection current and external cavity length. We firstly fixed an external cavity length with a certain value and then carefully adjust the attenuator to set C =4.6. After that, we applied different injection currents to the SL from large to small and recorded the injection current when undamped RO started to appear in the SMI signals. Then we set another external cavity length, and recorded the corresponding injection current by the same method. Nine pairs of injection currents and external cavity lengths were recorded. Note that we are able to use the variable optical attenuator to adjust the optical feedback strength to set C =4.6 for a given external cavity length thanks to the ability of continuously adjusting the transmission rate via rotation of the optical attenuator. Although C is related to many paramters, like linewith enhancement factor $\alpha $, laser internal cavity length, ect., it is still possible to set C with the required value, e.g., C =4.6 in this work. That is because C is directly proportional to feedback strength $\kappa $, and $\kappa $ is proportional to the power transmission rate of the attenuator. Once we obtained a reference value for C, we just need to adjust the power transmission rate of the attenuator to get C = 4.6 even though other parameters related to C are unknown but with fixed values. The value of C is determined as follows:

  • 1. For each external cavity length, adjust the injection current and get a typical stable moderate feedback SMI signal by observing the waveform.
  • 2. Use the method in [19] to estimate the value of C denoted by Cref and record the transmission rate of the optical attenuator denoted by Tref, which can be read from the dial of the attenuator.
  • 3. Adjust the transmission rate (T) of the attenuator and obtain the required feedback level factor C with the expression of $T = C{T_{ref}}/{C_{ref}}$, e.g., during the experiments, we obtained a typical SMI signal with ${C_{ref}}\textrm{ = }1.5$ and recorded the corresponding transmission rate with ${T_{ref}}\textrm{ = 10\%}$. Thus we can get the required transmission rate with $T \approx 31\mathrm{\%}$ for C = 4.6.
Figure 7 shows the experimental results. It can be found that the relationship between the injection current I and external cavity length L is an approximate hyperbolic curve, which is consistent with the theoretical results. As shown in the case of Fig. 7, when we choose a pair of L and I on the curve as the operation parameters, the SMI system with C =4.6 starts to exhibit undamped RO. When L and I are set with the values below the curve, undamped RO and different complex dynamics, e.g., period-one oscillation, quasi-periodical oscillation, or chaos may appear. In order to guarantee the SMI system is stable, L and I should be set with values above the curve. Although the curve in Fig. 7 is for moderate feedback regime, we are able to get a curve for operation conditions for a stable SMI system even in strong feedback regime by following the above experimental procedure. Note that, when we intend to determine the required transmission rate for C = 4.6, we need to estimate a value of Cref in advance, whereas the method we used has a uncertainty of about 8%, which may introduce some errors to the obtained optical feedback factor.

 figure: Fig. 7.

Fig. 7. The critical injection current and the external cavity length when undamped RO start to appear in SMI signals

Download Full Size | PDF

We took two pairs of L and I, i.e., point A (L = 12.0 cm, I = 45 mA) and point B (L = 20.0 cm, I = 60 mA) in Fig. 7 for the examples. Figure 8 presents the observed SMI signals with the L and I at points A and B respectively. Here, Fig. 8(a) is the controlling voltage applied on the PZT. Note that each 0.1 V controlling voltage applied on the PZT generates a displacement of 53 nm. The peak-peak amplitude of the PZT controlling voltage in Fig. 8(a) is around 3.8 V, which produces a displacement of about 2024nm. Figures 8(b)-(e) are the corresponding SMI signals, where (b) and (c) are for point A; (d) and (e) are for point B. In addition, Figs. 8(b) and (d) are with C = 1.6; (c) and (e) are with C = 4.6. Based on Fig. 7, we can predict that the SMI system will be unconditional stable with any C≤4.6 for point B. But for point A, we cannot ensure the system is stable or not in the whole moderate feedback regime. From Fig. 8, It can be found that the SMI signal is with typical form with C = 1.6, i.e., Fig. 8(b), but contains complex dynamics with C = 4.6, i.e., Fig. 8(c), at point A. When we set L and I with the values at point B, the SMI signals do not contain any complex dynamics no matter with C = 1.6 or C = 4.6, demonstrating that the SMI system is always stable without undamped RO in the moderate feedback regime. It is noteworthy that Figs. 8(b) and (d) are similar, which is because the system with C = 1.6 is stable when it operates at both point A and point B.

 figure: Fig. 8.

Fig. 8. Experimental SMI signals with different I and L, (a) the controlling voltage of PZT; (b) and (c) the corresponding SMI signals at point A; (d) and (e) the corresponding SMI signals at point B.

Download Full Size | PDF

4. Conclusion

In this work, we studied to suppress the undamped RO of SLs by controlling the system operation parameters in a laser SMI sensing system, especially in moderate feedback regime. The stability of a laser SMI system in a 3-parameter space of external cavity length, injection current and optical feedback factor was numerically analyzed through Hopf bifurcation diagram. The system is unstable with undamped RO when the feedback level factor is larger than that at the bifurcation point, which leads to the SMI signal waveforms and measurement performance of a laser SMI system are impacted by complex laser dynamics. Based on the stability analyses, we determined the system operation condition required for suppressing the undamped RO and derived an analytical expression for describing the relationship between the operation parameters, which shows that the relationship between the required injection current and external cavity is an approximate hyperbolic curve. We also derived an analytical expression for the required injection current and external cavity length, which is consistent with the numerical simulations. The SMI system is guaranteed to be stable with damped RO when it operates with the parameters above the curve. The theoretical and simulation results were verified by experiments. The experimental results showed that the built SMI system in moderate feedback regime can operate in steady state without undamped RO by setting proper parameters. The results in this work provide useful guidance to design a stable SMI sensing system without effects of complex laser dynamics for practical applications.

Funding

National Natural Science Foundation of China (62005234); Education Department of Hunan Province (20C1791).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. A. Uchida, Optical communication with chaotic lasers: applications of nonlinear dynamics and synchronization (John Wiley & Sons, 2012, Chaps. 2&4).

2. E. Lacot, R. Day, and F. Stoeckel, “Coherent laser detection by frequency-shifted optical feedback,” Phys. Rev. A 64(4), 043815 (2001). [CrossRef]  

3. I. Kanter, Y. Aviad, I. Reidler, E. Cohen, and M. Rosenbluh, “An optical ultrafast random bit generator,” Nat. Photonics 4(1), 58–61 (2010). [CrossRef]  

4. C.-H. Tseng, Y.-H. Hung, and S.-K. Hwang, “Frequency-modulated continuous-wave microwave generation using stabilized period-one nonlinear dynamics of semiconductor lasers,” Opt. Lett. 44(13), 3334–3337 (2019). [CrossRef]  

5. S. D. Cohen, A. Aragoneses, D. Rontani, M. Torrent, C. Masoller, and D. J. Gauthier, “Multidimensional subwavelength position sensing using a semiconductor laser with optical feedback,” Opt. Lett. 38(21), 4331–4334 (2013). [CrossRef]  

6. 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]  

7. G. A. Acket, D. Lenstra, A. Den Boef, and B. Verbeek, “The influence of feedback intensity on longitudinal mode properties and optical noise in index-guided semiconductor lasers,” IEEE J. Quantum Electron. 20(10), 1163–1169 (1984). [CrossRef]  

8. 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]  

9. K. Zhu, B. Guo, Y. Lu, S. Zhang, and Y. Tan, “Single-spot two-dimensional displacement measurement based on self-mixing interferometry,” Optica 4(7), 729–735 (2017). [CrossRef]  

10. J. I. Khan and U. Zabit, “Sparsity Promoting Frequency Sampling Method for Estimation of Displacement and Self-Mixing Interferometry Parameters,” IEEE Sens. J. 21(9), 10587–10594 (2021). [CrossRef]  

11. M. Norgia, G. Giuliani, and S. Donati, “Absolute distance measurement with improved accuracy using laser diode self-mixing interferometry in a closed loop,” IEEE Trans. Instrum. Meas. 56(5), 1894–1900 (2007). [CrossRef]  

12. Y. Zhao, C. Wang, Y. Zhao, D. Zhu, and L. Lu, “An All-Fiber Self-Mixing Range Finder With Tunable Fiber Ring Cavity Laser Source,” J. Lightwave Technol. 39(12), 4217–4224 (2021). [CrossRef]  

13. D. Guo and M. Wang, “Self-mixing interferometry based on a double-modulation technique for absolute distance measurement,” Appl. Opt. 46(9), 1486–1491 (2007). [CrossRef]  

14. D. Guo, H. Jiang, L. Shi, and M. Wang, “Laser self-mixing grating interferometer for MEMS accelerometer testing,” IEEE Photon. J. 10(1), 1–9 (2018). [CrossRef]  

15. L. Scalise, Y. Yu, G. Giuliani, G. Plantier, and T. Bosch, “Self-mixing laser diode velocimetry: application to vibration and velocity measurement,” IEEE Trans. Instrum. Meas. 53(1), 223–232 (2004). [CrossRef]  

16. X. Zhang, W. Gu, C. Jiang, B. Gao, and P. Chen, “Velocity measurement based on multiple self-mixing interference,” Appl. Opt. 56(23), 6709–6713 (2017). [CrossRef]  

17. 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 16(6), 928 (2016). [CrossRef]  

18. B. Wang, B. Liu, L. An, P. Tang, H. Ji, and Y. Mao, “Laser Self-Mixing Sensor for Simultaneous Measurement of Young's Modulus and Internal Friction,” Photonics 8(12), 550 (2021). [CrossRef]  

19. 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]  

20. B. Liu, Y. Ruan, Y. Yu, B. Wang, and L. An, “Influence of feedback optical phase on the relaxation oscillation frequency of a semiconductor laser and its application,” Opt. Express 29(3), 3163–3172 (2021). [CrossRef]  

21. Y. Zhao, R. Xiang, Z. Huang, J. Chen, and L. Lu, “Research on the Multi-longitudinal mode laser self-mixing static angle-measurement system using a right-angle prism,” Measurement 162, 107906 (2020). [CrossRef]  

22. C. Wang, X. Fan, G. You, H. Gui, H. Wang, J. Liu, B. Yu, and L. Lu, “Full-circle range and microradian resolution angle measurement using the orthogonal mirror self-mixing interferometry,” Opt. Express 26(8), 10371 (2018). [CrossRef]  

23. B. Liu, Y. Ruan, and Y. Yu, “All-Fiber Laser-Self-Mixing Sensor for Acoustic Emission Measurement,” J. Lightwave Technol. 39(12), 4062–4068 (2021). [CrossRef]  

24. Y. Wei, W. Huang, Z. Wei, J. Zhang, T. An, X. Wang, and H. Xu, “Double-path acquisition of pulse wave transit time and heartbeat using self-mixing interferometry,” Opt. Commun. 393, 178–184 (2017). [CrossRef]  

25. A. Arasanz, F. J. Azcona, S. Royo, A. Jha, and J. Pladellorens, “A new method for the acquisition of arterial pulse wave using self-mixing interferometry,” Opt. Laser Technol. 63, 98–104 (2014). [CrossRef]  

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

27. S. Donadello, V. Finazzi, A. G Demir, and B. Previtali, “Time-resolved quantification of plasma accumulation induced by multi-pulse laser ablation using self-mixing interferometry,” J. Phys. D: Appl. Phys. 53(49), 495201 (2020). [CrossRef]  

28. S. Donadello, A. G. Demir, and B. Previtali, “Probing multipulse laser ablation by means of self-mixing interferometry,” Appl. Opt. 57(25), 7232–7241 (2018). [CrossRef]  

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

30. M. Veng, J. Perchoux, and F. Bony, “Fringe Disappearance in Self-Mixing Interferometry Laser Sensors: Model and Application to the Absolute Distance Measurement Scheme,” IEEE Sens. J. 19(14), 5521–5528 (2019). [CrossRef]  

31. C. Bes, G. Plantier, and T. Bosch, “Displacement measurements using a self-mixing laser diode under moderate feedback,” IEEE Trans. Instrum. Meas. 55(4), 1101–1105 (2006). [CrossRef]  

32. 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]  

33. A. Magnani, A. Pesatori, and M. Norgia, “Real-Time Self-Mixing Interferometer for Long Distances,” IEEE Trans. Instrum. Meas. 63(7), 1804–1809 (2014). [CrossRef]  

34. K. Petermann, “External optical feedback phenomena in semiconductor lasers,” IEEE J. Select. Topics Quantum Electron. 1(2), 480–489 (1995). [CrossRef]  

35. S. Donati and M. T. Fathi, “Transition from short-to-long cavity and from self-mixing to chaos in a delayed optical feedback laser,” IEEE J. Quantum Electron. 48(10), 1352–1359 (2012). [CrossRef]  

36. 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]  

37. 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 16(9), 1546 (2016). [CrossRef]  

38. B. Liu, Y. Ruan, Y. Yu, Q. Guo, J. Xi, and J. Tong, “Modeling for optical feedback laser diode operating in period-one oscillation and its application,” Opt. Express 27(4), 4090–4104 (2019). [CrossRef]  

39. Y. Yu, J. Xi, J. F. Chicharo, and T. M. Bosch, “Optical Feedback Self-Mixing Interferometry With a Large Feedback Factor C: Behavior Studies,” IEEE J. Quantum Electron. 45(7), 840–848 (2009). [CrossRef]  

40. J. Mork, B. Tromborg, and J. Mark, “Chaos in semiconductor lasers with optical feedback: theory and experiment,” IEEE J. Quantum Electron. 28(1), 93–108 (1992). [CrossRef]  

41. M. Norgia and S. Donati, “A displacement-measuring instrument utilizing self-mixing interferometry,” IEEE Trans. Instrum. Meas. 52(6), 1765–1770 (2003). [CrossRef]  

42. J. I. Khan and U. Zabit, “Deformation Method of Self-Mixing Laser Sensor’s Feedback Phase for Estimation of Optical Feedback Coupling Factor and Displacement,” IEEE Sens. J. 21(6), 7490–7497 (2021). [CrossRef]  

43. V. Annovazzi-Lodi, S. Merlo, and M. Norgia, “Measurements on a micromachined silicon gyroscope by feedback interferometry,” IEEE/ASME Trans. Mechatron. 6(1), 1–6 (2001). [CrossRef]  

44. C.-M. Ri, C.-H. Kim, Y.-N. Oh, and S.-C. Kim, “Immediate estimation of feedback factor and linewidth enhancement factor from measured self-mixing signals under moderate or strong regime,” Meas. Sci. Technol. 31(6), 065204 (2020). [CrossRef]  

45. J. El Assad, T. Bosch, and G. Plantier, “Laser diode under strong feedback for mechatronics applications,” in Sensors, 2007 IEEE, (IEEE, 2007), 387–390.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

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 (8)

Fig. 1.
Fig. 1. Typical structure of a laser SMI sensing system
Fig. 2.
Fig. 2. Bifurcation diagram of a laser SMI system with I = 1.3Ith, L = 30.0 cm
Fig. 3.
Fig. 3. (a) Displacement of the target; (b), (c) and (d) the corresponding SMI signals with C = 2.80, C = 2.90 and C = 4.00, respectively
Fig. 4.
Fig. 4. The feedback level factor at bifurcation point, (a) described by external cavity length L and injection current I, (b) described by L with a fixed injection current I = 1.5Ith, (c) described by injection current with a fixed external cavity length L = 30.0 cm
Fig. 5.
Fig. 5. The critical injection current and external cavity length for damped RO in moderate feedback regime.
Fig. 6.
Fig. 6. Experimental setup
Fig. 7.
Fig. 7. The critical injection current and the external cavity length when undamped RO start to appear in SMI signals
Fig. 8.
Fig. 8. Experimental SMI signals with different I and L, (a) the controlling voltage of PZT; (b) and (c) the corresponding SMI signals at point A; (d) and (e) the corresponding SMI signals at point B.

Tables (1)

Tables Icon

Table 1. Typical values of parameters for simulations

Equations (10)

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

ϕ F ( τ ) = ϕ 0 ( τ ) Csin [ ϕ F ( τ ) + arctan α ] ,
P ( τ ) = P 0 [ 1 + m × cos ( ϕ F ( τ ) ] ,
d E ( t ) d t = 1 2 { G N [ N ( t ) N 0 ] [ 1 ε Γ E 2 ( t ) ] 1 / τ p } E ( t ) + κ τ i n E ( t τ ) cos [ ω 0 τ + ϕ ( t ) ϕ ( t τ ) ]   ,
d ϕ ( t ) d t = 1 2 α { G N [ N ( t ) N 0 ] [ 1 ε Γ E 2 ( t ) ] 1 / τ p } κ τ i n E ( t τ ) E ( t ) sin [ ω 0 τ + ϕ ( t ) ϕ ( t τ ) ] ,
d N ( t ) d t = I e V N ( t ) τ s G N [ N ( t ) N 0 ] E 2 ( t ) .
C c r i t i c a l = ( Ω ω R ) 2 L 1 + α 2 τ R c sin 2 ( Ω L / c ) [ 2 ( Ω / ω R ) 2 1 ] 2 + α 2 [ 2 ( Ω / ω R ) 2 1 ] 2 + α 2 .
ω R = ( I / I t h 1 ) ( 1 + G N N 0 τ p ) τ s τ p ,
Ω 2 ω R 2 = Ω τ R cot ( Ω τ 2 ) ,
1 τ R = 1 τ s + ( τ p + ε Γ G N ) ω R 2 .
L c [ ( ( I / I t h 1 ) ( N 0 G N τ p + 1 ) G N τ p τ s ) ( G N τ p + ε Γ ) + 1 τ s ]  =  4.6.
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.