Signaling on the Continuous Spectrum of
Nonlinear Optical Fiber
I MAN TAVAKKOLNIA *

AND

M AJID S AFARI

Institute for Digital Communications, School of Engineering, University of Edinburgh, EH9 3JL,
Edinburgh, UK

arXiv:1704.05537v2 [cs.IT] 21 Sep 2017

* i.tavakkolnia@ed.ac.uk

Abstract: This paper studies different signaling techniques on the continuous spectrum (CS) of
nonlinear optical fiber defined by nonlinear Fourier transform. Three different signaling techniques
are proposed and analyzed based on the statistics of the noise added to CS after propagation along
the nonlinear optical fiber. The proposed methods are compared in terms of error performance,
distance reach, and complexity. Furthermore, the effect of chromatic dispersion on the data
rate and noise in nonlinear spectral domain is investigated. It is demonstrated that, for a given
sequence of CS symbols, an optimal bandwidth (or symbol rate) can be determined so that the
temporal duration of the propagated signal at the end of the fiber is minimized. In effect, the
required guard interval between the subsequently transmitted data packets in time is minimized
and the effective data rate is significantly enhanced. Moreover, by selecting the proper signaling
method and design criteria a distance reach of 7100 km is reported by only singling on CS at a
rate of 9.6 Gbps.
© 2018 Optical Society of America.
OCIS codes: (060.1660) Coherent communications; (060.2330) Fiber optics communications; (070.4340) Nonlinear
optical signal processing.

References and links
1. J. C. Cartledge, F. P. Guiomar, F. R. Kschischang, G. Liga, and M. P. Yankov, “Digital signal processing for fiber
nonlinearities,” Opt. Express 25(3), 1916–1936 (2017).
2. R. Essiambre, G. Kramer, P. J. Winzer, G. J. Foschini, and B. Goebel, “Capacity limits of optical fiber networks,”
IEEE/OSA J. Lightw. Technol. 28(4), 662–701 (2010).
3. P. J. Winzer and D. T. Neilson, “From scaling disparities to integrated parallelism: A decathlon for a decade,”
IEEE/OSA J. Lightw. Technol. 35(5), 1099–1115 (2017).
4. M. Yousefi and F. Kschischang, “Information transmission using the nonlinear Fourier transform, part I: Mathematical
tools,” IEEE Trans. Inform. Theory 60(7), 4312–4328 (2014).
5. M. Yousefi and F. Kschischang, “Information transmission using the nonlinear Fourier transform, part II: Numerical
methods,” IEEE Trans. Inform. Theory 60(7), 4329–4345 (2014).
6. M. Yousefi and F. Kschischang, “Information transmission using the nonlinear Fourier transform, part III: Spectrum
modulation,” IEEE Trans. Inform. Theory 60(7), 4346–4369 (2014).
7. S. K. Turitsyn, J. E. Prilepsky, S. T. Le, S. Wahls, L. L. Frumin, M. Kamalian, and S. A. Derevyanko, “Nonlinear
Fourier transform for optical data processing and transmission: Advances and perspectives,” Optica 4(3), 307–322
(2017).
8. J. E. Prilepsky, S. A. Derevyanko, K. J. Blow, I. Gabitov, and S. K. Turitsyn, “Nonlinear inverse synthesis and
eigenvalue division multiplexing in optical fiber channels,” Phys. Rev. Lett. 113(1), 013901 (2014).
9. M. I. Yousefi and X. Yangzhang, “Nonlinear frequency-division multiplexing,” arXiv preprint arXiv:1603.04389
(2016).
10. I. Tavakkolnia and M. Safari, “Signalling over nonlinear fibre-optic channels by utilizing both solitonic and radiative
spectra,” in European Conference on Networks and Communications (EuCNC) (2015), pp. 103–107.
11. S. Hari, M. I. Yousefi, and F. R. Kschischang, “Multieigenvalue communication,” IEEE/OSA J. Lightw. Technol.
34(13), 3110–3117 (2016).
12. S. Hari and F. R. Kschischang, “Bi-directional algorithm for computing discrete spectral amplitudes in the NFT,”
IEEE/OSA J. Lightw. Technol. 34(15), 3529–3537 (2016).
13. H. Buelow, V. Aref, and W. Idler, “Transmission of waveforms determined by 7 eigenvalues with PSK-modulated
spectral amplitudes,” in European Conference on Optical Communication (2016), pp. 1–3.
14. H. Bülow, “Experimental demonstration of optical signal detection using nonlinear Fourier transform,” IEEE/OSA J.
Lightw. Technol. 33(7), 1433–1439 (2015).
15. Z. Dong, S. Hari, T. Gui, K. Zhong, M. I. Yousefi, C. Lu, P.-K. A. Wai, F. R. Kschischang, and A. P. T. Lau, “Nonlinear
frequency division multiplexed transmissions based on NFT,” IEEE Photon. Technol. Lett. 27(15), 1621–1623 (2015).

16. T. Gui, T. H. Chan, C. Lu, A. P. T. Lau, and P. K. A. Wai, “Alternative decoding methods for optical communications
based on nonlinear Fourier transform,” IEEE/OSA J. Lightw. Technol. 35(9),1542–1550 (2017).
17. N. A. Shevchenko, S. A. Derevyanko, J. E. Prilepsky, A. Alvarado, P. Bayvel, and S. K. Turitsyn, “On the capacity of
the noncentral Chi-channel with applications to soliton amplitude modulation,” arXiv preprint arXiv:1609.02318
(2016).
18. N. A. Shevchenko, J. E. Prilepsky, S. A. Derevyanko, A. Alvarado, P. Bayvel, and S. K. Turitsyn, “A lower bound on
the per soliton capacity of the nonlinear optical fibre channel,” in Information Theory Workshop (2015), pp. 104–108.
19. V. Aref, S. T. Le, and H. Buelow, “Demonstration of fully nonlinear spectrum modulated system in the highly
nonlinear optical transmission regime,” in European Conference on Optical Communication (2016), pp. 1–3.
20. S. Wahls and H. V. Poor, “Fast numerical nonlinear Fourier transforms,” IEEE Trans. Inform. Theory 61(12),
6957–6974 (2015).
21. S. Wahls and V. Vaibhav, “Fast inverse nonlinear Fourier transforms for continuous spectra of Zakharov-Shabat type,”
arXiv preprint arXiv:1607.01305 (2016).
22. S. Civelli, L. Barletti, and M. Secondini, “Numerical methods for the inverse nonlinear Fourier transform,” in
International Workshop on Digital Communications (2015), pp. 13–16.
23. M. Kamalian, J. E. Prilepsky, S. T. Le, and S. K. Turitsyn, “Periodic nonlinear Fourier transform for fiber-optic
communications, part F: theory and numerical methods,” Opt. Express 24(16), 18353–18369 (2016).
24. M. Kamalian, J. E. Prilepsky, S. T. Le, and S. K. Turitsyn, “Periodic nonlinear Fourier transform for fiber-optic
communications, part II: eigenvalue communication,” Opt. Express 24(16), 18370–18381 (2016).
25. S. A. Derevyanko, J. E. Prilepsky, and S. K. Turitsyn, “Capacity estimates for optical transmission based on the
nonlinear Fourier transform,” Nature Commun. 7, 12710 (2016).
26. I. Tavakkolnia and M. Safari, “Capacity analysis for the continuous spectrum of nonlinear optical fiber,” IEEE/OSA J.
Lightw. Technol. 35(11), 2086–2097 (2017).
27. X. Yangzhang, M. I. Yousefi, A. Alvarado, D. Lavery, and P. Bayvel, “Nonlinear frequency-division multiplexing in
the focusing regime,” in Optical Fiber Communication Conference (2017), p. Tu3D-1.
28. S. T. Le, J. E. Prilepsky, and S. K. Turitsyn, “Nonlinear inverse synthesis for high spectral efficiency transmission in
optical fibers,” Opt. Express 22(22), 26720–26741 (2014).
29. S. T. Le, J. E. Prilepsky, and S. K. Turitsyn, “Nonlinear inverse synthesis technique for optical links with lumped
amplification,” Opt. Express 23(7), 8317–8328 (2015).
30. S. T. Le, I. D. Philips, J. E. Prilepsky, P. Harper, A. D. Ellis, and S. K. Turitsyn, “Demonstration of nonlinear inverse
synthesis transmission over transoceanic distances,” IEEE/OSA J. Lightw. Technol. 34(10), 2459–2466 (2016).
31. S. T. Le, J. E. Prilepsky, P. Rosa, J. D. Ania-Castanon, and S. K. Turitsyn, “Nonlinear inverse synthesis for optical
links with distributed Raman amplification,” IEEE/OSA J. Lightw. Technol. 34(8), 1778–1786 (2016).
32. I. Tavakkolnia and M. Safari, “Dispersion pre-compensation for NFT-based optical fiber communication systems,” in
Conference on Lasers and Electro-Optics (2016), p. SM4F4.
33. M. J. Ablowitz, B. Prinari, and A. D. Trubatch, Discrete and continuous nonlinear Schrödinger systems (Cambridge
University, 2004).
34. G. P. Agrawal, Nonlinear fiber optics (Academic, 2013), 5th ed.
35. M. Sorokina, S. Sygletos, and S. Turitsyn, “Ripple distribution for nonlinear fiber-optic channels,” Opt. Express
25(3), 2228–2238 (2017).
36. A. S. Skidin, O. S. Sidelnikov, M. P. Fedoruk, and S. K. Turitsyn, “Mitigation of nonlinear transmission effects for
OFDM 16-QAM optical signal using adaptive modulation,” Opt. Express 24(26), 30296–30308 (2016).
37. A. A. Farid and S. Hranilovic, “Channel capacity and non-uniform signalling for free-space optical intensity channels,”
IEEE J. Sel. Areas Commun. 27(9), 1553–1563 (2009).
38. M. Bartlett, “The square root transformation in analysis of variance,” Supplement to the Journal of the Royal Statistical
Society pp. 68–78 (1936).
39. J. H. Curtiss, “On transformations used in the analysis of variance,” The Annals of Mathematical Statistics 14,
107–122 (1943).
40. M. Safari, “Efficient optical wireless communication in the presence of signal-dependent noise,” in International
Conference on Communication Workshop (ICCW) (2015), pp. 1387–1391.
41. A. Tsiatmas, F. M. Willems, and C. P. Baggen, “Square root approximation to the Poisson channel,” in International
Symposium on Information Theory (2013), pp. 1695–1699.
42. I. Tavakkolnia and M. Safari, “Variance normalizing transform for performance improvement in radio-over-fiber
systems,” in European Conference on Networks and Communications (2016), pp. 255-259.
43. P. R. Prucnal and B. E. Saleh, “Transformation of image-signal-dependent noise into image-signal-independent
noise,” Opt. lett. 6(7), 316–318 (1981).
44. J. Rice, Mathematical Statistics and Data Analysis (Duxbury, 1995), 2nd ed.
45. I. Tavakkolnia and M. Safari, “Efficient Signalling on the Continuous Spectrum of Nonlinear Optical Fibre,” in
Conference on Lasers and Electro-Optics/Europe and the European Quantum Electronics Conference (2017).
46. L. L. Frumin, A. A. Gelash, and S. K. Turitsyn, “New approaches to coding information using inverse scattering
transform, ” Phys. Rev. Lett. 118(22), 223901 (2017).
47. S. Dimitrov, S. Sinanovic, and H. Haas, “Clipping Noise in OFDM-Based Optical Wireless Communication Systems,”
IEEE Trans. Commun. 60 (4), 1072–1081 (2012).

1.

Introduction

Extensive efforts have been made to develop new techniques for higher achievable data rates in
optical fiber communication particularly over long-haul systems, where nonlinear effects are
significant [1–3]. Systems based on the nonlinear Fourier transform (NFT) have recently attracted
attention among researchers as an approach to overcome fiber nonlinearity and to achieve the true
capacity of optical fiber links [4–7]. In the absence of any perturbation (e.g., amplifier noise), the
nonlinear channel defined by nonlinear Schrödinger (NLS) equation, is transformed by NFT to
a linear memoryless channel. As a result, nonlinear frequency division multiplexing (NFDM)
was proposed [4–6], where, rather than on the time domain, the data is mapped on either of
the two available nonlinear spectra defined by NFT, namely, continuous spectrum (CS) [8–10]
and discrete spectrum (DS) [11–18]. Simultaneous signaling on both CS and DS has been also
recently investigated [10, 19]. Due to the complex nature of the problem, still many aspects of
NFDM systems are unknown and are being investigated. For instance, recently, efficient numerical
methods have been proposed for fast forward and backward NFT [20–22], and periodic NFT was
proposed in [23, 24].
Mapping the data on CS is an approach very similar to orthogonal frequency division
multiplexing (OFDM), but with nonlinear Fourier transform instead of ordinary linear Fourier
transform. Despite the interesting linearizing feature of NFT in ideal conditions, modeling the
system in the presence of perturbation is not trivial. In [10, 25, 26], it was demonstrated that the
CS channel is almost memoryless, yet the noise in nonlinear spectral domain was observed to
be non-Gaussian and signal-dependent. Nevertheless, remarkable capacity for CS channel was
predicted [25, 26], which promises the possibility of achieving the ultimate capacity of optical
fiber if all degrees of freedom (i.e., CS and DS) are exploited efficiently.
Different techniques were used for mapping data symbols on CS and detecting the received
symbols. In [9, 10, 26, 27], data is directly mapped on the CS and transformed to the time
domain using inverse nonlinear Fourier transform (INFT) before transmission over the fiber. At
the receiver side, NFT is applied to transfer back to the signal space in CS, where detection
is performed. In a series of other research papers [8, 28–31], nonlinear inverse synthesis was
proposed and investigated numerically and experimentally. In this method, the linear spectrum of
a predefined time domain signal (e.g., OFDM or Nyquist) rather than the signal itself is mapped
on CS. The important point is that each of the aforementioned methods manifests different
behavior and performance level, but this has not been investigated thoroughly in the literature.
An important characteristic of the CS channel in NFDM systems, similar to the traditional
optical fiber communication system, is the data rate limitation due to chromatic dispersion. In the
absence of DS, there is no mechanism, which can limit the induced dispersion before performing
NFT at the receiver. Thus, usually large temporal broadening is observed in NFT systems based
on CS signaling. Since interference of two neighboring data packets in time domain cannot be
removed after NFT operation, guard intervals between data packets in time domain are required to
avoid error due to the interference. For instance, a guard time of 17 ns, accounting for dispersion
induced memory plus a 20% margin, was used for transmitting a data packet of 2 ns long in [29].
In some studies [27–29], the effect of such guard intervals was not taken into account, when the
achievable data rate was calculated. However, this assumption only gives a reasonable estimate if
very long data packets are considered leading to a substantial increases in the complexity and the
delay at the receiver. In [32], we proposed a simple precompensation method that can decrease
the induced dispersion by 50%. Nevertheless, a shorter guard interval is still required and its
effect on the design of the NFDM system needs to be studied.
In this paper, first, the effect of dispersion on NFDM systems is investigated. It is demonstrated
that better performance and higher bit rate can be achieved by properly choosing the system
parameters and minimizing the effect of dispersion. Then, by considering the specific characteristics of the noise in nonlinear spectral domain, we propose three different signaling methods

for performance improvement and compare them with the conventional method in [6] in terms
of error performance, distance reach, and complexity. Taking into account the dispersion effect,
we can investigate data rates in bits per second rather than bits per symbol. In particular, it is
demonstrated that 9.6 Gbps can be transmitted over 7100 km using 26 GHz bandwidth.
2.

Preliminaries

In this section, the channel modeling and theoretical principles of NFDM systems are briefly
discussed [4–6, 8, 10, 26, 33].
2.1.

Optical fiber communication

Propagation of the complex envelope of a narrowband optical field in a standard single-mode
fiber can be described by the stochastic nonlinear Schrödinger (NLS) equation [34]. Assuming
the fiber loss to be perfectly compensated by ideal distributed amplification, the NLS equation
can be described as
j β2 ∂ 2 Q(T, l)
∂Q(T, l)
=−
+ jγQ(T, l) |Q(T, l)| 2 + N(T, l), 0 ≤ l ≤ L,
∂l
2
∂T 2

(1)

where Q(T, l), l, and T respectively
√ represent the complex envelope of the optical field, distance
in km, and time in seconds ( j = −1). N(T, l) represents the amplified spontaneous emission
(ASE) noise added by the amplifiers, which is a white Gaussian process with autocorrelation
E{N(T, l)N ∗ (T 0, l 0 )} =

N ASE
δ(T − T 0)δ(l − l 0),
L

(2)

where δ(.) is the Dirac delta function, and N ASE = αLhνs KT is the accumulated spectral density
of noise at fiber length L. Here, β2 , γ, α, hνs , and KT are respectively chromatic dispersion
coefficient, nonlinearity parameter, fiber loss, average photon energy at carrier frequency νs , and
phonon occupancy factor. Throughout this paper, we consider the focusing case β2 < 0 without
any kind of dispersion compensation. It is also assumed that the effect of higher order dispersions
can be ignored. The NLS equation can be normalized in the form of
j

∂q(t, z) ∂ 2 q(t, z)
=
+ 2|q(t, z)| 2 q(t, z) + n(t, z),
∂z
∂t 2

(3)

by normalization rules
q=

p

γLD Q,

z=

l
,
2LD

t=

T
,
T0

(4)

where the normalizing parameter T0 can be chosen independent of other parameters. Also, we
have LD = T02 /| β2 |. Note that the variables in the normalized NLS equation (3) are unitless.
2.2.

Nonlinear Fourier transform

Using inverse scattering method [33], NFT is defined, which transforms the time domain optical
M , {λ } M }, which evolve linearly along the fiber
signal into scattering data {ρ(λ, z), {Cm (z)}m=1
m m=1
in nonlinear spectral domain. The parameter M represents the number of eigenvalues. The
continuous spectrum (CS) ρ(λ, z) is defined on the real axis λ ∈ R, and the discrete spectrum
Cm (z) is defined on the upper half complex plane λm ∈ C+ . It can be shown that, in a noise-free
scenario, the eigenvalues are preserved along the fiber, and the evolution equations for continuous
and discrete spectra are expressed as [33]
λm (z) = λm (0) = λm,

m = 1, · · · , M,

(5)

2

Cm (z) = Cm (0)e−4jλm z ,
ρ(λ, z) = ρ(λ, 0)e

m = 1, · · · , M,
−4jλ2 z

(6)

.

(7)

It can be seen from (7) that for CS the evolution equation is only a phase shift since λ ∈ R.
The (forward) NFT is performed by solving the so-called Zakharov-Shabat eigenvalue problem,
and inverse nonlinear Fourier transform (INFT) can be performed by solving Gelfand-LevitanMarchenko (GLM) integral equations [33]. A number of other methods have been also proposed
for INFT [9, 21], and there exist several numerical methods in the literature for NFT and INFT
operations [5, 9, 10, 12, 21].
The block diagram of a typical NFDM system with direct mapping of data on CS is demonstrated
in Fig.1. In this paper, we consider only CS, which correspond to the dispersive radiation waves. In
M and {λ } M , which correspond to solitonic waves, are zero throughout
other words, {Cm (z)}m=1
m m=1
this paper. In spite of the fact that mapping data only on CS is similar to OFDM, still many
aspects of it are unknown and need to be investigated. For simplicity in notation, the dependence
to distance is dropped and shown as subscript. For instance, qL (t) = q(t, L) for fiber length
L. In Fig. 1, D0 is the vector of K symbols, which are driven from the binary data using a
specific modulation format. The symbols in D0 are then converted to an oversampled continuous
waveform ρ0 (λ) using a pulse shaping filter (e.g., raised-cosine or sinc pulses). Note that since
the integral equations of INFT are solved digitally in practice, oversampling is essential for
minimizing numerical errors. Then, INFT is applied to generate the time domain signal q0 (t).
The optical transmitter (Tx) and receiver (Rx) consist of all operations needed for launching the
signal into the fiber and detecting it as a digital signal, such as (de)normalization in (4), optical
filtering, and analog/digital conversion. Standard single mode fiber (SSMF) and ideal distributed
amplification are assumed. At the receiver, after NFT operation and removing the phase shift
introduced according to (7), the noisy CS signal ρ̃ L (λ) is obtained, from which the data symbols
can be recovered. Throughout the paper, we refer to the conventional signaling method depicted
in Fig. 1 [4–6, 9, 10, 26, 27] as "direct mapping on CS" and use it as a benchmark for performance
comparisons among proposed signaling techniques.

D0

Mapping

ρ0(λ)

DL

Sampling

ρL(λ)

Equalizer
eq. (7)

INFT

ρL(λ)

q0(t)

NFT

qL(t)

Tx

Q 0 ( T)
SSMF

On CS

Rx

QL(T)

Fig. 1. Block Diagram of a NFDM system with direct mapping on CS.
We use the Ablowitz method for NFT, and discretization of Marchenko equations and using
the symmetry of the resulted matrix to implement INFT operation [10]. Moreover, the split-step
Fourier method is applied for simulating the propagation of signal along the fiber. Throughout
this paper, T0 = 0.1 ns is chosen. More details about our simulation method can be found in [26].
2.3.

Propagation of the continuous spectrum in the presence of noise

In the low noise scenario, perturbation theory can be used to determine the characteristics of
noisy CS [25, 26]. The channel model in the nonlinear spectral domain can be expressed as [6]
ρ̃ L (λ) = ρ0 (λ) + η L (λ),

(8)

where ρ̃ L (λ) = ρ L (λ) exp(4 jλ2 z) refers to noisy CS signal after phase shift removal according to
(7), and η L (λ) is a complex generally non-Gaussian noise, the variance of which is dependent
on the initial signal as E{|η L (λ)| 2 } = f [| ρ0 (λ)|]. The signal-dependency of the noise in CS was
studied by asymptotic analysis and numerical simulation in [26], and by applying perturbation
theory in [25]. It is shown that f [| ρ0 (λ)|] can be approximated by a polynomial of order 4 [25,26].
This means that signals with higher amplitudes experience higher noise variance. Furthermore, the
real and imaginary parts of the noise in CS may not be assumed independent after propagation over
the nonlinear optical fiber unlike the conventional linear channels, and they can be modeled [26]
as E{<[η L (λ)]2 } = E{=[η L (λ)]2 } = f [| ρ0 (λ)|]/2 . However, in many practical scenarios, the
noise can be assumed white, i.e., E{η L (λ)η∗L (λ 0)} = f [| ρ0 (λ)|]δ(λ − λ 0). The linearizing effect
of NFT makes the channel model of a nonlinear optical fiber much less complex in CS as in (8)
than in time domain as in (1). Nonetheless, due to the special statistical behavior of the noise,
novel techniques are required for capacity estimation as we addressed in [26] and signaling,
which is the subject of this paper.
In this paper, we use complex modulation using sinc pulses as expressed in (9), where D0i is
the ith symbol from vector D0 .
ρ0 (λ) =

K
Õ
i=1

D0i sinc(

K
K 2i − 1
λ+ −
).
Λ
2
2

(9)

The width of CS signal in nonlinear spectral domain is called nonlinear spectral width and is
denoted by Λ. Note that the NFT operation at the receiver is only required at K symbol points
Λ
λi = − Λ2 + (2i − 1) 2K
.
For a given input symbol vector D0 , the question arises that in what rate should the symbols be
placed in nonlinear spectral domain. In other words, how should we choose the nonlinear spectral
width of the signal, Λ, and if there is any optimal value for it? Note that this question originally
arises since we are signaling in the nonlinear spectral domain, which does not have a linear
duality with the time domain unlike conventional OFDM systems. To answer this question, the
effects of chromatic dispersion on data rate and performance of NFDM systems are investigated
in section 3. Furthermore, due to the signal-dependency of the noise in this domain, conventional
signal processing techniques are no longer as efficient, and thus, new techniques are studied in
section 4.
3.

Dispersion effects and the optimum value of nonlinear spectral width

By mapping the data only on CS, which is the focus of the current paper, soliton formation is
suppressed. In the absence of DS (e.g. solitons), the mechanism of cancellation of dispersion
with nonlinearity no longer exists, and consequently excessive temporal broadening is expected.
Therefore, large guard bands are necessary to avoid mixing of neighboring signals [28–30]. This
reduces the effective data transmission rate in bits per second. In this section, the effects of
chromatic dispersion on NFDM communication systems are evaluated. Note that, throughout this
paper, "bandwidth" refers to the actual linear bandwidth of time domain signal Ql (T) rather than
the width of signal in the nonlinear spectrum CS. Furthermore, the temporal width and bandwidth
are calculated at the window, which consists 99.5% of signal energy. Using this definition, the
numerical error due to truncation of the time domain signal is negligible.
Chromatic dispersion is an important effect, which cannot be avoided in time domain in current
NFDM systems because, based on the inverse scattering theorem [33], the signal in time domain
needs to be separated from consecutive signals so that the NFT operation works properly and
linearizes the channel. The effect of dispersion is removed only after NFT operation (i.e., phase
shift in the CS as in (7)), and thus the broadening in time domain cannot be ignored. Therefore, it
is beneficial to know how the broadening can be minimized in order to maximize the effective

data rate. As an initial estimate on the dispersion effect, we consider the linear fiber channel as
in [34, Chapter 3.4], where broadening factor is derived for Gaussian pulses assuming a linear
fiber. It is shown that an optimum initial temporal pulse width (or pulse rate) exists that minimizes
the broadening effect of chromatic dispersion while the dispersion effect increases for a higher
bandwidth and fiber length. The same phenomena can be also observed for the NFDM system
when the data is modulated on the CS.
For a given input symbol vector D0 with fixed number of symbols K, the signal bandwidth is
directly related to the nonlinear spectral width Λ [26]. Similar to the argument about a single
Gaussian pulse in the linear regime, it is expected that, for a fixed K, Λ (or, in other words, the
symbol rate) can be adjusted to minimize the received temporal width of signal. The existence of
an optimum nonlinear spectral width Λ in NFDM systems can be shown through the asymptotic
analysis presented in [26]. Assuming the data is modulated on CS with the nonlinear spectral
width Λ, the initial temporal width of the signal in time domain (after INFT and at L = 0) can
g0 = T0 πK/Λ for low-amplitude signals by approximating INFT as a linear
be estimated as ∆T
Fourier transform. When the amplitude of the signal increases this temporal width increases and
cannot be accurately determined using the linear estimation any more. Defining the parameter r
as the ratio between the actual temporal width after INFT compared with the linear estimation
g0 ), the actual temporal width can be expressed as ∆T0 = rT0 πK/Λ, where
above (i.e., r = ∆T0 /∆T
r has been shown to be increasing as a function of the signal power in nonlinear spectral domain
and independent of Λ [26]. In asymptoticly long fibers, the received temporal width at length L
denoted by ∆TL can be expressed by [26]
rT0 πK
+ 4zT0 Λ,
(10)
Λ
where z is defined in (4).
p It can be readily seen from (10) that the received temporal width ∆TL can
be minimized at Λ = rπK/4z. Consequently, for a given sequence of K CS symbols containing
fixed amount of information, an optimal Λ (or equivalently symbol rate) can be selected to
minimize the required guard interval leading to higher packet rates and effectively higher data
rates in bits per second.
It can be also confirmed numerically that an optimal value exists for spectral width Λ (or
bandwidth BW) that minimizes ∆TL in practical (non-asymptotic) scenarios. Assuming that
K = 128 random QPSK modulated symbols are mapped on CS using sinc pulses as in (9), Fig. 2
shows the existence of the optimal bandwidth for different fiber lengths. The optimal bandwidth
decreases for longer fibers. Furthermore, it is confirmed by simulation that the initial width
(i.e., ∆TL at L = 0) decreases when the bandwidth increases, and r is constant for all values of
spectral width Λ. Signal power in nonlinear spectral domain used in this simulation corresponds
to r = 1.99. It is worth mentioning that, for instance, the optimal value of Λ for L = 2000 km and
r = 1.99, calculated based on the asymptotic formula (10), is obtained as 9.76, which is close to
the actual value obtained by simulation in Fig. 2.
It should be noted that, as stated in [26, Table IV], by increasing K, the effect of dispersion
can be reduced and higher overall data rates are achieved even when Λ remains fixed rather
than being adjusted to the optimal value for the new K. On the other hand, the computational
complexity increases for larger K, and thus K and Λ (or bandwidth) should be optimized based
on the trade-offs between performance metrics such as computational complexity and effective
data rate. Such optimization problems will be the subject of future research.
Another crucial consequence of minimizing the received temporal width by choosing the
optimal value for Λ is the noise reduction in nonlinear spectral domain, which may lead to higher
achievable data rates by reducing the number of errors. In NFDM systems, each single sample of
data at a nonlinear frequency λ receives signal and noise contributions from all the components
within the received temporal width of signal. Therefore, as stated in [25], a larger received
temporal width results in higher noise in nonlinear spectral domain. This is similar to ordinary
∆TL =

BW [GHz]
10

50

20

30

40

L = 2000 km
L = 1000 km
L = 3000 km
L = 0 km

40

∆TL [ns]

50

30
20
10
0

2

4

6

8

10

12

14

16

18

Λ

Fig. 2. The temporal width of the propagated signal after distance L for different BW
and Λ.

linear communication systems, in which the noise power in time domain is directly related to
the signal bandwidth. In effect, if the received temporal width, over which the NFT operation is
performed, shrinks, the effective noise added to CS in nonlinear spectral domain decreases.
Here, a simulation is performed for different values of Λ = 4, 8, 12, which shows that the largest
received temporal width (corresponding to the smallest Λ = 4) results in higher noise in nonlinear
spectral domain. Figure 3 illustrates this result by plotting the noise variance for different values
of Λ against the amplitude of the signal, which also shows the signal-dependency of the noise. The
results in Fig. 3 demonstrate that the noise induced in nonlinear spectral domain is directly related
to the received temporal width at any amplitude of CS. This figure is obtained by calculating the
noise variance for more than 2 × 106 symbols (2 × 107 samples) in different nonlinear frequencies.
Figure 3 suggests that choosing the optimum Λ not only minimizes the dispersion effect, but also
reduces the noise. Therefore, improved error performance and achievable data rate is expected.
0.5
Λ=4, BW=13.2 GHz, ∆TL=23.6 ns
Λ=8, BW=26.3 GHz, ∆TL=16 ns

E{|η|2}

0.4

Λ=12, BW=39.6 GHz, ∆TL=17.3 ns

0.3

0.2

0.1

0

0

0.2

0.4

0.6

0.8

1

1.2

1.4

1.6

|ρ0|

Fig. 3. The variance of the noise in nonlinear spectral domain for different Λ

Further reduction in noise variance of the nonlinear spectral domain is expected if a dispersion
precompensation method, similar to the technique proposed in [32], is used. In [32], the signal
2
is first phase shifted by e+2jλ z in CS. Then, it is transformed by INFT to time domain and
transmitted through the fiber. This reduces the effect of dispersion up to two times. Consequently,
the received temporal width and the required guard interval are decreased. Also, it is possible to
2
introduce a complete precompensation by performing INFT on ρ0 (λ)e+4jλ z at the transmitter
side. It is expected that this method further reduces the noise in nonlinear spectral domain,
but the required guard interval would remain equivalent to the case in which no dispersion
precompensation is used. Due to the added complexity, none of these precompensation techniques
are considered in the current paper.
4.

Signaling methods

In this section, three different methods are introduced for performance enhancement of the
continuous spectrum channel (8) [6, 10]. As a benchmark we also consider the conventional
method of “direct mapping on CS”, which refers to the block diagram shown in Fig. 1 and
channel model (8). In order to properly understand the effect of signal-dependent noise, a ring
constellation is used for modulation. It should be noted that even in conventional fiber optic
systems ring constellation is a popular choice [2, 35, 36]. Unless otherwise stated, K = 128
random symbols are mapped on CS using sinc pulses as in (9) for Λ = 8. The fiber length is
L = 2000 km and the signal bandwidth is 26 GHz.
4.1.

Nonuniform signaling

Using nonuniform levels for constellation diagram is an effective method for improving the
performance of a channel with signal-dependent noise [35, 37]. If the statistics are known, the
optimum levels can be found by solving an optimization problem such as minimization of BER.
For the CS channel model (8), unfortunately, a closed-form description of the noise statistics is
not known. Therefore, deriving and solving an optimization problem for this channel would be
cumbersome. Here, we propose another approach based on the variance normalizing transform
(VNT) as introduced in [38, 39] for optimizing the constellation levels for NFDM systems.
The VNT was used in [26] as a tool to approximate the capacity of CS channel (8). It has
also been used for efficient signaling and capacity approximation over shot-noise-limited optical
channels [40–42]. In general, the variance normalizing transform can be applied to a random
variable, where its variance varies as a function of its mean, to generate an approximately
Gaussian random variable with a variance independent of its mean. Considering (8) and removing
the dependence to λ for simplicity, we have random variable ρ̃ L , with mean ρ0 and variance
E{| ρ̃ L − ρ0 | 2 } = f (ρ0 ). VNT that normalizes this random variable is defined as
∫
1
T(u) =
du.
(11)
p
f (u)
The normalized (i.e., transformed) random variable Y = T( ρ̃ L ) has then the statistics of σY2 ' 1
and µY ' T(ρ0 ) for sufficiently large values of ρ0 . In [43], it was also shown that the probability
distribution of the normalized random variable tends to Gaussian distribution for a family of
originally non-Gaussian probability distribution. Moreover, it was demonstrated in [26] by
numerical simulation that the distribution of noise on CS becomes Gaussian after applying VNT.
Applying VNT as defined in (11) to the noisy signal ρ̃ L = ρ0 + η L , we will have
Y = T( ρ̃ L ) = T(ρ0 + η L ) ' T(ρ0 ) + ηT ,

(12)

where ηT can be well approximated as a zero mean Gaussian noise with unit variance independent
of the transformed signal T(ρ0 ). Consequently, a communication channel can be defined, in

which the signal is originally generated in a transformed domain and then mapped into the
original signal domain ρ0 using inverse VNT. After transmission through the channel (8) and
addition of signal-dependent noise η L , VNT is applied, and the output signal can be expressed as
Y = T(ρ0 )+ηT , which defines a conventional additive Gaussian noise channel. As in [26, Lemma],
the capacity of this transformed channel is equal to the capacity of the original channel (8) with
signal-dependent noise. Moreover, the optimal signaling techniques for the conventional additive
white Gaussian noise (AWGN) channels would be also optimum for the transformed channel (12).
The analysis above is valid for real signaling on CS and can be extended to complex signaling
only if the real and imaginary channels of CS are independent. However, as mentioned in section
2.3, when signal is modulated on CS, the real and imaginary channels are not independent [26].
Moreover, although the noise on the phase of the signal is not dependent to its mean (i.e., input
phase), but it depends on the input amplitude. Therefore, VNT with its form in (11) can not
be applied to the complex plane. However, if only real channel is used (={ρ0 } = 0), the noise
variance would only depend on the mean (i.e., input signal), and VNT would be applicable. As
in [26], the signaling can be performed in the transformed channel, but this requires the forward
and inverse VNT to be performed each time a signal is transmitted and received. Alternatively,
the nonuniform signal levels defined by VNT and the corresponding decision boundaries required
for detection at the receiver can be determined for a known channel at the beginning of the
communication and then the determined optimum levels can be used for mapping data on CS.
In this paper, first it is assumed that ={ρ0 } = 0, and uniform levels are considered for
T(ρ0 ). Then, inverse VNT is applied to obtain the nonuniform levels for ρ0 . Then, the complex
modulation format is defined on the rings with radius equal to the levels found. This makes our
approach a sub-optimal method for complex signaling on CS because of the dependence of real
and imaginary channels. Note that the number of constellation points on each ring can be chosen
based on the variance of the phase of the noise and the required performance. The variance
normalizing transform for the 2000 km link is shown in Fig. 4(a), where it is observed that T(u)
converges to 16.32 at u → ∞ implying that a peak amplitude constraint exists for the transformed
variable T(ρ0 ). In other words, the signal space is limited in the transformed domain because
the variance of noise depends on the signal amplitude by a polynomial of order larger than
2 [26, Theorem]. For demonstration purposes, four uniform levels for T(ρ0 ) are chosen between
0 and 16.1, and nonuniform levels of ρ0 are calculated using inverse VNT. The corresponding
rings and decision boundaries are shown in Fig. 4(b). It is expected that if nonuniform levels in
Fig. 4(b) are used, the error rate will be lower compared to uniform levels with the same average
power. This will be demonstrated in section 5.1.
6

20

4
10

T(u)

2
0

0
-2

-10

-4
-20
-10

-5

0

u

5

10

-6
-6

-4

(a)

-2

0

2

4

6

(b)

Fig. 4. (a) Variance normalizing transform for 2000 km (b) Nonuniform levels (solid)
and decision boundaries (dotted) derived from VNT.

4.2.

Direct mapping on CS and filtering

The second method, which can potentially improve the performance, is the application of
linear filtering at the receiver in the original signal space, i.e., continuous spectrum. This
method is similar to match filtering when detecting pulse modulation in the nonlinear spectral
domain [10, 27] and requires oversampling of ρ̃ L (λ) at the receiver. If INFT is applied to ρ̃ L (λ)
after removing the dispersion induced phase shift according to equation (7), the signal is squeezed
into the temporal width of the unpropagated signal q0 (t) because of dispersion compensation in
CS and while noise is still present over the uncompensated temporal width of qL (t). To get back
to the signal space, we can then apply an NFT to the time domain signal after windowing the
signal component distributed only within the unpropagated temporal width thereby removing
the excess noise. However, this method requires complex numerical computation because of
the extra NFT/INFT operations on the oversampled signal, and thus, we use ordinary Fourier
transform instead of NFT, which allows for linear filtering of the CS signal. In other words, linear
inverse Fast Fourier transform (IFFT) is applied to ρ̃ L (λ), the noise out of the signal’s window is
removed, and then forward linear Fast Fourier transform (FFT) is applied. It should be noted that
I/FFT should not be almost equivalent to I/NFT (i.e., asymptotic linear regime) to be effective in
noise cancellation, as significant amount of noise is observed outside the signal interval even
after IFT. In fact, employing FFT and NFT would result in different signal (and noise) power
distributions. For example, as explained in section 3, so-called low-amplitude long “tails” [8, 21]
typically appear after INFT while compact (finite duration) and much more uniform signal power
distribution is observed after linear Fourier transform operation. Therefore, the employment of
a simple constant-gain windowing operation can more effectively mitigate the noise effect for
linear filtering compared to the NFT-based filtering which leads to long tails.
Assume H(t) as the ideal filter with width equal to the width of IFFT{ρ0 (λ)}, then the output
can be described as
ρ̂ L (λ) = FFT{H(t)IFFT{ ρ̃ L (λ)}}.
(13)
For the simulation, an oversampled version of ρ̃(λ) by a factor of 20 is first generated using
NFT, and then ρ̂ L (λ) is obtained from equation (13). The probability distribution of noise are
demonstrated in Fig. 5 before and after filtering for signal amplitudes of ρ0 (λ) = 1, 2. It can
be seen that the variance of noise after filtering (Figs. 5(c) and 5(d)) is substantially reduced
compared to the original noise (Figs. 5(a) and 5(b)). Note that although the signal-dependency of
noise is not eliminated by filtering, its effect is reduced considerably.
0.06

0.06

0.06

0.06

0.04

0.04

0.04

0.04

0.02

0.02

0.02

0.02

0
-2

0

(a)

2

0
-2

0

(b)

2

0
-2

0

(c)

2

0
-2

0

(d)

Fig. 5. The histogram of noise on the CS signal (a)-(b) ρ̃ L (λ) for amplitudes ρ0 (λ) = 1, 2,
and (c)-(d) ρ̂ L (λ) (i.e., filtered signal) for amplitudes ρ0 (λ) = 1, 2.

2

4.3.

GLM-based signaling

The signal modulated on the continuous spectrum ρ(λ) is related to the time domain signal q(t)
by the Gelfand-Levitan-Marchenko (GLM) equation [33]
K(x, y) +

∫∞ ∫∞
x

K(x, s 0)F ∗ (s + y)F(s 0 + s)dsds 0 = F ∗ (x + y),

(14)

x

where

1
F(t) =
2π

∫

∞

ρ(λ)e jλt dλ

(15)

−∞

is the inverse Fourier transform (IFT) of ρ(λ). One can interpret that the GLM equation performs
a nonlinear mapping on the linear inverse Fourier transform of ρ(λ) (i.e. F(t)), so that the time
domain signal q(t) = −2K(t, t) only contains CS and no DS. This means that if GLM is applied to
any arbitrary input signal, the output signal would not include any soliton components. Therefore,
GLM can be used as a transformation, based on which a signaling method can be designed for
the CS channel.
In this section, a GLM-based signaling approach is proposed as shown in Fig. 6 [45]. The main
difference of the system in Fig. 6 compared to the direct mapping on CS (as in Fig. 1) is that the
inverse Fourier transform operation included in INFT operation is transferred from transmitter to
the receiver, and, effectively, data is mapped on F(t) (i.e., linear Fourier transform of CS) defined
as
K
Õ
K
K 2i − 1
F0 (t) =
D0i sinc( t + −
),
(16)
τ
2
2
i=1
where τ is the corresponding signal window in time domain determined based on system
constraints, such as bandwidth and power. Consequently, an IFT block is needed at the receiver
after NFT operation to return to the original signaling space, where the noisy version of (16) is
sampled at the center of pulses and are decoded to the original data. It should be also noted that
this GLM-based signaling using sinc pulses is in fact equivalent to the nonlinear inverse synthesis
method with Nyquist pulse shaping studied in [28–31]. However, either in (9) or (16), using sinc
pulse shapes are not necessary, and more efficient pulse shapes may be found by optimizing
the trade-off between performance and numerical error. A similar method was also proposed
independently in [46], in which data are mapped on the kernel of GLM equations for discrete and
continuous spectrum.
Using the law of large numbers, we can show that applying (inverse) Fourier transform on the
received CS diminishes the signal-dependency of the noise [40] while the probability distribution
of the transformed symbols tends to Gaussian according to the central limit theorem [44].
Therefore, conventional modulation and signal processing techniques should perform efficiently
if applied in the system illustrated in Fig. 6.

D0

Sampling

FL(t)

IFT

GLM

ρL(λ) Equalizer ρL(λ)
eq. (7)

q0(t)

NFT

Tx

qL(t)

Q0(T)
SSMF

DL

F0(t)

Mapping
on F(t)

Rx

QL(T)

Fig. 6. The block diagram of the GLM-based signaling system.

As discussed in Section 2.3, the variance of noise sample on CS, η L (λ), depends on the input
sample ρ0 (λ) = FT {F0 (t)}. Let κ denote the number of samples in nonlinear spectral domain,
h = Λ/κ be the distance between subsequent samples, and λk be the kth sample point. The noise
for GLM-based signaling is defined as ΓL (t) = FL (t) − F0 (t). Therefore, for a given input signal
F0 (t) = IFT{ρ0 (λ)} and by using the discretized version of (15), we have
 2 Õ
κ Õ
κ
h
2
E{<[ΓL (t)] } ≈
E{<[η L (λk )]<[η L (λ p )]} cos(λk t) cos(λ p t)
(17)
2π k=1 p=1
 2 κ κ
1 h ÕÕ
δ(k − p) f [| ρ0 (λk )|] cos(λ p t) cos(λk t)
(18)
=
2 2π k=1 p=1
 2 κ
1 h Õ
=
f [| ρ0 (λk )|] cos(λk t)2
(19)
2 2π k=1
κ

=

κ

hΛ 1 Õ
hΛ 1 Õ
f
[|
ρ
(λ
)|
+
f [| ρ0 (λk )|] cos(2λk t)
0
k
16π 2 κ k=1
16π 2 κ k=1

(20)

hΛ
E{ f [| ρ0 |]}.
(21)
16π 2
For the last equality, the law of large numbers can be used (at large κ) so that the first sum
term in (20) can be described as the expected value of f [| ρ0 |], and the second term tends to zero
because of the presence of the cosine term. A similar result can be demonstrated for =[ΓL (t)],
and it is observed from (21) that the variance of <[ΓL (t)] and =[ΓL (t)] at every sampling time
of t are independent of signal F0 (t) at that time. However, it is clear that the variance of noise
<[ΓL (t)] and =[ΓL (t)] depend on average value of f [| ρ0 |], indicating that higher energy (i.e.,
higher average power) results in larger noise variance. This will be investigated in terms of error
rate in Section 5.1. Also, according to the central limit theorem [44], the probability distribution
of ΓL (t) tends to Gaussian. This is similar to the effect observed in visible light communications
where the non-Gaussian clipping noise tends to a Gaussian one after Fourier transformation in
optical OFDM receivers [47].
Here, a simulation is performed in order to compare the statistics of the noise added to ρ0 (λ)
(i.e., direct mapping) and F0 (t) (i.e, GLM-based signaling) after fiber propagation and addition
of ASE noise. NFT is applied to the received signal at the output of the fiber, and also an
additional IFT operation is performed to obtain the noise on the signal F0 (t). The same simulation
parameters as previous sections are used here, and statistics of η L (λ) are compared with that of
ΓL (t). Note that oversampling (κ > K) leads to smaller noise variance as a result of averaging
effect of Fourier transform operation [40]. From (21), it can be seen that smaller h (larger number
of samples κ) results in noise reduction. However, in practical scenarios, noise samples become
correlated if they are so close to each other as shown in [26, Section V.D]. Thus, reducing h
is effective up to the point, for which the noise samples are correlated. In this simulation and
later in this paper h = 0.01 is chosen based on observations of noise correlation. In order to
investigate the effect of signal-dependency two additional parameters are defined, which represent
signal-to-noise ratio (SNR) associated with ρ(λ) and F(t) at a specific amplitude of x.
≈

SNRρ (|x|) =

|x| 2
,
E{|η L | 2 }| |ρ0 |=x

(22)

SNRF (|x|) =

|x| 2
.
E{|ΓL | 2 }| |F0 |=x

(23)

Figures 7(a) and 7(b) show the values for SNR parameters defined above. It can be seen clearly
in Figs. 7(a) and 7(b) that SNRF is increasing while SNRρ has a maximum. Also, the values

of SNRF are higher than SNRρ . Thus, better error performance is expected if data is mapped
on F0 (t). Furthermore, the probability distribution for ΓL (t) has a Gaussian-like distribution
with constant variance for all signal amplitudes, as shown in Fig. 7(d), while the distribution
of η L (λ) is non-Gaussian, as depicted in Fig. 7(c). The signal-independency of noise ΓL (t) will
be also shown numerically later in Fig. 8(d) of Section 5.1. We can therefore conclude that if
data is mapped on F0 (t) (i.e. GLM-based signaling), it is distorted by an approximately additive
signal-independent Gaussian noise as in the conventional AWGN channels where the underlying
uniform amplitude signaling remains efficient. On the other hand, in the direct mapping scheme,
the uniform distribution of the symbol amplitudes causes low SNR at higher amplitudes rendering
it essentially inefficient at high energy levels.
0.06

5

0.04

0.04

0

0.02

0.02

2

log(SNRρ)

10

log(SNRF)

0.06

3

1
0

-1

0

1 |ρ | 2

3

-5

0

0

(a)

0.2 |F | 0.4

0.6

0
-2

0

2

0

-0.02

0

0.02

0

(b)

(c)

(d)

Fig. 7. (a) SNRρ for different values of | ρ0 (λ)|. (b) SNRF for different values of |F0 (t)|.
(c) Distribution of noise η L (λ) for | ρ0 (λ)| = 2. (d) Distribution of noise ΓL (t) for all
values of |F0 (t)|.

5.

Numerical results and discussions

In this section, we present numerical results on the performance of different signaling techniques
discussed earlier. Two-ring constellations are used for modulation in order to demonstrate the
effect of dependence of noise to the signal amplitude. Unless otherwise stated, K = 128 symbols
are mapped on the nonlinear spectral width of Λ = 8, which are chosen randomly based on the
two-ring constellation. First, the methods described in previous section are compared regarding
symbol error rate (SER) only taking into account the error due to noise on the amplitude of two
particular symbols. Then, bit error rate (BER) performance is investigated based on different
constellation formats on two rings. The performance results are presented against the energy of
the CS signal denoted by E. The channel bandwidth is 26 GHz.
5.1.

Error rate performance

In the first simulation, SER is compared for different signaling methods assuming only two
symbols on the two rings with an equal phase of zero. This allows us to investigate the effect
of signal-dependency of noise. Figure 8 shows the received noisy symbols at the end of a
L = 2000 km fiber for different signaling methods with E = 2. Figure 8(a) shows the result for
the benchmark scheme, where the symbols are directly mapped on CS as explained in section 2.3.
It is observed that the variance of both real and imaginary parts (or amplitude and phase) of the
noise is dependent on the amplitude of the original transmitted symbol. Figure 8(b) demonstrates
the received symbol for the same energy but for nonuniform levels determined by VNT. In Fig.
8(a), ring radii are 0.8 and 1.6, and the midpoint decision boundary is at 1.2. However, in Fig.
8(b), ring radii are 0.575 and 1.685 with the decision boundary 0.975. Note that the statistics of
the noise in CS does not change by applying the VNT method, but the performance improvement
is achieved as a results of choosing optimal levels and decision boundary for the rings. For the
next two signaling methods in Figs. 8(c) and 8(d), the magnitude of both amplitude and phase

of the noise changes. As shown in Fig. 8(c), when the linear filtering is applied, as described
in section 4.2, the noise is significantly reduced. For this simulation, an oversampling of 20
CS samples per symbol is considered. Ring radii and decision boundaries are the same as Fig.
8(a). It should be noted that the the noise is still signal-dependent although significantly reduced.
Figure 8(d) depicts, the received symbols when GLM-based signaling is employed. The levels
are determined so that the energy E is the same as previous methods. It can be observed that the
noise statistics are almost uniform for two levels as predicted in section 4.3. Comparing Figs. 8(a)
and 8(d), it is observed how crucial it is to whether directly map digital data on CS or on the
GLM kernel as the signal would experience two completely different channels.
1

1

1

0.1

0.5

0.5

0.5

0.05

0

0

0

0

-0.5

-0.5

-0.5

-0.05

-1

0

1

2

3

-1

0

1

(a)

2

-1

3

0

1

(b)

2

3

-0.1

0

(c)

0.1

0.2

0.3

(d)

Fig. 8. Received noisy symbols at L = 2000 km based on (a) Direct mapping on CS, (b)
Nonuniform signaling, (c) Mapping on CS and filtering, and (d) GLM-based signaling.

Symbol error rate for two symbols considered in Fig. 8 are calculated and shown in Fig. 9 for
different signal energies. It is observed that all proposed methods result in improved performance
compared to the direct mapping on CS. The minimum SER is achieved at similar energies for
three out of the four methods, but for VNT-based nonuniform signaling the optimum point is
different. We conjecture that after E = 3.7 the VNT is not effective because the noise variance is
very large, and the large mean condition for the VNT is not satisfied. For the filtering method,
despite noise reduction, the signal-dependency of noise is not eliminated, and thus performance
is degraded after a specific energy. For the GLM-based signaling, although signal-dependency is
eliminated for a constant energy, the noise variance increases for signals with higher energies
based on (21), where higher error rate is expected.
100

Direct mapping on CS
Nonunifrom signaling
Mapping on CS and filtering
GLM-based signaling

SER

10-1

10-2

10-3

10-4

0

1

2

E

3

4

5

Fig. 9. SER for different signaling methods on the CS channel.
In Fig. 10, the BER performance for three different constellation designs based on the two-ring
pattern is depicted. In general, lower BER are achieved compared to direct mapping on CS,

when each of the proposed methods is employed, and therefore, higher achievable data rates are
expected. However, since NFT-based transmission only eliminates the signal-signal nonlinear
interactions, it is expected that signal-noise and noise-noise nonlinear interactions in time domain
would eventually limit the performance. Soft decision forward error correction (SDFEC) threshold
of 2 × 10−2 [16] is also shown in figures. For constellation (I) in Fig. 10(b), the minimum BER
is achieved with similar values for the filtering method and GLM-based signaling. Since the
symbols are well separated in each ring the behavior of BER for different methods is similar to
SER in Fig. 9. Nevertheless, when the outer ring is π/4 shifted in phase all the methods results
in better BER except the nonuniform signaling based on VNT. This is because the decision
boundary levels for nonuniform signaling are defined based on the one dimensional (i.e., real)
signal space considered in Section 4.1, which are not optimum for the complex constellation (II)
in Fig. 10(c). An optimal set of decision boundaries could be only determined based on maximum
likelihood requiring full analytical description of the noise statistics, which are not available or
through exhaustive numerical search. Note that, unlike constellation (II), the decision boundary
levels defined based on one dimensional signal space are nearly optimal for constellation (I) since
errors dominantly occur in the radial directions.
10-1

Direct mapping on CS
Nonunifrom signaling
Mapping on CS and filtering
GLM-based signaling

10

BER

(II)

(I)

SDFEC
-2

(III)

10-3

10-4

0

1

2

(a)

3

4

5

3

4

5

E

(b)
10-1

10-1
SDFEC
SDFEC

10-2

10-3

BER

BER

10-2

10-4

10-3

-5

10

10-6

0

1

2

3

4

5

10-4

0

1

2

E

E

(c)

(d)

Fig. 10. (a) Figure legends and constellation diagrams. (b)-(d) BER for different
signaling methods at L = 2000 km for constellation (I)-(III) respectively.

In Fig. 10(d), 8 symbols are chosen on each ring, and BER is obtained by simulation. For
constellation (III), error originated from the noise on the phase is dominant, which causes the
BER for nonuniform signaling not to be as low as other proposed methods. For the GLM-based
signaling, the slope of the performance degradation is higher compared to all other methods.

This can be explained by referring to equation (21). Based on this equation the noise variance
increases for higher signal energies with a polynomial of order 4, and this happens equally for all
symbol amplitudes. However, for other methods, the inner ring always experiences lower error
rate compared to the outer ring resulting in a slower increase of BER on average. It should be
noted that, as explained in Section 3, the method proposed in [32] can be used to further improve
the performance.
Note that, in Fig. 10, operating bandwidth is constant and BER is obtained for various energies,
but the effective data rates vary for different signaling methods and constellations. In fact, the
required temporal window for NFT operation without numerical error is different for each signal
energy E and signaling methods. This happens because of the effect of residual tail at the output
of INFT and the different effect of dispersion. For instance, at E = 2, taking into account the
dispersion effect, for Figs. 10(b) and 10(c) the effective bit rate is 21.9 Gpbs for first three methods
and 20.7 Gbps for GLM-based signaling, while for constellation in Fig. 10(d) these values are
respectively 29.2 and 27.6 Gpbs.
The corresponding launch power levels in Figs. 9 and 10 range from -17 dBm to -10dBm.
Not that, due to optimization performed for K = 128 in Section 3, the dispersion effect is
significantly reduced and in turn the nonlinear effect is intensified such that the CS channel with
power levels above can reach the high nonlinearity regimes where the BER performances start
declining. Assuming a larger K, the optimal bandwidth increases, which leads to higher launch
power; however, studying such a scenario would also require dealing with significantly more
time-consuming simulation of NFT and INFT operations.
It is worth mentioning that better performance for techniques presented here are obtained with
a cost of computational complexity. In order to compare the complexity of different methods, the
number of NFT operations at the receiver are taken into account. For direct mapping on CS and
nonuniform signaling, the number of NFT operations is equal to the number of symbols, which,
in our case, is K = 128. For GLM-based signaling, 6 × K samples of ρ̃ L (λ) (corresponding
to h = 0.01) are computed (768 NFT operations) for calculation of FL (t). Also for the linear
filtering method, as stated earlier, 20 × K samples of ρ̃ L (λ) are computed (2560 NFT operations),
and also additional FFT and IFFT operations are performed.
5.2.

Distance reach

Considering soft and hard decision FEC thresholds the distance reach are compared for different
methods in Fig. 11. For all methods shifted constellation (II) in Fig. 10 is used except for the
nonuniform signaling based on VNT, for which constellation (I) is selected. The input energy E
is chosen for the minimum achieved BER for each technique. Considering SDFEC and ignoring
the redundant bits due to FEC coding, the simulation results demonstrate that distances as long as
7100 km can be reached by applying the linear filtering technique. The effective bit rate at distance
7100 km is 9.6 Gbps. At second place, with GLM-based signaling, a distance reach of about
5900 km with effective bit rate 11.7 Gpbs is obtained. When the levels are optimized using VNT,
the achieved reach and effective bit rate are 5000 km and 12.7 Gbps. The corresponding values
are respectively 3000 km and 18.5 Gpbs for direct mapping on CS. Moreover, considering hard
decision forward error correction code (HDFEC) with 7% overhead with threshold of 3.8 × 10−3 ,
distance reaches of 1900, 3300, 4300, and 5300 are achieved respectively for direct mapping
on CS, nonuniform signaling, GLM-based signaling, and filtering. It should be noted that, for
simplicity and better presentation, E and Λ are considered to be fixed for all fiber lengths and
are chosen based on the best performance of a 2000 km link. However, for further performance
improvement, E and Λ should be optimized for each individual fiber length. This analysis requires
extensive numerical simulation and is out of scope of the current paper as marginal improvement
is expected.

10-1
SDFEC

BER

10-2

HDFEC

10-3
10-4

Direct mapping on CS
Nonuniform signaling
Mapping on CS and filtering
GLM-based signaling

10-5
10-6
1000

2000

3000

4000

5000

6000

7000

8000

L [km]

Fig. 11. Distance reach for different methods.
6.

Summary

Efficient signaling on CS was studied in this paper. First, the effect of chromatic dispersion on
NFDM systems was studied. It was demonstrated that the signal bandwidth can be chosen so
that the temporal signal broadening at the end of the fiber is minimized. By correctly selecting
the bandwidth, not only the data rate is maximized, but also the noise variance is minimized in
nonlinear spectral domain. Moreover, three different techniques were investigated for improving
error performance compared to direct mapping on CS. For instance, our analysis showed that
applying a linear filter on CS at the receiver significantly reduces the noise and improves the BER
performance so that 9.6 Gbps can be transfered over 7100 km of fiber by only mapping the data
on CS. The signal bandwidth was considered 26 GHz in our numerical results, however, data
rates beyond values presented here may potentially be achieved by exploiting higher bandwidth
and mapping data on DS as well.
Acknowledgment
The authors would like to thank Professor Frank R. Kschischang for insightful discussions on
simulation and signaling of the continuous spectrum channel.

