Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Frequency-resolved adaptive probabilistic shaping for DMT-modulated IM-DD optical interconnects

Open Access Open Access

Abstract

For decades, advanced modulation techniques have been proposed to increase the capacity for intensity-modulation and direct-detection (IM-DD) optical fiber interconnection systems. Typically, the frequency-resolved discrete multi-tone (DMT) modulation was proposed by loading modulations with different bit numbers to fit the channel’s frequency response. Capacity can thus be better improved through finer use of the signal-to-noise-ratio (SNR) distribution in the frequency domain. For conventional DMT, the constellations loaded on individual subcarriers are all equip probability distributed. In this work, we propose a probabilistically shaped DMT (PS-DMT) modulation with adaptively loaded entropies referring to channel frequency response for short-reach optical interconnects. Achievable information rate (AIR) improvements of PS-DMT with both Maxwell-Boltzmann and dyadic distribution are investigated based on generalized mutual information (GMI). Moreover, the proposed PS-DMT has been realized experimentally over a multimode optical link using vertical-cavity surface-emitting lasers (VCSELs) with 100-m-long multimode fiber (MMF) transmission. This method can significantly improve the signaling capacity since two significant benefits are simultaneously utilized: 1) the shaping gain of PS at limited SNR condition and 2) the frequency-resolved continuous entropy loading for better fitting to the channel frequency response. Improved capacity, in terms of AIR, can thus be expected for a practical channel when using PS-DMT. This method can potentially be extended to a wide range of application scenarios, including both multimode and single-mode IM-DD fiber-optic communications.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

In recent years, the continuously and explosively increasing amounts of data due to rising applications such as cloud computing/storage, 4G/5G wireless communication, artificial intelligence and so on have been boosting the demand for high-capacity optical communications, particularly for short-reach optical interconnects between/within data centers and high-performance computers. A wide variety of methods, including advanced modulation [1–6], error correction coding [7–12], and digital signal processing (DSP) [13–15], have been studied in attempts to improve the system capacity. For short-reach communication scenario, the intensity-modulation and direct-detection (IM-DD) solution provides great advantages of low power consumption and high density. While for most practical IM-DD channels, frequency response is usually non-flat even fluctuated in the frequency domain. That is caused by all possible contributions which influences channel frequency response, like limited bandwidth of optical and electrical devices and fiber dispersion. Discrete multi-tone (DMT) modulation was proposed as a way to address this problem, by loading modulations with different bit numbers on individual subcarriers with reference to the channel’s signal-to-noise-ratio (SNR) response [16,17]. However, for conventional DMT, the constellations of individual subcarriers are all equip probability distributed.

In fact, a uniform distribution with equip probability demonstrates high achievable information rate (AIR) at a high SNR. In noisy conditions, like most practical optical links, the optimal distributions with the highest capacities are not always uniform [18]. Particularly for DMT, the non-uniform distributions with the highest AIRs would also differ for different modulations of different subcarriers, according to their corresponding SNRs. Moreover, for the bit-loading process of DMT with targeted spectral efficiency, the optimal loaded entropies are continuous decimals rather than pure integers based on Fischer’s theory [16]. As a result, during the bit-loading process of conventional DMT, the stepped bit numbers (integers) cannot fit the channel response accurately. It requires extra power allocation to alter the subcarriers’ SNRs in return.

At the same time, probabilistic shaping (PS) method is widely utilized in coherent systems to adjust the distributions of the constellations, so that corresponding AIR can be improved effectively [19–24]. Very recently, an entropy loading scheme of DMT was reported to fit channel SNR response based on PS [25–27]. For [25], entropy-loaded DMT was realized by using an IQ modulator assisted with Kramers-Kronig receiver. In [26], entropy loading is proposed after water-filling power adaptation, to get rid of extra power reallocation. Corresponding investigations are done on a coherent system with maximum subcarrier number of 16. In [27], entropy-loaded DMT without power allocation is proposed by using Maxwell-Boltzmann (MB) distribution in visible light communication system, with AIR of 204 Mb/s. For IM-DD optical interconnection applications, there is still a research blank about the PS for frequency-resolved entropy loading, especially for achieving 100-Gbps per lane.

Here we demonstrate a frequency-resolved adaptive probabilistic shaping method refer to channel frequency response, for the 112-Gbps DMT-modulated IM-DD optical interconnection system. The continuously bit (in terms of entropy) loading is realized by adapted probability distributions, allowing for better fitting to channel frequency response with simultaneous shaping gain. Proof-of-concept studies have been carried out via both simulations and experiments. Performances of PS-DMT at a raw rate of 112 Gbps was investigated theoretically based on Maxwell-Boltzmann and dyadic distributions. The corresponding optical signaling was realized experimentally by direct modulation of a vertical-cavity surface-emitting laser (VCSEL) at 112 Gbps with 100-m multimode fiber (MMF) transmission successfully achieved.

2. Principle

For a practical optical channel, frequency responses always exhibit fading characteristic at high frequencies due to fiber dispersion, constrained bandwidth of optical and electrical devices. As a result, corresponding frequency response may roll down even fluctuate which resulting in non-flat response, as illustrated by the orange curve shown in Fig. 1(a). In this case, conventional DMT modulation maximizes capacity by loading paralleled bits on the frequency-resolved subcarriers with reference to the SNR-to-frequency response. Therefore, DMT modulation intuitively produces a higher capacity for such channels, than other formats, such as pulse-amplitude modulation (PAM) and carrier-less amplitude phase modulation. For DMT, the individual subcarriers are loaded using quadrature amplitude modulation (QAM) signals with different entropies. Typical AIR curves with respect to the SNR are briefly plotted in Fig. 1(b) for a non-uniform and the uniform probability distributions [19]. The QAM signals with the non-uniform probability distribution have a higher AIR in the low SNR region. Moreover, the non-uniform distribution can be adapted refer to different SNR.

 figure: Fig. 1

Fig. 1 AIR (in terms of GMI) gain by PS-DMT. (a) Frequency responses for ideal (green line) and practical channels (orange curve); (b) AIR with respect to SNR for the uniform distribution and a non-uniform distribution, along with the Shannon limit; (c) Corresponding AIR-frequency responses for the ideal (green) and practical (orange) channels. The solid curves indicate the Shannon limit. The dotted and dashed curves represent the AIR curves for the uniform distribution and adaptive non-uniform distribution, respectively.

Download Full Size | PDF

Consequently, for an ideal channel with a flat and high SNR response, AIR can be improved by loading with more bits per symbol. However, this is always a limited option during practical implementation. For DMT-modulated system, AIR-to-frequency curve can be obtained by known SNR and distribution. The AIR curves with respect to frequency can then be plotted as shown in Fig. 1(c). The green solid straight line represents Shannon limit for such a channel. The limited number of loaded bits leads to the AIR gap to Shannon limit, as indicated by the green dotted straight line. If using PS with the non-uniform probability distribution optimized for a low SNR, the AIR may be reduced because of such mismatched distribution for a high SNR, as indicated by the green dashed straight line. However, AIR performance is totally different for practical channels with constrained and fluctuating SNR responses, as illustrated by the orange curves shown in Fig. 1(c). When SNR is constraint, PS can offer considerable AIR gain, which more approaching to Shannon limit. The orange solid curve depicts the corresponding Shannon limit. Significant AIR gain can be obtained because of the shaping gain of the PS, particularly in the low SNR region, as indicated by the orange dotted and dashed curves.

In addition, for conventional DMT, the loaded bit numbers must be step integers to perform bit-to-constellation mapping. This is the reason why rounding-off is required during the conventional bit loading process. As a result, it requires extra power allocation of subcarriers to alter the corresponding SNRs, for adapting the channel response, as plotted in Fig. 2(a). When taking the digital subcarriers into the whole-link analysis, the adapted SNR-frequency response after ideal power reallocation should be flat for those subcarriers with same loaded bit numbers. However, such adaptation is always not precise enough for perfect channel fitting. The proposed PS-DMT method can get rid of power reallocation because continuous bit loading (in terms of entropy) can be realized without any reduction of net rate. Therefore, the curve of the effectively loaded bit can perfectly fit the channel response, with all subcarriers’ power unchanged, as shown in Fig. 2(b).

 figure: Fig. 2

Fig. 2 (a) Conventional DMT using discrete bit loading and power allocation; (b) PS-DMT using entropy loading.

Download Full Size | PDF

For conventional DMT, to achieve a balanced signal quality for all subcarriers of a DMT symbol with target spectrum efficiency, more bits are loaded on the subcarrier with the higher SNR, and fewer bits are loaded on those subcarriers with lower SNRs. The bit-loading process can be performed numerically based on Chow and Fischer algorithm [16]. The main reason why chose Fischer’s algorithm in this work, is to obtain the decimal entropies to be loaded on every subcarriers. The optimal loaded bit number Bi (which is always decimal) of the ith subcarrier can be calculated based on a known noise power for a DMT symbol, as described by Eq. (1) in [16], where Btotal is the total bit numbers to be loaded, M is the total number of subcarriers, and Ni is the noise power for the ith subcarrier assuming normalized signal power. Consequently, with known the SNR for a typical subcarrier, Ni can be directly obtained by calculated the ratio of adapted subcarrier’s power and corresponding SNR. For conventional DMT, the bit loading process includes channel SNR estimation, rounding-off function of Bi and power reallocation of individual subcarriers, as shown in Fig. 3.

 figure: Fig. 3

Fig. 3 Bit loading metrologies of conventional DMT and proposed PS-DMT schemes.

Download Full Size | PDF

Bi=Btotal+log2NiMlog2Ni

While for PS-DMT, the loaded entropies on the subcarriers can be arbitrarily decimal through careful selection of the non-uniform distributions. An adaptive approach to the optimal loaded bit number Bi of Eq. (1) can thus be achieved. In this way, entropies can be loaded very accurately with reference to the channel response. Additionally, the AIR gain produced by PS is significant when the SNR is limited to a low value. Benefitting from these two advantages, the proposed PS-DMT technique is expected to be an effective enabling method to approach the Shannon limit for SNR-constrained and fluctuating channels.

The proposed bit-loading methodology for PS-DMT is illustrated in the flowchart as shown in Fig. 3. A rounding-up function is performed for Bi to obtain standard constellations for QAM-N noted as Bi*, and N denotes the constellation number. The corresponding bit-to-symbol mapping is performed as N=2Bi*. Finally, the optimization problem must be solved to search for the distribution of the 1D PAM signals. The optimization object is to obtain the 1D probability distribution with entropy that is most close to Bi2. These distributions are subject to the Maxwell-Boltzmann (MB) equation P = e−μ|x| with a variable μ. Therefore, the optimization problem can be simplified as:

minμ(Plog2PBi2)s.t.P=e-μ|x|

Because in the case of Bi* = 3, two orthogonal components of QAM-8 are different: one is OOK, another is PAM-4. In this case, Eq. (2) with assuming the same modulation for two orthogonal parts, is no longer valid. As a result, Eq. (3) is stated specifically for the case of Bi* = 3, aiming to find the distribution P of the PAM-4 component:

minμ(Plog2P2)s.t.P=e-μ|x|

For quantitatively evaluating performance of proposed PS-DMT, there are several metrics like bit-error rate (BER) before FEC, Q factor, symbol-error rate (SER) and error vector magnitude (EVM). Very recently, AIR is proved as a powerful metric for coded-modulation systems, because it can indicate how much information can be reliably transmitted per symbol, and also exhibits inherent relationship to FEC [28]. In the case of bit-wise decoding, generalized mutual information (GMI) can represent AIR. Under this circumstance, normalized GMI can be used for predicting the code rate of optimal FEC for achieving error free [29,30]. Consequently, we use GMI and NGMI as the performance metrics in this work.

In fact, for IM-DD optical links, the inter-symbol interference (ISI) caused by limited bandwidth and fiber dispersion will induce non-Gaussian distortion. However, DMT separates its spectrum by several subcarriers, which can get rid of ISI. In this case, noise for every subcarrier can be approximated as Gaussian type. As a result, GMI value of the ith subcarrier can be calculated by Eq. (4) [31], in which, m is the bit number (equal to entropy) of a QAM symbol, i is the index of subcarrier, Bki is the kth bit of QAM symbol for the ith subcarrier, Lki is estimated log-likelihood ratio (LLR) of the kth bits for the ith subcarrier.

GMIi=k=1mI(Bki;Lki)

To evaluate the AIR performance of the proposed PS-DMT scheme, the optical signaling is modeled based on the additive white Gaussian noise (AWGN) channel with a constrained bandwidth. The attenuation at frequency domain is characterized using a Gaussian low-pass filter as described by Eq. (5), where SNRf = 0 is defined as the SNR value at DC, and BW3dB is the 3-dB bandwidth of the channel.

F(f)=SNRf=0e2.54f2BW3dB2

Based on known SNRf = 0 and BW3dB, the bit loading of the proposed PS-DMT can then be realized through solution of the optimization problem above. Typically, BW3dB is set at 25 GHz, with SNRf = 0 of 20 dB. With the aim of providing the raw rate of 112 Gbps, the bit loading results of both conventional DMT and the proposed PS-DMT scheme are investigated, as plotted in Fig. 4(a). The red circles are the loaded bit numbers for conventional DMT, which are all stepped integers. Moreover, we take the digital subcarriers into consideration for analyzing the whole-link SNR response. The reason why calculating the whole-link response with extra consideration of power-altered subcarriers, is for evaluating the channel adaptation of conventional DMT and proposed PS-DMT. Because loaded bit numbers for conventional DMT are integer only, power reallocation has to be performed to alter the subcarriers’ power (raw-tooth shape) to maintain same SNRs for those subcarriers with same loaded bit numbers, aiming to assure equal and optimal performance margin over all subcarriers according to the statement in Column 18, Page 2 of [17]. As a result of which, SNR response (including digital subcarriers) after power reallocation should be in the same shape to loaded bit numbers of subcarriers, which is flat for those subcarriers loaded by same bit numbers. The adapted channel response (including digital subcarriers) after power reallocation is shown as red line in Fig. 4(a). For PS-DMT, bit loading is performed by solving the optimization problem stated in the article. The entropies adapted to the channel response are obtained for all 127 subcarriers and plotted as the blue triangles in Fig. 4(a). The adapted distribution (MB) of the 60th subcarrier is shown in Fig. 4(b), with entropy equal to 4.975.

 figure: Fig. 4

Fig. 4 (a) Simulated channel SNR response (blue line), adapted channel response using power allocation (red line), and entropy loading results by PS-DMT; (b) The probability distribution of 60th subcarrier for PS-DMT.

Download Full Size | PDF

Under AWGN condition, theoretical GMI values are calculated for both conventional DMT and PS-DMT. For DMT signals, GMI calculation is performed over all the subcarriers. The GMI of a single subcarrier can be obtained by Eq. (4). The GMI values for all 127 subcarriers are plotted in Fig. 5(a), with BW3dB of 25 GHz and SNRf = 0 of 20 dB. It indicates that the AIR values of PS-DMT (MB distribution) are much higher than those of conventional DMT, at lower-order subcarriers (1st to 80th). While for higher-order subcarriers, shaped QAM signals has lower AIRs than uniform ones. In addition, the MB distributions require large-length block to perform PS coding. Here, we use Geometric Huffman coding (GHC) to match dyadic distributions to MB ones. GHC aims to obtain the optimized dyadic distribution, whose Kullback-Leibler distance to desired distribution is minimized. The detailed descriptions about GHC is given in [32]. To obtain dyadic distribution, a Huffman tree requires to be constructed with updating rule as Eq. (6), with assumption of p1 p2 p3 pm , and p' is the updated root during every interaction.

 figure: Fig. 5

Fig. 5 (a) Loaded entropy of conventional DMT (red circles), PS-DMT (MB) (blue triangles) and PS-DMT (dyadic) (black circles), and corresponding GMI; (b) The GMI-to-SNR curve for the 60th subcarrier.

Download Full Size | PDF

p'={pm1,2pm1pm,ififpm14pmpm1<4pm

Loaded entropies of PS-DMT (dyadic) are plotted as light gray dots in Fig. 5(a). Although the entropies of PS-DMT (dyadic) cannot fit channel response perfectly, positive AIR gain still can be obtained for lower-order subcarriers. Moreover, a phenomenon is observed that dyadic distribution performs better at low frequencies rather than MB distribution. We think that is because that the optimal entropies calculated by Eq. (1) exceed the Shannon capacity for those subcarriers. While, dyadic distribution has more probabilities gathering to the center of constellation. As a consequent, dyadic distribution has higher AIR than MB for those subcarriers. To investigate the reason of shaping gain more clearly, the GMI-to-SNR curve for the 60th subcarrier (QAM32) is plotted in Fig. 5(b). The dark rectangular marks the zone with positive shaping gain of PS. Without power reallocation, SNR at the 60th subcarrier is 17.46 dB. Adapted SNR is 17.06 dB after power reallocation. Thus the GMI gain of PS-DMT is 0.4 bit/symbol for the 60th subcarrier.

To further investigate the performance improvement of PS-DMT, total GMI value for a DMT symbol (contain 127 subcarriers) can be obtained by summing up GMI values of all subcarriers. Thus, the AIR improvement contours for PS-DMT (MB) and PS-DMT (dyadic) are plotted versus both SNR and 3-dB bandwidth in Figs. 6(a) and 6(b), respectively.

 figure: Fig. 6

Fig. 6 AIR gain versus varying SNR and bandwidth: (a), PS-DMT (MB); (b), PS-DMT (dyadic).

Download Full Size | PDF

Figure 6(a) indicates that the maximum AIR improvement over 189 bits/symbol can be obtained by using PS-DMT (MB) for bandwidths below 20 GHz, when SNR is lower than 15 dB. When bandwidth is over 20 GHz, the AIR improvement over 25 bits/symbol can be obtained due to the shaping gain of MB distributions. The improvement of data rate is not strictly monotone increasing with SNR. It can be explained that the positive shaping gain is achieved at certain SNRs, typically shown in Fig. 5(b). For dyadic distributions, the data rate improvement of 29.1 bits/symbol also can be achieved at a broad region with bandwidth over 25 GHz.

3. Experimental results and discussion

The experiment investigation of a 112-Gbps optical interconnection is carried out using a VCSEL-MMF link. The experimental setup is illustrated in Fig. 7 below. The gray box of Fig. 7 encloses the channel SNR estimation process, which is based on sending a prior-known multi-tone signal (QPSKs) to enable calculation of the SNR for every subcarrier. In the hardware part, the arbitrary waveform generator (AWG) converts the digital sequences into analogue signal to drive the VCSEL chip, with a DC source provided for tuning bias. The 3D alignment platform is for both electrical and optical coupling of the VCSEL. At the receiver, a photo detector (PD) is used for signal detecting, and a digital storage oscilloscope (DSO), is used to sample the detected signal for processing. The bit-loading process is realized using the Fischer algorithm (for the conventional DMT) and the proposed continuous bit-loading algorithm.

 figure: Fig. 7

Fig. 7 Experimental setup of VCSEL-MMF optical link.

Download Full Size | PDF

The sampling rate of the AWG is set at a fixed value of 54 Gsa/s. To achieve the target raw rate of 112 Gbps using 127 subcarriers, the total loaded bit number required for a DMT symbol is 112/54*(128*2 + 10) = 552. The term “128*2” represents the Hermitian transformation and “10” is the length of the cyclic prefix. Based on the measured SNR responses, bit-loading processes are performed on 127 subcarriers for both conventional DMT and PS-DMT. In the optical B2B link case, the bit-loading results for conventional DMT are shown as the red points in Fig. 8, while the bit-loading results for the proposed PS-DMT (MB) scheme are plotted as blue points. Entropies of dyadic shaped PS-DMT are marked by blue lines. The adapted SNR response by power reallocation is re-measured by sending multi-tone QPSK signals with reallocated power, plotted as pink line in Fig. 8. In fact, the adapted response is not precisely matched to loaded bit numbers for conventional DMT. Fortunately, PS-DMT can get rid of this deviation, by directly adapting entropies to channel response. The corresponding optimized distributions for 12th and 66th subcarriers are shown in Fig. 8(a), respectively. For MB distribution, a distribution matcher (DM) is required for coding uniform-distributed binary sequence to the sequence whose output symbol is non-uniformly distributed [33]. While for dyadic distributions, corresponding non-uniform distribution can be obtained directly by binary mapping. For an example, in order to obtain a symbol with probability of 2-L, L-length mapping has to be performed on binary sequence for this symbol. In practical case, such a variance-length coding may result in severe synchronization problem in the receiver. Insertion of ambiguity bits can ideal with this issue through maintaining same bit number for every symbol. The dyadic distributions are also inserted in Fig. 8(b). In optical B2B case, the entropies of MB and distribution for the 12th subcarrier are 6.613 and 6.5, respectively. For the 66th subcarrier, the entropies of MB and distribution are 5.3755 and 5.5. After 100-m OM3 fiber transmission, the bit-loading results and the optimized distributions are plotted as shown in Fig. 9. The entropies of MB and distribution for the 12th subcarrier are 6.9939 and 7, respectively. For the 66th subcarrier, the entropies of MB and distribution are 5.316 and 5.5. The channel fading becomes severe because of modal dispersion, and thus more bits are allocated on the lower-frequency subcarriers. These results show that continuous entropy loading provides a better fit to the SNR response when compared with the conventional DMT.

 figure: Fig. 8

Fig. 8 Experimental bit loading results for DMT and PS-DMT in the optical B2B case: (a) Bit-loading results for conventional DMT (red dots), PS-DMT (blue dots) and PS-DMT (dyadic) (blue line). (b) Shaped probability distributions of two typical subcarriers (12th and 66th subcarriers) for PS-DMT.

Download Full Size | PDF

 figure: Fig. 9

Fig. 9 Experimental bit loading results for DMT and PS-DMT after 100-m transmission: (a) Bit-loading results for conventional DMT (red dots), PS-DMT (blue dots) and PS-DMT (dyadic) (blue line). (b) Shaped probability distributions of two typical subcarriers (12th and 66th subcarriers) for PS-DMT.

Download Full Size | PDF

With an optical power of 3.5 dBm for the optical B2B case, the demodulated constellations for both DMT and PS-DMT at the receiver are plotted as shown in Fig. 10. In fact, the demodulated constellations will rotate if no synchronization is performed. Thus, we estimated the phase-to-frequency response during training process, for correcting the rotated constellations. Corrected constellations are drawn for two selected typical subcarriers (12nd and 66th). In addition, right side of Fig. 10 shows the constellations of these typical subcarriers after 100-m OM3 fiber transmission with a received optical power of 3.3 dBm. The improvement in the signal quality related to the PS shaping gain can be observed visually by the clearer constellations that appear after PS when compared with those obtained before PS. It can also be seen that the constellations become clearer after PS. This occurs because the Euclidean distances between the symbols are broadened, with more symbols being gathered at the center, when the average power is fixed. After 100-m OM3 fiber transmission, however, the channel-efficient bandwidth is reduced, and more bits are loaded on the lower-frequency subcarriers. As a result, the constellations are much noisier after the 100-m transmission. Additionally, the shaped constellations are obviously clearer, with more bits being allocated at their centers. Moreover, one may observe that the shaped constellations are not ideally circular, that is because constellations become distorted after modulation and fiber propagation.

 figure: Fig. 10

Fig. 10 Constellations of 12th and 66th subcarriers for optical B2B case (3.5-dBm received power) and after 100-m OM3 fiber transmission (3.3-dBm received power).

Download Full Size | PDF

Experimental GMI values for signals on individual subcarriers are calculated by Eq. (4). In detail, with assuming noise on every subcarrier as Gaussian type, LLR for the kth bits for the ith subcarrier can be estimated by Eq. (7), in which Xk1 means the real value of symbol whose kth bit is ‘1’, Yk1 is the imaginary value of symbol whose kth bit is ‘1’

Lk=log2e((xXk1)2σx2+(yYk1)2σy2)e((xXk0)2σx2+(yYk0)2σy2)

The GMI curve of the compensated signal has a similar shape to that of the bit-loading numbers, but with a reduced value because of noise. In addition, after PS, the GMI curve becomes much smoother, as shown in Fig. 11(a), and has a similar shape to the loaded entropies. The proposed continuous bit-loading algorithm can adapt the loaded entropies to fit the channel’s SNR response well, and thus the calculated GMI curve is much smoother for PS-DMT. Moreover, the theoretical GMI values of individual subcarriers are also plotted as dashed dots for comparison. The experimental results are in the similar shape to theoretical ones but smoother. After 100-m OM3 fiber transmission, the GMI curves for both conventional DMT and PS-DMT are plotted as shown in Fig. 11(b). Moreover, total GMI values for all subcarrier are summed up for calculating GMI of DMT symbol. In optical B2B case, GMI values for conventional DMT, PS-DMT (MB) and PS-DMT (dyadic) are 476.95 bits/symbol, 495.25 bits/symbol and 486.74 bits/symbol. Consequently, PS-DMT by using MB and dyadic distributions can offer 18.3 bits/symbol and 9.78 bits/symbol AIR gain, respectively.

 figure: Fig. 11

Fig. 11 Generalized Mutual information (GMI) values: (a) Conventional DMT and PS-DMT for the optical B2B case, with received optical power of 3.5 dBm; (b) Conventional DMT and PS-DMT after 100-m OM3 MMF transmission, with received optical power of 3.3 dBm.

Download Full Size | PDF

As stated before, NGMI can be used for predicting the post-FEC BER. Moreover, information can be reliably recovered only when  RcNGMI, in which Rc is the FEC code rate, according to [29,30]. As a result, the NGMI limit of 25%-overhead soft-decision FEC (25% SD-FEC) is 0.8 [29,30]. Consequently, in order to investigate the power sensitivity gain induced by PS-DMT, NGMI values are then calculated with varying received optical power values, as shown in Fig. 12. If there is no noise and signal distortion, NGMI value equals to 1. As the optical power decreases, NGMI also decreases rapidly because of the low SNR. In the optical B2B case, with an optical power of more than 1 dBm, the SNR is sufficient to keep the NGMI value stable. However, for optical powers below 1 dBm, any reduction in the SNR also reduces NGMI values of high-density constellations, which results in sharply reduced NGMI values. After 100-m OM3 fiber transmission, more bits are allocated to the lower-frequency subcarriers because the bandwidth is more constrained. In this case, when the optical power decreases, the higher-density constellations with larger bit-loading numbers suffer greater NGMI reduction than in the optical B2B case. To enable error-free signaling at 112 Gbps assuming 25% SD-FEC, power sensitivity gains of 0.7 dB and 0.64 dB can be obtained for the optical B2B case and the 100-m OM3 fiber transmission case, assisted by PS-DMT (MB). For PS-DMT (dyadic), power sensitivity gains are 0.44 dB and 0.29 dB for optical B2B and 100-m transmission.

 figure: Fig. 12

Fig. 12 Experimental NGMI values and data rate of reliable communication under various received optical powers.

Download Full Size | PDF

4. Conclusion

In this paper, an optical signaling method to improve the data rate of the DMT-modulated IM-DD optical interconnection links is proposed based on frequency-resolved adaptive probabilistic shaping. This method can significantly improve the signaling capacity since two significant benefits are simultaneously utilized: 1) the shaping gain of PS for carriers with limited signal-to-noise-ratio (SNR) and 2) the frequency-resolved continuous entropy loading for better fitting to the channel frequency response. Proof-of-concept investigations have been carried out via both simulations and experiments, for MB and dyadic distributions. AIR improvements of MB and dyadic distribution versus varying SNR and bandwidth were investigated theoretically. In addition, optical signaling was realized experimentally using a commercial VCSEL at 112 Gbps data rate with 100-m OM3 fiber transmission achieved. 0.64-dB and 0.29-dB power gains have been achieved for 100-m OM3 fiber transmissions, by using MB and dyadic shaping of DMT, respectively.

Funding

National Natural Science Foundation of China (NSFC) (61875124, 61875049, 61675128).

Disclosures

The authors declare that there are no conflicts of interest related to this article.

References

1. G. D. Forney and G. Ungerboeck, “Modulation and coding for linear Gaussian channels,” IEEE J. Sel. Areas Comm. 44, 2384–2415 (1998).

2. Y. Cai, J. X. Cai, A. Pilipetskii, G. Mohs, and N. S. Bergano, “Spectral efficiency limits of pre-filtered modulation formats,” Opt. Express 18(19), 20273–20281 (2010). [CrossRef]   [PubMed]  

3. J. M. Kahn and K. P. Ho, “Spectral efficiency limits and modulation/detection techniques for DWDM systems,” IEEE J. Sel. Top. Quantum Electron. 10(2), 259–272 (2004). [CrossRef]  

4. P. J. Winzer, “High-spectral-efficiency optical modulation formats,” J. Lightwave Technol. 30(24), 3824–3835 (2012). [CrossRef]  

5. A. Ellis, “Modulation formats which approach the Shannon limit,” in Optical Fiber Communication Conference (Optical Society of America, 2016), paper OMM4.

6. L. Sun, J. Du, and Z. He, “Multiband three-dimensional carrierless amplitude phase modulation for short reach optical communications,” J. Lightwave Technol. 34(13), 3103–3109 (2016). [CrossRef]  

7. A. J. Goldsmith and S. G. Chua, “Adaptive coded modulation for fading channels,” IEEE Trans. Commun. 46(5), 595–602 (1998). [CrossRef]  

8. S. Y. Chung, G. D. Forney, T. J. Richardson, and R. Urbanke, “On the design of low-density parity-check codes within 0.0045 dB of the Shannon limit,” IEEE Commun. Lett. 5(2), 58–60 (2001). [CrossRef]  

9. S. Ten Brink, “Rate one-half code for approaching the Shannon limit by 0.1 dB,” Electron. Lett. 36(15), 1293–1294 (2000). [CrossRef]  

10. J. Garcia-Frias and W. Zhong, “Approaching Shannon performance by iterative decoding of linear codes with low-density generator matrix,” IEEE Commun. Lett. 7(6), 266–268 (2003). [CrossRef]  

11. C. Berrou, A. Glavieux, and P. Thitimajshima, “Near Shannon limit error-correcting coding and decoding: Turbo-codes,” in Proceedings of International Conference on Communications (IEEE, 1993), pp. 1064–1070. [CrossRef]  

12. D. J. MacKay and R. M. Neal, “Near Shannon limit performance of low density parity check codes,” Electron. Lett. 32(18), 1645–1646 (1996). [CrossRef]  

13. B. Schmidt, A. Lowery, and J. Armstrong, “Experimental demonstrations of electronic dispersion compensation for long-haul transmission using direct-detection optical OFDM,” J. Lightwave Technol. 26(1), 196–203 (2008). [CrossRef]  

14. L. Zhang, X. Hong, X. Pang, O. Ozolins, A. Udalcovs, R. Schatz, C. Guo, J. Zhang, F. Nordwall, K. M. Engenhardt, U. Westergren, S. Popov, G. Jacobsen, S. Xiao, W. Hu, and J. Chen, “Nonlinearity-aware 200 Gbit/s DMT transmission for C-band short-reach optical interconnects with a single packaged electro-absorption modulated laser,” Opt. Lett. 43(2), 182–185 (2018). [CrossRef]   [PubMed]  

15. Y. Yang, G. Cao, K. Zhong, X. Zhou, Y. Yao, A. P. T. Lau, and C. Lu, “Fast polarization-state tracking scheme based on radius-directed linear Kalman filter,” Opt. Express 23(15), 19673–19680 (2015). [CrossRef]   [PubMed]  

16. R. F. Fischer and J. B. Huber, “A new loading algorithm for discrete multitone transmission,” in Global Telecommunications Conference (IEEE, 1996), pp. 724–728. [CrossRef]  

17. P. S. Chow, J. M. Cioffi, and J. A. Bingham, “A practical discrete multitone transceiver loading algorithm for data transmission over spectrally shaped channels,” IEEE Trans. Commun. 43(2/3/4), 773–775 (1995). [CrossRef]  

18. C. E. Shannon and W. Weaver, The Mathematical Theory of Communication (University of Illinois Press, 1963).

19. T. Fehenberger, A. Alvarado, G. Böcherer, and N. Hanik, “On probabilistic shaping of quadrature amplitude modulation for the nonlinear fiber channel,” J. Lightwave Technol. 34(21), 5063–5073 (2016). [CrossRef]  

20. M. P. Yankov, D. Zibar, K. J. Larsen, L. P. Christensen, and S. Forchhammer, “Constellation shaping for fiber-optic channels with QAM and high spectral efficiency,” IEEE Photonics Technol. Lett. 26(23), 2407–2410 (2014). [CrossRef]  

21. S. Chandrasekhar, B. Li, J. Cho, X. Chen, E. Burrows, G. Raybon, and P. Winzer, “High-spectral-efficiency transmission of PDM 256-QAM with parallel probabilistic shaping at record rate-reach trade-offs,” in Proceedings of European Conference on Optical Communications (IEEE, 2011), pp.1–3.

22. C. Pan and F. R. Kschischang, “Probabilistic 16-QAM shaping in WDM systems,” J. Lightwave Technol. 34(18), 4285–4292 (2016). [CrossRef]  

23. D. Semrau, T. Xu, N. A. Shevchenko, M. Paskov, A. Alvarado, R. I. Killey, and P. Bayvel, “Achievable information rates estimates in optically amplified transmission systems using nonlinearity compensation and probabilistic shaping,” Opt. Lett. 42(1), 121–124 (2017). [CrossRef]   [PubMed]  

24. Y. Zhu, A. Li, W. R. Peng, C. Kan, Z. Li, S. Chowdhury, and Y. Bai, “Spectrally-efficient single-carrier 400G transmission enabled by probabilistic shaping,” in Optical Fiber Communication Conference (Optical Society of America, 2017), paper M3C.1. [CrossRef]  

25. X. Chen, J. Cho, S. Chandrasekhar, P. Winzer, C. Antonelli, A. Mecozzi, and M. Shtaif, “Single-wavelength, single-polarization, single-photodiode kramers-kronig detection of 440-Gb/s entropy-loaded discrete multitone modulation transmitted over 100-km SSMF,” in Photonics Conference (IEEE, 2017), pp. 1–2. [CrossRef]  

26. C. Di and S. William, “Approaching the Capacity of Colored-SNR Optical Channels by Multicarrier Entropy Loading,” J. Lightwave Technol. 36(1), 68–78 (2018). [CrossRef]  

27. C. Xie, Z. Chen, S. Fu, W. Liu, Z. He, L. Deng, M. Tang, and D. Liu, “Achievable information rate enhancement of visible light communication using probabilistically shaped OFDM modulation,” Opt. Express 26(1), 367–375 (2018). [CrossRef]   [PubMed]  

28. A. Alvarado, T. Fehenberger, B. Chen, and F. M. Willems, “Achievable information rates for fiber optics: Applications and computations,” J. Lightwave Technol. 36(2), 424–439 (2018). [CrossRef]  

29. G. Böcherer, “On Joint Design of Probabilistic Shaping and Forward Error Correction for Optical Systems,” in Proceedings of the Optical Fiber Communication Conference (Optical Society of America, 2018), paper M4E.1. [CrossRef]  

30. J. Cho, L. Schmalen, and P. J. Winzer, “Normalized generalized mutual information as a forward error correction threshold for probabilistically shaped QAM,” in Proceedings of European Conference on Optical Communications (IEEE, 2017), pp. 1–3. [CrossRef]  

31. A. Alvarado, E. Agrell, D. Lavery, R. Maher, and P. Bayvel, “Replacing the soft-decision FEC limit paradigm in the design of optical communication systems,” J. Lightwave Technol. 33(20), 4338–4352 (2015). [CrossRef]  

32. G. Bocherer and R. Mathar, “Matching dyadic distributions to channels,” in Proceedings of the 2011 Data Compression Conference (AMC, 2011), pp. 23–32. [CrossRef]  

33. P. Schulte and G. Böcherer, “Constant composition distribution matching,” IEEE Trans. Inf. Theory 62(1), 430–434 (2016). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (12)

Fig. 1
Fig. 1 AIR (in terms of GMI) gain by PS-DMT. (a) Frequency responses for ideal (green line) and practical channels (orange curve); (b) AIR with respect to SNR for the uniform distribution and a non-uniform distribution, along with the Shannon limit; (c) Corresponding AIR-frequency responses for the ideal (green) and practical (orange) channels. The solid curves indicate the Shannon limit. The dotted and dashed curves represent the AIR curves for the uniform distribution and adaptive non-uniform distribution, respectively.
Fig. 2
Fig. 2 (a) Conventional DMT using discrete bit loading and power allocation; (b) PS-DMT using entropy loading.
Fig. 3
Fig. 3 Bit loading metrologies of conventional DMT and proposed PS-DMT schemes.
Fig. 4
Fig. 4 (a) Simulated channel SNR response (blue line), adapted channel response using power allocation (red line), and entropy loading results by PS-DMT; (b) The probability distribution of 60th subcarrier for PS-DMT.
Fig. 5
Fig. 5 (a) Loaded entropy of conventional DMT (red circles), PS-DMT (MB) (blue triangles) and PS-DMT (dyadic) (black circles), and corresponding GMI; (b) The GMI-to-SNR curve for the 60th subcarrier.
Fig. 6
Fig. 6 AIR gain versus varying SNR and bandwidth: (a), PS-DMT (MB); (b), PS-DMT (dyadic).
Fig. 7
Fig. 7 Experimental setup of VCSEL-MMF optical link.
Fig. 8
Fig. 8 Experimental bit loading results for DMT and PS-DMT in the optical B2B case: (a) Bit-loading results for conventional DMT (red dots), PS-DMT (blue dots) and PS-DMT (dyadic) (blue line). (b) Shaped probability distributions of two typical subcarriers (12th and 66th subcarriers) for PS-DMT.
Fig. 9
Fig. 9 Experimental bit loading results for DMT and PS-DMT after 100-m transmission: (a) Bit-loading results for conventional DMT (red dots), PS-DMT (blue dots) and PS-DMT (dyadic) (blue line). (b) Shaped probability distributions of two typical subcarriers (12th and 66th subcarriers) for PS-DMT.
Fig. 10
Fig. 10 Constellations of 12th and 66th subcarriers for optical B2B case (3.5-dBm received power) and after 100-m OM3 fiber transmission (3.3-dBm received power).
Fig. 11
Fig. 11 Generalized Mutual information (GMI) values: (a) Conventional DMT and PS-DMT for the optical B2B case, with received optical power of 3.5 dBm; (b) Conventional DMT and PS-DMT after 100-m OM3 MMF transmission, with received optical power of 3.3 dBm.
Fig. 12
Fig. 12 Experimental NGMI values and data rate of reliable communication under various received optical powers.

Equations (7)

Equations on this page are rendered with MathJax. Learn more.

B i = B total + log 2 N i M log 2 N i
min μ (P log 2 P B i 2 ) s.t. P= e -μ| x |
min μ (P log 2 P2) s.t. P= e -μ| x |
GM I i = k=1 m I( B k i ; L k i )
F(f)=SN R f=0 e 2.54 f 2 B W 3dB 2
p'={ p m1 , 2 p m1 p m , if if p m1 4 p m p m1 <4 p m
L k = log 2 e ( (x X k 1 ) 2 σ x 2 + (y Y k 1 ) 2 σ y 2 ) e ( (x X k 0 ) 2 σ x 2 + (y Y k 0 ) 2 σ y 2 )
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.