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

Impact of the carrier contribution factor in the self-coherent DC-value method

Open Access Open Access

Abstract

The carrier contribution factor (CCF) is used in the DC-value method to assure the minimum phase condition (MPC) at the transmitter. We analyze the impact of the CCF estimation accuracy at the receiver on the DC-value system performance. We found that the CCF estimation accuracy at the receiver should be within ±5% to exploit all advantages of the DC-value method. We propose an accurate method to estimate the CCF at the receiver that works effectively with both direct current (DC) and alternating current (AC) coupled photodetectors. We experimentally validate the proposed method employing a 24 GBaud 16QAM signal. The experimental results show that the minimum BER value can be approached within ±5% offset of the estimated CCF value almost independently of the carrier to signal power ratio (CSPR) value. The proposed method can also be used to reconstruct the DC component in other minimum phase signal self-coherent techniques, such as the Kramers-Kronig receiver.

© 2021 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Recently, several techniques have been proposed to extract amplitude and phase information from the optical field using direct detection (DD) schemes [13]. Among them, the DC-Value method offers upsampling free amplitude and phase recovery and works with low carrier to signal power ratio (CSPR) [35]. Nevertheless, the DC-Value performance is highly dependent on the receiver estimation of the carrier power that reaches the photodetector. Therefore, an accurate carrier contribution factor (CCF) estimation method is mandatory in order to exploit all advantages of the DC-Value method.

Single sideband (SSB) methods have been explored to retrieve full complex electric field of the optical signal using direct detection (DD) [1,6]. However, the performance of SSB transmission methods tend to be limited by the inherent signal-to-signal beating noise (SSBN) generated upon square-law detection [68]. To address this problem, more recently, a Kramers-Kronig (KK) receiver has been proposed that significantly alleviates the impact of SSBN and improves the system performance [1,9,10]. However, the nonlinear operations in the KK method demands higher sampling rates to accommodate spectral broadening [1113]. To address the high sampling rate problem, an upsampling free KK method was proposed that uses a Taylor expansion approximation of the associated nonlinear operations [2]. Nevertheless, the proposed approximations tend to require a higher CSPR, which results in a sensitivity penalty [2]. To address the aforementioned high sampling rate and high CSPR requirements of the conventional KK and upsampling free KK methods, respectively, a DC-Value method was recently proposed [3]. The DC-Value method provides an upsampling free phase reconstruction process at low CSPR [3]. The underlying idea of the DC-Value method is to iteratively impose the minimum phase condition (MPC) [3]. The signal reconstruction accuracy of the DC-Value method depends on the accurate estimation of CCF used in the MPC. The estimation of the CCF is simpler in the case of direct current coupled (DC-coupled) photodetectors. However, the DC component at the output of the photodetector causes unwanted bias offset to the subsequent circuits [14]. Also, the DC component degrades the signal quality when automatic gain control (AGC) is employed at the receiver [15,16]. The alternating current coupled (AC-coupled) photodetector alleviates the aforementioned problems and also allows to use full scale of the analog to digital convertor (ADC) that translates into improved signal to quantization noise ratio (SQNR). Therefore, the AC-coupled photodetector tends to become the preferred choice in the high-speed optical communication systems [16]. However, the AC-coupled photodetector makes the CCF estimation more challenging. Recently, a zero padded preamble based DC component recovery technique was proposed, which allows to recover CCF in AC-coupled photodetectors [16,17]. Nevertheless, it requires a preamble to estimate the CCF value [16,17]. Another alternative is to guess a constant bias and add it to the received signal. The guessed value can then be checked with the known CSPR value and be improved by successive iterations [18]. It means, however, that this method requires multiple iterations to estimate the DC component [18].

Here, we propose an accurate CCF estimation technique that does not require any iterative process and it is compatible with both DC and AC-coupled photodetectors. We also present an experimental validation of the proposed technique. We show that by applying this method the CCF estimation error falls below 1% for CSPR as low as 9 dB. Also, the experimental results confirm that the optimum BER value can be found within $\pm 5 \%$ offset of the estimated CCF value.

Besides the introduction, this paper comprises four more sections. In Section 2, we analyze the impact of CCF in the DC-Value method. In Section 3, we present a novel procedure to estimate the CCF. In Sections 4, we discuss experimental results obtained by the proposed method. In the last section, we present the major concluding remarks.

2. Impact of the carrier contribution factor

Consider the complex envelope of the incoming MPC optical signal as [3],

$$E(t) = E_{o} \; e^{i(2\pi f_{o}t+\phi)} + E_{s}(t),$$
where, $E_{o} \; e^{(i2\pi f_{o}t+\phi )}$ represents a carrier tone with an amplitude $E_{o}$, frequency $f_{o}$, and phase $\phi$ that coincides with the right edge of the information spectrum, and $E_{s}(t)$ is the complex information signal. When the field is detected using a DC-coupled photodetector, the following photocurrent will be generated,
$$I(t) = R \bigg(|E_{o}|^{2} + 2 \; \Re\{E_{s}(t) \; E_{o}^{*}e^{{-}i(2\pi f_{o}t+\phi)}\} + |E_{s}(t)|^{2} \bigg),$$
where, $R$ is the responsivity of the photodetector, and $\Re \{\cdot \}$ represents the real part of $\{\cdot \}$. The first term in Eq. (2) is the carrier-carrier beating DC component (square-root of the first term, $|E_{o}|$, is refereed as CCF), the second term is the desired carrier-signal beating, and the last term is the inherent SSBN. The working principle of the DC-Value method is based on iteratively imposing the SSB and DC-Value properties of the minimum phase signal, referred to as the MPC [3]. The MPC can be imposed on the Fourier transformed signal $\tilde{E'}_{n}(\omega )$, as it is done in [3], to attain the subsequent estimate of the minimum phase signal $\tilde{E}_n(\omega )$ in iteration $n$,
$$\tilde{E}_{n}(\omega) = \left\{\begin{array}{cl} 0, & \textrm{for}\; \mathrm{\omega} >0 \\ {N|E_{o}|}, & \textrm{for}\; \mathrm{\omega }= 0 \\ p\tilde{E'}_{n}(\omega), & \textrm{for}\; \mathrm{\omega} <0 \end{array}\right.$$
where $N$ represents the FFT length, and $p$ denotes the scaling factor, assuming the value $p = 2 \; \textrm{for} \; n=1 \; \textrm {and} \; {p = 1}$ $\textrm {for} \; n>1$.

2.1 Impact of the CCF

The normalized mean square error (NMSE) is an important parameter to characterize signal reconstruction accuracy in the DC-Value method. It is mathematically shown in [3] that the NMSE function between the known magnitude, $|E(t)|$, and the estimated magnitude, $|E_{n}(t)|$ (see Fig. 1),

$$\textrm{NMSE} = \frac{\int\limits_{-\infty}^{\infty} \Big | \; |\textrm{E(t)}| - |\textrm{E}_{\textrm{n}}(\textrm{t})| \; \Big |^{2} \textrm{dt}}{\int\limits_{-\infty}^{\infty} \Big | \textrm{E(t)} \Big | \textrm{dt} \int\limits_{-\infty}^{\infty} \Big | \textrm{E}_{\textrm{n}}(\textrm{t}) \Big | \textrm{dt}},$$
is a monotonically decreasing function after each iteration in the DC-Value method. The monotonically decreasing NMSE function shows the convergence of the DC-Value method, which leads to an accurate reconstruction of the minimum phase signal. An accurate reconstruction improves overall system performance in terms of BER and computation complexity as it requires less iterations in the DC-Value method, and propagates low distortions to the post-DSP. However, this result assumes that the receiver is able to retrieve the CCF value perfectly. In order to take into account the error in the CCF estimation at the receiver, the MPC can be re-written as,
$$\tilde{E}_n(\omega) = \begin{cases} \;\;\;0, & \textrm{for}\; \mathrm{\omega} > 0 \\ {N(|E_{o}| \pm \Delta E_{o})}, & \textrm{for} \; \mathrm{\omega} =0 \\ p\tilde{E'}_{n}(\omega), & \textrm{for} \; \mathrm{\omega } < 0 \end{cases}$$
where $\Delta E_{o}$ is the estimation error of the CCF. The estimation error $\Delta E_{o}$ can violate the MPC property and limits the signal reconstruction accuracy. Also, it should be noted that the signal reconstruction accuracy varies with the amount of CCF error estimation at the receiver end. Indeed, if the CCF estimation error is too high then no significant gain is obtained using the DC-Value method when compared with the other state of the art phase recovery methods. Indeed, this can justify the limited accuracy observed in the recently proposed experimental results for the time-domain implementation of the DC-Value method [19], where the DC-Value method provides a limited improvement when compared to KK methods.

 figure: Fig. 1.

Fig. 1. Schematic diagram of the NMSE calculation between available MPS magnitude, $|E(t)|$, and the DC-Value method estimated MPS magnitude, $|E_{n}(t)|$.

Download Full Size | PDF

Next, we quantify the impact of CCF estimation error $\Delta E_{o}$ in the signal reconstruction process using the DC-Value method. For the purpose of assessment of the CCF impact, let’s consider a 24 GBaud 16QAM DC-Value system with a CSPR of 12 dB. The electric field is detected using a single photodetector and full signal reconstructed using a DC-Value method algorithm. Results shown in Fig. 2 presents the NMSE between the known MPS magnitude, $|E(t)|$, and the DC-Value method estimated MPS magnitude, $|E_{n}(t)|$, as a function of $\Delta E_{o}$ for 20 iterations. In the signal reconstruction process, we varied the CCF error, $\Delta E_{o}$, from −25 to +25% and analyze its impact on the NMSE.

 figure: Fig. 2.

Fig. 2. Impact of CCF estimation error $\Delta E_{o}$ in signal reconstruction using the DC-Value method. Here, CSPR of 12 dB was considered in the numerical analysis.

Download Full Size | PDF

The analysis shows that for $\Delta E_{o} = 0 \%$, the NMSE converges to zero, as it was expected. We can also see that if $\Delta E_{o}$ increases then the NMSE becomes nearly independent of the number of iterations, preventing a proper convergence of the DC-Value method, thus resulting in a limited reconstruction accuracy. Results shown in Fig. 3 confirm that the NMSE is monotonically decreasing for the $\Delta E_{o} = 0\%$. For an estimation error greater than $\pm 2.5 \%$, a 17 dB NMSE penalty is observed in the numerical analysis for 5 iterations. Results presented in Fig. 4 show that for an SNR of 30 dB, a BER threshold of $10^{-3}$ can be assured after 10 iterations if CCF estimation error $\Delta E_{o}$ falls within approximately $\pm 5 \%$ range. Further, we have performed extensive simulation for other similar systems and the obtained results are analogous, i.e. the CCF must be estimated within a $\pm 5\%$ error to fully exploit the advantages of the DC-Value method. In the subsequent section, we show an accurate method to estimate CCF for the DC-Value method compatible with both DC and AC-coupled photodetectors.

 figure: Fig. 3.

Fig. 3. Impact of CCF estimation error $\Delta E_{o}$ on the NMSE as a function of iteration number.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. Impact of CCF estimation error $\Delta E_{o}$ on the BER as a function of SNR. For higher negative value of $\Delta E_{o}$ (usually less than −5%), the CCF value reduces significantly and the signal no longer follows the minimum phase condition in the phase reconstruction process which results in highly degraded system performance. Therefore, we have considered only positive $\Delta E_{o}$ values in the analysis to assess the performance.

Download Full Size | PDF

3. Carrier contribution factor estimation method

In a DC-coupled photodetector, the photocurrent generated upon direct-detection can be written from Eq. (2) as (considering responsivity $R=1$),

$$I(t) = |E_{o}|^{2} + 2 \; \Re\{E_{s}(t) \; E_{o}^{*}e^{{-}i(2\pi f_{o}t+\phi)}\} + |E_{s}(t)|^{2}.$$

The DC-coupled photocurrent contains nonzero DC component contributed from the carrier-carrier beating, $|E_{o}|^{2}$, and the signal-signal beating, $|E_{s}(t)|^{2}$. The carrier-signal beating does not contain any non-zero DC component as it has zero mean value. Thus, the DC component in Eq. (6) can be written as,

$$\langle I_{}(t) \rangle = |E_{o}(t)|^{2} + \langle |E_{s}(t)|^{2} \rangle,$$
therefore,
$$\langle I_{}(t) \rangle = P_{c} + \langle P_{s}(t) \rangle,$$
where $P_{c}=|E_{o}|^{2}$, $P_{s}(t)=|E_{s}(t)|^{2}$ and $\langle \cdot \rangle$ indicates an average operator. To make the notations simpler in the CCF estimation process, we can re-write Eq. (6) as,
$$I(t) = P_{c} + P_{s}(t) + 2\sqrt{ P_{c} P_{s}(t)} \; \textrm{cos}(\phi).$$

Notice that the photodetector sensitivity is assumed as $R=1$ for simplicity. In practice, the value of responsivity may differ from 1, nevertheless, the proposed CCF estimation method is insensitive to responsivity value.

In the case of AC-coupled photodetector, the DC component given by Eq. (8) gets filtered out by the DC block. Therefore, for the AC-coupled photodetector the current is given by,

$$I_{AC}(t)=I(t)-\langle I(t) \rangle.$$

From Eq. (8) and (9), we can write $I_{AC}(t)$ as,

$$I_{AC}(t) = P_{s}(t) + 2\sqrt{ P_{c} P_{s}(t)} \; \textrm{cos}(\phi) - \langle P_{s}(t) \rangle.$$

In case of sufficiently large CSPR (usually higher than 8 dB), we can assume that (refer to Fig. 5 for the approximation error),

$$\langle I_{AC}^{2}(t) \rangle \approx 2 P_{c} \langle P_{s}(t) \rangle.$$

Now considering $\langle P_{s}(t) \rangle = P_{c}/\textrm {CSPR}$ in the Eq. (12),

$$\langle I_{AC}^{2}(t) \rangle \approx 2 P_{c} \langle P_{s}(t) \rangle = 2 P_{c}^{2}/\textrm{CSPR}.$$

Therefore,

$$P_{c}^{2} = \langle I_{AC}^{2}(t) \rangle \dfrac{\textrm{CSPR}}{2}, $$
$$P_{c} = \sqrt{\langle I_{AC}^{2}(t) \rangle \dfrac{\textrm{CSPR}}{2}}.$$

 figure: Fig. 5.

Fig. 5. Approximation error in $\langle I^{2}_{AC}(t) \rangle$ when Eq. (12) is used.

Download Full Size | PDF

In a similar way, considering $P_{c} = \textrm {CSPR} \; \langle P_{s}(t) \rangle$ in Eq. (12), we obtain,

$$\langle P_{s}(t) \rangle = \sqrt{\dfrac{\langle I_{AC}^{2}(t) \rangle}{2 \; \textrm{CSPR}}}.$$

Therefore, from Eq. (11), (14) and (15), the photocurrent with its recovered DC component can be written as,

$$I(t) \simeq I_{AC}(t) + \sqrt{\langle I_{AC}^{2}(t) \rangle \dfrac{\textrm{CSPR}}{2}} (1+1/\textrm{CSPR}).$$

The CCF estimation error in percentage can be written as,

$$\% \Delta E_{o}= \frac{|E_{o}| - \hat{E}_{o}}{|E_{o}|} \times 100,$$
where $\hat {E}_{o}$ represents the estimated CCF from Eq. (14), and it can be written as,
$$\hat{E}_{o} = \sqrt{P_{c}} \simeq \bigg({{\langle I_{AC}^{2}(t) \rangle \dfrac{\textrm{CSPR}}{2}}}\bigg)^{\dfrac{1}{4}}.$$

Finally, the Eq. (18) can be used to formulate the MPC in the DC-Value method.

The results in Fig. 6 show the numerical analysis of CCF estimation accuracy of the purposed method as function of the CSPR. The analysis shows that the CSPR is not a relevant factor for the CCF estimation as the $\Delta E_{o}$ lies below 3.5% for the given CSPR range. The results in the Fig. 6 shows that less than 1% error can be ensured for the CSPR value as low as 8.5 dB for the given QPSK, 16QAM, and 64QAM modulation formats. The algorithmic flow shown in Fig. 7 presents the proposed method of DC component recovery and estimation of CCF. In practice, the known value of CSPR, and the photodetected current signal are supplied to the CCF estimation algorithm. First, the signal is passed through the DC remover and followed by a DC component recovery estimation (as per Eq. (16)) to calculate the photocurrent with its recovered DC component. Similarly, by using the known value of CSPR and output of the DC remover, the estimation of CCF can be carried out (as per Eq. (18)) before employing the DC-Value method.

 figure: Fig. 6.

Fig. 6. Numerical analysis of CCF estimation error $\Delta E_{o}$ in % using Eq. (17) with respect to CSPR for different modulation formats. The analysis shows that the CSPR is not a relevant factor for the CCF estimation as the $\Delta E_{o}$ lies below 3.5% for the given CSPR range.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. The algorithmic flow of DC component recovery and estimation of CCF compatible with both DC and AC-coupled photodetectors, i.e. $I_{AC}(t)$ or $I_{DC}(t)$ as an input.

Download Full Size | PDF

4. Experimental results

In this section, we experimentally validate the proposed CCF estimation method and its impact on the overall system performance. We start this section by describing the used experimental setup in detail.

4.1 DC-value transceiver experimental setup

4.1.1 Tx-DSP

At the transmitter side, see Fig. 8(a), a single-channel 30 Gbaud QPSK and 24 Gbaud 16QAM signals are generated using a root-raised cosine (RRC) pulse shaping filter with 0.1 roll-off factor, respectively. Following that, a complex tone $E_{o}e^{2 \pi f_{o}t}$, at frequency $f_{o}=16.5 \;\textrm{GHz}$ (considering 0.1 RRC roll-off factor, $\textrm{1}.1 \times \; 30\; \textrm{Gbaud}/2~=~16.5 \; \textrm{GHz}$), is added to the 30 Gbaud QPSK signal, $E_{s}(t)$, in the digital domain (similarly, complex tone at $f_{o}=13.2 \;\textrm{GHz}$ is added to the 24 Gbaud 16QAM signal). The digitally added complex tone coincides with the right edge of the information signal spectrum to generate an SSB information signal. The amplitude $E_{o}$ of the complex tone should be kept sufficiently high (usually higher than 9 to 10 dB) to satisfy the MPC condition upon photodetection. The desired CSPR is obtained by varying the amplitude $E_{o}$ of the complex tone.

 figure: Fig. 8.

Fig. 8. Experimental setup of a self-coherent transceiver employing DC-Value method. Here, virtual carrier added in the digital domain to reduce optical complexity. To reduce cost, only one optical amplifier employed (as a preamplifier) in the transmission link.

Download Full Size | PDF

4.1.2 Setup

After Tx-DSP, in our case, a Keysight M8194A arbitrary waveform generator (AWG) containing 120 GSa/s digital to analog converter (DAC) is used to generate the signal. It should be noted that no pre-equalization was employed at the transmitter end. The output of the DAC is modulated at $1550 \; \textrm{nm}$ using a single polarization IQ Mach-Zehnder modulator (IQ-MZM). The output of the IQ-MZM signal is directly launched to the standard single-mode fiber (SSMF). It should be noticed that the digitally added carrier tone could make the IQ-MZM bias control circuitry unstable due to its very high amplitude and make it difficult to lock IQ-MZM bias conditions. To circumvent this problem, first the bias conditions are locked using only the information signal without digitally added carrier tone and then after locking bias conditions, we can digitally add carrier tone to generate the SSB signal. Following SSMF, the optical signal is amplified by an Erbium-doped fiber amplifier (EDFA) and detected using a low-cost $40 \; \textrm{GHz}$ single photodetector without trans-impedance amplifier (TIA). The photodetector output is then sampled by a $100 \; \textrm{GSa}/s$ Tektronix real-time oscilloscope (RTO) with $33 \; \textrm{GHz}$ bandwidth. Finally, the experimental data captured by the RTO is processed offline using Matlab software.

4.1.3 Rx-DSP

In the RX-DSP, first, pre-equalization is carried to compensate the non-flat response of the photodetector. Following that, the proposed CCF estimation technique is applied to estimate the CCF and recover the DC component of the photodetected signal. Next, the signal is passed to the DC-Value method to recover the missing phase information. Notice that the pre-equalization at the receiver end includes the inversion of $S_{21}$ response of the photodetector. Moreover, the transmitter response is compensated by an adaptive equalizer after the signal reconstruction process. Also, it should be noted that no digital upsampling or downsampling is required in the DC-Value method.

In the DC-Value method, first, a square-root operation is carried out to obtain the magnitude $|E(t)|$ of the optical field and multiplied by a phase correction factor $e^{i\delta \theta _{n-1}(t)}$ which outputs the complex signal $E'_{n}(t)$. In the first iteration ($n=1$), the phase correction vector is assumed to be zero, i.e. $\delta \theta _0(t) = 0$. Next, the MPC is imposed (see Fig. 8(c) on the Fourier transformed signal $\tilde{E'}_{n}(\omega )$ to attain $\tilde{E}_n(\omega )$ as discussed in Eq. (3). The scaling factor $p$ used in the MPC greatly speeds up the convergence process (see Fig. 8(d)). It is worth noticing that the scaling factor $p$ improves the convergence speed because the input signal, $E'_n(t)$, for the first iteration consists of a real-valued amplitude signal $|E(t)|$, i.e. $\delta \theta _0(t) = 0$ [3]. The minimum phase condition imposed on the signal forces the positive frequency components to zero and it generates complex-valued signal $E_1(t)$ in the time domain (starts acquiring the phase information). As a consequence of forcing positive frequency components to zero, the amplitude of the signal $E_1(t)$ will have its real and imaginary parts scaled by a factor 0.5. Therefore, we set a scaling factor $p=2$ to adjust the amplitude of the information signal in the first iteration, which speeds up the convergence of the method (see Fig. 8(d)) [3]. It shows that the implementation of the scaling factor $p$ requires two less iterations to achieve the same accuracy without scaling factor. After imposing MPC, the IFFT of the $\tilde{E}_n(\omega )$ is computed to obtain the first estimate of the minimum phase signal $E_{1}(t)$. The phase correction vector corresponds to $E_{1}(t)$ and can be calculated as $e^{i \delta \theta _{1}(t)} = \frac {E_{1}(t)}{|E_{1}(t)|}$. This phase estimated is then supplied as an updated phase correction vector for the subsequent iteration ($n=2$). This process continues until the desired accuracy is reached. Also, the normalized mean squared error (NMSE) between the known magnitude $|E(t)|$ and the estimated magnitude $|E_{n}(t)|$ is monotonically decreasing after each iteration (see Fig. 8(d)), therefore, the reconstruction process converges as discussed earlier in Section 2. Afterwards, the reconstructed minimum phase signal $E_{n}(t)$ is passed through the DC remover and then upconverted by multiplying it with $e^{i2\pi f_{o}t}$. Finally, the chromatic dispersion compensation and equalization can be carried out before symbol recovery.

4.2 Performance analysis

In this subsection, we present the experimental analysis of the impact of the CCF (estimated by proposed method, and referred as $\hat {E}_{o}$) on the overall system BER performance in the DC-Value method. To analyze the impact of $\hat {E}_{o}$, we first calculate the $\hat {E}_{o}$ value as discussed previously in the Section 3., and add an offset error $\delta$ in the estimated $\hat {E}_{o}$ ranging from −25 to +25% with the step size of 0.55 % to assess its impact on the overall system performance. The $\hat {E}_{o,\textrm {offset}}$ can be given as,

$$\hat{E}_{o,\textrm{offset}} = \hat{E}_{o} \bigg(1 + \dfrac{\delta}{100}\bigg).$$

Results shown in Fig. 9 present the impact of $\delta$ on the overall BER performance for the different CSPR values for 24 GBaud 16QAM signal after 40 km SSMF. The dotted red line in Fig. 9 represents the estimated $\hat {E}_{o}$ value by the proposed CCF estimation method (with no added offset, i.e. $\delta =0 \%$). The surface between two dotted yellow lines represent the range between $\pm 5 \%$ offset. For CSPR values ranging from 7 to 15 dB, we change the offset value $\delta$ as discussed above and calculate the BER performance of the system. For each CSPR value, we find the lowest achievable BER in the given range of $\delta$ and it is represented as blue dots in the graph. Also, small green dots shown in the results represent valid BER values considering the soft decision forward error correction (SD-FEC) threshold of $2.4 \times 10^{-2}$. It shows that the value of BER becomes valid for the CSPR $\geq 9$ dB considering the same SD-FEC threshold. The major highlight of Fig. 9 is that the lowest BER lies between $\pm 5 \%$ of the estimated $\hat {E}_{o}$ value which shows the higher CCF estimation accuracy of the proposed method. Nevertheless, Fig. 9 shows that the estimation error $\delta$ increases for CSPR higher than $\sim$ 12 dB. This peculiar behaviour comes from the Mach-Zehnder modulator as it starts operating in a nonlinear regime, which effectively changes the CSPR at the transmitter end. The proposed method requires the knowledge of the CSPR at the receiver end to recover the DC component. Additionally, it should be noted that the actual CSPR may change due to the impairments caused by transmission and optoelectronic components. Therefore, starting from the known value of CSPR at the transmitter, the DC-Value method can be used as a tool to improve the CSPR estimation at the receiver by proper adjustments of the offset $\delta$ such that the NMSE becomes a monotonically decreasing function as we discussed earlier in Fig. 3.

 figure: Fig. 9.

Fig. 9. Performance analysis in terms of BER for the 24 GBaud 16QAM signal after 40 km SSMF as function of CSPR and estimated $\hat {E}_{o}$ offset $\delta$.

Download Full Size | PDF

Results shown in Fig. 10(a) presents the $\log_{10}(BER)$ as a function of $\delta$ for the 24 GBaud 16 QAM signal received after 40 km SSMF. The shaded gray area in Fig. 10(a) represents the $\pm 5 \%$ range of $\delta$. It is shown that the optimum BER value can be found in the close proximity of the estimated $\hat {E}_{o}$ value by the proposed CCF estimation method. Additionally, $\sim$ 0.2 dB penalty in terms of Q factor is observed for the $\delta \geq$ 10 % for CSPR of 11 dB. Similarly, the analysis of the received signal after 70 km SSMF shown in Fig. 10(b) represents the identical result like Fig. 10(a) such that the optimum BER values for the given CSPR lie within $\pm 5\%$ range of $\delta$. Another important observation is that after 13 dB CSPR value, an optimum value of $\log_{10}(BER)$ is found a little bit far from the $\hat {E}_{o}$ value. The possible reason behind this could be that the digitally added carrier tone makes the IQ-MZM bias controller unstable when a very high power tone (usually higher than 13 dB observed in our experiment) is inserted in the signal. In this situation, the actual CSPR value of the optical signal differs from the virtual carrier added digital signal. Also, notice that slightly higher CSPR can help reduce the phase jump [20].

 figure: Fig. 10.

Fig. 10. Impact of $\hat {E}_{o}$ estimation error on the system performance after (a) 40 km, and (b) 70 km. Here, virtual carrier assisted 24 GBaud 16 QAM signal used for the performance assessment. Results shows that the lowest BER can be assured in the DC-Value method when the $\hat {E}_{o}$ estimation falls below $\pm 5\%$.

Download Full Size | PDF

5. Conclusion

The CCF is an important parameter in the DC-Value method and it used to form the MPC for the DC-Value method algorithm. Using the knowledge of CSPR value, we proposed a novel CCF estimation method compatible with both DC and AC-coupled photodetector for the DC-Value self-coherent transceiver. We present the analysis and impact of the CCF in the DC-Value method and show that the CCF estimation error should be no more than $\pm 5\%$ to ensure better system performance. The numerical analysis shows that the CCF estimation is quite insensitive to the CSPR value. In practice, results shows that the optimum CSPR value usually lies around 12 to 13 dB and we can ensure less than 1% CCF estimation error by the proposed method. Also, we performed experimental analysis of 24 GBaud 16QAM signal transmission and employed proposed method to estimate the CCF.

Funding

Fundação para a Ciência e a Tecnologia (DSPMetroNet (POCI-01-0145-FEDER-029405), ORCIP (CENTRO-01- 0145-FEDER-022141), UIDB/50008/2020-UIDP/50008/2020 (action DigCORE)).

Acknowledgments

This work was supported in part by Fundação para a Ciência e a Tecnologia (FCT) through national funds, by the European Regional Development Fund (FEDER), through the Competitiveness and Internationalization Operational Programme (COMPETE 2020) of the Portugal 2020 framework, under the projects DSPMetroNet (POCI-01-0145-FEDER-029405), UIDB/50008/2020-UIDP/50008/2020 (action DigCORE), and ORCIP (CENTRO-01-0145-FEDER-022141). Fernando P. Guiomar acknowledges a fellowship from "la Caixa" Foundation (ID 100010434). The fellowship code is LCF/BQ/PR20/11770015. Marco Fernandes acknowledges a PhD fellowship from FCT. The fellowship code is 2020.07521.BD.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. A. Mecozzi, C. Antonelli, and M. Shtaif, “Kramers-Kronig Coherent Receiver,” Optica 3(11), 1220–1227 (2016). [CrossRef]  

2. T. Bo and H. Kim, “Towards Practical Kramers-Kronig Receiver: Resampling, Performance, and Implementation,” J. Lightwave Technol. 37, 461–469 (2018). [CrossRef]  

3. R. K. Patel, I. A. Alimi, N. J. Muga, and A. N. Pinto, “Optical Signal Phase Retrieval With Low Complexity DC-Value Method,” J. Lightwave Technol. 38(16), 4205–4212 (2020). [CrossRef]  

4. R. Patel, F. Guiomar, M. Fernandes, I. Alimi, P. Monteiro, N. Muga, and A. Pinto, “Virtual Carrier Assisted Self-Coherent Detection Employing DC-Value Method,” in 2021 Optical Fiber Communications Conference and Exhibition (OFC), (2021), pp. 1–3.

5. N. J. Muga, R. K. Patel, I. A. Alimi, N. A. Silva, and A. N. Pinto, “DSP Optimization for Simplified Coherent Receivers,” in 2020 22nd International Conference on Transparent Optical Networks (ICTON), (2020), pp. 1–4.

6. X. Chen, C. Antonelli, S. Chandrasekhar, G. Raybon, A. Mecozzi, M. Shtaif, and P. Winzer, “Kramers–kronig receivers for 100-km datacenter interconnects,” J. Lightwave Technol. 36(1), 79–89 (2018). [CrossRef]  

7. S. Randel, D. Pilori, S. Chandrasekhar, G. Raybon, and P. Winzer, “100-Gb/s Discrete-Multitone Transmission Over 80-km SSMF Using Single-Sideband Modulation with Novel Interference-Cancellation Scheme,” 2015 European Conference on Optical Communication (ECOC) pp. 1–3 (2015).

8. Z. Li, M. Erkilinç, S. Pachnicke, H. Griesser, B. Thomsen, P. Bayvel, and R. Killey, “Direct-Detection 16-QAM Nyquist-Shaped Subcarrier Modulation with SSBI Mitigation,” IEEE International Conference on Communications pp. 5204–5209 (2015).

9. A. Mecozzi, “Retrieving the Full Optical Response from Amplitude Data by Hilbert Transform,” Opt. Commun. 282(20), 4183–4187 (2009). [CrossRef]  

10. D. Cassioli and A. Mecozzi, “Minimum-phase impulse response channels,” IEEE Trans. Commun. 57(12), 3529–3532 (2009). [CrossRef]  

11. T. Bo, “Kramers-Kronig Receiver Without Digital Upsampling,” in 2018 Optical Fiber Communications Conference and Exposition (OFC), (2018), pp. 1–3.

12. T. Bo and H. Kim, “Kramers-Kronig Receiver Operable Without Digital Upsampling,” Opt. Express 26(11), 13810–13818 (2018). [CrossRef]  

13. N. J. Muga, R. K. Patel, I. A. Alimi, N. A. Silva, and A. N. Pinto, “Self-coherent optical detection for access and metro networks,” in 2019 21st International Conference on Transparent Optical Networks (ICTON), (2019), pp. 1–4.

14. K. Hara, S. Kimura, H. Nakamura, N. Yoshimoto, and H. Hadama, “New AC-Coupled Burst-Mode Optical Receiver Using Transient-Phenomena Cancellation Techniques for 10 Gbit/s-Class High-Speed TDM-PON Systems,” J. Lightwave Technol. 28(19), 2775–2782 (2010). [CrossRef]  

15. E. Sackinger, Broadband Circuits for Optical Fiber Communication (John Wiley & Sons, Inc., 2005).

16. T. Bo and H. Kim, “DC Component Recovery in Kramers-Kronig Receiver Utilizing AC-Coupled Photo-Detector,” J. Lightwave Technol. 38(16), 4307–4314 (2020). [CrossRef]  

17. S. An, Q. Zhu, J. Li, and Y. Su, “Modified KK Receiver with Accurate Field Reconstruction at Low CSPR Condition,” in 2019 Optical Fiber Communications Conference and Exhibition (OFC), (2019), pp. 1–3.

18. A. Mecozzi, C. Antonelli, and M. Shtaif, “Kramers-Kronig Receivers,” Adv. Opt. Photonics 11(3), 480–517 (2019). [CrossRef]  

19. W. Wang, D. Zou, Z. Li, Q. Sui, Z. Cao, C. Lu, F. Li, and Z. Li, “Optical Single Sideband Signal Reconstruction Based on Time-Domain Iteration,” J. Lightwave Technol. 39, 2319–2326 (2021). [CrossRef]  

20. S. T. Le, K. Schuh, M. Chagnon, F. Buchali, and H. Buelow, “1.6 Tb/s Virtual-Carrier Assisted WDM Direct Detection Transmission Over 1200 km,” J. Lightwave Technol. 37(2), 418–424 (2019). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

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

Fig. 1.
Fig. 1. Schematic diagram of the NMSE calculation between available MPS magnitude, $|E(t)|$, and the DC-Value method estimated MPS magnitude, $|E_{n}(t)|$.
Fig. 2.
Fig. 2. Impact of CCF estimation error $\Delta E_{o}$ in signal reconstruction using the DC-Value method. Here, CSPR of 12 dB was considered in the numerical analysis.
Fig. 3.
Fig. 3. Impact of CCF estimation error $\Delta E_{o}$ on the NMSE as a function of iteration number.
Fig. 4.
Fig. 4. Impact of CCF estimation error $\Delta E_{o}$ on the BER as a function of SNR. For higher negative value of $\Delta E_{o}$ (usually less than −5%), the CCF value reduces significantly and the signal no longer follows the minimum phase condition in the phase reconstruction process which results in highly degraded system performance. Therefore, we have considered only positive $\Delta E_{o}$ values in the analysis to assess the performance.
Fig. 5.
Fig. 5. Approximation error in $\langle I^{2}_{AC}(t) \rangle$ when Eq. (12) is used.
Fig. 6.
Fig. 6. Numerical analysis of CCF estimation error $\Delta E_{o}$ in % using Eq. (17) with respect to CSPR for different modulation formats. The analysis shows that the CSPR is not a relevant factor for the CCF estimation as the $\Delta E_{o}$ lies below 3.5% for the given CSPR range.
Fig. 7.
Fig. 7. The algorithmic flow of DC component recovery and estimation of CCF compatible with both DC and AC-coupled photodetectors, i.e. $I_{AC}(t)$ or $I_{DC}(t)$ as an input.
Fig. 8.
Fig. 8. Experimental setup of a self-coherent transceiver employing DC-Value method. Here, virtual carrier added in the digital domain to reduce optical complexity. To reduce cost, only one optical amplifier employed (as a preamplifier) in the transmission link.
Fig. 9.
Fig. 9. Performance analysis in terms of BER for the 24 GBaud 16QAM signal after 40 km SSMF as function of CSPR and estimated $\hat {E}_{o}$ offset $\delta$.
Fig. 10.
Fig. 10. Impact of $\hat {E}_{o}$ estimation error on the system performance after (a) 40 km, and (b) 70 km. Here, virtual carrier assisted 24 GBaud 16 QAM signal used for the performance assessment. Results shows that the lowest BER can be assured in the DC-Value method when the $\hat {E}_{o}$ estimation falls below $\pm 5\%$.

Equations (20)

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

E ( t ) = E o e i ( 2 π f o t + ϕ ) + E s ( t ) ,
I ( t ) = R ( | E o | 2 + 2 { E s ( t ) E o e i ( 2 π f o t + ϕ ) } + | E s ( t ) | 2 ) ,
E ~ n ( ω ) = { 0 , for ω > 0 N | E o | , for ω = 0 p E ~ n ( ω ) , for ω < 0
NMSE = | | E(t) | | E n ( t ) | | 2 dt | E(t) | dt | E n ( t ) | dt ,
E ~ n ( ω ) = { 0 , for ω > 0 N ( | E o | ± Δ E o ) , for ω = 0 p E ~ n ( ω ) , for ω < 0
I ( t ) = | E o | 2 + 2 { E s ( t ) E o e i ( 2 π f o t + ϕ ) } + | E s ( t ) | 2 .
I ( t ) = | E o ( t ) | 2 + | E s ( t ) | 2 ,
I ( t ) = P c + P s ( t ) ,
I ( t ) = P c + P s ( t ) + 2 P c P s ( t ) cos ( ϕ ) .
I A C ( t ) = I ( t ) I ( t ) .
I A C ( t ) = P s ( t ) + 2 P c P s ( t ) cos ( ϕ ) P s ( t ) .
I A C 2 ( t ) 2 P c P s ( t ) .
I A C 2 ( t ) 2 P c P s ( t ) = 2 P c 2 / CSPR .
P c 2 = I A C 2 ( t ) CSPR 2 ,
P c = I A C 2 ( t ) CSPR 2 .
P s ( t ) = I A C 2 ( t ) 2 CSPR .
I ( t ) I A C ( t ) + I A C 2 ( t ) CSPR 2 ( 1 + 1 / CSPR ) .
% Δ E o = | E o | E ^ o | E o | × 100 ,
E ^ o = P c ( I A C 2 ( t ) CSPR 2 ) 1 4 .
E ^ o , offset = E ^ o ( 1 + δ 100 ) .
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.