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

Performance analysis for OFDM-based multi-carrier continuous-variable quantum key distribution with an arbitrary modulation protocol

Open Access Open Access

Abstract

Multi-carrier continuous-variable quantum key distribution (CV-QKD) is considered to be a promising way to boost the secret key rate (SKR) over the existing single-carrier CV-QKD scheme. However, the extra excess noise induced in the imperfect multi-carrier quantum state preparation process of N subcarriers will limit the performance of the system. Here, a systematic modulation noise model is proposed for the multi-carrier CV-QKD based on the orthogonal frequency division multiplexing (OFDM). Subsequently, the performance of multi-carrier CV-QKD with arbitrary modulation protocol (e.g. QPSK, 256QAM and Gaussian modulation protocol) can be quantitatively evaluated by combining the security analysis method of the single-carrier CV-QKD. Under practical system parameters, the simulation results show that the SKR of the multi-carrier CV-QKD can still be significantly improved by increasing the carrier number N even with imperfect practical modulations. Specifically, the total SKR of multi-carrier CV-QKD can be optimized by carefully choosing N. The proposed model provides a feasible theoretical framework for the future multi-carrier CV-QKD experimental implementation.

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

1. Introduction

Quantum key distribution (QKD) allows to distribute an information-theoretically secure key between two legal communication parties through an insecure channel based on the fundamental principles of quantum mechanics [14]. Compared with the intensively developed discrete-variable QKD (DV-QKD) based on single-photon detection, continuous-variable QKD (CV-QKD) based on coherent detection is compatible with commercial off-the-shelf components and has the potential advantages of high secret key rate (SKR) in metropolitan area [59].

Currently, single-frequency optical carrier scheme is regularly used in most existing CV-QKD setups, and it has been classified into two schemes according to the modulation method of the coherent state, i.e., Gaussian modulation coherent state (GMCS) [1012] and discrete modulation coherent state (DMCS) [1316]. For the single-carrier GMCS CVQKD, the SKR has achieved a sub-Mbps level over 100 km [17]. Meanwhile, the SKR of the single-carrier DMCS CVQKD has reached a sub-Gbps level within metropolitan area [18,19]. Specifically, the SKRs of these single-carrier schemes can be improved by increasing the system repetition rate and optimizing the excess noise. However, the repetition rate of the state-of-art CV-QKD system has reached 5 GHz [18], where the technical challenge lies in that it requires faster data acquisition cards with high linearity, wider bandwidth homodyne detectors with low noise and faster post-processing with high complexity [2024]. Meanwhile, the excess noise of CV-QKD system has been controlled to a reasonably low level at the order of 0.001∼0.01 [18,25,26]. To further improve the SKR, the multi-carrier CV-QKD scheme has been proposed to distribute multiplexing independent secret keys encoded on N subcarriers within a single fiber channel [27,28]. At present, two multi-carrier CV-QKD schemes have been proposed. One needs multiple CV-QKD transmitters, which is basically the same as deploying multiple CV-QKD systems [29]. The other is based on one transmitter with N independent radio-frequency (RF) oscillators while it is difficult to realize due to the complicated electrical structure [30]. In fact, the multi-carrier CV-QKD scheme can be simply realized with one transmitter and one receiver by employing the orthogonal-frequency-division-multiplexing (OFDM) method [3133], which is efficient and economical for practical deployment. Moreover, the SKR can be intuitively improved by increasing the carrier number N for an OFDM-based multi-carrier CV-QKD. However, extra excess noise will be introduced in the imperfect multi-carrier quantum state preparation process, which will deteriorate the system performance for larger N. So, to effectively improve the system performance, there is a trade-off between the carrier number N and the extra excess noise. The extra excess noise is mainly the modulation noise that originated from the IQ imbalance and intermodulation distortion. In [31], the IQ imbalance noise has been modeled and analyzed for the QPSK multi-carrier CV-QKD scheme. Nevertheless, there is still a lack of systematic modulation noise model for multi-carrier CV-QKD with arbitrary modulation protocol (e.g. QPSK, 256QAM or Gaussian protocol), which is seriously needed to evaluate the performance of the scheme.

In this paper, a systematic modulation noise model is proposed for the OFDM-based multi-carrier CV-QKD with arbitrary modulation protocol. Based on the proposed modulation noise model, the SKRs of the OFDM-based multi-carrier CV-QKD system with QPSK, 256QAM and Gaussian modulation protocols are quantitatively evaluated by combining the security analysis method of the single-carrier CV-QKD. Specifically, the modulation noise from IQ imbalance and third-order intermodulation distortion are quantitatively modeled and analyzed, which limits the performance of the system with N subcarriers. Under practical system parameters, the simulation results show that the OFDM-based multi-carrier CV-QKD system can achieve much higher asymptotic SKR by increasing N even with imperfect modulation compared with the single-carrier CV-QKD scheme. Moreover, a maximum total SKR for multi-carrier CV-QKD can be obtained by carefully optimizing N. Therefore, our work offers a promising way to further improve the SKR performance of quantum secure network.

2. Description of the multi-carrier CV-QKD scheme

The working process of the OFDM-based multi-carrier CV-QKD scheme can be described as follows. (i) Quantum state preparation: Alice modulates N independent random keys on a N mode coherent state as based on QPSK, 256QAM or Gaussian modulation CV-QKD protocol, and sends the prepared quantum state to Bob through the quantum channel. (ii) Quantum state measurement: Bob measures the received quantum state with coherent receiver (e.g. homodyne or heterodyne detection), and gets the measurement results $\vec{y} = \{ {y_k}|k = 1,2, \ldots ,N\} $ based on digital demodulation. Alice and Bob repeat steps (i) and (ii) for n times. (iii) Post-processing: Alice and Bob perform sifting (if only Bob uses homodyne detection), parameter estimation, reverse reconciliation and privacy amplification independently for N modes of encoded quantum states. The total asymptotic SKR between Alice and Bob with N subcarriers can be written as:

$$R = \sum\limits_{k = 1}^N {{R_k}} ,\quad {R_k} = {\beta _k}I({{A_k}:{B_k}} )- \chi ({{B_k}:E} )$$
where Rk is the asymptotic SKR of k-th mode, I(Ak:Bk) is the Shannon mutual entropy between Alice and Bob for k-th mode, χ(Bk:E) is the quantum mutual entropy between Bob and Eve for k-th mode, and ${\beta _k}$ is the reconciliation efficiency for k-th mode.

In the following, we describe detailly the OFDM-based multi-carrier CV-QKD scheme with local local oscillator (LLO), as shown in Fig. 1. At Alice’s site, a continuous optical carrier with central frequency fA that generated from Alice laser is split into two branches. In the upper branch, the optical carrier is modulated in an IQ modulator by OFDM signal Is and Qs generated from an OFDM generator. The processing routine of the OFDM generator is demonstrated in Fig. 2. Firstly, a serial of true random bit stream with high bit rate originated from quantum random number generator (QRNG) is converted to N independent parallel bit streams by a serial-to-parallel conversion (S/P). Secondly, the N parallel bit streams are mapped with arbitrary modulation distribution based on the chosen CV-QKD protocol, such as QPSK, 256QAM or Gaussian distribution. After inverse fast Fourier transform (IFFT), the processed parallel data is converted to the desired OFDM signal by a parallel-to-serial conversion (P/S), given by

$${I_s} = \sum\limits_{k = 1}^N {{I_k}\cos ({2\pi {f_k}t} )- } \sum\limits_{k = 1}^N {{Q_k}\sin ({2\pi {f_k}t} )}$$
$${Q_s} = \sum\limits_{k = 1}^N {{I_k}\sin ({2\pi {f_k}t} )+ } \sum\limits_{k = 1}^N {{Q_k}\cos ({2\pi {f_k}t} )}$$
where N is the total number of the generated subcarriers, fk = kΔf corresponds to the frequency of k-th subcarrier and Δf denotes the frequency spacing between N subcarriers. It is noteworthy from Fig. 2 that the cyclic prefix (CP) and training sequence (TS) are appended into the OFDM signal in processing routine for avoiding mutual crosstalk between subcarriers in fiber channel and achieving the data-aided time domain equalization at Bob’s site.

 figure: Fig. 1.

Fig. 1. Schematic diagram of the OFDM-based multi-carrier CV-QKD with local local oscillator. LD: laser diode; OFDM: orthogonal frequency division multiplexing; MZM: Mach-Zehnder modulator; PBS: polarization beam splitter; VOA, variable optical attenuator; PBC: polarization beam combiner; SMF: single mode fiber; PSA: polarization synthesis analyzer; OC: optical coupler; BPD: balanced photo-detector; DSP: digital signal processing

Download Full Size | PDF

 figure: Fig. 2.

Fig. 2. The data processing of the OFDM generator. S/P: serial-to-parallel conversion; TS: training sequence; IFFT: inverse fast Fourier transform; P/S: parallel-to-serial conversion; CP: cyclic prefix; DAC: digital-to-analog converter.

Download Full Size | PDF

In the upper branch at Alice’s site in Fig. 1, the OFDM-based optical signal after IQ modulation is attenuated by a variable optical attenuator (VOA) to generate the encoded multi-carrier quantum state

$$|{{X_{sig}} + j{P_{sig}}} \rangle ={\otimes} _{k = 1}^N|{{X_k} + j{P_k} + \Delta {X_k} + j\Delta {P_k}} \rangle$$
where Xk and Pk are the expected modulated quadrature components for the k-th subcarrier. ΔXk and ΔPk are the extra modulation noise on Xk and Pk, respectively, which mainly come from IQ imbalance and intermodulation distortion in the multi-carrier IQ modulation process. The detailed theoretical model of the modulation process is shown in next section. In the lower branch, the optical carrier is sent into a frequency shifter (e.g. MZM with carrier suppression double sideband modulation) with fs for generating the desired pilot tone with frequency fA-fs [11]. The generated pilot tone is used for eliminating the fast-drift phase noise in the LLO-CV-QKD system. Significantly, the intense pilot tone and the weak quantum signal are respectively generated in two different optical paths, which is beneficial to improve the preparation accuracy of the multi-carrier quantum state in the case of finite DAC quantization bits and modulation extinction ratio (ER). Finally, the multi-carrier quantum signal and the pilot tone are co-transmitted into a single mode fiber (SMF) channel with different frequency bands and orthogonal polarization states for avoiding the crosstalk between them.

At Bob’s site, the received quantum signal and pilot tone are separated by a polarization beam splitter (PBS). In order to separate the multi-carrier quantum signal and pilot tone efficiently, a polarization synthesis analyzer (PSA) is used for correcting the polarization deterioration resulted from the fiber channel disturbance. The separated quantum signal and pilot tone are independently coupled into two balanced photo-detectors (BPDs) with their corresponding LLO signals of frequency fB supplied by the Bob laser. The detected quantum signal and pilot tone are analog-to-digital converted for the subsequent DSP, respectively. As is shown in Fig. 3, the frequency fsfAB is firstly efficiently estimated from the pilot tone, with which the central frequency ΔfAB of the OFDM-based quantum signal can be obtained. Secondly, the OFDM-based quantum signal ysig and the pilot tone ypilot are bandpass filtered, down-converted and matched filtered for recovering their quadrature components. Thirdly, the fast-drift phase noise of the OFDM-based quantum signal can be compensated based on the pilot-assisted channel equalization method, which is recently applied in optical fiber LLO-CV-QKD system [34]. After removing CP, the serial OFDM-based quantum signal is converted to N parallel signals and then transformed to N independent subcarriers by FFT. Subsequently, the phase noises of N independent subcarriers are further compensated based on data-aided time domain equalization method, respectively. Finally, a high-efficient post-processing setup for the multi-carrier CV-QKD is required to extract N independent secure keys. For distilling the final key based on the security analysis method of the single-carrier CV-QKD, the excess noise of the designed multi-carrier CV-QKD system, mainly including laser intensity noise, DAC quantization noise, modulation noise, fiber channel crosstalk noise, and phase noise, should be controlled to be a reasonably low level. Moreover, the shot noise for each subcarrier should be calibrated in real time at Bob’s site [35].

 figure: Fig. 3.

Fig. 3. The DSP routine of the OFDM-based multi-carrier LLO-CV-QKD scheme. CP: cyclic prefix; S/P: serial-to- parallel conversion; FFT: fast Fourier transform.

Download Full Size | PDF

3. Modulation noise model of the multi-carrier CV-QKD

In the OFDM-based multi-carrier modulation process, the preparation accuracy of the multi-carrier quantum state will suffer from the IQ imbalance and intermodulation distortion. Therefore, a detailed modulation noise model of the multi-carrier CV-QKD scheme should be derived to effectively evaluate the SKR.

The OFDM-based multi-carrier quantum signal after the IQ modulation can be written as:

$${E_{sig}}(t )= 2{A_{sig}}{e^{j2\pi {f_A}t}}\left\{ {{G_1}\sin \left[ {\sum\limits_{k = 1}^N {{\gamma_k}\cos ({2\pi {f_k}t + {\varphi_k}} )} } \right] + j{G_2}\sin \left[ {\sum\limits_{k = 1}^N {{\gamma_k}\sin ({2\pi {f_k}t + {\varphi_k}} )} } \right]} \right\}$$
with
$${\gamma _k} = {\mu _k}\sqrt {{I_k}^2 + {Q_k}^2}$$
$$\cos {\varphi _k} = \frac{{{I_k}}}{{\sqrt {{I_k}^2 + {Q_k}^2} }}$$
where the IQ imbalance factor G1 = [1 + κkexp(k)]/2 and G2 = [1 + κkexp(-k)]/2, and κk and θk are the gain imbalance and quadrature skew of the k-th subcarrier modulation, respectively. µk = Vkπ/Vπ/2 is the modulation index of the k-th subcarrier. Vk and Vπ are the driving amplitude and the half-wave voltage of the employed I/Q modulator, respectively, and they should be chosen as a reasonably lower value for guaranteeing the small signal condition and the multi-carrier quantum state with enough signal-to-noise ratio (SNR) in practical experiment. Asig is the amplitude of the quantum signal. Significantly, under IQ balance (κk = 1 and θk = 0) and no intermodulation distortion, Eq. (4) can be simplified as:
$${E_{sig}}(t )= 2{A_{sig}}{e^{j2\pi {f_A}t}}\sum\limits_{k = 1}^N {{\gamma _k}{e^{j({2\pi {f_k}t + {\varphi_k}} )}}} = 2{A_{sig}}{e^{j2\pi {f_A}t + j2\pi {f_k}t}} \otimes _{k = 1}^N|{{I_k} + j{Q_k}} \rangle$$

Comparing Eq. (3) and Eq. (6), one can see that the expected modulated quadrature Xk and Pk of the k-th subcarrier are determined by the random variables Ik and Qk with arbitrary modulation distribution.

In order to derive modulation noise model for ΔXk and ΔPk shown in Eq. (3), we apply the Jacobi-Anger expansion [36] into Eq. (4), which can be expanded with the pth-order Bessel function of the first kind Jp(·) (p = 0, ± 1, ± 2, ± 3, …) as follows:

$${E_{sig}}(t )= 2{A_{sig}}{e^{j2\pi {f_A}t}}\left\{ \begin{array}{l} {G_1}\sum\limits_{{p_1} ={-} \infty }^{ + \infty } {{J_{{p_1}}}({{\gamma_1}} )\cdots \sum\limits_{{p_k} ={-} \infty }^{ + \infty } {{J_{{p_k}}}({{\gamma_k}} )} } \cdots \sum\limits_{{p_N} ={-} \infty }^{ + \infty } {({{\gamma_{_N}}} ){J_{{p_N}}}} \cdot \\ \sin \left[ \begin{array}{l} \left( {2\pi {f_1}t{p_1} + {p_1}{\varphi_1} + {p_1}\frac{\pi }{2}} \right) +\cdots{+} \left( {2\pi {f_k}t{p_k} + {p_k}{\varphi_{_k}} + {p_k}\frac{\pi }{2}} \right)\\ +\cdots{+} \left( {2\pi {f_N}t{p_N} + {p_N}{\varphi_{_N}} + {p_N}\frac{\pi }{2}} \right) \end{array} \right]\\ + j{G_2}\sum\limits_{{p_1} ={-} \infty }^{ + \infty } {{J_{{p_1}}}({{\gamma_1}} )\sum\limits_{{p_k} ={-} \infty }^{ + \infty } {{J_{{p_k}}}({{\gamma_2}} )} } \cdots \sum\limits_{{p_N} ={-} \infty }^{ + \infty } {{J_{{p_N}}}({{\gamma_{_N}}} )} \cdot \\ \sin \left[ \begin{array}{l} ({2\pi {f_1}t{p_1} + {p_1}{\varphi_1}} )+\cdots{+} ({2\pi {f_k}t{p_k} + {p_k}{\varphi_{_k}}} )\\ +\cdots{+} ({2\pi {f_N}t{p_N} + {p_N}{\varphi_{_N}}} )\end{array} \right] \end{array} \right\}$$

From Eq. (7), we can quantify the frequency component at f0 + fk for expressing the k-th subcarrier in consideration of the IQ imbalance and the third-order intermodulation distortion, which is given by

$$\begin{array}{l} {E_{sig(k)}}(t )= 2{A_{sig}}{e^{j\pi {f_A}t}} \cdot \\ \left\{ \begin{array}{l} {J_1}({{\gamma_k}} )[{({{G_1} + {G_2}} ){e^{j({2\pi {f_k}t + {\varphi_k}} )}} + ({{G_1} - {G_2}} ){e^{ - j({2\pi {f_k}t + {\varphi_k}} )}}} ]\\ - \sum\limits_{,n = 1\atop { \ne n \ne k\atop \scriptstyle2m + n = k}}^N {{J_2}({{\gamma_m}} ){J_1}({{\gamma_n}} )[{({{G_1} + {G_2}} ){e^{ - j({2\pi {f_k}t + 2{\varphi_m} + {\varphi_n}} )}} + ({{G_1} - {G_2}} ){e^{j({2\pi {f_k}t + 2{\varphi_m} + {\varphi_n}} )}}} ]} \\ - \sum\limits_{,s = 1\atop { \ne n \ne k\atop \scriptstyle2m - n = k}}^N {{J_2}({{\gamma_m}} ){J_1}({{\gamma_n}} )[{({{G_1} + {G_2}} ){e^{j({2\pi {f_k}t + 2{\varphi_m} - {\varphi_n}} )}} + ({{G_1} - {G_2}} ){e^{ - j({2\pi {f_k}t + 2{\varphi_m} - {\varphi_n}} )}}} ]} \\ - \sum\limits_{,s,l = 1\atop { \ne n \ne l \ne k\atop + n + l = k}}^N {{J_1}({{\gamma_m}} ){J_1}({{\gamma_n}} ){J_1}({{\gamma_l}} )[{({{G_1} + {G_2}} ){e^{ - j({2\pi {f_k}t + {\varphi_m} + {\varphi_n} + {\varphi_l}} )}} + ({{G_1} - {G_2}} ){e^{j({2\pi {f_k}t + {\varphi_m} + {\varphi_n} + {\varphi_l}} )}}} ]} \\ - \sum\limits_{,n,l = 1\atop { \ne n \ne l \ne k\atop + n - l = k}}^N {{J_1}({{\gamma_m}} ){J_1}({{\gamma_n}} ){J_1}({{\gamma_l}} )[{({{G_1} + {G_2}} ){e^{j({2\pi {f_k}t + {\varphi_m} + {\varphi_n} - {\varphi_l}} )}} + ({{G_1} - {G_2}} ){e^{ - j({2\pi {f_k}t + {\varphi_m} + {\varphi_n} - {\varphi_l}} )}}} ]} \\ - \sum\limits_{,n,l = 1\atop { \ne n \ne l \ne k\atop - n - l = k}}^N {{J_1}({{\gamma_m}} ){J_1}({{\gamma_n}} ){J_1}({{\gamma_l}} )[{({{G_1} + {G_2}} ){e^{ - j({2\pi {f_k}t + {\varphi_m} - {\varphi_n} - {\varphi_l}} )}} + ({{G_1} - {G_2}} ){e^{j({2\pi {f_k}t + {\varphi_m} - {\varphi_n} - {\varphi_l}} )}}} ]} \end{array} \right\} \end{array}$$

Under small signal condition of J1(x)≈x/2 and J2(x)≈x2/4, the quadrature Xsig(k) of the k-th subcarrier can be obtained from Eq. (8) and expressed as:

$$\begin{array}{l} {X_{sig(k )}} = {X_k} + \Delta {X_k} = 2{A_{sig}}\textrm{real} ({{\gamma_k}{e^{j{\varphi_k}}}} )\\ + 2{A_{sig}}real \left\{ \begin{array}{l} \frac{{{\gamma_k}}}{2}[{({{G_1} + {G_2} - 2} ){e^{j{\varphi_k}}} + ({{G_1} - {G_2}} ){e^{ - j{\varphi_k}}}} ]\\ - \frac{{{M_1}({N,k} )\gamma_m^2{\gamma_n}}}{8}[{({{G_1} + {G_2}} ){e^{ - j({2{\varphi_m} + {\varphi_n}} )}} + ({{G_1} - {G_2}} ){e^{j({2{\varphi_m} + {\varphi_n}} )}}} ]\\ - \frac{{{M_2}({N,k} )\gamma_m^2{\gamma_n}}}{8}[{({{G_1} + {G_2}} ){e^{j({2{\varphi_m} - {\varphi_n}} )}} + ({{G_1} - {G_2}} ){e^{ - j({2{\varphi_m} - {\varphi_n}} )}}} ]\\ - \frac{{{W_1}({N,k} ){\gamma_m}{\gamma_n}{\gamma_l}}}{8}[{({{G_1} + {G_2}} ){e^{ - j({{\varphi_m} + {\varphi_n} + {\varphi_l}} )}} + ({{G_1} - {G_2}} ){e^{j({{\varphi_m} + {\varphi_n} + {\varphi_l}} )}}} ]\\ - \frac{{{W_2}({N,k} ){\gamma_m}{\gamma_n}{\gamma_l}}}{8}[{({{G_1} + {G_2}} ){e^{j({{\varphi_m} + {\varphi_n} - {\varphi_l}} )}} + ({{G_1} - {G_2}} ){e^{ - j({{\varphi_m} + {\varphi_n} - {\varphi_l}} )}}} ]\\ - \frac{{{W_3}({N,k} ){\gamma_m}{\gamma_n}{\gamma_l}}}{8}[{({{G_1} + {G_2}} ){e^{ - j({{\varphi_m} - {\varphi_n} - {\varphi_l}} )}} + ({{G_1} - {G_2}} ){e^{j({{\varphi_m} - {\varphi_n} - {\varphi_l}} )}}} ]\end{array} \right\} \end{array}$$
where M1(N, k) and M2(N, k) are the combinatorial numbers of the condition 2m + n = k and 2m-n = k, respectively. W1(N, k), W2(N, k) and W3(N, k) are the combinatorial numbers of the condition m + n + l = k, m + n-l = k and m-n-l = k, respectively. In Eq. (9), Xk is the k-th subcarrier component under the condition of IQ balance (κk = 1 and θk = 0) and no intermodulation distortion [M1(N, k)=M2(N, k)=W1(N, k)=W2(N, k)=W3(N, k) = 0], which can be expressed as
$${X_k} = 2{A_{sig}}{\mu _k}{I_k}$$

In Eq. (9), ΔXk is the extra modulation noise introduced from the IQ imbalance and third-order intermodulation distortion, which includes three parts

$$\Delta {X_k} = \Delta {X_{k1}} + \Delta {X_{k2}} + \Delta {X_{k3}}$$
with
$$\Delta {X_{k1}} = {A_{sig}}{\mu _k}[{({{\kappa_k}\cos {\theta_k} - 1} ){I_k} + ({{\kappa_k}\sin {\theta_k}} ){Q_k}} ]$$
$$\Delta {X_{k2}} = \frac{1}{4}{A_{sig}}\mu _m^2{\mu _n}\left\{ \begin{array}{l} ({1 + {\kappa_k}\cos {\theta_k}} )\left\{ \begin{array}{l} [{{M_1}({N,k} )+ {M_2}({N,k} )} ]Q_m^2{I_n}\\ + 2[{{M_1}({N,k} )- {M_2}({N,k} )} ]{I_m}{Q_m}{Q_n}\\ - [{{M_1}({N,k} )+ {M_2}({N,k} )} ]I_m^2{I_n} \end{array} \right\}\\ + ({{\kappa_k}\sin {\theta_k}} )\left[ \begin{array}{l} [{{M_1}({N,k} )+ {M_2}({N,k} )} ]I_m^2{Q_n}\\ + 2[{{M_1}({N,k} )- {M_2}({N,k} )} ]{I_m}{Q_m}{I_n}\\ - [{{M_1}({N,k} )+ {M_2}({N,k} )} ]Q_m^2{Q_n} \end{array} \right] \end{array} \right\}$$
$$\Delta {X_{k3}} = \frac{1}{4}{A_{sig}}{\mu _m}{\mu _n}{\mu _l}\left\{ \begin{array}{l} ({1 + {\kappa_k}\cos {\theta_k}} )\left\{ \begin{array}{l} [{{W_1}({N,k} )- {W_2}({N,k} )+ {W_3}({N,k} )} ]{I_m}{Q_n}{Q_l}\\ + [{{W_1}({N,k} )- {W_2}({N,k} )- {W_3}({N,k} )} ]{Q_m}{I_n}{Q_l}\\ + [{{W_1}({N,k} )+ {W_2}({N,k} )- {W_3}({N,k} )} ]{Q_m}{Q_n}{I_l}\\ - [{{W_1}({N,k} )+ {W_2}({N,k} )+ {W_3}({N,k} )} ]{I_m}{I_n}{I_l} \end{array} \right\}\\ + ({{\kappa_k}\sin {\theta_k}} )\left\{ \begin{array}{l} [{{W_1}({N,k} )+ {W_2}({N,k} )- {W_3}({N,k} )} ]{I_m}{I_n}{Q_l}\\ + [{{W_1}({N,k} )- {W_2}({N,k} )- {W_3}({N,k} )} ]{I_m}{Q_n}{I_l}\\ + [{{W_1}({N,k} )- {W_2}({N,k} )+ {W_3}({N,k} )} ]{Q_m}{I_n}{I_l}\\ - [{{W_1}({N,k} )+ {W_2}({N,k} )+ {W_3}({N,k} )} ]{Q_m}{Q_n}{Q_l} \end{array} \right\} \end{array} \right\}$$
where ΔXk1 is modulation noise resulted from the IQ imbalance. The modulation noise ΔXk2 and ΔXk3 are originated from the third-order intermodulation with the frequency relationships (2fm-fn = fk and 2fm + fn = fk) and (fm + fn + fl = fk, fm + fn-fl = fk and fm-fn-fl = fk), respectively.

It is obvious that the quadrature Ig and Qg (g = k, m, n, l) of the g-th subcarrier are mutually random and independent. One can define the variance 〈Ig2〉=〈Qg2〉=σ12 and 〈Ig4〉=〈Qg4〉=σ22. According to Eq. (11), the variances of ΔXkcan be expressed as

$$\left\langle {\Delta X_k^2} \right\rangle = \left\langle {\Delta X_{k1}^2} \right\rangle + \left\langle {\Delta X_{k2}^2} \right\rangle + \left\langle {\Delta X_{k3}^2} \right\rangle$$
with
$$\left\langle {\Delta X_{k1}^2} \right\rangle = A_{sig}^2\mu _k^2\sigma _1^2({\kappa_k^2 + 1 - 2{\kappa_k}\cos {\theta_k}} )$$
$$\left\langle {\Delta X_{k2}^2} \right\rangle = \frac{{A_{sig}^2\mu _m^4\mu _n^2}}{8}({1 + 2{\kappa_k}\cos {\theta_k} + {\kappa_k}^2} )\left\{ \begin{array}{l} {[{{M_1}({N,k} )+ {M_2}({N,k} )} ]^2}\sigma_2^2\sigma_1^2\\ + 2{[{{M_1}({N,k} )- {M_2}({N,k} )} ]^2}\sigma_1^6 \end{array} \right\}$$
$$\left\langle {\Delta X_{k3}^2} \right\rangle = \frac{{A_{sig}^2\mu _m^2\mu _n^2\mu _l^2}}{4}({1 + 2{\kappa_k}\cos {\theta_k} + \kappa_k^2} )[{W_1^2({N,k} )+ W_2^2({N,k} )+ W_3^2({N,k} )} ]\sigma _1^6$$

At the same time, µk, µm, µn and µl are the modulation indexes of the k, m, n, l-th subcarrier modulation, respectively, which are approximately equal to each other in the multi-carrier modulation process. So, we can use µk to replace µm, µn and µl in Eqs. (14b) and (14c). Moreover, according to 〈Ik2〉=〈Qk2〉=σ12 and the variance of Eq. (10), the expected modulation variance of the k-th subcarrier can be expressed as

$${V_A} = \left\langle {X_k^2} \right\rangle = 4A_{sig}^2\mu _k^2\left\langle {I_k^2} \right\rangle = 4A_{sig}^2\mu _k^2\sigma _1^2$$
where the modulation variance VA is normalized in shot noise unit (SNU). From Eq. (15), it is obvious that the modulation index µk is directly related to the modulation variance VA of the multi-carrier quantum state. Based on Eqs. (13), (14) and (15), we can simplify the modulation noise of the k-th subcarrier with the modulation variance VA, given by
$$\begin{array}{l} {\varepsilon _{mod}}(k) = \left\langle {\Delta X_k^2} \right\rangle = \frac{{{V_A}}}{4}({\kappa_k^2 + 1 - 2{\kappa_k}\cos {\theta_k}} )\\ + \frac{{{V_A}\mu _k^4}}{{32}}({1 + 2{\kappa_k}\cos {\theta_k} + {\kappa_k}^2} )\left\{ \begin{array}{l} {[{{M_1}({N,k} )+ {M_2}({N,k} )} ]^2}\sigma_2^2\\ + 2\left\{ \begin{array}{l} {[{{M_1}({N,k} )- {M_2}({N,k} )} ]^2}\\ \textrm{ + }W_1^2({N,k} )+ W_2^2({N,k} )+ W_3^2({N,k} )\end{array} \right\}\sigma_1^4 \end{array} \right\} \end{array}$$

One can see from Eq. (16) that the modulation noise of the k-th subcarrier in the multi-carrier CV-QKD system is related with the modulation index µk, the IQ imbalance factors (κk and θk) and the third-order intermodulation factors [M1(N, k), M2(N, k), W1(N, k), W2(N, k) and W3(N, k)]. Moreover, the modulation noise is also determined by the variance σ12 and σ22 related with different modulation distributions (e.g. QPSK, 256QAM or Gaussian modulation) in the quantum state preparation process. One can see from Eq. (16) that the modulation noise of the single-carrier CVQKD can be also obtained, which mainly comes from the IQ imbalance noise where it has no the intermodulation distortion, i.e., [M1(N, k)=M2(N, k)=W1(N, k)=W2(N, k)=W3(N, k) = 0].

4. Performance evaluation of multi-carrier CV-QKD

It is well known that the modulation noise is estimated as a part of excess noise at Alice’s site [37]. In this case, the extra modulation noise εmod(k) of the k-th subcarrier in the multi-carrier CV-QKD will cause an increment of the excess noise compared to a single-carrier CV-QKD, and the total excess noise εmulti(k) of the k-th subcarrier can be expressed as

$${\varepsilon _{multi}}(k )= {\varepsilon _{mod}}(k )+ {\varepsilon _{single}}$$
where εsingle is the excess noise of a single-carrier CV-QKD scheme without the modulation noise. Based on the asymptotic security analysis method of the single-carrier CV-QKD scheme, the asymptotic SKR performances of the multi-carrier CV-QKD with QPSK, 256 QAM and Gaussian modulation can be evaluated with modulation noise εmod(k) simulated from the proposed modulation noise model. It is worth noting that the asymptotic SKR that is the maximum achievable SKR can be directly reflect the performance of the designed CV-QKD system.

4.1 QPSK modulation protocol

For the QPSK modulation multi-carrier CV-QKD scheme, the variance σ12 and σ22 are determined to be 1 because the random variables Ig and Qg (g = k, m, n, l) follow mutually independent binomial distribution (-1, 1). Under the total carrier number N = 10, 25, 40 and 60, the gain imbalance κk = 0.98 and the quadrature skew θk = π/50, the modulation noise ratio ${\varepsilon _{mod}}(k )/A_{sig}^2$ of the k-th subcarrier at different modulation index µk are calculated according to Eq. (16) and shown in Fig. 4. It is seen from Fig. 4 that ${\varepsilon _{mod}}(k )/A_{sig}^2$ of the k-th subcarrier is an increasing function of modulation index µk. Therefore, it is better to choose a lower modulation index µkto reduce the modulation noise especially for the larger total carrier number N in practical experiment, e.g. the modulation noise ratio for N = 60 is close to 0 when µk = 0.01.

 figure: Fig. 4.

Fig. 4. Simulated εmod(k)/A2sig of the k-th subcarrier at different modulation index µk with N = 10, 25, 40 and 60, respectively.

Download Full Size | PDF

With system parameters µk = 0.01, κk = 0.98 and θk = π/50, the worst modulation noise εmod(kworst) of the kworst-th subcarrier in different N carriers is numerically obtained as shown in Fig. 5. According to our previous QPSK single-carrier CV-QKD experimental results [18], the single-carrier excess noise εsingle can be controlled to be about 0.007 shot noise unit (SNU) with optimized modulation variance VA = 0.45 SNU. Suppose the k-th subcarrier share the same excess noise performance as the single-carrier CV-QKD scheme, except with the extra modulation noise induced in the multi-carrier quantum state prepration process. Then, the excess noise εmulti(kworst) of the kworst-th subcarrier in the QPSK multi-carrier CV-QKD scheme can be calculated based on Eq. (16), which is shown in Fig. 5. The excess noise εmulti(kworst) of the kworst-th subcarrier gradually increases with the total carrier number N, which should be set to below 60 for avoiding null key rate at 25 km transmission distance.

 figure: Fig. 5.

Fig. 5. Calculated noise of the kworst-th subcarrier at different total carrier number N in the QPSK multi-carrier CV-QKD. The blue dash dot line denotes the worst modulation noise εmod(kworst), the magenta solid line corresponds to the excess noise εmulti(kworst) and the red dash line represents the null key rate threshold at 25 km.

Download Full Size | PDF

Under the modulation variance VA = 0.45 SNU, the reconciliation efficiency βk = 95% (k = 1, 2, …, N), and the excess noise εmulti(kworst) calculated by Eq. (16), the single SKR ${R_{{k_{\textrm{worst}}}}}$ of the kworst-th subcarrier is evaluated by using the semidefinite programming (SDP) security analysis method [38,39], as shown in Fig. 6(a). In fact, the single SKR can be further improved by employing the security analysis method in [40] for QPSK CV-QKD protocol. Note that the modulation noise εmod(kworst) of the kworst-th subcarrier is the worst in N subcarriers. So, the total SKR of the N multi-carrier CV-QKD system with extra modulation noise can be obtained as $R \ge N{R_{{k_{\textrm{worst}}}}}$ according to Eq. (1), and the total SKR $N{R_{{k_{\textrm{worst}}}}}$ is shown in Fig. 6(b). From Fig. 6(b) that the total SKR of the multi-carrier CV-QKD can be greatly improved by increasing N. Moreover, the multi-carrier SKR gain is defined as $N{R_{{k_{\textrm{worst}}}}}/{R_{single}}$ and shown in Fig. 7, where Rsingle is the SKR of the single-carrier CV-QKD system without modulation noise. It is noteworthy that the multi-carrier SKR gain drastically drops with the increase of total carrier number N in long distance transmission due to the extra modulation noise. Besides, the optimal carrier number N can be numerically calculated based on the proposed modulation noise model as shown in Fig. 8, where the optimal N = 61, 56 and 43 for L = 5 km, 10 km and 25 km, respectively. In Fig. 8, the maximum N at different distances is obtained when SKR ${R_{{k_{\textrm{worst}}}}} = 0$, which reflects the null SKR threshold of the kworst-th subcarrier.

 figure: Fig. 6.

Fig. 6. The SKR of the QPSK multi-carrier CV-QKD with N = 10, 25, 40 and 60, respectively. (a) Single SKR of the kworst-th subcarrier, (b) total SKR of the QPSK multi-carrier CV-QKD.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. Calculated multi-carrier SKR gains for the QPSK multi-carrier CV-QKD with N = 10, 25, 40 and 60, respectively.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. Calculated QPSK multi-carrier SKR gains at different total carrier number N with transmission distance L = 5 km, 10 km, 25 km, respectively.

Download Full Size | PDF

4.2 256QAM modulation protocol

The SKR performance of the multi-carrier CV-QKD with 256QAM discrete modulation protocol is evaluated in a similar way based on the proposed modulation noise model with system parameters µk = 0.01, κk = 0.98 and θk = π/50. For 256QAM modulation, the random variables Ig and Qg (g = k, m, n, l) obey mutually independent discrete Maxwell-Boltzmann distribution, and the corresponding probability density function is given by [19]

$${\rho _{{x_i}}} = \frac{{\exp ({ - \nu x_i^2} )}}{{\sum\limits_{i = 1}^{16} {\exp ({ - \nu x_i^2} )} }}$$
where xi is normalized within (-1, 1). ν is a positive free parameter whose optimal value is about 0.04. Therefore, the variance σ12 of the quadrature Ig and the variance σ22 of the quadrature Ig2 can be obtained as, respectively
$$\sigma _1^2 = \sum\limits_{i = 1}^{16} {{{[{{x_i} - E({{x_i}} )} ]}^2}{\rho _{{x_i}}}}$$
$$\sigma _2^2 = \sum\limits_{i = 1}^{16} {{{[{x_i^2 - E({x_i^2} )} ]}^2}\rho _{{x_i}}^2}$$
with σ12 = 0.37 and σ22 = 0.11 applied in Eq. (16). So, the single SKR of the kworst-th subcarrier in 256QAM multi-carrier CV-QKD can be evaluated by using the improved SDP method in [39], as is shown in Fig. 9(a). Moreover, the modulation variance VA = 5 SNU, the reconciliation efficiency βk = 95%, the quantum efficiency η=0.56, the electronic noise υele = 0.15 SNU and the excess noise εsingle = 0.03 SNU are chosen in the following simulations referring to our previous 256QAM single-carrier CV-QKD experiment [19]. Note that the single-carrier excess noise εsingle = 0.03 SNU has excluded the IQ imbalance noise by DSP in the single-carrier CV-QKD experiment [19]. Furthermore, the corresponding total SKR and the multi-carrier SKR gain are respectively given in Fig. 9(b) and Fig. 10 for verifying the performance of 256QAM multi-carrier CV-QKD. Besides, the multi-carrier SKR gains with different N are calculated in Fig. 11, where the optimal N = 96, 85 and 72 for L = 25 km, 50 km and 100 km, respectively.

 figure: Fig. 9.

Fig. 9. The SKR of the 256QAM multi-carrier CV-QKD with N = 10, 40, 80 and 120, respectively. (a) Single SKR of the kworst-th subcarrier, (b) total SKR of the 256QAM multi-carrier CV-QKD.

Download Full Size | PDF

 figure: Fig. 10.

Fig. 10. Calculated multi-carrier SKR gains for the 256QAM multi-carrier CV-QKD with N = 10, 40, 80 and 120, respectively.

Download Full Size | PDF

 figure: Fig. 11.

Fig. 11. Calculated 256QAM multi-carrier SKR gains at different total carrier number N with transmission distance L = 25 km, 50 km and 100 km, respectively.

Download Full Size | PDF

4.3 Gaussian modulation protocol

We also estimate the SKR performance of the multi-carrier CV-QKD with Gaussian modulation protocol based on the proposed modulation noise model with system parameters µk = 0.01, κk = 0.98 and θk = π/50. The random variables Ig and Qg (g = k, m, n, l) are mutually independent and follow the normal Gaussian distribution $N\left( {0,\,\sigma _1^2 } \right)$ with amplitude normalization (-1, 1), whose with probability density function is given by

$$f(x )= \frac{1}{{\sqrt {2\pi } {\sigma _1}}}{e^{ - \frac{{{x^2}}}{{2\sigma _1^2}}}}$$
and the corresponding random variables Ig2 and Qg2 follow the Gaussian distribution N12, 2σ14), which can be verified as
$$E({{x^2}} )= \int_{ - \infty }^{ + \infty } {{x^2}} f(x )dx = \sigma _1^2$$
$$\sigma _2^2\textrm{ = }E({{x^4}} )- {[{E({{x^2}} )} ]^2} = \int_{ - \infty }^{ + \infty } {{x^4}} f(x )dx - {\sigma _1}^4 = 2{\sigma _1}^4$$

Therefore, the variance σ12 = 1/9 and σ22 = 2/81 can be obtained when the standard normal distribution N(0, σ2 = 1) is normalized to (-1, 1) from the amplitude section (-3σ, 3σ). Similarly, applying σ12 and σ22 in Eq. (16), the single SKR of the kworst-th subcarrier for the Gaussian modulated multi-carrier CV-QKD can be evaluated based on the no-switch protocol security analysis method [12], as shown in Fig. 12(a). In the SKR calculation, the modulation variance VA = 5 SNU, the quantum efficiency η=0.56, the electronic noise υele = 0.1 SNU, the reconciliation efficiency βk = 95% and the excess noise εsingle = 0.03 SNU are chosen as in [11]. Similarly, the corresponding total SKR and multi-carrier SKR gain are given in Fig. 12(b) and Fig. 13, respectively, which also verifies that the SKR can be greatly increased based on the multi-carrier CV-QKD scheme. Moreover, the multi-carrier SKR gains at different N are demonstrated Fig. 14, where the optimal N = 155, 130 and 110 for L = 50 km, 100 km and 150 km, respectively.

 figure: Fig. 12.

Fig. 12. The SKR of the Gaussian modulated multi-carrier CV-QKD with N = 10, 80, 120 and 160, respectively. (a) Single SKR of the kworst-th subcarrier, (b) total SKR of the Gaussian modulated multi-carrier CV-QKD.

Download Full Size | PDF

 figure: Fig. 13.

Fig. 13. Calculated multi-carrier SKR gains for Gaussiam modulated multi-carrier CV-QKD with N = 10, 80, 120 and 160, respectively

Download Full Size | PDF

 figure: Fig. 14.

Fig. 14. Calculated Gaussian modulated multi-carrier SKR gains at different total carrier number N with transmission distance L = 50 km, 100 km and 150 km, respectively.

Download Full Size | PDF

5. Conclusion

We have demonstrated a systematic modulation noise model for OFDM-based multi-carrier CV-QKD with arbitrary modulation protocol. Specifically, the extra modulation noise including the IQ imbalance noise and the third-order intermodulation noise in the multi-carrier quantum state preparation is theoretically modeled and analyzed, which limits the SKR performance of system with total carrier number N. Under considering the extra modulation noise, the SKR performance of the OFDM-based multi-carrier LLO-CV-QKD system with QPSK, 256QAM and Gaussian modulation protocol can be quantitatively evaluated by combining the security analysis method of the single-carrier CV-QKD system. Under practical system parameters, our simulation results show that the OFDM-based multi-carrier CV-QKD scheme can greatly improve the asymptotic SKR by increasing N even with imperfect modulations. Moreover, a maximum total SKR can be obtained by carefully optimizing N. Note that the IQ imbalance noise and the third-order intermodulation noise can be further reduced by optimizing the related parameters and designing the corresponding compensation algorithm at Bob’s site for achieving a better SKR performance in the practical multi-carrier CVQKD system. Compared with the single-carrier CV-QKD scheme, our work offers a promising way to significantly improve the performance of CV-QKD. In the future, other noises from the channel and receiver in practical experiment should be further considered to evaluate the performance of multi-carrier CV-QKD system, such as the fiber channel crosstalk noise between these subcarriers.

Funding

National Key Research and Development Program of China (2020YFA0309704); Sichuan Science and Technology Program (2021YJ0313, 2022YFG0330, 2022ZYD0118); National Natural Science Foundation of China (61901425, 62101516, 62171418, 62201530, U19A2076); Technology Innovation and Development Foundation of China Cyber Security (JSCX2021JC001); Chengdu Major Science and Technology Innovation Program (2021-YF08-00040-GX); Chengdu Key Research and Development Support Program (2021-YF05-02430-GX, 2021-YF09-00116-GX); Foundation of Science and Technology on Communication Security Laboratory (61421030402012111); Major Project of the Department of Science and Technology of Sichuan (2022ZDZX0009).

Disclosures

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

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. F. Grosshans, G. V. Assche, J. Wenger, R. Brouri, N. J. Cerf, and P. Grangier, “Quantum key distribution using Gaussian-modulated coherent states,” Nature 421(6920), 238–241 (2003). [CrossRef]  

2. S. Pirandola, “Composable security for continuous variable quantum key distribution: Trust levels and practical key rates in wired and wireless networks,” Phys. Rev. Res. 3(4), 043014 (2021). [CrossRef]  

3. F. H. Xu, X. F. Ma, Q. Zhang, H. K. Lo, and J. W. Pan, “Secure quantum key distribution with realistic devices,” Rev. Mod. Phys. 92(2), 025002 (2020). [CrossRef]  

4. M. Pittaluga, M. Minder, M. Lucamarini, M. Sanzarol, R. Woodward, M. J. Li, Z. L. Yuan, and A. J. Shields, “600 km repeater-like quantum communications with dual-band stabilization,” Nat. Photonics 15(7), 530–535 (2021). [CrossRef]  

5. H. Guo, Z. Y. Li, S. Yu, and Y. C. Zhang, “Toward practical quantum key distribution using telecom components,” Fundam. Res. 1(1), 96–98 (2021). [CrossRef]  

6. T. Matsuura, K. Maeda, T. Sasaki, and M. Koashi, “Finite-size security of continuous-variable quantum key distribution with digital signal processing,” Nat. Commun. 12(1), 252 (2021). [CrossRef]  

7. G. Zhang, J. Y. Haw, H. Cai, F. Xu, S. M. Assad, J. F. Fitzsimons, X. Zhou, Y. Zhang, S. Yu, J. Wu, W. Ser, L. C. Kwek, and A. Q. Liu, “An integrated silicon photonic chip platform for continuous-variable quantum key distribution,” Nat. Photonics 13(12), 839–842 (2019). [CrossRef]  

8. N. Jain, H. Chin, H. Mani, C. Lupo, D. S. Nikolic, A. Kordts, S. Pirandola, T. B. Pedersen, M. Kolb, B. Ömer, C. Pacher, T. Gehring, and U. L. Andersen, “Practical continuous-variable quantum key distribution with composable security,” Nat. Commun. 13(1), 4740 (2022). [CrossRef]  

9. Y. Tian, P. Wang, J. Q. Liu, S. N. Du, W. Y. Liu, Z. G. Lu, X. Y. Wang, and Y. M. Li, “Experimental demonstration of continuous-variable measurement-device-independent quantum key distribution over optical fiber,” Optica 9(5), 492–500 (2022). [CrossRef]  

10. S. J. Ren, S. Yang, A. Wonfor, I. White, and R. Penty, “Demonstration of high-speed and low-complexity continuous variable quantum key distribution system with local local oscillator,” Sci. Rep. 11(1), 9454 (2021). [CrossRef]  

11. H. Wang, Y. D. Pi, W. Huang, Y. Li, Y. Shao, J. Yang, J. L. Liu, C. L. Zhang, Y. C. Zhang, and B. J. Xu, “High-speed Gaussian-modulated continuous-variable quantum key distribution with a local local oscillator based on pilot-tone-assisted phase compensation,” Opt. Express 28(22), 32882–32893 (2020). [CrossRef]  

12. T. Wang, P. Huang, Y. M. Zhou, W. Q. Liu, and G. H. Zeng, “Pilot-multiplexed continuous-variable quantum key distribution with a real local oscillator,” Phys. Rev. A 97(1), 012310 (2018). [CrossRef]  

13. S. Kleis, M. Rueckmann, and C. G. Schaeffer, “Continuous variable quantum key distribution with a real local oscillator using simultaneous pilot signals,” Opt. Lett. 42(8), 1588–1591 (2017). [CrossRef]  

14. F. Laudenbach, B. Schrenk, C. Pacher, M. Hentschel, C. F. Fung, F. Karinou, A. Poppe, M. Peev, and H. Hübel, “Pilot-assisted intradyne reception for high-speed continuous-variable quantum key distribution with true local oscillator,” Quantum 3, 193193 (2019). [CrossRef]  

15. M. Li and M. Cvijetic, “Continuous-variable quantum key distribution with self-reference detection and discrete modulation,” IEEE J. Quantum Electron. 54(5), 1–8 (2018). [CrossRef]  

16. Z. Qu, I. B. Djordjevic, and M. A. Neifeld, “RF-subcarrier-assisted four-state continuous-variable QKD based on coherent detection,” Opt. Lett. 41(23), 5507–5510 (2016). [CrossRef]  

17. Y. D. Pi, H. Wang, Y. Pan, Y. Shao, Y. Li, Y. C. Zhang, W. Huang, and B. J. Xu, “Sub-Mbps key-rate continuous-variable quantum key distribution with local-local-oscillator over 100 km fiber,” arXivarXiv:2212.11534 (2022).

18. H. Wang, Y. Li, Y. D. Pi, Y. Pan, Y. Shao, L. Ma, Y. C. Zhang, J. Yang, T. Zhang, W. Huang, and B. J. Xu, “Sub-Gbps key rate four-state continuous-variable quantum key distribution within metropolitan area,” Commun. Phys. 5(1), 162 (2022). [CrossRef]  

19. Y. Pan, H. Wang, Y. Shao, Y. D. Pi, Y. Li, B. Liu, W. Huang, and B. J. Xu, “Experimental demonstration of high-rate discrete-modulated continuous-variable quantum key distribution system,” Opt. Lett. 47(13), 3307–3310 (2022). [CrossRef]  

20. F. Roumestan, A. Ghazisaeidi, J. Renaudier, T. V. Vidarte, E. Diamanti, and P. Grangier, “High-Rate Continuous variable quantum key distribution based on probabilistically shaped 64 and 256-QAM,” European Conference on Optical Communication, Th2G3, 1–4 (2021).

21. X. X. Zhang, Y. C. Zhang, Z. Y. Li, S. Yu, and H. Guo, “1.2 GHz balanced homodyne detector for continuous-variable quantum information technology,” IEEE Photonics J. 10(5), 1–10 (2018). [CrossRef]  

22. Y. Li, X. F. Zhang, Y. Li, B. J. Xu, L. Ma, J. Yang, and W. Huang, “High-throughput GPU layered decoder of multi-edge type low density parity check codes in continuous-variable quantum key distribution systems,” Sci. Rep. 10(1), 14561 (2020). [CrossRef]  

23. S. S. Yang, Z. G. Lu, and Y. M. Li, “High-speed post-processing in continuous-variable quantum key distribution based on FPGA implementation,” J. Lightwave Technol. 38(15), 3935–3941 (2020). [CrossRef]  

24. C. Zhou, Y. Li, L. Ma, J. Yang, W. Huang, H. Wang, Y. J. Luo, F. C. M. Lau, Y. Li, and B. J. Xu, “High-throughput decoder of quasi-cyclic LDPC codes with limited precision for continuous-variable quantum key distribution systems,” arXivarXiv:2207.01860 (2022).

25. Y. C. Zhang, Z. Y. Chen, S. Pirandola, X. Y. Wang, C. Zhou, B. J. Chu, Y. J. Zhao, B. J. Xu, S. Yu, and H. Guo, “Long-distance continuous-variable quantum key distribution over 202.81 km of fiber,” Phys. Rev. Lett. 125(1), 010502 (2020). [CrossRef]  

26. J. Lodewyck, T. Debuisschert, R. Tualle-Brouri, and P. Grangier, “Controlling excess noise in fiber optics continuous variables quantum key distribution,” Phys. Rev. A 72(5), 050303 (2005). [CrossRef]  

27. T. A. Eriksson, R. S. Luís, B. J. Puttnam, G. Rademacher, M. Fujiwara, Y. Awaji, H. Furukawa, N. Wada, M. Takeoka, and M. Sasaki, “Wavelength division multiplexing of 194 continuous variable quantum key distribution channels,” J. Lightwave Technol. 38(8), 2214–2218 (2020). [CrossRef]  

28. L. Gyongyosi, “Multi-carrier continuous-variable quantum key distribution,” Theor. Comput. Sci. 816, 67–95 (2020). [CrossRef]  

29. Y. J. Wang, Y. Y. Mao, W. T. Huang, D. Huang, and Y. Guo, “Optical frequency comb-based multichannel parallel continuous-variable quantum key distribution,” Opt. Express 27(18), 25314–25329 (2019). [CrossRef]  

30. J. Fang, P. Huang, and G. H. Zeng, “Multi-channel parallel continuous variable quantum key distribution with Gaussian modulation,” Phys. Rev. A 89(2), 022315 (2014). [CrossRef]  

31. H. Zhang, Y. Mao, D. Huang, J. W. Li, L. Zhang, and Y. Guo, “Security analysis of orthogonal-frequency-division-multiplexing-based continuous-variable quantum key distribution with imperfect modulation,” Phys. Rev. A 97(5), 052328 (2018). [CrossRef]  

32. S. Bahrani, M. Razavi, and J. A. Salehi, “Orthogonal frequency-division multiplexed quantum key distribution,” J. Lightwave Technol. 33(23), 4687–4698 (2015). [CrossRef]  

33. W. Zhao, Q. Liao, D. Huang, and Y. Guo, “Performance analysis of the satellite-to-ground continuous-variable quantum key distribution with orthogonal frequency division multiplexed modulation,” Quantum Inf. Process 18(1), 39 (2019). [CrossRef]  

34. C. Roberto, “Phase-noise limitations in continuous-variable quantum key distribution with homodyne detection,” Phys. Rev. A 95(2), 022315 (2017). [CrossRef]  

35. Y. C. Zhang, Y. D. Huang, Z. Y. Chen, Z. Y. Li, S. Yu, and H. Guo, “One-time shot-noise unit calibration method for continuous-variable quantum key distribution,” Phys. Rev. Appl. 13(2), 024058 (2020). [CrossRef]  

36. H. Wang, S. J. Zhang, X. H. Zou, Y. L. Zhang, S. Liu, and Y. Liu, “Bias-Independent and self-calibrated electrical method for microwave characterization of dual-parallel Mach-Zehnder modulators based on two-tone and bias-swing modulation,” IEEE Trans. Microwave Theory Techn. 65(7), 2636–2643 (2017). [CrossRef]  

37. F. Laudenbach, C. Pacher, C. F. Fung, A. Poppe, M. Peev, B. Schrenk, M. Hentschel, P. Walther, and H. Hübel, “Continuous-variable quantum key distribution with Gaussian modulation-the theory of practical implementations,” Adv. Quantum Technol. 1(1), 1800011 (2018). [CrossRef]  

38. S. Ghorai, P. Grangier, E. Diamanti, and A. Leverrier, “Asymptotic security of continuous-variable quantum key distribution with a discrete modulation,” Phys. Rev. X 9(2), 021059 (2019). [CrossRef]  

39. A. Denys, P. Brown, and A. Leverrier, “Explicit asymptotic secret key rate of continuous-variable quantum key distribution with an arbitrary modulation,” Quantum 5, 540540 (2021). [CrossRef]  

40. J. Lin, T. Upadhyaya, and N. Lütkenhaus, “Asymptotic security analysis of discrete-modulated continuous-variable quantum key distribution,” Phys. Rev. X 9(4), 041064 (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 (14)

Fig. 1.
Fig. 1. Schematic diagram of the OFDM-based multi-carrier CV-QKD with local local oscillator. LD: laser diode; OFDM: orthogonal frequency division multiplexing; MZM: Mach-Zehnder modulator; PBS: polarization beam splitter; VOA, variable optical attenuator; PBC: polarization beam combiner; SMF: single mode fiber; PSA: polarization synthesis analyzer; OC: optical coupler; BPD: balanced photo-detector; DSP: digital signal processing
Fig. 2.
Fig. 2. The data processing of the OFDM generator. S/P: serial-to-parallel conversion; TS: training sequence; IFFT: inverse fast Fourier transform; P/S: parallel-to-serial conversion; CP: cyclic prefix; DAC: digital-to-analog converter.
Fig. 3.
Fig. 3. The DSP routine of the OFDM-based multi-carrier LLO-CV-QKD scheme. CP: cyclic prefix; S/P: serial-to- parallel conversion; FFT: fast Fourier transform.
Fig. 4.
Fig. 4. Simulated εmod(k)/A2sig of the k-th subcarrier at different modulation index µk with N = 10, 25, 40 and 60, respectively.
Fig. 5.
Fig. 5. Calculated noise of the kworst-th subcarrier at different total carrier number N in the QPSK multi-carrier CV-QKD. The blue dash dot line denotes the worst modulation noise εmod(kworst), the magenta solid line corresponds to the excess noise εmulti(kworst) and the red dash line represents the null key rate threshold at 25 km.
Fig. 6.
Fig. 6. The SKR of the QPSK multi-carrier CV-QKD with N = 10, 25, 40 and 60, respectively. (a) Single SKR of the kworst-th subcarrier, (b) total SKR of the QPSK multi-carrier CV-QKD.
Fig. 7.
Fig. 7. Calculated multi-carrier SKR gains for the QPSK multi-carrier CV-QKD with N = 10, 25, 40 and 60, respectively.
Fig. 8.
Fig. 8. Calculated QPSK multi-carrier SKR gains at different total carrier number N with transmission distance L = 5 km, 10 km, 25 km, respectively.
Fig. 9.
Fig. 9. The SKR of the 256QAM multi-carrier CV-QKD with N = 10, 40, 80 and 120, respectively. (a) Single SKR of the kworst-th subcarrier, (b) total SKR of the 256QAM multi-carrier CV-QKD.
Fig. 10.
Fig. 10. Calculated multi-carrier SKR gains for the 256QAM multi-carrier CV-QKD with N = 10, 40, 80 and 120, respectively.
Fig. 11.
Fig. 11. Calculated 256QAM multi-carrier SKR gains at different total carrier number N with transmission distance L = 25 km, 50 km and 100 km, respectively.
Fig. 12.
Fig. 12. The SKR of the Gaussian modulated multi-carrier CV-QKD with N = 10, 80, 120 and 160, respectively. (a) Single SKR of the kworst-th subcarrier, (b) total SKR of the Gaussian modulated multi-carrier CV-QKD.
Fig. 13.
Fig. 13. Calculated multi-carrier SKR gains for Gaussiam modulated multi-carrier CV-QKD with N = 10, 80, 120 and 160, respectively
Fig. 14.
Fig. 14. Calculated Gaussian modulated multi-carrier SKR gains at different total carrier number N with transmission distance L = 50 km, 100 km and 150 km, respectively.

Equations (29)

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

R = k = 1 N R k , R k = β k I ( A k : B k ) χ ( B k : E )
I s = k = 1 N I k cos ( 2 π f k t ) k = 1 N Q k sin ( 2 π f k t )
Q s = k = 1 N I k sin ( 2 π f k t ) + k = 1 N Q k cos ( 2 π f k t )
| X s i g + j P s i g = k = 1 N | X k + j P k + Δ X k + j Δ P k
E s i g ( t ) = 2 A s i g e j 2 π f A t { G 1 sin [ k = 1 N γ k cos ( 2 π f k t + φ k ) ] + j G 2 sin [ k = 1 N γ k sin ( 2 π f k t + φ k ) ] }
γ k = μ k I k 2 + Q k 2
cos φ k = I k I k 2 + Q k 2
E s i g ( t ) = 2 A s i g e j 2 π f A t k = 1 N γ k e j ( 2 π f k t + φ k ) = 2 A s i g e j 2 π f A t + j 2 π f k t k = 1 N | I k + j Q k
E s i g ( t ) = 2 A s i g e j 2 π f A t { G 1 p 1 = + J p 1 ( γ 1 ) p k = + J p k ( γ k ) p N = + ( γ N ) J p N sin [ ( 2 π f 1 t p 1 + p 1 φ 1 + p 1 π 2 ) + + ( 2 π f k t p k + p k φ k + p k π 2 ) + + ( 2 π f N t p N + p N φ N + p N π 2 ) ] + j G 2 p 1 = + J p 1 ( γ 1 ) p k = + J p k ( γ 2 ) p N = + J p N ( γ N ) sin [ ( 2 π f 1 t p 1 + p 1 φ 1 ) + + ( 2 π f k t p k + p k φ k ) + + ( 2 π f N t p N + p N φ N ) ] }
E s i g ( k ) ( t ) = 2 A s i g e j π f A t { J 1 ( γ k ) [ ( G 1 + G 2 ) e j ( 2 π f k t + φ k ) + ( G 1 G 2 ) e j ( 2 π f k t + φ k ) ] , n = 1 n k 2 m + n = k N J 2 ( γ m ) J 1 ( γ n ) [ ( G 1 + G 2 ) e j ( 2 π f k t + 2 φ m + φ n ) + ( G 1 G 2 ) e j ( 2 π f k t + 2 φ m + φ n ) ] , s = 1 n k 2 m n = k N J 2 ( γ m ) J 1 ( γ n ) [ ( G 1 + G 2 ) e j ( 2 π f k t + 2 φ m φ n ) + ( G 1 G 2 ) e j ( 2 π f k t + 2 φ m φ n ) ] , s , l = 1 n l k + n + l = k N J 1 ( γ m ) J 1 ( γ n ) J 1 ( γ l ) [ ( G 1 + G 2 ) e j ( 2 π f k t + φ m + φ n + φ l ) + ( G 1 G 2 ) e j ( 2 π f k t + φ m + φ n + φ l ) ] , n , l = 1 n l k + n l = k N J 1 ( γ m ) J 1 ( γ n ) J 1 ( γ l ) [ ( G 1 + G 2 ) e j ( 2 π f k t + φ m + φ n φ l ) + ( G 1 G 2 ) e j ( 2 π f k t + φ m + φ n φ l ) ] , n , l = 1 n l k n l = k N J 1 ( γ m ) J 1 ( γ n ) J 1 ( γ l ) [ ( G 1 + G 2 ) e j ( 2 π f k t + φ m φ n φ l ) + ( G 1 G 2 ) e j ( 2 π f k t + φ m φ n φ l ) ] }
X s i g ( k ) = X k + Δ X k = 2 A s i g real ( γ k e j φ k ) + 2 A s i g r e a l { γ k 2 [ ( G 1 + G 2 2 ) e j φ k + ( G 1 G 2 ) e j φ k ] M 1 ( N , k ) γ m 2 γ n 8 [ ( G 1 + G 2 ) e j ( 2 φ m + φ n ) + ( G 1 G 2 ) e j ( 2 φ m + φ n ) ] M 2 ( N , k ) γ m 2 γ n 8 [ ( G 1 + G 2 ) e j ( 2 φ m φ n ) + ( G 1 G 2 ) e j ( 2 φ m φ n ) ] W 1 ( N , k ) γ m γ n γ l 8 [ ( G 1 + G 2 ) e j ( φ m + φ n + φ l ) + ( G 1 G 2 ) e j ( φ m + φ n + φ l ) ] W 2 ( N , k ) γ m γ n γ l 8 [ ( G 1 + G 2 ) e j ( φ m + φ n φ l ) + ( G 1 G 2 ) e j ( φ m + φ n φ l ) ] W 3 ( N , k ) γ m γ n γ l 8 [ ( G 1 + G 2 ) e j ( φ m φ n φ l ) + ( G 1 G 2 ) e j ( φ m φ n φ l ) ] }
X k = 2 A s i g μ k I k
Δ X k = Δ X k 1 + Δ X k 2 + Δ X k 3
Δ X k 1 = A s i g μ k [ ( κ k cos θ k 1 ) I k + ( κ k sin θ k ) Q k ]
Δ X k 2 = 1 4 A s i g μ m 2 μ n { ( 1 + κ k cos θ k ) { [ M 1 ( N , k ) + M 2 ( N , k ) ] Q m 2 I n + 2 [ M 1 ( N , k ) M 2 ( N , k ) ] I m Q m Q n [ M 1 ( N , k ) + M 2 ( N , k ) ] I m 2 I n } + ( κ k sin θ k ) [ [ M 1 ( N , k ) + M 2 ( N , k ) ] I m 2 Q n + 2 [ M 1 ( N , k ) M 2 ( N , k ) ] I m Q m I n [ M 1 ( N , k ) + M 2 ( N , k ) ] Q m 2 Q n ] }
Δ X k 3 = 1 4 A s i g μ m μ n μ l { ( 1 + κ k cos θ k ) { [ W 1 ( N , k ) W 2 ( N , k ) + W 3 ( N , k ) ] I m Q n Q l + [ W 1 ( N , k ) W 2 ( N , k ) W 3 ( N , k ) ] Q m I n Q l + [ W 1 ( N , k ) + W 2 ( N , k ) W 3 ( N , k ) ] Q m Q n I l [ W 1 ( N , k ) + W 2 ( N , k ) + W 3 ( N , k ) ] I m I n I l } + ( κ k sin θ k ) { [ W 1 ( N , k ) + W 2 ( N , k ) W 3 ( N , k ) ] I m I n Q l + [ W 1 ( N , k ) W 2 ( N , k ) W 3 ( N , k ) ] I m Q n I l + [ W 1 ( N , k ) W 2 ( N , k ) + W 3 ( N , k ) ] Q m I n I l [ W 1 ( N , k ) + W 2 ( N , k ) + W 3 ( N , k ) ] Q m Q n Q l } }
Δ X k 2 = Δ X k 1 2 + Δ X k 2 2 + Δ X k 3 2
Δ X k 1 2 = A s i g 2 μ k 2 σ 1 2 ( κ k 2 + 1 2 κ k cos θ k )
Δ X k 2 2 = A s i g 2 μ m 4 μ n 2 8 ( 1 + 2 κ k cos θ k + κ k 2 ) { [ M 1 ( N , k ) + M 2 ( N , k ) ] 2 σ 2 2 σ 1 2 + 2 [ M 1 ( N , k ) M 2 ( N , k ) ] 2 σ 1 6 }
Δ X k 3 2 = A s i g 2 μ m 2 μ n 2 μ l 2 4 ( 1 + 2 κ k cos θ k + κ k 2 ) [ W 1 2 ( N , k ) + W 2 2 ( N , k ) + W 3 2 ( N , k ) ] σ 1 6
V A = X k 2 = 4 A s i g 2 μ k 2 I k 2 = 4 A s i g 2 μ k 2 σ 1 2
ε m o d ( k ) = Δ X k 2 = V A 4 ( κ k 2 + 1 2 κ k cos θ k ) + V A μ k 4 32 ( 1 + 2 κ k cos θ k + κ k 2 ) { [ M 1 ( N , k ) + M 2 ( N , k ) ] 2 σ 2 2 + 2 { [ M 1 ( N , k ) M 2 ( N , k ) ] 2  +  W 1 2 ( N , k ) + W 2 2 ( N , k ) + W 3 2 ( N , k ) } σ 1 4 }
ε m u l t i ( k ) = ε m o d ( k ) + ε s i n g l e
ρ x i = exp ( ν x i 2 ) i = 1 16 exp ( ν x i 2 )
σ 1 2 = i = 1 16 [ x i E ( x i ) ] 2 ρ x i
σ 2 2 = i = 1 16 [ x i 2 E ( x i 2 ) ] 2 ρ x i 2
f ( x ) = 1 2 π σ 1 e x 2 2 σ 1 2
E ( x 2 ) = + x 2 f ( x ) d x = σ 1 2
σ 2 2  =  E ( x 4 ) [ E ( x 2 ) ] 2 = + x 4 f ( x ) d x σ 1 4 = 2 σ 1 4
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.