Experimental realization of Fermi-Pasta-Ulam-Tsingou recurrence in a long-haul optical fiber transmission system

The integrable nonlinear Schrödinger equation (NLSE) is a fundamental model of nonlinear science which also has important consequences in engineering. The powerful framework of the periodic inverse scattering transform (IST) provides a description of the nonlinear phenomena modulational instability and Fermi-Pasta-Ulam-Tsingou (FPUT) recurrence in terms of exact solutions. It associates the complex nonlinear dynamics with invariant nonlinear spectral degrees of freedom that may be used to encode information. While optical fiber is an ideal testing ground of its predictions, maintaining integrability over sufficiently long distances to observe recurrence, as well as synthesizing and measuring the field in both amplitude and phase on the picosecond timescales of typical experiments is challenging. Here we report on the experimental realization of FPUT recurrence in terms of an exact space-time-periodic solution of the integrable NLSE in a testbed for optical communication experiments. The complex-valued initial condition is constructed by means of the finite-gap integration method, modulated onto the optical carrier driven by an arbitrary waveform generator and launched into a recirculating fiber loop with periodic amplification. The measurement with an intradyne coherent receiver after a predetermined number of revolutions provides a non-invasive full-field characterization of the space-time dynamics. The recurrent space-time evolution is in close agreement with theoretical predictions over a distance of 9000 km. Nonlinear spectral analysis reveals an invariant nonlinear spectrum. The space-time scale exceeds that of previous experiments on FPUT recurrence in fiber by three orders of magnitude.

The NLSE is an important exactly solvable model for the study of nonlinear phenomena. An example is modulational instability (MI) 1 , an exponential amplification of periodic random fluctuations at the expense of a pump wave that has been suggested as a possible mechanism for the generation of rogue waves 2 . The reversal of this process can give rise to repeated cycles of growth and decay, which constitute a realization of FPUT recurrence 3,4 . In the framework of the IST, these phenomena find a description in terms of exact solutions 5,6 associated with conserved nonlinear spectral degrees of freedom. From an engineering perspective, the prospect of encoding information in the invariant nonlinear spectrum is of high interest for optical communication systems, which today are limited by nonlinear interference 7 .
Various predictions of the underlying analytical NLSE theory have been observed in optical fiber experiments, including solitons 8 , Akhmediev breathers 9 and their collisions 10 , the Peregrine soliton 11 , and the Kuznetsov-Ma soliton 12 . Such experiments are not without challenges. They are typically conducted at average signal powers up to a few Watts 10,[12][13][14][15][16][17] . The dynamics take place over distances of several hundred meters up to few kilometers and on the picosecond scale. At such timescales, the generation of arbitrary initial conditions is difficult and has been approximated by amplitude modulation based on dual-frequency excitation 9,18 , by beating of two narrow-linewidth lasers to create a low-frequency modulation 11 , or excitation of superpositions of complex exponentials with tuned relative phases and amplitudes. The latter can be obtained from an optical frequency comb shaped with a programmable optical filter 10 . The observation of the spatial dynamics has been achieved with fiber cut-back experiments 12 . Simultaneous observation of amplitude and phase information can be realized by frequency-resolved optical gating (FROG) 11 [22][23][24] . They establish a deep mathematical link between the NLSE and the theory of compact Riemann surfaces and algebraic geometry 25 and had important influence on mathematical physics since their discovery more than 40 years ago 26 . Finite-gap solutions derive their name from the structure of their underlying nonlinear spectrum. They are defined through the IST in terms of a finite number of nonlinear spectral parameters. An important characteristic is that the main IST spectrum remains invariant under propagation of the solution governed by the NLSE -a property which bears the potential to eliminate nonlinear interference limiting data rates in modern fiber optic communication systems 7 .
Finite-gap solutions have provided analytical insight into the study of MI 6 . Here we construct solutions of the NLSE which exhibit periodic growth-decay cycles, as a hallmark of FPUT recurrence. In its dimensionless form, the focusing NLSE reads where ζ, τ are normalized space and time coordinates. In the IST, the field ψ ζ τ ( , ) plays the role of a scattering potential. For periodic potentials, the main IST spectrum is a discrete set of points in the complex plane for which the Bloch functions of the scattering problem are (anti-)periodic. A solution ψ ζ τ ( , ) is called finite-gap or finite-band if it can be described in terms of a finite number of spectral components. The main spectrum consists of N complex-conjugate pairs λ j , λ j , = … j N 1 , and remains invariant when ψ ζ τ ( , ) evolves according to the NLSE. The dynamics are encoded in the N auxiliary spectrum variables μ ζ τ ( , ) j . They evolve on the two-sheeted Riemann surface of genus = − g N 1 of an algebraic curve 5,23,25 The solution of the Jacobi inversion problem on the hyperelliptic Riemann surface leads to exact solutions in terms of g-dimensional Riemann θ-functions in the form  Figure 1a shows the branch cuts between pairs of eigenvalues λ j , λ j for = g 2 where the function λ λ = w P ( ) ( ) changes sign. To make w(λ) single-valued, con-www.nature.com/scientificreports www.nature.com/scientificreports/ sider two sheets of the same complex plane on which it takes a definite value. The projective plane { } is topologically equivalent to a sphere. If we cut the spheres and glue them together at their three branch cuts we obtain two spheres connected by three tubes, or equivalently, a torus with two holes: a genus-2 surface (Fig. 1b). There are 2g distinct closed curves, which can neither be deformed continuously into each other, nor shrunk to zero. We can imagine the g a-cycles going around the handles in Fig. 1b the short way, while the b-cycles go around the long way. They give rise to the g individual phases in the Riemann θ-functions. Most of the parameters in Eq. 2, in particular the g-dimensional vectors k and ω which determine the space and time period of the solution, are given by integrals over holomorphic differentials over the basis of cycles (see Methods and Fig. 1a).
Finite-gap solutions can hence be characterized by the genus of their underlying Riemann surface. Genus-0 solutions are plane waves, while genus-1 solutions are soliton-like and preserve the shape of the amplitude. Waveforms with the simplest nonlinear spectrum that exhibit nontrivial space-time dependence are genus-2. The Akhmediev breather, the Kuznetsov-Ma soliton and the Peregrine soliton are in fact limits of periodic genus-2 solutions which are known as solitons on finite background (SFB), and for which the period goes to infinity in space, time, or both 27 . Figure 2 shows two examples of finite-gap solutions. The genus-2 solution in a develops high peaks surrounded by troughs of minimal intensity (dark blue regions) starting from an almost flat initial condition and clearly shows out-of-phase recurrent behavior 16 . It has the symmetric spectrum − ± ± ± a bi ci a bi { , , } , which corresponds to an algebraic curve with two involutions (complex conjugation and negation) and implies the existence of double coverings of two genus-1 surfaces by the genus-2 surface. This property allows the linear combination of the cycles on the genus-2 surface in terms of cycles on two genus-1 surfaces with integer coefficients. As a consequence, wavevectors and frequencies become integer multiples of each other ( = k k 2 ) so that the solution is exactly periodic in time (ω = 0 0 ) and quasi-periodic (periodic up to a phase) in space 28 . The argument ωτ ζ + k in Eq. 2 remains invariant under a simultaneous shift of half a period in both space and time ζ ζ → + ζ p /2, τ τ → + τ p /2. Hence the waveform essentially repeats after a distance of half a spatial period, which defines the typical spatial scale of a growth-decay cycle. The genus-3 solution in Fig. 2b illustrates that finite-gap solutions can have rather intricate space-time dependence. The individual periods π ω =  www.nature.com/scientificreports www.nature.com/scientificreports/ quencies are not necessarily commensurate, so that not every finite-gap solution is periodic. In general, finding a spectrum for which the solution has a desired period is a difficult problem 25 . The nonlinear wave of the experiment. We focus on a periodic genus-2 solution for the experiment. The fiber link with periodic amplification can be approximated by an effective loss-less form of the integrable NLSE (1), provided the typical spatial scale of the solution is large compared to the span length L (see Methods). The equivalent integrable system has the same length as the experimental fiber link, but is governed by a renormalized nonlinearity parameter γ eff . From (2) we obtain a solution to the dimensionful NLSE through the substitutions is the effective nonlinearity parameter and Note that T 0 can be chosen freely to set the timescale of the solution. We haven chosen T 0 to obtain a time period = p 1 T ns to sample the waveform with reasonably high accuracy. In order to guarantee sufficiently slow evolution of the waveform we adjust the spectrum to give a spatial half-period that is a large approximate multiple of the span length  p L /2 z . By further choosing a spectrum with a large imaginary part and closely spaced points, we can obtain a significant peak power compared to the background 28 The resulting finite gap-solution has a recurrence length (half-period) of = p /2 5760 z km (~77 fiber spans) and average power of 2.1 mW. The space-time dependence of the power and phase of the waveform as obtained from Eq. 2 and scaled to the experimental units are shown in Fig. 3a,c, respectively. The point of maximal temporal compression occurs at 2880 km (~38 spans) with a peak power of 13.9 mW, roughly six times larger than that of the initial condition. A second peak reaches its maximum after 8640 km.

Measurements.
We performed measurements with the experimental setup shown in Fig. 4. A telecom-grade narrow linewidth distributed feedback fiber laser (1550 nm) is modulated by a pair of Mach-Zehnder modulators in an interferometric structure relatively phase-shifted by π/2 to transfer the complex electrical signal to the optical domain. The modulators are driven by an arbitrary waveform generator (AWG) which allows us to generate the initial condition in both amplitude and phase with a high sampling rate of 64 Gs/s per second, or 64 samples www.nature.com/scientificreports www.nature.com/scientificreports/ per period. An EDFA amplifies the signal to the desired theoretical average power at the entrance of the fiber. Losses due to fiber attenuation and components are compensated by the EDFA inside the loop. A third EDFA after the loop ensures a significant electrical amplitude at the oscilloscope. A balanced coherent receiver is used to separate real and imaginary parts of the waveform before it is sampled and digitized by a real-time oscilloscope (80 Gs/s). All further data processing, in particular downsampling to 64 Gs/s, timing synchronization, scaling and phase retrieval is done off-line in software (see Methods).
We inject multiple signal bursts consisting of a timing synchronization symbol and 5000 periods of the waveform into the fiber loop (first acousto-optical modulator (AOM) closed and second open to stop propagation of a previously recirculating signal). The number of bursts is determined by the roundtrip time (360 μs) of the signal inside the loop. When the loop is filled, the setting of the AOMs is reversed, blocking the signal coming from the modulator and letting the signal recirculate inside the loop. Half of the signal is coupled out of the loop via a 50:50 coupler at each revolution. An electric pulse generator synchronizes the operation of the AOMs and the oscilloscope to record the signal after a preset number of rotations in the loop. We vary the number of spans between 0 (back-to-back configuration) and 120, covering a total distance of 9000 km, corresponding roughly to the distance between Paris and Los Angeles. Amplitude and phase of the propagating wave are thus recorded with a spatial resolution of 75 km. Figure 3d is a false-color plot of the experimentally measured space-time evolution of the power. It has been obtained after timing synchronization, scaling and averaging over 25 samples of 5 consecutive periods for each distance (see Methods section for details). We can see good agreement with the theoretical waveform shown in Fig. 3a. In particular we see the buildup of the first and second peak and surrounding troughs of low intensity. Figure 3e shows a single realization of a single period. Depending on the realization of the amplified spontaneous emission (ASE) noise from repeated signal regeneration, the measurement can deviate from the theoretical prediction. As a result the second peak in Fig. 3d is smeared out in the averaging procedure.

Space-time evolution.
The exact solution of Fig. 3a has high symmetry (inversion symmetry ω ω ↔ − in frequency), which is apparently broken in the experimental data. In particular one can see a'tail' emerging from the first peak which crosses over into the second. The troughs can be seen to be located slightly asymmetrically around the peak. We emphasize that this qualitative deviation is not due to the fact that the system is only approximately integrable due to attenuation, because such deviation cannot break the symmetry. Deviations from the theoretical prediction occur where temporal compression and associated spectral broadening set in and are due to the tunable optical bandpass filter (OBPF). The Gaussian-shaped OBPF inside the loop has a measured full width at half maximum (FWHM) of 147 GHz. Assuming a misalignment of the filter with respect to the center frequency by 5 GHz, one obtains a small, but measurable asymmetry in the spectrum. The AOMs further cause a small drift of the signal in frequency in each pass through the fiber loop. Figure 3b shows a simulation based on the NLSE for the effective model, including a Gaussian-shaped filter of measured width, which is assumed to be offset by 5 GHz. The result qualitatively reproduces the features observed in the experiment.
In Figs. 5 and 6 we show slices of the experimental data at characteristic distances. Figure 5a shows the power of the waveform measured after passing through the setup of Fig. 4 without the fiber loop. It illustrates that we can precisely target the initial condition. 99% of the spectrum shown in b is confined in a bandwidth of 2 GHz. We compare the experimental data with the theoretical result provided by Eq. 2 scaled to the experimentally measured power and to split-step simulation of the NLSE (see Methods). All NLSE simulations include the effect of fiber attenuation and the optical bandpass filter inside the recirculating loop. Figure 5c,d show data at the point of maximal temporal compression in excellent agreement with the theoretical prediction. The spectrum has broadened by roughly a factor of 5 (99% of spectrum in 10 GHz bandwidth). The comb decays following a geometric progression with a maximum in the first sideband, showing that the energy is distributed over an infinite number of modes. It leads to the triangular shape that is characteristic of many parametrically driven systems 29 and which has been derived analytically for a particular exact periodic solution of the NLSE 30 . The triangular comb has been experimentally observed for the Akhmediev breather 9 and Kuznetsov-Ma soliton 12 . After this point the process is reversed, leading to FPUT recurrence.
Excellent agreement is also observed for the phase data shown in Fig. 6a,b, reconstructed from averaged measurements of the real and imaginary parts of the waveform. The phase difference between the peak and the background at the point of maximal temporal compression shown in b is very close to π, which can also be seen from www.nature.com/scientificreports www.nature.com/scientificreports/ Fig. 3c,f. In Fig. 6c we plot the difference between phases on two slices of fixed time, one traversing the peak ( = .
T 0 76 ns) and one traversing the background ( = . T 1 26 ns, see arrows in Fig. 3f) as a function of distance. Over a spatial half-period the relative phase difference between the peak and the background accumulates to roughly 1.2π. For the Peregrine soliton a shift of 2π has been observed 15 . Because of the periodicity of the solution we do not observe this limit.

Nonlinear spectral analysis.
We have shown that the waveform qualitatively evolves according to the theoretical prediction. In order to unambiguously demonstrate that the propagating wave is a finite-gap solution, we perform nonlinear spectral analysis of the real and imaginary components of the received data.
Nonlinear spectral analysis has previously been used to characterize SFB solutions and to identify similar structures in noise seeded MI 27,31 . Since SFB solutions correspond to degenerate cases of finite-gap solutions with an infinite period, the relevant peak structures have to be cut and periodized to obtain a non-degenerate spectrum, which, as a consequence, is not exactly constant even if the dynamics were ideally integrable. In our case we would expect to exactly recover the main spectrum of Fig. 7 if the dynamics are integrable (and without spectrum distortions). In NLSE simulations of the non-integrable system including attenuation we observe small oscillations of the main spectrum with distance, which however preserves its original shape and symmetry as expected.
We compute the nonlinear spectrum from the waveform via a fast nonlinear Fourier transform algorithm 32 implemented in the open-source FNFT package 33 . The algorithm is applied to one period of the raw data at a time. The resulting complex spectrum values are collected in a histogram. To obtain sufficient statistics we use five bursts of 5000 periods. The result is shown in the top row of Fig. 7 at a few characteristic distances (see arrows in Fig. 3d). At 975 km (15 spans) one can clearly distinguish three clouds of spectral points centered around the theoretical values as expected for a genus-2 finite-gap solution. The spread of the spectrum values is due to noise which primarily stems from the ASE of the optical amplifiers. Additional eigenvalues at small imaginary part arise because the periodic waveform affected by noise deviates from a pure genus-2 solution. Further similar artifacts occur along the real axis (not shown). As the distance increases, the noise distributions grow as ASE noise is added due to repeated amplification. While three eigenvalues can still be distinguished at 1950 km (26 spans), the distributions have merged after 3000 km (40 spans).
In the bottom row of Fig. 7 we show data obtained by applying the algorithm to averages over 25 periods to reduce the noise. We note that the procedures of averaging and nonlinear Fourier transform approximately commute, suggesting that noise is a small perturbation. The noise in the resulting spectra is considerably reduced. After 975 km (and in back-to-back configuration, not shown) the result almost perfectly agrees with the theoretical spectrum, which shows that departure from integrability is a small effect. Deviations of the spectrum values from the theoretical prediction start to occur around 1950 km, where spectral broadening sets in and the optical filter distorts the spectrum. Our NLSE simulations qualitatively reproduce the observed deviations from the ideal nonlinear spectrum as a result of the skewing of the Fourier spectrum through the optical filter. At 3000 km and beyond, the point clouds corresponding to two of the three spectral points have merged, but the distributions retain two discernible peaks of high intensity in each case. These results unambiguously confirm the propagation of a genus-2 finite-gap solution in our experiment. www.nature.com/scientificreports www.nature.com/scientificreports/

Discussion
We have observed the space-time evolution of a periodic finite-gap solution of the integrable NLSE in close agreement with theoretical predictions. The experiment has been conducted in a testbed for optical transmission experiments built from commercially available standard components. Modern coherent telecommunication equipment is designed for modulation and detection of both amplitude and phase with high sampling rates. Hence we were able to precisely target the initial condition and measure the real and imaginary parts of the complex envelope of the electrical field with high accuracy. Together with the recirculation loop, the setup permits a non-destructive full-field characterization of the spacetime dependence of the propagating wave. This allowed us to experimentally verify a key prediction of the finite-gap integration theory: through nonlinear spectral analysis we showed that the nonlinear main spectrum is preserved over transcontinental distances. The prospect of encoding information in the invariant nonlinear spectrum is of high interest for applications in optical communication.
We further found that despite the presence of ASE noise from repeated amplification, the evolution of the optical field is accurately described by the integrable NLSE. Even though ASE noise becomes a limiting factor at large distances, the remarkable agreement between experiment and theory shows that such a setup provides highly controlled conditions for observing nonlinear wave dynamics.  Fig. 3f for fixed times, one passing through the peak and one between two peaks, showing the accumulated phase of the peak relative to its background.

Figure 7.
Histograms of nonlinear main spectra at characteristic distances. The nonlinear main spectra computed based on the experimental raw data for the waveform (top) and using the averaged data (bottom). Noise accumulates with distance due to repeated amplification, leading to distributions around the theoretical prediction (red dots) and artifacts with small imaginary part along the real axis. In the spectra obtained from the averaged data, the three distinct spectral points of the genus-2 finite-gap solution can be distinguished for all distances. Note that while for 3000 and 4050 km the point clouds have merged, they contain two discernible peaks of high intensity. Departure from the theoretical prediction is due to the optical filter, see text. www.nature.com/scientificreports www.nature.com/scientificreports/ The observed periodic compression dynamics constitutes a realization of FPUT recurrence, which we have observed for close to two growth decay cycles. The results are in line with recent observations of the space-time dependence of FPUT recurrence 16,17,34,35 . In particular, our results bear remarkable qualitative resemblance to the very recent ones of ref. 17 , while exceeding them in terms of scale in space and time by three orders of magnitude. At these scales the NLSE is an excellent model of the nonlinear dynamics. In their experiment, loss is compensated by a Raman pump. A small number of waves (pump, signal, idler and harmonics) is launched into the fiber and the time-domain signal is reconstructed by means of the Fourier transform of up to five waves. The propagating and reconstructed waves therefore correspond to time-periodic fields similar to the one in our experiment.
Here we even achieve quantitative agreement with the underlying theory.
We have observed FPUT recurrence over more than one cycle. Previously, approximate recurrence has been observed over three complete cycles, albeit in a photorefractive KLTN crystal 35 . In order to observe multiple recurrence cycles with our setup, the number of spans can be increased, or the spatial period of the waveform may be reduced for fixed fiber length. Both lead to deviations from integrability: due to increasing noise for the first option and a departure from the validity of the effective model for the second. For the latter, the effect can be mitigated in two ways. By either decreasing the span length, again decreasing the signal to noise-ratio, or by increasing the separation between nonlinear spectrum points, which reduces the dynamic range of the signal power. Distributed Raman amplification has the potential to provide a flat gain to further increase the achievable distances.
While our experiment has been conducted neglecting polarization effects, the two polarization components of single-mode fibers and the various modes in multi-mode fibers are described by the integrable Manakov system 36,37 . This opens possibilities to design analogous experiments on generalizations of finite-gap solutions to the Manakov system 38 or for the study of vector solitons.

Methods
Inverse scattering transform and finite-gap integration. The IST is a well-established method to solve certain classes of integrable equations. It was first discovered for the Korteweg-de Vries equation 39 and later by Zhakharov and Shabat 40 for the NLSE. In the IST, the NLSE is represented as the compatibility condition ∂ ∂ Φ = ∂ ∂ Φ τ ζ ζ τ of two linear differential equations for the scattering data ζ τ λ Φ( , ; ), which holds true if and only if ψ ζ τ ( , ) satisfies the NLSE. For given initial condition ψ ζ τ = ( 0, ), obtaining the scattering solution ζ τ λ Φ = ( 0, ; ) is known as the Zakharov-Shabat scattering problem. The spatial evolution of the scattering data is governed by  and allows one to determine the solution of the NLSE at any ζ by solving the inverse scattering problem, that is, by recovering the potential from the scattering data. The main spectrum of a periodic finite-gap solution is defined by the 2N discrete values λ j which occur in complex conjugate pairs and for which the scattering solution, i.e. the Bloch functions, are (anti-)periodic in τ: ζ τ λ ζ τ λ Φ =±Φ + T ( , ; ) ( , ; ). Solutions of the NLSE can be constructed by solving a Riemann-Hilbert problem 41 , or the finite-gap integration method, an algebro-geometric approach 6,25 . Here we adopt the latter method to construct a periodic solution (2) by computing its parameters in terms of integrals over the hyperelliptic Riemann surface defined by the main spectrum. To this end, we define a canonical homology basis , } i j of a-and b-cycles on the two-sheeted Riemann surface such that each a-cycle crosses one b-cycle exactly once and otherwise crosses no other a-or b-cycle, and the same holds for the b-cycles. This is expressed in terms of the intersection numbers as The procedure can be automatized and for example yields the homology basis shown in Fig. 1a. In terms of the basis of holomorphic differentials , the invertible g × g matrices of a-and b-periods are defined in terms of the integrals We evaluate these integrals numerically. From these we obtain the frequencies ω , and the period matrix τ = − A B 1 . In terms of the normalized differ- They involve the initial condition for the auxiliary variables. However their difference Not every initial condition of the μ j corresponds to a finite-gap solution. Since we are not interested in the precise evolution of the μ j , we exploit the fact that any choice of δ − yields a valid solution as long as its imaginary part is zero 6 . Hence we can set δ ≡ − 0 j and compute δ + j using the formula for their difference. For a finite-gap solution half of the 2N real degrees of freedom in the μ j are constrained so that this choice completely determines the auxiliary spectrum. The remaining parameters K, k 0 and ω 0 are determined 25  www.nature.com/scientificreports www.nature.com/scientificreports/ crosses a branch cut once to change sheets, while avoiding the holomorphic basis of differentials. see Fig. 1a. The Riemann θ-functions finally are evaluated numerically for given argument by truncating the series according to a precision criterion.
Experiment and data preprocessing. To map the spatiotemporal evolution of the waveform we perform independent measurements for each distance (number of turns in the recirculating loop). The AWG periodically generates sequences of 5000 periods of periodic waveform (5 μs) separated by a 42 ns guard interval and preceded by a Schmidl-Cox synchronization symbol 42 (42 ns) separated from the waveform by a 21 ns guard interval. We employ a telecommunication-grade distributed feedback fiber laser with linewidth ~100 Hz operating at 1550 nm. The optical IQ modulator (Fujitsu 32GB DP-QPSK) is driven by an arbitrary waveform generator with a sampling rate of 64 Gs/s, bandwidth 25 GHz, and effective number of bits (ENOB) of 5.5. Out-of-band spectral components from the modulator are filtered using a rectangular optical bandpass filter (OBPF, Yenista XTA-50). To remove any DC components of the envelope (the low-frequency limit of electrical drivers of the MZ is 35 kHz), we digitally shift the spectrum by f = 5 GHz. Insertion losses of the acousto-optic modulators (AOM) and couplers are compensated by an Erbium-doped Fiber amplifier (EDFA). At system setup we open the recirculating loop and adjust the power setting of the first EDFA to obtain the desired theoretical average power at the entrance of the fiber. With closed loop and circulating signal, the EDFA inside the loop is adjusted based on the output of the low-bandwidth photodetector (PD) signal to compensate fiber attenuation and losses due to components. The recirculating fiber loop itself consists of 75 km of standard single-mode fiber (G.652) and a Gaussian-shaped OBPF (measured FWHM 147 GHz) to suppress out-of-band noise. We minimize power fluctuations in the loop due to polarization dependent loss through a polarization controller (PC). The power is monitored by a low-bandwidth photodetector. An electronic pulse generator is used to control the number of roundtrips inside the loop through the AOMs and to synchronize the oscilloscope. A second PC gathers the signal in a single polarization component. The output of the intradyne phase diversity coherent receiver is proportional to the amplitude of the in-phase and quadrature components of the signal. They are independently digitized using an 80Gs/s oscilloscope (Keysight DSOZ504A, ENOB 6 per sample). Each pass through an AOM downshifts the center frequency by 40 MHz. On the receiver side, the 5 GHz frequency shift and any residual carrier phase offset caused by the AOMs are compensated by shifting the center frequency of the spectrum to 0 Ghz (in digital domain). The signal exhibits a peak at multiples of 1 GHz due to its periodicity, which can be used as a reference. The signal in digital domain is further downsampled to 64Gs/s. An approximate relative timing to locate the relevant data is obtained through the Schmidl-Cox synchronization algorithm. Due to noise in the transmission, the timing can be off by up to 6 samples 42 (64 samples per period). The data are processed independently for each distance. To avoid edge effects, 100 periods on each end are discarded. The remaining 4800 periods are subdivided into groups of 25 × 5 periods. For each group of the 40 groups we average over the 25 realizations per group to obtain averaged waveforms of 5 periods length. Since our setup only detects relative phase, absolute phase for each distance is recovered by aligning the phase at the maximum amplitude at the respective distance to the corresponding theoretical value. Sample-accurate timing synchronization is finally achieved by aligning the phase of the averaged data (within a window of one period) to the phase of the theoretical waveform at a fixed time (T = 4.875 ns). The same timing correction is applied to real and imaginary parts of the averaged and unaveraged data. . Both are significantly larger than the span length, but are shorter than the fiber link. The scales are of the same order of magnitude and cannot be separated. Dispersion and nonlinearity act together in governing the evolution of the waveform.
We perform simulations of the NLSE (5) using a standard split-step integration scheme 1 with adaptive step size controlled by a maximum phase change of 0.05°. Since the implementation is based on FFT and implicitly assumes periodic boundary conditions, it is sufficient to simulate a single period. To quantify the departure from integrability, we performed NLSE simulations of the fiber link including attenuation. We simulate the system consisting of a fixed number of identical fiber spans including the Gaussian-shaped filter with measured FWHM and with lumped amplification. We neglect ASE noise to emphasize small deviations from integrability arising from attenuation and the filter. We further assume a chromatic dispersion coefficient β = − . − 21 5 ps km 2 2 1 , nonlinearity parameter γ = . ⋅ − − 1 3 W km 1 1 and attenuation α = 0.2 dB/km (L eff = 21 km). The simulations show a relative deviation from the theoretical result for the integrable system of less than 1% on average, with a maximum difference of 9% reached at the smallest amplitude. Idealized NLSE simulations of the effective model neglecting attenuation yield essentially the same results as the theoretical prediction by Eq. (2), confirming the high accuracy of the split-step Fourier integration method.