1. Introduction
Investigation of mode-locked laser noises is important both for studying delicate laser dynamics and for optimizing laser performance. Theoretically, by introducing proper noise sources into the Master equation model for mode-locked lasers, different laser noise properties, such as the relative intensity noises (RIN), as well as the timing jitter noises, can be calculated based on soliton perturbation theory [
1,
2,
3,
4,
5,
6] or stochastic numerical simulation [
7,
8,
9,
10,
11]. In the literature, both passive and active mode-locked lasers have been studied extensively, and the agreement between theoretical and experimental results are reasonable for the considered single-pulse mode-locking cases [
1,
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17]. It is well known that some mode-locked lasers can be operated under a bound-pulse mode-locking state, in which two or more adjacent bound-pulses are formed inside the laser cavity [
18,
19,
20,
21,
22,
23,
24,
25]. Because of the sophisticated balancing effects among different optical linear/nonlinear processes, bound-pulse mode-locking states can exhibit more interesting appearances. The noise properties of bound-pulse mode-locked lasers have also been less investigated, but remain of particular research interest. Some preliminary experimental results have been reported in our early work [
26]. Due to the more complicated pulse shape of the bound-pulses, analytic soliton perturbation theory is no longer applicable. The stochastic numerical simulation is applicable in principle, but in practice will be extremely time consuming, due to the required computation efforts. To overcome these difficulties, in the present work we develop a new theoretical method based on linearized backpropagation to calculate the laser RIN noise spectra of bound-soliton mode-locking. The linearized backpropagation approach was originally developed for modeling quantum nonlinear optical pulse propagation problems [
27,
28,
29]. Here we generalize the method to calculate the laser noise spectra, as well as the total noise variance, for the first time. Due to the deterministic computational nature of the approach, the required computation is much less than that of the stochastic simulation method. Since no particular pulse shape is assumed, the developed new approach is a general method applicable to the considered bound-pulse mode-locking case. As a real example, we investigate an environmentally stable 10 GHz hybrid mode-locked Er-doped fiber laser with active phase modulation. The laser is experimentally built with the interesting property that its mode-locking state can transmit from the single-soliton to the bound-soliton mode-locking state by simply increasing pump power [
25,
26]. The RIN noises of the laser are experimentally measured under both the bound- and single-soliton states for comparison. The observed experimental tendency agrees reasonably with the theoretical prediction based on the linearized backpropagation approach. One interesting finding is that the RIN noise of the bound-soliton state can be lower than that of the single-soliton state by following the laser power scaling tendency. The results also demonstrate that the developed theoretical method is very well-suited for studying general mode-locked laser noise problems.
2. Materials and Methods
Our experimental mode-locked fiber laser configuration is constructed as a sigma-shape cavity which is able to achieve environmentally-stable operation, as indicated in
Figure 1 [
30,
31,
32,
33,
34]. The cavity is composed of a polarization-maintaining (PM) reflection loop, and a non-PM linear fiber section with a Farady rotator mirror for automatic birefringence cancellation. In the PM-loop, a LiNbO
3 phase modulator is inserted to provide an active mode-locking mechanism, and the phase modulation frequency is set at 10 GHz so that the laser can be harmonically mode-locked under a high repetition rate. In the non-PM linear section, a piece of 4 m dispersion-shifted high nonlinearity fiber (DS-HNLF) is used to enhance the cavity Kerr nonlinearity, and a piece of bent dispersion-shifted fiber (DSF) is inserted to adjust the lasing wavelength through its wavelength-dependent loss. The Er-doped gain laser is bi-directionally pumped so as to achieve enough lasing power for observing the bound-soliton mode-locking. Notice that in our mode-locked fiber laser, there is no equivalent saturable absorption effect [
33,
34]. The bound-soliton state can be formed in the higher pump power levels in comparison to the single-soliton state [
25]. The RIN noises of both states are measured with a fast photo-detector and a radio frequency (RF) spectrum analyzer (Agilent E4440A, Santa Clara, CA, USA) for detailed comparison.
Our theoretical model is based on the Master equation, as shown in Equation (1), which assumes all the cavity parameters are averagely distributed. The physical meanings of the symbols are described as follow. : complex pulse envelope, : small time scale, : large time scale in terms of the cavity roundtrip number, : small signal gain, : gain saturation energy, : linear loss, : cavity dispersion, optical filtering, : equivalent nonlinear saturable absorption (=0 for the studied laser), : self-phase modulation (SPM), active phase modulation depth, : modulation frequency, and : noise source.
The backpropagation method starts from the linearization of the above nonlinear equation. We assume a small perturbation field around the mean field solution and substitute into Equation (1) to derive the linearized equation for . The result is shown in Equation (2).
One can rewrite Equation (2) in a compact form, as in Equation (3). Note that the operator S now denotes a special operator that operates both on and as defined in Equation (2).
Corresponding to the above linearized perturbation equation, we then define an adjoint equation such that the inner-product conservation between two systems can be ensured. By introducing the adjoint operator of
S according to the following Equation (4),
the adjoint equation is defined by Equation (5) as below:
Here, the inner product of two functions is defined as . The complete form of the adjoint equation is given in Equation (6), where represents the mean field pulse energy.
One can easily prove that the solutions of the linearized equation and its adjoint will satisfy:
Equation (7) is main result of the backpropagation method, and is very physically meaningful. The final noise parameter defined by the projection in the left-hand side of Equation (7) comes from two contributions: the initial noise (the first term in the right-hand side) and the integrated added noises (the second term in the right-hand side). Thus, one can obtain information about the final noise parameter after numerically solving the adjoint system backwardly (from T to 0) with the initial condition given by , which is the projection function for the noise parameter to be analyzed. By choosing a suitable projection function, the inner product in the left-hand side of Equation (7) can correspond to the perturbation of the chosen pulse parameter. From the definition of the pulse energy , for calculating the pulse energy noise of the optical pulse, the projection function is simply 2 times the mean field, as shown in Equation (8).
The noise variance
and the noise auto-correlation function
can also be derived after knowing the statistical properties (i.e., white noises) of the noise source term
. In the following, we will simply assume the noise is delta-function correlated, as given by Equation (9) [
1,
2,
3,
4], where
is the steady state gain coefficient of the laser and
is the noise enhancement factor compared to the ideal quantum noise level.
The final formula for the intensity noise variance and the noise auto-correlation function are listed in Equations (10) and (11) below. Here we have assumed the contribution from the initial noise can be ignored, and only the integrated added noises are taken into account. For the mode-locked laser problems considered in the present work, this is a valid assumption.
where in Equation (11),
and
are the two backpropagation solutions of the adjoint system corresponding to
and
respectively.
The general calculation procedures can be summarized in 3 steps. First of all, we forwardly propagate the nonlinear Equation (1) until the steady state of the mean field is reached. Secondly, we backwardly propagate the linear adjoint Equation (6) with the initial condition given by
at the end point of the first step. Thirdly, we then utilize the Equations (10) and (11) to determine the noise variance and the noise auto-correlation function with the solution of the adjoint equation. The noise spectrum can then be computed by Fourier transforming the noise auto-correlation function [
1,
4]. The finite difference Crank-Nicholson method is used for propagating all the nonlinear and linearized equations. The whole calculation is deterministic, and thus, the required computation efforts are much reduced when compared to those required for the stochastic simulation method.
3. Results
The typical output characteristics of the experimental fiber laser are summarized in
Table 1 and illustrated in
Figure 2. The laser can be operated under the single-soliton state at lower power levels, and under the bound-soliton state at higher power levels. In the middle of the transition, the laser is not very stable, mainly due to the co-existence of the two states. The measured pulse widths are in the ps level for both states, and the maximum output averaged power is near 60 mW.
Figure 2a shows the auto-correlation trace for the single-soliton state. The red dashed line is the fitting curve with an assumed hyperbolic secant pulse shape. The insets of the sub-figures indicate the corresponding optical spectra. In the inset of
Figure 2b, one can observe an interference dip located at the center of the spectrum, which indicates that the two bound-pulses are with π phase difference.
Figure 2c,d show the RF spectral measurement results through a high speed photo-detector for both the single- and bound-soliton states. The repetition rate is around 10 GHz and the laser is with a good super mode suppression ratio (SMSR).
Figure 3 shows the measurement data for the laser intensity noise spectra with a span from 1 kHz to 500 kHz, and a resolution bandwidth of 10 Hz. In the measurement range from 1 kHz to 500 kHz, one can avoid the dominant pump-induced amplitude noises near direct current (DC) so as to more clearly observe the intrinsic pulse intensity fluctuations of the mode-locked fiber laser [
13,
14,
15]. The measurement results approach the background floor at frequencies higher than 500 kHz. One can observe a flat intensity noise spectrum for both states from 1 kHz to the turning point near 100 kHz. The noise of the bound-soliton state is lower than that of the single-soliton state. The averaged noise power levels for the single- and bound- soliton states within the span from 1 kHz to 30 kHz are −96 dBm and −101 dBm respectively, and the total integrated power within the span from 1 kHz to 500 kHz are −36.3 dBm and −38.3 dBm. One can observe some oscillation peaks near the turning point for both states, and the flat bandwidth of the bound-soliton state is wider than that of the single-soliton state. One possible cause for these peaks is the relaxation oscillation effect caused by the gain dynamics inside the cavity [
6,
35]. However, as we will see in the theoretical results below, parts of these resonance peaks may also be attributed to the mode-locking dynamics. Note that, experimentally, we cannot operate the single- and bound-soliton states under the same power level, because the laser becomes unstable due to the co-existence of the two states. Therefore, for direct experimental comparison, the measurement results for both states are performed with the same averaged optical power incident on the photo-detector through the use of an optical attenuator.
For numerically calculating the intensity noises, the physical normalization units and the various cavity parameters used in our simulation are given below. The normalization unit of
is 0.5 ps, and the unit of
is the cavity roundtrip time, which is about 200 ns, as determined from the cavity length. The normalization unit of the pule energy is 10 pJ, which is estimated from the experimental results. The group velocity dispersion parameter
is estimated from the experimental condition. We set the SPM parameter at
according to the soliton assumption, and finely adjust it to make the numerical and experimental results better matched. The saturable absorption parameter
is set to 0 because there was no passive saturable absorption in our laser due to the sigma-type cavity. The phase modulation frequency is set around 10 GHz, with a modulation depth around 1. The cavity linear gain is set to be variable from
= 2.5 to 5 to accord with the lower and higher power level cases according to the relation
[
25]. The noise enhancement factor
is set to be 10. Other parameters are also reasonably determined so that the simulation results can be closely matched with the experimental results. The initial condition of the single-soliton case is a symmetrically single secant-shape pulse, and the initial condition of the bound-soliton case is two anti-symmetrical secant-shape pulses with
phase difference. The total propagation distance in the simulation is 1000 roundtrips with the step size
dT = 0.005. The small-scale time step is
dt = 0.01 in the normalization units, and the total simulation window is 25 ps. We have checked that these step sizes are small enough to ensure the accuracy of the calculation.
Figure 4 shows the steady state propagation solutions of the single- and bound-soliton states in the case of
= 2.5 and 5. The pulse width for the single-soliton is around 2.41 ps, while the pulse width and timing separation of the bound-soliton state are around 1.72 and 4.58 ps. We gradually increase
from 2.5 to 5 to investigate the pulse properties under different pulse energies.
Figure 5a indicates the simulated pulse energy for both states versus
. Basically, a linear relation can be obtained to match with the relation
.
Figure 5b is the pulse width and bound-pulse timing separation as a function of
. The pulse width and timing separation are both decreased when
is increased, because of larger nonlinearity. The pulse width reduction is extremely sensitive in the single-soliton state; this is because the single-soliton is with a larger peak power to induce larger nonlinearity.
Figure 5c shows the required steady-state lasing gain for the considered mode-locking states by analyzing the Master equation as an eigenvalue problem [
25,
36]. The negative real part of the solved eigenvalue exactly corresponds to the required steady-state lasing gain for the corresponding mode-locking state in the presence of all the cavity effects. Physically, the state with the minimum steady-state lasing gain is the most stable solution that will emerge as the final steady state of the system. This is because once the state emerges, other states will see net losses, and thus, decay away. This is why the steady state lasing gain can be used as a meaningful physical quantity for determining the stability of the lasing states. Nevertheless, if there are several states with enough close steady-state lasing gains, the perturbations in the laser can cause state jumping, and the differential losses may be too small to make the perturbations decay quickly. Under this situation the laser will become unstable. One can see that the required steady state lasing gains for the single- and bound-soliton states are very close in the overlap region marked by the blue circle. The laser operation within this region is expected to be not very stable experimentally, due to the co-existence of two states with close steady state lasing gains, since small fluctuations of the laser will easily lead to state jumping under this condition. This explains why, experimentally, we can only operate the laser under the single-soliton state in the lower power level, and under the bound-soliton state in the higher power level.
Figure 6a,b are the theoretical single sideband RIN spectra, derived by taking the Fourier transform of the noise auto-correlation functions calculated by the linearized backpropagation method for both the single- and bound-soliton states respectively. The main difference of the noise spectra between the two states is the level of the flat spectral distribution from 1 kHz to the turning point near 100 kHz. The bound-soliton state has a lower noise level, which is in good agreement with the observed experimental tendency. Moreover, the flat bandwidth of the bound-soliton state is wider than that of the single-soliton state, i.e., about tens of kHz. This again agrees with the experiment. Some resonance peaks can also be resolved near the turning point for both states, which can be attributed to the coupling effects of different pulse parameters [
4]. Since the relaxation oscillation effects of the Er-doped fiber lasers may also show up in this range of frequencies, one may not be able to resolve these resonance peaks experimentally.
Figure 6c shows the integrated noise variance within the span from 1 kHz to 2.5 MHz versus different values of
. The integrated RIN noise variance for the single- and bound-soliton states are
and
at
= 2.5 and
= 5 respectively, which corresponds to the energy fluctuations of 0.087% and 0.052% respectively. Both the single- and bound-soliton RIN variances decrease when the pulse energy gets larger. The bound-soliton states have smaller RIN noises compared to the single-soliton states of the same laser by following the similar lasing power scaling tendency as shown in
Figure 6c.