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

Analytical and experimental performance evaluation of an integrated Si-photonic balanced coherent receiver in a colorless scenario

Open Access Open Access

Abstract

We study analytically and experimentally the performance limits of a Si-photonic (SiP) balanced coherent receiver (CRx) co-packaged with transimpedance amplifiers (TIAs) in a colorless WDM scheme. Firstly, the CRx architecture is depicted and characterization results are presented. Secondly, an analytical expression for the signal-to-noise ratio (SNR) at the CRx output is rigorously developed and various noise sources in the context of colorless reception are outlined. Thirdly, we study experimentally the system-level CRx performance in colorless reception of 16 × 112 Gbps PDM-QPSK WDM channels. Using a 15.5 dBm local oscillator (LO) power, error free transmissions over 4800 and 4160 km at received powers of −3 and −21 dBm per channel, respectively, were achieved in a fully colorless and preamplifierless reception. Next, a set of measurements on one of the center WDM channels is performed where the LO power, received signal power, distance, and number of channels presented to the CRx are swept to evaluate the performance limits of colorless reception. Results reveal that the LO beating with optical noise incoming with the signal is a dominant noise source regardless of received signal power. In the high received signal power regime (~0 dBm/channel), the self-beat noise from out-of-band (OOB) channels is an additional major noise source especially for small LO-to-signal power ratio, short reach and large number of OOB channels. For example, at a received signal power of 0 dBm/channel after 1600 km transmission, the SNR difference between the fully filtered and colorless scenarios, where 1 and 16 channels are passed to the CRx respectively, grows from 0.5 to 3.3 dB as the LO power changes from 12 to 0 dBm. For low received power (~-12 dBm/channel), the effect of OOB channels becomes minor while the receiver shot and thermal noises become more significant. We identify the common mode rejection ratio (CMRR) and sensitivity as the two important CRx specifications that impact the performance at high and low received signal power regimes, respectively. Finally, an excellent match between experimental and analytical SNRs is proven after the derived SNR model is fitted to the experimental data in a least-squares sense. The model is then used to predict that the CRx can operate colorlessly for a fully populated WDM spectrum with 80 channels provided that the LO-to-signal power ratio is properly set.

© 2014 Optical Society of America

1. Introduction

The incessantly growing capacity demand in optical transport has been spurring extensive research efforts in the realm of coherent optical technology [110]. One key enabling component of this technology is the coherent optical receiver (CRx) that linearly maps the in-phase (I) and quadrature (Q) components of the received light in both orthogonal polarizations into the electrical domain by beating the received field with that of a local oscillator (LO) laser [1]. In addition, a CRx offers higher sensitivity compared to direct-detection [2] and enables the exploitation of M-ary quadrature amplitude modulation (M-QAM) formats consequently providing the needed larger spectral efficiency [2,3]. When incorporated with high speed analog-to-digital converters (ADCs), it also enables mitigating transmission impairments as well as compensating for hardware imperfections using digital signal processing (DSP) [1,5,6].

By incorporating coherent technology along with frequency division multiplexing (FDM) on multiple optical carriers, up to 400G coherent transceivers have been recently commercialized [7]. However, it still remains a challenge to meet the relentless growing network traffic demand [11]. Technological solutions such as space division multiplexing (SDM) that utilizes multiple fiber modes / cores have the potential to multiply the capacity per fiber [8]. Although being heavily investigated throughout the research community, SDM requires a potential replacement of the field installed fiber and developing novel enabling components. This may hinder near future field deployments of SDM. Another promising technological solution that is probably more suited for the near future is the concept of an agile/flexible software-defined optical network [9,10]. In the envisioned agile network, flexible coherent transceivers, that are capable of adapting to dynamic traffic needs and thus maximizing the network capacity, are key elements that need to be deployed at ingress/egress nodes. These flexible transceivers need to be smart, data-rate and modulation format agile and impairment aware/tolerant. In addition, the notion of agility necessitates the ability of the transceiver to transmit / receive at any optical frequency in order to adapt to the dynamic nature of the network. Since external cavity lasers (ECLs) fully tunable over the C or L band are commercially available [12], they can be incorporated as transmitter and LO lasers and hence, in principle, enable the transceiver to operate at any wavelength. At the transmitter side, tuning the laser simply varies the transmission wavelength without any observable change in the transmitted signal quality provided that the electro-optics performs similarly across the tuning range. However on the receiver side, the situation is more complex. A CRx can in principle operate in a colorless fashion by utilizing the LO laser for channel selection. Such a scheme mandates removing the DEMUX filter before the receiver which is beneficial on one hand as it reduces cost and footprint besides offering the desired transceiver flexibility, but on the other hand, it may degrade the performance depending on various parameters e.g., the number of channels co-incident at the receiver with the channel of interest, the LO power, etc. Colorless reception in the context of coherent optical communication was studied from a theoretical standpoint more than two decades ago [13]; however, the absence of commercial narrow linewidth lasers and high speed integrated circuits (ICs) prevented any practical implementation. Recently, there have been some reports on experimental colorless coherent reception by utilizing either a conventional CRx front-end followed by single-ended detection [14], which suffers from limited maximum channel count and poor dynamic range due to the inability to reject the self-beat common mode noise, or a more complicated receiver front-end followed by single-ended detection to enable rejecting the self-beat noise [15]. Later in [16], Zhang et al. presented design guidelines and performance limits of colorless reception employing balanced detection based on simulation and analytical predictions. More recently in [17], the authors demonstrated experimentally the use of an integrated PLC-based balanced coherent receiver for colorless reception of one 100G channel among 80 WDM channels.

The recent tremendous progress of the photonic integration technology allows for cost effective and volume production of coherent receivers employing balanced photodetectors which are more suited for colorless operation. Recently, integrated dual-polarization phase-diversity balanced coherent receiver photonic integrated circuits (PICs) have been demonstrated in both Si [18,19] and InP [20,21]. Si has advantages over InP as a host material since it can be made into larger wafers allowing for lower cost chips [22]. In addition, the world of analog ICs now offers linear and differential transimpedance amplifiers (TIAs) with variable gain control and bandwidths commensurate with the baud rates at which a coherent transceiver operates. Recently, researchers reported TIAs packaged and wire bonded to the PIC in one module [19,20,23]. The TIAs enhance the receiver dynamic range and provide the target voltage swing required by the analog-to-digital converters (ADCs). In addition, TIAs accompanied with the high sensitivity of coherent reception potentially enable preamplifierless operation in which the optical preamplifier prior to the receiver is removed and essentially replaced by the electrical TIAs thus saving footprint and power consumption. Recently in our work in [24], we evaluated the sensitivity limits of the Si-photonic integrated balanced CRx presented in [19] when employed in colorless and preamplifierless reception of single channel 112 Gbps PDM-QPSK and 224 Gbps PDM-16QAM. Also, we have shown preliminary results of colorless reception of 16 × 112 Gbps PDM-QPSK WDM channels.

In this paper, we extend our work in [24] and study from theoretical and experimental standpoints the performance limits of the Si-photonic CRx when employed in a WDM colorless configuration. Firstly, we show the building blocks and introduce the key characterization results of the CRx assembly. Secondly, we develop rigorously an analytical expression of the SNR at the CRx output in the case of colorless reception taking into consideration all sources of noise namely, LO-ON, SIG-SIG, SIG-ON, ON-ON and LO-LO beatings as well as the receiver shot and thermal noises, where SIG stands for signal and ON stands for optical noise incoming with the signal. As will be later explained in details, the ON incoming with the signal contains contributions of both the amplified spontaneous emission (ASE) noise from in-line optical amplifiers and the fiber nonlinearity (NL) induced noise, which can be both modeled, according to [2527], as additive white Gaussian noise (AWGN) in the context of dispersion unmanaged transmission (UT). Compared to [16] where an analytical SNR expression was also given, the derivation performed in this manuscript, which is based on evaluating the autocorrelation functions and the power spectral densities of all terms in the differential photocurrent, is more rigorous. In addition, we propose an analytical expression for an effective common mode rejection ratio (CMRR) for the SIG-SIG beating term which takes into account the frequency dependent nature of the power imbalance and the skew mismatch between the P and N photodiodes of the CRx weighted by the power spectral density of the received signal intensity. This parameter was only arbitrarily defined in [16]. Furthermore, another distinguishing difference from the model in [16] is the inclusion of the impact of fiber nonlinearity as excess AWGN in addition to the conventional ASE noise from in-line EDFAs. Thirdly, we study experimentally the performance of the CRx in a fully colorless and preamplifierless reception of 16 × 112 Gbps PDM-QPSK WDM channels. Using a maximum LO power of 15.5 dBm, we evaluate the BER of each channel at various transmission distances at received signal powers of −3 and −21dBm per channel. Next, we restrict our measurements to one of the center channels of the WDM spectrum and perform an extensive set of measurements where we sweep the LO power, the received signal power, the distance and the number of channels presented to the CRx. This set of measurements aims at evaluating the performance limits of the CRx when used for colorless reception. All trends observed in the experimental results match those expected from simulations in [16]. In fact, results reveal that in the regime of high received signal power, the residual SIG-SIG beat noise from the out-of-band (OOB) channels degrades the performance especially for lower LO to signal power ratios and this degradation becomes relatively less pronounced as the transmission distance increases. In the regime of low received signal power, the shot noise and receiver thermal noise degrade the performance especially for very low LO powers and short transmission distances. Two key CRx specifications that dictate the achievable dynamic operating range are highlighted namely the responsivity and the CMRR, which impact the performance in the regimes of low and high received signal powers, respectively. Finally, we employ our experimental SNR results to fit the proposed analytical model to our experimental data using the least-squares criterion. An excellent match between experimental and analytical SNRs is found. We finally use the model to extrapolate the SNR curves beyond our experimental capabilities and predict the SNR when the number of channels is larger than 16. We conclude that the CRx allows for colorless selection of the test channel out of a fully populated WDM spectrum with 80 channels while achieving error-free operation after 4800 km transmission provided that the LO and received signal powers are properly set.

This paper is organized as follows. In section 2, the architecture of the SiP CRx assembly is presented and the main device-level characterization results are introduced. Section 3, supported by derivations in appendices A and B and the list of mathematical symbols used in Table 1, presents an analytical model for the SNR at the output of the CRx where various noise sources that impact the SNR in a colorless reception scheme are outlined. In section 4, we depict the system-level WDM experimental results obtained using the CRx on 16 × 112 Gbps PDM-QPSK channels where the impact of LO power, received signal power, number of channels and transmission distance on the performance of colorless reception are studied experimentally and compared to theory presented in section 3. Finally, concluding remarks are given in section 5.

Tables Icon

Table 1. List of Symbols and Notation

2. Si-photonic integrated balanced coherent receiver design and characterization

The coherent receiver is based on silicon photonics (SiP) technology which enables a very compact design. Figure 1(a) gives a view of the assembly. Strip waveguides of 220 nm (height) × 500 nm (width) are used to convey the light in the optical chip. A single-mode and a PM fiber whose end facets are polished at an angle close to 45° are used to couple the light from the data-carrying signal and LO respectively into the SiP chip. Polarization diversity is obtained via a 2D surface grating that couples orthogonal linear polarization states from the data-carrying signal in two separate waveguides. For the LO, a 1D grating coupler is used, followed by a 1 × 2 multimode interference (MMI) splitter. The light from each polarization of the incoming signal is combined to the light from the split LO using 2 × 4 MMI couplers. The light from the two 2 × 4 MMI couplers is coupled out of the chip using eight 1D grating coupler and sent to two 1 × 4 photodiode arrays which are flip chip mounted onto the 6 mm × 6 mm SiP chip. The photodiodes are used in balanced pairs and connected to TIAs. The SiP chip and TIAs are fixed onto a ceramic substrate. RF lines connect the AC-coupled differential outputs of the TIAs to eight pins on the CRx package in a GSGSG configuration. For the sake of the present experiments, an RF interface board is used to provide access to a set of four single-ended outputs (XI+, XQ+, YI+ and YQ+) on K connectors. The assembled CRx presents a small-signal bandwidth of 22 GHz and a quadrature error of less than 3°. The effective responsivity at 192.9 THz averaged over all 8 photodiodes (PDs) is 0.025 and 0.029 A/W for the SIG and LO ports, respectively. These effective responsivities include (besides the PD responsivities) the 3 dB intrinsic loss of the polarization beam splitter (PBS) and the 6 dB intrinsic loss of the 2 × 4 MMI coupler in addition to their excess loss as well as coupling losses in and out of the SiP chip. Figure 1(b) shows the overall small signal frequency response including the photodiodes, the TIA and the RF lines on the ceramic. The overshoot in the response is due to the TIA. Finally, Figs. 2(a) and 2(b) show the CMRR versus frequency curves for the LO and SIG ports, respectively measured at the four output ports at 192.9 THz according to the procedure described in [28]. We notice that the CRx provides a DC CMRR better than −14 and −18 dBe for the LO and SIG ports, respectively.

 figure: Fig. 1

Fig. 1 (a) Picture of the PIC of the CRx showing the main buliding blocks, (b) Overall frequency response of the Si-photonic CRx assembly.

Download Full Size | PDF

 figure: Fig. 2

Fig. 2 CMRR versus frequency curves measured at 192.9 THz at the four output RF ports of the CRx when light is applied at the (a) signal port, (b) LO port.

Download Full Size | PDF

3. Theoretical performance evaluation of a colorless coherent receiver

After laying out the CRx chip architecture, the performance of the CRx when employed for colorless coherent reception is evaluated analytically in this section. We present the various noise sources that dictate the performance of colorless CRx operation, highlighting the key CRx specifications that impact the performance. An analytical expression for the SNR at the CRx output is developed. The SNR analytical model derived will be employed later in the manuscript to verify the experimental results presented in forthcoming sections.

Although an analytical SNR expression for the colorless scenario when employing coherent balanced detection has been presented in [16], a more rigorous and complete approach is adopted in this manuscript. Without loss of generality, we consider one of the four CRx output ports {XI,XQ,YI,YQ} and derive an expression of the autocorrelation function (ACF) (and then the power spectral density PSD) of the differential photocurrent after imperfect balanced detection. The desired signal term as well as various noise terms in the final differential photocurrent are identified and their respective ACFs (and PSDs) are derived methodically in terms of input fields to the CRx. Compared to [16], the expressions we obtain hereafter are more complete and the systematic approach adopted throughout the derivation leads to an expression for an “effective CMRR” which is a system parameter that quantifies the ability of balanced detection to reject the SIG-SIG beating, taking into account the frequency dependent nature of the power imbalance and the skew mismatch between the P and N photodiodes of the CRx weighted by the PSD of the received signal intensity. This quantity termed “effective CMRR” was proposed and only arbitrarily defined in [16]. Furthermore, another distinguishing difference from the work in [16] is that the optical noise incoming with the received signal assumed in our model contains the impact of not only the ASE noise from in-line EDFAs but also the induced noise due to fiber nonlinearity where both noise sources can be modeled as AWGN in the context of dispersion unmanaged transmission [2527].

The procedure used to derive the ACFs and PSDs for various terms in the differential photocurrent follows the one outlined in [29]. We start by writing the photocurrents at the outputs of the P and N photodiodes (PDs), denoted by iP(t) and iN(t), in terms of the real electric fields of the signal (all WDM channels and their accompanying optical noise) and LO as

iP(t)=[i=1NchRSIG,iP(ESIG,i(t)+EON,i(t))+RLO,sP(ELO(t))]2+ishP+ithP2RSIG,sPRLO,sP[ESIG,s(t)ELO(t)+EON,s(t)ELO(t)]+i=1NchRSIG,iP[ESIG,i2(t)+EON,i2(t)+2ESIG,i(t)EON,i(t)]+RLO,sPELO2(t)+ishP+ithP,
iN(t)=[i=1NchRSIG,iN(ESIG,i(tδ)+EON,i(tδ))RLO,sN(ELO(tδ))]2+ishN+ithN2RSIG,sNRLO,sN[ESIG,s(tδ)ELO(tδ)+EON,s(tδ)ELO(tδ)]+i=1NchRSIG,iN[ESIG,i2(tδ)+EON,i2(tδ)+2ESIG,i(tδ)EON,i(tδ)]+RLO,sNELO2(tδ)+ishN+ithN,
where ESIG,i, EON,i, ELO are the real-valued electric fields of the ith signal out of Nch WDM signals presented to the CRx in colorless scheme, optical noise accompanying the ith WDM signal, and local oscillator (LO) all measured at the input to the CRx, respectively. Without loss of generality, we assume that channel s{1,2,...,Nch} is the desired channel in the WDM spectrum. Also, we assume that the LO frequency is tuned such that the frequency offset between the LO and the center frequency of channel s is much less than its baseband bandwidth, i.e. intradyne reception is used. Thus, the LO beating with OOB channels will lie outside the CRx bandwidth and hence was eliminated in the second equality of Eqs. (1) and (2). Also, we assume that EON,i is the total incoming optical noise with the ith WDM channel which contains both the ASE noise from in-line optical amplifiers and the noise induced due to fiber NL. According to [2527], the NL-induced noise in case of dispersion uncompensated transmission can be modeled accurately as AWGN similar to ASE noise and hence, both additive noises can be combined in one total optical noise field which finally results in a worse optical signal-to-noise ratio (OSNR) governing the system performance [27]. According to [27], this modified OSNR at the CRx input is defined as OSNR=PSIG/PON where PON=PASE+PNL is the total optical noise power due to both ASE and NL contributions. As will be shown in subsection 4.2.3, both the ASE noise and the NL-induced noise scale linearly with transmission distance (assuming the NL-induced noise contributions from multiple spans are approximately summed incoherently [25,27]) and will be indistinguishable in our experimental data since the launch power was kept constant throughout the measured data. In addition, RSIG,iP,RSIG,iN are the effective responsivities at the wavelength of the ith channel measured at the P and N PDs to the optical power applied to the signal port. In addition, RLO,sP,RLO,sN are the effective responsivities measured at the P and N PDs to the LO power applied at the wavelength of channel s. It is also noteworthy that the effective responsivities defined above take into account (besides the photodiode responsivities) the 3 dB intrinsic loss of the polarization beam splitter (PBS) and the 6 dB intrinsic loss of the 2 × 4 MMI coupler in addition to their excess loss as well as coupling losses in and out of the SiP chip. In addition, δis the time skew between the P and N paths measured from the output of the MMI couplers up to current subtraction which occurs inside the TIA. Finally, ishP,ishN are the shot noise currents at the P and N photodiode outputs and ithP,ithN are the input-referred thermal noise currents at the P and N inputs of the TIA, respectively. In Eqs. (1) and (2), we also neglect the SIG-SIG beating between each WDM channel with its neighbors and consider only the beating between each signal with itself which lies within the CRx bandwidth. Although a portion of the spectral components from the beating between two adjacent WDM channels may fall inside the baseband CRx bandwidth, this portion is negligible compared to the SIG-SIG beating of each channel with itself especially after being lowpass filtered by the CRx frequency response.

Similar to [29], we express the real-valued fields in Eqs. (1) and (2) in terms of baseband complex envelops (denoted by boldface symbols throughout the manuscript) as

ESIG,i(t)=Re{E˜SIG,i(t)ej2πνSIG,it},ELO(t)=Re{E˜LO(t)ej2πνLOt},EON,i(t)=Re{E˜ON,i(t)ej2πνSIG,it},
where E˜k(t) is the baseband complex envelope representing the real bandpass field Ek(t) and ej2πνktis the complex-valued carrier at the optical frequency υk. Now, one can use the relation
E1(t)E2(t)=Re{E˜1(t)ej2πν1t}Re{E˜2(t)ej2πν2t}=12Re{E˜1(t)E˜2*(t)ej2π(ν1ν2)t}+12Re{E˜1(t)E˜2(t)ej2π(ν1+ν2)t},
to re-write Eqs. (1) and (2) noting that the terms at the sum of optical frequencies will be filtered out, yielding
iP(t)=2RSIG,sPRLO,sP12Re{E˜SIG,s(t)E˜LO*(t)ej2πΔνt+E˜ON,s(t)E˜LO*(t)ej2πΔνt}+12i=1NchRSIG,iP[|E˜SIG,i(t)|2+|E˜ON,i(t)|2+2Re{E˜SIG,i(t)E˜ON,i*(t)}]+RLO,sP12|E˜LO(t)|2+ishP+ithP,
iN(t)=2RSIG,sNRLO,sN12Re{E˜SIG,s(tδ)E˜LO*(tδ)ej2πΔν(tδ)+E˜ON,s(tδ)E˜LO*(tδ)ej2πΔν(tδ)}+12i=1NchRSIG,iN[|E˜SIG,i(tδ)|2+|E˜ON,i(tδ)|2+2Re{E˜SIG,i(tδ)E˜ON,i*(tδ)}]+RLO,sN12|E˜LO(tδ)|2+ishN+ithN,
where Δν=νSIG,sνLO is the frequency offset between the center frequency of channel s and the LO. Next, the difference between the P and N photocurrents in Eqs. (5) and (6), which represents the final differential photocurrent, can be written as
Δi(t)=iP(t)hP(t)iN(t)hN(t)=ΔiLOSIG(t)+ΔiLOON(t)+ΔiSIGSIG(t)+ΔiONON(t)+ΔiSIGON(t)+ΔiLOLO(t)+Δish(t)+Δith(t),
where denotes convolution, hP(t) and hN(t) are the impulse responses of the P and N paths from the PD outputs up to current subtraction in the TIAs, i.e. they include the impact of limited bandwidth of PDs, TIA and wire bonding between them. Ideally, hP(t) and hN(t) should be the same but the slight difference between the frequency responses should be taken into account since it will impact the frequency dependence of the common mode rejection ratio of the CRx. In Eq. (7), we separate the individual contributions in the final photocurrent where the first term is the LO-SIG beating term which is the desired signal term resulting from subtracting the first terms of Eqs. (5) and (6) after convolving each of them with hP(t) and hN(t). The next terms in Eq. (7) represent all noise contributions in the final differential photocurrent due to LO-ON, SIG-SIG, ON-ON, SIG-ON, LO-LO beatings as well as shot and thermal noises which also result from subtracting all terms subsequent to the LO-SIG term in Eqs. (5) and (6) after convolving with the P and N impulse responses.

In case of ideal balanced detection, i.e. RSIG,iP=RSIG,iN, RLO,sP=RLO,sN, hP(t)=hN(t) and δ=0, the SIG-SIG, ON-ON, SIG-ON and LO-LO noise terms in Eq. (7) cancel out since these terms are common to both the P and N PDs and will be entirely rejected. In this case, the LO-ON beating as well as the thermal and shot noises become the only noise sources contributing to the final differential photocurrent. However in reality, imperfect balanced detection due to frequency dependent power imbalance and skew mismatch between the P and N photocurrents will lead to residual noise power due to imperfect common mode rejection. For the rest of the derivation, we assume imperfect balanced detection and derive the time-averaged ACF of Δi(t) in Eq. (7), which is denoted by Γ¯Δi(τ), as follows

Γ¯Δi(τ)=Δi(t)Δi(tτ)¯=Γ¯ΔiLOSIG(τ)+Γ¯ΔiLOON(τ)+Γ¯ΔiSIGSIG(τ)+Γ¯ΔiONON(τ)+Γ¯ΔiSIGON(τ)+Γ¯ΔiLOLO(τ)+Γ¯Δish(τ)+Γ¯Δith(τ),
where denotes ensemble (statistical) averaging and the bar denotes time averaging. In deriving Eq. (8), all the cross correlation terms between currents from various sources can be shown to cancel out by the use of statistical independence between different fields as well as the fact that E˜(t)=0 for all fields. This finally means that the ACF of the final photocurrent is the sum of the ACFs of all individual terms in Eq. (7). The next step is to find the individual ACFs, i.e. each term in Eq. (8), in terms of the input fields as well as the CRx parameters. After that, we make use of the Wiener–Khintchine theorem for non-stationary signals [30] which states that the power spectral density (PSD) of a non-stationary signal is the Fourier transform (FT) of its time-averaged ACF Γ¯E(τ) which can be applied to Eq. (8), yielding
SΔi(f)=FT {Γ¯Δi(τ)}=SΔiLOSIG(f)+SΔiLOON(f)+SΔiSIGSIG(f)+SΔiONON(f)+SΔiSIGON(f)+SΔiLOLO(f)+SΔish(f)+SΔith(f)SΔiLOSIG(f)+SΔiLOON(f)+SΔiSIGSIG(f)+SΔish(f)+SΔith(f),
where in the last equality, we neglect the ON-ON, SIG-ON and LO-LO contributions to the overall photocurrent difference [31]. The argument behind neglecting those three terms is that they were found to be experimentally negligible compared to the SIG-SIG beating. In fact, this experimental observation is intuitive since the SIG-SIG, ON-ON and SIG-ON beatings are proportional toPSIG2, PON2 and PSIGPON, respectively, where PSIG,PON are the average optical powers of the received signal and the total optical noise real-valued fields within the bandwidth of one channel. Given that the OSNR level (defined similar to [27] as the ratio between the optical signal power to the total optical noise power including both ASE and NL-induced noises) at the CRx input at which error-free operation is achieved is relatively high (e.g. > 12 dB in 0.1 nm for 28 Gbaud PDM-QPSK), we can notice that ON-ON and SIG-ON beatings are negligible compared to SIG-SIG beating. Regarding the LO-LO beating, SΔiLOLO(f) is proportional to PLO2RIN(f), where PLO is the average optical power of the LO field and RIN(f) is the PSD of the relative intensity noise of the LO laser. Since the LO used in experiments had a RIN of around −145 dB/Hz and the maximum LO power used throughout experiments is 15.5 dBm, the LO RIN contribution to the final photocurrent was experimentally found to be negligible.

Using Eq. (9), we can evaluate the SNR achieved after offline DSP if we assume that the receiver side DSP will mitigate transmission impairments by employing various algorithms, some of which involving filtering of the final sampled photocurrent. If we assume only linear equalization within the DSP, i.e. no NL mitigation algorithms, and lump all filtering performed in DSP in the transfer function HrxDSP(f), one can calculate the final SNR after DSP as the total power of the useful LO-SIG term after DSP filtering divided by the total AC power from various noise sources after filtering, which yields

SNR=SΔiLOSIG(f)|HrxDSP(f)|2df[SΔiLOON(f)+SΔiSIGSIG(f)+SΔish(f)+SΔith(f)]AC|HrxDSP(f)|2df=PΔiLOSIGσΔiLOON2+σΔiSIGSIG2+σΔish2+σΔith2,
where []AC denotes keeping only the AC portion of each PSD involved in the integration. To further pursue the derivation, we should evaluate the PSDs involved in the integrations in Eq. (10) in terms of the input fields to the CRx as well as the CRx parameters. The rigorous way to perform such a task is to find the time-averaged ACFs of each term from which the Fourier transformation leads to the PSDs in Eq. (10). In appendix A, we provide the complete derivation of all ACFs and PSDs from which the powers required to evaluate the SNR in Eq. (10) are evaluated. If we substitute Eqs. (29), (32), (41), (42) and (43) derived in appendix A into Eq. (10), we reach the following final analytical expression for the SNR of one of the four electrical photocurrents
SNR=2PLOPSIG2c1PLOPON+c2(σΔish2+σΔith2)+CMRR¯SIGNchβPSIG2,
where c1, c2 and CMRR¯SIG are obtained as follows
c1=Bch/2Bch/2|HP(f)RSIG,sPRLO,sP+HN(f)ej2πfδRSIG,sNRLO,sN|2|HrxDSP(f)|2df/Bch|HP(f)RSIG,sPRLO,sP+HN(f)ej2πfδRSIG,sNRLO,sN|2|G(f)|2|HrxDSP(f)|2df/|G(f)|2df,
c2=1|HP(f)RSIG,sPRLO,sP+HN(f)ej2πfδRSIG,sNRLO,sN|2|G(f)|2|HrxDSP(f)|2df/|G(f)|2df,
CMRR¯SIG=|HP(f)RSIGPHN(f)ej2πfδRSIGN|2S|E˜SIG|2AC(f)|HrxDSP(f)|2df/S|E˜SIG|2AC(f)|Hav(f)|2|HrxDSP(f)|2df|HP(f)RSIGPRLOP+HN(f)ej2πfδRSIGNRLON|2|G(f)|2|HrxDSP(f)|2df/|G(f)|2df,
and σΔish2,σΔith2 are the variances (AC power) of the shot noise and thermal noise contributions in the differential photocurrent which are expressed by Eqs. (42) and (43) in appendix A.

Also, the scaling factor β is defined in appendix A as the ratio between the AC power from the filtered photocurrent generated by the field intensity of one channel and its DC electrical power. In appendix B, a more thorough study of β is performed to assess its dependence on the residual chromatic dispersion (CD) and polarization orientation of the received signal. It should be noted that Eq. (11) was proved by considering only one of the four outputs provided by the CRx; however, the derived SNR expression has the same form for all ports except for using the parameters of each port, e.g. responsivities and frequency responses. Hence, an overall average SNR can be calculated by averaging the SNRs calculated at the four ports.

The quantity CMRR¯SIG denotes the so-called “effective CMRR” of the SIG-SIG beating which is a system parameter (not a CRx parameter) that quantifies the ability of imperfect balanced detection to reject the common-mode SIG-SIG beating while taking into account the CRx parameters, namely, the effective responsivities, frequency responses and time skew for the P and N paths up to current subtraction inside the TIAs, as well as the system parameters, namely, the residual CD, polarization orientation and pulse shape which are imbedded in the PSD of the intensity S|E˜SIG|2AC(f) (see Eq. (40) in appendix A). In fact, the effective CMRR definition given above contains the quantity |HP(f)RSIGPHN(f)ej2πfδRSIGN|2 which quantifies the CRx ability to reject the common mode SIG-SIG beating incoming from the SIG port. Since this rejection ability is frequency dependent, the effective CMRR expression also weighs this quantity using the PSD of the common mode quantity being rejected, i.e. the field intensity, as represented by S|E˜SIG|2AC(f) in the CMRR¯SIGexpression. Dividing by the denominator in the effective CMRR expression means that CMRR¯SIG is a measure of how strong the power of the residual common mode SIG-SIG beating is compared to that of the differential LO-SIG beating. This explains why the denominator is proportional to the product of responsivities of the SIG and LO ports whereas the numerator is proportional to the square of the responsivity of the SIG port. Finally, we should note that the above expression for the effective CMRR was obtained in a systematic and rigorous manner in the derivation included in appendix A. This differs from [16] where the authors proposed the interesting concept of an effective CMRR, being a single averaged value applied to the common mode beating and defined quite arbitrarily as the CMRR value at half the Nyquist frequency.

4. WDM experiment, results and discussion

In this section, we evaluate experimentally the system-level performance of the SiP CRx module whose chip architecture was described in section 2. First, we introduce the experimental setup used to assess the CRx performance using 16 × 112 Gbps PDM-QPSK WDM channels. Subsequently, we present the experimental results, starting with the BER for all 16 channels at different distances and received signal powers. Afterwards, we present a series of extensive measurements obtained for one of the center WDM channels where the LO power, received signal power, distance, and number of channels presented to the CRx are swept. A thorough discussion of the results is given, focusing on the major noise sources affecting the performance of the colorless operation of the CRx for different operating parameters in view of the analytical SNR expression derived in section 3. Finally, we use the experimental SNR data to fit the analytical SNR expression derived in section 3. The fitted SNR model is then used to predict the SNR at other operating parameters beyond our experimental capabilities, e.g. larger number of WDM channels.

4.1. Experimental setup

Figure 3 shows the experimental setup used in various WDM measurements presented in the following subsections. CW light from 16 distributed feedback (DFB) lasers (50 GHz spaced and centered at 192.9 THz) are multiplexed using a 50 GHz arrayed waveguide grating WDM multiplexer and then bulk modulated by means of a QAM transmitter. The QAM transmitter comprises a single polarization LiNbO3 double-nested Mach-Zehnder modulator and two high-speed RF amplifiers to drive the two quadratures. Binary I and Q electrical signals are fed to the QAM transmitter from two MICRAM DACs driven by two Xilinx FPGAs at 28 Gbaud to realize 16 × 112 Gbps PDM-QPSK WDM channels on a 50 GHz grid. Then, a polarization beam splitter (PBS), variable optical delay line (VODL) and polarization beam combiner (PBC) are employed for PDM emulation with a decorrelation delay between the two orthogonal polarizations equivalent to 700 symbols at 28 Gbaud. After optical amplification, odd and even WDM channels are decorrelated by ~7.5 ns using a 50 GHz interleaver / de-interleaver. The resulting optical signal was attenuated using a variable optical attenuator (VOA) to set the launch power at −2 dBm per channel which was experimentally found to be the optimum launch power that achieves the right balance between fiber nonlinearity and ASE noise from in-line EDFAs. At this launch power, which was kept constant throughout the experiment, both ASE and NL-induced noises are non-negligible and hence, they were both included in the SNR model presented in section 3. Next, the resulting signal was launched into a recirculating loop containing four spans of 80 km of SMF-28e+ fiber, each span being followed by an erbium-doped fiber amplifier (EDFA) having a noise figure of 5 dB. In addition, an optical WaveShaper from Finisar is inserted after the second in-line EDFA for gain flattening.

 figure: Fig. 3

Fig. 3 Experimental setup for WDM experiments, (IL: Interleaver, DFB laser: Distributed Feedback laser, AWG: Arrayed Waveguide Grating, PBS: Polarization Beam Splitter, PBC: Polarization Beam Combiner, VODL: Variable Optical Delay Line, ECL: External Cavity Laser, OSA: Optical Spectrum Analyzer, VOA: Variable Optical Attenuator).

Download Full Size | PDF

At the output of the loop, no noise loading is performed hence the OSNR is only dictated by the ASE noise from in-line EDFAs around the loop as well as the fiber NL-induced noise. We used a 99:1 coupler to monitor the WDM spectrum on an optical spectrum analyzer (OSA) and to measure the OSNR. After 1600, 3200 and 4800 km of transmission, i.e. 5, 10 and 15 loops respectively, the measured OSNRs, assuming a 0.1 nm reference bandwidth, were 21.4, 18.9 and 17.6 dB, respectively. After the 99:1 coupler, we distinguish between two reception configurations: 1) a fully colorless configuration in which the attenuated signal carrying all 16 channels is fed directly to the signal port of the CRx module at a total power level PSIG,tot ( = NchPSIG) controlled by the VOA, i.e. elements in the grey shaded box are skipped, 2) another configuration in which the WDM signal is first amplified to a constant total power of 22 dBm, subsequently filtered using a wavelength & bandwidth tunable filter, and then attenuated using a VOA to set the total received power level (PSIG,tot) before finally entering the CRx, i.e. elements in the grey shaded box are employed. Inserting the EDFA prior to the CRx allows us to sweep PSIG,tot, and subsequently the received power per channel (PSIG), to levels larger than the received power from the loop. As will be observed in the next subsection, the sweep range of PSIG extends up to 5 dBm per channel in order to evaluate the impact of colorless reception on the dynamic range of the CRx. After the EDFA, the wavelength & bandwidth tunable filter had its center frequency constantly set to that of channel 9 (with channel 1 having the smallest frequency) and its bandwidth swept to vary the number of channels presented to the CRx denoted by Nch. Using this filter, we can vary Nch between 1 and 16, where the extreme cases correspond to fully filtered and fully colorless reception schemes, respectively.

In reception configuration (1), the frequency of the LO laser is swept to the center frequency of the channel desired to be selected out of the 16 available channels. However in configuration (2), the LO frequency is constantly set to channel 9 (channel of interest in this configuration) and the bandwidth of the tunable filter is set to control the number of out-of-band (OOB) WDM neighbors (Nch - 1) co-incident on the CRx with channel 9. In both configurations, the VOA sets the received power level per channel (PSIG) as desired prior to the CRx. At all received power levels, the voltage swing provided by the TIAs at the output of the CRx is maintained at a constant 320 mVpp (single-ended). Finally, the four electrical signals are sampled at 80 GSa/s using one Agilent DSOX96204Q real-time scope and stored for offline processing using a myriad of standard DSP algorithms to mitigate various impairments.

Offline processing starts by IQ power imbalance compensation and quadrature error correction using Gram-Schmidt orthogonalization, resampling to 2 samples per symbol, chromatic dispersion (CD) compensation in the frequency domain [1,5,6] and frequency offset removal using the periodogram method [32]. Polarization is then demultiplexed using a 31-tap butterfly equalizer whose taps are updated using the constant modulus algorithm (CMA) with a step size of 1 × 10−3 [1]. Phase noise was mitigated by a decision-directed phase locked loop (DD-PLL) [33]. After hard decision on a symbol-by-symbol basis, the achieved electrical SNR was evaluated and BER was finally calculated.

4.2 Experimental results and discussion

In this subsection, we present all experimental results obtained using the setup depicted in Fig. 3. In subsection 4.2.1, we show the results of reception configuration (1) in which fully colorless and preamplifierless strategy is employed. Next, subsection 4.2.2 presents all results collected on channel 9 employing reception configuration (2) where the EDFA and the tunable wavelength & bandwidth filter in the grey shaded box in Fig. 3 are inserted prior to the CRx. Using this configuration, various operating parameters such as the received signal power per channel (PSIG), LO power (PLO), number of channels (Nch) incident on the CRx, and distance (L) are swept. Results collected from channel 9 include the SNR measured after offline DSP while sweeping the aforementioned parameters to evaluate the impact of each parameter on the performance of colorless CRx operation and to identify the dynamic range of the CRx for various number of channels. Based on the results, we also identify the key specifications of the CRx that are crucial for colorless operation. In subsection 4.2.3, we use the experimental SNR data at various operating parameters to fit the analytical SNR expression derived in section 3. The fitted SNR model is then used to predict the SNR at other operating parameters beyond our experimental capabilities, e.g. larger Nch.

4.2.1 Configuration (1): BER for all channels at various distances using fully colorless and preamplifierless reception at either −3 dBm/ch or −21 dBm/ch

BERs of all 16 WDM channels were calculated and are plotted in Fig. 4 for 3200, 4480 and 4800 km transmission at two different PSIG (−3 and −21 dBm per channel). All results in Fig. 4 are collected in the fully colorless and preamplifierless reception termed above by reception configuration (1). The LO power (PLO) was set to 15.5 dBm and its frequency was tuned to select each one of the 16 channels before data was acquired. Insets at the right of Fig. 4 also show the WDM spectra in back-to-back and after 3200 km transmission obtained using a 0.05 nm resolution bandwidth on the OSA. Throughout the manuscript, we assume 3.8 × 10−3 forward error correction (FEC) threshold, which is typical for the hard-decision FEC algorithm with 7% overhead defined in [34]. Then, if we carefully examine the BERs of all 16 channels, we can claim a 4800 km reach for all channels at PSIG of −3 dBm which slightly decreases to 4160 km at −21 dBm per channel. As will be explained in upcoming subsections, the good performance (reach) achieved by the CRx for fully colorless and preamplfierless reception with very low received signal power of −21 dBm per channel is enabled by moderate CMRR and responsivity possessed by the CRx.

 figure: Fig. 4

Fig. 4 BER for all WDM channels using fully colorless and preamplifierless reception at different transmission distances and at per channel received signal power of −3 and −21 dBm. Insets show the received WDM spectra in back-to-back and after 3200 km transmission where we observe a fairly flat spectrum due to the gain flattening waveshaper in the loop.

Download Full Size | PDF

4.2.2 Configuration (2): SNR of channel 9 when using a varying number of channels

Having evaluated the CRx colorless / preampliferless operation for each channel across the spectrum, we now focus on channel 9 specifically. All measurements on channel 9 presented in the upcoming subsections, which involve sweeping various experimental parameters such as PLO, PSIG, Nch and L, will use the SNR calculated after offline DSP instead of BER as the performance metric. The reason for using SNRs instead of BERs is that SNRs allow us to assess the performance quantitatively in back-to-back and short transmission distance scenarios when BERs are usually below 10−5. Evaluating BERs in this low-error regime would require acquiring and storing a number of symbols exceeding the memory of the real-time scope used. Also, later in this section, utilizing the experimental SNRs allows us to perform direct comparison with analytically predicted SNRs based on the derived SNR expression in section 3.

Figure 5 shows the relationship between experimentally evaluated BER and SNR for channel 9. In Fig. 5(a), we plot the calculated BER on a double log scale versus the transmission distance whereas Fig. 5(b) shows the corresponding SNR evaluated at these distances. Using both curves, we plot in Fig. 5(c) the BER versus SNR relationship found experimentally. We find that a 9.2 dB SNR corresponds to the 3.8 × 10−3 FEC threshold which is 0.67 dB higher than the theoretical SNR of 8.53 dB at this BER level for 28 Gbaud QPSK assuming an ideal additive white Gaussian noise (AWGN) channel [35]. By means of Fig. 5(c), all SNR results that will be presented in the forthcoming subsections can be mapped into BER values in a straightforward manner.

 figure: Fig. 5

Fig. 5 Measurement results of channel 9: (a) BER versus transmission distance, (b) SNR evaluated after offline DSP versus distance, (c) BER versus SNR using both experimental results and theory where the experimental SNR corresponding to FEC threshold is marked.

Download Full Size | PDF

4.2.2.1 Impact of received signal power and number of channels presented to the CRx on the SNR of channel 9 at various distances and LO powers

From this point forward, the entire measurements are performed only for channel 9 in the WDM spectrum. In Figs. 6(a)6(f), we show the SNR versus PSIG at different LO powers of 15.5, 12, 9, 6, 3 and 0 dBm, respectively. In each subplot, we plot the SNR curves with different colors corresponding to various number of channels (Nch) presented to the CRx where Nch is set to either 1, 5, 11 or 16 channels, ranging from fully filtered to colorless reception respectively. Also in each subplot, SNR curves measured at various transmission distances of 0, 1600, 3200 and 4800 km are plotted together where the top group of four curves always corresponds to the back-to-back case whereas the bottom set of curves represents the 4800 km transmission case. By examining Figs. 6(a)6(f), we notice that the sweep range of PSIG starts at −24 dBm for PLO = 15.5 and 12 dBm (see Figs. 6(a) and 6(b)), whereas for PLO = 9, 6, 3 and 0 dBm the sweep range of PSIG had to be reduced to start at −21, −18, −15 and −12 dBm, respectively (see Figs. 6(c)6(f))). The reason for increasing the lower end of the PSIG sweep range for PLO ≤ 9 dBm is to be able to achieve the fixed target voltage swing of 320 mVpp after the TIAs at the CRx output. For that reason, at each PLO ≤ 9 dBm, the minimal value of the PSIG sweep range was fixed to the least power level at which the TIAs could provide the 320 mVpp target swing when their respective gain was set to maximum (5 kΩ differential transimpedance Zt). Since for a fixed target voltage swing, the required Zt for the TIAs is proportional to PLOPSIG, we notice in Figs. 6(c)6(f) that as PLO is decreased by 3 dB going from one subfigure to the next, the starting point of the PSIG sweep range (at maximum Zt) increases by 3 dB supporting the above proportionality.

 figure: Fig. 6

Fig. 6 SNR of channel 9 versus received signal power per channel PSIG after 0, 1600, 3200 and 4800 km transmission with varying number of channels presented to the CRx at an LO power of: (a) 15.5 dBm, (b) 12 dBm, (c) 9 dBm, (d) 6 dBm, (e) 3 dBm, (f) 0 dBm. In each subfigure, different colors represent different number of channels Nch according to the legend. Also, SNR evaluated after 0, 1600, 3200 and 4800 km transmission are shown in each subfigure where the top set of four curves in each subfigure is evaluated in back-to-back and the bottom set is after 4800 km transmission.

Download Full Size | PDF

Looking at the SNR variation with PSIG throughout Figs. 6(a)6(f), we notice that regardless of PLO, L and Nch, the SNR improves as PSIG is increased to an extent after which the SNR starts to decrease if PSIG is further increased. This means that an optimum PSIG exists for each PLO, L, and Nch. Furthermore, if a certain target required SNR at a certain PLO, L and Nch is to be achieved, one can specify a dynamic range for PSIG within which the target SNR is at least guaranteed. It is noticeable that the dynamic range of PSIG at a given L and PLO decreases as Nch is increased from 1 to 16 showing the impact of colorless reception. Finally, the dynamic range of PSIG at a given L and Nch is also found to decrease when PLO is decreased.

The two reasons for the SNR degradation at the low and high ends of PSIG are distinct. At high received signal power (PSIG ~0 dBm), the SNR degradation observed with a further increase of PSIG is due to the residual SIG-SIG beat noise which is expressed by the last term in the denominator of Eq. (11). This SIG-SIG beat noise is proportional to NchPSIG2 which explains why the SNR degrades as PSIG is increased and why this degradation is more pronounced as more channels are presented to the CRx. In fact, we can re-write Eq. (11) by neglecting the thermal and shot noises at high PSIG and dividing by PLOPSIG, which yields

SNR=22c1OSNR1+CMRRSIG¯NchβLSR1,

where LSR is the LO to signal power ratio at the CRx input. We can see from Eq. (12) that for a fixed OSNR, i.e. fixed distance, the SNR depends on the effective CMRR defined in section 3, Nch and LSR. If more channels are introduced to the CRx, the SNR degrades due to larger SIG-SIG beating from the OOB neighbors. Also, the impact of the SIG-SIG beat noise can be greatly reduced if high enough LSR is employed, which is also apparent in Fig. 6 by observing the impact of PLO on the SNR difference between the cases of Nch = 1 and Nch = 16, i.e. the SNR penalty due to colorless reception compared to filtered reception denoted by ΔSNR1,16, at a given L and at PSIG = 0 dBm.

Since we are now interested in the high end of PSIG and studying the impact of Nch on SNR, we re-plot in Fig. 7 a portion of the results shown in Fig. 6 by focusing only on the SNR at PSIG = 0 dBm. Figures 7(a)7(f) show the SNR after 0, 1600, 3200 and 4800 km transmission versus Nch at LO powers from 15.5 dBm down to 0 dBm, all at PSIG = 0 dBm. This exactly corresponds to LSRs ranging from 15.5 dB to 0 dB. By examining Fig. 7 carefully, we notice that for LSR ≥ 9 dB, i.e. PLO ≥ 9 dBm, ΔSNR1,16 in Figs. 7(a)7(c) is found to be < 1.2 dB at any of the various distances travelled, which is a reasonably small degradation given the transceiver flexibility and hardware savings gained in an entirely colorless reception scheme. In fact, this 1.2 dB SNR penalty due to colorless reception is only pronounced in back-to-back at PLO = 9 dBm whereas ΔSNR1,16 is less after transmission. This can be explained in light of the approximate SNR expression in Eq. (12) by noting that the OSNR decreases as L increases and hence the SIG-SIG beat term becomes less significant relative to the first term. On the other hand for LSR < 9 dB, ΔSNR1,16 in Figs. 7(d)7(f) increases significantly where it becomes most apparent in Fig. 7(f) when LSR = 0 dB. In view of Eq. (12), it should be finally noted that the SNR degradation due to colorless reception can be greatly reduced by developing coherent receivers with improved CMRRs which results in a better rejection for the SIG-SIG beating from the OOB channels.

 figure: Fig. 7

Fig. 7 SNR of channel 9 versus Nch at PSIG = 0 dBm after 0, 1600, 3200 and 4800 km transmission at PLO of (a) 15.5 dBm, (b) 12 dBm, (c) 9 dBm, (d) 6 dBm, (e) 3 dBm, (f) 0 dBm.

Download Full Size | PDF

Looking back at Fig. 6 for low received signal powers (PSIG ≤ −12 dBm), we notice that the SNR degrades quickly as PSIG decreases and that the SNR is independent of Nch. Referring back to the general SNR expression in Eq. (11), we re-write the expression by neglecting the SIG-SIG beat term for low PSIG and dividing by PSIGPLO, which yields

SNR=22c1OSNR1+c24eRLOavΔfPSIG+c2iTIA2ΔfPLOPSIG,
where we neglected the signal contribution to the shot noise and kept only the LO shot noise term. In view of Eq. (13), we can easily explain the SNR trends at the low end of PSIG in Fig. 6. First, it is clear that as PSIG increases while PLO is fixed, the SNR increases due to the second and third terms (shot and thermal noises) in the denominator of Eq. (13) being smaller which matches our observation in Fig. 6. In addition, we also observe in Fig. 6 that the impact of changing PSIG is less pronounced as L increases which can be explained by the decrease in the OSNR with L, thus making the first term in the denominator of Eq. (13) larger compared to the second and third ones. Our next observation based on Eq. (13) is that increasing PLO for a fixed PSIG should also improve the SNR which goes along with the experimental results in Fig. 6. Comparing the impact of PSIG and PLO on the SNR in the low received signal power regime, we can observe in Fig. 6 that increasing PSIG by 3 dB for a fixed PLO results in a larger SNR improvement compared to increasing PLO by 3 dB for a fixed PSIG. As an example, if PSIG is increased from −21 to - 18 dBm at PLO = 9 dBm, the SNR improves by 1.68 dB whereas the improvement is 1.26 dB if PLO increases from 9 to 12 dBm while keeping PSIG at −21 dBm. The observation that PSIG has a slightly larger impact on the SNR compared to PLO can also be explained in light of the denominator of Eq. (13) by noting that PSIG affects both the second and third terms compared to PLO impacting only one term. Next, we notice that the SNR in Eq. (13) is totally independent of Nch at low PSIG which agrees with Fig. 6. In order to further make this independence clear, we re-plot in Fig. 8 a portion of the results in Fig. 6 which are collected at PSIG = −12 dBm. Figures 8(a)8(f) show the SNR after 0, 1600, 3200 and 4800 km transmission versus Nch at LO powers from 15.5 dBm down to 0 dBm, all at low PSIG of −12 dBm. In Fig. 8, it is evident that Nch has no impact on the SNR at either any distance or any PLO at PSIG = −12 dBm.

 figure: Fig. 8

Fig. 8 SNR of channel 9 versus Nch at PSIG = −12 dBm after 0, 1600, 3200 and 4800 km transmission at PLO of (a) 15.5 dBm, (b) 12 dBm, (c) 9 dBm, (d) 6 dBm, (e) 3 dBm, (f) 0 dBm.

Download Full Size | PDF

Finally, it should be also noted that the CRx responsivity is a key specification in this low PSIG regime because according to Eq. (13) if all parameters stay the same while RSIG and RLO are increased, the SNR will consequently increase (note that c2(RSIGRLO)1). Physically speaking, larger responsivities translate into larger photocurrents from the PDs relative to the fundamental CRx noise. Although the CRx used in our experiments possess a moderate responsivity due to relatively high coupling losses into the SiP chip, it still allows for good performance at low PSIG. For example when PLO = 12 dBm, the SNR after 1600-4800 km transmission degrades by < 1 dB when PSIG decreases from 0 to −18 dBm while providing the same output voltage swing using the TIA gain.

4.2.2.2 Impact of LO power on the SNR of channel 9 at various distances and number of channels presented to the CRx for both regimes of high and low received signal power

In this subsection, we explicitly show the impact of the LO power on the SNR of channel 9. In Fig. 9, we re-plot a portion of the results shown in Fig. 6 where we show the SNR versus PLO at various distances and number of channels presented to the CRx for both regimes of high and low received signal powers (PSIG = 0 and −15 dBm, respectively). For the PSIG = 0 dBm case in Fig. 9(a), we observe how significant it is to choose a high LSR which greatly reduces the SNR penalty due to colorless reception. For example, changing the LSR from 0 to 12 dB reduces ΔSNR1,16 from 1.5 to 4.7 dB (for distances travelled from 4800 to 0 km) to 0.1-0.65 dB. This agrees with Eq. (12) as discussed in the previous subsection. Moving to Fig. 9(b) when PSIG = −15 dBm, we notice that increasing PLO improves the SNR and that the SNR is independent of Nch which agrees with Eq. (13). As mentioned in section 3, the LO-LO beat noise was not observed throughout all experimental data because the RIN of our lasers is below −145 dB/Hz. However, it is noteworthy that as PLO is increased beyond 15.5 dBm (maximum output power from our ECLs) at PSIG = −15 dBm (see Fig. 9(b)), the LO-LO beat noise should eventually lead to SNR degradation.

 figure: Fig. 9

Fig. 9 SNR of channel 9 versus PLO after 0, 1600, 3200 and 4800 km transmission at various Nch and at PSIG of: (a) 0 dBm, (b) −15 dBm.

Download Full Size | PDF

4.2.2.3 Two dimensional (2-D) plots of SNR of channel 9 versus LO power and received signal power at different distances in case of 1 (filtered) and 16 channels (colorless)

In this subsection, we re-present the experimental data of Fig. 6 corresponding to the extreme cases of Nch = 1, i.e. filtered operation and Nch = 16, i.e. fully colorless operation. Here, we present the SNR data as two dimensional (2-D) plots versus both PSIG and PLO as the two independent variables being swept. Figure 10 depicts these 2-D plots where different SNR values are encoded into different colors according to the color bars beside each subfigure (red representing higher SNR whereas blue representing smaller SNR). In all Figs. 10(a)10(h), the small x’s mark the positions of experimental data points across the 2-D space (PSIG, PLO). All values between those points are evaluated using linear interpolation using the two nearest experimental data points resulting in continuous and smooth color-coded surfaces. The four rows of subfigures in Fig. 10 represent data after 0, 1600, 3200, and 4800 km transmission, respectively, whereas the left column of subfigures represent the filtered case (Nch = 1) and the right column represent the colorless case (Nch = 16). We should also notice that each pair of subfigures in any row share the same SNR color bar to ease the comparison within each row between the filtered and colorless scenarios. Also, we observe in all subfigures a white space at the bottom left corner where no surface is plotted. This range of (PSIG, PLO) has no experimental data since the target voltage swing of 320 mVpp could not be achieved within this range when both PSIG and PLO are small.

 figure: Fig. 10

Fig. 10 Color coded 2-D plots of SNR versus both PSIG and PLO for (a) back-to-back and Nch = 1, (b) back-to-back and Nch = 16, (c) 1600 km transmission and Nch = 1, (d) 1600 km transmission and Nch = 16, (e) 3200 km transmission and Nch = 1, (f) 3200 km transmission and Nch = 16, (g) 4800 km transmission and Nch = 1, and (h) 4800 km transmission and Nch = 16.

Download Full Size | PDF

Also, we plot in each subfigure (if applicable) a black contour line corresponding to the parameters (PSIG, PLO) that achieve an SNR of 9.2 dB (FEC threshold). If we now compare Figs. 10(g) and 10(h) in the last row, we notice that the operating range where error-free operation is achieved, i.e. when SNR > 9.2 dB which is the area above the contour, is tighter in case of Nch = 16 compared to Nch = 1 due to the residual SIG-SIG self beat noise impact in the case of colorless reception which degrades the SNR as PSIG increases and PLO decreases. Within the operating range, we notice that for a fixed PLO, there is a dynamic range for PSIG within which error-free operation is possible and that this dynamic range gets wider as PLO increases. Generally, the SNR degradation observed in the left edge of the surfaces in Figs. 10(a)10(h) is due to the fundamental receiver shot and thermal noises. On the other hand, we notice a SNR degradation in the bottom right corner (where the LSR is smallest) of all subfigures in the right column representing the colorless scenario, due to the residual SIG-SIG self-beat noise due to neighboring WDM channels in this case. In fact, one should also expect SNR degradation near the top right corner of all subfigures due to saturation effects when both PSIG and PLO are large. In this case, the TIAs may get overloaded leading to nonlinear distortion of the signal which in turn degrades the SNR. From a system design perspective, the TIA overload problem is governed by two important parameters namely the maximum single-ended DC current (IDC) and the peak-to-peak differential AC current (IACppd), which are calculated similar to [16] as follows

IDC=RLOPLO+NchRSIGPSIG,
IACppd=8RLOPLOPAPRRSIGPSIG,
where PAPR is the peak-to-average power ratio of the photocurrent of one WDM channel which depends on the residual CD and polarization orientation [16]. If those two currents calculated in Eqs. (14) and (15) are guaranteed to be below the maximum ratings of the TIA, TIA saturation will be avoided and linear operation will be achieved. In our experiment, the TIAs had maximum ratings of 4 mA and 3 mA for IDC and IACppd, respectively. Regarding DC saturation, if we use the effective responsivities of our CRx which are 0.025 and 0.029 A/W for the signal and LO ports, respectively, we can infer from Eq. (14) that for the experimental point with the largest PLO, PSIG and Nch in Fig. 10, i.e. when PLO = 15.5 dBm, PSIG = 5 dBm and Nch = 16, IDC will be 2.27 mA which is below 4 mA, meaning that no DC saturation has occurred in our experiment and that the results in Fig. 10 are not affected by the TIA DC overload. On the other hand, the AC overload problem affects a portion of our experimental data. If we use a PAPR of 10 dB [36] and the CRx responsivities defined above, we can use Eq. (15) to evaluate (PLO, PSIG) that corresponds to IACppd = 3 mA which yields the magenta lines plotted in Figs. 10(a)10(h) near the top right corner of each subfigure. This means that the experimental points above these contours were affected by TIA AC saturation whereas all points lying below these lines were not impacted by AC saturation and the TIAs were in the linear regime of operation. Finally, it is noteworthy that the two currents calculated according to Eqs. (14) and (15) are not the only key parameters that should be considered when assessing TIA saturation. In fact, one should also consider the common mode AC current (which is dependent on Nch). One can actually evaluate the single-ended AC current, i.e. before subtraction that takes place inside the TIA, which should include both the common mode and differential AC contributions. However, the TIA specification regarding the maximum common mode AC current was not available and hence, was not considered.

4.2.3 Verification of experimental results using the analytical SNR model

In this subsection, the goal is to verify the match between the experimental data presented in the preceding subsections and the behavior predicted by the analytical SNR model given by Eq. (11). In order to do so, the analytical SNR model in Eq. (11) is fitted to the experimental data. We first re-write Eq. (11) after re-expressing the denominator which represents the various noise sources as follows

SNR=PLOPSIGα1PLOPSIG+α2NloopsPLOPSIG+α3+α4(PLO+NchPSIG)+α5NchPSIG2,
where α1,α1,...,α5 are the coefficients that need to be estimated in order to fit the model in Eq. (16) to the experimental data. In fact, the first two terms in the denominator of Eq. (16) represent the LO beating with the noise coming along with the received signal which was represented by the first term (LO-ON beating) in the denominator of Eq. (11). Here, we split this term into two terms representing the LO beating with the noise accompanying the received signal in back-to-back and transmission scenarios, respectively. The reason to add a noise accompanying the signal in the back-to-back case is to take into account the noise due to the transmitter. This noise from the transmitter side impacts the signal even without transmission and originates from the booster EDFA, RF drivers of the IQ modulator, etc. Hence, the α1 term in Eq. (16) is independent of the transmission distance. On the other hand, the α2 term represents the LO beating with the incoming optical noise due to both ASE noise of the in-line EDFAs and NL-induced noise after dispersion uncompensated transmission. According to [27], as the number of WDM channels increases and the channel spacing gets smaller, i.e. more dense spectrum, the NL-induced noise accumulation along the link is approximately linear with the number of spans which corresponds to summing incoherently the powers of NL noise contributions from various spans. In our experiment, we always had 16 WDM channels with 50 GHz spacing propagating along the link which results in NL accumulation approximately proportional to Nspans1.06 [27], which is very close to being linear. Hence, both the ASE and NL-induced AWGNs can be safely assumed to scale linearly with the number of spans, and consequently, with the number of loops in the recirculating loop,Nloops. Next, by carefully comparing the first noise term in Eq. (11) with the first two noise terms in Eq. (16), we notice that the first noise term in Eq. (11) is expressed in terms of the product PLOPON whereas in Eq. (16), we make use of the proportionality between PON and PSIG for a fixed OSNR and write both the α1 and α2 terms in terms of PLOPSIG instead, i.e. the inverses of the OSNRs in back-to-back and after one loop transmission are imbedded into α1 and α2, respectively. We should note that the launch power, which affects PON by changing the power of the NL-induced noise contribution, was kept constant throughout the experiment. In fact, PSIG in Eq. (16) denotes the received signal power (not the launch power) that was swept using the VOA prior to the CRx after NL accumulation already happened along the fiber. In addition, it is noteworthy that the number of WDM channels propagating along the link in our experiment was also kept constant at 16, whereas the number of channels presented to the CRx, denoted by Nch in Eq. (16), was swept using the tunable bandwidth optical filter at the end of the link prior to the CRx. Hence, the NL-induced noise along the link in our experiment was constant. Next, the α3 term represents the thermal noise which is assumed independent of PLO and PSIG. Then, the α4 term is due to the shot noise including both the LO and SIG contributions. It is noteworthy that both α3 and α4 depend on the responsivity of the CRx beside various other parameters. Finally, the α5 term represents the residual SIG-SIG beating and hence is proportional to NchPSIG2 where α5 includes the impact of the effective CMRR of the signal port and the scaling factor β.

The next step is to estimate α1,α1,...,α5 in Eq. (16) using the experimental data plotted in Fig. 6. For the purpose of fitting the model, we use only the experimental data obtained after transmission corresponding to Nloops equal to 5, 10 and 15 because at these distances, the scaling factor β saturates to a constant value of 0.55 (see appendix B) and hence α5 can be assumed constant. Since we have many SNR data points at various operating parameters namely, PLO,PSIG,Nch and Nloops, estimating the coefficients α1,α1,...,α5 requires solving an over-determined system. In order to do so, we use the linear least-squares method where the following two vectors and matrix are formed

α=[α1α5],   y=[PLO1PSIG1/SNR1PLOmPSIGm/SNRm],   X=[PLO1PSIG1PLOmPSIGmNloops1PLO1PSIG1NloopsmPLOmPSIGm11PLO1+Nch1PSIG1PLOm+NchmPSIGmNch1PSIG12NchmPSIGm2],
where m is the number of data points. Using the least-squares method, an estimate α^ for the coefficients α is computed such that αα^2 is minimized and will be on the form [37]
α^=(XTX)1XTy,
where x denotes the norm of a vector x.

By plugging the estimated coefficients from Eq. (17) into the SNR model of Eq. (16), we now have a known analytical SNR model that can be employed to either verify the experimental data or evaluate the SNR at other operating conditions, e.g. another Nch. Figure 11 shows a comparison between the experimental SNR data and the evaluated SNR data from the fitted analytical model where the SNR is plotted against PSIG for various transmission distances. Figures 11(a) and 11(b) present the data when PLO is set to 12 dBm and Nch is 1 and 16, respectively whereas Figs. 11(c) and 11(d) present the data when PLO is set to 3 dBm and Nch is 1 and 16, respectively. In obtaining Fig. 11, we performed two least-squares fits of the SNR model: one for the back-to-back case and one for the transmission case, since α5 is expected to be different in both scenarios. The reason for the expected difference of α5 between the two scenarios is that α5 includes the effect of the scaling factor β which is dependent on the amount of CD (see appendix B for more details on the scaling factor β) and hence fitting the SNR model is done for the back-to-back data and transmitted data separately. Indeed, α5 found from fitting the model to the transmitted data was 1.74 times larger than the fitted α5 using the back-to-back data. This ratio matches well the expected value based on the theory in appendix B (see Fig. 13) if we use β of 0.55 for the transmitted data and β of 0.316 for the back-to-back data assuming the worst case polarization orientation of the received signal (45° relative to the PBS of the CRx). Furthermore, we should anticipate that α5 is the same for all transmitted data since β saturates quickly to a value of 0.55 when the residual CD is more than 3 ns/nm (see Fig. 13 in appendix B). Looking back at Fig. 11, we notice a good agreement between the SNR curves derived from the model compared to the experimental curves which verifies both the validity of the proposed model in Eq. (16) and the accuracy of the experimental data.

 figure: Fig. 11

Fig. 11 SNR versus PSIG comparing both experimental data points with the points evaluated from the fitted SNR analytical model after 0, 1600, 3200 and 4800 km at the following operating conditions: (a) PLO = 12 dBm and Nch = 1, (b) PLO = 12 dBm and Nch = 16, (c) PLO = 3 dBm and Nch = 1, and (d) PLO = 3 dBm and Nch = 16.

Download Full Size | PDF

Besides the agreement between the model and the experimental SNR data, we can also use the estimated value of α5 (in case of transmitted data) to calculate CMRR¯SIG assuming β=0.55 which can be easily done by comparing the last term in the denominators of Eqs. (11) and (16). We found that CMRR¯SIG obtained from the fitted coefficient is −19 dB. In order to verify this value, we used the analytical definition of CMRR¯SIG presented below Eq. (11) in section 3 in addition to the CMRR versus frequency curves obtained via characterization of the CRx which were depicted in Fig. 2(a) of section 2. The value found through the analytical expression of CMRR¯SIG was −17.5 dB which is in good agreement with the one obtained from fitting.

Next, we make use of the SNR model found in Eq. (16) to predict the SNR performance if the same CRx is used in operating conditions different from those of our experiment. For example, we use the model to predict the performance when the number of WDM channels co-incident on the CRx is larger than 16. In Fig. 12, we plot the SNR versus Nch after 0, 1600, 3200 and 4800 km transmission for both the model and experiment. It can be observed that the sweep range of Nch is limited to 16 in case of the experimental data whereas we employ the analytical model to predict the SNR if the CRx is used for colorless reception when Nch is up to 80. In Figs. 12(a) and 12(b), a high PLO of 12 dBm is chosen and PSIG is set to 0 and −12 dBm, respectively whereas in Figs. 12(c) and 12(d) a low PLO of 3 dBm is chosen and PSIG is set to 0 and −12 dBm, respectively. If we compare Figs. 12(a)12(d), we notice that the impact of increasing Nch on the SNR degradation is the most in Fig. 12(c) where the residual SIG-SIG beat noise has the largest impact for small PLO and high PSIG, i.e. small LSR. For example, the SNR after 1600 km transmission degrades by nearly 6 dB as Nch changes from 1 to 80. On the other hand, we observe no impact from varying Nch on the SNR performance in Fig. 12(b) when PLO = 12 dBm and PSIG = −12 dBm. In fact, the model allows us to predict error free operation (SNR > 9.2 dB) after 4800 km transmission when Nch = 80. In order to verify that our SNR predictions when Nch = 80 are correct from the TIA saturation perspective, we also use Eqs. (14) and (15) to verify if IDC and IACppd are going to be below the TIA ratings. Indeed, if we plug the largest PSIG and PLO used in Fig. 12 into Eqs. (14) and (15), i.e. 0 dBm and 12 dBm respectively, and use Nch = 80, we obtain IDC = 2.43 and IACppd = 2.69 mA which are below the TIA maximum ratings. We finally conclude that all SNR points drawn from the model in Fig. 12 are within the linear regime of operation of the TIAs.

 figure: Fig. 12

Fig. 12 SNR versus Nch comparing the experimental data points with the points evaluated from the fitted SNR analytical model after 0, 1600, 3200 and 4800 km at the following operating conditions: (a) PLO = 12 dBm and PSIG = 0 dBm, (b) PLO = 12 dBm and PSIG = −12 dBm, (c) PLO = 3 dBm and PSIG = 0 dBm, and (d) PLO = 3 dBm and PSIG = −12 dBm.

Download Full Size | PDF

5. Conclusion

In this paper, we present a comprehensive theoretical and experimental study of the performance of a SiP integrated balanced CRx when employed for colorless reception in a WDM setup. The layout of the CRx chip is depicted and characterization results of its parameters are shown. Then, an analytical SNR expression of the output photocurrent is derived rigorously while pinpointing noise sources in a colorless scheme. Next, experimental results are presented where 16 × 112 Gbps PDM-QPSK WDM channels are used. We study the impact of varying system parameters on the performance of colorless reception. Results show that the LO beating with incoming optical noise including ASE and NL-induced noises is a dominant noise source regardless of the received signal power. In addition, we identify the SIG-SIG beating due to OOB WDM channels as an additional important noise source in the regime of high received signal power. On the other hand, thermal and shot noises have a significant impact on SNR in the regime of low received signal power. In addition, we conclude that the CMRR and sensitivity are two important CRx specifications that impact the performance at high and low received signal power regimes, respectively. Finally, an excellent match between experimental and analytical SNRs is proven. After fitting the SNR model, we conclude that the SiP CRx is capable of colorless operation with a WDM spectrum comprising 80 channels as long as the LO-to-signal power ratio is properly set.

6. Appendix A: Derivation of ACFs, PSDs and powers of various terms in the differential photocurrent

In this appendix, we derive the time-averaged ACFs of all terms in the total differential photocurrent given by Eq. (7) from which the PSDs are evaluated through Fourier transformation (FT). Then, the PSDs derived are integrated over frequency to calculate the powers required to evaluate the SNR according to Eq. (10). We divide the appendix into four subsections which are laid out for different beating terms of the differential photocurrent given in Eq. (7).

6.1 Derivation of ACF, PSD and power of the LO-SIG and LO-ON beating terms in the differential photocurrent

We start by neglecting the effect of the impulse responses of hP(t) and hN(t). However, we will later modify the final result derived to include their impact. The time averaged ACF of the LO-SIG term in the differential photocurrent in Eq. (7) can be calculated starting from the definition of the ACF of ΔiLOSIG(t), obtained from the difference between the LO-SIG terms in the P and N photocurrents in Eqs. (5) and (6), yielding

Γ¯ΔiLOSIG(τ)=ΔiLOSIG(t)ΔiLOSIG(tτ)¯=RSIG,sPRLO,sPRe{E˜SIG,s(t)E˜LO*(t)ej2πΔνt}Re{E˜SIG,s(tτ)E˜LO*(tτ)ej2πΔν(tτ)}¯+RSIG,sNRLO,sNRe{E˜SIG,s(tδ)E˜LO*(tδ)ej2πΔν(tδ)}Re{E˜SIG,s(tτδ)E˜LO*(tτδ)ej2πΔν(tτδ)}¯+RSIG,sPRLO,sPRSIG,sNRLO,sNRe{E˜SIG,s(t)E˜LO*(t)ej2πΔνt}Re{E˜SIG,s(tτδ)E˜LO*(tτδ)ej2πΔν(tτδ)}¯+RSIG,sPRLO,sPRSIG,sNRLO,sNRe{E˜SIG,s(tδ)E˜LO*(tδ)ej2πΔν(tδ)}Re{E˜SIG,s(tτ)E˜LO*(tτ)ej2πΔν(tτ)}¯.
Using the relation in Eq. (4) to express the product of two real parts and also realizing that the terms having ej2πΔν(2tτ) cancel after time averaging as in [29], we arrive at
Γ¯ΔiLOSIG(τ)=12RSIG,sPRLO,sPRe{E˜SIG,s(t)E˜SIG,s*(tτ)E˜LO*(t)E˜LO(tτ)¯ej2πΔντ}+12RSIG,sNRLO,sNRe{E˜SIG,s(tδ)E˜SIG,s*(tτδ)E˜LO*(tδ)E˜LO(tτδ)¯ej2πΔντ}+12RSIG,sPRLO,sPRSIG,sNRLO,sNRe{E˜SIG,s(t)E˜SIG,s*(tτδ)E˜LO*(t)E˜LO(tτδ)¯ej2πΔν(τ+δ)}+12RSIG,sPRLO,sPRSIG,sNRLO,sNRe{E˜SIG,s(tδ)E˜SIG,s*(tτ)E˜LO*(tδ)E˜LO(tτ)¯ej2πΔν(τδ)}.
Then, if use the statistical independence between the LO and signal fields, we obtain
Γ¯ΔiLOSIG(τ)=12RSIG,sPRLO,sPRe{Γ¯E˜SIG(τ)ΓE˜LO*(τ)ej2πΔντ}+12RSIG,sNRLO,sNRe{Γ¯E˜SIG(τ)ΓE˜LO*(τ)ej2πΔντ}+12RSIG,sPRLO,sPRSIG,sNRLO,sNRe{Γ¯E˜SIG(τ+δ)ΓE˜LO*(τ+δ)ej2πΔν(τ+δ)}+12RSIG,sPRLO,sPRSIG,sNRLO,sNRe{Γ¯E˜SIG(τδ)ΓE˜LO*(τδ)ej2πΔν(τδ)}.
The first term in Eq. (19) follows directly from the first term in Eq. (18) after using the statistical independence between the LO and signal fields. However, the second, third and fourth terms are each obtained as follows
E˜SIG,s(tδ)E˜SIG,s*(tτδ)E˜LO*(tδ)E˜LO(tτδ)¯=E˜SIG,s(tδ)E˜SIG,s*(tτδ)¯E˜LO*(tδ)E˜LO(tτδ)=ΓE˜SIG(tδ,τ)¯ΓE˜LO*(τ)=Γ¯E˜SIG(τ)ΓE˜LO*(τ),
E˜SIG,s(t)E˜SIG,s*(tτδ)E˜LO*(t)E˜LO(tτδ)¯=E˜SIG,s(t)E˜SIG,s*(tτδ)¯E˜LO*(t)E˜LO(tτδ)=ΓE˜SIG(t,τ+δ)¯ΓE˜LO*(τ+δ)=Γ¯E˜SIG(τ+δ)ΓE˜LO*(τ+δ),
E˜SIG,s(tδ)E˜SIG,s*(tτ)E˜LO*(tδ)E˜LO(tτ)¯=E˜SIG,s(tδ)E˜SIG,s*(tτ)¯E˜LO*(tδ)E˜LO(tτ)=ΓE˜SIG(tδ,τδ)¯ΓE˜LO*(τδ)=Γ¯E˜SIG(τδ)ΓE˜LO*(τδ),
where the LO field is assumed to be stationary (unmodulated), whereas the signal field is cyclostationary with its time-dependent and lag-dependent (two-argument) ACF being periodic with a period equal to the symbol duration. After time averaging which involves omitting the time variable, it can be shown that ΓE˜SIG(tδ,τ)¯,ΓE˜SIG(t,τ+δ)¯, and ΓE˜SIG(tδ,τδ)¯ will finally lead to, Γ¯E˜SIG(τ+δ) and Γ¯E˜SIG(τδ), respectively.

In Eq. (19), if we neglect the lineshape of the LO, i.e. SE˜LO(f)=2PLOδ(fνLO), then ΓE˜LO(τ)=2PLO. The reason for multiplying PLO by 2 is that the total power of the complex envelope is twice that of the real field [29] since its respective ACF is given by ΓELO(τ)=12Re{ΓE˜LO(τ)ej2πντ}. Now if we substitute for ΓE˜LO(τ) in Eq. (19), we obtain

Γ¯ΔiLOSIG(τ)=PLO(RSIG,sPRLO,sP+RSIG,sNRLO,sN)Re{Γ¯E˜SIG(τ)ej2πΔντ}+PLORSIG,sPRLO,sPRSIG,sNRLO,sN.[Re{Γ¯E˜SIG(τ+δ)ej2πΔν(τ+δ)+Γ¯E˜SIG(τδ)ej2πΔν(τδ)}].
Then, if we take the FT of ACF, we get the PSD of the LO-SIG term as
SΔiLOSIG(f)=FT{Γ¯ΔiLOSIG(τ)}=PLO|RSIG,sPRLO,sP+ej2πfδRSIG,sNRLO,sN|212[SE˜SIG(fΔν)+SE˜SIG((f+Δν))].
In obtaining the expression in Eq. (21), we initially neglected the impulse responses of the P and N paths from the PD outputs up to current subtraction in the TIA which were denoted by hP(t) and hN(t). It can be shown that the effects of those two impulse responses can be taken into consideration in the ACF of the LO-SIG as follows
Γ¯ΔiLOSIG(τ)=PLORSIG,sPRLO,sP(Re{Γ¯E˜SIG(τ)ej2πΔντ}ΓhP(τ))+PLORSIG,sNRLO,sN(Re{Γ¯E˜SIG(τ)ej2πΔντ}ΓhN(τ))+PLORSIG,sPRLO,sPRSIG,sNRLO,sN(Re{Γ¯E˜SIG(τ+δ)ej2πΔν(τ+δ)}ΓhPhN(τ))+PLORSIG,sPRLO,sPRSIG,sNRLO,sN(Re{Γ¯E˜SIG(τδ)ej2πΔν(τδ)}ΓhNhP(τ)),
where ΓhP(τ), ΓhP(τ), ΓhPhN(τ) and ΓhNhP(τ) are the deterministic ACFs and cross-correlations of the impulse responses hP(t) and hN(t)that can be calculated as follows
ΓhP(τ)=hP(τ)hP*(τ),  ΓhN(τ)=hN(τ)hN*(τ),ΓhPhN(τ)=hP(τ)hN*(τ),  ΓhNhP(τ)=hN(τ)hP*(τ).
Taking the FT of the ACF in Eq. (22) leads to the PSD of the LO-SIG beating term as
SΔiLOSIG(f)=PLO|HP(f)RSIG,sPRLO,sP+ej2πfδHN(f)RSIG,sNRLO,sN|212[SE˜SIG(fΔν)+SE˜SIG((f+Δν))],
where HP(f) and HN(f) are the FT pairs of hP(t) and hN(t), respectively.

Next, we evaluate the PSD of the received signal field SE˜SIG(f), which appears in Eq. (23). In fact, this PSD is related to the pulse shape and modulation format used at the transmitter. In fact, we can write the received signal field as

E˜SIG(t)=n=(aXn+bYn)h(tnT),
where h(t)=g(t)hCD(t) is the overall impulse response including the transmitted pulse shape g(t) and CD impulse response hCD(t). Also, Xn and Yn are the nth information symbols on the two polarization states, i.e. Xn,Yn12{(1+j),(1+j),(1j),(1j)} in the case of QPSK. In addition, a and b are the matrix entries of the 2 x 2 Jones matrix (J) representing the polarization rotation due to the misalignment between the state of polarization of the received signal and the axis of the CRx PBS, i.e.
J=[abb*a*].
Using Eq. (24), we find the time averaged ACF of E˜SIG(t), denoted by Γ¯E˜SIG(τ), as follows
Γ¯E˜SIG(τ)=E˜SIG(t)E˜SIG*(tτ)¯=|a|2n=m=XnXm*h(tnT)h*(tτmT)¯+|b|2n=m=YnYm*h(tnT)h*(tτmT)¯+ab*n=m=XnYm*h(tnT)h*(tτmT)¯+a*bn=m=Xn*Ymh(tnT)h*(tτmT)¯.
We know that for QPSK modulated symbols on the same polarization, we have
XnXm*=YnYm*={1           n=m0          nm,
and for independent QPSK symbols on different polarizations, we have
XnYm*=Xn*Ym=0.
Hence, we can substitute with the statistical averages evaluated above into Eq. (26) to derive the ACF of the signal field as follows
Γ¯E˜SIG(τ)=(|a|2+|b|2)n=h(tnT)h*(tτnT)¯=1Tn=T/2T/2h(tnT)h*(tτnT)dt=1Th(t)h*(tτ)dt=1TΓh(τ),
where Γh(τ)=h(t)h*(tτ)dt is the deterministic ACF of the overall pulse shape and where we used the fact that |a|2+|b|2=1 since J is unitary. The integral in the second line of Eq. (27) comes from the time averaging operator f(t)¯ being 1/TT/2T/2f(t)dt provided that f(t) is periodic with a period T. Hence, the PSD of the signal field SE˜SIG(f) required in Eq. (23) is given by the FT of the ACF as follows
SE˜SIG(f)=1T|H(f)|2=1T|G(f)|2,
where H(f) and G(f) are the FT pairs of h(t) and g(t), respectively. In the second equality of Eq. (28), we used the fact that |HCD(f)|2=1. Now if we back substitute in Eq. (23) after neglecting Δν (since it will be corrected in DSP afterwards) and integrate over frequency, we get the power of the desired LO-SIG photocurrent term (numerator of the SNR in Eq. (10))
PΔiLOSIG=PLO|HP(f)RSIG,sPRLO,sP+HN(f)ej2πfδRSIG,sNRLO,sN|21T|G(f)|2|HrxDSP(f)|2df=2PLOPSIG|HP(f)RSIG,sPRLO,sP+HN(f)ej2πfδRSIG,sNRLO,sN|2|G(f)|2|HrxDSP(f)|2df|G(f)|2df,
where we employ the fact that the integral of the PSD of the complex envelope of the received signal field in Eq. (28) gives 2PSIG.

Regarding the LO-ON beating term, its ACF and PSD can be derived in a similar way as that used to obtain those quantities for the LO-ON term. In fact, the PSD of the LO-ASE beating will have a similar form to the one derived for the LO-SIG term which was given by Eq. (23), and can be expressed as

SΔiLOON(f)=PLO|HP(f)RSIG,sPRLO,sP+HN(f)ej2πfδRSIG,sNRLO,sN|212[SE˜ON(fΔν)+SE˜ON((f+Δν))],
where SE˜ON(f) is the PSD of the complex envelope of the ASE field within one channel, which can be assumed as that of a narrowband filtered white noise given by
SE˜ON(f)={2No         |f|Bch/20              |f|>Bch/2,
where Nois the PSD level of the real-valued ASE noise field in W/Hz and Bch is the optical bandwidth of one WDM channel. In Eq. (31), we made use of the assumption that the NL-induced noise portion inside E˜ON can be assumed locally white, i.e. it has an approximately flat spectrum within one WDM channel [27]. Hence, the total optical noise field E˜ON can be assumed white. Next, we can easily find the variance (or AC electrical power) of the LO-ON beating photocurrent term in the denominator of the SNR in Eq. (10) as follows

σΔiLOON2=2PLOPONBch/2Bch/2|HP(f)RSIG,sPRLO,sP+HN(f)ej2πfδRSIG,sNRLO,sN|2|HrxDSP(f)|2dfBch.

6.2 Derivation of ACF, PSD and power of the residual SIG-SIG beating term in the differential photocurrent

Similar to subsection 6.1 in this appendix, we begin by neglecting the effect of the impulse responses hP(t), hN(t) in Eq. (7) and later take them into consideration by modifying the final PSD expression obtained. We start by writing the definition of the time averaged ACF of ΔiSIGSIG(t) obtained from the difference between the SIG-SIG beating terms in the P and N photocurrents in Eqs. (5) and (6), yielding

Γ¯ΔiSIGSIG(τ)=ΔiSIGSIG(t)ΔiSIGSIG(tτ)¯=14[i=1Nchj=1NchRSIG,iPRSIG,jP|E˜SIG,i(t)|2|E˜SIG,j(tτ)|2¯+i=1Nchj=1NchRSIG,iNRSIG,jN|E˜SIG,i(tδ)|2|E˜SIG,j(tτδ)|2¯i=1Nchj=1NchRSIG,iPRSIG,jN|E˜SIG,i(t)|2|E˜SIG,j(tτδ)|2¯i=1Nchj=1NchRSIG,iNRSIG,jP|E˜SIG,i(tδ)|2|E˜SIG,j(tτ)|2¯].
In the above equation, all summation terms where ij will involve the intensities of different channels which can be assumed independent, i.e. the expectation of the product can be evaluated as the product of expectations. Hence,
Γ¯ΔiSIGSIG(τ)=i=1Nchj=1jiNch(RSIG,iPRSIG,jP+RSIG,iNRSIG,jN2RSIG,iPRSIG,jN)PSIG,iPSIG,j+14i=1NchRSIG,iP2|E˜SIG,i(t)|2|E˜SIG,i(tτ)|2¯+14i=1NchRSIG,iN2|E˜SIG,i(tδ)|2|E˜SIG,i(tτδ)|2¯14i=1NchRSIG,iPRSIG,iN|E˜SIG,i(t)|2|E˜SIG,i(tτδ)|2¯14i=1NchRSIG,iPRSIG,iN|E˜SIG,i(tδ)|2|E˜SIG,i(tτ)|2¯=i=1Nchj=1jiNch(RSIG,iPRSIG,jP+RSIG,iNRSIG,jN2RSIG,iPRSIG,jN)PSIG,iPSIG,j+14i=1NchRSIG,iP2Γ¯|E˜SIG,i|2(τ)+14i=1NchRSIG,iN2Γ¯|E˜SIG,i|2(τ)14i=1NchRSIG,iPRSIG,iNΓ¯|E˜SIG,i|2(τ+δ)14i=1NchRSIG,iPRSIG,iNΓ¯|E˜SIG,i|2(τδ),
where Γ¯|E˜SIG,i|2(τ) is the time averaged ACF of the squared envelope of the electric field of the ith WDM signal which will be derived later. Then, the PSD of the SIG-SIG beating term in the differential photocurrent can be evaluated by taking the FT of its ACF given in Eq. (33) and written as

SΔiSIGSIG(f)=δ(f)i=1Nchj=1jiNch(RSIG,iPRSIG,jP+RSIG,iNRSIG,jN2RSIG,iPRSIG,jN)PSIG,iPSIG,j+14i=1Nch|RSIG,iPej2πfδRSIG,iN|2S|E˜SIG,i|2(f).

Now, we can take into consideration the impact of the impulse responses hP(t) and hN(t) in a way similar to that described prior to Eq. (23), which yields

SΔiSIGSIG(f)=δ(f)i=1Nchj=1jiNch(RSIG,iPRSIG,jP+RSIG,iNRSIG,jN2RSIG,iPRSIG,jN)PSIG,iPSIG,j+14i=1Nch|HP(f)RSIG,iPHN(f)ej2πfδRSIG,iN|2S|E˜SIG,i|2(f),
where the first term in Eq. (35) represents a DC portion that is not of interest while calculating the residual SIG-SIG beat noise power σΔiSIGSIG2in Eq. (10). The second term in Eq. (35) represents the sum of the PSDs of the field intensities of all channels scaled by the quantity |HP(f)RSIG,iPHN(f)ej2πfδRSIG,iN|2 which represents the ability of the CRx to reject the common mode on both the P and N PDs. It is evident that this quantity is frequency dependent due to both the time skew and unequal frequency responses of the P and N paths.

Our next step is to derivethe PSD of the field intensity of one channel S|E˜SIG,i|2(f) required in Eq. (35). Without loss of generality, we will drop the subscript i and focus on the PSD of the squared envelope of one of the WDM channels S|E˜SIG|2(f). We begin by evaluating the squared envelope of the received signal field given by Eq. (24) and then finding the time averaged ACF of the squared field envelope Γ¯|E˜SIG|2(τ) as

Γ¯|E˜SIG|2(τ)=E˜SIG(t)E˜SIG*(t)E˜SIG(tτ)E˜SIG*(tτ)¯=n=m=l=k=(aXn+bYn)(aXm+bYm)*(aXl+bYl)(aXk+bYk)*h(tnT)h*(tmT)h(tτlT)h*(tτkT)¯,
where the statistical average of the product of four brackets can be evaluated as follows
(aXn+bYn)(aXm+bYm)*(aXl+bYl)(aXk+bYk)*=|a|4XnXm*XlXk*+|b|4YnYm*YlYk*+|a|2|b|2XnXm*YlYk*+|a|2|b|2YnYm*XlXk*+|a|2|b|2YnXm*XlYk*+|a|2|b|2XnYm*YlXk*+terms with vanishing statistical averages,
while noting that
XnXm*XlXk*=YnYm*YlYk*={1           n=m&l=k1          n=k&m=l&nm0          otherwise,
XnXm*YlYk*=YnYm*XlXk*={1           n=m&l=k0          otherwise,
YnXm*XlYk*=XnYm*YlXk*={1           n=k&l=m0          otherwise.
Back substituting into Eq. (36) yields
Γ¯|E˜SIG|2(τ)=(|a|4+|b|4)[n=l=|h(tnT)|2|h(tτlT)|2¯+n=m=nmh(tnT)h*(tmT)h(tτmT)h*(tτnT)¯]+2|a|2|b|2[n=l=|h(tnT)|2|h(tτlT)|2¯+n=m=h(tnT)h*(tmT)h(tτmT)h*(tτnT)¯]=1T(|a|4+|b|4+2|a|2|b|2)l=Γ|h|2(τlT)+1T(|a|4+|b|4)k=k0Ch(τ,k)+1T2|a|2|b|2k=Ch(τ,k),
where Γ|h|2(τ)=|h(t)|2|h(tτ)|2dt is the deterministic ACF of the squared envelope of the overall pulse shape and Ch(τ,k)=h(t)h*(tτ)h*(tkT)h(tτkT)dt. In order to obtain the last line in Eq. (37), we perform some lengthy algebraic manipulations on the two double summations. The first double summation is manipulated as
n=l=|h(tnT)|2|h(tτlT)|2¯=1Tn=l=T/2T/2|h(tnT)|2|h(tτlT)|2dt=1Tn=T/2T/2|h(tnT)|2l=|h(tτlT)|2dt     ,lett-nT=z=1Tn=nTT/2nT+T/2|h(z)|2l=|h(zτ(ln)T)|2dz=1Tl=|h(z)|2|h(zτlT)|2dz=1Tl=Γ|h|2(τlT),          where     Γ|h|2(τ)=|h(t)|2|h(tτ)|2dt
and the second double summation is evaluated similarly as
n=m=nmh(tnT)h*(tmT)h(tτmT)h*(tτnT)¯=1Tn=m=nmT/2T/2h(tnT)h*(tmT)h(tτmT)h*(tτnT)dt=1Tn=T/2T/2h(tnT)h*(tτnT)[m=nmh*(tmT)h(tτmT)]dt   ,let  t-nT=z=1Tn=nTT/2nT+T/2h(z)h*(zτ)[m=nmh*(z(mn)T)h(zτ(mn)T)]dz=1Tk=k0h(t)h*(tτ)h*(tkT)h(tτkT)dt=1Tk=k0Ch(τ,k),         where  Ch(τ,k)= h(t)h*(tτ)h*(tkT)h(tτkT)dt   
Next, applying the FT to the ACF obtained in Eq. (37) yields the desired PSD of the squared envelope (or intensity) of one channel S|E˜SIG|2(f) on the form
S|E˜SIG|2(f)=1T2|FT{h(t)h*(t)}|2l=δ(flT)+1Tk=|FT{h(t)h*(tkT)}|21T(|a|4+|b|4)|FT{h(t)h*(t)}|2=1T2|H(f)H*(f)|2l=δ(flT)+1Tk=|H(f)(ej2πfkTH*(f))|21T(|a|4+|b|4)|H(f)H*(f)|2,
where in obtaining Eq. (38), we first proved and then employed the following FT pairs
FT{Γ|h|2(τ)}=|FT{h(t)h*(t)}|2=|H(f)H*(f)|2,
FT{Ch(τ,k)}= |FT{h(t)h*(tkT)}|2=|H(f)(ej2πfkTH*(f))|2.
By carefully examining Eq. (38), we notice that the first term of S|E˜SIG|2(f) is a summation of Dirac delta functions at discrete frequencies of l/T where l in an integer. The second and third terms represent the continuous portion of the PSD which depends on the pulse shape, residual CD and polarization orientation. Next, we make use of the derived PSD of the intensity of one channel in Eq. (38) to evaluate the PSD of the residual SIG-SIG beating contribution in the differential photocurrent given by Eq. (35) while keeping only the relevant AC terms required for the SNR evaluation
SΔiSIGSIGAC(f)=14i=1Nch|HP(f)RSIG,iPHN(f)ej2πfδRSIG,iN|2S|E˜SIG,i|2AC(f),
where the PSD of the intensity of the ith channel is written based on the form given by Eq. (38) in terms of the pulse shape and polarization orientation of the ith channel as follows
S|E˜SIG,i|2AC(f)=1Tk=|FT{hi(t)hi*(tkT)}|21T(|ai|4+|bi|4)|FT{hi(t)hi*(t)}|2,
where the subscripts i in Eq. (40) take into account the possibility of having different pulse shapes, CD and polarization orientation of different WDM channels. Finally if we integrate the PSD in Eq. (39) over frequency, we obtain the variance (AC power) of the residual SIG-SIG beating contribution in the differential photocurrent (required in Eq. (10) to evaluate the SNR) as follows
σΔiSIGSIG2=14i=1Nch|HP(f)RSIG,iPHN(f)ej2πfδRSIG,iN|2S|E˜SIG,i|2AC(f)|HrxDSP(f)|2df=14Nchσ|E˜SIG,se|22|HP(f)RSIGPHN(f)ej2πfδRSIGN|2S|E˜SIG|2AC(f)|HrxDSP(f)|2dfS|E˜SIG|2AC(f)|Hav(f)|2|HrxDSP(f)|2df,
where in the second equality, we assumed the same responsivity for different WDM channels as well as equal intensity PSDs for various channels, i.e. S|E˜SIG,i|2AC(f) is the same for any i which implies assuming the same CD, pulse shape and polarization orientation for different WDM channels [38]. Also, σ|E˜SIG,se|22=S|E˜SIG|2AC(f)|Hav(f)|2|HrxDSP(f)|2df is the variance (or AC electrical power) of the photocurrent generated by the field intensity of one channel if detected with a single-ended PD having a unity responsivity and frequency response equal to the average frequency responses of the P and N PDs followed by DSP filtering. It is also noticeable that 14σ|E˜SIG,se|22=σ|ESIG,se|22 because the power of the complex field is twice that of the real field. In addition, it can be shown that σ|ESIG,se|22=βPSIG2, where β is a fractional scaling factor defined as the ratio between the AC power of the filtered photocurrent generated by the field intensity of one channel and its DC electrical power (which is the squared DC power of the field itself assuming a unit responsivity). In fact, β was shown by simulations in [14] to be dependent on the residual CD and polarization orientation (see appendix B for further details on the scaling factor β).

6.3 Power of the thermal and shot noise contributions in the photocurrent difference

Since shot noises generated from P and N PDs are uncorrelated, we can write the variance of their photocurrent difference Δish(t)=ishP(t)ishN(t) assuming, for simplicity, flat responsivity across the WDM spectrum and equal received signal powers for all WDM channels as follows

σΔish2=2e[(RLON+RLOP)PLO+Nch(RSIGN+RSIGP)PSIG]Δf4e[RLOavPLO+NchRSIGavPSIG]Δf,
where Δf is the effective BW that depends on the overall frequency response of the receiver including PDs, TIAs, and DSP filtering as well as the transmitted pulse shape. In the second equality, we use the averaged LO and SIG responsivities over both the P and N PDs. Next, we also calculate the variance of the final photocurrent due to thermal noise as follows
σΔith2=(iTIA)2Δf,
where iTIA is the input-referred noise current density which is an important specification for the TIAs, often provided in pA/Hz.

6.4 Expressions of PSDs of the LO-LO, SIG-ON and ON-ON beating contributions in the final differential photocurrent

Although we neglected the LO-LO, SIG-ON and ON-ON beating contributions when evaluating the PSD of the differential photocurrent in Eq. (9) since they were found negligible in our experiment, we present hereafter the final expressions of the PSDs of these three beat noise sources for the sake of completeness.

The PSDs of the differential photocurrents ΔiLOLO(t), ΔiSIGON(t), and ΔiONON(t) can be obtained using the same procedure employed for other beating terms, that is by first evaluating the ACF of the photocurrent which, after a FT leads to the PSD. Below are the final expressions of the PSDs of the three beating contributions assuming the same responsivity for all WDM channels

SΔiLOLO(f)=|HP(f)RLOPHN(f)ej2πfδRLON|2PLO2(δ(f)+RIN(f)),
SΔiSIGON(f)=14|HP(f)RSIGPHN(f)ej2πfδRSIGN|2Nch(SE˜SIG(f)SE˜ON(f)+SE˜SIG(f)SE˜ON(f)),
SΔiONON(f)=|HP(f)RSIGPHN(f)ej2πfδRSIGN|2((NoBo)2δ(f)+No2Botri(fBo)),
where RIN(f) is the spectrum of the relative intensity noise of the LO laser defined as
RIN(f)=FT{δPLO(t)δPLO(tτ)/PLO2},
where δPLO(t) denotes the random intensity fluctuations of the LO laser optical power around its average power PLO. Also, the triangular function in Eq. (46) is defined as
tri(fBo)={1|f|Bo            |f|Bo0                    |f|>Bo,
where Bo is the optical bandwidth of the WDM spectrum.

7. Appendix B: Dependence of the scaling factor of the SIG-SIG beating photocurrent on the residual CD and polarization orientation

In this appendix, we study in detail the scaling factor β presented in section 3. It is a fractional scaling factor equal to the ratio of the variance (or AC power) of the filtered photocurrent generated by the field intensity of one WDM channel relative to the DC electrical power of the photocurrent (which is the squared average optical power of the field itself assuming unit responsivity). In fact, it can be expressed as follows

β=14S|E˜SIG|2AC(f)|Hav(f)|2|HrxDSP(f)|2dfPSIG2,
where S|E˜SIG|2AC(f) is the AC portion of the PSD of the intensity of the complex envelope of the signal field defined in Eq. (38), Hav(f) is the average frequency response of the P and N ports of the CRx, HrxDSP(f) is the frequency response representing all filtering performed at the receiver side DSP, and PSIG is the average optical power of the field of one channel. By looking back at Eq. (38), we notice that S|E˜SIG|2AC(f) depends on the pulse shape and the residual CD (included in h(t)) as well as on the polarization orientation of the received signal relative to the PBS axis of the CRx (represented by a and b). Hence, the scaling factor β will in turn be dependent on the above parameters. In addition, β also tracks the overlap between the intensity spectrum S|E˜SIG|2AC(f) and the overall transfer function of the receiver side (Hav(f)and HrxDSP(f)). In fact, the intensity spectrum is broader than the LO-SIG beating spectrum (desired term) which can be explained in view of Eq. (38) as follows: S|E˜SIG|2AC(f) is proportional to k=|H(f)(ej2πfkTH*(f))|2 meaning that the intensity spectrum is proportional to the squared magnitude of the convolution of the frequency response of the received field with phase shifted versions of itself and hence, the bandwidth of the intensity spectrum is expected to be double that of the received field. Thus, almost one half of the power of the intensity spectrum is rejected by filtering from the CRx itself and from the receiver DSP since this filtering is designated to pass only the desired signal spectrum.

In order to evaluate the impact of residual CD and polarization orientation on the scaling factor β, simulations are carried using a 28 Gbaud PDM-QPSK signal having 105 symbols. First, the symbols were pulse shaped using NRZ pulses at 16 samples/symbol. The NRZ pulse shape used in the simulation produces the QPSK intensity eyediagram shown in Fig. 13(a) which matches quite well the experimental intensity eyediagram obtained by direct detection after the single polarization IQ modulator (see Fig. 13(b)). After pulse shaping, the waveform was transmitted through various CD amounts ranging from 0 to 5 ns/nm. After CD was applied, the polarization of the waveform was rotated by some Jones angle (θ) by multiplying by the Jones matrix J defined in Eq. (25) where a=cos(θ) and b=sin(θ), in order to simulate the misalignment between the state of polarization of the received signal and the axes of the CRx PBS. Then, the resulting waveform was convolved with the overall receiver’s impulse response (assumed to be the matched filter in simulations). Finally, the variance (AC power) of the squared envelope (intensity) of the resulting filtered waveform is calculated and divided by its DC power resulting in the scaling factor β.

 figure: Fig. 13

Fig. 13 (a) Simulated intensity eyediagram of the single polarization QPSK signal using NRZ pulse shape, (b) experimental intensity eyediagram obtained after the single polarization IQ modulator, (c) scaling factor (β) versus the amount of residual CD for various polarization rotation angles (θ) obtained from simulations and theory on a 28 Gbaud PDM-QPSK signal.

Download Full Size | PDF

Figure 13(c) plots β as evaluated from the simulation described above versus the amount of CD in ns/nm for various polarization rotation angles θ. For the purpose of comparison, we also plot in Fig. 13(c) β versus the amount of CD at different polarization orientations as calculated from our theoretical expression in Eq. (48) with the application of Eq. (38) (derived in subsection 6.2) that expresses the PSD of the intensity in terms of the pulse shape, CD and polarization orientation. An excellent agreement is observed between the β values calculated from simulation and theory which verifies our analytical expressions. It is also noteworthy that the theoretical expressions derived can be employed to calculate β for any other pulse shape by plugging the impulse response of the pulse shape in Eq. (38) and then using Eq. (38) into Eq. (48). In addition, β can be evaluated for other modulation formats in a similar fashion to the one used for PDM-QPSK, however, a part of the derivation in subsection 6.2 has to be re-worked for the desired modulation format, specifically, the statistical average performed on the information symbols in Eq. (36), which was evaluated assuming PDM-QPSK modulation, has to be re-calculated for the new format.

Looking back to Fig. 13(c), we notice that the impact of the polarization orientation on β is significant for CD values less than 2 ns/nm (which corresponds to ~117 km propagation over a fiber with dispersion parameter of 17 ps/(nm.km)). Within this CD range, β is at its maximum when θ=π/4 since for this polarization rotation angle, maximum crosstalk (and beating) between the two orthogonal polarizations of the received signal after passing through the PBS of the CRx results. Hence, this maximizes the peak-to-average ratio of the electrical signals at the output of the CRx which in turn increase the AC power of the photocurrent and finally increases β. On the other hand when the residual CD is larger than 2 ns/nm, the impact of θ is less significant. In fact as the residual CD increases beyond 3 ns/nm, β increases until it saturates at a maximum value of around 0.55 regardless of polarization.

References and links

1. S. J. Savory, “Digital coherent optical receivers: algorithms and subsystems,” IEEE J. Sel. Top. Quantum Electron. 16, 1164–1179 (2010).

2. J. Kahn, “Modulation and detection techniques for optical communication systems,” in Proceedings of Optical Amplifiers and Their Applications/Coherent Optical Technologies and Applications, Technical Digest (CD) (Optical Society of America, 2006), paper CThC1.

3. K. Roberts, D. Beckett, D. Boertjes, J. Berthold, and C. Laperle, “100G and beyond with digital coherent signal processing,” IEEE Commun. Mag. 48(7), 62–69 (2010). [CrossRef]  

4. C. Laperle, “Advances in high-speed ADCs, DACs, and DSP for optical transceivers,” in Proceedings of OFC (17–21 March 2013).

5. S. J. Savory, “Digital filters for coherent optical receivers,” Opt. Express 16(2), 804–817 (2008). [CrossRef]   [PubMed]  

6. E. Ip and J. M. Kahn, “Fiber impairment compensation using coherent detection and digital signal processing,” J. Lightwave Technol. 28, 502–519 (2010).

7. K. Roberts, “100G and beyond,” in Proceedings of 17th Opto-Electronics and Communications Conference (OECC) (2012), pp. 144–145. [CrossRef]  

8. R. Essiambre, R. Ryf, N. Fontaine, and S. Randel, “Breakthroughs in photonics 2012: Space-division multiplexing in multimode and multicore fibers for high-capacity optical communication,” IEEE Photonics J. 5(2), 0701307 (2013). [CrossRef]  

9. K. Roberts and C. Laperle, “Flexible transceivers,” in Proceedings of European Conference and Exhibition on Optical Communication, OSA Technical Digest (Optical Society of America, 2012), paper We.3.A.3. [CrossRef]  

10. D. Plant, Q. Zhuge, M. Morsy-Osman, M. Chagnon, X. Xu, and M. Qiu, “Flexible transceivers using adaptive digital signal processing for single carrier and OFDM,” in Proceedings of Optical Fiber Communication Conference (OFC) (17–21 March 2013).

11. Cisco white paper, “Cisco Visual Networking Index: Forecast and methodology, 2012–2017” (Cisco, 2013), http://www.cisco.com/en/US/solutions/collateral/ns341/ns525/ns537/ns705/ns827/white_paper_c11-481360.pdf.

12. Emcore product web page, “ITLA (Integrable Tunable Laser Assembly),” http://www.emcore.com/fiber-optics/telecom-tunable-dwdm/itla-integrable-tunable-laser-assembly/.

13. L. Kazovsky, “Multichannel coherent optical communications systems,” J. Lightwave Technol. 5, 1095–1102 (1987).

14. L. Nelson, S. Woodward, S. Foo, M. Moyer, D. Beckett, M. OSullivan, and P. Magill, “Detection of a single 40 Gb/s polarization-multiplexed QPSK channel with a real-time intradyne receiver in the presence of multiple coincident WDM channels,” J. Lightwave Technol. 28, 2933–2943 (2010).

15. C. Xie, P. J. Winzer, G. Raybon, A. H. Gnauck, B. Zhu, T. Geisler, and B. Edvold, “Colorless coherent receiver using 3x3 coupler hybrids and single-ended detection,” Opt. Express 20(2), 1164–1171 (2012). [CrossRef]   [PubMed]  

16. B. Zhang, C. Malouin, and T. J. Schmidt, “Towards full band colorless reception with coherent balanced receivers,” Opt. Express 20(9), 10339–10352 (2012). [CrossRef]   [PubMed]  

17. L. Nelson, X. Zhou, R. Isaac, Y. Lin, J.Chon, and W. Way, “Colorless reception of a single 100Gb/s channel from 80 coincident channels via an intradyne coherent receiver,” in Proceedings of IEEE Photonics Conference (2012), paper TuE4.

18. C. Doerr, P. Winzer, Y. Chen, S. Chandrasekhar, M. Rasras, L. Chen, T. Liow, K. Ang, and G. Lo, “Monolithic polarization and phase diversity coherent receiver in silicon,” J. Lightwave Technol. 28, 520–525 (2010).

19. Y. Painchaud, M. Pelletier, M. Poulin, F. Pelletier, C. Latrasse, G. Robidoux, S. Savard, J. Gagné, V. Trudel, M. Picard, P. Poulin, P. Sirois, F. D'Amours, D. Asselin, S. Paquet, C. Paquet, M. Cyr, M. Guy, M. Morsy-Osman, Q. Zhuge, X. Xu, M. Chagnon, and D. Plant, “Ultra-compact coherent receiver based on hybrid integration on silicon,” in Proceedings OFC (2013), paper OM2J.

20. A. Matiss, R. Ludwig, J. Fischer, L. Molle, C. Schubert, C. Leonhardt, H. Bach, R. Kunkel, and A. Umbach, “Novel integrated coherent receiver module for 100G serial transmission,” in Proc. OFC (2010), paper PDPB3. [CrossRef]  

21. C. Doerr, L. Zhang, P. Winzer, N. Weimann, V. Houtsma, T. Hu, N. Sauer, L. Buhl, D. Neilson, S. Chandrasekhar, and Y. Chen, “Monolithic InP dual-polarization and dual-quadrature coherent receiver,” IEEE Photonics Technol. Lett. 23, 694–696 (2011).

22. D. Liang and J. E. Bowers, “Photonic integration: Si or InP substrates?” Electron. Lett. 45(12), 578–581 (2009). [CrossRef]  

23. C. R. Doerr, L. L. Buhl, Y. Baeyens, R. Aroca, S. Chandrasekhar, X. Liu, L. Chen, and Y. K. Chen, “Packaged monolithic silicon 112-Gb/s coherent receiver,” IEEE Photonics Technol. Lett. 23, 762–764 (2011).

24. M. Morsy-Osman, M. Chagnon, X. Xu, Q. Zhuge, M. Poulin, Y. Painchaud, M. Pelletier, C. Paquet, and D. V. Plant, “Colorless and preamplifierless reception using an integrated Si-photonic coherent receiver,” IEEE Photonics Technol. Lett. 25, 1027–1030 (2013).

25. A. Carena, V. Curri, G. Bosco, P. Poggiolini, and F. Forghieri, “Modeling of the impact of nonlinear propagation effects in uncompensated optical coherent transmission links,” J. Lightwave Technol. 30, 1524–1539 (2012).

26. P. Poggiolini, A. Carena, V. Curri, G. Bosco, and F. Forghieri, “Analytical modeling of nonlinear propagation in uncompensated optical transmission links,” IEEE Photonics Technol. Lett. 23, 742–744 (2011).

27. P. Poggiolini, “The GN model of non-linear propagation in uncompensated coherent optical systems,” J. Lightwave Technol. 30, 3857–3879 (2012).

28. Y. Painchaud, M. Poulin, M. Morin, and M. Têtu, “Performance of balanced detection in a coherent receiver,” Opt. Express 17(5), 3659–3672 (2009). [CrossRef]   [PubMed]  

29. M. Nazarathy, W. Sorin, D. Baney, and S. Newton, “Spectral analysis of optical mixing measurements,” J. Lightwave Technol. 7, 1083–1096 (1989).

30. Although both LO and signal fields are stationary in their CW form, a modulated signal field is non-stationary rendering the final photocurrent difference Δi non-stationary. Hence, we apply both ensemble and time averaging to obtain a single argument time-averaged ACF from which the PSD is obtained using the Wiener–Khintchine theorem.

31. Despite that we neglect the ON-ON, SIG-ON and LO-LO terms in the final SNR expression derived, we provide the final expressions of the PSDs of each of those three noise contributions in appendix A (see subsection 6.4) for completeness.

32. P. Ciblat and M. Ghogho, “Blind NLLS carrier frequency-offset estimation for QAM, PSK, and PAM modulations: Performance at low SNR,” IEEE Trans. Commun. 54, 1725–1730 (2006).

33. I. Fatadin, D. Ives, and S. J. Savory, “Blind equalization and carrier phase recovery in a 16-QAM optical coherent system,” J. Lightwave Technol. 27, 3042–3049 (2009).

34. ITU-T Recommendation G.975.1 Appendix I.9 (2004).

35. J. Proakis, Digital Communications (McGraw Hill, 2001), Chap. 5.

36. Although the PAPR depends on the accumulated CD, polarization orientation and number of channels, we use a PAPR value of 10 dB which represents an average value evaluated from the different raw waveforms captured by the real-time scope at various distances, number of channels, etc.

37. S. Haykin, Adaptive Filter Theory (Prentice Hall, 1996), Chap. 11.

38. Although assuming the same pulse shape, CD and polarization orientation is not true in practice, the discussion in appendix B regarding the scaling factor β will show that after a short transmission distance, the PSDs of the intensities of various WDM channels (and hence the power of the SIG-SIG beating) should eventually be the same regardless of the pulse shape, polarization orientation and CD.

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 (13)

Fig. 1
Fig. 1 (a) Picture of the PIC of the CRx showing the main buliding blocks, (b) Overall frequency response of the Si-photonic CRx assembly.
Fig. 2
Fig. 2 CMRR versus frequency curves measured at 192.9 THz at the four output RF ports of the CRx when light is applied at the (a) signal port, (b) LO port.
Fig. 3
Fig. 3 Experimental setup for WDM experiments, (IL: Interleaver, DFB laser: Distributed Feedback laser, AWG: Arrayed Waveguide Grating, PBS: Polarization Beam Splitter, PBC: Polarization Beam Combiner, VODL: Variable Optical Delay Line, ECL: External Cavity Laser, OSA: Optical Spectrum Analyzer, VOA: Variable Optical Attenuator).
Fig. 4
Fig. 4 BER for all WDM channels using fully colorless and preamplifierless reception at different transmission distances and at per channel received signal power of −3 and −21 dBm. Insets show the received WDM spectra in back-to-back and after 3200 km transmission where we observe a fairly flat spectrum due to the gain flattening waveshaper in the loop.
Fig. 5
Fig. 5 Measurement results of channel 9: (a) BER versus transmission distance, (b) SNR evaluated after offline DSP versus distance, (c) BER versus SNR using both experimental results and theory where the experimental SNR corresponding to FEC threshold is marked.
Fig. 6
Fig. 6 SNR of channel 9 versus received signal power per channel PSIG after 0, 1600, 3200 and 4800 km transmission with varying number of channels presented to the CRx at an LO power of: (a) 15.5 dBm, (b) 12 dBm, (c) 9 dBm, (d) 6 dBm, (e) 3 dBm, (f) 0 dBm. In each subfigure, different colors represent different number of channels Nch according to the legend. Also, SNR evaluated after 0, 1600, 3200 and 4800 km transmission are shown in each subfigure where the top set of four curves in each subfigure is evaluated in back-to-back and the bottom set is after 4800 km transmission.
Fig. 7
Fig. 7 SNR of channel 9 versus Nch at PSIG = 0 dBm after 0, 1600, 3200 and 4800 km transmission at PLO of (a) 15.5 dBm, (b) 12 dBm, (c) 9 dBm, (d) 6 dBm, (e) 3 dBm, (f) 0 dBm.
Fig. 8
Fig. 8 SNR of channel 9 versus Nch at PSIG = −12 dBm after 0, 1600, 3200 and 4800 km transmission at PLO of (a) 15.5 dBm, (b) 12 dBm, (c) 9 dBm, (d) 6 dBm, (e) 3 dBm, (f) 0 dBm.
Fig. 9
Fig. 9 SNR of channel 9 versus PLO after 0, 1600, 3200 and 4800 km transmission at various Nch and at PSIG of: (a) 0 dBm, (b) −15 dBm.
Fig. 10
Fig. 10 Color coded 2-D plots of SNR versus both PSIG and PLO for (a) back-to-back and Nch = 1, (b) back-to-back and Nch = 16, (c) 1600 km transmission and Nch = 1, (d) 1600 km transmission and Nch = 16, (e) 3200 km transmission and Nch = 1, (f) 3200 km transmission and Nch = 16, (g) 4800 km transmission and Nch = 1, and (h) 4800 km transmission and Nch = 16.
Fig. 11
Fig. 11 SNR versus PSIG comparing both experimental data points with the points evaluated from the fitted SNR analytical model after 0, 1600, 3200 and 4800 km at the following operating conditions: (a) PLO = 12 dBm and Nch = 1, (b) PLO = 12 dBm and Nch = 16, (c) PLO = 3 dBm and Nch = 1, and (d) PLO = 3 dBm and Nch = 16.
Fig. 12
Fig. 12 SNR versus Nch comparing the experimental data points with the points evaluated from the fitted SNR analytical model after 0, 1600, 3200 and 4800 km at the following operating conditions: (a) PLO = 12 dBm and PSIG = 0 dBm, (b) PLO = 12 dBm and PSIG = −12 dBm, (c) PLO = 3 dBm and PSIG = 0 dBm, and (d) PLO = 3 dBm and PSIG = −12 dBm.
Fig. 13
Fig. 13 (a) Simulated intensity eyediagram of the single polarization QPSK signal using NRZ pulse shape, (b) experimental intensity eyediagram obtained after the single polarization IQ modulator, (c) scaling factor (β) versus the amount of residual CD for various polarization rotation angles (θ) obtained from simulations and theory on a 28 Gbaud PDM-QPSK signal.

Tables (1)

Tables Icon

Table 1 List of Symbols and Notation

Equations (69)

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

i P ( t ) = [ i = 1 N c h R S I G , i P ( E S I G , i ( t ) + E O N , i ( t ) ) + R L O , s P ( E L O ( t ) ) ] 2 + i s h P + i t h P 2 R S I G , s P R L O , s P [ E S I G , s ( t ) E L O ( t ) + E O N , s ( t ) E L O ( t ) ] + i = 1 N c h R S I G , i P [ E S I G , i 2 ( t ) + E O N , i 2 ( t ) + 2 E S I G , i ( t ) E O N , i ( t ) ] + R L O , s P E L O 2 ( t ) + i s h P + i t h P ,
i N ( t ) = [ i = 1 N c h R S I G , i N ( E S I G , i ( t δ ) + E O N , i ( t δ ) ) R L O , s N ( E L O ( t δ ) ) ] 2 + i s h N + i t h N 2 R S I G , s N R L O , s N [ E S I G , s ( t δ ) E L O ( t δ ) + E O N , s ( t δ ) E L O ( t δ ) ] + i = 1 N c h R S I G , i N [ E S I G , i 2 ( t δ ) + E O N , i 2 ( t δ ) + 2 E S I G , i ( t δ ) E O N , i ( t δ ) ] + R L O , s N E L O 2 ( t δ ) + i s h N + i t h N ,
E S I G , i ( t ) = R e { E ˜ S I G , i ( t ) e j 2 π ν S I G , i t } , E L O ( t ) = R e { E ˜ L O ( t ) e j 2 π ν L O t } , E O N , i ( t ) = R e { E ˜ O N , i ( t ) e j 2 π ν S I G , i t } ,
E 1 ( t ) E 2 ( t ) = R e { E ˜ 1 ( t ) e j 2 π ν 1 t } R e { E ˜ 2 ( t ) e j 2 π ν 2 t } = 1 2 R e { E ˜ 1 ( t ) E ˜ 2 * ( t ) e j 2 π ( ν 1 ν 2 ) t } + 1 2 R e { E ˜ 1 ( t ) E ˜ 2 ( t ) e j 2 π ( ν 1 + ν 2 ) t } ,
i P ( t ) = 2 R S I G , s P R L O , s P 1 2 R e { E ˜ S I G , s ( t ) E ˜ L O * ( t ) e j 2 π Δ ν t + E ˜ O N , s ( t ) E ˜ L O * ( t ) e j 2 π Δ ν t } + 1 2 i = 1 N c h R S I G , i P [ | E ˜ S I G , i ( t ) | 2 + | E ˜ O N , i ( t ) | 2 + 2 R e { E ˜ S I G , i ( t ) E ˜ O N , i * ( t ) } ] + R L O , s P 1 2 | E ˜ L O ( t ) | 2 + i s h P + i t h P ,
i N ( t ) = 2 R S I G , s N R L O , s N 1 2 R e { E ˜ S I G , s ( t δ ) E ˜ L O * ( t δ ) e j 2 π Δ ν ( t δ ) + E ˜ O N , s ( t δ ) E ˜ L O * ( t δ ) e j 2 π Δ ν ( t δ ) } + 1 2 i = 1 N c h R S I G , i N [ | E ˜ S I G , i ( t δ ) | 2 + | E ˜ O N , i ( t δ ) | 2 + 2 R e { E ˜ S I G , i ( t δ ) E ˜ O N , i * ( t δ ) } ] + R L O , s N 1 2 | E ˜ L O ( t δ ) | 2 + i s h N + i t h N ,
Δ i ( t ) = i P ( t ) h P ( t ) i N ( t ) h N ( t ) = Δ i L O S I G ( t ) + Δ i L O O N ( t ) + Δ i S I G S I G ( t ) + Δ i O N O N ( t ) + Δ i S I G O N ( t ) + Δ i L O L O ( t ) + Δ i s h ( t ) + Δ i t h ( t ) ,
Γ ¯ Δ i ( τ ) = Δ i ( t ) Δ i ( t τ ) ¯ = Γ ¯ Δ i L O S I G ( τ ) + Γ ¯ Δ i L O O N ( τ ) + Γ ¯ Δ i S I G S I G ( τ ) + Γ ¯ Δ i O N O N ( τ ) + Γ ¯ Δ i S I G O N ( τ ) + Γ ¯ Δ i L O L O ( τ ) + Γ ¯ Δ i s h ( τ ) + Γ ¯ Δ i t h ( τ ) ,
S Δ i ( f ) = FT  { Γ ¯ Δ i ( τ ) } = S Δ i L O S I G ( f ) + S Δ i L O O N ( f ) + S Δ i S I G S I G ( f ) + S Δ i O N O N ( f ) + S Δ i S I G O N ( f ) + S Δ i L O L O ( f ) + S Δ i s h ( f ) + S Δ i t h ( f ) S Δ i L O S I G ( f ) + S Δ i L O O N ( f ) + S Δ i S I G S I G ( f ) + S Δ i s h ( f ) + S Δ i t h ( f ) ,
S N R = S Δ i L O S I G ( f ) | H r x D S P ( f ) | 2 d f [ S Δ i L O O N ( f ) + S Δ i S I G S I G ( f ) + S Δ i s h ( f ) + S Δ i t h ( f ) ] A C | H r x D S P ( f ) | 2 d f = P Δ i L O S I G σ Δ i L O O N 2 + σ Δ i S I G S I G 2 + σ Δ i s h 2 + σ Δ i t h 2 ,
S N R = 2 P L O P S I G 2 c 1 P L O P O N + c 2 ( σ Δ i s h 2 + σ Δ i t h 2 ) + C M R R ¯ S I G N c h β P S I G 2 ,
c 1 = B c h / 2 B c h / 2 | H P ( f ) R S I G , s P R L O , s P + H N ( f ) e j 2 π f δ R S I G , s N R L O , s N | 2 | H r x D S P ( f ) | 2 d f / B c h | H P ( f ) R S I G , s P R L O , s P + H N ( f ) e j 2 π f δ R S I G , s N R L O , s N | 2 | G ( f ) | 2 | H r x D S P ( f ) | 2 d f / | G ( f ) | 2 d f ,
c 2 = 1 | H P ( f ) R S I G , s P R L O , s P + H N ( f ) e j 2 π f δ R S I G , s N R L O , s N | 2 | G ( f ) | 2 | H r x D S P ( f ) | 2 d f / | G ( f ) | 2 d f ,
C M R R ¯ S I G = | H P ( f ) R S I G P H N ( f ) e j 2 π f δ R S I G N | 2 S | E ˜ S I G | 2 A C ( f ) | H r x D S P ( f ) | 2 d f / S | E ˜ S I G | 2 A C ( f ) | H a v ( f ) | 2 | H r x D S P ( f ) | 2 d f | H P ( f ) R S I G P R L O P + H N ( f ) e j 2 π f δ R S I G N R L O N | 2 | G ( f ) | 2 | H r x D S P ( f ) | 2 d f / | G ( f ) | 2 d f ,
S N R = 2 2 c 1 O S N R 1 + C M R R S I G ¯ N c h β L S R 1 ,
S N R = 2 2 c 1 O S N R 1 + c 2 4 e R L O a v Δ f P S I G + c 2 i T I A 2 Δ f P L O P S I G ,
I D C = R L O P L O + N c h R S I G P S I G ,
I A C p p d = 8 R L O P L O P A P R R S I G P S I G ,
S N R = P L O P S I G α 1 P L O P S I G + α 2 N l o o p s P L O P S I G + α 3 + α 4 ( P L O + N c h P S I G ) + α 5 N c h P S I G 2 ,
α = [ α 1 α 5 ] ,     y = [ P L O 1 P S I G 1 / S N R 1 P L O m P S I G m / S N R m ] ,     X = [ P L O 1 P S I G 1 P L O m P S I G m N l o o p s 1 P L O 1 P S I G 1 N l o o p s m P L O m P S I G m 1 1 P L O 1 + N c h 1 P S I G 1 P L O m + N c h m P S I G m N c h 1 P S I G 1 2 N c h m P S I G m 2 ] ,
α ^ = ( X T X ) 1 X T y ,
Γ ¯ Δ i L O S I G ( τ ) = Δ i L O S I G ( t ) Δ i L O S I G ( t τ ) ¯ = R S I G , s P R L O , s P Re { E ˜ S I G , s ( t ) E ˜ L O * ( t ) e j 2 π Δ ν t } Re { E ˜ S I G , s ( t τ ) E ˜ L O * ( t τ ) e j 2 π Δ ν ( t τ ) } ¯ + R S I G , s N R L O , s N Re { E ˜ S I G , s ( t δ ) E ˜ L O * ( t δ ) e j 2 π Δ ν ( t δ ) } Re { E ˜ S I G , s ( t τ δ ) E ˜ L O * ( t τ δ ) e j 2 π Δ ν ( t τ δ ) } ¯ + R S I G , s P R L O , s P R S I G , s N R L O , s N Re { E ˜ S I G , s ( t ) E ˜ L O * ( t ) e j 2 π Δ ν t } Re { E ˜ S I G , s ( t τ δ ) E ˜ L O * ( t τ δ ) e j 2 π Δ ν ( t τ δ ) } ¯ + R S I G , s P R L O , s P R S I G , s N R L O , s N Re { E ˜ S I G , s ( t δ ) E ˜ L O * ( t δ ) e j 2 π Δ ν ( t δ ) } Re { E ˜ S I G , s ( t τ ) E ˜ L O * ( t τ ) e j 2 π Δ ν ( t τ ) } ¯ .
Γ ¯ Δ i L O S I G ( τ ) = 1 2 R S I G , s P R L O , s P Re { E ˜ S I G , s ( t ) E ˜ S I G , s * ( t τ ) E ˜ L O * ( t ) E ˜ L O ( t τ ) ¯ e j 2 π Δ ν τ } + 1 2 R S I G , s N R L O , s N Re { E ˜ S I G , s ( t δ ) E ˜ S I G , s * ( t τ δ ) E ˜ L O * ( t δ ) E ˜ L O ( t τ δ ) ¯ e j 2 π Δ ν τ } + 1 2 R S I G , s P R L O , s P R S I G , s N R L O , s N Re { E ˜ S I G , s ( t ) E ˜ S I G , s * ( t τ δ ) E ˜ L O * ( t ) E ˜ L O ( t τ δ ) ¯ e j 2 π Δ ν ( τ + δ ) } + 1 2 R S I G , s P R L O , s P R S I G , s N R L O , s N Re { E ˜ S I G , s ( t δ ) E ˜ S I G , s * ( t τ ) E ˜ L O * ( t δ ) E ˜ L O ( t τ ) ¯ e j 2 π Δ ν ( τ δ ) } .
Γ ¯ Δ i L O S I G ( τ ) = 1 2 R S I G , s P R L O , s P R e { Γ ¯ E ˜ S I G ( τ ) Γ E ˜ L O * ( τ ) e j 2 π Δ ν τ } + 1 2 R S I G , s N R L O , s N R e { Γ ¯ E ˜ S I G ( τ ) Γ E ˜ L O * ( τ ) e j 2 π Δ ν τ } + 1 2 R S I G , s P R L O , s P R S I G , s N R L O , s N R e { Γ ¯ E ˜ S I G ( τ + δ ) Γ E ˜ L O * ( τ + δ ) e j 2 π Δ ν ( τ + δ ) } + 1 2 R S I G , s P R L O , s P R S I G , s N R L O , s N R e { Γ ¯ E ˜ S I G ( τ δ ) Γ E ˜ L O * ( τ δ ) e j 2 π Δ ν ( τ δ ) } .
E ˜ S I G , s ( t δ ) E ˜ S I G , s * ( t τ δ ) E ˜ L O * ( t δ ) E ˜ L O ( t τ δ ) ¯ = E ˜ S I G , s ( t δ ) E ˜ S I G , s * ( t τ δ ) ¯ E ˜ L O * ( t δ ) E ˜ L O ( t τ δ ) = Γ E ˜ S I G ( t δ , τ ) ¯ Γ E ˜ L O * ( τ ) = Γ ¯ E ˜ S I G ( τ ) Γ E ˜ L O * ( τ ) ,
E ˜ S I G , s ( t ) E ˜ S I G , s * ( t τ δ ) E ˜ L O * ( t ) E ˜ L O ( t τ δ ) ¯ = E ˜ S I G , s ( t ) E ˜ S I G , s * ( t τ δ ) ¯ E ˜ L O * ( t ) E ˜ L O ( t τ δ ) = Γ E ˜ S I G ( t , τ + δ ) ¯ Γ E ˜ L O * ( τ + δ ) = Γ ¯ E ˜ S I G ( τ + δ ) Γ E ˜ L O * ( τ + δ ) ,
E ˜ S I G , s ( t δ ) E ˜ S I G , s * ( t τ ) E ˜ L O * ( t δ ) E ˜ L O ( t τ ) ¯ = E ˜ S I G , s ( t δ ) E ˜ S I G , s * ( t τ ) ¯ E ˜ L O * ( t δ ) E ˜ L O ( t τ ) = Γ E ˜ S I G ( t δ , τ δ ) ¯ Γ E ˜ L O * ( τ δ ) = Γ ¯ E ˜ S I G ( τ δ ) Γ E ˜ L O * ( τ δ ) ,
Γ ¯ Δ i L O S I G ( τ ) = P L O ( R S I G , s P R L O , s P + R S I G , s N R L O , s N ) R e { Γ ¯ E ˜ S I G ( τ ) e j 2 π Δ ν τ } + P L O R S I G , s P R L O , s P R S I G , s N R L O , s N . [ R e { Γ ¯ E ˜ S I G ( τ + δ ) e j 2 π Δ ν ( τ + δ ) + Γ ¯ E ˜ S I G ( τ δ ) e j 2 π Δ ν ( τ δ ) } ] .
S Δ i L O S I G ( f ) = FT { Γ ¯ Δ i L O S I G ( τ ) } = P L O | R S I G , s P R L O , s P + e j 2 π f δ R S I G , s N R L O , s N | 2 1 2 [ S E ˜ S I G ( f Δ ν ) + S E ˜ S I G ( ( f + Δ ν ) ) ] .
Γ ¯ Δ i L O S I G ( τ ) = P L O R S I G , s P R L O , s P ( R e { Γ ¯ E ˜ S I G ( τ ) e j 2 π Δ ν τ } Γ h P ( τ ) ) + P L O R S I G , s N R L O , s N ( R e { Γ ¯ E ˜ S I G ( τ ) e j 2 π Δ ν τ } Γ h N ( τ ) ) + P L O R S I G , s P R L O , s P R S I G , s N R L O , s N ( R e { Γ ¯ E ˜ S I G ( τ + δ ) e j 2 π Δ ν ( τ + δ ) } Γ h P h N ( τ ) ) + P L O R S I G , s P R L O , s P R S I G , s N R L O , s N ( R e { Γ ¯ E ˜ S I G ( τ δ ) e j 2 π Δ ν ( τ δ ) } Γ h N h P ( τ ) ) ,
Γ h P ( τ ) = h P ( τ ) h P * ( τ ) ,    Γ h N ( τ ) = h N ( τ ) h N * ( τ ) , Γ h P h N ( τ ) = h P ( τ ) h N * ( τ ) ,    Γ h N h P ( τ ) = h N ( τ ) h P * ( τ ) .
S Δ i L O S I G ( f ) = P L O | H P ( f ) R S I G , s P R L O , s P + e j 2 π f δ H N ( f ) R S I G , s N R L O , s N | 2 1 2 [ S E ˜ S I G ( f Δ ν ) + S E ˜ S I G ( ( f + Δ ν ) ) ] ,
E ˜ S I G ( t ) = n = ( a X n + b Y n ) h ( t n T ) ,
J = [ a b b * a * ] .
Γ ¯ E ˜ S I G ( τ ) = E ˜ S I G ( t ) E ˜ S I G * ( t τ ) ¯ = | a | 2 n = m = X n X m * h ( t n T ) h * ( t τ m T ) ¯ + | b | 2 n = m = Y n Y m * h ( t n T ) h * ( t τ m T ) ¯ + a b * n = m = X n Y m * h ( t n T ) h * ( t τ m T ) ¯ + a * b n = m = X n * Y m h ( t n T ) h * ( t τ m T ) ¯ .
X n X m * = Y n Y m * = { 1             n = m 0            n m ,
X n Y m * = X n * Y m = 0.
Γ ¯ E ˜ S I G ( τ ) = ( | a | 2 + | b | 2 ) n = h ( t n T ) h * ( t τ n T ) ¯ = 1 T n = T / 2 T / 2 h ( t n T ) h * ( t τ n T ) d t = 1 T h ( t ) h * ( t τ ) d t = 1 T Γ h ( τ ) ,
S E ˜ S I G ( f ) = 1 T | H ( f ) | 2 = 1 T | G ( f ) | 2 ,
P Δ i L O S I G = P L O | H P ( f ) R S I G , s P R L O , s P + H N ( f ) e j 2 π f δ R S I G , s N R L O , s N | 2 1 T | G ( f ) | 2 | H r x D S P ( f ) | 2 d f = 2 P L O P S I G | H P ( f ) R S I G , s P R L O , s P + H N ( f ) e j 2 π f δ R S I G , s N R L O , s N | 2 | G ( f ) | 2 | H r x D S P ( f ) | 2 d f | G ( f ) | 2 d f ,
S Δ i L O O N ( f ) = P L O | H P ( f ) R S I G , s P R L O , s P + H N ( f ) e j 2 π f δ R S I G , s N R L O , s N | 2 1 2 [ S E ˜ O N ( f Δ ν ) + S E ˜ O N ( ( f + Δ ν ) ) ] ,
S E ˜ O N ( f ) = { 2 N o           | f | B c h / 2 0                | f | > B c h / 2 ,
σ Δ i L O O N 2 = 2 P L O P O N B c h / 2 B c h / 2 | H P ( f ) R S I G , s P R L O , s P + H N ( f ) e j 2 π f δ R S I G , s N R L O , s N | 2 | H r x D S P ( f ) | 2 d f B c h .
Γ ¯ Δ i S I G S I G ( τ ) = Δ i S I G S I G ( t ) Δ i S I G S I G ( t τ ) ¯ = 1 4 [ i = 1 N c h j = 1 N c h R S I G , i P R S I G , j P | E ˜ S I G , i ( t ) | 2 | E ˜ S I G , j ( t τ ) | 2 ¯ + i = 1 N c h j = 1 N c h R S I G , i N R S I G , j N | E ˜ S I G , i ( t δ ) | 2 | E ˜ S I G , j ( t τ δ ) | 2 ¯ i = 1 N c h j = 1 N c h R S I G , i P R S I G , j N | E ˜ S I G , i ( t ) | 2 | E ˜ S I G , j ( t τ δ ) | 2 ¯ i = 1 N c h j = 1 N c h R S I G , i N R S I G , j P | E ˜ S I G , i ( t δ ) | 2 | E ˜ S I G , j ( t τ ) | 2 ¯ ] .
Γ ¯ Δ i S I G S I G ( τ ) = i = 1 N c h j = 1 j i N c h ( R S I G , i P R S I G , j P + R S I G , i N R S I G , j N 2 R S I G , i P R S I G , j N ) P S I G , i P S I G , j + 1 4 i = 1 N c h R S I G , i P 2 | E ˜ S I G , i ( t ) | 2 | E ˜ S I G , i ( t τ ) | 2 ¯ + 1 4 i = 1 N c h R S I G , i N 2 | E ˜ S I G , i ( t δ ) | 2 | E ˜ S I G , i ( t τ δ ) | 2 ¯ 1 4 i = 1 N c h R S I G , i P R S I G , i N | E ˜ S I G , i ( t ) | 2 | E ˜ S I G , i ( t τ δ ) | 2 ¯ 1 4 i = 1 N c h R S I G , i P R S I G , i N | E ˜ S I G , i ( t δ ) | 2 | E ˜ S I G , i ( t τ ) | 2 ¯ = i = 1 N c h j = 1 j i N c h ( R S I G , i P R S I G , j P + R S I G , i N R S I G , j N 2 R S I G , i P R S I G , j N ) P S I G , i P S I G , j + 1 4 i = 1 N c h R S I G , i P 2 Γ ¯ | E ˜ S I G , i | 2 ( τ ) + 1 4 i = 1 N c h R S I G , i N 2 Γ ¯ | E ˜ S I G , i | 2 ( τ ) 1 4 i = 1 N c h R S I G , i P R S I G , i N Γ ¯ | E ˜ S I G , i | 2 ( τ + δ ) 1 4 i = 1 N c h R S I G , i P R S I G , i N Γ ¯ | E ˜ S I G , i | 2 ( τ δ ) ,
S Δ i S I G S I G ( f ) = δ ( f ) i = 1 N c h j = 1 j i N c h ( R S I G , i P R S I G , j P + R S I G , i N R S I G , j N 2 R S I G , i P R S I G , j N ) P S I G , i P S I G , j + 1 4 i = 1 N c h | R S I G , i P e j 2 π f δ R S I G , i N | 2 S | E ˜ S I G , i | 2 ( f ) .
S Δ i S I G S I G ( f ) = δ ( f ) i = 1 N c h j = 1 j i N c h ( R S I G , i P R S I G , j P + R S I G , i N R S I G , j N 2 R S I G , i P R S I G , j N ) P S I G , i P S I G , j + 1 4 i = 1 N c h | H P ( f ) R S I G , i P H N ( f ) e j 2 π f δ R S I G , i N | 2 S | E ˜ S I G , i | 2 ( f ) ,
Γ ¯ | E ˜ S I G | 2 ( τ ) = E ˜ S I G ( t ) E ˜ S I G * ( t ) E ˜ S I G ( t τ ) E ˜ S I G * ( t τ ) ¯ = n = m = l = k = ( a X n + b Y n ) ( a X m + b Y m ) * ( a X l + b Y l ) ( a X k + b Y k ) * h ( t n T ) h * ( t m T ) h ( t τ l T ) h * ( t τ k T ) ¯ ,
( a X n + b Y n ) ( a X m + b Y m ) * ( a X l + b Y l ) ( a X k + b Y k ) * = | a | 4 X n X m * X l X k * + | b | 4 Y n Y m * Y l Y k * + | a | 2 | b | 2 X n X m * Y l Y k * + | a | 2 | b | 2 Y n Y m * X l X k * + | a | 2 | b | 2 Y n X m * X l Y k * + | a | 2 | b | 2 X n Y m * Y l X k * + terms with vanishing statistical averages,
X n X m * X l X k * = Y n Y m * Y l Y k * = { 1             n = m & l = k 1            n = k & m = l & n m 0           otherwise ,
X n X m * Y l Y k * = Y n Y m * X l X k * = { 1             n = m & l = k 0           otherwise ,
Y n X m * X l Y k * = X n Y m * Y l X k * = { 1             n = k & l = m 0           otherwise .
Γ ¯ | E ˜ S I G | 2 ( τ ) = ( | a | 4 + | b | 4 ) [ n = l = | h ( t n T ) | 2 | h ( t τ l T ) | 2 ¯ + n = m = n m h ( t n T ) h * ( t m T ) h ( t τ m T ) h * ( t τ n T ) ¯ ] + 2 | a | 2 | b | 2 [ n = l = | h ( t n T ) | 2 | h ( t τ l T ) | 2 ¯ + n = m = h ( t n T ) h * ( t m T ) h ( t τ m T ) h * ( t τ n T ) ¯ ] = 1 T ( | a | 4 + | b | 4 + 2 | a | 2 | b | 2 ) l = Γ | h | 2 ( τ l T ) + 1 T ( | a | 4 + | b | 4 ) k = k 0 C h ( τ , k ) + 1 T 2 | a | 2 | b | 2 k = C h ( τ , k ) ,
n = l = | h ( t n T ) | 2 | h ( t τ l T ) | 2 ¯ = 1 T n = l = T / 2 T / 2 | h ( t n T ) | 2 | h ( t τ l T ) | 2 d t = 1 T n = T / 2 T / 2 | h ( t n T ) | 2 l = | h ( t τ l T ) | 2 d t      ,let t - n T = z = 1 T n = n T T / 2 n T + T / 2 | h ( z ) | 2 l = | h ( z τ ( l n ) T ) | 2 d z = 1 T l = | h ( z ) | 2 | h ( z τ l T ) | 2 d z = 1 T l = Γ | h | 2 ( τ l T ) ,           where      Γ | h | 2 ( τ ) = | h ( t ) | 2 | h ( t τ ) | 2 d t
n = m = n m h ( t n T ) h * ( t m T ) h ( t τ m T ) h * ( t τ n T ) ¯ = 1 T n = m = n m T / 2 T / 2 h ( t n T ) h * ( t m T ) h ( t τ m T ) h * ( t τ n T ) d t = 1 T n = T / 2 T / 2 h ( t n T ) h * ( t τ n T ) [ m = n m h * ( t m T ) h ( t τ m T ) ] d t    ,let   t - n T = z = 1 T n = n T T / 2 n T + T / 2 h ( z ) h * ( z τ ) [ m = n m h * ( z ( m n ) T ) h ( z τ ( m n ) T ) ] d z = 1 T k = k 0 h ( t ) h * ( t τ ) h * ( t k T ) h ( t τ k T ) d t = 1 T k = k 0 C h ( τ , k ) ,          where   C h ( τ , k ) =   h ( t ) h * ( t τ ) h * ( t k T ) h ( t τ k T ) d t    
S | E ˜ S I G | 2 ( f ) = 1 T 2 | F T { h ( t ) h * ( t ) } | 2 l = δ ( f l T ) + 1 T k = | F T { h ( t ) h * ( t k T ) } | 2 1 T ( | a | 4 + | b | 4 ) | F T { h ( t ) h * ( t ) } | 2 = 1 T 2 | H ( f ) H * ( f ) | 2 l = δ ( f l T ) + 1 T k = | H ( f ) ( e j 2 π f k T H * ( f ) ) | 2 1 T ( | a | 4 + | b | 4 ) | H ( f ) H * ( f ) | 2 ,
FT { Γ | h | 2 ( τ ) } = | FT { h ( t ) h * ( t ) } | 2 = | H ( f ) H * ( f ) | 2 ,
FT { C h ( τ , k ) } =   | FT { h ( t ) h * ( t k T ) } | 2 = | H ( f ) ( e j 2 π f k T H * ( f ) ) | 2 .
S Δ i S I G S I G A C ( f ) = 1 4 i = 1 N c h | H P ( f ) R S I G , i P H N ( f ) e j 2 π f δ R S I G , i N | 2 S | E ˜ S I G , i | 2 A C ( f ) ,
S | E ˜ S I G , i | 2 A C ( f ) = 1 T k = | F T { h i ( t ) h i * ( t k T ) } | 2 1 T ( | a i | 4 + | b i | 4 ) | F T { h i ( t ) h i * ( t ) } | 2 ,
σ Δ i S I G S I G 2 = 1 4 i = 1 N c h | H P ( f ) R S I G , i P H N ( f ) e j 2 π f δ R S I G , i N | 2 S | E ˜ S I G , i | 2 A C ( f ) | H r x D S P ( f ) | 2 d f = 1 4 N c h σ | E ˜ S I G , s e | 2 2 | H P ( f ) R S I G P H N ( f ) e j 2 π f δ R S I G N | 2 S | E ˜ S I G | 2 A C ( f ) | H r x D S P ( f ) | 2 d f S | E ˜ S I G | 2 A C ( f ) | H a v ( f ) | 2 | H r x D S P ( f ) | 2 d f ,
σ Δ i s h 2 = 2 e [ ( R L O N + R L O P ) P L O + N c h ( R S I G N + R S I G P ) P S I G ] Δ f 4 e [ R L O a v P L O + N c h R S I G a v P S I G ] Δ f ,
σ Δ i t h 2 = ( i T I A ) 2 Δ f ,
S Δ i L O L O ( f ) = | H P ( f ) R L O P H N ( f ) e j 2 π f δ R L O N | 2 P L O 2 ( δ ( f ) + R I N ( f ) ) ,
S Δ i S I G O N ( f ) = 1 4 | H P ( f ) R S I G P H N ( f ) e j 2 π f δ R S I G N | 2 N c h ( S E ˜ S I G ( f ) S E ˜ O N ( f ) + S E ˜ S I G ( f ) S E ˜ O N ( f ) ) ,
S Δ i O N O N ( f ) = | H P ( f ) R S I G P H N ( f ) e j 2 π f δ R S I G N | 2 ( ( N o B o ) 2 δ ( f ) + N o 2 B o t r i ( f B o ) ) ,
R I N ( f ) = FT { δ P L O ( t ) δ P L O ( t τ ) / P L O 2 } ,
t r i ( f B o ) = { 1 | f | B o              | f | B o 0                      | f | > B o ,
β = 1 4 S | E ˜ S I G | 2 A C ( f ) | H a v ( f ) | 2 | H r x D S P ( f ) | 2 d f P S I G 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.