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

Paired-subcarrier equalization for phase noise and transmitter IQ skew in digital subcarrier multiplexing system

Open Access Open Access

Abstract

The digital subcarrier multiplexing (DSCM) transmission scheme is expected for future ultra-large baud rate transmission. However, the phase noise and transmitter (Tx) IQ skew tolerance are decreased due to the narrow sub-band transmission and conjugated interference from symmetric subcarrier. In this paper, we propose a paired-subcarrier equalization scheme to jointly mitigate the phase noise and Tx IQ skew. We use a phase locking loop (PLL) embedded 4 × 4 MIMO equalizer to simultaneously realize polarization demultiplexing, phase noise and Tx IQ skew compensation. The 4 × 4 MIMO can deal with the paired-subcarrier interference in the DSCM transmission. Besides, since the inner subcarrier suffers smaller interference from its symmetric subcarrier, we estimate the phase noise by inner subcarriers and share the phase noise information with other subcarriers to reduce the overall complexity. Through simulations of 100-GBaud 64-QAM DSCM coherent optical fiber transmission consisting of eight 12.5-Gbaud subcarriers and experiment of 10-GBaud four-subcarriers PM-16QAM transmission, we find that the PLL embedded equalizer for DSCM scheme exhibits better skew and phase noise compensation ability compared with other equalizers. Additionally, we compare the performance of single-carrier and DSCM schemes with the proposed equalizers in simulation. The influence of phase noise and Tx IQ skew on DSCM transmission can be largely relaxed.

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

1. Introduction

With the rapid development of the Internet, mobile broadband and data centers, the requirement of network transmission bandwidth has been explosive growth, promoting the development of ultra-large capacity optical transmission systems. Beyond 100 Tb/s transmissions in the C + L band of single-mode fiber have been demonstrated [13], where the transmission symbol rate is expected to be 100 Gbaud per channel. The single-carrier transmission scheme has achieved huge commercial success these years, the digital signal processing (DSP) algorithms are mature and have many advantages in signal recovery, e.g. dispersion compensation, carrier recovery and fast-tracking of the state of polarization (SOP). However, as the application of ultra-large bandwidth and higher spectral-efficiency modulated signal, the communication scheme and various linear impairments have to be reexamined. The complexity of dispersion compensation is greatly increased due to the ultra-high baud rate [4]. Besides, the equalization-enhanced phase noise (EEPN) impairment imposes a tighter constraint on the receiver laser phase noise for coherent systems with higher modulation format and large electronic-compensated dispersion [5]. In wideband single-channel transmission, the polarization mode dispersion (PMD) parameter is dependent on frequency, which implies that we have to cope with increasingly complicated PMD impairments in wideband channels [6]. In addition, the self-phase modulation penalty increases with the baud rate [7]. Therefore, DSCM has been attracting more attention to enable flexible data rate adaptivity and provide robustness against various linear impairments, such as dispersion compensation, EEPN and frequency-dependent PMD effect [8]. The cumulated dispersion compensation increases with the square of the baud rate and the variance of EEPN is inversely proportional to the symbol rate. This simplifies the chromatic dispersion compensation (CDC) and improves the EEPN tolerance. In the narrowband sub-channel transmission, the PMDs at different frequencies are strongly correlated, therefore, the traditional 2 × 2 MIMO equalizer is still feasible for polarization demultiplexing.

The DSCM is expected in the evolution of ultra-wideband signal transmission. However, when there exists IQ mismatch and skew, subcarriers suffer interference from their “symmetric” subcarriers [9]. Besides, the IQ skew introduces frequency-dependent interference which is larger at higher frequency bands [10,11]. Currently, numerous analog or digital domain calibration methods for the IQ skew estimation of coherent optical transceivers have been proposed for single-carrier transmissions [1215]. However, the residual IQ skew after calibration or from some dynamic effects still needs to be mitigated, especially for the application of higher-level modulation format signals. Several adaptive equalization schemes have demonstrated their effectiveness to mitigate the dynamic receiver (Rx) IQ skew [1622]. But in fact, the skew of the transmitter is more difficult to deal with since it impairs the whole DSP chain compared with Rx-side impairments [23]. Some algorithms for transmitter (Tx) IQ skew equalization have been proposed [2426], where the 4 × 2 widely-linear (WL) MIMO equalizer has demonstrated its effectiveness to compensate the Tx IQ skew. However, these schemes are ineffective in DSCM transmission because of the conjugate interference resulting from the “symmetric” subcarrier. The 4 × 4 WL equalizer is useful to compensate for the Rx skew in the DSCM system, however, its performance is affected by the phase noise for Tx IQ skew compensation [27].

Since DSCM transmissions have longer symbol periods per subcarrier when compared to the equivalent-rate single-carrier transmissions, there is a performance degradation if a conventional per-subcarrier laser phase noise compensation method is employed [28,29]. In general, the carrier phase estimation (CPE) can be operated in data-aided or blind manner. By inserting pilot symbols, the phase noise can be directly estimated from the phase difference between the transmitted and received symbols. The extension of pilot-based joint CPE has already been proposed in [30], but the additional overhead reduces the net transmission rate. Blind CPE approaches include techniques such as the Viterbi-Viterbi (VV) algorithm and BPS, which provide an overhead-free scheme at the cost of increased complexity and modulation format dependence [29,31]. In addition, the training symbols used for de-polarization can similarly be used for digital PLL to achieve phase noise compensation. Since the accuracy of phase noise estimation affects the performance of Tx IQ skew compensation, we need to use a better phase noise estimation scheme that can mitigate the effect of Tx IQ skew. When the PLL is used, the carrier phase is updated by using the complex signals after MIMO equalizer and Tx impairment compensation filter. So carrier phase recovery is not contaminated by Tx impairments [26]. However, the phase offset increases at higher frequency bands due to the Tx IQ skew. It is difficult for PLL to obtain accurate phase tracking. Moreover, the PLL implemented on each subcarrier for the DSCM transmission increases the overall complexity.

In this paper, we propose a paired-subcarrier equalization scheme that uses the PLL embedded 4 × 4 equalizer to jointly compensate for laser phase noise and Tx IQ skew in DSCM transmission systems. We use training symbols from each pair of subcarriers as the input of 4 × 4 equalizer since the Tx IQ skew for each subcarrier suffers interference from the conjugate components of the “symmetric” subcarrier. To alleviate the complexity, we use a shared phase noise estimation scheme. We first estimate the phase noise by the inner paired subcarriers to reduce the frequency-dependent influence from the Tx IQ skew. This helps us have a more accurate phase noise estimation. Then, we share the estimated phase noise with other paired subcarriers. We conduct numerical simulations of 100-GBaud 64-QAM DSCM coherent optical fiber transmission that consists of eight 12.5-Gbaud subcarriers and experiment of 10-GBaud four-subcarriers PM-16QAM transmission. The simulation results show that the PLL embedded 4 × 4 equalizer can simultaneously compensate for the laser phase noise and Tx IQ skew. The scheme of shared phase noise estimation exhibits better skew tolerance and can reduce the complexity. Additionally, compared with the single-carrier system, it shows the same performance under varied optical signal-to-noise ratio (OSNR) and better performance on the tolerance of EEPN. The proposed equalization method can relax the phase noise and Tx IQ skew limitations for DSCM transmission schemes in ultra-large capacity transmission applications.

2. Principle of the proposed equalization

The IQ skew results from the misalignment of in-phase- and quadrature-component, e.g. the mismatch of the high-speed cable, group delay of DAC/ADC, drivers, etc. We use $x(t )\textrm{ } = \textrm{ }I(t )\textrm{ } + \textrm{ }jQ(t )$ to express the transmitted signal, so the distorted signal due to IQ skew can be expressed as, $y(t )$

$$\left[ {\begin{array}{{c}} {{I_t}(t )}\\ {{Q_t}(t )} \end{array}} \right] = \left[ {\begin{array}{{cc}} 1&0\\ 0&{\delta ({t - \tau } )} \end{array}} \right] \otimes \left[ {\begin{array}{{c}} {I(t )}\\ {Q(t )} \end{array}} \right]$$
where $y(t )= {I_t}(t )+ j{Q_t}(t )$ is the received signal containing the IQ skew and τ is the skew. The received signal can be rewritten in terms of the transmitted signal and its conjugate as,
$$y(t )= \frac{{x(t )+ {x^ \ast }(t )}}{2} + j\left( {\delta ({t - \tau } )\otimes \left( { - \frac{{j({x(t )- {x^ \ast }(t )} )}}{2}} \right)} \right) = {k_1}(t )\otimes x(t )+ {k_2}(t )\otimes {x^ \ast }(t )$$
where ${k_1} = \frac{{1 + \delta ({t - \tau } )}}{2},{k_2} = \frac{{1 - \delta ({t - \tau } )}}{2}$, “*” and ${\otimes}$ mean conjugation and convolution respectively. Because the 2 × 2 equalizer structure enforces joint filtering of the in-phase and quadrature signal components, four complex-valued independent filters are unable to compensate for any imbalance between the in-phase and quadrature components [32]. Therefore, the augmented WL complex signal processing is required in the receiver. Conjugation which is applied on the time domain signal corresponds to a conjugation and spectral sign flip in the frequency domain. So the expression in the frequency domain of received signal by Fourier transform of Eq. (2) can be written as,
$$Y(w )= \frac{1}{2}[{X(w )({1 + {e^{ - jw\tau }}} )+ {X^ \ast }({ - w} )({1 - {e^{ - jw\tau }}} )} ]$$
Where w is the signal’s angular frequency, $X(w )$ and $Y(w )$ are the Fourier transform of $x(t )$ and $y(t )$, respectively. Assuming that $x(t )$ is a DSCM signal, which is expressed as,
$$x(t )= \sum\limits_{n = 1}^N {{x_n}(t ){e^{j2\pi {f_n}t}}}$$
where ${x_n}(t )$ and ${f_n}$ are the baseband signal and center frequency for each subcarrier, respectively. The N is the number of subcarriers. The Eq. (3) can be rewritten as,
$$\begin{array}{l} Y(w )= \frac{1}{2}\sum\limits_{n = 1}^N {({{X_n}({w - 2\pi {f_n}} )({1 + {e^{ - jw\tau }}} )+ {X_n}^ \ast ({ - w - 2\pi {f_n}} )({1 - {e^{ - jw\tau }}} )} )} \\ \textrm{ } = \frac{1}{2}\sum\limits_{n = 1}^N {({{X_n}({w - 2\pi {f_n}} )({1 + \cos w\tau - j\sin w\tau } )+ {X_n}^ \ast ({ - w - 2\pi {f_n}} )({1 - \cos w\tau + j\sin w\tau } )} )} \end{array}$$

In Eq. (5), each subcarrier experiences frequency up-conversion. The received signal after down-conversion consists of the useful signal and the interference from its “symmetric” subcarrier. These paired interfering subcarriers follow the relationship $n + m = N + 1$, where m is the index of “symmetric” subcarrier-n. There is no interference from other subcarriers except for the conjugate subcarrier.

Figure 1 illustrates the Tx IQ skew effect of a DSCM signal with eight subcarriers in the frequency domain. The amplitude of the spectrum is flat because the spectrum is the superposition of subcarriers and their conjugate interference. It is worth noting that the penalty caused by a certain skew varies periodically with the frequency of subcarriers [9,10]. This can be explained by the varied with $w\tau$ in Eq. (5). As the $w\tau$ increases, $1 - \cos w\tau$ is increased for $0 \le w\tau \le \pi$ and decreased for $\pi \le w\tau \le 2\pi$. Since the alias of a pair of symmetric subcarriers exists in the spectrum, the 2 × 2 real-valued equalizer used in single-carrier system after carrier recovery is unable to compensate for the Tx IQ skew. Thus, we have to use the 2 × 2 WL complex-value equalizer to compensate for the Tx IQ skew after depolarization or the 4 × 4 WL complex-valued equalizer to achieve joint depolarization and Tx IQ skew compensation. To facilitate explanation, we focus on two generic frequency-symmetric subcarriers in the following discussions. We use sub-n and sub-m, chosen such that ${f_n} ={-} {f_m}$.

 figure: Fig. 1.

Fig. 1. Illustration of Tx IQ skew effect in a DSCM signal with eight subcarriers.

Download Full Size | PDF

Figure 2 shows four equalization configurations. Configuration I is a conventional de-polarization and phase noise compensation scheme [8]. The de-polarization is realized through a 2 × 2 strictly linear (SL) complex-valued equalizer for each subcarrier. Then, pilot symbols are used to estimate the laser phase noise. Configuration II is also a conventional de-polarization and phase noise compensation as Configuration I, and 2 × 2 WL complex-valued equalizer which has the same performance with 4 × 4 real-valued equalizer is used to compensate the Tx IQ skew [9]. Moreover, two 2 × 2 WL complex-valued equalizers can be substituted by a 4 × 4 WL equalizer and the phase noise is compensated independently [33]. They have similar performance, we attribute these two algorithms as one kind. Configuration III uses the PLL to track the phase noise and realizes the joint equalization to mitigate the Tx IQ skew [26]. Configuration IV is the proposed PLL embedded 4 × 4 WL equalizer for each pair of subcarriers. Since the existence of Tx IQ skew causes penalties for de-polarization and phase noise compensation, errors are accumulated in the following compensation. So we propose a PLL embedded 4 × 4 WL equalizer to jointly perform de-polarization, phase noise and Tx IQ skew compensation to reduce the error accumulation problem. The finite impulse response (FIR) structure of the equalizer is,

$$\begin{array}{c} \left[ {\begin{array}{{c}} {{X_n}}\\ {{Y_n}} \end{array}} \right] = \left[ {\begin{array}{{cc}} {{h_{{x_n}{x_n}}}}&{{h_{{x_n}{y_n}}}}\\ {{h_{{y_n}{x_n}}}}&{{h_{{y_n}{y_n}}}} \end{array}} \right]\left[ {\begin{array}{{c}} {{x_n}{e^{ - i{\theta_n}}}}\\ {{y_n}{e^{ - i{\theta_n}}}} \end{array}} \right] + \left[ {\begin{array}{{cc}} {{h_{{x_n}x_m^\ast }}}&{{h_{{x_n}y_m^\ast }}}\\ {{h_{{y_n}x_m^\ast }}}&{{h_{{y_n}y_m^\ast }}} \end{array}} \right]\left[ {\begin{array}{{c}} {x_m^\ast {e^{i{\theta_m}}}}\\ {y_m^\ast {e^{i{\theta_m}}}} \end{array}} \right]\\ \left[ {\begin{array}{{c}} {{X_m}}\\ {{Y_m}} \end{array}} \right] = \left[ {\begin{array}{{cc}} {{h_{{x_m}{x_m}}}}&{{h_{{x_m}{y_m}}}}\\ {{h_{{y_m}{x_m}}}}&{{h_{{y_m}{y_m}}}} \end{array}} \right]\left[ {\begin{array}{{c}} {{x_m}{e^{ - i{\theta_m}}}}\\ {{y_m}{e^{ - i{\theta_m}}}} \end{array}} \right] + \left[ {\begin{array}{{cc}} {{h_{{x_m}x_n^\ast }}}&{{h_{{x_m}y_n^\ast }}}\\ {{h_{{y_m}x_n^\ast }}}&{{h_{{y_m}y_n^\ast }}} \end{array}} \right]\left[ {\begin{array}{{c}} {x_n^\ast {e^{i{\theta_n}}}}\\ {y_n^\ast {e^{i{\theta_n}}}} \end{array}} \right] \end{array}$$
where ${\theta _{^n}}$ and ${\theta _{^m}}$ are the estimated phase noise of the n-th and m-th subcarrier, respectively. For simple description, the m-th subcarrier is taken as an example, and the n-th subcarrier only needs to convert m and n in the following Eq. (7)–(9) according to the relationship $n + m = N + 1$. The error signals ${e_{{x_m}}}$ and ${e_{{y_m}}}$ using the least mean square (LMS) algorithm are expressed as,
$${e_{{x_m}}}(k )= T{x_m}(k )- {X_m}(k ),{e_{{y_m}}}(k )= T{y_m}(k )- {Y_m}(k )$$
where $T{x_m}$ and $T{y_m}$ are the training symbols of the m-th subcarrier. Then we can define the cost function using the LMS algorithm as,
$${J_{{x_m}}}(h )= {({{e_{{x_m}}}(k )} )^2},{J_{{y_m}}}(h )= {({{e_{{y_m}}}(k )} )^2}$$

 figure: Fig. 2.

Fig. 2. The structure of four equalization configurations. Configuration I: 2 × 2 equalizer for de-polarization, phase noise compensation based pilot symbols; Configuration II: 2 × 2 equalizer for de-polarization, phase noise compensation based pilot symbols and Tx IQ skew compensation for each subcarrier by 2 × 2 equalizer independent; Configuration III: 2 × 2 equalizer for de-polarization and Tx skew compensation for each subcarrier by 2 × 2 equalizer combined with PLL; Configuration IV: PLL embedded 4 × 4 equalizer for joint equalization of polarization, phase noise and Tx IQ skew.

Download Full Size | PDF

Correspondingly, the updated equations of the filter coefficients are obtained by the gradient descent method to minimize the cost function,

$$\begin{array}{l} {h_{{x_m}{x_m}}}({k + 1} )= {h_{{x_m}{x_m}}}(k )+ \mu {e_{{x_m}}}x_m^\ast (k ){e^{i{\theta _m}}},{h_{{x_m}{y_m}}}({k + 1} )= {h_{{x_m}{y_m}}}(k )+ \mu {e_{{x_m}}}y_m^\ast (k ){e^{i{\theta _m}}}\\ {h_{{x_m}x_n^\ast }}({k + 1} )= {h_{{x_m}x_n^\ast }}(k )+ \mu {e_{{x_m}}}{x_n}(k ){e^{ - i{\theta _n}}},{h_{{x_m}y_n^\ast }}({k + 1} )= {h_{{x_m}y_n^\ast }}(k )+ \mu {e_{{x_m}}}{y_n}(k ){e^{ - i{\theta _n}}}\\ {h_{{y_m}{x_m}}}({k + 1} )= {h_{{y_m}{x_m}}}(k )+ \mu {e_{{y_m}}}x_m^\ast (k ){e^{i{\theta _m}}},{h_{{y_m}{y_m}}}({k + 1} )= {h_{{y_m}{y_m}}}(k )+ \mu {e_{{y_m}}}y_m^\ast (k ){e^{i{\theta _m}}}\\ {h_{{y_m}x_n^\ast }}({k + 1} )= {h_{{y_m}x_n^\ast }}(k )+ \mu {e_{{y_m}}}{x_n}(k ){e^{ - i{\theta _n}}},{h_{{y_m}y_n^\ast }}({k + 1} )= {h_{{y_m}y_n^\ast }}(k )+ \mu {e_{{y_m}}}{y_n}(k ){e^{ - i{\theta _n}}} \end{array}$$

Since the Tx IQ skew leads to the phase offset increasing at higher frequencies, it is difficult for PLL to realize accurate phase tracking. Besides, the phase noise of subcarriers is approximately equal when the transmission distance is short. So, we use a shared architecture to realize the phase noise compensation, which can improve the anti-skew ability and reduce the complexity of DSP. First, we estimate the phase noise of the inner subcarriers using PLL, which is less affected by the skew, so the estimated phase noise is more accurate. Then, this phase noise is directly shared with other pairs of subcarriers. The walk-off effect is induced by the fact that different frequency components of the optical pulse have different group velocities, which leads to the time delay between subcarriers. So the phase noise for different subcarriers has a relative delay [34]. For longer transmission distances, we can increase the number of reference subcarriers for phase noise estimation to alleviate the influence of the walk-off effect.

3. Simulation setup and discussions

As shown in Fig. 3, we carry out numerical simulations by VPItransmissionMakerTM9.1 of 100-Gbaud PDM-64-QAM digital subcarrier multiplexing transmission consisting of eight subcarriers to investigate the performance of our proposed joint equalizer. At the Tx, the roll-off factor of the RRC spectral shaping is 0.1. And we also digitally emulate the Tx IQ skew. Then the transmitted signals are fed into a dual-polarization IQ modulator (DP-IQM). An external cavity laser (ECL) with a central wavelength of 1550 nm and a linewidth of 100 kHz is used as the optical carrier. After electrical-to-optical (E/O) conversion, an erbium-doped fiber amplifier (EDFA) is used for launched optical power (LOP) adjusting. The signal is transmitted over 100km-SSMF with an attenuation of 0.2 dB/km and a dispersion coefficient of 16.0 ps/nm/km. The loss of the SSMF is fully compensated by the EDFA. An optical AWGN source loads noise to the optical signal transmitted through the fiber in order to adjust optical signal’s OSNR. At the receiver, the linewidth and optical power of the local oscillator laser are 100 KHz and 10 dBm, respectively. The frequency offset (FO) is set to 200 MHz. Then, a coherent receiver is utilized for optical-to-electrical (O/E) conversion. Figure 3(b) shows the received spectrum with eight sub-bands. Figure 3 (c) depicts the DSP flow. After the FO compensation, the subcarriers are demultiplexed through down-conversion and matched filtering. Next, we compensate the CD for each subcarrier. Then, the proposed pair-subcarrier joint equalizer is used to realize the depolarization, carrier phase recovery and Tx skew compensation. Finally, the BER is measured by the direct counting method. The Q-factor is calculated from the measured BER $Q[{dB} ]= 20{\log _{10}}\left( {\sqrt 2 erf{c^{ - 1}}({2 \times BER} )} \right)$.

 figure: Fig. 3.

Fig. 3. (a)Simulation setup. PBS: polarization beam splitter, PBC: polarization beam coupler, SSMF: standard single mode fiber, EDFA: Erbium-doped optical fiber amplifier. (b)The spectrum of DSCM after the effect of Tx IQ skew (c) Tx DSP and Rx DSP.

Download Full Size | PDF

3.1 Performance of the proposed algorithm in DSCM transmission

We first investigate the performance of the proposed equalizer in the DSCM system under optical back-to-back without phase noise. The equalization configurations are similar to that in Section 2, Configuration I only uses one 2 × 2 SL for each subcarrier. Configurations II and III have the same structure that one 2 × 2 SL equalizer is used for each subcarrier and two 2 × 2 WL equalizers are used for each pair of subcarriers for Tx IQ skew compensation since the phase noise is not considered. So we use Configuration II&III in Fig. 4(a) to represent it. Configuration IV uses the 4 × 4 WL equalizer for each pair of subcarriers. Figure 4(a) shows the Q-factor versus Tx IQ skew using different equalization configurations without phase noise. In Fig. 4(a), the performance of the four configurations is the same when the Tx IQ skew is 0ps because the depolarization in each configuration is not disturbed by skew and phase noise. We find that the penalty without equalizer is up to 9.4 dB within the case of 3-ps Tx IQ skew and only the Configuration IV can fully compensate the Tx IQ skew. Compared with the case without equalizer, Configuration I only reduces 1.4 dB penalty when the Tx IQ skew is 3 ps. This is because the 2 × 2 SL complex-valued equalizer cannot deal with the Tx IQ skew and the small gain comes from equalizing the divergent constellations caused by the Tx IQ skew. In addition, the Q-penalty is only 0.26 dB for Configuration II&III. It effectively compensates for most of the IQ-skew impairment, and the residual penalty comes from the insufficient equalization accuracy caused by the mutual influence of Tx IQ skew and polarization. Overall, the 4 × 4 WL equalizer is the best scheme to deal with the Tx IQ skew.

 figure: Fig. 4.

Fig. 4. Q-factor Vs Tx IQ skew for different equalization configurations in DSCM system (a) optical back-to-back without phase noise (b) transmission over 100 km SSMF with phase noise.

Download Full Size | PDF

Figure 4(b) depicts the Q-factor versus Tx IQ skew after 100-km fiber transmission with phase noise. We compare the four equalization configurations defined in Fig. 2. The pilot insertion rate is 1/32 for all configurations. In Fig. 4(b), the difference of Q-factors at Tx IQ skew of 0ps is dependent on the phase noise tolerance of different configurations. We find that Configuration IV can almost completely compensate for Tx IQ skew while Configurations II and III have 1.50-dB and 1.15-dB Q-penalty at Tx IQ skew of 3ps, respectively. The penalty is due to the presence of Tx IQ skew which results in the inaccuracy of fully compensating for the polarization effect and phase noise. Therefore, the PLL embedded 4 × 4 WL equalizer can simultaneously achieve polarization demultiplexing, phase noise and Tx skew compensation that is more suitable for the large capacity and high-level modulated transmission systems.

Figure 5(a) shows the Q-factor versus Tx IQ skew with the two schemes of phase noise estimation. 1) independent mode: each subcarrier uses PLL to track the phase noise, 2) shared mode: only a pair of subcarriers use the PLL to track the phase noise, and other subcarriers use the estimated phase noise directly. The difference is that the phase noise is estimated in independent or shared mode, while they both use the same 4 × 4 equalizer architecture. The “Sub. m&n” or “Sub. m&n in shared mode” means the performance of each paired subcarriers in independent or shared mode, where the shared phase noise is estimated by Sub. 4&5. The Q-factor of Fig. 5 (a) is calculated from the average BER of both polarization states for each paired subcarrier. The Q-factor in all other figures is calculated from the average BER of both polarization states for all subcarriers. In Fig. 5(a), the inner subcarriers show better performance while the outer subcarriers sub-1&8 have over 0.23-dB Q-penalty than the inner subcarriers sub-4&5 with 3-ps Tx skew using independent phase noise compensation. This result is also consistent with Fig. 1, where the impairment of outer subcarrier is aggravated by Tx IQ skew, leading to a decrease in the accuracy of phase noise estimation. Therefore, when we share the phase noise estimation results from the inner two subcarriers with the minimum interference to the rest of the subcarriers, we can obtain overall better skew tolerance.

 figure: Fig. 5.

Fig. 5. (a) Q-factor Vs Tx IQ skew of PLL embedded 4 × 4 equalizer in independent and shared mode (b) Q-factor versus transmission distance of the overall performance of DSCM system with the PLL embedded 4 × 4 equalizer in independent mode and shared mode with different reference subcarriers.

Download Full Size | PDF

Figure 5(b) shows the Q-factor versus transmission distance of the overall performance of DSCM system with the PLL embedded 4 × 4 equalizer in independent mode and shared mode with different reference subcarriers. The Tx IQ skew is 3 ps. In Fig. 5 (b), we find that the performance of Configuration IV in shared mode with Sub.4&5 as reference subcarriers degrades faster compared with Configuration IV in independent mode. Since there is a relative delay in the received phase of each subcarrier because of the existence of the walk-off effect, the phase offset is larger for outer subcarriers when we choose Sub.4&5 as the reference subcarriers to estimate the phase noise. In this case, we should use the other inner subcarrier pairs as reference subcarriers to estimate the phase noise or increase the number of reference subcarriers for phase noise estimation to alleviate the influence of the walk-off effect. To maintain a low-complexity compensation, we use Sub.2&7 as reference subcarriers for comparison. As the dash-dotted lines shown in Fig. 5(b), the overall performance with Sub.2&7 as reference subcarriers shows about 0.4-dB performance gain compared with Sub.4&5 as reference subcarriers in 3000-km transmission. Besides, the shared mode with both Sub.4&5 and Sub.2&7 as reference subcarriers shows better performance compared with only one paired reference subcarrier. The shared mode with both Sub. 4&5 and Sub. 2&7 is achieved through sharing the phase estimated from Sub. 4&5 to Sub. 3&6 and Sub. 2&7 to Sub. 1&8, respectively. Moreover, the shared mode with Sub.4&5 as reference subcarriers has a better performance compared with the independent mode when the transmission distance is within 150 km. So, the phase noise compensation in shared mode is more suitable for short-reach applications. For long-haul transmissions, we can use independent mode or more reference subcarriers.

3.2 Performance comparison of single-carrier and DSCM transmission

We also compare the transmission performance of DSCM system using the proposed joint equalizer and single-carrier system using the PLL embedded 4 × 2 WL equalizer [26]. The symbol rate and modulation format of the singe-carrier system are also 100GBaud and 64-QAM, respectively. Moreover, we also verify the DSCM system performance with independent and shared phase noise estimation modes, respectively. Figure 6 shows the Q-factor versus OSNR with or without 3-ps Tx IQ skew after transmission over 100-km SSMF. The performance of DSCM system with independent phase estimation or shared phase estimation mode embedded 4 × 4 equalizer is almost the same to that of the single-carrier system.

 figure: Fig. 6.

Fig. 6. Q-factor Vs OSNR for comparison of single-carrier and DSCM systems with/without 3ps Tx skew after transmission over 100 km SSMF.

Download Full Size | PDF

Then, we compare the influence of EEPN on the two transmission systems. Figure 7 depicts the Q-factor of single-carrier and DSCM transmission systems with and without 3-ps Tx skew as a function of transmission distance. From Fig. 7, we can see that the performance of single-carrier system decreases sharply while the DSCM system has a penalty within 0.2 dB as transmission distance increases. This indicates that the DSCM system has a better tolerance to EEPN than the single-carrier system.

 figure: Fig. 7.

Fig. 7. Q-factor Vs transmission distance of single-carrier and DSCM systems with/without 3ps Tx IQ skew.

Download Full Size | PDF

4. Experiment setup and results

We experimentally evaluate the effectiveness of the proposed method by testing a 40-GBaud PM-16QAM transmission consisting of four subcarriers. This mismatch of data rate with previous simulations is resulted from the limited band-width of transceivers. The 40-GBaud PM-16QAM transmission can also be used to verify the performance of the proposed algorithm. Figure 8 shows the experimental setup and offline DSP. On the Tx side, the 215 16-QAM symbols are pulse shaped by a root-raised cosine (RRC) filter with a roll-off factor of 0.1 for each subcarrier. The Tx IQ skew is digitally emulated in Tx DSP. Since the Tx IQ skew of the two polarizations is basically different in the real system, the Tx IQ skew of digital emulation in X-polarization is 2ps larger than that in Y polarization. They are uploaded to an 8-bit, four-channel, 120-GSa/s arbitrary waveform generator (AWG) to generate the driving signals. A tunable 100-KHz linewidth (ECL) operating at 1550 nm generates the optical signal at a power of 11 dBm. Then, a commercially available LiNbO3 dual-polarization IQ modulator with a 3-dB bandwidth of 23 GHz is used for electro-optic conversion.

 figure: Fig. 8.

Fig. 8. Experimental setup and offline DSP.

Download Full Size | PDF

During the transmission, an erbium-doped fiber amplifier (EDFA) with a noise figure of 5 dB is used to set the launch power. The signal is transmitted over 100-km-SSMF with an attenuation of 0.2 dB/km and a dispersion coefficient of 16.8 ps/nm/km. After compensating for the fiber loss through the post-EDFA, a variable optical attenuator (VOA) is used to control the received optical power (ROP). The LO laser is an ECL with 100-KHz linewidth operating at 1549.947 nm with a power of 9.5 dBm. After coherent detection and electrical amplification, the signal is digitized by a digital storage oscilloscope (DSO) operating at 256 GSa/s for offline DSP.

In the offline DSP, the Rx IQ skew is first compensated by a fixed Rx frontend filter [26]. Then the signal is resampled to 2 Sa/sym. The remaining DSP processes include frame synchronization, FO compensation, subcarrier demultiplexing and matched filter. The FO estimation is based on the fourth power of the inner subcarrier. Then we compensate for the FO for the DSCM signal in the frequency domain. After the CDC for each subcarrier, we compare the four equalizers as shown in Fig. 2. After demodulation, the BER is measured by directly counting and the Q-factor is calculated from the BER.

Figure 9 depicts the Q-factor versus Tx Y-IQ skew after 100-km fiber transmission for four equalizer configurations. Note that the Q-factor is averaged over two polarization signals. In Fig. 9, Configuration IV, i.e., the PLL embedded 4 × 4 WL equalizer, can almost completely compensate for the Tx skew and outperform other equalizers. Configuration I is used to measure the impact of Tx IQ skew, which linearly degrades the Q-factor. There is a large gap between Configuration I to others when the Tx Y-IQ skew is 0ps. Because the X-polarization is affected with a 2-ps Tx X-IQ skew. Compared with configuration IV, configurations II and III have 0.60-dB and 1.54-dB Q-penalty at Tx Y-IQ skew of 10ps, respectively. The trend of experiment results agrees well with simulation results. The proposed joint equalizer can simultaneously achieve polarization demultiplexing, phase noise and Tx IQ skew compensation that is suitable for large-capacity transmission systems.

 figure: Fig. 9.

Fig. 9. Q-factor Vs Tx Y-IQ skew for different algorithms in DSCM system transmission over 100 km SSMF.

Download Full Size | PDF

5. Conclusion

In this paper, we have proposed a PLL-embedded paired-subcarrier equalization scheme for the joint equalization of polarization, phase noise and Tx IQ skew for DSCM transmission systems. The proposed PLL embedded 4 × 4 equalizer can deal with the paired-subcarrier interference in the DSCM transmission, so it has superior performance. Moreover, the comparison between the single-carrier and the DSCM system using the proposed equalizer illustrates that the influence of phase noise and Tx IQ skew on DSCM transmission can be largely relaxed.

Funding

National Natural Science Foundation of China (62111530150, U22A2086); Science and Technology Commission of Shanghai Municipality (SKLSFO2021-01); Fundamental Research Funds for the Central Universities (ZYGX2019J008, ZYGX2020ZB043).

Acknowledgments

This work was supported by National Science Foundation of China (NSFC) (U22A2086 and 62111530150), STCSM (SKLSFO2021-01), Fundamental Research Funds for the Central Universities (ZYGX2020ZB043 and ZYGX2019J008).

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. D. Qian, M.-F. Huang, E. Ip, Y.-K. Huang, Y. Shao, J. Hu, and T. Wang, “High Capacity/Spectral Efficiency 101.7-Tb/s WDM Transmission Using PDM-128QAM-OFDM Over 165-km SSMF Within C- and L-Bands,” J. Lightwave Technol. 30(10), 1540–1548 (2012). [CrossRef]  

2. A. Sano, T. Kobayashi, S. Yamanaka, A. Matsuura, H. Kawakami, Y. Miyamoto, K. Ishihara, and H. Masuda, “102.3-Tb/s (224 × 548-Gb/s) C- and Extended L-band All-Raman Transmission over 240 km Using PDM-64QAM Single Carrier FDM with Digital Pilot Tone,” in Optical Fiber Communication Conference, OSA Technical Digest (Optica Publishing Group, 2012), paper PDP5C.3.

3. M. Luo, C. Li, Q. Yang, Z. He, J. Xu, Z. Zhang, and S. Yu, “100.3-Tb/s(375×267.27-Gb/s) C- and L-band Transmission over 80-km SSMF Using DFT-S OFDM 128-QAM,” in Asia Communications and Photonics Conference 2014, OSA Technical Digest (online) (Optica Publishing Group, 2014), paper AF4B.1.

4. B. Spinnler, “Equalizer design and complexity for digital coherent receivers,” IEEE J. Select. Topics Quantum Electron. 16(5), 1180–1192 (2010). [CrossRef]  

5. W. Shieh and K. Ho, “Equalization-enhanced phase noise for coherent-detection systems using electronic digital signal processing,” Opt. Express 16(20), 15718–15727 (2008). [CrossRef]  

6. N. Cui, X. Zhang, Q. Zhang, and X. Zhang, “Narrow- or wide-band channel for a high baud rate fiber communication system: a judgment based on a temporal and spectral evolution PMD model,” Opt. Express 29(23), 38497–38511 (2021). [CrossRef]  

7. “The Ultimate Guide to Nyquist Subcarriers,” 2019. [Online] Available: https://www.infinera.com/wp-content/uploads/The-Ultimate-Guide-to-Nyquist-Subcarriers-0208-WP-RevA-0719.pdf.

8. H. Sun, M. Torbatian, M. Karimi, et al., “800 G DSP ASIC Design Using Probabilistic Shaping and Digital Sub-Carrier Multiplexing,” J. Lightwave Technol. 38(17), 4744–4756 (2020). [CrossRef]  

9. Z. Zhai, M. Fu, L. Liu, H. Jiang, H. Ren, L. Yi, W. Hu, and Q. Zhuge, “Transmitter IQ Mismatch Compensation and Monitoring for Digital Subcarrier-multiplexing Systems,” in Asia Communications and Photonics Conference/International Conference on Information Photonics and Optical Communications 2020 (ACP/IPOC), OSA Technical Digest (Optica Publishing Group, 2020), paper M4A.319.

10. G. Bosco, S. M. Bilal, A. Nespola, P. Poggiolini, and F. Forghieri, “Impact of the Transmitter IQ-Skew in Multi-Subcarrier Coherent Optical Systems,” in Optical Fiber Communication Conference, OSA Technical Digest (online) (Optica Publishing Group, 2016), paper W4A.5.

11. B. Baeuerle, A. Josten, R. Bonjour, D. Hillerkuss, and J. Leuthold, “Effect of Transmitter Impairments on Nyquist-FDM Signals with Increasing Sub-band Granularity,” in Advanced Photonics 2016 (IPR, NOMA, Sensors, Networks, SPPCom, SOF), OSA Technical Digest (online) (Optica Publishing Group, 2016), paper SpW3F.4.

12. Y. Fan, X. Su, H. Chen, J. Liang, Z. Tao, H. Nakashima, and T. Hoshlda, “Experimental verification of IQ imbalance monitor for high order modulated transceivers,” in Proc. Eur. Conf. Opt. Commun., Rome, Italy, 2018, pp. 23–27.

13. N. Liu, C. Ju, C. Li, Z. Chi, Z. Zhang, and X. Chen, “In-service blind transceiver calibration for long-haul non-dispersion managed systems,” Opt. Quantum Electron. 52(1), 37 (2020). [CrossRef]  

14. Y. Fan, Y. Jiang, J. Liang, Z. Tao, H. Nakashima, and T. Hoshida, “Transceiver IQ imperfection monitor by digital signal processing in coherent receiver,” in Proc. Optoelectron. Commun. Conf. Int. Conf. Photon. Switching Comput., Fukuoka, Japan, 2019, pp. 7–11.

15. D. L. Longquan Dai, C. Huang, H. Li, Y. Yang, H. Song, M. Cheng, Q. Yang, M. Tang, D. Liu, and L. Deng, “Experimental Demonstration of Simultaneously Precise Tx and Rx Skew Calibration for Coherent Optical Transceiver,” J. Lightwave Technol. 40(4), 1043–1054 (2022). [CrossRef]  

16. M. Paskov, D. Lavery, and S. J. Savory, “Blind Equalization of Receiver In-Phase/Quadrature Skew in the Presence of Nyquist Filtering,” IEEE Photonics Technol. Lett. 25(24), 2446–2449 (2013). [CrossRef]  

17. M. S. Faruk and K. Kikuchi, “Compensation for In-Phase/Quadrature Imbalance in Coherent-Receiver Front End for Optical Quadrature Amplitude Modulation,” IEEE Photonics J. 5(2), 7800110 (2013). [CrossRef]  

18. E. P. Silva and D. Zibar, “Widely Linear Equalization for IQ Imbalance and Skew Compensation in Optical Coherent Receivers,” J. Lightwave Technol. 34(15), 3577–3586 (2016). [CrossRef]  

19. R. Rios-Müller, J. Renaudier, and G. Charlet, “Blind Receiver Skew Compensation and Estimation for Long-Haul Non-Dispersion Managed Systems Using Adaptive Equalizer,” J. Lightwave Technol. 33(7), 1315–1318 (2015). [CrossRef]  

20. J. Liang and W. Wang, “Receiver IQ Imbalance and Skew Compensation By Frequency Domain 4 × 2MIMO Equalizer,” in Conference on Lasers and Electro-Optics, OSA Technical Digest (Optica Publishing Group, 2021), paper JTu3A.129.

21. . J. Liang, W. Wang, M. Xiang, and Y. Cai, “Receiver IQ Imbalance and Skew Compensation By Frequency Domain Widely Linear Equalizer,” in Asia Communications and Photonics Conference/International Conference on Information Photonics and Optical Communications 2020 (ACP/IPOC), OSA Technical Digest (Optica Publishing Group, 2020), paper T4B.4.

22. J. Liang, Z. Feng, Z. Wang, L. Zhang, Y. Wang, Y. Tu, and W. Wang, “Receiver IQ Imbalance And Skew Compensation for High order Modulation Formats By Frequency Domain 4×2 MIMO,” in Asia Communications and Photonics Conference 2021, Technical Digest Series (Optica Publishing Group, 2021), paper T3G.5.

23. G. Di Rosa, R. Emmerich, M. Sena, J. K. Fischer, C. Schubert, and A. Richter, “Impact of wavelength-dependent I/Q imbalances of standard c-band transceivers in rate-adaptive multiband systems,” in Proc. Eur. Conf. Opt. Commun, 2021, pp. 1–4.

24. C. Fludger and T. Kupfer, “Transmitter impairment mitigation and monitoring for high baud-rate, high order modulation systems,” in Proc. Eur. Conf. Opt. Commun., Düsseldorf, Germany, 2016, pp. 1–3.

25. P. Skvortcov, C. Sanchez-Costa, I. Phillips, and W. Forysiak, “Receiver DSP highly tolerant to transmitter IQ impairments,” in Optical Fiber Communication Conference (OFC) 2019, OSA Technical Digest (Optica Publishing Group, 2019), paper Th1D.2.

26. M. Sato, M. Arikawa, H. Noguchi, J. Matsui, J. ‘i. Abe, and E. Le Taillandier de Gabory, “Mitigation of Transmitter Impairment with 4×2 WL MIMO Equalizer Embedding Preliminary CPR,” in Optical Fiber Communication Conference (OFC)2022, paper M1H.5.

27. E. P. D. Silva and D. Zibar, “Widely Linear Blind Adaptive Equalization for Transmitter IQ-Imbalance/Skew Compensation in Multicarrier Systems,” ECOC 2016; 42nd European Conference on Optical Communication, Dusseldorf, Germany, 2016, pp. 1–3.

28. M. S. Neves, A. Lorences-Riesgo, C. S. Martins, S. Mumtaz, Y. Frignac, G. Charlet, P. P. Monteiro, S. Dris, and F. P. Guiomar, “Leveraging Dispersion-Aware Phase Recovery for Long-Haul Digital Multi-Carrier Transmission: An Experimental Demonstration,” J. Lightwave Technol. 40(16), 5432–5439 (2022). [CrossRef]  

29. S. M. Bilal, C. Fludger, and G. Bosco, “Carrier phase estimation in multisubcarrier coherent optical systems,” IEEE Photonics Technol. Lett. 28(19), 2090–2093 (2016). [CrossRef]  

30. M. S. Neves, A. Carena, A. Nespola, P. P. Monteiro, and F. P. Guiomar, “Joint Carrier-Phase Estimation for Digital Subcarrier Multiplexing Systems With Symbol-Rate Optimization,” J. Lightwave Technol. 39(20), 6403–6412 (2021). [CrossRef]  

31. D. Martinez, Y. Mori, H. Hasegawa, and K. Sato, “Novel Subcarrier Multiplexing and Subcarrier-Synchronous Phase Estimation Tolerant to Laser Phase Noise,” in Advanced Photonics 2017 (IPR, NOMA, Sensors, Networks, SPPCom, PS), OSA Technical Digest (online) (Optica Publishing Group, 2017), paper SpW1F.5.

32. M. Saifuddin Faruk and S. J. Savory, “Digital Signal Processing for Coherent Transceivers Employing Multilevel Formats,” J. Lightwave Technol. 35(5), 1125–1141 (2017). [CrossRef]  

33. W. Tong, J. Zhang, M. Zhu, X. Liu, Y. Wei, B. Hua, Z. Xin, M. Lei, Y. Cai, Y. Zou, L. Tian, and G. Zhao, “Experimental Comparison of Single-Carrier and Digital Subcarrier Multiplexing Transmissions in a W-Band 200 Gb/s Fiber-Wireless System Considering Transmitter IQ Imbalance and Skew Mitigation,” IEEE Photonics J. 15(4), 1–8 (2023). [CrossRef]  

34. M. S. Neves, P. P. Monteiro, and F. P. Guiomar, “Enhanced Phase Estimation for Long-Haul Multi-Carrier Systems Using a Dual-Reference Subcarrier Approach,” J. Lightwave Technol. 39(9), 2714–2724 (2021). [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 (9)

Fig. 1.
Fig. 1. Illustration of Tx IQ skew effect in a DSCM signal with eight subcarriers.
Fig. 2.
Fig. 2. The structure of four equalization configurations. Configuration I: 2 × 2 equalizer for de-polarization, phase noise compensation based pilot symbols; Configuration II: 2 × 2 equalizer for de-polarization, phase noise compensation based pilot symbols and Tx IQ skew compensation for each subcarrier by 2 × 2 equalizer independent; Configuration III: 2 × 2 equalizer for de-polarization and Tx skew compensation for each subcarrier by 2 × 2 equalizer combined with PLL; Configuration IV: PLL embedded 4 × 4 equalizer for joint equalization of polarization, phase noise and Tx IQ skew.
Fig. 3.
Fig. 3. (a)Simulation setup. PBS: polarization beam splitter, PBC: polarization beam coupler, SSMF: standard single mode fiber, EDFA: Erbium-doped optical fiber amplifier. (b)The spectrum of DSCM after the effect of Tx IQ skew (c) Tx DSP and Rx DSP.
Fig. 4.
Fig. 4. Q-factor Vs Tx IQ skew for different equalization configurations in DSCM system (a) optical back-to-back without phase noise (b) transmission over 100 km SSMF with phase noise.
Fig. 5.
Fig. 5. (a) Q-factor Vs Tx IQ skew of PLL embedded 4 × 4 equalizer in independent and shared mode (b) Q-factor versus transmission distance of the overall performance of DSCM system with the PLL embedded 4 × 4 equalizer in independent mode and shared mode with different reference subcarriers.
Fig. 6.
Fig. 6. Q-factor Vs OSNR for comparison of single-carrier and DSCM systems with/without 3ps Tx skew after transmission over 100 km SSMF.
Fig. 7.
Fig. 7. Q-factor Vs transmission distance of single-carrier and DSCM systems with/without 3ps Tx IQ skew.
Fig. 8.
Fig. 8. Experimental setup and offline DSP.
Fig. 9.
Fig. 9. Q-factor Vs Tx Y-IQ skew for different algorithms in DSCM system transmission over 100 km SSMF.

Equations (9)

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

[ I t ( t ) Q t ( t ) ] = [ 1 0 0 δ ( t τ ) ] [ I ( t ) Q ( t ) ]
y ( t ) = x ( t ) + x ( t ) 2 + j ( δ ( t τ ) ( j ( x ( t ) x ( t ) ) 2 ) ) = k 1 ( t ) x ( t ) + k 2 ( t ) x ( t )
Y ( w ) = 1 2 [ X ( w ) ( 1 + e j w τ ) + X ( w ) ( 1 e j w τ ) ]
x ( t ) = n = 1 N x n ( t ) e j 2 π f n t
Y ( w ) = 1 2 n = 1 N ( X n ( w 2 π f n ) ( 1 + e j w τ ) + X n ( w 2 π f n ) ( 1 e j w τ ) )   = 1 2 n = 1 N ( X n ( w 2 π f n ) ( 1 + cos w τ j sin w τ ) + X n ( w 2 π f n ) ( 1 cos w τ + j sin w τ ) )
[ X n Y n ] = [ h x n x n h x n y n h y n x n h y n y n ] [ x n e i θ n y n e i θ n ] + [ h x n x m h x n y m h y n x m h y n y m ] [ x m e i θ m y m e i θ m ] [ X m Y m ] = [ h x m x m h x m y m h y m x m h y m y m ] [ x m e i θ m y m e i θ m ] + [ h x m x n h x m y n h y m x n h y m y n ] [ x n e i θ n y n e i θ n ]
e x m ( k ) = T x m ( k ) X m ( k ) , e y m ( k ) = T y m ( k ) Y m ( k )
J x m ( h ) = ( e x m ( k ) ) 2 , J y m ( h ) = ( e y m ( k ) ) 2
h x m x m ( k + 1 ) = h x m x m ( k ) + μ e x m x m ( k ) e i θ m , h x m y m ( k + 1 ) = h x m y m ( k ) + μ e x m y m ( k ) e i θ m h x m x n ( k + 1 ) = h x m x n ( k ) + μ e x m x n ( k ) e i θ n , h x m y n ( k + 1 ) = h x m y n ( k ) + μ e x m y n ( k ) e i θ n h y m x m ( k + 1 ) = h y m x m ( k ) + μ e y m x m ( k ) e i θ m , h y m y m ( k + 1 ) = h y m y m ( k ) + μ e y m y m ( k ) e i θ m h y m x n ( k + 1 ) = h y m x n ( k ) + μ e y m x n ( k ) e i θ n , h y m y n ( k + 1 ) = h y m y n ( k ) + μ e y m y n ( k ) e i θ n
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.