Abstract
We report recent results obtained with a novel optical fiber experimental setup based on a heterodyne optical time-domain reflectometer in the context of FPU recurrence process. Moreover, we actively compensate the dissipation of the system. We show that we can observe several FPU recurrences by monitoring the power and relative phase evolutions of the main discrete frequency components involved in the process.
1 Introduction
The Fermi-Pasta-Ulam (FPU) recurrence process describes the ability of a multimodal nonlinear system to come back to its initial state after its energy had been distributed over its modes during the evolution. It had been discovered in the 50s by Fermi and co-workers who were studying the dynamics of a chain of nonlinearly coupled oscillators [1]. The problem has been extensively studied and led to the development of a theory allowing to explain the formation of solitons by Zabusky and Kruskal in 1965 [2]. Since then, the FPU recurrence process has been especially investigated for its role in chaos [3] and for its relation with the nonlinear stage of modulation instability (MI), in the framework of nonlinear systems which can be described by the nonlinear Schrödinger equation (NLSE). In particular, it led to the development of breathers theory [4], which plays an important role in the description of rogue waves formation [5–10]. Experimentally, FPU had been observed in several fields of physics such as hydrodynamics [11] and electric transmission lines [12].
The first observation of FPU recurrence in optics has been achieved in 2001 by Van Simaeys et al. [13]. The main issue was to monitor the longitudinal evolution of the optical power along the fiber to confirm that the system indeed came back to its initial state. Van Simaeys et al. [13] succeeded via an indirect measurement method. These results attracted a lot of attention from the nonlinear scientific community with the publication of a short editorial in the prestigious Nature journal [14]. Indeed, optical fibers constitute a fantastic test bed for investigations of such complex nonlinear processes as they offer an easy tuning of initial conditions, an accurate control of fiber parameters and an accurate characterization of the output field with sub-ten femtosecond resolution. For those reasons, this pioneering result has motivated numerous studies of FPU, breathers or rogue wave formation in fiber optics [8, 10, 15–21]. This non-exhaustive list of works highlights the strong interest for this experimental system. However, there are two main limitations in fiber optics systems. Firstly, the monitoring of the longitudinal evolution of the field is not straightforward while it is crucial in FPU recurrence investigations to unambiguously prove that the systems evolved before coming back to its initial state. Van Simaens et al. exploited the invariant properties of the NLSE. They used a fiber of fixed length and performed indirect distributed measurements by varying the optical power instead of the distance z. This method was also used for the first observation of Peregrine soliton [16]. However, varying also changes the MI frequency which scales in so the frequency detuning of the seed has also to be varied to keep a normalized frequency detuning constant, which is not really convenient. Another set of solutions based on cut-back techniques had been reported. They are either fully destructive by cutting the fiber in several pieces [15] or by using fibers of different lengths [19] and repeating the experiments by copying the electric field characteristics to the input of the following fiber. Another approach to perform distributed measurements is the use of scattering processes, either elastic (Rayleigh) or inelastic (Brillouin), which presents the advantage of being non-invasive and offers a good resolution [22]. As an example, Brillouin optical time-domain reflectometry (OTDR) has been used [23, 24] to characterize fiber optical parametric amplifiers. Concerning OTDR based on Rayleigh scattering, it enabled the distributed measurement of one FPU recurrence via the distributed characterization of the pump and low-orders sidebands intensity when the process is coherently seeded [25] or the pump and the amplified noise intensity when it is spontaneous [26]. These options are very interesting but have the drawback to only measure power evolutions. We will see below that the knowledge of the whole electric field is crucial to accurately characterize the FPU recurrence process. The second limitation is due to dissipation, despite the ultra-low attenuation of optical fibers with the latest record of 0.1419 dB/km [27]. This is the reason why in previous experiments on FPU recurrence, most recordings had been limited to one or one and a half period of recurrence [13, 15, 23, 25, 28]. The capability to observe more than one recurrence period is of great interest to study in detail the complex dynamics of the FPU process appearing at the second compression point [29]. There is a need to develop an experimental optical fiber system enabling the monitoring of the longitudinal evolution of the electric field (phase and amplitude) without loss. Note that multiple recurrences had been observed in bulk optical systems [30], which can be considered as lossless untill a certain point, but initial conditions are trickier to tune as compared to fiber optics.
In this paper we present a review of recent results obtained with a novel optical fiber experimental setup allowing a non-invasive distributed characterization of the amplitude and phase of the main frequency components of a laser pulse train based on a heterodyne optical time domain reflectometer (HOTDR) [31]. Moreover, we implemented an active loss compensation scheme by exploiting the Raman amplification effect to get an almost fully transparent fiber [31]. Note that other groups took benefit of a recirculating loop to observe two recurrences as it was originally done in telecommunications applications to emulate long-haul fiber systems [32].
The paper is organized as follows. In the first part, we present the experimental setup. We provide detailed information on the HOTDR technique and on the active loss compensation compared to already published papers. Then we review various striking results that we recently published [33–38]. We first show a detailed characterization of the FPU recurrence, highlighting the symmetry breaking of the process [33, 35], demonstrate a record of 4 recurrences in an ultra-low loss fiber [37] and reveal the observation of A- and B-type solutions of the NLSE [36]. Then we show that a spatio-temporal evolution of the whole electrical field can be extracted from HOTDR recordings [34]. We illustrate this with the fine characterization of first and second order breathers with sub-ps resolution.
2 Basic Concepts
MI results from the amplification of weak sidebands by a strong pump wave. A numerical example is provided in Figure 1. As initial conditions, we used a strong pump ( mW), surrounded by two symmetric weak sidebands (20 dB below) located at GHz. This frequency shift corresponds to the maximum MI gain (Figure 1A). These three waves are launched in the anomalous region of an optical fiber. In the early steps, the sidebands are amplified at the expense of the pump and rapidly additional side bands resulting from multiple cascaded four wave mixing processes are generated during the propagation. At about 5 km, more than twelve of these Fourier modes are generated (Figure 1B) and the pump is maximally depleted. This broadest spectrum corresponds to maximally compressed pulses in the temporal domain (Figure 1C). Then, due to the linear and nonlinear phase acquired by these waves, the energy flow is reversed back into the pump. At about 9 km, the system presents characteristics which are almost identical to initial conditions.
FIGURE 1
This behaviour is a manifestation of the FPU recurrence process. By further propagating in the fiber, the process can be repeated several times. In experiments we will aim at monitoring the evolution of the pump and sideband power evolution as well as their relative phase to characterize in detail the FPU recurrence process.
3 Experimental Setup
The full experimental setup of the HOTDR we have developed is shown in Figure 2A. It is relatively complex and cumbersome as seen in the laboratory picture in Figure 2B. In the following, we will explain each element step by step.
FIGURE 2
3.1 Initial Conditions
The initial conditions we use to observe the FPU process consist of a pump and a pair of sidebands. To generate such a three-wave input, we first modulate the intensity of the CW beam delivered by laser 1 in order to get 50 ns square-shaped pulses (Figure 3A) at a 9.8 KHz (102 ms) repetition rate (it has to be decreased when fiber length exceeds ). To do so, an arbitrary waveform generator produces a clock at 9.8 kHz to trigger a pulse generator that delivers the RF signal driving the intensity modulator (IM1). Then, the optical pulses pass through a phase modulator (PM) to generate a frequency comb with a 35 GHz line spacing. A programmable optical filter (waveshaper) with a bandwidth of 10 GHz allows us to tailor the frequency comb to select the pump, signal and idler waves and to control their relative phases. Figure 3B shows the spectrum at the output of the waveshaper when we do not apply any filters (blue line) and when a filter is applied to get a symmetric three waves input (red line). The unwanted harmonics are filtered out. Following the waveshaper, the pulses are amplified via an erbium-doped fiber amplifier (EDFA 1) to reach the high peak power required to trigger the FPU recurrence process (usually around 450 mW). Then the signal goes through an acousto-optic modulator with a high extinction ratio (typically higher than 50 dB) to reduce the CW background between pulses in order to mitigate the stimulated Brillouin scattering (SBS). Finally, a filter with a 1 nm bandwidth (Filter 1) removes the excess amplified spontaneous emission. The 4 GHz OSA resolution, allows us to get a signal to noise ratio of about 50 dB for the pump and 40 dB for the signal and idler (see Figure 3B). Finally, short squared pulses composed of a strong pump surrounded by weaker signal and idler waves with adjustable amplitude and relative phase are generated.
FIGURE 3
3.2 Local Oscillator and Phase-Locking
In order to perform phase measurements via heterodyning, the local oscillator (laser 2) has to be initially phase-locked with laser 1. To do so, we used the method described in ref. 39 based on laser difference-frequency. A 5 GHz bandwidth photodiode receives the beat note of laser 1 and 2 which is then mixed with a 800 MHz RF reference signal provided by a stable synthesizer. In order to have a beat note close to 800 MHz, the frequency of laser 2 can be coarsely tuned via its power supply. The resulting intermediate frequency goes then through a PID controller (with a response time less than 15 ns) which drives the phase of the local oscillator. The efficiency of the phase-locking can be monitored by looking at the beat signal of laser 1 and 2 on a signal analyser as shown in Figures 4A,B, which display the RF spectrum (with a 100 Hz resolution) of the beat signal when the locking is turned off and on, respectively. Once the lasers are phase-locked, the locking can last up to several hours. As we intend to perform measurements not only on the pump component but also on the signal/idler one, the laser 2 is modulated in intensity by IM2 to create sidebands spaced by 35 GHz so that each frequency component of interest has its own local oscillator. The RF signal driving IM2 is delivered by the same microwave source that is driving the previous phase modulator in order to ensure a fixed phase relation between them. The local oscillator is then amplified by a semiconductor optical amplifier (SOA) and the amplitude of its three main components is roughly equalized by tuning the DC bias voltage of IM2 to get similar local oscillators, as can be seen in Figure 4C, which shows the spectrum of the local oscillator after the SOA.
FIGURE 4
3.3 Active Loss Compensation
Linear losses have a strong influence on the FPU recurrence as illustrated in Figure 5. These figures display the pump and signal power evolutions for an initial signal to pump ratio of dB and for (Figure 5A) or (Figure 5B) (with and the signal and pump phase respectively). We chose these specific value because they trigger the system in its two main operating regimes described by inner () or outer ) trajectories in the phase plane. In both cases, one can see that dissipation strongly influences the dynamic of the process. It reduces the amplitude of both pump and signal waves (Figure 5B) or completely disable the formation of the second recurrence (Figure 5A). The impact of the attenuation on the recurrent process is more or less pronounced depending on the initial relative phase, which is intrinsically linked to the complex nonlinear dynamics of the system. It is thus necessary to compensate losses to get an almost fully transparent optical fiber to observe more than two recurrences in the two main operating mode of the system. With perfect loss compensation, the dynamics of the process is ruled by the pure NLSE and then direct comparison between experiments and theory will be achievable. Hence, we introduced a Raman amplification scheme. This idea, borrowed from the telecommunication field where distributed Raman amplification is commonly used [40], has been recently implemented successfully in a handful of nonlinear fiber optics experiments [31, 34, 41–43]. The Raman laser source is located at about 13 THz from the signal to benefit from the maximum Raman gain amplification. It counter-propagates in order to minimize the relative intensity noise transfer compared to the co-propagative case [44]. The Raman pump power is set empirically as follows. We chose the configuration where and increased the Raman pump power until the level of signal sideband at the second peak of conversion was similar to the first one. This allows to get the best overall compensation of the loss along the fiber length. It lead to very similar input and output signal power values. The efficiency of the loss compensation will be illustrated hereafter, with experimental results in excellent agreements with numerics from pure NLSE.
FIGURE 5
3.4 Fading Suppression Technique and Signal Processing
In OTDR, it is known that a random noise in amplitude and phase is superimposed on the backscattered signal [45]. This noise, which is detrimental as it leads to jagged measurements, finds its origins in the random state of polarization of the backscattered light and a speckle-like phenomenon due to huge number of scattered waves involved in the process and/or in the thermo-mechanical fluctuations of the propagating medium. In order to avoid it, many techniques have been developed such as averaging over a huge number of backscattered signal and polarization scrambling [46, 47]. However, none of these techniques have proven to be effective and compatible with accurate distributed phase measurements in our setup. Thus, we have developed a novel method which allows us to remove effectively this additional source of noise. We made two assumptions: 1) the fading phenomenon is a purely linear effect and 2) the thermo-mechanical fluctuations have a characteristic time of a few milliseconds. Therefore, we propose the following post-processing treatment to suppress the fading. We first launch a strong 50 ns pulse into the fiber which will experience both linear and nonlinear effects. Then, 102 μs later, a weak pulse, attenuated by 13 dB, follows experiencing mostly linear effects. The attenuation is obtained via the acousto-optic modulator which is driven at half-rate clock (4.9 KHz) so that one pulse over two is attenuated as it is displayed in the schematic in Figure 6A. Note that the 102 μs delay between two consecutive pulses is long enough to avoid any overlapping between two backscattered signals but very short compared to the response time of the thermo-mechanical fluctuations which ensures that the linear fading effects experienced by both backscattered waves (from the strong and weak pulses) are strongly correlated. The output spectra at the input and output of the fiber of weak and strong pulses are displayed in Figures 6B,C. We observe that for strong pulses (magenta lines), the spectrum is composed of 11 waves in this typical example corresponding to the broadest spectrum achievable by increasing the pump power to get the first maximum compression point at the output of the fiber. In the opposite case, we observe a spectrum made of 3 main waves (cyan lines) plus weak harmonics 25 dB lower than signal/idler wave. We can thus consider that these weak pulses mainly experience linear effects. Note that we analyze each state of polarization of backscattered signals independently (which are later recombined in post-processing) thanks to the use of polarization beam splitters, thus avoiding the need of polarization scrambling. Once the heterodyned signals are logged with the four channels of the oscilloscope (two channels for each polarization states of the pump and two for those of the signal), we proceed to the following data processing. We will refer to each traces observed on the oscilloscope as with z being the distance along the fiber (obtained simply by converting the time of flight of the backscattered wave into a distance), can either be (for the nonlinear pulse) or (for the reference), is either 1 or 2 for each polarization state and refers to the or . The ratios and are calculated (the procedure is similar for the signal sideband) in order to remove all linear contributions and to cancel the fading effect in amplitude and phase1. Then, we used a short time fast Fourier transform (FFT) to extract phase and intensity value for each recordings. All demodulated traces are then averaged over at least one hundred shots to remove the noise and filtered afterward, by two types of filters: a median filter and a Savitzky-Golay filter.
FIGURE 6
The performances of our measurement technique are illustrated in Figure 7. Figures 7A,E show a single-shot of the pump and signal power evolutions (normalized to their respective maxima) and their relative phase without calibration with the reference (weak pulse). One can distinguish the exchange of energy between the pump and the signal but the measurement is very noisy and the relative phase evolution does not agree at all with the one expected. Figures 7B,F show the same measurements but this time with the calibration (fading removed). Although the measurements are still very noisy, one can notice almost two recurrences by looking at the power evolutions and a specific relative phase evolution with jumps around 2 km and 6 km which are approximately the positions where the direction of the power flow between the pump and the sideband reverses. After averaging over one hundred shots, we observe a significant reduction of the noise as displayed in Figures 7C,G. We then obtain clean traces after applying a Savitzky–Golay filter as shown in Figures 7D,H. All results presented in this paper passed this data processing procedure. From this experimental setup and this data processing, we can record clean phase and intensity evolutions of pump and signal waves along the fiber length and to almost perfectly compensate attenuation. In the following, we are going to present a selection of results that we obtained with our HOTDR system, which were not accessible with other setups.
FIGURE 7
4 Fermi-Pasta-Ulam Recurrence Observation
4.1 Evolution as a Function of the Initial Relative Phase
We first present the impact of the initial relative phase value on the dynamics of the FPU recurrence process in Figure 8. The experimental parameters are detailed in Figure 8 caption. Firstly, we looked for the initial phase value allowing the system to evolve as close as possible to the separatrix (homoclinic orbit). This curve marks the border between the two main regimes of the system which consists of the inner or outer trajectories in the phase plane [29, 33, 48]. Theoretically, it corresponds to an input relative phase between the pump and the signal from the truncated three wave model [48]. Note that, for weak input signals, the separatrix corresponds to an Akhmediev breather [4]. By finally tuning the relative phase in experiments, the closest results we got were obtained for (Figures 8A,B), very close from the theoretical value of .
FIGURE 8
A single cycle of recurrence can be seen corresponding to a single loop in the right half side of the phase plane, as expected. Note that, a slight change of makes the system switch o the inner or outer trajectories [33]. Then, we set the input relative phase to . Signal and pump power evolutions are shown in Figure 8C and the corresponding trajectory in the phase plane in Figure 8D (solid lines). As seen, the system experiences two recurrences. The pump extrema coincide with the signal ones, as expected from the theory [28, 48]. The curve in the phase plane occupies the right and left sides, showing a symmetry breaking of the FPU recurrence process due to the crossing of the semi-plane [33]. These results are in excellent agreement with numerical simulations depicted in dashed lines. Next, the relative initial phase is set to to trigger the symmetric case. It is illustrated in Figures 8E,F. Two recurrences are still visible, but do not experience π shift as can be seen in Figure 8F, where only the right plane is occupied.
The global evolution is summarized in Figure 8G. It shows the evolution of the position of the first and second maximum compression points as a function of the initial relative phase value . We compared our experimental results (crosses), with theoretical predictions from Santini and Grinevitch [49] (solid lines) and numerics from the NLSE (dashed lines). The region for which the dynamics of the system corresponds to an outer (inner) trajectory is colored in shaded green (violet) areas, the borders correspond to the initial phase value to trigger the separatrix of the system.
As we can see, the positions of the first and the second maximum compression points diverge both for . This is because this particular phase corresponds to the stable manifold of the separatrix, which entails asymptotic conversion from the input sidebands to the pump. Conversely, the position of the first maximum compression point is minimum for , which entails the most rapid growth of the sidebands to the apex. Note that the position of the second maximum compression point still diverges for this phase value since backconversion still occurs asymptotically along the separatrix. The phase being the critical parameter of the system, influences the transition from inner to outer trajectory across this phase. The relatively strong modulation used in our experiment and NLSE simulations (the sidebands contains about 15% of the total power) induces a shift of the critical phase which is here . The measured values of the position of the first and second recurrence are in excellent agreement with the numerical simulations and theoretical predictions [49]. Finally, our system allowed us to finely characterize the dynamics of FPU recurrence in a nonlinear system modeled by the pure NLSE and to observe two maximum compression points. More specifically, we evidenced the two main operating regimes of the system characterized by inner or outer trajectories as well as the regime very close to the one modeled by the homoclinic orbit. Detailed information on the symmetry breaking of the FPU recurrence process can be found in refs. 33 and 35.
4.2 Record Observation of 4 Recurrences in an Ultra-Low Loss Optical Fiber
In the previous section, we reported the observation of two recurrences thanks to an active compensation of the loss using a counter-propagating Raman pump. Nonetheless, much remains to be done, since at present it is not clear how far one can go in observing a higher number of recurrences. This is a preliminary step for further studies on the onset of thermalization [50] for instance. The recirculating loop architecture reported by Goosens et al. [32] seems to be an other option, but reported results are still limited to two recurrences. Here we show that by improving our setup it is possible to observe more than 4 recurrences. To do so, we used an ultra-low loss optical fiber developed for ultra-long-haul submarine telecommunication systems [27] with 0.147 dB/km at 1550 nm (see details in ref. 37).
Figure 9 shows an example of recordings. We set the input phase to . This leads to phase-shifted recurrences [48] (outer trajectories). In Figure 9A, the power evolutions along the fiber length are represented in blue solid line for the pump and red solid line for the signal, respectively. The phase evolution is shown in Figure 9B. Numerics are represented in dashed lines. As seen, four recurrences have been observed, that is a record achievement [37]. These four recurrences correspond to ten nonlinear lengths in quasi-lossless optical fiber. The excellent agreement between experiments and simulations proves that the active compensation is not far from being perfect. Note that we also reported the observation of more than four recurrences in ref. 37 by triggering the regime characterized by inner trajectories (). However, the system does not spatially evolve in a purely periodic fashion, but is rather following orbits in the infinite dimensional space that can be closer to a higher-order solution of the nonlinear Schrödinger equation [51]. It indicates that the input sideband power was a bit too high to remain close to the same trajectory. These results have been made possible due to the lowering of losses at 1550 nm by 0.05 dB/km from one side but also at the Raman pump wavelength (1450 nm) with an improvement of 0.08 dB/km. These results pave the way to study further complex problems such as the role of multiple unstable sideband pairs and the route to the thermalization using longer fiber spans.
FIGURE 9
4.3 Observation of Type-A and B Solutions of the Nonlinear Schrödinger Equation
Akhmediev breathers are exact solutions of NLSE that allow to predict for a single stage of growth to the apex and the asymptotic return to the initial state [4]. They do not allow to predict the formation of multiple FPU recurrences. Approximate solutions have been reported recently by means of finite-gap theory and asymptotic expansions to describe multiple recurrences [49]. Exact solutions have been derived by Akhmediev et al. and recently revisited to extract their Fourier coefficients [38]. There are two types of first-order solutions, termed type-A and type-B, according to their location in the infinite-dimensional phase space with respect to the homoclinic orbit corresponding to the Akhmediev breather trajectory [51, 52]. Here we show that they can be observed experimentally with our experimental setup for the first time. We triggered both type-A and type-B solutions by using three initial waves and setting their relative phases to or 0, respectively [36].
First we set to trigger B-type solutions. Experimental results are displayed in Figures 10A–C in solid lines. Two maximum compression points are observed. The corresponding trajectory is limited in the right half side of the phase plane in Figure 10C as expected from theory [4] for B-type solutions. These results are in excellent agreement with analytical solutions (Eqs. 13–16 in ref. 4 whose parameters were determined through the perfect match condition between the initial three central waves and the experimental three-wave input). A-type solutions are then triggered by setting the initial phase to . Experimental results are represented in Figures 10D–F. This time, the first maximum point appears in the right half side of the phase plane and the second one in the left half side. This π shift is characteristic of A-type solutions [4]. These results are also confirmed by exact solutions (dashed lines, Eqs. 13–16 in ref. 4). Detailed information can be found in ref. 36. These results highlight the fact that this system might be a fantastic test bed to observe more complex solutions, such as the rogue waves on the doubly periodic background recently derived in ref. 53.
FIGURE 10
4.4 Spatio-Temporal Characterization in Intensity and Phase
4.4.1 Principle
We saw that the measurement the phase and the intensity evolutions of different laser lines in the spectral domain was insightful and allowed the observation of new physical phenomena. However, getting the spatio-temporal evolution of the electric field would be more meaningful, providing an observation of FPU recurrence in space. This requires phase and intensity measurements with sub-picosecond resolution. Standard techniques such as frequency-resolved optical gating [54] and spectral phase interferometry for direct-electric field reconstruction [55] (a review can be found in ref. 56) provide these information with an excellent temporal resolution. Recently, new methods allowing real-time full-field characterization have been developed, based in particular on advanced time lens [21] and/or standard time lens combined with a dispersive Fourier transform technique [57] with typically sub-hundred’s femtosecond resolution. However, these techniques [21, 54–58] are limited to localized measurements. Note that Xu [59] reported recently the longitudinal full-field characterization of Peregrine-like structures using a combination of temporal and spectral measurements and a reconstruction algorithm but this technique suffers from the same limitations as the cut-back one. Based on our HOTDR setup, we implemented a new method to perform fast and non-invasive full-field characterization of time-periodic pulses (i.e., with discrete line spectra) along the whole length of an optical fiber [34]. The idea is to calculate the inverse Fourier transform of the data recorded with the HOTDR setup to get the temporal characteristics of the electric field. In most investigations presented in this paper, we measured a truncated spectrum (three main spectral components, i. e., the pump and signal/idler) since it was enough to get a clear description of the process. However, we remind that at maximum compression points, more than eleven spectral components appear in the spectrum in our configurations [33]. We estimated from numerics that the truncation lowers the accuracy of the recordings, specifically in regions where sharp transitions exist. By calculating the error in the sense of the mean square difference from numerics, we estimated that the error is about by considering 3 waves, drops down to with 5 waves and to less than with 7 waves [34]. For practical reasons, we limit our recordings to 5 waves, but there is no fundamental limitations to make it with a larger number and to obtain almost perfect spatio-temporal recordings with less than 1 error [34].
4.4.2 First Order Breathers
We show here the spatio-temporal evolution of the FPU recurrence process presented in previous paragraphs. The detection of the five waves requires launching five waves at the fiber input in order to improve the signal-to-noize ratio in the heterodyne detection stage. So we launched two additional waves compared to previous experiments, but we checked numerically that they are so weak (20 dB lower than the pump), that they do not affect the dynamics of the system. We only record the evolution of the pump, signal and its first harmonic as displayed in Figures 11A,B. It shows the longitudinal power (Figure 11A) and phase evolutions (Figure 11B) in solid lines for experiments. The system exhibits two recurrences, in good agreement with numerical simulations depicted in dashed lines. We then calculate the spatio-temporal evolution of the field via the inverse Fourier transform from the five waves. The intensity plot is displayed in Figure 11C and exhibits two maximum compression points at about 1.5 km and 6 km, respectively, which are in phase, as expected (here ). This is in good overall agreement with the results obtained from numerical simulations accounting for the full spectrum and displayed in Figure 11D. The experimental spatio-temporal 2D-map shows broader pattern maximum compression points due to the truncation [34].
FIGURE 11
The spatio-temporal evolution of the phase of the field is displayed in Figures 11E,F. The experimental evolution is in very good agreement with numerical simulations, the impact of the truncation being less pronounced on the phase evolution. Noteworthy, in the temporal case (Figure 11G), we observe a phase jump close to π between the center of the pulses and their wings and equal to π in simulations which is one of the characteristic features of Akhmediev breathers and Peregrine soliton [8]. More details can be found in ref. 34 where we also reported the spatio-temporal evolution of the FPU recurrence when the symmetry is broken. Finally, we demonstrated that our HOTDR system enables non-invasive distributed measurements of the electric field in the temporal domain with a resolution lying in the picosecond scale with a 20 m spatial resolution.
4.4.3 Second Order Breathers
Second order breathers is a family of higher order solutions of the NLSE [60]. They have already been observed in fiber optics systems by Frisquet et al. [18]. They are known to experience sharper edges and higher maximum powers compared to first-order solutions (e.g., the maximal normalized intensity of a first order breather is 9 while it is 25 for a second-order one). Thus, they are of great interest because they could constitute good prototypes to model extremely rare and powerful rogue waves. We will show that it is possible to record the spatio-temporal evolution of their intensity and phase with our experimental setup. The solution of the field envelope of a second-order breather is reported in ref. 60. We use the same normalized parameters in the following. In order to design the experimental initial conditions to trigger second-order breathers, we started from the analytical expression of these breathers [60] and truncated this solution to the 5 central components. We then use these 5 waves to trigger the formation of second order breathers.
We checked numerically (not shown here) that these solutions are very robust and while being excited by a truncated initial condition, they remain quasi-unchanged. We set the phase and amplitude in experiments by using the waveshaper (Figure 2).
As an example, we choose a solution for which signal and harmonics experience the same linear MI gain: (and ). In physical units, this corresponds to and in normalized units from ref. 60. For a total power of mW, s/m and W/m this corresponds to GHz and GHz. With these parameters, the maximum compression point of this second order breather is predicted at km. The Fourier transform of this solution [] is calculated in order to get the input spectrum. This spectrum is displayed in red in Figure 12A. As we can see, it is composed of many sidebands with only four falling within the MI gain band (the MI gain curve is superimposed in blue): the signal and idler at frequencies and their first harmonics at . These bands have indeed experience similar parametric gain. We extracted the phase and amplitude value of the 5 main components and use them as initial conditions in experiments. The input initial spectrum is shown in Figure 12B. Experimental results are depicted in Figures 13C,D,G,H. The evolution of the pump, signal and harmonic power (Figure 13C) are in excellent agreement with numerical simulations (solid lines in Figure 13A), triggered with the truncated 5 waves solution, as well as the analytical solution (dashed black lines in Figure 13A). The relative phase evolution is shown in Figure 13B (numerics and analytical solutions) and in Figure 13D (experiments). The agreement between experiments, numerics and analytical predictions is excellent too. It means that while initial conditions does not correspond exactly to second order solutions, it is still possible to trigger this solution. By calculating the inverse Fourier transform of these 5 waves it is possible to plot the spatio-temporal evolution of the power (Figure 13G) and the phase (Figure 13H). As seen, a sharp and powerful second-order breather can be observed. We obtained a pretty good agreement with numerical simulations (Figures 13E,G). The main disagreement at the maximum compression point is due to the truncation. We show that using the HOTDR setup and then calculating the inverse Fourier transform, it is possible to obtain the spatio-temporal evolution of the electric field of a second-order breather. An even better agreement would be possible by using more waves.
FIGURE 12
FIGURE 13
5 Conclusion
We have reviewed recent results obtained with a novel optical fiber experimental setup based on a heterodyne optical time-domain reflectometer allowing to monitor the power and relative phase evolution of the main discrete frequency components of a periodic laser pulse train. Moreover, we actively compensated fiber losses with a conter-propagating Raman pump leading to an almost fully transparent optical fiber which can be accurately modeled by the NLSE. We provided details about the experimental setup, notably about the multi-heterodyne detection technique, the phase-locking between the laser and the local oscillator, the solution we developed to suppress the fading effect, and the data processing. Then, we presented a selection of results recently obtained with this system in the context of the FPU recurrence process. It enabled the observation of several FPU recurrences for the first time, the study of their dynamics, and a direct comparison with analytical predictions from the pure NLSE. More precisely, the symmetry breaking of the FPU process has been observed by varying initial conditions, in excellent agreement with analytical predictions. Moreover, we showed by recording first harmonics in addition to the signal, idler, and pump waves it is possible to obtain the spatial evolution of the electric field of time-periodic pulses with a sub-ps resolution. We illustrated the performance of this technique by showing the power and phase evolution along the fiber of first and second-order breathers. All these results have been confirmed by numerical simulations, and remarkably with analytical solutions thanks to the efficient loss compensation scheme we developed. This original experimental setup combining dissipation compensation and non-invasive characterization of the electric field of a periodic laser pulse train is a powerful tool in nonlinear fiber optics experiments that could be implemented to investigate other nonlinear effects.
Statements
Author contributions
CN and GV performed experiments and numerical simulations. PS and AM designed the experimental setup with AK and GM. MD provided optimized optical fibers. NA, MC, and ST developed theoretical/numerical tools. All authors analyzed results and participated in the writing of the paper.
Funding
This work was partly supported by the Agence Nationale de la Recherche through the Equipex Fibres optiques pour les hauts flux (FLUX)through the “Programme Investissements d’Avenir”, by the Ministry of Higher Education and Research, Hauts de France Council and European Regional Development Fund (ERDF) through the Contrat de Projets Etat-Region (CPER Photonics for Society, P4S), HEAFISY project and I-SITE through the FUHNKC, VERIFICO and EXAT projects and H2020 Marie Skodowska-Curie Actions (MSCA)(713694) and MEFISTA.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2021.637812/full#supplementary-material.
Footnotes
1.^Unfortunately, this removed also the linear phase due to the group velocity dispersion acquired during the propagation, which cannot be neglected as it directly impacts the dynamics of the system. The following phase term is then added later on to the measured pump-signal phase difference to get .
References
1.
FermiEPastaPUlamSTsingouM. Studies of the nonlinear problems. In: Technical report. New Mexico: Los Alamos Scientific Lab. (1955).
2.
ZabuskyNJKruskalMD. Interaction of "solitons" in a collisionless plasma and the recurrence of initial states. Phys Rev Lett (1965) 15(6):240. 10.1103/physrevlett.15.240
3.
ZabuskyNJ. Fermi-Pasta-Ulam, solitons and the fabric of nonlinear and computational science: history, synergetics, and visiometrics. Chaos (2005) 15(1):15102. 10.1063/1.1861554
4.
AkhmedievNNEleonskiiVMKulaginNE. Exact first-order solutions of the nonlinear Schrödinger equation. Theor Math Phys (1987) 72(2):809–18. 10.1007/bf01017105
5.
SolliDRRopersCKoonathPJalaliB. Optical rogue waves. Nature (2007) 450(7172):1054–57. 10.1038/nature06402
6.
AkhmedievNSoto-CrespoJMAnkiewiczA. Extreme waves that appear from nowhere: on the nature of rogue waves. Phys Lett (2009) 373(25):2137–45. 10.1016/j.physleta.2009.04.023
7.
AkhmedievNDudleyJMSolliDRTuritsynSK. Recent progress in investigating optical rogue waves. J Optic (2013) 15(6):060201. 10.1088/2040-8978/15/6/060201
8.
DudleyJMDiasFErkintaloMGentyG. Instabilities, breathers and rogue waves in optics. Nat Photon (2014) 8(10):755–64. 10.1038/nphoton.2014.220
9.
AkhmedievNKiblerBBaronioFBelićMZhongWPZhangYet alRoadmap on optical rogue waves and extreme events. J Opt (2016) 18(6):063001. 10.1088/2040-8978/18/6/063001
10.
DudleyJMGentyGMussotAChabchoubADiasF. Rogue waves and analogies in optics and oceanography. Nat Rev Phys1(11):675–89. 10.1038/s42254-019-0100-0
11.
LakeBMYuenHCRungaldierHFergusonWE. Nonlinear deep-water waves: theory and experiment. part 2. evolution of a continuous wave train. J Fluid Mech (1977) 83(1):49–74. 10.1017/S0022112077001037
12.
JägerD. Soliton propagation along periodic-loaded transmission line. Appl Phys (1978) 16(1):35–8. 10.1007/bf00931418
13.
Van SimaeysGEmplitPHaeltermanM. Experimental demonstration of the Fermi-Pasta-Ulam recurrence in a modulationally unstable optical wave. Phys Rev Lett (2001) 87(3):033902. 10.1103/PhysRevLett.87.033902
14.
AkhmedievNN. Nonlinear physics. Déjà vu in optics. Nature (2001) 413(6853):267–8. 10.1038/35095154
15.
MussotAKudlinskiADroquesMSzriftgiserNAkhmedievN. Fermi-Pasta-Ulam recurrence in nonlinear fiber optics: the role of reversible and irreversible losses. Phys Rev X (2014) 4(1):011054. 10.1103/physrevx.4.011054
16.
KiblerBFatomeJFinotCMillotGDiasFGentyGet alThe Peregrine soliton in nonlinear fibre optics. Nat Phys (2010) 6(10):790–5. 10.1038/nphys1740
17.
KiblerBFatomeJFinotCMillotGGentyGWetzelBet alObservation of Kuznetsov-Ma soliton dynamics in optical fibre. Sci Rep (2012) 2:463. 10.1038/srep00463
18.
FrisquetBKiblerBMillotG. Collision of Akhmediev breathers in nonlinear fiber optics. Phys Rev X (2013) 3(4):041032. 10.1103/physrevx.3.041032
19.
HammaniKWetzelBKiblerBFatomeJFinotCMillotGet alSpectral dynamics of modulation instability described using Akhmediev breather theory. Opt Lett (2011) 36(11):2140–2. 10.1364/OL.36.002140
20.
ErkintaloMHammaniKKiblerBFinotCAkhmedievNDudleyJMet alHigher-order modulation instability in nonlinear fiber optics. Phys Rev Lett (2011) 107(25):253901. 10.1103/PhysRevLett.107.253901
21.
TikanABielawskiSSzwajCRandouxSSuretP. Single-shot measurement of phase and amplitude by using a heterodyne time-lens system and ultrafast digital time-holography. Nat Photon (2018) 12(4):228–34. 10.1038/s41566-018-0113-8
22.
ThévenazL. Review and progress in distributed fiber sensing. In: Optical fiber sensors. Cancun, Mexico: Optical Society of America (2006). ThC1.
23.
ArmandVAlasiaDLantzEMaillotteHThévenazLGonzalez-HerraezMet alBrillouin optical time-domain analysis of fiber-optic parametric amplifiers. IEEE Photon Technol Lett (2007) 19(3):179–81. 10.1109/LPT.2006.890039
24.
AlishahiFVedadiAShoaieMASotoMADenisovAMehranyKet alPower evolution along phase-sensitive parametric amplifiers: an experimental survey. Opt Lett (2014) 39(21):6114–7. 10.1364/OL.39.006114
25.
HuXChenWLuYYuZChenMMengZ. Distributed measurement of fermi-pasta-ulam recurrence in optical fibers. IEEE Photon Technol Lett (2018) 30(1):47–50. 10.1109/lpt.2017.2773615
26.
MartinsHFMartin-LopezSCorrederaPSalgadoPFrazãoOGonzález-HerráezM. Modulation instability-induced fading in phase-sensitive optical time-domain reflectometry. Opt Lett (2013) 38(6):872–4. 10.1364/OL.38.000872
27.
TamuraYSakumaHMoritaKSuzukiMYamamotoYShimadaKet alLowest-ever 0.1419-db/km loss optical fiber. In: Optical fiber communication Conference, 2017 March 19-23, Los Angeles, CA. Optical Society of America (2017). Th5D–1.
28.
BendahmaneAMussotAKudlinskiASzriftgiser,PConfortiMWabnitzSet alOptimal frequency conversion in the nonlinear stage of modulation instability. Optic Express (2015) 23(24):30861–71. 10.1364/OE.23.030861
29.
MoonHT. Homoclinic crossings and pattern selection. Phys Rev Lett (1990) 64(4):412. 10.1103/PhysRevLett.64.412
30.
PierangeliDFlamminiMZhangLMarcucciGAgranatAJGrinevichPGet alObservation of Fermi-Pasta-Ulam-Tsingou recurrence and its exact dynamics. Phys Rev X (2018) 8(4):041017. 10.1103/physrevx.8.041017
31.
MussotANaveauCConfortiMKudlinskiACopieFSzriftgiserPet alFibre multi-wave mixing combs reveal the broken symmetry of Fermi-Pasta-Ulam recurrence. Nat Photon (2018) 12(5):303–8. 10.1038/s41566-018-0136-1
32.
GoossensJWHafermannHJaouënY. Experimental realization of fermi-pasta-ulam-tsingou recurrence in a long-haul optical fiber transmission system. Sci Rep (2019) 9(1):1–11. 10.1038/s41598-019-54825-4
33.
MussotAConfortiMTrilloSCopieFAlexandreK. Modulation instability in dispersion oscillating fibers. Adv Opt Photon (2018) 10(1):1–42. 10.1364/AOP.10.000001
34.
NaveauCSzriftgiser,PKudlinskiAConfortiMTrilloSMussotA. Full-field characterization of breather dynamics over the whole length of an optical fiber. Opt Lett (2019) 44(4):763–6. 10.1364/OL.44.000763
35.
NaveauCSzriftgiserPKudlinskiAConfortiMTrilloSMussotA. Experimental characterization of recurrences and separatrix crossing in modulation instability. Optic Lett (2019) 44(22):5426–9. 10.1364/ofc.2018.th3d.5
36.
VanderhaegenGSzriftgiserPNaveauCKudlinskiAConfortiMTrilloSet alObservation of doubly periodic solutions of the nonlinear schrödinger equation in optical fibers. Opt Lett45(13):3757–60. 10.1364/OL.394604
37.
VanderhaegenGSzriftgiserPKudlinskiAConfortiMTrilloSDroquesMet al2020. Observation of four fermi-pasta-ulam-tsingou recurrences in an ultra-low-loss optical fiber, Opt Express. 28(12):17773–81. 10.1364/OE.391560
38.
ConfortiMMussotAKudlinskiATrilloSAkhmedievN. Doubly periodic solutions of the focusing nonlinear Schrödinger equation: recurrence, period doubling, and amplification outside the conventional modulation-instability band. Phys Rev A (2020) 101:023843. 10.1103/physreva.101.023843
39.
FriederichFSchurichtGDeningerALisonFSpickermannGHaring BolívarPet alPhase-locking of the beat signal of two distributed-feedback diode lasers to oscillators working in the mhz to thz range. Optic Express (2010) 18(8):8621–9. 10.1364/OE.18.008621
40.
HeadleyCAgrawalG. Raman amplification in fiber optical communication systems. 1st ed. Academic Press (2005). 392.
41.
Juan Diego Ania-CastañónEllinghamTJIbbotsonRChenXZhangLTuritsynSK. Ultralong raman fiber lasers as virtually lossless optical media. Phys Rev Lett (2006) 96:023902. 10.1103/physrevlett.96.023902
42.
XuGConfortiMKudlinskiAMussotATrilloS. Dispersive dam-break flow of a photon fluid. Phys Rev Lett (2017) 118(25):254101. 10.1103/PhysRevLett.118.254101
43.
KraychAESuretPGennadyERandouxS. Nonlinear evolution of the locally induced modulational instability in fiber optics. Phys Rev Lett (2019) 122(5):054101. 10.1103/physrevlett.122.054101
44.
AgrawalGP. Nonlinear fiber optics. 5th ed. Academic Press (2012). 648.
45.
HealeyP. Fading in heterodyne OTDR. Electron Lett (1984) 20(1):30–2. 10.1049/el:19840022
46.
IzumitaHFurukawaSI.KoyamadaYSankawaI. Fading noise reduction in coherent OTDR. IEEE Photon Technol Lett (1992) 4(2):201–3. 10.1109/68.122361
47.
YangGFanXWangBLiuQHeZ. Polarization fading elimination in phase-extracted OTDR for distributed fiber-optic vibration sensing. In: 2016 21st OptoElectronics and Communications Conference (OECC) held jointly with 2016 International Conference on Photonics in Switching (PS), 2016 July 3-7. Niigata, Japan. IEEE (2016). pages 1–3.
48.
TrilloSWabnitzS. Dynamics of the nonlinear modulational instability in optical fibers. Opt Lett (1991) 16(13):986–8. 10.1364/ol.16.000986
49.
GrinevichPGSantiniPM. The exact rogue wave recurrence in the NLS periodic setting via matched asymptotic expansions, for 1 and 2 unstable modes. Phys Lett (2018) 382(14):973–9. 10.1016/j.physleta.2018.02.014
50.
WabnitzSAkhmedievN. Efficient modulation frequency doubling by induced modulation instability. Opt Commun (2010) 283(6):1152–4. 10.1016/j.optcom.2009.11.030
51.
KimmounOHsuHCBrangerHLiMSChenYYKharifCet alModulation instability and phase-shifted Fermi-Pasta-Ulam recurrence. Sci Rep (2016) 6:28516. 10.1038/srep28516
52.
AkhmedievNNKorneevVI. Modulation instability and periodic solutions of the nonlinear Schrödinger equation. Theor Math Phys (1986) 69(2):1089–93. 10.1007/bf01037866
53.
ChenJPelinovskyDEWhiteRE. Rogue waves on the double-periodic background in the focusing nonlinear Schrödinger equation. Phys Rev E (Nov 2019) 100:052219. 10.1103/PhysRevE.100.052219
54.
RickT. Frequency-resolved optical gating: the measurement of ultrashort laser pulses. Springer, Boston, MA: Springer Science & Business Media (2012).
55.
IaconisCWalmsleyIA. Spectral phase interferometry for direct electric-field reconstruction of ultrashort optical pulses. Opt Lett (1998) 23(10):792–4. 10.1364/ol.23.000792
56.
WalmsleyIADorrerC. Characterization of ultrashort electromagnetic pulses. Adv Optic Photon (2009) 1(2):308–437. 10.1364/aop.1.000308
57.
RyczkowskiPNärhiMBilletCMerollaJMGentyGDudleyJM. Real-time full-field characterization of transient dissipative soliton dynamics in a mode-locked laser. Nat Photon (2018) 12(4):221–7. 10.1038/s41566-018-0106-7
58.
ChengLGodaK. The complete optical oscilloscope. Nat Photon (2018) 12(4):190–1. 10.1038/s41566-018-0141-4
59.
XuGHammaniKChabchoubADudleyJMKiblerBFinotC, 2019. Phase evolution of peregrine-like breathers in optics and hydrodynamics, Phys Rev E. 99(1):012207. 10.1103/PhysRevE.99.012207
60.
KedzioraDJAnkiewiczAAkhmedievN. Second-order nonlinear Schrödinger equation breather solutions in the degenerate and rogue wave limits. Phys Rev E - Stat Nonlinear Soft Matter Phys (2012) 85(6):066601. 10.1103/PhysRevE.85.066601
61.
DevineNAnkiewiczAGentyGDudleyJMAkhmedievN. Recurrence phase shift in Fermi-Pasta-Ulam nonlinear dynamics. Phys Lett (2011) 375(46):4158–61. 10.1016/j.physleta.2011.10.006
Summary
Keywords
Fermi Pasta Ulam recurrence, optical fibers, modulation instability, breathers and rogue waves, heterodyne, Rayleig
Citation
Naveau C, Vanderhaegen G, Szriftgiser P, Martinelli G, Droques M, Kudlinski A, Conforti M, Trillo S, Akhmediev N and Mussot A (2021) Heterodyne Optical Time Domain Reflectometer Combined With Active Loss Compensation: A Practical Tool for Investigating Fermi Pasta Ulam Recurrence Process and Breathers Dynamics in Optical Fibers. Front. Phys. 9:637812. doi: 10.3389/fphy.2021.637812
Received
04 December 2020
Accepted
12 January 2021
Published
15 April 2021
Volume
9 - 2021
Edited by
Bertrand Kibler, UMR6303 Laboratoire Interdisciplinary Carnot De Bourgogne (ICB), France
Reviewed by
Christophe Finot, Université De Bourgogne, France
Miguel Gonzalez-Herraez, University of Alcalá, Spain
Updates
Copyright
© 2021 Naveau, Vanderhaegen, Szriftgiser, Martinelli, Droques, Kudlinski, Conforti, Trillo, Akhmediev and Mussot.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Arnaud Mussot, arnaud.mussot@univ-lille.fr
This article was submitted to Optics and Photonics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.