Using the integrable nonlinear Schrödinger equation (NLSE) as a channel model, we describe the application of nonlinear spectral management for effective mitigation of all nonlinear distortions induced by the fiber Kerr effect. Our approach is a modification and substantial development of the so-called “eigenvalue communication” idea first presented in A. Hasegawa, T. Nyu, J. Lightwave Technol. 11, 395 (1993). The key feature of the nonlinear Fourier transform (inverse scattering transform) method is that for the NLSE, any input signal can be decomposed into the so-called scattering data (nonlinear spectrum), which evolve in a trivial manner, similar to the evolution of Fourier components in linear equations. We consider here a practically important weakly nonlinear transmission regime and propose a general method of the effective encoding/modulation of the nonlinear spectrum: The machinery of our approach is based on the recursive Fourier-type integration of the input profile and, thus, can be considered for electronic or all-optical implementations. We also present a novel concept of nonlinear spectral pre-compensation, or in other terms, an effective nonlinear spectral pre-equalization. The proposed general technique is then illustrated through particular analytical results available for the transmission of a segment of the orthogonal frequency division multiplexing (OFDM) formatted pattern, and through WDM input based on Gaussian pulses. Finally, the robustness of the method against the amplifier spontaneous emission is demonstrated, and the general numerical complexity of the nonlinear spectrum usage is discussed.
© 2013 Optical Society of America
Signal transmission down the fiber-optic line is affected simultaneously by several physical phenomena in fibers: dispersion, signal corruption due to the amplifier spontaneous emission, fiber nonlinearity, and other effects [1,2]. The fiber Kerr nonlinearity is the important feature distinguishing the fiber-channel transmission from the radio-type (open air or space) channel : The dependence of the transmitting media properties on the signal intensity results in a number of detrimental effects produced by the coupling and nonlinear intertwining between the signal degrees of freedom, accompanied by the influence of noise and dispersion. Thus, the data encoded into the signal propagating through the nonlinear channel become corrupt, leading to loss of the transmitted information and decrease in the spectral efficiency of transmission [3–15]. For the suppression of nonlinearity-induced distortions, a number of techniques have been proposed and studied, e.g. digital backpropagation , receiver-based digital signal processing [17–20], pilot-based nonlinearity compensator , optical  and twin-wave  phase conjugation, fractional Fourier-based pre-compensation , to mention several recent advancements, and many other types of pre-compensation, in-line and post-compensation methods [19,23]. In this work we describe the application of a method (developed a long time ago in the theory of nonlinear evolutionary equations [1, 25–29]) that allows effective linearization of some classes of nonlinear communication channels, to the problem of optical communications.
An important subclass of nonlinear transmission channels is governed by the so-called integrable evolutionary equations [1, 25–29]. One of the particular manifestations of this integrability property is that given the initial conditions (in the considered context, the waveform of the input signal), we can propagate the signal to a distance L in three steps, which have direct analogies with the same stages in the consideration of linear problems (although the implementation of these steps differs significantly from the linear case). (i) The first step is the mapping of the input profile to the spectral domain. For a linear channel, this operation corresponds to the ordinary forward Fourier transform (FT), and according to the already established terminology, we will call this stage the forward nonlinear Fourier transform (NFT), by analogue with the linear situation. For nonlinear propagation, this stage involves solving the specific direct scattering problem associated with an integrable equation and produces a set of scattering data, where the particular quantities (continuous spectrum and a set of localized complex eigenvalues, if the latter exists) are then associated with orthogonal nonlinear “normal modes”. (ii) The next step is the propagation of the initial spectral distribution (again, a continuous spectrum and complex eigenvalues) to distance L: Here, in both linear and nonlinear cases, the spectrum evolves according to the linear dispersion law, and this operation is usually quite trivial, amounting to the phase rotation of the corresponding decoupled spectral component. We would like to stress that even if the original problem is nonlinear, the nonlinear spectral components evolve without interaction with each other in full analogy with the linear problem, which explains the term “linearization”. (iii) The last stage is the recovery of the solution profile in the space-time domain: It is the backward FT in the linear case. For the nonlinear integrable problem, the backward NFT amounts to the solution of the so-called Gelfand-Levitan-Marchenko equations [1, 25–29], and this step accomplishes the finding of a solution (signal profile) at distance L. These three stages which provide the solution of a nonlinear equation at distance L, constitute the method of the inverse scattering transform (IST) [1, 25–29]. Another important aspect of the NFT [1, 25–30] is that, due to effective linearization of the channel, all nonlinear propagation effects such as self-phase modulation, cross phase modulation, and four-wave mixing are effectively removed owing to the orthogonality property of the nonlinear spectrum (NS) components (involving both discrete (soliton) and continuous spectrum) , thus removing the very reason for the spectral efficiency decay at high signal power. In engineering terms, the NFT method can be treated as effective nonlinear orthogonal frequency division multiplexing . In practical terms, the usage of NS is a novel approach to combat the nonlinear distortions, being an alternative/addition to the commonly used engineering techniques (see e.g. [16, 17, 19, 20, 22–24] and references therein).
The evolution of the slow-varying optical field envelope in a passive optical single-mode fiber is described by the nonlinear Schrödinger equation (NLSE) [1, 2, 32, 33], which introduces the effects of the mutual counteraction of dispersion and nonlinearity. In the absence of losses and external actions, the NLSE is integrable, and so the explicit form for the three aforementioned operations (forward/backward NFT and the spectrum evolution law) is known . Since we can decompose the input signal into analogues of orthogonal Fourier modes by performing the forward NFT operation, the dynamics of the signal through the NLSE channel can be completely described by the evolution of these orthogonal modes [1, 25–29], and then, the informational content encoded at the transmitter can be potentially recovered at the receiver without distortions.
The quest for the optimal design of a nonlinear transmission channel and usage of the nonlinearity in a ‘constructive’ way has occupied researchers for a long time. This includes soliton transmission techniques, dispersion management (e.g. dispersion-managed soliton) for mitigation of nonlinear impairments, optical regeneration, and other approaches. The idea of manipulating the robust NS of NLSE for the purposes of optical transmission was already put forward nearly 20 years ago  and has been intensively scrutinized recently in . Notably, in the original study , only the soliton part of the nonlinear spectrum was considered (and these ideas were advanced further in ) and no practical fiber optical implementation scheme was suggested. The solitons are described by the localized (complex) eigenvalues of the corresponding spectral problem defining NS. These eigenvalues are the conserved quantities: Their particular magnitudes emerge directly from the parameters of the input profile (encoded signal) and do not change during the propagation. Therefore, these eigenvalues (or even the total number of discrete eigenvalues) can be utilized for encoding and transferring the information, owing to their conservation. However, the usage of pure solitons for transmission purposes has a number of challenges, e.g., the soliton collision problems in WDM systems, low symbol rates etc. [2, 32], although this direction has been revived recently [35, 36] through plugging in the multilevel modulation formats and using coherent transmission with the subsequent effective increase in the transmission rate.
We would like to point out that the direct implementation of the transmission based on the NS brings about significant technical challenges. First, there is no fully developed electronic or optical scheme that could perform forward/backward NFT without the computer processing of the signal. Second, the properties of the NS can be significantly different from its linear counterpart, and therefore, if one deals with the “conventional” linear formats, no general method for encoding and recovery of the information exists: Each particular configuration of the NS emerging from the input signal should be studied from scratch, as the NS is a nonlinear functional of the input. Even the number of analytically solvable cases for the forward NFT is very limited . Thus, the principal possibility of NLSE channel linearization should not create any wrong impression, digressing one from the hurdles and challenges in the practical implementation of this technique.
The major goal of the current paper is to outline some solutions for the aforementioned problems. First, in our study, we withhold the use of solitons (so that we do not deal with the discrete part of the NS), thus disposing of the problems associated with the soliton WDM collisions, etc. In this paper, we propose to utilize the continuous part of the NS for information encoding and transmission. As mentioned earlier, the advantage of using NS in the integrable NLSE channel is that the nonlinear modes of the continuous NS part are still orthogonal and evolve according to the linear dispersion law: The NS coding is free from the “coupling” nonlinearity influence allowing virtually undistorted transmission of initial data provided that the linear dispersion is compensated at the endpoint.
Further, we describe a general method of calculating the continuous NS when the average power of the input signal (localized within the boundaries of a considered time interval) is not very high and nonlinearity can be treated as a perturbation. We show that in this case, the forward NFT can be performed iteratively and that the NS can be computed with a desirable precision in terms of expansion with respect to the input signal power by implementing the Fourier-type recursive integrations, involving only the input signal itself. Note that these operations can already be suggested for the electronic and/or all-optical implementation.
Finally, based on the expansions of the NS, we are able to cancel the “nonlinear” additions of the NS and, thus, can equalize the nonlinear and linear spectra by means of nonlinear pre-distortion of the input signal. The equalization of the linear and nonlinear spectra after the nonlinear pre-distortion paves the way for the utilization of “linear” methods, such as ordinary FT for information encoding and decoding. This approach is then studied in detail for illustrative purposes in application to a fragment of an OFDM input pattern with the phase-shift keying (PSK) encoding of OFDM coefficients. It is shown that pre-distortion allows one to obtain an almost undistorted initial signal profile at the receiver after winding out the accumulated linear dispersion in the NS. In the case of the OFDM input, the pre-distortion amounts to the extraction of the lowest-order terms associated with the inter-carrier (IC) and inter-slot (IS) interference effects.
Our paper is organized as follows. First, in Section 2, we introduce the NLSE and the associated normalizations, determining the scales of the entering parameters, and provide the details of the forward NFT, explaining the relevant quantities characterizing the NS (i.e., we present the elements of the direct scattering problem associated with the NLSE). Section 3 deals with the general theory of NS when the amplitude of the input signal is sufficiently low: We present the expression of the NS in terms of the signal amplitude series. Subsection 3.2 explains how these results can be geared for the effective equalization of the linear and nonlinear spectra, used further for our nonlinear pre-distortion method. In Section 4, we consider the NFT for the OFDM-type inputs. First, in Subsection 4.1, we remind the basics of the OFDM format. After that, in Subsection 4.2, we derive the analytic form of the NS for a single OFDM tone with a finite base extent, and in Subsection 4.3, we mention some general properties of NS associated with the OFDM-type input. Section 5 deals with the direct application of our spectra equalization method. First, we consider the pre-distortion of a single OFDM tone (Subsection 5.1). Then, Subsection 5.2 contains the results for the pre-compensation of a sequence of OFDM-coded slots. Section 6 contains the extension of our approach and discusses how to deal with the situation where the recursive spectra equalization is not effective. First, in Subsection 6.1 we consider precompensation of the input in the form of Gaussian-based WDM pulse and show, that our method is applicable in this case as well. In Subsection 6.2 we analyze the robustness of our method against the amplifier spontaneous emission (ASE). Finally, in Subsection 6.3 we present the estimation of the numerical complexity of the methods based on NFT as compared to digital back-propagation and outline the way of how to further extend the ideas of NFT-based transmission. Our findings are summarized in the Conclusion.
Throughout the paper, we adopt the following conventions for the linear forward/backward FT:
2. The model, normalizations, and NFT basics
2.1. NLSE and normalizationsEq. (2)) and hence deal with the so-called focusing type of NLSE. We take the typical value of β2 = −22 ps2/km, and the higher-order dispersion terms are not considered here. The instantaneous Kerr nonlinearity coefficient γ is expressed through the nonlinear part of refractive index n2 and an effective mode area Aeff : γ = n2ω0/cAeff, with c being the vacuum speed of light and ω0 = 2πν0 being the carrier frequency of the envelope q(t, z); the typical value is γ = 1.27 (W·km)−1.
We normalize time in Eq. (2) to the characteristic time related to an input signal Ts, which can be, e.g., the extent of the RZ signal or the duration of a single OFDM slot, and then use the effective z-scale associated with . Then, we measure the power of the input in units of γZs and normalize the signal amplitude. The summary of normalizations is as follows
Further, we typically assume, without loss of generality, the following values: Ts ≈ 300ps, giving the characteristic z normalization scale as Zs ≈ 4000km. For instance, for the typical input peak power of P0 ∼ 0.1mW, the characteristic measure of the effective nonlinearity is ε = (γZsP0)1/2 ≈ 0.5, and for P0 = 1mW, ε ≈ 2.5
The normalized NLSE (2) is rewritten as
2.2. Forward NFT as a part of the IST method
The forward NFT operation for the NLSE (4) requires solutions of the so-called Zakharov-Shabat spectral problem (ZSSP) [1, 25–29], which corresponds to the scattering problem for a non-Hermitian (in the case of anomalous dispersion) Dirac-type system of equations for two auxiliary functions ϕ1,2(t), with the NLSE input pulse profile q(0, t) ≡ q(t) serving as an effective potential entering the equations1, 25, 29]).
At the left end t → −∞ we fix the “initial” condition for the incident wave scattered by the potential q(t) to have the so-called Jost solution Φ⃗(t, ζ) = [ϕ1(t, ζ), ϕ2(t, ζ)]T:Eq. (5) is then defined as Eq. (5) (if these are present). The quantity r(ξ) from Eqs. (7), (8) defined for real ξ (“frequency”) plays the role of the continuous nonlinear spectral distribution, while the quantities associated with discrete eigenvalues describe solitonic degrees of freedom and do not have analogues in the linear problems. In our case, we deal with the situation where no solitons are present, and therefore, the ZSSP (5) does not contain a discrete eigenspectrum:
The backward NFT maps the scattering data Σ onto the space-time domain and requires the solution of the Gelfand-Levitan-Marchenko (GLM) system of equations [1, 25, 29] associated with a particular Σ. However, for the method explained in the beginning of the paper we need only the forward NFT, so that we postpone the discussion of backward NFT to Section 6.3.Eq. (7), i.e. via solving the ZSSP (5) where the input profile serves as a potential. From Eq. (10), one can see that the NS obeys the linear dispersion law of NLSE (4) if one associates the linear frequency ω with the quantity ξ as ξ = −ω/2. Indeed, from the IST theory, it is known  that asymptotically in the linear limit |q(t)| → 0, the following formula is valid: Eq. (1). In view of Eq. (11), it is useful to define the spectral function N(ω) associated with forward NFT via the reflection coefficient: (12) coincides with its linear counterpart. Further, the term NS implies the spectral function N(ω) defined in Eq. (12), rather than the reflection coefficient r(ξ) itself.
3. Expansions for the low signal amplitude ε ≪ 1 and the concept of spectra equalization
In this section we first, considering the amplitude of the input signal |q(t)| ∼ ε to be small, derive the expansion for the reflection coefficient r(ξ) in terms of ε. Note that the results given in Subsection 3.1 below are general and do not imply any specific modulation format. In fact, the only limitation of the method is the relative smallness of the input signal power, given by ε. The smallness of ε allows one to obtain the desired expression using perturbation theory, similar to that used in  in the context of the fiber Bragg gratings design. After that, for the sake of illustration and outlining how the method can be applied to specific formats, we successively consider the application of the results obtained to a single tone of OFDM, a single slot of OFDM, and several OFDM slots serving as the NLSE input.
3.1. NFT in the case of a weak nonlinearity (small enough input signal power)
First, one makes the transformation to the slow varying functions φ1,2 as: φ1,2 = e∓iξtϕ1,2. In terms of φ1,2, the expression for r(ξ) (7) now changes toEqs. (5), we arrive at the system (14) as Eq. (15) into Eq. (16) and vise versa recursively to obtain the following expansions:
The expression for r(ξ) takes the form r ≈ r0(ξ)+r1(ξ), where r0 ∼ ε and r1 ∼ ε3 are given as follows (the next term in this expansion is ∼ ε5):Eq. (10) (2) itself rather than spectral domain) may also be used for calculation of the NS corrections.
We note, that of course one could also in principle achieve the equalization result without resorting to IST by linearizing NLSE (4) directly in the time domain. Yet when considering the effects of other perturbations on signal propagation (like e.g. ASE noise) the NS-based approach is the best way to go forward. This is due to the very mathematical structure of integrable NLS : The nonlinear spectrum plays the role of so called adiabatic invariants and as such represents a valid basis for any perturbative expansion (as it does not develop secular terms). The NS based approach also provides the correct anzatz for the large distance asymptotes of the NLSE solutions .
3.2. Linear and nonlinear spectra equalization using signal pre-distortion
Suppose that at the input z = 0, we apply pre-distortion s(t) ∼ ε3 to the initial signal waveform q(t):Eq. (20) and, thus, the term R1(L, ξ) from the spectral density at the end point z = L in Eq. (21), by using the additional pre-processing given by s(t). Looking at the form of r(ξ) expansion, we see that when a small quantity s(t) ∼ ε3 is added to the input signal, then we gain a correction rs(ξ) ∼ ε3 to the expression for r1(ξ), see Eq. (20): (12), we can now obtain the ordinary Fourier spectrum S(ω) for our correction s(t) as Eq. (25), we restore the profile of s(t) in the time domain. Thus, for the predistorted signal qs(t) = q(t) + s(t) with the FT of s(t) given by Eq. (25), the addition to the nonlinear spectrum ∼ ε3 disappears altogether and the NS associated with qs(t) coincides with the linear spectrum of initial q(t) up to the terms ∼ ε5.
The flowchart of the pre-distortion scheme and the signal recovery from the pre-processed NS at distance z = L is given in Fig. 1.
Note that by considering the higher-order nonlinear pre-distortions to q(t), we can cancel the higher terms in the expansion of NS. At the same time, we emphasize that the recursive Fourier-type integration in Eq. (20), used to obtain r(ξ), can be potentially implemented through the electronic or all-optical realization in order to calculate the NS associated with a particular given signal. One of the novel and significant points of our approach is that we are able, by means of precompensation described above, to translate the information encoded in the traditional way in time domain to the nonlinear spectral domain without using any sophisticated special new formats, and control the accuracy of the data mapping. Finally, the transmission itself is effectively performed through the NS evolution.
4. An illustration of the NFT application: OFDM-modulated input signal
4.1. OFDM basics
Coherent optical OFDM (CO-OFDM) has recently become a popular transmission technique owing to its robustness against chromatic dispersion, polarization mode dispersion, and practicality of implementation [19–22,39,40]. OFDM is a multi-carrier transmission format where a data stream is carried with many lower-rate sub-carriers (tones):
The linear spectrum of the OFDM-modulated signal (26) is a comb of sinc-like shapes:
In this work, for illustration purpose, we employ the widely used QPSK encoding of the coefficients cαk [19, 36, 40–42] (without losing the generality of the analytical results), with four allowed phase levels used for encoding the informational content: This means that the absolute value of cαk is the same for each coefficient, |cαk| = c = const., and the phase of each cαk takes four discrete values from the set:
4.2. NFT in the case of a single OFDM tone
Let us start with the idealized but analytically solvable case when we have only one OFDM slot of length T (extending, say, from −T/2 to T/2) with a single encoded coefficient, see Eq. (26), as the input for the NLSE channel. This means that, at z = 0, we have the following profile q(0, t) = q(t)(5). The solution for Eqs. (5) for the potential given by Eq. (29) can obtained analytically using the results for the rectangular pulse [43, 44]: Eq. (7): Eq. (31) is
4.3. Correspondence between the NS and linear spectrum, soliton creation threshold, and general remarks on the NS associated with the OFDM input
In this particular case, one can check directly the relations (11), (12) by seeking the limiting expression of a single-carrier NS when |ck| → 0. Indeed, in the limit of a small ck from Eq. (12) using (31) or (32), one immediately gets Eq. (33).
In , it was demonstrated that the singularities of r(ξ) arising on the real axis ξ due to the increase of the amplitude/length of the signal correspond to the appearance of bound (discrete) states in the ZSSP (5), i.e., bright optical solitons. Such a creation process has a threshold character with respect to both pulse width and amplitude. In the case of a single-carrier OFDM pulse, the amplitude threshold for |ck| needed for creating the first soliton can be easily obtained from Eq. (32): The expression for r(ξ), Eq. (31), has a singularity when the denominator vanishes, which happens if ψk = 0 and ξkT = π/2. These two conditions give29]) for a soliton nucleation, although the latter is strictly valid only for single-hump pulses with a constant phase  and is usually violated in the case of chirped inputs [45, 47], where the creation of solitons is significantly suppressed (the latter fact was utilized in the concept of a soliton-based discriminator ). In the context of coherent optical communication, due to highly mixed phases of the optical signal resulting from signal overlap from different time slots or from the use of many sub-carriers in OFDM, the phase of an optical signal becomes effectively random  and the soliton formation is dramatically suppressed . Note that such a phase mixture together with the signal considered as noise (for technical purposes) can be enhanced by initial pre-dispersion technique . Finally, the value of ξth gives the position of appearance of the first bound state on the real axis and is related to the velocity of the emerging soliton: vsol = Ωk.
Now, let us compare the nonlinear and linear spectra of an OFDM tone, i.e., Eq. (31) and the corresponding NS, Eqs. (12) and (33). The comparison is made in Fig. 2, where we plotted an NS (red) and a linear spectrum (black) for different values of the input amplitude. In the linear case, the main spectral peak corresponds to the single OFDM frequency Ωk, and the sinc-shape of the spectrum is the consequence of the rectangular pulse-shaping function in the time domain. The NFT spectrum (32) after transformation (12) maintains a similar shape—it is also symmetric with respect to the main peak Ωk regardless of the value of ck. The only property that changes drastically is the scaling of the peak with the amplitude. In the linear case it scales as |ck|2, but in the nonlinear case, it actually scales as tan2(|ck| T), becoming infinite at the threshold |ck|th. So, for the subthreshold ck, the NS deviates only insignificantly from the linear spectrum, see Fig. 2(a). Then, when we approach |ck|th, we observe an appearance of a pronounced peak indicating the tendency for soliton nucleation, see Fig. 2(b). For the values of ck beyond the threshold, we already observe the significant structural changes of NS as compared to the linear spectrum, in Fig. 2(c).
For the general form of the OFDM input (26) inserted as a potential in the ZSSP (5), it is impossible to find the analytical expression for the reflection coefficient r(ξ) and, therefore, for the NS (12), as well as for the soliton creation thresholds in terms of the signal extent and amplitude. Moreover, when random values of coefficients entering Eq. (26) are used, each particular realization produces its own NS profile and threshold parameters, and hence, one has to resort to numerical methods  for finding the form of NS [see Fig. 3(a)]. However, the conclusions reached above remain qualitatively valid: In the case of low power of the pulse (i.e. below the soliton creation threshold), the NS should have the structure similar to the linear spectrum of the signal (see Fig. 3) with only quantitative deviations, the measure of which is stipulated by how close to the threshold values we are. This fact is used in the next section for performing effective equalization of the NS and linear spectrum by extracting these small deviations between them.
Note that for the chirped and random profiles, the soliton nucleation is strongly suppressed by both the chirp intensity (rapid phase variations)  and the incoherence in the input signal [44, 45, 48–50], and one can expect the mitigation of the threshold values (34) produced by a single tone in the case of a multitone randomly coded signal.
In Fig. 3 we also present the results for the dynamics of the linear and nonlinear spectra associated with a fragment of the OFDM input (26): We took 3 slots of the 10-subcarrier OFDM with the pseudo-random QPSK encoding of coefficients for each subcarrier. The extent of each OFDM slot is T = 1, and the amplitude of each coefficient is |cαk| = 0.25. The case considered produced the subthreshold pattern, i.e., no solitons were nucleated from the input used in our plot. The propagation distance was z = 1.2 in normalized units, which corresponds to L ≈ 5000km for the normalization values described at the end of Subsection 2.1. First, we calculated the NS N(ω), red line in Figs. 3(a) and 3(b), propagated the input pattern to the distance z = 1.2 and again calculated N(L, ω), blue line in Figs. 3(a) and 3(b) (for the latter case, we removed the linear dispersion). One can see that the NS does not change, which is in accordance with the NS dynamics formula (10): By winding out the linear dispersion corresponding to the distance L, one obtains the exactly initial spectral distribution N(0, ω) up to the numerical errors associated with the NLSE integration and calculation of NS. In contrast, the profile of Q(L, ω), green line in Figs. 3(c), 3(d), visibly differs from Q(0, ω), black line in Figs. 3(c), 3(d), even when we remove the dispersion, because of the nonlinear distortions of the signal during the propagation. Some results for the NS attributed to several other types of input profiles are given in .
5. Implementation of the nonlinear pre-distortion
Let us now consider how the general pre-distortion scheme considered above works in a particular case of the OFDM-modulated pulses.
5.1. Pre-distortion of a single OFDM subcarrier in a single time slot
First, we consider the pre-distortion for the input in the form of a single OFDM subcarrier in a single slot given by Eq. (29) having the finite extent with |ck| ∼ ε (recall that Ωk = 2π (k − 1)/T) [Fig. 4(a)]. Then, we implement the pre-distortion of the NS associated with such an input using the results given in the previous section. Since, in this case, we have the exact analytical expression for the reflection coefficient, Eq. (31), we can simultaneously test the validity of the r(ξ) expansion given in Subsection 3.1. Indeed, one can expand the exact formula (31) in terms of the pulse amplitude ck to ensure that the general expansion formulas, Eqs. (19) and (20), are correct.
First, for such an input, the first-order NS term is given byEq. (21))
For the nonlinear addition ∼ ε3, we have(37) at ξ = −Ωk/2 is (εT)3, and therefore, we see that the expansion takes place in terms of the pulse area εT; so, it is actually valid if εT is small.
Now, we can calculate the pre-distortion correction in the NS domain using Eqs. (24) and the explicit expression (37), and the corresponding correction in the time domain is then obtained by using the backward FT of Eq. (25). The resulting function s(t) is given in Fig. 4(b) for the case T = 1, ε = 0.5, where we used the OFDM subcarrier with k = 2, i.e., with Ω2 = 2π [see the profile of q(t) in Fig. 4(a)]. Interestingly, the resulting profile of s(t) is asymmetric with respect to time axis origin, in contrast to the obvious symmetry of the input pulse.
In Fig. 4(c), we present the comparison of the absolute errors, |q(L, t) − q(0, t)|, for the following three cases.
- First, we numerically calculated the NS N(ω) for the input given by Eq. (29), pre-distorted it using Eq. (37), and performed the backward FT to obtain the new preprocessed waveform profile qs(0, t) in the time domain. Then, we computed the absolute error |qs(0, t) − q(t)| (red dots, almost indistinguishable on the figures): As expected, the error is small . The sharp peaks at the borders of the pulse base are due to the aliasing caused by the abrupt variation of the profile q(0, t).
- Then, we propagated the pre-distorted profile qs(0, t) to the normalized distance z = 1 (corresponding to L = 4000km for the normalizations given in Subsection 2.1), thus receiving qs(L, t). This quantity was used again to compute the NS at distance L, N(L, ω). After that, we removed linear dispersion from the resulting NS using Eq. (10) and performed the backward FT to obtain q̃s(0, t). The absolute error |q̃s(0, t) − q(0, t)| in this case is plotted in Fig. 4, represented by blue points. We see that the error almost precisely coincided with the initial error (expected error, red dots), as it should be, with very slight deviations occurring because of numerical errors.
- In the third run, we launched the pure initial profile q(t) without any pre-distortion and then simply removed the linear dispersion after the signal passed the distance z = 1, getting the nonlinearity-distorted quantity q̃(0, t). The black points indicate the absolute error |q̃(0, t) − q(t)| obtained without the pre-distortion.
5.2. Nonlinear pre-distortion of a general finite OFDM-modulated pattern
Now, consider a finite sequence of OFDM slots, each having length T. The input pulse has the form
Again, the leading order term in the r(ξ) expansion can be straightforwardly obtained by generalizing Eq. (35):
In the next order, we have different contributions:(43) is nonsingular for any ξ. The second term in Eq. (41) is the contribution from the inter-slot (IS) interference, the expression for each partial entry is given by Eq. (40). The two last terms are responsible for the mixed IC/IS interference effects, and the expressions for and are given by
Now, using Eqs. (41) and (25), we can calculate the precompensation profile s(t). The absolute errors associated with the propagation of the precompensated OFDM versus the error induced by the nonlinear distortions of the same OFDM sequence are summarized in Fig. 5. We can see that the error is relatively low (notice the aliasing contribution at the ends of each slot) in spite of the extended total duration of the pulse (for the simulation in Fig. 5, we used 3 OFDM slots; so, the duration is now 3T). This is because the rapid abrupt changes of the OFDM-coded pattern still arrest soliton formation and the NS remains similar to the linear spectrum.
We also checked the preprocessing error for the longer OFDM sequence and found that the precompensation worked relatively well up to 10 OFDM slots with normalized parameters: slot duration T = 1 and |cαk| = 0.1.
6. Other input profiles, noisy propagation and general comments on NS modulation beyond the perturbative approach
6.1. Precompensation of other profiles: WDM input with Gaussian base
We have so far considered only OFDM-type input, where the form of NS expansion (20) and the corresponding precompensation profile (23), (25), can be obtained analytically. However, the method described in Section 3 is general and can be applied to arbitrary input forms, with the only requirement that the input pulse does not generate singularities in the NS expansion, i.e., as before, we are far from soliton creation threshold. To give an example of such a situation, we consider a finite WDM sequence based on Gaussian pulses. The input pulse in this case has the form(46) the explicit analytical form of the precompensation profile (25) cannot be obtained, and thus we resort to the numerical procedure : first, we calculate numerically the NS N(ω) associated with q(t) in the form (46), and then the spectral profile of the pre-compensation correction S(ω) is given by S(ω) = N(ω) −Q(ω), where Q(ω) is the linear FT of Eq. (46). The descired pre-compensation profile in the time-domain is given by the backward FT of S(ω), and, same as it was done in previous sections, one can evaluate the errors associated with the propagation of precompensated input (cf. Fig. 5) and compare them with the errors arising in the propagation without precompensation.
The absolute errors calculated for the propagation of the precompensated Gaussian-based WDM versus the error induced by the nonlinear distortions of the same symbol sequence as in Fig. 5 are summarized in Fig. 6. In the simulation we used 3 slots with the QPSK encoding of cαk, the absolute value of each coefficient c = 0.15 (cf. Fig. 5); the extent of each Gaussian from (46) is ρ = 0.3 (in normalized units), the (normalized) distance between the adjacent Gaussians (slot duration) τα+1 − τα = 1. We see that the error now is even lower than that for the QPSK-OFDM with the same parameters, Fig. 5 (the aliasing errors are absent at all since the Gaussians are smooth, in contrast to the rectangular OFDM bases). This is because the Gaussians with ρ = 0.3 contain less energy compared to the OFDM with the same parameters, and the remaining terms from the expansions given in Section 3, responsible for the error of the precompensated pulse (i.e. for the difference between the NS and FT), are effectively smaller.
To end up, we see that the spectral equalization method works well for different kinds and power of input profiles provided that the conditions of the validity of expansions from Section 3 are fulfilled.
6.2. Channel with noise
Let us study the robustness of our method against the amplifier spontaneous emission (ASE). In our case we employ the model with distributed Raman amplification, and consider stochastic NLSE  with the additive symmetric white Gaussian noise (AWGN) term η (t, z):7]:
Now, to study how the noise influences the interrelation of errors associated with the pre-compensated and undistorted signal, we take the example input from Subsection 5.2: 3 slots of 10-mode QPSK-OFDM with the same parameters as were used in Fig. 5, and launch the undistorted and precompensated inputs into noisy channel (48). As there possibly exists a very small probability that for certain critical values of signal amplitude the addition of noise can cause creation of solitonic components in the received signal, one must chose the signal amplitude even weaker than the singularity condition; however, this situation is, generally, highly unlikely as the addition of noise usually arrests the creation of new solitons rather than helps it [44, 45, 48]. However we note that, when performing the Monte-Carlo simulations of signal propagation in noisy model, Eq. (48), and evaluating average errors over several runs (i.e. over several realizations of the AWGN η(z, t)), each randomly-coded set of OFDM coefficients cαk used in a certain run produces its own particular input profile q(t) and, hence, its own respective output error picture. Because of this the average local absolute error, like the one depicted in Fig. 5, is not informative, as its profile strongly depends on the particular realization of cαk and, thus, disguises the ASE-induced distortions themselves. Therefore we take the integral absolute error for each particular realization of cαk as the measure of signal corruption (here the corruption occurs already due to both noise influence and nonlinearity). For the 3 OFDM slots it means that we take the integral of local absolute error (the error definition is given in Subsection 5.1) over the full signal extent, i.e. from −1.5 to 1.5 in normalized time, evaluating it for precompensated and undistorted inputs separately (see Fig. 5). Next we evaluate the average of this quantity over several runs as a function of i) noise intensity Γ and ii) propagation length L.
The results of our simulations are shown in Fig. 7, where the averaging over 50 realizations was made for each point and we plotted separate graphs for L = 2000km, Fig. 7(a), and L = 4000km, Fig. 7(b). One can observe that even for the noisy channel the transmission based on the precompensated signal outperforms the ordinary transmission (the linear dispersion at the end point was, as before, removed): the average integral absolute error associated with the precompensated inputs (red dots) is lower than that for the non-preprocessed inputs (black squares) for every value of the noise intensity Γ studied. We notice, however, that as far as we operate in the quasilinear regime, the transmission is effectively noise-dominated, so that the nonlinearity induced absolute error, compensated by the preprocessing, is several times lower than the noise-induced error even for the lowest value of noise intensity used, 2Γ = 10−21. In spite of this, the usage of precompensation obviously contributes towards the improvement of the transmission and decrease of distortions in the noisy channel, and our method is robust against the ASE-induced corruption.
6.3. General remarks on the NS-based transmission
In our current study, due to the validity requirements of our spectra equalization method, Section 3, we have addressed the case of quasilinear transmission, where the power of the input is sufficiently low and soliton formation is arrested. However, we want to emphasize that even if the nonlinearity strength and pattern length are such that a simple perturbation-based approach of Section 3 fails, the general idea of the transmission based on NS can still have a significant practical advantage over split-step back-propagation (BP) compensation schemes . Indeed, in any split step back-propagation scheme the number of operations needed for one elementary step scales as N logN (where N is the resolution in time domain). This, however, is compounded by the need to adjust the coordinate step to compensate for the nonlinear phase shift (this is particularly important for the high pulse power), so that the total number of operation scales as M · N · logN, where M ≫ 1 is the number of coordinate steps that grows linearly with distance and signal power. In the IST scheme the propagation of the spectral data is the trivial phase rotation so M = 1 always, independent of nonlinearity or initial power. As for the direct and inverse NFT, the former requires solving the ZS problem (5) (or a first order Ricatti equation for the scattering coefficient, see e.g. ): well-developed methods based on e.g. the piece-wise linear approximation (peeling) of the input q(t)  usually take order of N2 operations. Other methods, based on the reformulation of ZS problem (5) as a matrix eigenvalue problem (see e.g. the part II of ), can require a larger number of operations but provide some advantages in the accuracy or more efficiency in the search for complex ZS eigenvalues (solitons). Note that the processing schemes utilizing the NFT have already been proposed [52, 53], although those rather dealt with complex ZS eigenvalues responsible for solitons.
The inverse NFT can be solved using ∼ N2 operations, by employing fast Gelfand-Levitan-Marchenko solvers, see e.g.  and references therein. So the computational speed-up of the NFT over conventional split step BP is ∼ (M logN)/N operations. To make a real world estimate we have chosen the latest state-of-the-art massively parallel simulations recently reported in  with N = 211 and M = 8000. Even with the reported impressive speed-up of 30 due to highly efficient parallel algorithm used in the above paper, the potential speed-up of the fully numerical NFT scheme is still 8000 × 11/(30 × 211) ≈ 1.5, i.e. even in the quasilinear regime and for a relatively small number of z-steps considered in , a simple serial numerical NFT-based scheme performs at least as good as a massively parallel direct propagation algorithm. And at higher values of nonlinearity where the required number of coordinate steps, M, is about to grow, the potential effectiveness of the NFT-based modelling schemes is likely to increase. What is most important here is the fact that both algorithms for direct and inverse NFT can be modified to include the appearance of discrete spectra at no extra computational cost. This certainly proves the viability of NFT even beyond the domain of validity of perturbative iteration scheme proposed in Section 3.
Finally we notice that the matrix methods for both the ZS problem , and for the GLM equations  can be effectively reduced to the Toeplitz matrix inversions. There already exist “superfast” stable algorithms allowing one to perform the Toeplitz matrix inversion by using only ∼ N logN operations [56,57]. Therefore, in view of the arguments given above, we believe that the NFT-based processing and transmission methods can be highly competitive, no less effective than ordinary BP, and even seriously outperform the traditional split-step BP processing.
One of the ideas that allows one to use an encoded continuous NS at larger powers, i.e. when the iterative approach of Section 3 does not rapidly converge failing to produce the correct approximation of NS in a small number of steps, can be the method of nonlinear inverse syntheses. The strategy of this method is similar to that widely used in the fiber Bragg gratings design : one starts from the given (encoded) NS profile in the spectral domain (and corresponding reflection coefficient r(ξ), Eq. (7)) and then reconstructs the corresponding profile in time domain by employing the backward NFT, i.e. by solving the corresponding GLM equation. For the focusing NLSE in the absence of discrete ZS spectrum (solitons), the GLM equation (in normalized units) has the form [26–28]:equation (50) for K(t, t′), the sought profile in time domain is eventually recovered as q(t) = −2K(t, t′)|t′→t. After the particular form of q(t) corresponding to the encoded NS has been found, one launches this profile into the NLSE. At the receiver we read the output field q(t, L), perform the forward NFT, thus getting the NS N(ω; L), wind out the accumulated dispersion, and recover the information encoded without nonlinearity-induced distortions. The detailed analysis of this approach is beyond the scope of this paper and will be published elsewhere.
In this paper, we have examined an application of the inverse scattering transform (nonlinear Fourier transform) to digital signal processing in the coherent communication channel modeled by the integrable NLSE channel. We have demonstrated that the nonlinear spectrum (NS) associated with this channel can be used for the improvement of the quality of long-haul information transmission. Without loss of generality, we have illustrated the use of the NS in the case of the OFDM-coded signal. The advantage of using the NS is that the encoded spectral content undergoes just a trivial transformation, similar to the dispersive evolution of spectral component in the linear problem. We would like to stress that the linear spectrum in a nonlinear channel is distorted and its evolution differs drastically from the evolution in a linear dispersive channel. In contrast, the change in the spectral phase of the NS can be easily eliminated at the receiver side, and thus, the informational content can be straightforwardly restored.
To explain the method, we first calculated analytically the NS in the case of a single tone of OFDM and studied the threshold for the nucleation of solitons. We demonstrated that the threshold in the considered case coincided with that for a pure rectangular pulse of the same energy. For the general situation of several OFDM slots serving as an NLSE input, we showed that NS is qualitatively close to the linear signal spectrum when the energy of the signal is below the soliton-creation threshold: In this situation, the NS can be used for information encoding in the same way as we do with the linear spectrum. The important nontrivial point here is a calculated nonlinear correction to the linear spectrum.
The approach considered in our paper allows one to perform consecutive elimination of the nonlinearity-induced interference effects appearing during the transmission through a nonlinear channel, when the input energy of the single information-carrying pulse is not very high. When the input energy of the pulses block is above the soliton threshold, a more sophisticated modulation can be needed, e.g. separation of small individual information-bearing blocks, such that the shape-sensitive soliton creation threshold is not reached; but the straightest way in this case is to employ the inverse syntheses method briefly described at the end of Subsection 6.3.
Using the expansion of the NS, we propose the nonlinear pre-distortion technique, which allows one to equalize the NS with the linear spectrum of the undistorted input signal. Then, the informational content encoded on the NS can be recovered at the receiver with standard linear methods after the linear accumulated dispersion has been compensated. We found that our technique worked well even for a relatively long OFDM sequence, provided that the input pulse does not nucleate a soliton. We have also demonstrated that our approach is general and works well also for the Gaussian-based WDM input, where the NS expansion cannot be computed analytically. Further, we analyzed the robustness of our approach against the noise associated with amplifier spontaneous emission and revealed that even though the quasilinear regime is noise-dominated, the average error associated with the pre-compensated inputs is always lower than that for non-processed transmission. We hope that our work represents an important first step in the development of efficient digital signal techniques based on the inverse scattering transform capable of eliminating nonlinear distortion of the optical signal in the NLSE channel via pre-or post-processing using the integrability of the NLSE (2). The forward NFT itself consists in the solution of just the linear system – the ZS spectral problem (5). We note that the squares of ZS eigenfunctions (Jost solutions) constitute an orthogonal basis (5), and this paves the way for various perturbative approaches for the consideration of noise influence and other nearly-integrable (weakly perturbed NLSE) cases. We anticipate that with the further development of practical implementation of the proposed approaches, certain nonlinear fiber communication channels can be effectively linearized. The numerical complexity of the transmission based on NS and, generally, of the NFT-based methods, can be of the same order of even lower than that for traditional slit-step back propagation.
We would like to thank Keith Blow and Ildar Gabitov for valuable discussions. The work was supported by the EPSRC project UNLOC (Unlocking the Capacity of Optical Communications) EP/J017582/1. The support of the Russian Ministry of Education and Science, European Research Council and Marie Curie IRSES program is also acknowledged.
References and links
1. A. Hasegawa and Y. Kodama, Solitons in Optical Communications (Oxford University Press, 1995).
2. G. P. Agrawal, Nonlinear Fiber Optics, 4th ed. (Elsevier, 2006).
4. A. Mecozzi, “Limits to the long haul coherent transmission set by the Kerr nonlinearity and noise of in-line amplifiers,” J. Lightwave Technol. 12, 1993–2000 (1994). [CrossRef]
6. E. B. Desurvire, “Capacity demand and technology challenges for lightwave systems in the next two decades,” J. Lightwave Technol. 24, 4697–4710 (2006). [CrossRef]
8. A. D. Ellis, J. Zhao, and D. Cotter, “Approaching the non-linear Shannon limit,” J. Lightwave Technol. 28, 423–433 (2010). [CrossRef]
10. R. I. Killey and C. Behrens, “Shannon’s theory in nonlinear systems,” J. Mod. Opt. 58, 1–10 (2011). [CrossRef]
11. A. Mecozzi and R.-J. Essiambre, “Nonlinear Shannon limit in pseudo-linear coherent systems,” J. Lightwave Technol. 30, 2011–2024 (2012). [CrossRef]
13. see Impact of nonlinearities on fiber-optic communication systems, Ed. S. Kumar, (Springer, 2011). [CrossRef]
14. A. Splett, C. Kurtzke, and K. Petermann, “Ultimate transmission capacity of amplified optical fiber communication systems taking into account fiber nonlinearities,” in Tech. Digest of European Conference on Optical Communication , 1993, paper MoC2.4.
15. D. Rafique and A. D. Ellis, “Impact of signal-ASE four-wave mixing on the effectiveness of digital back-propagation in 112 Gb/s PM-QPSK systems,” Opt. Express 19, 3449–3454 (2011). [CrossRef] [PubMed]
16. E. Ip and J. Kahn, “Compensation of dispersion and nonlinear impairments using digital backpropagation,” J. Lightwave Technol. 26, 3416–3425 (2008). [CrossRef]
17. D. S. Millar, S. Makovejs, C. Behrens, S. Hellerbrand, R. I. Killey, P. Bayvel, and S. J. Savory, “Mitigation of fiber nonlinearity using a digital coherent receiver,” IEEE J. Sel. Top. Quantum Electron. 16, 1217–1226 (2010). [CrossRef]
18. M. Nazarathy, J. Khurgin, R. Weidenfeld, Y. Meiman, P. Cho, R. Noe, I. Shpantzer, and V. Karagodsky, “Phased-array cancellation of nonlinear FWM in coherent OFDM dispersive multi-span links,” Opt. Express 15, 15777–15810 (2008). [CrossRef]
19. S. Shieh and I. Djordjevic, OFDM for Optical Communications (Academic Press, 2010).
21. L. B. Y. Du and A. J. Lowery, “Pilot-based XPM nonlinearity compensator for CO-OFDM systems,” Opt. Express 19, B862–B869 (2011). [CrossRef]
23. X. Liu, A. R. Chraplyvy, P. J. Winzer, R. W. Tkach, and S. Chandrasekhar, “Phase-conjugated twin waves for communication beyond the Kerr nonlinearity limit,” Nat. Photonics 7, 560–568 (2013). [CrossRef]
24. H. Cheng, W. Li, Y. Fan, Z. Zhang, S. Yu, and Z. Yang, “A novel fiber nonlinearity suppression method in DWDM optical fiber transmission systems with an all-optical pre-distortion module,” Opt. Comm. 290, 152–157 (2013). [CrossRef]
25. V. E. Zakharov and A. B. Shabat, “Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media,” Soviet Physics-JETP 34, 62–69 (1972).
26. M. J. Ablowitz, D. J. Kaup, A. C. Newell, and H. Segur, “The inverse scattering transform-Fourier analysis for nonlinear problems,” Stud. Appl. Math. 53, 249–315 (1974).
27. A. C. Newell, Solitons in Mathematics and Physics (SIAM, 1985). [CrossRef]
28. V. E. Zakharov, S. V. Manakov, S. P. Novikov, and L. P. Pitaevskii, Theory of Solitons. The Inverse Scattering Method(Colsultants Bureau, 1984).
29. M. J. Ablowitz and H. Segur, Solitons and the Inverse Scattering Transform (SIAM, 1981). [CrossRef]
30. M. I. Yousefi and F. R. Kschischang, “Information transmission using the nonlinear Fourier transform, Part I: Mathematical tools,” http://arxiv.org/abs/1202.3653, submitted to IEEE Trans. Inf. Theory, Feb. 2012; “Part II: Numerical methods,” Apr. 2012, online: http://arxiv.org/abs/1204.0830, submitted to IEEE Trans. Inf. Theory, Feb. 2012; “Part III: Spectrum modulation,” Feb. 2013, online: http://arxiv.org/abs/1302.2875, submitted to IEEE Trans. Inf. Theory, Feb. 2013.
31. D. J. Kaup, “Closure of the squared Zakharov-Shabat eigenstates,” J. Math. Anal. Appl. 54, 849–864 (1976). [CrossRef]
32. L. F. Mollenauer and J. P. Gordon, Solitons in Optical Fibers: Fundamentals and Applications (Academic Press, 2006).
33. Yu. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic Press, 2003).
34. A. Hasegawa and T. Nyu, “Eigenvalue communication,” J. Lightwave Technol. 11, 395–399 (1993). [CrossRef]
35. J. E. Prilepsky, S. A. Derevyanko, and S. K. Turitsyn, “Lattice approach to the dynamics of phase-coded soliton trains,” J. Phys. A: Math. Theor. 45, 025202 (2012). [CrossRef]
37. J. Satsuma and N. Yadjima, “Initial value problems of one-dimensional self-modulation of nonlinear waves in dispersive media,” Progr. Theor. Phys. Suppl. 55, 284–306 (1974). [CrossRef]
38. R. Feced, M. N. Zervas, and M. A. Muriel, “An efficient inverse scattering algorithm for the design of nonuniform fiber Bragg gratings,” IEEE J. Quant. Electron. 35, 1105–1115 (1999). [CrossRef]
39. W. Shieh, X. Yi, Y. Ma, and Q. Yang, “Coherent optical OFDM: has its time come? [Invited],” J. Opt. Networking 7, 234–255 (2008). [CrossRef]
40. R. Bouziane, P. Milder, R. Koutsoyannis, Y. Benlachtar, J. C. Hoe, M. Puschel, M. Glick, and R. I. Kille, “Design studies for ASIC implementations of 28 GS/s optical QPSK- and 16-QAM-OFDM transceivers,” Opt. Express 19, 20857–20864 (2011). [CrossRef] [PubMed]
41. S. Herbst, S. Bayer, H. Wernz, and H. Griesser, “21-GHz single-band OFDM transmitter with QPSK modulated subcarriers,” in Optical Fiber Communication Conference and Exposition (OFC/NFOEC), 2011 and the National Fiber Optic Engineers Conference, IEEE, 2011, paper OMS3.
42. A.H. Gnauk and P. J. Winzer, “Optical phase-shift-keyed transmission,” J. Lightwave Technol. 23, 115–130 (2005). [CrossRef]
43. S.V. Manakov, “On nonlinear Fraunhofer diffraction,” Soviet Physics-JETP 38, 693–696 (1974).
44. S. A. Derevyanko and J. E. Prilepsky, “Soliton generation from randomly modulated return-to-zero pulses,” Opt. Comm. 281, 5439–5443 (2008). [CrossRef]
45. S. A. Derevyanko and J. E. Prilepsky, “Random input problem for the nonlinear Schrodinger equation,” Phys. Rev. E 78, 046610 (2008). [CrossRef]
46. M. Klaus and J. K. Shaw, “Purely imaginary eigenvalues of Zakharov-Shabat systems,” Phys. Rev. E 65, 036607 (2002). [CrossRef]
47. J. E. Prilepsky, S. A. Derevyanko, and S. K. Turitsyn, “Conversion of a chirped Gaussian pulse to a soliton or a bound multisoliton state in quasilossless and lossy optical fiber spans,” J. Opt. Soc. Am. B 24, 1254–1261 (2007). [CrossRef]
48. S. K. Turitsyn and S. A. Derevyanko, “Soliton-based discriminator of noncoherent optical pulses,” Phys. Rev. A 78, 063819 (2008). [CrossRef]
49. S. Vergeles and S. K. Turitsyn, “Optical rogue waves in telecommunication data streams,” Phys. Rev. A 83, 061801(R) (2011). [CrossRef]
51. G. Boffetta and A. R. Osborne, “Computation of the direct scattering transform for the nonlinear Schroedinger equation,” J. Comput. Phys. 102, 252–264 (1992). [CrossRef]
52. S. Oda, A. Maruta, and K. Kitayama, “All-Optical Quantization Scheme Based on Fiber Nonlinearity,” IEEE Photon. Technol. Lett. 16, 587–589 (2004). [CrossRef]
53. H. Terauchi and A. Maruta, “Eigenvalue Modulated Optical Transmission System Based on Digital Coherent Technology,” in 18th OptoElectronics and Communications Conference held jointly with 2013 International Conference on Photonics in Switching (OECC/PS), IEICE, 2013, paper WR2-5.
54. O. V. Belai, L. L. Frumin, E. V. Podivilov, and D. A. Shapiro, “Efficient numerical method of the fiber Bragg grating synthesis,” J. Opt. Soc. Am. B 24, 1451–1457 (2007). [CrossRef]
55. A. O. Korotkevich and P. M. Lushnikov, “Proof-of-concept implementation of the massively parallel algorithm for simulation of dispersion-managed WDM optical fiber systems,” Opt. Lett. 36, 1851–1853 (2011). [CrossRef] [PubMed]
56. M. Van Barel, G. Heinig, and P. Kravanja, “A stabilized superfast solver for nonsymmetric Toeplitz systems,” SIAM J. Matrix Anal. Appl. 23, 494–510 (2001). [CrossRef]
57. S. Chandrasekaran, M. Gu, X. Sun, J. Xia, and J. Zhu, “A superfast algorithm for Toeplitz systems of linear equations,” SIAM J. Matrix. Anal. Appl. 29, 1247–1266 (2007). [CrossRef]