<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.1 plus MathML 2.0//EN" "http://www.w3.org/TR/MathML2/dtd/xhtml-math11-f.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" xmlns:m="http://www.w3.org/1998/Math/MathML" xml:lang="en">
<head><meta http-equiv="Content-Type" content="text/xml; charset=UTF-8"/><title>Secure Coherent-state Quantum Key Distribution Protocols with Efficient Reconciliation</title><link rel="stylesheet" type="text/css" href="document.css"/></head>

<body><h1>Secure Coherent-state Quantum Key Distribution Protocols with Efficient Reconciliation</h1><p class="authors"><span class="author"><a href="mailto:gvanassc@AT@ulb.ac.be">G. Van Assche</a></span><sup>1</sup>, <span class="author">S. Iblisdir</span><sup>1,2</sup> and <span class="author">N. J. Cerf</span><sup>1</sup></p><div class="institutions"><p class="institution"><sup>1</sup> QuIC, Ecole Polytechnique, Université Libre de Bruxelles, CP 165/59, 1050 Brussels, Belgium</p><p class="institution"><sup>2</sup> GAP-Optique, University of Geneva, 20 rue de l'Ecole-de-Médecine, CH-1211 Genève, Switzerland</p></div><h2 id="id2247311">Table of Contents</h2><div class="toc"><p class="toc1"><a href="#id2247587">1 Introduction</a></p><p class="toc1"><a href="#id2248168">2 From Entanglement Purification to Secret Key Distillation</a></p><p class="toc2"><a href="#id2248191">2.1 Binary CSS codes</a></p><p class="toc2"><a href="#id2249156">2.2 Quantum key distribution based on entanglement purification</a></p><p class="toc2"><a href="#id2249788">2.3 Prepare-and-measure quantum key distribution</a></p><p class="toc1"><a href="#id2250009">3 Error Rates Estimation using Tomography</a></p><p class="toc2"><a href="#id2250289">3.1 Estimating phase errors in the average state</a></p><p class="toc2"><a href="#id2251294">3.2 Estimating phase errors using coherent states and homodyne detection</a></p><p class="toc1"><a href="#id2253094">4 Encoding of Multiple QuBits in an Oscillator</a></p><p class="toc2"><a href="#id2253132">4.1 Sliced error correction with invertible mappings</a></p><p class="toc2"><a href="#id2255560">4.2 Quantum sliced error correction</a></p><p class="toc3"><a href="#id2255966">4.2.1 The mappings <span class="eq"><m:math><m:mi mathvariant="script">QS</m:mi></m:math></span> and <span class="eq"><m:math><m:mi mathvariant="script">QE</m:mi></m:math></span></a></p><p class="toc3"><a href="#id2259027">4.2.2 Phase coherence</a></p><p class="toc3"><a href="#id2259600">4.2.3 Construction of <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> and <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover></m:math></span></a></p><p class="toc1"><a href="#id2263219">5 The Attenuation Channel</a></p><p class="toc1"><a href="#id2264357">6 Asymptotic Behavior</a></p><p class="toc2"><a href="#id2264429">6.1 Direct reconciliation</a></p><p class="toc2"><a href="#id2268521">6.2 Reverse reconciliation</a></p><p class="toc1"><a href="#id2269346">7 Conclusion</a></p><p class="toc1"><a href="#id2247160">Acknowledgments</a></p><p class="toc1"><a href="#id2247175">Bibliography</a></p></div><h2>Abstract</h2><div class="abstract">
<p>We study the equivalence between a realistic quantum key distribution protocol using coherent states and homodyne detection and a formal entanglement purification protocol. Maximally-entangled qubit pairs that one can extract in the formal protocol correspond to secret key bits in the realistic protocol. More specifically, we define a qubit encoding scheme that allows the formal protocol to produce more than one entangled qubit pair per coherent state, or equivalently for the realistic protocol, more than one secret key bit. The entanglement parameters are estimated using quantum tomography. We analyze the properties of the encoding scheme and investigate its application to the important case of the attenuation channel.</p>
</div><h2 id="id2247587">1 Introduction</h2><p>Quantum key distribution (QKD), also called quantum cryptography, allows two parties, Alice and Bob, to share a secret key that can be used for encrypting messages using a classical cipher, e.g., the one-time pad. The main interest of such a key distribution scheme is that any eavesdropping is, in principle, detectable as the laws of quantum mechanics imply that measuring a quantum state generally disturbs it.</p><p>The resources required by QKD comprise a source of non-orthogonal quantum states on Alice's side, a quantum channel conveying these states to Bob, a measuring apparatus on Bob's side, and a (public) authenticated classical channel between Alice and Bob. In addition to being used to generate a secret key, the quantum channel is subject to probing by the legitimate parties, so as to determine how many secret bits can be generated.</p><p>Most interest in QKD has been devoted to protocols involving (an approximation to) a single-photon source on Alice's side and a single-photon detector on Bob's side (see <span class="bib-cite">[<a href="#id2247184">1</a>]</span> and the references therein). However, protocols involving quantum continuous variables have been considered with an increasing interest <span class="bib-cite">[<a href="#id2269712">2</a>, <a href="#id2269752">3</a>, <a href="#id2269800">4</a>, <a href="#id2269853">5</a>]</span>. Of special importance are coherent-state protocols <span class="bib-cite">[<a href="#id2269906">6</a>, <a href="#id2269953">7</a>]</span>. The quantum source at Alice's side randomly generates coherent states of a light mode with Gaussian-distributed quadratures, and Bob's measurements are homodyne measurements. These protocols seem to allow for facilitated implementations and higher secret-key generation rates than the protocols involving single-photon sources <span class="bib-cite">[<a href="#id2269953">7</a>]</span>.</p><p>Consequently, there is an increasing interest for studying the security of coherent state protocols under general classes of attacks. Individual Gaussian attacks are considered in <span class="bib-cite">[<a href="#id2269906">6</a>, <a href="#id2269953">7</a>]</span>, and are found to be optimal in the more general class of finite-width non-Gaussian incoherent attacks <span class="bib-cite">[<a href="#id2270025">8</a>]</span>. Individually-probed collective attacks are also considered  in <span class="bib-cite">[<a href="#id2270073">9</a>, <a href="#id2270108">10</a>]</span>. The recent techniques of <span class="bib-cite">[<a href="#id2270149">11</a>, <a href="#id2270197">12</a>]</span> do not make any assumptions on the eavesdropper's technology and are also considered in <span class="bib-cite">[<a href="#id2270073">9</a>, <a href="#id2270108">10</a>]</span> for coherent state protocols, although giving lower secret key rates.</p><p>In this paper, we study the security of a prepare-and-measure QKD protocol <span class="bib-cite">[<a href="#id2269906">6</a>, <a href="#id2269953">7</a>]</span> by establishing its equivalence to an entanglement purification (EP) protocol, which produces maximally-entangled qubit pairs. A maximally-entangled qubit pair is by definition completely factored from its environment, and thus the values obtained by measuring each side are fully correlated and secret. The equivalent prepare-and-measure QKD protocol also enjoys this property. This particular technique thus allows one to relieve from any assumptions on the eavesdropper's strategy and was used in <span class="bib-cite">[<a href="#id2270238">13</a>]</span> to assess the security of the BB84 protocol and in <span class="bib-cite">[<a href="#id2269752">3</a>]</span> for a squeezed state protocol. More recently, this technique was extended to the case of coherent state protocols <span class="bib-cite">[<a href="#id2270285">14</a>]</span>.</p><p>To show the equivalence between a QKD protocol and of an EP protocol, one has to explicitly take into account the secret key distillation, that is, the techniques used to make Alice's and Bob's keys equal (reconciliation) and fully secret (privacy amplification). In <span class="bib-cite">[<a href="#id2270238">13</a>]</span>, the EP protocol uses CSS quantum codes <span class="bib-cite">[<a href="#id2270327">15</a>, <a href="#id2270373">16</a>]</span>, which are equivalent in QKD to reconciliation with syndromes of binary linear codes and privacy amplification by multiplication with a parity-check matrix. In contrast to BB84, the modulation of coherent states in the protocol we consider here is <em>continuous</em>, therefore producing continuous key elements from which to extract a secret key. Reconciliation of a Gaussian-distributed key was studied in <span class="bib-cite">[<a href="#id2270414">17</a>]</span>, and a generic protocol called sliced error correction was designed so as to distill a <em>binary</em> key.</p><p>In contrast to <span class="bib-cite">[<a href="#id2270285">14</a>]</span>, the EP protocol is constructed in such a way that it is equivalent to a QKD protocol with sliced error correction for reconciliation. The advantage is the higher secret key rate and the better resistance to attenuation that one can achieve. In particular, more than one maximally-entangled pair (or secret key bit) can be produced per coherent state. Furthermore, thanks to its generality, the asymptotic efficiency of the EP protocol inherits to some extent from the asymptotic efficiency of the classical reconciliation protocol.</p><p>The paper is organized as follows. First, in Sec. <span class="ref"><a href="#id2248168">2</a></span>, we macroscopically describe the formal EP protocol and its equivalent QKD protocol that are used throughout the paper. Then, in Sec. <span class="ref"><a href="#id2250009">3</a></span>, we show how the channel can be probed so as to determine the number of secret key bits that Alice and Bob can generate. The encoding of qubits, that is, the generalization of sliced error correction to EP, is described in Sec. <span class="ref"><a href="#id2253094">4</a></span>. Then, Sec. <span class="ref"><a href="#id2263219">5</a></span> deals with the important particular case of an attenuation channel. Finally, the asymptotic properties of the qubit encoding scheme are detailed in Sec. <span class="ref"><a href="#id2264357">6</a></span>.</p><h2 id="id2248168">2 From Entanglement Purification to Secret Key Distillation</h2><p>After we review the case of EP using CSS codes and its equivalence to BB84, we give a high-level description of a QKD protocol based on EP. We consider this protocol as formal, that is, we do not expect a physical implementation of it. Instead, we propose a prepare-and-measure QKD protocol, derived from the formal one, which also encompasses error correction and privacy amplification.</p><h3 id="id2248191">2.1 Binary CSS codes</h3><p>In the case of BB84, the CSS codes can readily be used to establish the equivalence between an EP protocol and a QKD protocol <span class="bib-cite">[<a href="#id2270238">13</a>]</span>. Since we will use CSS codes as an ingredient for the EP and QKD protocols below, let us briefly review their properties.</p><p>Starting from the Einstein-Podolski-Rosen (EPR) state
<table class="equation" id="id2248218"><tr><td class="eq"><m:math><m:mrow><m:mo>|</m:mo><m:msup><m:mi>φ</m:mi><m:mo>+</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mo>=</m:mo><m:msup><m:mn>2</m:mn><m:mrow><m:mo>-</m:mo><m:mn>1</m:mn><m:mo>/</m:mo><m:mn>2</m:mn></m:mrow></m:msup><m:mo>⁢</m:mo><m:mo stretchy="false">(</m:mo><m:mrow><m:mo>|</m:mo><m:mn>0</m:mn><m:mn>0</m:mn><m:mo>〉</m:mo></m:mrow><m:mo>+</m:mo><m:mrow><m:mo>|</m:mo><m:mn>1</m:mn><m:mn>1</m:mn><m:mo>〉</m:mo></m:mrow><m:mo stretchy="false">)</m:mo><m:mtext>,</m:mtext></m:math></td><td class="number">(2.1)</td></tr></table>
Alice keeps half of the state and sends the other half to Bob. His part may undergo a bit error (<span class="eq"><m:math><m:mrow><m:mo>|</m:mo><m:msup><m:mi>φ</m:mi><m:mo>+</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mo>→</m:mo><m:mrow><m:mo>|</m:mo><m:msup><m:mi>ψ</m:mi><m:mo>+</m:mo></m:msup><m:mo>〉</m:mo></m:mrow></m:math></span>), phase error (<span class="eq"><m:math><m:mrow><m:mo>|</m:mo><m:msup><m:mi>φ</m:mi><m:mo>+</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mo>→</m:mo><m:mrow><m:mo>|</m:mo><m:msup><m:mi>φ</m:mi><m:mo>-</m:mo></m:msup><m:mo>〉</m:mo></m:mrow></m:math></span>) or both errors (<span class="eq"><m:math><m:mrow><m:mo>|</m:mo><m:msup><m:mi>φ</m:mi><m:mo>+</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mo>→</m:mo><m:mrow><m:mo>|</m:mo><m:msup><m:mi>ψ</m:mi><m:mo>-</m:mo></m:msup><m:mo>〉</m:mo></m:mrow></m:math></span>), with <span class="eq"><m:math><m:mrow><m:mo>|</m:mo><m:msup><m:mi>φ</m:mi><m:mo>-</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mo>=</m:mo><m:msup><m:mn>2</m:mn><m:mrow><m:mo>-</m:mo><m:mn>1</m:mn><m:mo>/</m:mo><m:mn>2</m:mn></m:mrow></m:msup><m:mo>⁢</m:mo><m:mo stretchy="false">(</m:mo><m:mrow><m:mo>|</m:mo><m:mn>0</m:mn><m:mn>0</m:mn><m:mo>〉</m:mo></m:mrow><m:mo>-</m:mo><m:mrow><m:mo>|</m:mo><m:mn>1</m:mn><m:mn>1</m:mn><m:mo>〉</m:mo></m:mrow><m:mo stretchy="false">)</m:mo></m:math></span> and <span class="eq"><m:math><m:mrow><m:mo>|</m:mo><m:msup><m:mi>ψ</m:mi><m:mo>±</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mo>=</m:mo><m:msup><m:mn>2</m:mn><m:mrow><m:mo>-</m:mo><m:mn>1</m:mn><m:mo>/</m:mo><m:mn>2</m:mn></m:mrow></m:msup><m:mo>⁢</m:mo><m:mo stretchy="false">(</m:mo><m:mrow><m:mo>|</m:mo><m:mn>0</m:mn><m:mn>1</m:mn><m:mo>〉</m:mo></m:mrow><m:mo>±</m:mo><m:mrow><m:mo>|</m:mo><m:mn>1</m:mn><m:mn>0</m:mn><m:mo>〉</m:mo></m:mrow><m:mo stretchy="false">)</m:mo></m:math></span>. Given that not too many such errors occurs, Alice and Bob can obtain, from many instances of such a transmitted state, a smaller number of EPR pairs using only local operations and classical communications (LOCC). One way to do this is to use CSS codes.</p><p>Let <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>1</m:mn></m:msub></m:math></span> and <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>2</m:mn></m:msub></m:math></span> be two binary error correcting codes of <span class="eq"><m:math><m:mi>n</m:mi></m:math></span> bits (i.e., <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>1</m:mn></m:msub></m:math></span> and <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>2</m:mn></m:msub></m:math></span> are vector spaces of <span class="eq"><m:math><m:msubsup><m:mi mathvariant="bold">F</m:mi><m:mn>2</m:mn><m:mi>n</m:mi></m:msubsup></m:math></span>) with parity check matrices <span class="eq"><m:math><m:msub><m:mi>H</m:mi><m:mn>1</m:mn></m:msub></m:math></span> and <span class="eq"><m:math><m:msub><m:mi>H</m:mi><m:mn>2</m:mn></m:msub></m:math></span>, resp. They are chosen such that <span class="eq"><m:math><m:mo stretchy="false">{</m:mo><m:mn>0</m:mn><m:mo stretchy="false">}</m:mo><m:mo>⊂</m:mo><m:msub><m:mi>C</m:mi><m:mn>2</m:mn></m:msub><m:mo>⊂</m:mo><m:msub><m:mi>C</m:mi><m:mn>1</m:mn></m:msub><m:mo>⊂</m:mo><m:msubsup><m:mi mathvariant="bold">F</m:mi><m:mn>2</m:mn><m:mi>n</m:mi></m:msubsup></m:math></span>. A CSS code is a <span class="eq"><m:math><m:mi>k</m:mi></m:math></span>-dimensional subspace of <span class="eq"><m:math><m:msup><m:mi mathvariant="script">H</m:mi><m:mi>n</m:mi></m:msup></m:math></span>, the Hilbert space of <span class="eq"><m:math><m:mi>n</m:mi></m:math></span> qubits, with <span class="eq"><m:math><m:mi>k</m:mi><m:mo>=</m:mo><m:mi>dim</m:mi><m:mo>⁡</m:mo><m:msub><m:mi>C</m:mi><m:mn>1</m:mn></m:msub><m:mo>-</m:mo><m:mi>dim</m:mi><m:mo>⁡</m:mo><m:msub><m:mi>C</m:mi><m:mn>2</m:mn></m:msub></m:math></span> <span class="bib-cite">[<a href="#id2270327">15</a>, <a href="#id2270373">16</a>]</span>. The code <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>1</m:mn></m:msub></m:math></span> allows to correct bit errors, while <span class="eq"><m:math><m:msubsup><m:mi>C</m:mi><m:mn>2</m:mn><m:mo>⊥</m:mo></m:msubsup></m:math></span> (the dual code of <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>2</m:mn></m:msub></m:math></span>) allows to correct phase errors—one important property of the CSS codes is to be able to correct bit errors and phase errors independently.</p><p>For entanglement purification, Alice and Bob must compare their syndromes, both for bit errors and phase errors. The relative syndrome determines the correction that Bob must apply to align his qubits to Alice's. Translating this into the BB84 protocol, one can show <span class="bib-cite">[<a href="#id2270238">13</a>]</span> that the relative syndrome for bit errors in the EP protocol is equal to the relative syndrome for bit errors that Alice and Bob would have reconciled in the BB84 protocol. So, reconciliation can be done using the <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>1</m:mn></m:msub></m:math></span> code. Phase errors of the EP protocol do not have such a direct equivalent in the BB84 protocol: The prepare-and-measure protocol works as if Alice measured her part of the state in the <span class="eq"><m:math><m:mo stretchy="false">{</m:mo><m:mo>|</m:mo><m:mn>0</m:mn><m:mo>〉</m:mo><m:mo>,</m:mo><m:mo>|</m:mo><m:mn>1</m:mn><m:mo>〉</m:mo><m:mo stretchy="false">}</m:mo></m:math></span> basis, thereby discarding information on the phase. However, one does not really need to correct the phase errors in the BB84 protocol. Instead, if <span class="eq"><m:math><m:msubsup><m:mi>C</m:mi><m:mn>2</m:mn><m:mo>⊥</m:mo></m:msubsup></m:math></span> would be able to correct them in the EP protocol, the syndrome of <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>2</m:mn></m:msub></m:math></span> in <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>1</m:mn></m:msub></m:math></span> of Alice and Bob's bit string turns out to be a valid secret key in the prepare-and-measure protocol. Stated otherwise, <span class="eq"><m:math><m:msub><m:mi>H</m:mi><m:mn>1</m:mn></m:msub></m:math></span> determines the syndrome Alice has to send to Bob to perform reconciliation, while <span class="eq"><m:math><m:msub><m:mi>H</m:mi><m:mn>2</m:mn></m:msub></m:math></span> determines the way the final key is computed for privacy amplification.</p><p>Overall, the number of secret key bits is thus <span class="eq"><m:math><m:mi>k</m:mi><m:mo>=</m:mo><m:mi>dim</m:mi><m:mo>⁡</m:mo><m:msub><m:mi>C</m:mi><m:mn>1</m:mn></m:msub><m:mo>-</m:mo><m:mi>dim</m:mi><m:mo>⁡</m:mo><m:msub><m:mi>C</m:mi><m:mn>2</m:mn></m:msub></m:math></span>, provided that <span class="eq"><m:math><m:msub><m:mi>C</m:mi><m:mn>1</m:mn></m:msub></m:math></span> (resp. <span class="eq"><m:math><m:msubsup><m:mi>C</m:mi><m:mn>2</m:mn><m:mo>⊥</m:mo></m:msubsup></m:math></span>) is small enough to correct all the bit (resp. phase) errors. When considering asymptotically large block sizes, the CSS codes can produce
<table class="equation" id="id2248893"><tr><td class="eq"><m:math><m:mi>k</m:mi><m:mo>=</m:mo><m:mi>r</m:mi><m:mo>⁢</m:mo><m:mi>n</m:mi><m:mo>→</m:mo><m:mi>n</m:mi><m:mo>⁢</m:mo><m:mo stretchy="false">(</m:mo><m:mn>1</m:mn><m:mo>-</m:mo><m:mi>h</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>e</m:mi><m:mtext>b</m:mtext></m:msup><m:mo stretchy="false">)</m:mo><m:mo>-</m:mo><m:mi>h</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>e</m:mi><m:mtext>p</m:mtext></m:msup><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:mi>R</m:mi><m:mo>⁢</m:mo><m:mi>n</m:mi></m:math></td><td class="number">(2.2)</td></tr></table>
EPR pairs or secret key bits, with <span class="eq"><m:math><m:msup><m:mi>e</m:mi><m:mtext>b</m:mtext></m:msup></m:math></span> (resp. <span class="eq"><m:math><m:msup><m:mi>e</m:mi><m:mtext>p</m:mtext></m:msup></m:math></span>) the bit (resp. phase) error rate and <span class="eq"><m:math><m:mi>h</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:mo>-</m:mo><m:mi>p</m:mi><m:mo>⁢</m:mo><m:msub><m:mi>log</m:mi><m:mn>2</m:mn></m:msub><m:mo>⁡</m:mo><m:mi>p</m:mi><m:mo>-</m:mo><m:mo stretchy="false">(</m:mo><m:mn>1</m:mn><m:mo>-</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mo>⁢</m:mo><m:msub><m:mi>log</m:mi><m:mn>2</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mn>1</m:mn><m:mo>-</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> <span class="bib-cite">[<a href="#id2270238">13</a>]</span>.</p><p>We conclude this section by noting that the bit error rate <span class="eq"><m:math><m:msup><m:mi>e</m:mi><m:mtext>b</m:mtext></m:msup></m:math></span> determines the number of bits revealed by reconciliation (asymptotically <span class="eq"><m:math><m:mi>h</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>e</m:mi><m:mtext>b</m:mtext></m:msup><m:mo stretchy="false">)</m:mo></m:math></span>), whereas the phase error rate <span class="eq"><m:math><m:msup><m:mi>e</m:mi><m:mtext>p</m:mtext></m:msup></m:math></span> determines the number of bits discarded by privacy amplification due to eavesdropping (asymptotically <span class="eq"><m:math><m:mi>h</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>e</m:mi><m:mtext>p</m:mtext></m:msup><m:mo stretchy="false">)</m:mo></m:math></span>).</p><h3 id="id2249156">2.2 Quantum key distribution based on entanglement purification</h3><p>In BB84, the modulation of qubits can be transposed as if Alice prepares a <span class="eq"><m:math><m:mo>|</m:mo><m:msup><m:mi>φ</m:mi><m:mo>+</m:mo></m:msup><m:mo>〉</m:mo></m:math></span> state and measures her part. In the case of the QKD protocol with Gaussian-modulated coherent states, the formal state that Alice prepares is of course different, as it must reduce to the proper modulation when Alice measures her part. We define the formal state as:</p><table class="equation" id="id2249189"><tr><td class="eq"><m:math><m:mo>|</m:mo><m:mi>Ψ</m:mi><m:mo>〉</m:mo><m:mo>=</m:mo><m:mo>∫</m:mo><m:mo>ⅆ</m:mo><m:mi>x</m:mi><m:mo>ⅆ</m:mo><m:mi>p</m:mi><m:mspace width="0.222em"/><m:mi>g</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>,</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:msub><m:mo>⊗</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mi>p</m:mi><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:msub><m:mo>⊗</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo>〉</m:mo></m:mrow><m:mi mathvariant="sans-serif">b</m:mi></m:msub><m:mtext>,</m:mtext></m:math></td><td class="number">(2.3)</td></tr></table><p>where <span class="eq"><m:math><m:mi>g</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>,</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> denotes a bi-variate Gaussian distribution <span class="eq"><m:math><m:mrow><m:mi>g</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>,</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>=</m:mo><m:msqrt><m:msub><m:mi>G</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>⁢</m:mo><m:msub><m:mi>G</m:mi><m:mn>2</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:msqrt></m:math></span>. The kets <span class="eq"><m:math><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>〉</m:mo></m:math></span>, <span class="eq"><m:math><m:mo>|</m:mo><m:mi>p</m:mi><m:mo>〉</m:mo></m:math></span>, <span class="eq"><m:math><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo>〉</m:mo></m:math></span> are shorthand notations for respectively a <span class="eq"><m:math><m:mi mathvariant="bold">x</m:mi></m:math></span>-quadrature eigenstate with eigenvalue <span class="eq"><m:math><m:mi>x</m:mi></m:math></span>, a <span class="eq"><m:math><m:mi mathvariant="bold">p</m:mi></m:math></span>-quadrature eigenstate with eigenvalue <span class="eq"><m:math><m:mi>p</m:mi></m:math></span> and a coherent state whose <span class="eq"><m:math><m:mi mathvariant="bold">x</m:mi></m:math></span> mean value equals <span class="eq"><m:math><m:mi>x</m:mi></m:math></span> and whose <span class="eq"><m:math><m:mi mathvariant="bold">p</m:mi></m:math></span> mean value equals <span class="eq"><m:math><m:mi>p</m:mi></m:math></span>. The subscripts <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:math></span>, <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:math></span> (resp. <span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span>) denote that the system is lying on Alice's side (resp. Bob's side).</p><p>The state (<span class="ref"><a href="#id2249189">2.3</a></span>) does not have a direct physical meaning. In particular, the systems <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:math></span> and <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:math></span> must be understood as classical pointers, e.g., resulting from the (formal) homodyne measurement of an EPR state as studied in <span class="bib-cite">[<a href="#id2270466">18</a>]</span>.</p><p>In the entanglement purification picture, the <span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span> part of the system is sent to Bob (and possibly attacked by Eve) and the <span class="eq"><m:math><m:mi mathvariant="sans-serif">a</m:mi></m:math></span> part stays at Alice's station. If Alice measures <span class="eq"><m:math><m:mi mathvariant="bold">x</m:mi></m:math></span> in <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:math></span> and <span class="eq"><m:math><m:mi mathvariant="bold">p</m:mi></m:math></span> in <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:math></span>, the state is projected as if Alice sent Bob a coherent state centered on <span class="eq"><m:math><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi></m:math></span>.</p><p>Let us now describe the EP protocol, which reduces to the prepare-and-measure QKD protocol described further.</p><ul>
<li>Alice creates <span class="eq"><m:math><m:mi>l</m:mi><m:mo>+</m:mo><m:mi>n</m:mi></m:math></span> copies of the state <span class="eq"><m:math><m:mo>|</m:mo><m:mi>Ψ</m:mi><m:mo>〉</m:mo></m:math></span>, of which she sends the <span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span> part to Bob.</li>

<li>Bob acknowledges reception of the states.</li>
<li>Out of the <span class="eq"><m:math><m:mi>l</m:mi><m:mo>+</m:mo><m:mi>n</m:mi></m:math></span> states, <span class="eq"><m:math><m:mi>n</m:mi></m:math></span> will serve for estimation purposes. These states are chosen randomly and uniformly by Alice, who informs Bob about their positions.</li>
<li>For the remaining <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> states, Alice and Bob perform entanglement purification, so as to produce <span class="eq"><m:math><m:mi>r</m:mi><m:mo>⁢</m:mo><m:mi>l</m:mi></m:math></span> (<span class="eq"><m:math><m:mn>0</m:mn><m:mo>≤</m:mo><m:mi>r</m:mi><m:mo>≤</m:mo><m:mn>1</m:mn></m:math></span>) states very close to <span class="eq"><m:math><m:mo>|</m:mo><m:msup><m:mi>φ</m:mi><m:mo>+</m:mo></m:msup><m:mo>〉</m:mo></m:math></span>. Measured in the computational bases, the produced states yield <span class="eq"><m:math><m:mi>r</m:mi><m:mo>⁢</m:mo><m:mi>l</m:mi></m:math></span> secret bits on both Alice's and Bob's sides.</li>
</ul><p>The details of the EP procedure, which uses CSS codes as an ingredient, are given in Sec. <span class="ref"><a href="#id2253094">4</a></span>, while the estimation is detailed in Sec. <span class="ref"><a href="#id2250009">3</a></span>.</p><h3 id="id2249788">2.3 Prepare-and-measure quantum key distribution</h3><p>By virtually measuring the <span class="eq"><m:math><m:mi mathvariant="sans-serif">a</m:mi></m:math></span> part of the state <span class="eq"><m:math><m:mo>|</m:mo><m:mi>Ψ</m:mi><m:mo>〉</m:mo></m:math></span>, the protocol above reduces to the following one.</p><ul>
<li>Alice modulates <span class="eq"><m:math><m:mi>l</m:mi><m:mo>+</m:mo><m:mi>n</m:mi></m:math></span> coherent states <span class="eq"><m:math><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo>〉</m:mo></m:math></span> that she sends to Bob. The choice of the values of <span class="eq"><m:math><m:mi>x</m:mi></m:math></span> and <span class="eq"><m:math><m:mi>p</m:mi></m:math></span> follow the distribution <span class="eq"><m:math><m:msup><m:mrow><m:mo>|</m:mo><m:mi>g</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>,</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mo>|</m:mo></m:mrow><m:mn>2</m:mn></m:msup><m:mo>=</m:mo><m:msub><m:mi>G</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>⁢</m:mo><m:msub><m:mi>G</m:mi><m:mn>2</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>. </li>
<li>Bob acknowledges reception of the states.</li>
<li>Out of the <span class="eq"><m:math><m:mi>l</m:mi><m:mo>+</m:mo><m:mi>n</m:mi></m:math></span> states, <span class="eq"><m:math><m:mi>n</m:mi></m:math></span> will serve for estimation purposes. These states are chosen randomly and uniformly by Alice, who informs Bob about their positions.</li>
<li>For the remaining <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> states, Bob measures <span class="eq"><m:math><m:mi mathvariant="bold">x</m:mi></m:math></span>. Alice and Bob perform secret key distillation (reconciliation and privacy amplification), so as to produce <span class="eq"><m:math><m:mi>r</m:mi><m:mo>⁢</m:mo><m:mi>l</m:mi></m:math></span> secret bits.</li>
</ul><p>The reconciliation and privacy amplification procedures are based on classical error correcting codes, which derive from the CSS codes used in the formal EP protocol.</p><h2 id="id2250009">3 Error Rates Estimation using Tomography</h2><p>In QKD protocols derived from EP, an important step is to show how one can infer the bit and phase error rates of the samples that compose the key. A fraction of the samples sent by Alice to Bob are sacrificed so as to serve as test samples. By randomly choosing them within the stream of data, they are statistically representative of the whole stream.</p><p>In <span class="bib-cite">[<a href="#id2269752">3</a>, <a href="#id2270238">13</a>]</span>, one can simply make measurements and directly count the number of bit and phase errors from the results. This is possible since Bob's apparatus can measure both bit and phase values. In <span class="bib-cite">[<a href="#id2270285">14</a>]</span>, however, it is not possible to measure directly phase errors. Yet some data post-processing can be applied on measurements so as to infer the number of phase errors in the stream of data. In this section, we wish to show that we can extend this to more general (and more efficient) encodings of qubits (in the EP picture) or bits (in the derived QKD protocol).</p><p>The encoding of bits will be described in a further section—for the moment, the qubit pair system (i.e., one among possibly several ones), which Alice and Bob will process using CSS codes, is abstractly represented by its Pauli operators acting in <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:math></span>: <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">Z</m:mi><m:mi mathvariant="sans-serif">s</m:mi></m:msub></m:math></span> (phase flip) and <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">s</m:mi></m:msub></m:math></span> (bit flip), and in
<span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span>: <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">Z</m:mi><m:mi mathvariant="sans-serif">e</m:mi></m:msub></m:math></span> and <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">e</m:mi></m:msub></m:math></span>. (The subscripts <span class="eq"><m:math><m:mi mathvariant="sans-serif">s</m:mi></m:math></span> and <span class="eq"><m:math><m:mi mathvariant="sans-serif">e</m:mi></m:math></span> stand for slice and estimator, resp., to follow the convention of the following sections.)
The bit errors are assumed to be easy to determine, that is, <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">Z</m:mi><m:mi mathvariant="sans-serif">s</m:mi></m:msub></m:math></span> has a diagonal expansion in <span class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:msub><m:mo>〈</m:mo><m:mi>x</m:mi><m:mo>|</m:mo></m:math></span>, and <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">Z</m:mi><m:mi mathvariant="sans-serif">e</m:mi></m:msub></m:math></span> can directly be determined by a single homodyne measurement on <span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span>. This ensures, in the derived prepare-and-measure QKD protocol, that Alice knows the bit value she sent, and Bob can determine the received bit value. A measurement of the observable <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">s</m:mi></m:msub><m:msub><m:mi mathvariant="bold">I</m:mi><m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:mrow></m:msub><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">e</m:mi></m:msub></m:math></span> associated to the phase error rate, however, cannot be implemented by a single homodyne measurement on <span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span>. Therefore, we
have to invoke quantum tomography with a quorum of operators <span class="bib-cite">[<a href="#id2270532">19</a>]</span> to get an estimate of the phase error rate.
</p><h3 id="id2250289">3.1 Estimating phase errors in the average state</h3><p>In the EP picture, let <span class="eq"><m:math><m:msup><m:mi>ρ</m:mi><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>n</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msup></m:math></span> be the state of the <span class="eq"><m:math><m:mi>n</m:mi></m:math></span> samples used for estimation of the phase error rate (i.e., <span class="eq"><m:math><m:mi>n</m:mi></m:math></span> instances of the <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub><m:mi mathvariant="sans-serif">b</m:mi></m:math></span> system). To count the number of phase errors in a set of <span class="eq"><m:math><m:mi>n</m:mi></m:math></span> samples, one needs to measure
<span class="eq"><m:math><m:mi mathvariant="bold">O</m:mi><m:mo>=</m:mo><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">s</m:mi></m:msub><m:msub><m:mi mathvariant="bold">I</m:mi><m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:mrow></m:msub><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">e</m:mi></m:msub></m:math></span> on the <span class="eq"><m:math><m:mi>n</m:mi></m:math></span> samples and sum the results (with <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">I</m:mi><m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:mrow></m:msub></m:math></span> the identity in the system <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:math></span>).
This is equivalent to measuring <span class="eq"><m:math><m:msup><m:mi mathvariant="bold">O</m:mi><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>n</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msup><m:mo>=</m:mo><m:msub><m:mo>∑</m:mo><m:mi>i</m:mi></m:msub><m:msubsup><m:mi mathvariant="bold">I</m:mi><m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub><m:mi mathvariant="sans-serif">b</m:mi></m:mrow><m:mrow><m:mo>⊗</m:mo><m:mi>i</m:mi><m:mo>-</m:mo><m:mn>1</m:mn></m:mrow></m:msubsup><m:mo>⊗</m:mo><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">s</m:mi></m:msub><m:msub><m:mi mathvariant="bold">I</m:mi><m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:mrow></m:msub><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">e</m:mi></m:msub><m:mo>⊗</m:mo><m:msubsup><m:mi mathvariant="bold">I</m:mi><m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub><m:mi mathvariant="sans-serif">b</m:mi></m:mrow><m:mrow><m:mo>⊗</m:mo><m:mi>n</m:mi><m:mo>-</m:mo><m:mi>i</m:mi></m:mrow></m:msubsup></m:math></span>. If the true phase error probability in the <span class="eq"><m:math><m:mi>n</m:mi><m:mo>+</m:mo><m:mi>l</m:mi></m:math></span> samples is <span class="eq"><m:math><m:msup><m:mi>e</m:mi><m:mtext>p</m:mtext></m:msup></m:math></span>, the error variance is <span class="eq"><m:math><m:msubsup><m:mi>σ</m:mi><m:mn>1</m:mn><m:mn>2</m:mn></m:msubsup><m:mo>=</m:mo><m:mn>2</m:mn><m:mo>⁢</m:mo><m:msup><m:mi>e</m:mi><m:mtext>p</m:mtext></m:msup><m:mo stretchy="false">(</m:mo><m:mn>1</m:mn><m:mo>-</m:mo><m:msup><m:mi>e</m:mi><m:mtext>p</m:mtext></m:msup><m:mo stretchy="false">)</m:mo><m:mo>/</m:mo><m:mi>n</m:mi></m:math></span>, and thus the probability of making an estimation error of more than <span class="eq"><m:math><m:mi>Δ</m:mi></m:math></span> is <span class="bib-cite">[<a href="#id2269752">3</a>, <a href="#id2270238">13</a>]</span> asymptotically <span class="eq"><m:math><m:mi>exp</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mo>-</m:mo><m:msup><m:mi>Δ</m:mi><m:mn>2</m:mn></m:msup><m:mi>n</m:mi><m:mo>/</m:mo><m:mn>4</m:mn><m:msup><m:mi>e</m:mi><m:mtext>p</m:mtext></m:msup><m:mo stretchy="false">(</m:mo><m:mn>1</m:mn><m:mo>-</m:mo><m:msup><m:mi>e</m:mi><m:mtext>p</m:mtext></m:msup><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo></m:math></span>. It is easy to see that
<table class="equation" id="id2250724"><tr><td class="eq"><m:math><m:mi>Tr</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi mathvariant="bold">O</m:mi><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>n</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msup><m:mo>⁢</m:mo><m:msup><m:mi>ρ</m:mi><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>n</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msup><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:mi>n</m:mi><m:mo>⁢</m:mo><m:mi>Tr</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi mathvariant="bold">O</m:mi><m:mo>⁢</m:mo><m:mi>ρ</m:mi><m:mo stretchy="false">)</m:mo><m:mtext>,</m:mtext></m:math></td><td class="number">(3.1)</td></tr></table>
where <span class="eq"><m:math><m:mi>ρ</m:mi><m:mo>=</m:mo><m:msup><m:mi>n</m:mi><m:mrow><m:mo>-</m:mo><m:mn>1</m:mn></m:mrow></m:msup><m:mo>⁢</m:mo><m:msub><m:mo>∑</m:mo><m:mi>i</m:mi></m:msub><m:msub><m:mi>Tr</m:mi><m:mrow><m:mtext>All</m:mtext><m:mo>∖</m:mo><m:mo stretchy="false">{</m:mo><m:mi>i</m:mi><m:mo stretchy="false">}</m:mo></m:mrow></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>ρ</m:mi><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>n</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msup><m:mo stretchy="false">)</m:mo></m:math></span> is the density matrix of the average state measured. So, we can estimate the number of phase errors using the average state, even if eavesdropping is joint (<span class="eq"><m:math><m:msup><m:mi>ρ</m:mi><m:mrow><m:mo>⊗</m:mo><m:mi>n</m:mi></m:mrow></m:msup><m:mo>≠</m:mo><m:msup><m:mi>ρ</m:mi><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>n</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msup></m:math></span>).
</p><p>If the measurement of
<span class="eq"><m:math><m:mi mathvariant="bold">O</m:mi><m:mo>=</m:mo><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">s</m:mi></m:msub><m:msub><m:mi mathvariant="bold">I</m:mi><m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:mrow></m:msub><m:msub><m:mi mathvariant="bold">X</m:mi><m:mi mathvariant="sans-serif">e</m:mi></m:msub></m:math></span> cannot be made directly, one instead looks for a quorum of operators <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">Q</m:mi><m:mi>λ</m:mi></m:msub></m:math></span> such that
<span class="eq"><m:math><m:mi mathvariant="bold">O</m:mi><m:mo>=</m:mo><m:mo>∫</m:mo><m:mo>ⅆ</m:mo><m:mi>λ</m:mi><m:mo>⁢</m:mo><m:mi>o</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>λ</m:mi><m:mo stretchy="false">)</m:mo><m:mo>⁢</m:mo><m:msub><m:mi mathvariant="bold">Q</m:mi><m:mi>λ</m:mi></m:msub></m:math></span>;
estimating <span class="eq"><m:math><m:mo>〈</m:mo><m:mi mathvariant="bold">O</m:mi><m:mo>〉</m:mo></m:math></span> comes down to measuring several times
<span class="eq"><m:math><m:msub><m:mi mathvariant="bold">Q</m:mi><m:mi>λ</m:mi></m:msub></m:math></span> for values of <span class="eq"><m:math><m:mi>λ</m:mi></m:math></span> chosen randomly and independently of each other, and averaging the results weighted by <span class="eq"><m:math><m:mi>o</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>λ</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>: <span class="eq"><m:math><m:mi mathvariant="bold">O</m:mi><m:mo>≈</m:mo><m:msub><m:mo>∑</m:mo><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mi>o</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>λ</m:mi><m:mi>i</m:mi></m:msub><m:mo stretchy="false">)</m:mo><m:mo>⁢</m:mo><m:msub><m:mi mathvariant="bold">Q</m:mi><m:mrow><m:msub><m:mi>λ</m:mi><m:mi>i</m:mi></m:msub></m:mrow></m:msub></m:math></span> <span class="bib-cite">[<a href="#id2270532">19</a>]</span>.
If the values of <span class="eq"><m:math><m:mi>λ</m:mi></m:math></span> are chosen independently of the sample index on which <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">Q</m:mi><m:mi>λ</m:mi></m:msub></m:math></span> is applied, we get unbiased results, as <span class="eq"><m:math><m:mi>Tr</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi mathvariant="bold">O</m:mi><m:mo>⁢</m:mo><m:mi>ρ</m:mi><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:msub><m:mi>E</m:mi><m:mi>λ</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">[</m:mo><m:mi>Tr</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi mathvariant="bold">Q</m:mi><m:mi>λ</m:mi></m:msub><m:mo>⁢</m:mo><m:mi>ρ</m:mi><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">]</m:mo></m:math></span>, with <span class="eq"><m:math><m:mi>E</m:mi></m:math></span> the expectation. Of course, the estimation of <span class="eq"><m:math><m:mi>Tr</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi mathvariant="bold">O</m:mi><m:mo>⁢</m:mo><m:mi>ρ</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> with a quorum cannot be perfect and an estimation variance <span class="eq"><m:math><m:msubsup><m:mi>σ</m:mi><m:mn>2</m:mn><m:mn>2</m:mn></m:msubsup></m:math></span> must also be considered and added, <span class="eq"><m:math><m:msup><m:mi>σ</m:mi><m:mn>2</m:mn></m:msup><m:mo>=</m:mo><m:msubsup><m:mi>σ</m:mi><m:mn>1</m:mn><m:mn>2</m:mn></m:msubsup><m:mo>+</m:mo><m:msubsup><m:mi>σ</m:mi><m:mn>2</m:mn><m:mn>2</m:mn></m:msubsup></m:math></span>.
</p><h3 id="id2251294">3.2 Estimating phase errors using coherent states and homodyne detection</h3><p>We now explain how the phase error rate can be estimated, in principle, using coherent states modulated in both quadratures and homodyne detection in all quadratures.</p><p>
It is clear that the knowledge of matrix elements of the average state <span class="eq"><m:math><m:mi>ρ</m:mi></m:math></span> gives the knowledge of <span class="eq"><m:math><m:mo>〈</m:mo><m:mi mathvariant="bold">O</m:mi><m:mo>〉</m:mo></m:math></span>. Let <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mn>0</m:mn></m:msub><m:mo>=</m:mo><m:mo>|</m:mo><m:mi>Ψ</m:mi><m:mo>〉</m:mo><m:mo>〈</m:mo><m:mi>Ψ</m:mi><m:mo>|</m:mo></m:math></span> be the state that Alice and Bob would share if the transmission was perfect. Since the <span class="eq"><m:math><m:mi mathvariant="sans-serif">a</m:mi></m:math></span> part of the system stays at Alice's station, we only need to learn about how the <span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span> part of the system is affected. In the prepare-and-measure picture, let <span class="eq"><m:math><m:mi>T</m:mi></m:math></span> be the completely positive (CP) map that maps the states sent by Alice onto the states received by Bob,
<span class="eq"><m:math><m:mo stretchy="false">(</m:mo><m:mi>Id</m:mi><m:mo>⊗</m:mo><m:mi>T</m:mi><m:mo stretchy="false">)</m:mo><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>ρ</m:mi><m:mn>0</m:mn></m:msub><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:mi>ρ</m:mi></m:math></span>.
In particular, let the coherent state <span class="eq"><m:math><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo>〉</m:mo><m:mo>〈</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo>|</m:mo></m:math></span> be mapped onto <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> and the (pseudo-)position state <span class="eq"><m:math><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>〉</m:mo><m:mo>〈</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>|</m:mo></m:math></span> be mapped onto <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>,</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo stretchy="false">)</m:mo></m:math></span>. The functions <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> and <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>,</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo stretchy="false">)</m:mo></m:math></span> are related by the following identity:
<table class="equation" id="id2251605"><tr><td class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mo>∝</m:mo><m:mo>∫</m:mo><m:mo>ⅆ</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>ⅆ</m:mo><m:msup><m:mi>x</m:mi><m:mrow><m:mo>′</m:mo><m:mo>′</m:mo></m:mrow></m:msup><m:msup><m:mi>e</m:mi><m:mrow><m:mo>-</m:mo><m:msup><m:mrow><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>-</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:mn>2</m:mn></m:msup><m:mo>/</m:mo><m:mn>4</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub><m:mo>-</m:mo><m:msup><m:mrow><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mrow><m:mo>′</m:mo><m:mo>′</m:mo></m:mrow></m:msup><m:mo>-</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:mn>2</m:mn></m:msup><m:mo>/</m:mo><m:mn>4</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub></m:mrow></m:msup><m:msup><m:mi>e</m:mi><m:mrow><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>-</m:mo><m:msup><m:mi>x</m:mi><m:mrow><m:mo>′</m:mo><m:mo>′</m:mo></m:mrow></m:msup><m:mo stretchy="false">)</m:mo><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo>/</m:mo><m:mn>2</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub></m:mrow></m:msup><m:mo>⁢</m:mo><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:msup><m:mi>x</m:mi><m:mrow><m:mo>′</m:mo><m:mo>′</m:mo></m:mrow></m:msup><m:mo stretchy="false">)</m:mo><m:mtext>,</m:mtext></m:math></td><td class="number">(3.2)</td></tr></table>
with <span class="eq"><m:math><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub></m:math></span> the variance of the vacuum fluctuations.
By setting <span class="eq"><m:math><m:mi>D</m:mi><m:mo>=</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>-</m:mo><m:msup><m:mi>x</m:mi><m:mrow><m:mo>′</m:mo><m:mo>′</m:mo></m:mrow></m:msup></m:math></span> and <span class="eq"><m:math><m:mi>S</m:mi><m:mo>=</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>+</m:mo><m:msup><m:mi>x</m:mi><m:mrow><m:mo>′</m:mo><m:mo>′</m:mo></m:mrow></m:msup><m:mo>-</m:mo><m:mn>2</m:mn><m:mo>⁢</m:mo><m:mi>x</m:mi></m:math></span>, we get:
<table class="equation" id="id2251905"><tr><td class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mo>∝</m:mo><m:mo>∫</m:mo><m:mo>ⅆ</m:mo><m:mi>D</m:mi><m:mo>ⅆ</m:mo><m:mi>S</m:mi><m:msup><m:mi>e</m:mi><m:mrow><m:mo>-</m:mo><m:msup><m:mi>S</m:mi><m:mn>2</m:mn></m:msup><m:mo>/</m:mo><m:mn>8</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub><m:mo>-</m:mo><m:msup><m:mi>D</m:mi><m:mn>2</m:mn></m:msup><m:mo>/</m:mo><m:mn>8</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>D</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo>/</m:mo><m:mn>2</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub></m:mrow></m:msup><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>S</m:mi><m:mo>+</m:mo><m:mi>D</m:mi><m:mo>,</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>S</m:mi><m:mo>-</m:mo><m:mi>D</m:mi><m:mo stretchy="false">)</m:mo><m:mtext>,</m:mtext></m:math></td><td class="number">(3.3)</td></tr></table>
which shows that <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>,</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo stretchy="false">)</m:mo></m:math></span> is integrated with an invertible kernel (Gaussian convolution in <span class="eq"><m:math><m:mi>S</m:mi></m:math></span>, multiplication by
<span class="eq"><m:math><m:msup><m:mi>e</m:mi><m:mrow><m:mo>-</m:mo><m:msup><m:mi>D</m:mi><m:mn>2</m:mn></m:msup><m:mo>/</m:mo><m:mn>8</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub></m:mrow></m:msup></m:math></span>, and Fourier transform in <span class="eq"><m:math><m:mi>D</m:mi></m:math></span>). So in principle, any different CP map <span class="eq"><m:math><m:msup><m:mi>T</m:mi><m:mo>′</m:mo></m:msup><m:mo>≠</m:mo><m:mi>T</m:mi></m:math></span> implies a different effect on coherent states, <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mo>≠</m:mo><m:msub><m:mi>ρ</m:mi><m:msup><m:mi>T</m:mi><m:mo>′</m:mo></m:msup></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>. The modulation of coherent states in both quadratures is thus crucial for this implication being possible.</p><p>By inspecting Eq. (<span class="ref"><a href="#id2251905">3.3</a></span>), it seems that due to the factors <span class="eq"><m:math><m:msup><m:mi>e</m:mi><m:mrow><m:mo>-</m:mo><m:msup><m:mi>S</m:mi><m:mn>2</m:mn></m:msup><m:mo>/</m:mo><m:mn>8</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub></m:mrow></m:msup></m:math></span> and <span class="eq"><m:math><m:msup><m:mi>e</m:mi><m:mrow><m:mo>-</m:mo><m:msup><m:mi>D</m:mi><m:mn>2</m:mn></m:msup><m:mo>/</m:mo><m:mn>8</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub></m:mrow></m:msup></m:math></span>, two different CP-maps <span class="eq"><m:math><m:mi>T</m:mi></m:math></span> and <span class="eq"><m:math><m:msup><m:mi>T</m:mi><m:mo>′</m:mo></m:msup></m:math></span> may make <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> and <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:msup><m:mi>T</m:mi><m:mo>′</m:mo></m:msup></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> only vanishingly different. It thus seems unlikely that Eq. (<span class="ref"><a href="#id2251905">3.3</a></span>) should allow us to extract the coefficients
<span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>S</m:mi><m:mo>+</m:mo><m:mi>D</m:mi><m:mo>,</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>S</m:mi><m:mo>-</m:mo><m:mi>D</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>. However, assuming that <span class="eq"><m:math><m:mi>T</m:mi></m:math></span> depends only on a finite number of parameters, a variation of these parameters
will induce a measurable variation of <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>. We will now discuss why it is reasonable to make such an assumption.</p><p>Due to the finite variance of the modulation of coherent states, the probability of emission of a large number of photons vanishes—this intuitively indicates that we only need to consider the description of <span class="eq"><m:math><m:mi>T</m:mi></m:math></span> for a bounded number of emitted photons. More precisely, one can consider the emission of <span class="eq"><m:math><m:mi>w</m:mi></m:math></span> joint copies of the state <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mrow><m:mn>0</m:mn><m:mi mathvariant="sans-serif">b</m:mi></m:mrow></m:msub><m:mo>=</m:mo><m:msub><m:mi>Tr</m:mi><m:mi mathvariant="sans-serif">a</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>ρ</m:mi><m:mn>0</m:mn></m:msub><m:mo stretchy="false">)</m:mo></m:math></span>. For <span class="eq"><m:math><m:mi>w</m:mi></m:math></span> sufficiently large <span class="eq"><m:math><m:msubsup><m:mi>ρ</m:mi><m:mrow><m:mn>0</m:mn><m:mi mathvariant="sans-serif">b</m:mi></m:mrow><m:mrow><m:mo>⊗</m:mo><m:mi>w</m:mi></m:mrow></m:msubsup></m:math></span> can be represented in the typical subspace <span class="eq"><m:math><m:msub><m:mi>Γ</m:mi><m:mi>δ</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>ρ</m:mi><m:mrow><m:mn>0</m:mn><m:mi mathvariant="sans-serif">b</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span> of dimension not greater than <span class="eq"><m:math><m:msup><m:mn>2</m:mn><m:mrow><m:mi>w</m:mi><m:mo>⁢</m:mo><m:mo stretchy="false">(</m:mo><m:mi>H</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>ρ</m:mi><m:mrow><m:mn>0</m:mn><m:mi mathvariant="sans-serif">b</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo><m:mo>+</m:mo><m:mi>δ</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msup></m:math></span>, for any <span class="eq"><m:math><m:mi>δ</m:mi><m:mo>&gt;</m:mo><m:mn>0</m:mn></m:math></span> <span class="bib-cite">[<a href="#id2270585">20</a>]</span>, where <span class="eq"><m:math><m:mi>H</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>ρ</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> is the Von Neumann entropy of a state <span class="eq"><m:math><m:mi>ρ</m:mi></m:math></span>. The probability mass of <span class="eq"><m:math><m:msubsup><m:mi>ρ</m:mi><m:mrow><m:mn>0</m:mn><m:mi mathvariant="sans-serif">b</m:mi></m:mrow><m:mrow><m:mo>⊗</m:mo><m:mi>w</m:mi></m:mrow></m:msubsup></m:math></span> outside the typical subspace can be made arbitrarily small and does not depend on the eavesdropping strategy. This means that the support for the input of <span class="eq"><m:math><m:mi>T</m:mi></m:math></span> has finite dimension, up to an arbitrarily small deviation.</p><p>The number of photons received by Bob can also be upper bounded. Alice and Bob can first assume that no more than <span class="eq"><m:math><m:msub><m:mi>n</m:mi><m:mtext>max</m:mtext></m:msub></m:math></span> photons are received. This fact may depend on a malicious eavesdropper, so Bob has to do hypothesis testing. The test comes down to estimating <span class="eq"><m:math><m:mo>〈</m:mo><m:mi>Π</m:mi><m:mo>〉</m:mo></m:math></span> with <span class="eq"><m:math><m:mi>Π</m:mi><m:mo>=</m:mo><m:msub><m:mo>∑</m:mo><m:mrow><m:mi>n</m:mi><m:mo>&gt;</m:mo><m:msub><m:mi>n</m:mi><m:mtext>max</m:mtext></m:msub></m:mrow></m:msub><m:mo>|</m:mo><m:mi>n</m:mi><m:mo>〉</m:mo><m:mo>〈</m:mo><m:mi>n</m:mi><m:mo>|</m:mo></m:math></span>. If the threshold is well chosen so that <span class="eq"><m:math><m:mi>n</m:mi><m:mo>&gt;</m:mo><m:msub><m:mi>n</m:mi><m:mtext>max</m:mtext></m:msub></m:math></span> never occurs, we can apply the central limit theorem and upper bound the probability that <span class="eq"><m:math><m:mo>〈</m:mo><m:mi>Π</m:mi><m:mo>〉</m:mo><m:mo>&gt;</m:mo><m:mi>ε</m:mi></m:math></span> for any chosen <span class="eq"><m:math><m:mi>ε</m:mi><m:mo>&gt;</m:mo><m:mn>0</m:mn></m:math></span>. The positivity of the density matrices implies that the off-diagonal coefficients are also bounded. We can thus now express <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> as <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:msub><m:mo>∑</m:mo><m:mrow><m:mi>n</m:mi><m:mo>,</m:mo><m:msup><m:mi>n</m:mi><m:mo>′</m:mo></m:msup><m:mo>≤</m:mo><m:msub><m:mi>n</m:mi><m:mtext>max</m:mtext></m:msub></m:mrow></m:msub><m:msub><m:mi>ρ</m:mi><m:mi>T</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo>,</m:mo><m:mi>n</m:mi><m:mo>,</m:mo><m:msup><m:mi>n</m:mi><m:mo>′</m:mo></m:msup><m:mo stretchy="false">)</m:mo><m:mo>⁢</m:mo><m:mo>|</m:mo><m:mi>n</m:mi><m:mo>〉</m:mo><m:mo>〈</m:mo><m:msup><m:mi>n</m:mi><m:mo>′</m:mo></m:msup><m:mo>|</m:mo></m:math></span>. Note that the test can be implemented either by explicitly measuring the intensity of the beam (therefore requiring an additional photodetector) or by exploiting the correlation between the high intensity of the beam and the high absolute values obtained when doing homodyne measurements in all directions.</p><p>Finally, the estimation of the coefficient of <span class="eq"><m:math><m:mo>|</m:mo><m:mi>n</m:mi><m:mo>〉</m:mo><m:mo>〈</m:mo><m:msup><m:mi>n</m:mi><m:mo>′</m:mo></m:msup><m:mo>|</m:mo></m:math></span> can be done with arbitrarily small statistical error using homodyne detection in all directions <span class="bib-cite">[<a href="#id2270532">19</a>, <a href="#id2270632">21</a>]</span>. This is achieved by considering the quorum of operators <span class="eq"><m:math><m:msub><m:mrow><m:mo stretchy="false">(</m:mo><m:msub><m:mi mathvariant="bold">x</m:mi><m:mi>θ</m:mi></m:msub><m:mo stretchy="false">)</m:mo></m:mrow><m:mrow><m:mn>0</m:mn><m:mo>≤</m:mo><m:mi>θ</m:mi><m:mo>&lt;</m:mo><m:mn>2</m:mn><m:mo>⁢</m:mo><m:mi>π</m:mi></m:mrow></m:msub></m:math></span>, where <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">x</m:mi><m:mi>θ</m:mi></m:msub><m:mo>=</m:mo><m:mi>cos</m:mi><m:mo>⁡</m:mo><m:mi>θ</m:mi><m:mspace width="0.222em"/><m:mi mathvariant="bold">x</m:mi><m:mo>+</m:mo><m:mi>sin</m:mi><m:mo>⁡</m:mo><m:mi>θ</m:mi><m:mspace width="0.222em"/><m:mi mathvariant="bold">p</m:mi></m:math></span> denotes the amplitude of the quadrature in direction <span class="eq"><m:math><m:mi>θ</m:mi></m:math></span>. Considering a finite combination of arbitrarily small statistical errors on parameters also gives arbitrarily small overall statistical error on the phase error rate.

</p><h2 id="id2253094">4 Encoding of Multiple QuBits in an Oscillator</h2><p>Reconciliation and privacy amplification are integral parts of the prepare-and-measure protocols derived from entanglement purification protocols. In our case, we wish to derive a prepare-and-measure protocol with sliced error correction (SEC) <span class="bib-cite">[<a href="#id2270414">17</a>]</span> as reconciliation, which allows us to obtain a higher secret key rate and a better resistance to losses than in <span class="bib-cite">[<a href="#id2270285">14</a>]</span>. We therefore need to describe an entanglement purification procedure that reduces to SEC when the corresponding prepare-and-measure protocol is derived. An overview of SEC is proposed next.</p><h3 id="id2253132">4.1 Sliced error correction with invertible mappings</h3><p>We here recall the main principles of SEC in a form that is slightly different from the presentation in <span class="bib-cite">[<a href="#id2270414">17</a>]</span>. To suit our needs, we here describe SEC in terms of invertible functions giving the slices and the estimators—the invertibility property will be required when we generalize SEC to entanglement purification. Also, from the generality of <span class="bib-cite">[<a href="#id2270414">17</a>]</span>, two parameters are fixed here: Binary error correction is operated by sending syndromes of classical linear error-correcting codes (ECC), and we momentarily restrict ourselves to the case of scalar values.</p><p>Suppose Alice and Bob have <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> pairs of correlated random variables <span class="eq"><m:math><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mn>1</m:mn></m:msub><m:mo>,</m:mo><m:msubsup><m:mi>X</m:mi><m:mn>1</m:mn><m:mo>′</m:mo></m:msubsup><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:mo>…</m:mo><m:mo>,</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mi>l</m:mi></m:msub><m:mo>,</m:mo><m:msubsup><m:mi>X</m:mi><m:mi>l</m:mi><m:mo>′</m:mo></m:msubsup><m:mo stretchy="false">)</m:mo></m:math></span>, with <span class="eq"><m:math><m:msub><m:mi>X</m:mi><m:mi>i</m:mi></m:msub><m:mo>,</m:mo><m:msubsup><m:mi>X</m:mi><m:mi>i</m:mi><m:mo>′</m:mo></m:msubsup><m:mo>∈</m:mo><m:mi mathvariant="bold">R</m:mi></m:math></span>, <span class="eq"><m:math><m:mi>i</m:mi><m:mo>=</m:mo><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:math></span>, from which they intend to extract common bits.</p><p>First, Alice wishes to convert each of her variables <span class="eq"><m:math><m:mi>X</m:mi></m:math></span> into <span class="eq"><m:math><m:mi>m</m:mi></m:math></span> bits and thereby defines <span class="eq"><m:math><m:mi>m</m:mi></m:math></span> binary functions: <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:mspace width="0.222em"/><m:mo>…</m:mo><m:mo>,</m:mo><m:msub><m:mi>S</m:mi><m:mi>m</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>. To make the mapping invertible, she also defines a function <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> such that mapping from <span class="eq"><m:math><m:mi>X</m:mi></m:math></span> to the vector
<span class="eq"><m:math><m:mo stretchy="false">(</m:mo><m:mrow><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>,</m:mo><m:mrow><m:msub><m:mi>S</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:mo stretchy="false">)</m:mo></m:math></span> is bijective.
As a convention, the range of <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> is <span class="eq"><m:math><m:mo stretchy="false">[</m:mo><m:mn>0</m:mn><m:mo>;</m:mo><m:mn>1</m:mn><m:mo stretchy="false">]</m:mo></m:math></span>. The mapping from <span class="eq"><m:math><m:mi mathvariant="bold">R</m:mi></m:math></span> to <span class="eq"><m:math><m:mrow><m:mo stretchy="false">[</m:mo><m:mn>0</m:mn><m:mo>;</m:mo><m:mn>1</m:mn><m:mo stretchy="false">]</m:mo></m:mrow><m:mo>×</m:mo><m:msup><m:mrow><m:mo stretchy="false">{</m:mo><m:mn>0</m:mn><m:mo>,</m:mo><m:mn>1</m:mn><m:mo stretchy="false">}</m:mo></m:mrow><m:mi>m</m:mi></m:msup></m:math></span>: <table class="equation" id="id2253543"><tr><td class="eq"><m:math><m:mi>x</m:mi><m:mo>→</m:mo><m:mo stretchy="false">(</m:mo><m:mrow><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>,</m:mo><m:mrow><m:msub><m:mi>S</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:mo stretchy="false">)</m:mo></m:math></td><td class="number">(4.1)</td></tr></table>
is collectively denoted as <span class="eq"><m:math><m:mi mathvariant="script">S</m:mi></m:math></span>.</p><p>Concretely, the functions <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> implicitly cut the real line into intervals (see <span class="bib-cite">[<a href="#id2270414">17</a>]</span> for more details), whereas <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> indicates where to find <span class="eq"><m:math><m:mi>X</m:mi></m:math></span> within a given interval.</p><p>Then, we can assemble the bits produced by the <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> random variables <span class="eq"><m:math><m:msub><m:mi>X</m:mi><m:mn>1</m:mn></m:msub><m:mo>…</m:mo><m:msub><m:mi>X</m:mi><m:mi>l</m:mi></m:msub></m:math></span> into <span class="eq"><m:math><m:mi>m</m:mi></m:math></span> <span class="eq"><m:math><m:mi>l</m:mi></m:math></span>-bit vectors. To each bit vector ("slice") <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:mrow><m:mo stretchy="false">(</m:mo><m:msub><m:mi>S</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mn>1</m:mn></m:msub><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:mo>…</m:mo><m:mo>,</m:mo><m:msub><m:mi>S</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mi>l</m:mi></m:msub><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo></m:mrow></m:math></span>, is associated an ECC that Alice and Bob agreed upon. To proceed with the correction, Alice sends the syndrome
<span class="eq"><m:math><m:msubsup><m:mi>ξ</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo>=</m:mo><m:msubsup><m:mi>H</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo>⁢</m:mo><m:msub><m:mi>S</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span> to Bob over the public authenticated channel, where <span class="eq"><m:math><m:msubsup><m:mi>H</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup></m:math></span> is the <span class="eq"><m:math><m:msubsup><m:mi>l</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo>×</m:mo><m:mi>l</m:mi></m:math></span> parity check matrix of the ECC associated to slice <span class="eq"><m:math><m:mi>i</m:mi></m:math></span>. Alice also sends <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span>.</p><p>
Bob would like to recover
<span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span> from his knowledge of <span class="eq"><m:math><m:msubsup><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span>, <span class="eq"><m:math><m:msubsup><m:mi>ξ</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mtext>b</m:mtext></m:msubsup></m:math></span> and <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span>.
To do so, he also converts each of his variables <span class="eq"><m:math><m:msubsup><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span> into <span class="eq"><m:math><m:mi>m</m:mi></m:math></span> bits, but he does so in a consecutive manner. He tries to produce bits that are best correlated to Alice's and takes advantage of the corrected bits of slices <span class="eq"><m:math><m:mi>j</m:mi><m:mo>&lt;</m:mo><m:mi>i</m:mi></m:math></span> before trying to estimate the bits of slice <span class="eq"><m:math><m:mi>i</m:mi></m:math></span>.
In particular, to produce bits that are best correlated to Alice's first slice <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span>, he uses a function <span class="eq"><m:math><m:msub><m:mi>E</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>X</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo></m:math></span>, which gives his best estimate on Alice's corresponding bit <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> given the known correlations between <span class="eq"><m:math><m:mi>X</m:mi></m:math></span> and <span class="eq"><m:math><m:msup><m:mi>X</m:mi><m:mo>′</m:mo></m:msup></m:math></span>.
By applying the function <span class="eq"><m:math><m:msub><m:mi>E</m:mi><m:mn>1</m:mn></m:msub></m:math></span> on all the variables <span class="eq"><m:math><m:msubsup><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span> and <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span>, Bob is able to construct a string of <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> bits that is equal to Alice's up to some error rate <span class="eq"><m:math><m:msubsup><m:mi>e</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup></m:math></span>.
Given the knowledge of <span class="eq"><m:math><m:msubsup><m:mi>ξ</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup></m:math></span> and assuming the adequacy of the ECC, Bob has enough information to determine <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span> with high probability.
Then, for slice <span class="eq"><m:math><m:mi>i</m:mi><m:mo>&gt;</m:mo><m:mn>1</m:mn></m:math></span>, he estimates <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span> using the estimator function <span class="eq"><m:math><m:msub><m:mi>E</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>X</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:msub><m:mi>β</m:mi><m:mn>1</m:mn></m:msub><m:mo>,</m:mo><m:mo>…</m:mo><m:mo>,</m:mo><m:msub><m:mi>β</m:mi><m:mrow><m:mi>i</m:mi><m:mo>-</m:mo><m:mn>1</m:mn></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span>, where <span class="eq"><m:math><m:msub><m:mi>β</m:mi><m:mi>j</m:mi></m:msub></m:math></span> is the random variable indicating Bob's knowledge of <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mi>j</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>, so that <span class="eq"><m:math><m:msub><m:mi>β</m:mi><m:mi>j</m:mi></m:msub><m:mo>=</m:mo><m:msub><m:mi>S</m:mi><m:mi>j</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> with arbitrarily high probability. (Note that the estimators can also be written as jointly working on <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> samples at once: <span class="eq"><m:math><m:msub><m:mi>E</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msubsup><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:msubsup><m:mi>ξ</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup><m:mo>,</m:mo><m:mo>…</m:mo><m:mo>,</m:mo><m:msubsup><m:mi>ξ</m:mi><m:mrow><m:mi>i</m:mi><m:mo>-</m:mo><m:mn>1</m:mn></m:mrow><m:mtext>b</m:mtext></m:msubsup><m:mo stretchy="false">)</m:mo></m:math></span>,
but we will preferably use the previous notation for its simplicity since, besides the ECC decoding, all the operations are done on each variable <span class="eq"><m:math><m:mi>X</m:mi></m:math></span> or <span class="eq"><m:math><m:msup><m:mi>X</m:mi><m:mo>′</m:mo></m:msup></m:math></span> independently.)
</p><p>
We also need a supplementary function to ensure that the process on Bob's side is described using bijective functions: <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>X</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:msub><m:mi>β</m:mi><m:mn>1</m:mn></m:msub><m:mo>,</m:mo><m:mo>…</m:mo><m:mo>,</m:mo><m:msub><m:mi>β</m:mi><m:mi>m</m:mi></m:msub><m:mo stretchy="false">)</m:mo></m:math></span> (or jointly <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msubsup><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:msubsup><m:mi>ξ</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup><m:mo>,</m:mo><m:mo>…</m:mo><m:mo>,</m:mo><m:msubsup><m:mi>ξ</m:mi><m:mi>m</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo stretchy="false">)</m:mo></m:math></span>). As a convention, the range of <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo></m:math></span> is <span class="eq"><m:math><m:mo stretchy="false">[</m:mo><m:mn>0</m:mn><m:mo>;</m:mo><m:mn>1</m:mn><m:mo stretchy="false">]</m:mo></m:math></span>. <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> is chosen so that the mapping <span class="eq"><m:math><m:mi mathvariant="script">E</m:mi></m:math></span> defined below is invertible,
<table class="equation" id="id2254845"><tr><td class="eq"><m:math><m:mi mathvariant="script">E</m:mi><m:mtext>: </m:mtext><m:mrow><m:mo stretchy="false">(</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:mrow><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:mrow><m:mo>,</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>→</m:mo><m:mo stretchy="false">(</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo>,</m:mo><m:mrow><m:msub><m:mi>E</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>,</m:mo><m:mo>…</m:mo><m:mo>,</m:mo><m:mrow><m:msub><m:mi>E</m:mi><m:mi>m</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi><m:mo>-</m:mo><m:mn>1</m:mn></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>,</m:mo><m:mrow><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo stretchy="false">)</m:mo></m:mrow><m:mo stretchy="false">)</m:mo><m:mtext>.</m:mtext></m:math></td><td class="number">(4.2)</td></tr></table></p><p>Similarly to <span class="eq"><m:math><m:mi mathvariant="script">S</m:mi></m:math></span>, the functions <span class="eq"><m:math><m:msub><m:mi>E</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow></m:msub></m:math></span> of <span class="eq"><m:math><m:mi mathvariant="script">E</m:mi></m:math></span> cut the real line into intervals. However, these intervals are adapted as a function of the information sent by Alice, so as to estimate Alice's bits more reliably. Like for <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover></m:math></span>, the function <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> indicates where to find <span class="eq"><m:math><m:msup><m:mi>X</m:mi><m:mo>′</m:mo></m:msup></m:math></span> within an interval.</p><p>
The mapping <span class="eq"><m:math><m:mi mathvariant="script">S</m:mi></m:math></span> summarizes Alice's process of conversion of her real values <span class="eq"><m:math><m:mi>X</m:mi></m:math></span> into <span class="eq"><m:math><m:mi>m</m:mi></m:math></span> bits (plus a continuous component). The mapping <span class="eq"><m:math><m:mi mathvariant="script">E</m:mi></m:math></span> represents the bits (and a continuous component) produced by Bob from his real values <span class="eq"><m:math><m:msup><m:mi>X</m:mi><m:mo>′</m:mo></m:msup></m:math></span> and his knowledge of <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> and of the syndromes
<span class="eq"><m:math><m:msubsup><m:mi>ξ</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mtext>b</m:mtext></m:msubsup></m:math></span>. The bits produced by the functions <span class="eq"><m:math><m:msub><m:mi>E</m:mi><m:mi>i</m:mi></m:msub></m:math></span> are not yet corrected by the ECC, even though they take as input the corrected values of the previous slices <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mi>j</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>, <span class="eq"><m:math><m:mi>j</m:mi><m:mo>&lt;</m:mo><m:mi>i</m:mi></m:math></span>.
The description of the mapping <span class="eq"><m:math><m:mi mathvariant="script">E</m:mi></m:math></span> with the bits prior to ECC correction allows us to easily express the bit error rate between Alice's slices and Bob's estimators and thereby to deduce the size of the parity matrices of the ECCs needed for the binary correction. Simply, we define <span class="eq"><m:math><m:msubsup><m:mi>e</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo>=</m:mo><m:mi>Pr</m:mi><m:mo stretchy="false">[</m:mo><m:mrow><m:msub><m:mi>S</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo><m:mo>≠</m:mo><m:msub><m:mi>E</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>X</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:msub><m:mi>S</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>i</m:mi><m:mo>-</m:mo><m:mn>1</m:mn></m:mrow></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo></m:mrow><m:mo stretchy="false">]</m:mo></m:math></span>. As the block size <span class="eq"><m:math><m:mi>l</m:mi><m:mo>→</m:mo><m:mo>∞</m:mo></m:math></span>, there exist ECCs with size <span class="eq"><m:math><m:msubsup><m:mi>l</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo>→</m:mo><m:mi>l</m:mi><m:mo>⁢</m:mo><m:mi>h</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msubsup><m:mi>e</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo stretchy="false">)</m:mo></m:math></span> and arbitrarily low probability of decoding error. The number of common (but not necessarily secret) bits produced by SEC is therefore asymptotically equal to <span class="eq"><m:math><m:mi>H</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>S</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo><m:mo>-</m:mo><m:mrow><m:msubsup><m:mo>∑</m:mo><m:mrow><m:mi>i</m:mi><m:mo>=</m:mo><m:mn>1</m:mn></m:mrow><m:mi>m</m:mi></m:msubsup><m:mi>h</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msubsup><m:mi>e</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo stretchy="false">)</m:mo></m:mrow></m:math></span> per sample <span class="bib-cite">[<a href="#id2270414">17</a>]</span>.
</p><p>The generalization of the SEC to a quantum entanglement purification protocol is examined next.</p><h3 id="id2255560">4.2 Quantum sliced error correction</h3><p>From classical binary error correcting codes, one can construct CSS quantum codes and use them to extract EPR pairs from noisy qubit pairs. We will now show that, similarly, from SEC, it is possible to construct an encoding and decoding procedure, which when applied on entangled quantum oscillator systems, also allows to extract pure EPR pairs. Such a purification protocol is formal, as it would of course be very difficult to implement in practice.</p><p>The purification uses a few quantum registers, which we now list. Alice's system <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:math></span> is split into <span class="eq"><m:math><m:mi>m</m:mi></m:math></span> qubit systems <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow></m:msub></m:math></span> and a continuous register <span class="eq"><m:math><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:math></span>. On Bob's side, the system <span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span> is split into <span class="eq"><m:math><m:mi>m</m:mi></m:math></span> qubit systems <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow></m:msub></m:math></span> and a continuous register <span class="eq"><m:math><m:mover accent="true"><m:mi mathvariant="sans-serif">e</m:mi><m:mo>¯</m:mo></m:mover></m:math></span>. He also needs <span class="eq"><m:math><m:mi>m</m:mi></m:math></span> qubit registers <span class="eq"><m:math><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span> for temporary storage. All these registers must of course be understood per exchanged sample: As Alice generates <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> copies of the state <span class="eq"><m:math><m:mo>|</m:mo><m:mi>Ψ</m:mi><m:mo>〉</m:mo></m:math></span>, the legitimate parties use <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> instances of the registers listed above.</p><p>The usual bit-flip and phase-flip operators <span class="eq"><m:math><m:mi mathvariant="bold">X</m:mi></m:math></span> and <span class="eq"><m:math><m:mi mathvariant="bold">Z</m:mi></m:math></span>, resp., can be defined as acting on a specific qubit register among the systems <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">s</m:mi><m:mi>i</m:mi></m:msub></m:math></span> and <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mi>i</m:mi></m:msub></m:math></span>. E.g., <span class="eq"><m:math><m:msub><m:mi mathvariant="bold">Z</m:mi><m:msub><m:mi mathvariant="sans-serif">s</m:mi><m:mi>i</m:mi></m:msub></m:msub></m:math></span> is defined as acting on <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">s</m:mi><m:mi>i</m:mi></m:msub></m:math></span> only. These operators are used by Alice and Bob to construct the CSS codes that produce entangled qubits, which are in turn used to produce EPR pairs in the registers <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">s</m:mi><m:mi>i</m:mi></m:msub><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mi>i</m:mi></m:msub></m:math></span> for <span class="eq"><m:math><m:mi>i</m:mi><m:mo>=</m:mo><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:math></span>. Since each CSS code operates in its own register pair, the action of one does not interfere with the action of the other. It is thus possible to extract
more than one EPR pair <span class="eq"><m:math><m:mo>|</m:mo><m:msup><m:mi>φ</m:mi><m:mo>+</m:mo></m:msup><m:mo>〉</m:mo></m:math></span> per state <span class="eq"><m:math><m:mo>|</m:mo><m:mi>Ψ</m:mi><m:mo>〉</m:mo></m:math></span>. If asymptotically efficient binary codes are used, the rate or EPR pairs produced is <span class="eq"><m:math><m:mi>R</m:mi><m:mo>=</m:mo><m:msub><m:mo>∑</m:mo><m:mi>i</m:mi></m:msub><m:mo stretchy="false">(</m:mo><m:mn>1</m:mn><m:mo>-</m:mo><m:mi>h</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msubsup><m:mi>e</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo stretchy="false">)</m:mo><m:mo>-</m:mo><m:mi>h</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msubsup><m:mi>e</m:mi><m:mi>i</m:mi><m:mtext>p</m:mtext></m:msubsup><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo></m:math></span>, where <span class="eq"><m:math><m:msubsup><m:mi>e</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup></m:math></span> (resp. <span class="eq"><m:math><m:msubsup><m:mi>e</m:mi><m:mi>i</m:mi><m:mtext>p</m:mtext></m:msubsup></m:math></span>) indicates the bit error rate (resp. the phase error rate) <span class="bib-cite">[<a href="#id2270238">13</a>]</span>.</p><p>The process that defines the content of the registers is described next.</p><h4 id="id2255966">4.2.1 The mappings <span class="eq"><m:math><m:mi mathvariant="script">QS</m:mi></m:math></span> and <span class="eq"><m:math><m:mi mathvariant="script">QE</m:mi></m:math></span></h4><p>First, we define the unitary transformation <span class="eq"><m:math><m:mi mathvariant="script">QS</m:mi></m:math></span>: <span class="eq"><m:math><m:msup><m:mi>L</m:mi><m:mn>2</m:mn></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi mathvariant="bold">R</m:mi><m:mo stretchy="false">)</m:mo><m:mo>→</m:mo><m:msup><m:mi>L</m:mi><m:mn>2</m:mn></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mo stretchy="false">[</m:mo><m:mn>0</m:mn><m:mo>;</m:mo><m:mn>1</m:mn><m:mo stretchy="false">]</m:mo><m:mo stretchy="false">)</m:mo><m:mo>⊗</m:mo><m:msup><m:mi mathvariant="script">H</m:mi><m:mrow><m:mo>⊗</m:mo><m:mi>m</m:mi></m:mrow></m:msup></m:math></span> by its application on the basis of quadrature eigenstates:
<table class="equation" id="id2256075"><tr><td class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:msub><m:mo>→</m:mo><m:mrow><m:mi>σ</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>⁢</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:mo>⊗</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>S</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn></m:msub></m:msub><m:mo>⊗</m:mo><m:mo>…</m:mo><m:mo>⊗</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>S</m:mi><m:mi>m</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">s</m:mi><m:mi>m</m:mi></m:msub></m:msub><m:mtext>.</m:mtext></m:math></td><td class="number">(4.3)</td></tr></table>
The states <span class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:msub></m:math></span>, <span class="eq"><m:math><m:mn>0</m:mn><m:mo>≤</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>≤</m:mo><m:mn>1</m:mn></m:math></span>, form an orthogonal basis of <span class="eq"><m:math><m:msup><m:mi>L</m:mi><m:mn>2</m:mn></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mo stretchy="false">[</m:mo><m:mn>0</m:mn><m:mo>;</m:mo><m:mn>1</m:mn><m:mo stretchy="false">]</m:mo><m:mo stretchy="false">)</m:mo></m:math></span>, <span class="eq"><m:math><m:mi>σ</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:msup><m:mrow><m:mo stretchy="false">(</m:mo><m:msub><m:mo>ⅆ</m:mo><m:mi>x</m:mi></m:msub><m:mo>⁡</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo></m:mrow><m:mrow><m:mo>-</m:mo><m:mn>1</m:mn><m:mo>/</m:mo><m:mn>2</m:mn></m:mrow></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo></m:math></span> is a normalization function, and <span class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>s</m:mi><m:mi>i</m:mi></m:msub><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">s</m:mi><m:mi>i</m:mi></m:msub></m:msub></m:math></span>, <span class="eq"><m:math><m:msub><m:mi>s</m:mi><m:mi>i</m:mi></m:msub><m:mo>∈</m:mo><m:mo stretchy="false">{</m:mo><m:mn>0</m:mn><m:mo>,</m:mo><m:mn>1</m:mn><m:mo stretchy="false">}</m:mo></m:math></span>, denotes the canonical basis of <span class="eq"><m:math><m:mi mathvariant="script">H</m:mi></m:math></span>, the Hilbert space of a qubit. As a convention, the system <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">s</m:mi><m:mi>i</m:mi></m:msub></m:math></span> is called slice <span class="eq"><m:math><m:mi>i</m:mi></m:math></span>. The transformation <span class="eq"><m:math><m:mi mathvariant="script">QS</m:mi></m:math></span> is depicted in Fig. <span class="ref"><a href="#id2256485">4.1</a></span>.</p><p id="id2256485" class="figure"><img src="QS.png" type="image/png"/></p><p class="caption">Figure 4.1: Schematic description of <span class="eq"><m:math><m:mi mathvariant="script">QS</m:mi></m:math></span></p><p>
For each slice <span class="eq"><m:math><m:mi>i</m:mi></m:math></span>, Alice and Bob agree on a CSS code, defined by its parity matrices <span class="eq"><m:math><m:msubsup><m:mi>H</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup></m:math></span> for bit error correction and <span class="eq"><m:math><m:msubsup><m:mi>H</m:mi><m:mi>i</m:mi><m:mtext>p</m:mtext></m:msubsup></m:math></span> for phase error correction. For the entanglement purification, let us assume that Alice computes the syndromes of the CSS code with a quantum circuit. For each slice, she produces <span class="eq"><m:math><m:msubsup><m:mi>l</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup></m:math></span> qubits in the state <span class="eq"><m:math><m:mo>|</m:mo><m:msubsup><m:mi>ξ</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup><m:mo>〉</m:mo></m:math></span> and <span class="eq"><m:math><m:msubsup><m:mi>l</m:mi><m:mi>i</m:mi><m:mtext>p</m:mtext></m:msubsup></m:math></span> qubits in the state <span class="eq"><m:math><m:mo>|</m:mo><m:msubsup><m:mi>ξ</m:mi><m:mi>i</m:mi><m:mtext>p</m:mtext></m:msubsup><m:mo>〉</m:mo></m:math></span> that she sends to Bob over a perfect quantum channel, so that the syndromes are received without any distortion. In the entanglement purification picture, the syndromes can be transmitted over a non-perfect channel if they are encoded using appropriate error correcting codes. Also, after reduction to a prepare-and-measure protocol, this perfect transmission is actually done over the public authenticated channel. Alice also sends the <span class="eq"><m:math><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> system to Bob.
</p><p>
Then, the slice estimators are defined as the unitary transformation <span class="eq"><m:math><m:mi mathvariant="script">QE</m:mi></m:math></span> from <span class="eq"><m:math><m:mrow><m:msup><m:mi>L</m:mi><m:mn>2</m:mn></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mo stretchy="false">[</m:mo><m:mn>0</m:mn><m:mo>;</m:mo><m:mn>1</m:mn><m:mo stretchy="false">]</m:mo><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>⊗</m:mo><m:msup><m:mi mathvariant="script">H</m:mi><m:mrow><m:mo>⊗</m:mo><m:mi>m</m:mi></m:mrow></m:msup><m:mo>⊗</m:mo><m:mrow><m:msup><m:mi>L</m:mi><m:mn>2</m:mn></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi mathvariant="bold">R</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:math></span> to <span class="eq"><m:math><m:mrow><m:msup><m:mi>L</m:mi><m:mn>2</m:mn></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mo stretchy="false">[</m:mo><m:mn>0</m:mn><m:mo>;</m:mo><m:mn>1</m:mn><m:mo stretchy="false">]</m:mo><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>⊗</m:mo><m:msup><m:mi mathvariant="script">H</m:mi><m:mrow><m:mo>⊗</m:mo><m:mi>m</m:mi></m:mrow></m:msup><m:mo>⊗</m:mo><m:msup><m:mi mathvariant="script">H</m:mi><m:mrow><m:mo>⊗</m:mo><m:mi>m</m:mi></m:mrow></m:msup><m:mo>⊗</m:mo><m:mrow><m:msup><m:mi>L</m:mi><m:mn>2</m:mn></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mo stretchy="false">[</m:mo><m:mn>0</m:mn><m:mo>;</m:mo><m:mn>1</m:mn><m:mo stretchy="false">]</m:mo><m:mo stretchy="false">)</m:mo></m:mrow></m:math></span>:
<table class="equation" id="id2256815"><tr><td class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:mo>⁢</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:msub><m:mo>⁢</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mi mathvariant="sans-serif">b</m:mi></m:msub><m:mo>→</m:mo><m:mrow><m:mi>ε</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo stretchy="false">)</m:mo></m:mrow><m:mo>⁢</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:mo>⁢</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:msub><m:mo>⁢</m:mo>

<m:msubsup><m:mo largeop="true">⊗</m:mo><m:mrow><m:mi>i</m:mi><m:mo>=</m:mo><m:mn>1</m:mn></m:mrow><m:mi>m</m:mi></m:msubsup>
<m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>E</m:mi><m:mi>i</m:mi></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>i</m:mi><m:mo>-</m:mo><m:mn>1</m:mn></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mi>i</m:mi></m:msub></m:msub><m:mo>⁢</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">e</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:mtext>,</m:mtext></m:math></td><td class="number">(4.4)</td></tr></table>
where <span class="eq"><m:math><m:mi>ε</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:msup><m:mrow><m:mo stretchy="false">(</m:mo><m:msub><m:mo>∂</m:mo><m:mrow><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup></m:mrow></m:msub><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo></m:mrow><m:mrow><m:mo>-</m:mo><m:mn>1</m:mn><m:mo>/</m:mo><m:mn>2</m:mn></m:mrow></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup><m:mo stretchy="false">)</m:mo></m:math></span> is a normalization function; <span class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mi mathvariant="sans-serif">b</m:mi></m:msub></m:math></span> is a quadrature eigenstate with <span class="eq"><m:math><m:mi mathvariant="bold">x</m:mi></m:math></span>-eigenvalue <span class="eq"><m:math><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup></m:math></span>; <span class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>e</m:mi><m:mi>i</m:mi></m:msub><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mi>i</m:mi></m:msub></m:msub></m:math></span>, <span class="eq"><m:math><m:msub><m:mi>e</m:mi><m:mi>i</m:mi></m:msub><m:mo>∈</m:mo><m:mo stretchy="false">{</m:mo><m:mn>0</m:mn><m:mo>,</m:mo><m:mn>1</m:mn><m:mo stretchy="false">}</m:mo></m:math></span>, denotes the canonical basis of <span class="eq"><m:math><m:mi mathvariant="script">H</m:mi></m:math></span>; <span class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>e</m:mi><m:mo>¯</m:mo></m:mover><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">e</m:mi><m:mo>¯</m:mo></m:mover></m:msub></m:math></span>, <span class="eq"><m:math><m:mn>0</m:mn><m:mo>≤</m:mo><m:mover accent="true"><m:mi>e</m:mi><m:mo>¯</m:mo></m:mover><m:mo>≤</m:mo><m:mn>1</m:mn></m:math></span>, form an orthogonal basis of <span class="eq"><m:math><m:msup><m:mi>L</m:mi><m:mn>2</m:mn></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mo stretchy="false">[</m:mo><m:mn>0</m:mn><m:mo>;</m:mo><m:mn>1</m:mn><m:mo stretchy="false">]</m:mo><m:mo stretchy="false">)</m:mo></m:math></span>. As the classical mapping <span class="eq"><m:math><m:mi mathvariant="script">E</m:mi></m:math></span> is invertible, <span class="eq"><m:math><m:mi mathvariant="script">QE</m:mi></m:math></span> is unitary with the appropriate normalization function <span class="eq"><m:math><m:mi>ε</m:mi></m:math></span>.
This mapping is defined to act on individual states, with the slice values <span class="eq"><m:math><m:msubsup><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span> as input in the system <span class="eq"><m:math><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span>, whose purpose is actually to hold Bob's sequentially corrected bit values <span class="eq"><m:math><m:msub><m:mi>β</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow></m:msub></m:math></span>. The complete transformation jointly involving <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> systems would be fairly heavy to describe. Only the ECC correction needs to be described jointly, and assuming it is correctly sized (i.e., <span class="eq"><m:math><m:msubsup><m:mi>l</m:mi><m:mi>i</m:mi><m:mtext>b</m:mtext></m:msubsup></m:math></span> are large enough), Bob has enough information to reconstruct Alice's bit values. Let us now sketch how the system <span class="eq"><m:math><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span> is constructed.
</p><p id="id2257605" class="figure"><img src="QE.png" type="image/png"/></p><p class="caption">Figure 4.2: Schematic description of <span class="eq"><m:math><m:mi mathvariant="script">QE</m:mi></m:math></span> and the use of the systems <span class="eq"><m:math><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span></p><p>
Assume that Bob first calculates, using a quantum circuit, the first slice estimator (classically: <span class="eq"><m:math><m:msub><m:mi>E</m:mi><m:mn>1</m:mn></m:msub><m:mo stretchy="false">(</m:mo><m:msup><m:mi>X</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo></m:math></span>), which does not depend on any syndrome. That is, he applies the following mapping, defined on the bases of <span class="eq"><m:math><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> and <span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span>: <span class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mi mathvariant="sans-serif">b</m:mi></m:msub><m:mo>→</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>E</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mn>1</m:mn></m:msub></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msub><m:mover accent="true"><m:mi mathvariant="sans-serif">e</m:mi><m:mo>¯</m:mo></m:mover><m:mn>1</m:mn></m:msub></m:msub></m:math></span> (up to normalization), where the function <span class="eq"><m:math><m:msub><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mn>1</m:mn></m:msub></m:math></span> is needed only to make the mapping unitary. From the <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> qubits in the <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> systems <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mn>1</m:mn></m:msub></m:math></span> and the syndrome sent by Alice <span class="eq"><m:math><m:mo>|</m:mo><m:msubsup><m:mi>ξ</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup><m:mo>〉</m:mo></m:math></span>, there exists a quantum circuit that calculates the relative syndrome of Alice's and Bob's bits, that is a superposition of the classical quantities
<span class="eq"><m:math><m:msubsup><m:mi>ξ</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup><m:mo>⊕</m:mo><m:msubsup><m:mi>H</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup><m:msub><m:mi>E</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msub><m:mi>X</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>l</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span>. From this, a quantum circuit calculates the coset leader of the syndrome, that is (a superposition of) the most probable difference vector between Alice's and Bob's qubits. An extra <span class="eq"><m:math><m:mi>l</m:mi><m:mo>-</m:mo><m:msubsup><m:mi>l</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup></m:math></span> blank qubits are needed for this operations; we assume they are all initialized to <span class="eq"><m:math><m:mo>|</m:mo><m:mn>0</m:mn><m:mo>〉</m:mo></m:math></span>:
<table class="equation" id="id2258045"><tr><td class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:msubsup><m:mi>H</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup><m:mo>⁢</m:mo><m:mo stretchy="false">(</m:mo><m:msubsup><m:mi>s</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup><m:mo>⊕</m:mo><m:msubsup><m:mi>e</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn><m:mrow><m:mo>′</m:mo><m:mo stretchy="false">(</m:mo><m:msubsup><m:mi>l</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:mn>0</m:mn><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn><m:mrow><m:mo>′</m:mo><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo>-</m:mo><m:msubsup><m:mi>l</m:mi><m:mn>1</m:mn><m:mtext>b</m:mtext></m:msubsup><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup></m:msub><m:mo>→</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:msubsup><m:mi>s</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup><m:mo>⊕</m:mo><m:msubsup><m:mi>e</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn><m:mrow><m:mo>′</m:mo><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup></m:msub><m:mtext>.</m:mtext></m:math></td><td class="number">(4.5)</td></tr></table>
Then, using a controlled-not operation between Bob's bits (control) and the difference vector (target), we produce <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> qubits containing the same bit values as Alice's, with an arbitrarily large probability:
<table class="equation" id="id2258260"><tr><td class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:msubsup><m:mi>e</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">e</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msubsup><m:mi>s</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup><m:mo>⊕</m:mo><m:msubsup><m:mi>e</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn><m:mrow><m:mo>′</m:mo><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup></m:msub><m:mo>→</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:msubsup><m:mi>e</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">e</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msubsup><m:mi>s</m:mi><m:mn>1</m:mn><m:mrow><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn><m:mrow><m:mo>′</m:mo><m:mo stretchy="false">(</m:mo><m:mi>l</m:mi><m:mo stretchy="false">)</m:mo></m:mrow></m:msubsup></m:msub><m:mtext>.</m:mtext></m:math></td><td class="number">(4.6)</td></tr></table>
This is how the <span class="eq"><m:math><m:mi>l</m:mi></m:math></span> systems <span class="eq"><m:math><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn><m:mo>′</m:mo></m:msubsup></m:math></span> are created.</p><p>Following this approach for the next slices, we can define: <span class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>s</m:mi><m:mn>1</m:mn></m:msub><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn><m:mo>′</m:mo></m:msubsup></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>E</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mn>1</m:mn></m:msub></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msub><m:mover accent="true"><m:mi mathvariant="sans-serif">e</m:mi><m:mo>¯</m:mo></m:mover><m:mn>1</m:mn></m:msub></m:msub><m:mo>→</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>s</m:mi><m:mn>1</m:mn></m:msub><m:mo>〉</m:mo></m:mrow><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn><m:mo>′</m:mo></m:msubsup></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>E</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mn>1</m:mn></m:msub></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mi>E</m:mi><m:mn>2</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msub><m:mi>s</m:mi><m:mn>1</m:mn></m:msub><m:mo stretchy="false">)</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">e</m:mi><m:mn>2</m:mn></m:msub></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msub><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mn>2</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msub><m:mi>s</m:mi><m:mn>1</m:mn></m:msub><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:msub><m:mover accent="true"><m:mi mathvariant="sans-serif">e</m:mi><m:mo>¯</m:mo></m:mover><m:mn>2</m:mn></m:msub></m:msub></m:math></span>, and reasonably assume that the bit value given in <span class="eq"><m:math><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mn>1</m:mn><m:mo>′</m:mo></m:msubsup></m:math></span> is equal to Alice's <span class="eq"><m:math><m:msub><m:mi>S</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>X</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>. This reasoning can be applied iteratively, so as to fill the system <span class="eq"><m:math><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span> with all the corrected bit values, and with an extra step to set <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:msub><m:mi>s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow></m:msub><m:mo stretchy="false">)</m:mo></m:math></span> in <span class="eq"><m:math><m:mover accent="true"><m:mi mathvariant="sans-serif">e</m:mi><m:mo>¯</m:mo></m:mover></m:math></span>.</p><p>As a last step, Bob can revert the ECC decoding operations and come back to the situation where he has blank qubits in <span class="eq"><m:math><m:msubsup><m:mi mathvariant="sans-serif">s</m:mi><m:mrow><m:mn>1</m:mn><m:mo>…</m:mo><m:mi>m</m:mi></m:mrow><m:mo>′</m:mo></m:msubsup></m:math></span> as depicted in Fig. <span class="ref"><a href="#id2257605">4.2</a></span>.</p><p>Finally, the qubits produced by <span class="eq"><m:math><m:mi mathvariant="script">QE</m:mi></m:math></span> can be transformed into EPR pairs using the CSS codes and the syndromes Alice sent to Bob.</p><h4 id="id2259027">4.2.2 Phase coherence</h4><p>Neither the unitary transformation <span class="eq"><m:math><m:mi mathvariant="script">QS</m:mi></m:math></span> nor <span class="eq"><m:math><m:mi mathvariant="script">QE</m:mi></m:math></span> take into account the modulation of the coherent state in the <span class="eq"><m:math><m:mi mathvariant="bold">p</m:mi></m:math></span>-quadrature. By ignoring what happens in the <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:math></span> system of Eq. (<span class="ref"><a href="#id2249189">2.3</a></span>), the reduced system <span class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub><m:mi mathvariant="sans-serif">b</m:mi></m:mrow></m:msub></m:math></span> lacks phase coherence:
<table class="equation" id="id2259105"><tr><td class="eq"><m:math><m:msub><m:mi>ρ</m:mi><m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub><m:mi mathvariant="sans-serif">b</m:mi></m:mrow></m:msub><m:mo>=</m:mo><m:mo>∫</m:mo><m:mo>ⅆ</m:mo><m:mi>x</m:mi><m:mo>ⅆ</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>ⅆ</m:mo><m:mi>p</m:mi><m:msqrt><m:mrow><m:msub><m:mi>G</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>⁢</m:mo><m:msub><m:mi>G</m:mi><m:mn>1</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo stretchy="false">)</m:mo></m:mrow></m:msqrt><m:mo>⁢</m:mo><m:msub><m:mi>G</m:mi><m:mn>2</m:mn></m:msub><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mo>⁢</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:msub><m:mrow><m:mo>〈</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>|</m:mo></m:mrow><m:mi>D</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mo>⁡</m:mo><m:msub><m:mrow><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mn>0</m:mn><m:mo>〉</m:mo></m:mrow><m:mi mathvariant="sans-serif">b</m:mi></m:msub><m:mrow><m:mo>〈</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>+</m:mo><m:mi>i</m:mi><m:mo>⁢</m:mo><m:mn>0</m:mn><m:mo>|</m:mo></m:mrow><m:mo>⁡</m:mo><m:msup><m:mi>D</m:mi><m:mo>†</m:mo></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>i</m:mi><m:mi>p</m:mi><m:mo stretchy="false">)</m:mo><m:mtext>.</m:mtext></m:math></td><td class="number">(4.7)</td></tr></table>
To remedy this, we assume that Alice also sends the <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:math></span> system to Bob, just like she does for the <span class="eq"><m:math><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> system and the syndromes, since the modulation in the <span class="eq"><m:math><m:mi mathvariant="bold">p</m:mi></m:math></span>-quadrature is independent of the key. Bob can take it into account before applying <span class="eq"><m:math><m:mi mathvariant="script">QE</m:mi></m:math></span>, by displacing his state along the <span class="eq"><m:math><m:mi mathvariant="bold">p</m:mi></m:math></span>-quadrature in order to bring it on the <span class="eq"><m:math><m:mi mathvariant="bold">x</m:mi></m:math></span>-axis.</p><p>Actually, we could formally include this <span class="eq"><m:math><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:math></span>-dependent operation in the <span class="eq"><m:math><m:mi mathvariant="script">QE</m:mi></m:math></span> mapping, by adding <span class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:mi>p</m:mi><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>2</m:mn></m:msub></m:msub></m:math></span> to its input and output (unmodified) and by multiplying by a factor of the form <span class="eq"><m:math><m:msup><m:mi>e</m:mi><m:mrow><m:mi>i</m:mi><m:mo>⁢</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>⁢</m:mo><m:mi>p</m:mi><m:mo>/</m:mo><m:mn>4</m:mn><m:mo>⁢</m:mo><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub></m:mrow></m:msup></m:math></span> in Eq. (<span class="ref"><a href="#id2256815">4.4</a></span>), with <span class="eq"><m:math><m:msub><m:mi>N</m:mi><m:mn>0</m:mn></m:msub></m:math></span> the vacuum fluctuations. For notation simplicity, however, we mention it here without explicitly writing it.</p><p>Also, for the simplicity of the notation in the next section, we can assume without loss of generality that the coefficients of <span class="eq"><m:math><m:mo>|</m:mo><m:mi>Ψ</m:mi><m:mo>〉</m:mo></m:math></span> in the <span class="eq"><m:math><m:mi mathvariant="bold">x</m:mi></m:math></span>-basis of <span class="eq"><m:math><m:mi mathvariant="sans-serif">b</m:mi></m:math></span> are real, after adjustment by Bob as a function of <span class="eq"><m:math><m:mi>p</m:mi></m:math></span>.</p><h4 id="id2259600">4.2.3 Construction of <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> and <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover></m:math></span></h4><p>Let us now make explicit the construction of the functions <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> and <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover></m:math></span>. First assume, for simplicity, that we have only one slice (<span class="eq"><m:math><m:mi>m</m:mi><m:mo>=</m:mo><m:mn>1</m:mn></m:math></span>)—for this we do not write the slice index as a subscript. The mapping has thus the following form:
<table class="equation" id="id2259674"><tr><td class="eq"><m:math><m:msub><m:mrow><m:mo>|</m:mo><m:mi>x</m:mi><m:mo>〉</m:mo></m:mrow><m:msub><m:mi mathvariant="sans-serif">a</m:mi><m:mn>1</m:mn></m:msub></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>〉</m:mo></m:mrow><m:mi mathvariant="sans-serif">b</m:mi></m:msub><m:mo>→</m:mo><m:mrow><m:mi>σ</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo></m:mrow><m:msub><m:mrow><m:mo>|</m:mo><m:mi>S</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:mi mathvariant="sans-serif">s</m:mi></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:mrow><m:mi>ε</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:mi>S</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo></m:mrow><m:msub><m:mrow><m:mo>|</m:mo><m:mi>E</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:mi mathvariant="sans-serif">e</m:mi></m:msub><m:msub><m:mrow><m:mo>|</m:mo><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>,</m:mo><m:mi>S</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo stretchy="false">)</m:mo><m:mo>〉</m:mo></m:mrow><m:mover accent="true"><m:mi mathvariant="sans-serif">e</m:mi><m:mo>¯</m:mo></m:mover></m:msub><m:mtext>,</m:mtext></m:math></td><td class="number">(4.8)</td></tr></table>
where <span class="eq"><m:math><m:mi>σ</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:msup><m:mrow><m:mo stretchy="false">(</m:mo><m:msub><m:mo>ⅆ</m:mo><m:mi>x</m:mi></m:msub><m:mo>⁡</m:mo><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo></m:mrow><m:mrow><m:mo>-</m:mo><m:mn>1</m:mn><m:mo>/</m:mo><m:mn>2</m:mn></m:mrow></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:mi>x</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>,
<span class="eq"><m:math><m:mi>ε</m:mi><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:mi>s</m:mi><m:mo stretchy="false">)</m:mo><m:mo>=</m:mo><m:msup><m:mrow><m:mo stretchy="false">(</m:mo><m:msub><m:mo>∂</m:mo><m:mrow><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup></m:mrow></m:msub><m:mo>⁡</m:mo><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover><m:mo stretchy="false">)</m:mo></m:mrow><m:mrow><m:mo>-</m:mo><m:mn>1</m:mn><m:mo>/</m:mo><m:mn>2</m:mn></m:mrow></m:msup><m:mo>⁡</m:mo><m:mo stretchy="false">(</m:mo><m:msup><m:mi>x</m:mi><m:mo>′</m:mo></m:msup><m:mo>,</m:mo><m:mover accent="true"><m:mi>s</m:mi><m:mo>¯</m:mo></m:mover><m:mo>,</m:mo><m:mi>s</m:mi><m:mo stretchy="false">)</m:mo></m:math></span>, and <span class="eq"><m:math><m:mover accent="true"><m:mi>S</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> and <span class="eq"><m:math><m:mover accent="true"><m:mi>E</m:mi><m:mo>¯</m:mo></m:mover></m:math></span> range between <span class="eq"><m:math><m:mn>0</m:mn></m:math></span> and <span class="eq"><m:math><m:mn>1</m:mn></m:math></span>.
</p><p>Let us take some state <span class="eq"><m:math><m:mi>ρ</m:mi></m:math></span> of the systems <span class="eq"><m:math><m:mi mathvariant="sans-serif">s</m:mi><m:mover accent="true"><m:mi mathvariant="sans-serif">s</m:mi><m:mo>¯</m:mo></m:mover><m:mi mathvariant="sans-serif">e</m:mi>