Atmospheric continuous-variable quantum communication

We present a quantum communication experiment conducted over a point-to-point free-space link of 1.6 km in urban conditions. We study atmospheric influences on the capability of the link to act as a continuous-variable (CV) quantum channel. Continuous polarization states (that contain the signal encoding as well as a local oscillator in the same spatial mode) are prepared and sent over the link in a polarization multiplexed setting. Both signal and local oscillator undergo the same atmospheric fluctuations. These are intrinsically auto-compensated which removes detrimental influences on the interferometric visibility. At the receiver, we measure the Q-function and interpret the data using the framework of effective entanglement. We compare different state amplitudes and alphabets (two-state and four-state) and determine their optimal working points with respect to the distributed effective entanglement. Based on the high entanglement transmission rates achieved, our system indicates the high potential of atmospheric links in the field of CV QKD.


Introduction
Quantum communication refers to the distribution of quantum states between two parties via a quantum channel. With regard to this, it is crucial that this quantum channel preserves the quantum properties of the distributed states. The most common channel implementation are optical fibers and the free space. The latter offers great flexibility in terms of infrastructure establishment and links to moving objects are also feasible, see e.g. [1]. For a review of representative free-space quantum communication experiments see [2].
Quantum key distribution (QKD) [3,4] is probably the most practical branch of quantum communication and concerns the establishment of a secret key jointly between two legitimate parties, Alice and Bob. As the security is based on the laws of quantum mechanics, in principle information theoretic-security can be achieved [4].
Free-space QKD over a real atmospheric channel was first demonstrated in 1996 [5]. Since then, a number of prepare-and-measure as well as entanglement based schemes have been implemented in free space; the longest distance so far achieved on Earth is 144 km [6,7]. Nowadays, even quantum communication between earth and space is being conceived [8,9,10,11,12,13,14]. All of the systems so far referred to here have one major aspect in common: they are based on discrete quantum variables and use single photon threshold ("click") detectors, which involves spatial, spectral and/or temporal filtering in order to reduce background noise (see e.g. [15,16] for the first point-to-point demonstrations of free-space QKD in daylight).
As already shown by Bennett [17], any two non-orthogonal quantum states suffice to ensure secure key distribution. This paved the way for continuous-variable protocols (for a review see [18,19]), based on a different approach: performing homodyne measurements on weak coherent states with the help of a bright local oscillator (LO). Generally, homodyne detectors offer immunity to stray light and the PIN photodiodes used therin own higher quantum efficiencies than avalanche photodiodes. Initially, a discrete modulation of Gaussian states was proposed [20,21,22], followed by Gaussian modulation [23,24] shortly thereafter.
The homodyning technique is well established in classical optical communication including satellite based links [25]. There, the signal states are at high intensities, their overlap is negligible, and generating a LO locally at the receiver is appropriate. CV QKD uses strongly overlapping and thus almost indistinguishable states, which requires high detection efficiencies as well as a high interferometric visibility. We fulfilled this demand by means of a specially developed protocol using the polarization degree of freedom to multiplex signal and LO in the same spatial channel mode already at the sender [26,27]. A successful pilot project proved the feasibility of CV QKD via a 100 m free-space link on the roof of our institute [28,29]. We now succeeded to established a point-to-point free-space link connecting two buildings at a distance of 1.6 km in an urban environment within the city of Erlangen. A two-dimensional modulation of the signal states and their simultaneous detection along two conjugate quadratures (as initially shown in [26] and [30]) was implemented, using the aforementioned polarization multiplexing. For the moment, we focus on the discrete modulation of two or four signal states, but in principle, a variety of alphabets such as Gaussian modulation [23,24] or ring-type alphabets [31] can be realized with our system. Note that due to the inherent auto-compensation of phase fluctuations, a freely drifting phase, such as might be the case with fiber channels [32], can be avoided. Effective entanglement (EE) [33,34] is a powerful measure in describing the quantum correlations and their preservation during channel transmission and a necessary precondition for QKD. We verify and quantify EE in terms of minimal negativity [36,37] according to the method introduced in [33,34,35,36]. In our free-space scenario, special consideration has to be given to atmospheric turbulence that not only affects the signal states but also leads to fluctuations of the LO and to undesired classical excess noise in general. Such effects on quantum continuous-variable states in a turbulent atmosphere have been studied recently from a theoretical [38,39,40] as well as an experimental [28,29,41] point of view. Furthermore, there have also been ideas suggesting to exploit the fading channel properties in terms of improved quantum state propagation [42,43,44,45,38,39,46,47].

Stokes operators
In quantum mechanics, polarization can be described by the quantum Stokes operators, which are the quantum counterparts of the classical Stokes parameters [48]. The Stokes operators are introduced and defined, for example, in [49] and read aŝ in terms of the creation and annihilation operatorsâ † andâ of the respective polarization modes. The arrows in brackets display the operational definitions of the Stokes parameters as intensity differences of polarization types. In our experiment, where the LO is circularly polarized and the signal states are measured by a homodyne detection of the S 1/2 -Stokes operators (see Figure 1), this results in the uncertainty relation: In the case of coherent states, the equality holds and the variances ofŜ 1 andŜ 2 are equal. Polarization excess noise results in increased variances. The co-propagation of signal and LO ensures a perfect spatial interference between the two, enabling high detection efficiency without any additional interference stabilization, as channel induced phase fluctuations are auto-compensated. For our free-space system, the co-propagation has further advantageous side effects: the LO acts as a spatial and spectral filter, such that only those photons that are mode-matched with it will result in a significant detector signal. In contrast to what is the case in single photon experiments, here there is no need for spatial or spectral filtering, and background light (that is not local oscillator (S3 ) + signal (S )

Experimental Setup
The optical setup of our QKD system is shown in Figure 1. At the sender, a grating stabilized continuous wave diode laser (Toptica DLL 100) at a wavelength of 809 nm is spatially mode cleaned and guided to an alignable breadboard (see upper left part in Figure 1) via a single mode polarization maintaining fiber. A small portion is then split off to enable monitoring of the shot-noise limitation of the laser output states in a selfhomodyning setup. With the remaining portion, we prepare the circularly polarized LO using a sequence of two polarizing beam splitters (PBS), a half wave plate (HWP) and a quarter wave plate (QWP). Two electro-optical modulators (Thorlabs, EO-AM-NR-C1, 600-900 nm, bandwidth 100 MHz), that sandwich a HWP for 45 • polarization rotation, are voltage driven by two synchronized arbitrary waveform generators (Agilent 33250A) and serve for the generation of signal states using the Pockels effect. At a repetition rate of 3.05 MHz, Gaussian shaped voltage pulses in the order of only few mV vary the crystal's birefringence. They cause a time dependent phase shift, diminishing the amount of circular (LO) polarization and introducing weak signal components. The signal energy is taken from the LO, which, due to the weak modulation, remains essentially unaffected. Each signal pulse is followed by a zero voltage period of the same length, acting as the vacuum reference. After 263 signal pulses, the modulation voltage is increased for exactly one pulse that later serves as the trigger for a clock recovery in the receiving process. The state generation process is verified by a double homodyne detection scheme realized as a simultaneous measurement of the S 1/2 -Stokes parameters. The two outputs of a non polarizing, symmetric beam splitter are equipped with a combination of a quarter and a half wave plate (to compensate for polarization offsets and to adjust the S 1/2 basis), a polarizing beam splitter and a homemade detector, respectively. The required linearity of the detection system is verified by an attenuation measurement of the unmodulated LO. The detectors measure the electronic difference current of two photo diodes (the overall detection efficiency including optical losses and the diodes' quantum efficiencies is 0.84 ± 0.02.) and amplify the translated voltage. The signal obtained is then electronically high-pass filtered (Minicircuits BLK-89-S+, 100 kHz) ‡ to compensate for long term drifts, and measured by an oscilloscope at a sampling rate of 250 MHz. One time slot, containing signal and vacuum reference, consists of 82 samples; the signal state -or corresponding vacuum -is determined by integrating over 41 samples. To prevent any influence of the trigger pulse on the quantum signals, we disregard the first 92 out of 264 slots, leading to an effective state modulation rate of 2.22 MHz. After the state preparation, the signal and LO beam is spatially expanded to 4 cm in diameter and sent to the receiver through the 1.6 km free-space link. Atmospheric turbulence causes the beam to diverge more than expected by the diffraction limit, leading to an increased beam width after channel propagation. This effect can be attributed to both, wavefront distortions and "beam wandering" or spatial jitter. The latter is visualized in the inset in Figure 2 showing three examples of shots of the spatial beam profile arriving at Bob. At the receiver, a telescope with a front aperture of 150 mm is needed to capture most of the arriving beam intensity. After reduction of the beam diameter and splitting off a small part (< 4%) for transmission and LO monitoring, the signal states are detected in the same way as described above for Alice (see Figure 1). The channel induces intensity fluctuations, Figure 2 shows the statistical distribution of the channel transmission T , monitored in a direct detection simultaneously with the quantum signal. The width of one of the 35 bins (∆T ) is 0.9%, T mean is ≈ 76.1%. The triangles indicate those bins that had sufficiently large sample sizes and were thus used for further processing. They contain 92 % of events. These measurements were made on a rather clear and calm night (September 4 th 2013). We want to point out that daylight operation of our homodyne system is easily feasible without any spectral or spatial filters. However, we chose to measure at night, as atmospheric turbulence is increased during the day, which would directly lead to a broader transmission distribution at a lower mean value.
The linearity of Bob's detection setup was equally well verified by an attenuation ‡ Note that we verified that the influence of this filter on the signal pulses at a repetition rate of 3.05 MHz is negligible. measurement of the unmodulated LO at intensities down to zero. Here, the overall detection efficiency (not including the monitoring split-off) is 0.83 ± 0.02. For a typical received LO power of 14 mW, the electronic noise is more than 10 dB below the shotnoise level. The fading channel, however, has further impact on Bob's detection. Most important is the fluctuation of the shot-noise level, that is given by the unmodulated LO. Hence the signal modulation scheme is performed such that each signal slot is followed directly by its corresponding vacuum reference. We apply again a high-pass filter to avoid slow drifts in the offset of the signal, as the dynamic range of the oscilloscope's A/D conversion is limited.

Results
We use the classical transmission value to sort the collected data according to the prevailing transmission bins (see [41,39]). This technique allows us to directly compare the measured Stokes parameter variances to the corresponding shot-noise variance. For an alphabet of four states distributed symmetrically (|+α , |+iα , |−α , |−iα ), the variances of 15000 states and corresponding vacua are calculated for each individual transmission bin. We discard the bins that contain too few samples; the remaining data have small statistical errors. The results are shown in Figure 3, averaged over four states and individually for the S 1 /S 2 direction. The error is estimated by one standard deviation, the expected linearity of the vacuum states is shown by the linear Only those bins containing statistically sufficient numbers of samples (as indicated by the triangles in Figure 2) are considered. We averaged over four states, measured simultaneously along the S 1 and S 2 direction. The error is estimated conservatively by one standard deviation. The measurements of the vacuum are consistent with a linear fit, as theoretically expected.
fit. The (mean) excess noise of 0.01 shot-noise units (SNU) of the signal states compared to vacuum is mainly of technical origin, and is already introduced in Alice's signal generation step. This is due to the stronger fluctuating voltage pulses applied to the EOMs compared to the unmodulated vacuum slots. Note that the signal amplitudes at Alice have been chosen near 1, as this was determined as an optimal value for the four state alphabet in terms of negativity guaranteeing the maximal possible rate of effective entanglement. This will be explained in more detail in the next section.
For the largest occupied transmission bin, the Q-function [50] of four signal states received at Bob was measured; each signal state dataset contained 270000 samples. Figure 4 shows on the left hand side each state plotted individually as well as the combined Q-function of the mixed state (right). The height of the Q-function is an indicator of mixedness of the depicted state [27]. The maximal value of 1 π can be achieve by a pure state, as e.g. the vacuum state (not shown here). The individual states (left) are almost pure, the peak height of the mixed state is distinctly less (≈ 0.14), as theoretically expected.
The signal amplitudes, which have been attenuated by the channel transmission of ≈ 81.2%, are: 0.88, -0.87 (S 1 observable) and 0.92, -0.92 (S 2 observable). The slight asymmetry is already introduced at Alice, mainly by the fact that the same voltages are applied to both EOMs, the transmission of which is ≈ 95%. As the first EOM modulates along the S 1 axis, these signals are slightly attenuated.
Effective entanglement verification and quantification The use of entanglement verification describing correlations of a bipartite state shared between Alice and Bob as a first necessary precondition has been proposed in the context of QKD protocols [51,33]. This is especially useful in the absence of a full security proof considering realistic amounts of excess noise as well as finite size effects for discrete modulation CV QKD protocols. (Still, promising advances in this direction have been made within the last years [52,53,54]). The method was extended to the Stokes operators ( [34], experimental realization in [27,45]), where monitoring of the LO can help identify the impact of an adversary. Following [36,37], we use the negativity [55,56,57] as an entanglement measure that can be minimized over all bipartite states consistent with the information available. By means of semi-definite programming [58], we find optimal working points with respect to the initial state overlap and measured excess noise for our continuous-variable system. Formally, this treatment applies to quadrature operators, not Stokes operators. As in [37], we obtain quadrature measurements from Stokes measurements by considering a specific trusted form for the local oscillator. As shown in Figure 5, we compare two alphabets using discrete modulation of two versus four signal states. The four-state alphabet clearly outperforms the two-state counterpart in terms of acceptable excess noise and higher negativity bounds. The optimal amplitudes also depend greatly on the choice of alphabets. For these plots, the transmission was assumed to be fixed at 63 % (referring to its mean value, including detector losses).
Unlike the fixed channel transmission of the experimental fiber based setup in [37], we have to cope with the naturally fluctuating transmission. For the four-state alphabet, we now perform the minimization for each of the transmission sub-channels individually, at a fixed and optimal amplitude near 1. The result is shown in Figure 6. At a channel transmission close to 100% (or 83% with detection efficiency considered), obviously, the entanglement transfer would be optimal. For the estimation of an overall entanglement transfer rate, we use the concept of logarithmic negativity [59], as it is additive. Figure 7 shows the individual rates for each sub-channel. We sum over all of them, weighted according to their frequency of occurance. This leads to an average minimal entanglement transfer rate of 2.17 M log-neg units/s. We also performed the entanglement minimization while taking into account experimental error bars. Using standard error propagation techniques, these measurement uncertainties lead to looser constraint sets for the minimization (this is the same approach taken in [37]). We optimized at 1σ, 2σ and 3σ confidence levels, each giving progressively lower overall entanglement distribution rates.  Figure 6. Effective entanglement minimization performed individually for each transmission sub-channel, for the four-state alphabet at a fixed and optimal amplitude near 1. The negativity bounds in the shaded region were calculated theoretically by assuming all variances equal and thus only serve as a reference for the slightly varying experimental values (see inserted mean variances). At a channel transmission close to 100% (or 83% with detection efficiency considered), the entanglement transfer would be optimal.

Conclusion
We experimentally demonstrated the preservation of quantum properties through the verification of EE for an atmospheric point-to-point link for CV quantum communication in an urban environment. We compare a two-state and four-state alphabet of discretely modulated signal states. The four-state alphabet clearly outperforms the two-states in terms of potentially achievable negativity and tolerable excess noise. Thus we determined the optimal working points for the latter and calculated the rate of distributed entanglement, also taking into account the channel induced drifting shotnoise variance. Note that even though we focused on discrete modulation so far, in principle Gaussian modulation [23,24] can also be easily realized with our system. The shown results clearly indicate the strong potential of free-space channels for CV QKD systems. In addition, they pave the way for further free-space CV quantum communication experiments such as CV entanglement-based QKD ( [60] and see [41] for a demonstration of the free-space distribution of squeezed states) or the use of CV quantum states for atmospheric sensing.   Figure 7. Minimal rate of distributed entanglement per subchannel. By using the logarithmic negativity [59], we can sum over all sub-channels, weighted according to their frequency of occurance. This leads to an overall entanglement transfer rate of 2.17 M log-neg units/s. By expanding the constraint regions in the minimization process, we gain 1σ, 2σ and 3σ intervals, with respectively lower average entanglement distribution rates.