Modeling for optical feedback laser diode operating in period-one oscillation and its application

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. Disciplines Engineering | Science and Technology Studies Publication Details B. Liu, Y. Ruan, Y. Yu, Q. Guo, J. Xi & J. Tong, "Modeling for optical feedback laser diode operating in period-one oscillation and its application," Optics Express, vol. 27, (4) pp. 4090-4104, 2019. Authors Bin Liu, Yuxi Ruan, Yanguang Yu, Qinghua Guo, Jiangtao Xi, and Jun Tong This journal article is available at Research Online: https://ro.uow.edu.au/eispapers1/2442 Modeling for optical feedback laser diode operating in period-one oscillation and its application BIN LIU, YUXI RUAN, YANGUANG YU,* QINGHUA GUO, JIANGTAO XI, AND JUN TONG School of Electrical, Computer and Telecommunications Engineering, University of Wollongong, Wollongong 2522, NSW, Australia yanguang@uow.edu.au 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. 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


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][4][5][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][19][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: 1+ / in C κτ α τ , where κ is the optical feedback strength, τ and in τ is the laser roundtrip time in the external and internal cavity respectively, and α 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.

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, ( ) t φ is the phase of the electric field, and ( ) N t is the carrier density. 2 ( ) E t is usually considered as the laser intensity or power, denoted by 2 ( ) ( ) where  is the modal gain per unit time.In the following derivation, the non-linear gain is ignored [18].There are 3 controllable parameters for a specific OFLD, i.e. injection current ( I ), feedback strength ( κ ) and external cavity laser round-trip time (τ ).The physical meanings and typical values of the other symbols appearing in Eqs. ( 1)-( 3) are listed in Table 1, which are adopted from [19,24].Starting from the L-K equations, we aim to get an approximate analytical expression for OFLD-POO signal (i.e., 2 ( ) E 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 ( ) t φ can be written in the form of below: ). E t and ( ) N t respectively.We will need to express the mean values and the modulation amplitudes as well as the initial phases using the parameters related to the LD and its external cavity listed in Table 1.
Regarding E and N , they can be determined by solving Eqs. ( 1 As Eq. ( 10) is an identical equation, letting the corresponding terms in both sides are equal, we have: Interestingly, it can be found that Eq. ( 13) is the same as Eq. ( 7).In a similar way, substituting ( ) E t , ( ) N t and ( ) t φ in Eq. ( 1) by Eqs. ( 4), ( 5) and ( 6), we can get: After expansion and reorganization of Eq. ( 14), we have: where 0 ( ) Δ and 1 ( ) Δ are the 0 th -order and 1 st -order Bessel function respectively.Note that higher harmonic components of r ω are neglected [18].Then, we have: After expansion and reorganization of Eq. ( 17), we have: As Eq. ( 18) is an identical equation as well, letting the corresponding terms in both sides are equal, we have: Inserting Eq. ( 12) into Eq.( 19) yields: Similarly, substituting ( ) E t , ( ) N t and ( ) t φ in Eq. ( 2) with Eqs. ( 4), ( 5) and ( 6), we get: sin .
Combining Eqs. ( 16) and ( 8), we have: With the Taylor expansion of the Bessel functions, we have the approximations: 20), ( 22), ( 23) and ( 24), we get: Inserting Eq. ( 27) into Eq.( 11) leads to Eq. ( 28) can be simplified as: where 2  E is shown in Eq. ( 7).Inserting Eq. ( 8) into Eq.( 7), in the case of 0.01 κ < , it can be rewritten as [1]: where 2 0 S E is the laser intensity of the LD without optical feedback, which is determined by the injection current.Inserting Eqs. ( 7) and ( 8) into Eq.( 25), we have , we obtain: 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 , where c is the vacuum light speed) [2,[11][12][13] as well as r ω [24].We rewrite Eq. (34) as below to include displacement sensing: where, 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 ( 0 / 2 λ ) [2].For a varying external target, the external cavity length can be expressed as where 0 L is the initial external cavity length, and 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 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. r ω contained in Eqs. ( 35) and (37), as in our previous work [24], which shows that r ω changes with displacement L Δ in the form of a saw-tooth waveform with a period of 0 / 2 λ when the OFLD system is in periodone oscillation.Thus, we have ) can be expressed as: where N is an integer, 0 ( )   .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 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 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 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: , 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 20 K > .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.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].

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: = .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 0 / 2 λ .As described in Fig. 1, there is a certain linear relationship L Δ , 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 .Thus, , 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 ) 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 RO f and displacement L Δ .  .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.634 GHz and 2.595 GHz respectively.From Fig. 7, it can be seen the variation range of RO f in each period is from 2.579 GHz to 2.645 GHz .Thus, we can get  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 0 / 2 λ .In this experiment, the frequency resolution of the oscilloscope in the experiments is 62.5 KHz and the variation of RO f 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 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.

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.
photon life time

ω
is the RO frequency of the LD without optical feedback, which is determined by the injection current.Inserting Eqs.(32) and (33), into Eq.(31) and considering the case of 0.01 κ <

rωr
is the angular RO frequency of the solitary laser diode which is determined by injection current density, - r offset ω is the offset value which is determined by the initial external cavity length and feedback level, and 0 operation conditions of an LD are fixed.Equation (37) deserves a close examination.Since both r periodical with same period of 0 / 2 λ for a varying external cavity length, we have:

λ.
Figure 1(a) shows an example of the relationship between the angular RO frequency and displacement of external target when 0 =24 threshold injection current.Figure 1(b) shows the waveform of external cavity with the same operation conditions.

Fig. 1 .
Fig. 1.(a) Relationship between RO frequency and displacement of external cavity, (b) waveform of

Figure 2
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 0 3λ , the initial external cavity length 0 =24 L cm , and the injection current 1.3 th I I = .Note that 23 r ω τ π ≈ 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 2.8 C = and 4.5 C = ; 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 r

Figure 2 (
a) is the displacement; Figs.2(b) and 2(c) show the numerically simulation results of L-K equations with 22(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 r ω and L Δ by

Fig. 2 .
Fig. 2. Verification of OFLD-POO signal expressed in Eq. (40) by the numerical results from L-K equations.(a) displacement of the external target; (b) and (c) results from L-K equations with C = 2.8 and C = 4.5 respectively; (d) and (e) results from Eq. (40) when C = 2.8 and C = 4.5 respectively.
by fringe counting with 0 / 2 λ resolution by using the envelop in OFLD-POO signal, and f L Δ is measured by the RO frequency for the fractional displacement within 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.

.r ω and 2 rωω ω and 2 rω
between r ω and the displacement within each fringe.As in Eq. (38), the maximum and minimum RO frequency within 0 We denote the RO angular frequency at the start and end point in a fractional fringe by 1 ω respectively.Then the displacement corresponding to fractional fringes can be obtained as: is the RO angular frequency at the start point of the displacement, for the last fractional fringe, is the RO angular frequency at the end point of the displacement.Because max r ω and min r ω are fixed for a LD with certain operation conditions, the displacement corresponding to fractional fringes can be determined by only measuring the RO frequency of the OFLD when displacement starts and ends.Note that the first and last fringes are always considered as fractional fringes.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 0 =3 L λ Δ 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 is the pulse train number).It can be seen from Fig.4(d) that there are 6 need to obtain the fractional displacement determined by the first fractional fringe and last fractional fringe denoted by
50 I mA = ,and 0 18 L cm = in this work.Before staring the measurement, we firstly need to get the relationship between RO frequency ( r ω ) and external target displacement ( L Δ ) when the external cavity length changing by 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 (
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