1. Introduction
The vast majority of the global Internet traffic is conveyed through fiber-optical networks, which form the backbone of our information society. To cope with the growing data demand, the fiber-optical networks have evolved from regenerated direct-detection systems to coherent wavelength division multiplexing (WDM) ones. Newly emerging bandwidth-hungry services, like Internet-of-Things (IoT) applications and cloud processing, require even higher data rates. Motivated by this ever-growing demand, an increasing attention has been devoted in recent years to the analysis of the capacity of the fiber-optical channel.
Finding the capacity of the fiber-optical channel that is governed by the stochastic nonlinear Schrödinger (NLS) equation ([
1], Equation (
1)), which captures the effects of Kerr nonlinearity, chromatic dispersion, and amplification noise, remains an open problem. An information-theoretic analysis of the NLS channel is cumbersome because of a complicated signal–noise interaction caused by the interplay between the nonlinearity and the dispersion [
2]. In general, capacity analyses of optical fibers are performed either by considering simplified channels, or by evaluating mismatched decoding lower bounds [
3] via simulations (see [
4] and ([
5], Sec. I) for excellent literature reviews). Lower bounds based on the mismatch-decoding framework go to zero after reaching a maximum (see, for example, [
2,
6,
7,
8,
9]). Capacity lower bounds with a similar behavior are also reported in [
10]. In [
11], it has been shown that the maximum value of a capacity lower bound can be increased by increasing fiber dispersion, which mitigates the effects of nonlinearity. To establish a capacity upper bound, Kramer et al. [
12] used the split-step Fourier (SSF) method, which is a standard approach to solve the NLS equation numerically ([
13], Sec. 2.4.1), to derive a discrete-time channel model. They proved that the capacity of this discrete-time model is upper-bounded by that of an equivalent additive white Gaussian noise (AWGN) channel. In contrast to the available lower bounds, which fall to zero or saturate at high powers, this upper bound, which is the only one available for a realistic fiber channel model, grows unboundedly.
Since the information-theoretic analysis of the NLS channel is difficult, to approximate capacity, one can resort to simplified models, a number of which have been studied in the literature (see [
14] and references therein for a recent review). Two approaches to obtain such models are to use the regular perturbation or the logarithmic perturbation methods. In the former, the effects of nonlinearity are captured by an additive perturbative term [
15,
16]. This approach yields a discrete-time channel with input–output relation
([
14], Equation (
5)), where
and
are the transmitted and the received symbols, respectively;
is the amplification noise; and
is the perturbative nonlinear distortion. This model holds under the simplifying assumption that both the nonlinearity and the signal–noise interaction are weak, which is reasonable only at low power.
Regular perturbative fiber-optical channel models, with or without memory, have been extensively investigated in the literature. In [
17], a first-order perturbative model for WDM systems with arbitrary filtering and sampling demodulation, and coherent detection is proposed. The accuracy of the model is assessed by comparing the value of a mismatch-decoding lower bound, which is derived analytically based on the perturbative model, with simulation results over a realistic fiber-optical channel. A good agreement at all power levels is observed. The capacity of a perturbative multiple-access channel is studied in [
18]. It is shown that the nonlinear crosstalk between channels does not affect the capacity region when the information from all the channels is optimally used at each detector. However, if joint processing is not possible (it is typically computationally demanding [
19]), the channel capacity is limited by the inter-channel distortion.
Another class of simplified models, which are equivalent to the regular perturbative ones up to a first-order linearization, is that of logarithmic perturbative models, where the nonlinear distortion term
is modeled as a phase shift. This yields a discrete-time channel with input–output relation
([
14], Equation (
7)). In [
5], a single-span optical channel model for a two-user WDM transmission system is developed from a coupled NLS equation, neglecting the dispersion effects within the WDM bands. The channel model in [
5] resembles the perturbative logarithmic models. The authors study the capacity region of this channel in the high-power regime. It is shown that the capacity pre-log pair (1,1) is achievable, where the capacity pre-log is defined as the asymptotic limit of
for
, where
P is the input power and
is the capacity.
Despite the fact that the aforementioned simplified channels are valid in the low-power regime, these models are often used also in the moderate- and high-power regimes. Currently, it is unclear to what extent the simplifications used to obtain these models influence the capacity at high powers. To find out, we study the capacity of two single-channel memoryless perturbative models, namely, a regular perturbative channel (RPC), and a logarithmic perturbative channel (LPC). To assess accuracy of these two perturbative models, we also investigate the per-sample capacity of a memoryless NLS channel (MNC).
To enable an information-theoretic analysis of the fiber-optical channel, we deploy two common assumptions on the channel model. First, the dispersion is set to zero and second, a sampling receiver is used to obtain discrete-time models from continuous-time channels. These two assumptions were first applied to the NLS equation in [
1] to obtain an analytically tractable channel model. This channel model was developed also in [
20,
21,
22] using different methods. In this paper, we refer to this model as MNC.
In [
21], a lower bound on the per-sample capacity of the memoryless NLS channel is derived, which proves that the capacity goes to infinity with power. In [
22], the capacity of the same channel is evaluated numerically. Furthermore, it is shown that the capacity pre-log is
. Approximations of the capacity and optimal input distribution in the intermediate power range are derived in [
23,
24]. These results are extended to a channel with a more realistic receiver than the sampling one in [
25]. The only known nonasymptotic upper bound on the capacity of this channel is
(bits per channel use) [
12], where
is the signal-to-noise ratio. This upper bound holds also for the general case of nonzero dispersion.
The novel contributions of this paper are as follows. First, we tightly bound the capacity of the RPC model and prove that its capacity pre-log is 3. Second, the capacity of the LPC is readily shown to be the same as that of an AWGN channel with the same input and noise power. Hence, the capacity pre-log of the LPC is 1. Third, we establish a novel upper bound on the capacity of the MNC (first presented in the conference version of this manuscript [
26]). Our upper bound improves the previously known upper bound [
12] on the capacity of this channel significantly and, together with the proposed lower bound, allows one to characterize the capacity of the MNC accurately.
Although all three models represent the same physical optical channel, their capacities behave very differently in the high-power regime. This result highlights the profound impact of the simplifying assumptions on the capacity at high powers, and indicates that care should be taken in translating the results obtained based on these models into guidelines for system design.
The rest of this paper is organized as follows. In
Section 2, we introduce the three channel models. In
Section 3, we present upper and lower bounds on the capacity of these channels and establish the capacity pre-log of the perturbative models. Numerical results are provided in
Section 4. We conclude the paper in
Section 5. The proofs of all theorems are given in the appendices.
Notation: Random quantities are denoted by boldface letters. We use to denote the complex zero-mean circularly symmetric Gaussian distribution with variance . We write , , and x to denote the real part, the absolute value, and the phase of a complex number x. All logarithms are in base two. The mutual information between two random variables and is denoted by . The entropy and differential entropy are denoted by and , respectively. Finally, we use * for the convolution operator.
2. Channel Models
The fiber-optical channel is well-modeled by the NLS equation, which describes the propagation of a complex baseband electromagnetic field through a lossy single-mode fiber as
Here,
is the complex baseband signal at time
t and location
z. The parameter
is the nonlinear coefficient,
is the group-velocity dispersion parameter,
is the attenuation constant,
is the gain profile of the amplifier, and
is the Gaussian amplification noise, which is bandlimited because of the inline channel filters. The third term on the left-hand side of (
1) is responsible for the channel memory and the fourth term for the channel nonlinearity.
To compensate for the fiber losses, two types of signal amplification can be deployed, namely, distributed and lumped amplification. The former method compensates for the fiber loss continuously along the fiber, whereas the latter method boosts the signal power by dividing the fiber into several spans and using an optical amplifier at the end of each span. With distributed amplification, which we focus on in this paper, the noise can be described by the autocorrelation function [
2]
Here,
is the spontaneous emission factor,
h is Planck’s constant, and
is the optical carrier frequency. In addition,
is the Dirac delta function and
, where
is the noise bandwidth. In this paper, we shall focus on the ideal distributed-amplification case
.
We use a sampling receiver to go from continuous-time channels to discrete-time ones. A comprehensive description of the sampling receiver and of the induced discrete-time channel is provided in ([
22], Section III). Here, we review some crucial elements of this description for completeness. Assume that a signal
, which is band-limited to
hertz, is transmitted through a zero-dispersion NLS channel ((1) with
) in the time interval
. Because of nonlinearity, the bandwidth of the received signal
may be larger than that of
. To avoid signal distortion by the inline filters, we assume that
is set such that
is band-limited to
hertz for
. Since
, assuming
, both the transmitted and the received signal can be represented by
equispaced samples. The transmitter encodes data into subsets of these samples of cardinality
, referred to as the principal samples. At the receiver, demodulation is performed by sampling
at instances corresponding to the principal samples. This results in
parallel independent discrete-time channels that have the same input–output relation.
The sampling receiver has a number of shortcomings [
27] and using it should be considered a simplification. The resulting discrete-time model is used extensively in the literature (see, for example, [
1,
20,
21,
22,
28,
29]), since it makes analytical calculation possible. In this paper, we apply the sampling receiver not only to the memoryless NLS channel but also to the memoryless perturbative models.
Next, we review two perturbative channel models that are used in the literature to approximate the solution of the NLS Equation (
1). Among the multiple variations of perturbative models available in the literature, we use the ones proposed in [
30]. For both perturbative models, first continuous-time dispersive models are introduced, and then memoryless discrete-time channels are developed by assuming that
and by using a sampling receiver. Finally, we introduce the MNC model, which is derived from (
1) under the two above-mentioned assumptions.
Regular perturbative channel (RPC): Let
be the solution of the linear noiseless NLS equation (Equation (
1) with
and
). It can be computed as
, where
and
denotes the inverse Fourier transform. In the regular perturbation method, the output of the noiseless NLS channel (Equation (
1) with
) is approximated as
Here,
L is the fiber length and
is the nonlinear perturbation term. If now the model is expanded to include amplification noise as an additive noise component, neglecting signal–noise interactions, then the accumulated amplification noise
can be added to the signal at the receiver to obtain the channel model ([
14], Equation (
5))
The first-order approximation of
is ([
30], Equation (
13))
where the convolution is over the time variable. (Using higher-order nonlinear terms improves the accuracy of the regular perturbative channels. However, in this paper, we focus only on the channel model based on the first-order approximation, which is commonly used in the literature.) Neglecting dispersion (i.e., setting
), we have
and
. Using this in (
6), and then substituting (
6) into (
5), we obtain
Finally, by deploying sampling receiver, we obtain from (
7) the discrete-time channel model
Here,
,
is the total noise power, and
We refer to (
8) as the RPC.
Logarithmic perturbative channel (LPC): Another method for approximating the solution of the NLS Equation (
1) is to use logarithmic perturbation. With this method, the output signal is approximated as ([
14], Equation (
7))
where
is the same noise term as in (
5)–(
4). The first-order approximation of
is ([
30], Equation (19))
Under the zero-dispersion assumption (
), we have
and
. Using this in (
12), and then substituting (
12) into (
11), we obtain
Finally, by sampling the output signal, the discrete-time channel
is obtained, where
,
is given in (
9), and
is defined in (
10). We note that the channels (
8) and (
14) are equal up to a first-order linearization, which is accurate in the low-power regime. Furthermore, one may also obtain the model in (
13) by solving (
1) for
,
, and
and by adding the noise at the receiver.
Memoryless NLS Channel (MNC): Here, we shall study the underlying NLS channel in (
1) under the assumptions that
and that a sampling receiver is used to obtain a discrete-time channel. Let
and
be the amplitude and the phase of a transmitted symbol
, and let
and
be those of the received samples
. The discrete-time channel input–output relation can be described by the conditional probability density function (pdf) ([
20], Ch. 5) (see also ([
28], Sec. II))
The conditional pdf
and the Fourier coefficients
in (
15) are given by
Here,
denotes the
mth order modified Bessel function of the first kind, and
The complex square root in (
18) is a two-valued function, but both choices give the same values of
and
.
In the next section, we study the capacity of the channel models given in (
8), (
14), and (
15). Since all of these models are memoryless, their capacities under a power constraint
P are given by
where the supremum is over all complex probability distributions of
that satisfy the average-power constraint
4. Numerical Examples
In
Figure 1, we evaluate the bounds derived in
Section 3 for a fiber-optical channel whose parameters are listed in
Table 1. (The channel parameters are the same as in ([
22], Table I).) Using (
10), we obtain
.
As can be seen from
Figure 1, the capacity of the RPC is tightly bounded between the upper bound
in (
25) and the lower bound
in (
22). Furthermore, one can observe that although the alternative upper bound
in (
26) is loose at low powers, it becomes tight in the moderate- and high-power regimes.
We also plot the upper bound
on the capacity of the MNC. It can be seen that
improves substantially on the upper bound given in [
12], i.e., the capacity of the corresponding AWGN channel (
33) (which coincides with
). As a lower bound on the MNC capacity, we propose the mutual information in (
20) with an input
with uniform phase and amplitude
following a chi distribution with
k degrees of freedom. Specifically, we set
where
denotes the gamma function. The parameter
k is optimized for each power. (Due to the computational complexity, we only considered
k values from
to
in steps of
.) We calculated the bound numerically and include it in
Figure 1 (referred to as max–chi lower bound). We also include two lower bounds corresponding to
(with half-Gaussian amplitude distribution, first presented in [
22]) and
(with Rayleigh-distributed amplitude, or equivalently, a complex Gaussian input
, first presented in [
26]). The max–chi lower bound coincides with these two lower bounds at asymptotically low and high power, and improves slightly thereon at intermediate powers (around 0 dBm), similarly to the numerical bound in [
32]. Specifically, at asymptotically low powers,
(Gaussian lower bound) is optimal. This is expected, since the channel is essentially linear at low powers. At high powers, on the other hand, the optimal
k value approaches 1 (half-Gaussian lower bound), which is consistent with [
22], where it has been shown that half-Gaussian amplitude distribution is capacity-achieving for the MNC in the high-power regime. Based on our numerical evaluations, we observed that
maximizes the max–chi lower bound (among the set of considered values of
k) in the power range
dBm. Finally, in
Figure 1, we plot the lower bound based on the input distribution ([
23], Equation (45)). As can be seen, based on our numerical evaluation, this lower bound almost coincides with the max-chi lower bound at low powers (
) and improves on it in the intermediate power range (
dBm); however, it is suboptimal at high powers (
dBm).
Figure 1 suggests that
experiences changes in slope at about 0 and 30 dBm (corresponding to the inflection points at about
dBm and 20 dBm). To explain this behavior, we evaluate the phase and the amplitude components of the half-Gaussian lower bound. Specifically, we split the mutual information into two parts as
The first term in (39) is the amplitude component and the second term is the phase component of the mutual information. These two components are evaluated for the half-Gaussian amplitude distribution and plotted in
Figure 1. It can be seen from
Figure 1 that the amplitude component is monotonically increasing with power while the phase component goes to zero with power after reaching a maximum. Indeed, at high powers, the phase of the received signal becomes uniformly distributed over
and independent of the transmitted signal ([
33], Lem. 5). By adding these two components, one obtains a capacity lower bound that changes concavity at two points. The reduction of the capacity slope at intermediate powers is consistent with [
23,
24], where it is shown that the capacity grows according to
in this regime.
As a final observation, we note that diverges from at about dBm, whereas diverges from at about dBm. Since the MNC describes the nondispersive NLS channel more accurately than the other two channels, this result suggests that the perturbative models are grossly inaccurate in the high-power regime.