Experimental Demonstration of Dual Polarization Nonlinear
Frequency Division Multiplexed Optical Transmission System
S. Gaiarin(1) , A. M. Perego(2) , E. P. da Silva(1) , F. Da Ros(1) , and D. Zibar(1)
(1)
(2)

DTU Fotonik, Technical University of Denmark, Lyngby, 2800, Denmark simga@fotonik.dtu.dk
AIPT, Aston University, Birmingham B4 7E7 UK, peregoa@aston.ac.uk

Abstract Multi-eigenvalues transmission with information encoded simultaneously in both orthogonal
polarizations is experimentally demonstrated. Performance below the HD-FEC limit is demonstrated for
8-bits/symbol 1-GBd signals after transmission up to 207 km of SSMF.

arXiv:1708.00350v1 [] 1 Aug 2017

Introduction
The increasing demand for capacity and the limitations in achievable information rates imposed by
fiber nonlinearities to the conventional coherent
transmission systems have pushed part of the research community to investigate new techniques
to encode and transmit information over the fiber
channel. Starting from the inverse scattering
transform theory, the concept of nonlinear frequency division multiplexed (NFDM) system has
been introduced few years ago 1 .The first investigations focused on modulating only the position of
a few discrete eigenvalues and transmitting over
short distances. Further extensive research enabled to demonstrate NFDM systems using multiple degrees of freedom (eigenvalue, spectral amplitude) transmitting through thousands of kilometers 2,3 .
Nevertheless, NFDM systems still suffer from
many challenges which have not yet allowed to
match the spectral efficiencies achievable by the
technologically mature linear coherent transmission systems. One of the main limitations is
the inability of taking advantage of both polarizations supported by the standard single mode fiber
(SMF), which can potentially double the spectral
efficiency of such systems. A first preliminary attempt in this direction has been reported by the
numerical analysis of 4 .
In this work, exploiting a rigorous mathematical technique, we report, to the best of our
knowledge, the first experimental demonstration
of NFDM transmission independently modulating
the two polarizations. The information data is encoded in the phases of the four nonlinear Fourier
(NF) coefficients b(λ) corresponding to two discrete eigenvalues, using quadrature phase shift
keying (QPSK). The system is able to reach a
transmission distance of 207 km and 166 km at
the hard-decision forward error correction (HDFEC, BER=3.8×10−3 ) threshold, over a standard
link with EDFA amplification of 41.5 km and 83 km
SMF fiber spans, respectively.

Eigenvalue
Im{λ}

Pol. X
Im{b1(λ)}

Pol. Y
Im{b2(λ)}

0.6j

0.3j

Re{b1(λ)}

Re{b2(λ)}

Re{λ}

Fig. 1: Ideal eigenvalues and corresponding spectral
amplitude constellations diagrams.

Mathematical framework
The signal propagation in single mode fibers
(SMFs) can be modeled by the so called Manakov system (MS) which describes the evolution
of the pair of complex field amplitudes q1 and q2
corresponding to the two polarizations of the field
envelope . The MS for the focusing regime in normalized form reads:
∂q1
=
∂z
∂q2
=
∂z


∂ 2 q1
+ i2q1 |q1 |2 + |q2 |2 ,
∂t2

∂ 2 q2
i 2 + i2q2 |q1 |2 + |q2 |2 ,
∂t
i

(1)

where z and t are the normalized space and time
coordinates respectively. The MS in Eqs. (1) is
integrable by the inverse scattering transform 5 ,
also called nonlinear Fourier transform (NFT). It is
therefore possible to associate a Lax pair of operators to the MS and define an eigenvalue problem
with eigenvalues λ and eigenvectors ψ. The time
evolution differential equation for ψ , induced by
the eigenvalue problem, is
  −iλ
ψ1

∂t ψ2  = −q1∗
ψ3
−q ∗
2

q1
iλ
0

q2

 
ψ1

0  ψ2  .
ψ3
iλ

(2)

By integrating Eq. (2) using appropriate boundary conditions, the eigenvalues λ and the
NF coefficients {a(λ), b1 (λ), b2 (λ)} can be obtained. Compared to the single-polarization case
{a(λ), b1 (λ)}), the MS has an additional NF coefficient, b2 (λ). This new degree of freedom can
be used in a NFDM communication system to en-

Eigen. 2

Transmission link

Coherent
frontend

Rx DSP

Q

SMF
EDFA 40km
FL

RF Amp.

BPD

DSO 80 GSa/s

SMF
EDFA 40km

I

BPD

PC
DP-IQ

BPD

Q

Pol. diverse hybrid

I

BPD

AWG 64 GS/s, 20 GHz

Signal predistortion

PRBS generation and
NFT mapping
INFT

Eigen. 1

Eigen. 2

Eigen. 1

Optical Carrier
IQ Modulation

Tx DSP

Fiber laser (Tx/Rx-LO)

Fig. 2: Experimental setup with transmitter and receiver DSP chain. Four constellations (two per polarization) associated to the
two eigenvalues are shown at the transmitter side (top left) and after 373.5 km transmission with 41.5-km spans (top right).

code twice as many information bits.
The time domain waveform associated to a
set of NF coefficients can be synthesized from
a null field by using the Darboux transformation
(DT) for the Manakov system. Such a powerful
technique allows to add iteratively new eigenvalues to the NF discrete spectrum by updating the
time domain signals correspondingly. The generic
solution ψ of (2) is transformed by the DT as
ψ̃ = (λ0 I3 − G0 )ψ, where G0 = ΨM0 Ψ−1 ,9 with


ψ1 ψ2∗
ψ3∗
0 .
Ψ = ψ2 −ψ1∗
(3)
ψ3
0
−ψ1∗
The ψj are the components of the eigenvalue
problem solution that are being transformed,
M0 = diag(λ0 , λ∗0 , λ∗0 ) and λ0 is the new eigenvalue we want to add to the spectrum. The signals are changed by the DT in the following way:
q̃j = qj + 2i(λ∗0 − λ0 )

1+

u∗j
P2

k=1

|uk |2

, i = 1, 2 (4)

being q̃j the modified component and
uj = ψj+1 /ψ1 , where j = 1, 2 indicates one
of the two polarizations.
Experimental setup
In order to demonstrate the advantage provided
by the new degree of freedom b2 (λ), the dualpolarization NFDM signaling depicted in Fig. 1 is
implemented by choosing two reference eigenvalues λ1 = 0.3i and λ2 = 0.6i.
Data symbols are generated at 1 GBd by mapping the bits from a pseudo random bit sequence
of length 211 − 1 to the phases of each of the
NFT coefficient pairs {b1 (λi ), b2 (λi )} according to
a QPSK constellation for b1,2 (λ1 ) and to a QPSK
constellation rotated by π/4 for b1,2 (λ2 ), (Fig 1).
The data so encoded are experimentally transmitted through the setup of Fig. 2. Inverse NFT

(INFT) is used to generate time waveforms from
the symbol {b1 (λ1 ), b2 (λ1 ), b1 (λ2 ), b2 (λ2 )} using
an iterative algorithm based on the DT (see previous section). The digital signal is pre-distorted
to compensate for the Mach-Zehnder modulator
nonlinear response and the waveforms are uploaded to an arbitrary waveform generator (AWG,
20-GHz analog bandwidth and 64-GSa/s) which
drives a dual polarization I/Q modulator. The optical input of the modulator is a fiber laser (FL) with
∼1 kHz linewidth, which is amplified after modulation and transmitted to a fiber link of multiple
spans of SMF with dispersion D = 17.5 ps/nm·km,
nonlinear coefficient γ = 1.25 W−1 km−1 and attenuation α = 0.195 dB/km.
At the receiver, the polarization demultiplexing
is performed manually through a polarization controller as no blind polarization demultiplexing algorithms for the NFT are yet available. The signal
is then detected in homodyne configuration by using a standard coherent receiver (33-GHz analog
bandwidth, 80 GSa/s). The acquired digital signal
is rescaled to the ideal power of the signal generated by the INFT and low-pass filtered to remove
the out of band noise. FFT-based coarse clock
recovery is used in order to align the pulse to the
center of the processing time window. The signal is then passed to the NFT block which first
locates the eigenvalues λi (zeros of the NFT coefficient a(λi )) and then computes b1,2 (λi ). In order
to minimize numerical errors, the NFT coefficients

Fig. 3: BER as a function of the OSNR in back-to-back.

10 -1
Pol. 1
Pol. 2

BER

10 -2
10

HD-FEC

HD-FEC

HD-FEC

-3

10 -4

40km span
80km span

Pol. 1
Pol. 2

40

120

200

280

360 80

240

400

560

720 40

200

360

520

680

Distance [Km]
Distance [Km]
Distance [Km]
Fig. 4: BER as a function of the transmission distance for each individual eigenvalue for 41.5 km (a) and 83 km (b) spans.
(c) Comparison of the average BER versus transmission distance between links of the two different span lengths.

are computed by integrating (2) using a generalized version of the forward-backward trapezoidal
algorithm 7 . The standard blind phase search algorithm is used independently on each of the constellations to perform carrier phase estimation.
Experimental results
The proposed system is first tested in a back-toback (BTB) configuration by removing the transmission link and varying the optical signal to noise
ratio (OSNR) at the receiver input through noise
loading. The measured signal bit error rate (BER)
is shown in Fig. 3. The BER averaged over the
four constellations saturates toward a value of
5×10−4 as the OSNR is increased. The main
contribution of errors comes from the two constellations of λ2 . As shown in Fig. 2 (top left) at
maximum OSNR, both constellations related to
λ2 are more noisy and distorted at the transmitter due to the higher peak-to-average power ratio of the λ2 -pulse portion, and the limited digital
to analog converter resolution. Additionally, imperfections in the digital signal processing (DSP)
chain at the receiver can introduce further degradation. For example, sub-optimum timing recovery due to noise, results in distortions of the received constellations, and an increment in BER.
Such a transceiver-induced penalty could be improved by more sophisticated DSP schemes.
The system is then tested by transmitting over
a link of SMF fibers spans, first with spans of
41.5 km and then with spans of 83 km. In
Fig. 4 (a) and (b) the BER curves for the four different b(λ) are shown as a function of the transmission distance for the 41.5 km and 83 km spans
case respectively. As in Fig. 3, the BER is sensibly higher for b1,2 (λ2 ) as already hinted by the received constellations after a 373.5-km transmission (Fig. 2). The constellations associated to λ2
are more degraded while those associated to λ1
are still well defined. The performance of b1 (λi )
and b2 (λi ), which are associated to the two polarizations, are instead similar for the same eigenvalue i. In Fig. 4 (c) we compare the average BER
for the two span lengths. Theoretically, worse performance may be expected for the longer spans
(further away from lossless transmission). The
results of Fig. 4 (c) instead, do not show a rele-

vant difference and the performance are similar
in the two cases at the same transmission distance. This is probably due to the performance
being mainly limited by transceivers impairments,
so that the losses play a minor role. Transmission
distances up to 207 km for the configuration using
41.5 km long spans and up to 166 km for the one
using spans of 83 km are achievable with a BER
below the HD-FEC threshold of 3.8×10−3 .
Conclusions
The results presented in the paper demonstrate
that polarization multiplexing in NFDM systems is
feasible and can lead to a potential increase of the
spectral efficiency achieved up to now. The performance of the two polarizations fields reported
are similar, while there is a clear difference in
BER between constellations associated to different eigenvalues.
Acknowledgements
This work is supported by the Marie Curie Actions
through ICONE Project (no. 608099). We thank M.
Kamalian for the stimulating discussions.

References
[1] Y. I. Mansoor, F. Kschischang, ”Information Transmission Using the Nonlinear Fourier Transform, Part
I:MathematicalTools,” IEEE Trans on Inf. Theory, Vol.
60(7), p. 4312-4328 (2014).
[2] S. K. Turitsyn et al. “Nonlinear Fourier transform for optical data processing and transmission: advances and perspectives,” Optica, Vol. 4(3), p. 307 (2017).
[3] V. Aref, S. T. Le, and H. Buelow, “Demonstration of Fully
Nonlinear Spectrum Modulated System in the Highly Nonlinear Optical Transmission Regime,” ECOC, Th.3.B.2,
Düsseldorf (2016).
[4] A. Maruta, Y. Matsuda, “Polarization division multiplexed
optical eigenvalue modulation”, PS, Florence, 2015
[5] S. V. Manakov, “On the theory of two-dimensional stationary self-focussing of electromagnetic waves”, Sov. Phys.JETP Vol. 38, no.2 p. 248 (1974).
[6] O. C. Wright, “The Darboux transformation of some Manakov systems,” Applied Mathematics Letters Vol.16, p. 647
(2003).
[7] V. Aref, “Control and detection of discrete spectral amplitudes in nonlinear fourier transform,” arXiv:1605.06328v1
(2016).

