Arxiv:1903.00012V1 [Quant-Ph] 28 Feb 2019 Brational Modes of Trapped Ions [11]
Total Page:16
File Type:pdf, Size:1020Kb
All-Gaussian universality and fault tolerance with the Gottesman-Kitaev-Preskill code Ben Q. Baragiola,1 Giacomo Pantaleoni,1 Rafael N. Alexander,2 Angela Karanjai,3 and Nicolas C. Menicucci1 1Centre for Quantum Computation and Communication Technology, School of Science, RMIT University, Melbourne, Victoria, Australia 2Center for Quantum Information and Control, Department of Physics and Astronomy, University of New Mexico, Albuquerque, USA 3Centre for Engineered Quantum Systems, School of Physics, The University of Sydney, Sydney, Australia (Dated: Monday 4th March, 2019) The Gottesman-Kitaev-Preskill (GKP) encoding of a qubit within an oscillator is particularly ap- pealing for fault-tolerant quantum computing with bosons because Gaussian operations on encoded Pauli eigenstates enable Clifford quantum computing with error correction. We show that applying GKP error correction to Gaussian input states, such as vacuum, produces distillable magic states, achieving universality without additional non-Gaussian elements. Fault tolerance is possible with sufficient squeezing and low enough external noise. Thus, Gaussian operations are sufficient for fault-tolerant, universal quantum computing given a supply of GKP-encoded Pauli eigenstates. Introduction.|The promise of a quantum computer 21] due to the built-in redundancy afforded by their lies in its ability to dramatically outpace classical com- infinite-dimensional Hilbert space. High precision con- puters for certain tasks [1]. The celebrated Threshold trollability of optical-cavity [22, 23] and vibrational [11] Theorem [2] proves that this feature survives even in the modes further enhances their appeal. presence of (low enough) noise|a property called fault Using a bosonic code, one may define logical-Clifford tolerance, which is based on quantum error correction [3]. QC, comprising encoded Pauli eigenstates and logical- Computation using operations restricted to Pauli- Clifford operations, which allows error correction at the eigenstate preparation, Clifford transformations, and encoded-qubit level. This, too, is an efficiently simu- Pauli measurements|henceforth referred to as Clifford lable subtheory and thus requires a logical-non-Clifford quantum computing (QC)|provides all the necessary element for universality. tools for quantum error correction but cannot outperform Our main result is that a magic state for logical- classical computation since it is efficiently simulable on a Clifford QC (using a particular bosonic code) can be classical computer [4]. Universal quantum computation found within Gaussian QC. Thus, the union of these two requires supplementing Clifford QC by a non-Clifford re- simulable subtheories is universal and|with low enough source|that is, a preparation, gate or measurement that physical noise|fault tolerant. is not an element of Clifford QC. Often this resource is a non-Pauli eigenstate, referred to as a magic state. The The bosonic code that enables this is the Gottesman- union of Clifford QC and a supply of magic states is uni- Kitaev-Preskill (GKP) encoding of a qubit into an os- versal for quantum computing and can be made fault cillator [15]. This is the only bosonic code that allows tolerant if the physical noise level is low enough [5]. logical-Clifford QC and CV-level error correction to be implemented using Gaussian QC along with a supply of The continuous-variable (CV) analog of Clifford QC logical-Pauli eigenstates, which are non-Gaussian.1 Un- is Gaussian QC, which includes Gaussian state prepara- til now, GKP logical-Clifford QC has been elevated to tion, Gaussian operations (i.e., Hamiltonians quadratic universal quantum computation through non-Gaussian ina; ^ a^y), and homodyne detection. CV systems arise gates (cubic phase gate) or the preparation of a logi- naturally in many quantum architectures, including op- cal magic state (logical Hadamard eigenstate or cubic tical modes [6{9], microwave-cavity modes [10], and vi- phase state) [15], all three of which require additional arXiv:1903.00012v1 [quant-ph] 28 Feb 2019 brational modes of trapped ions [11]. Gaussian QC lends non-Gaussian machinery beyond the logical-Pauli eigen- itself to optics because the nonlinearities required are lim- state. ited and low order and because homodyne detection is very high efficiency. Gaussian QC is also efficiently sim- In what follows, we show how to produce a distill- ulable by a classical computer [12] and therefore requires able [5] GKP magic state using GKP error correction on a a non-Gaussian resource (preparation, gate, or measure- thermal state (vacuum or finite temperature), along with ment) to elevate it to universal QC [13]. a complete analysis of the success probability of prepar- Fault tolerance requires discrete quantum information. ing a high-quality magic state from any given thermal Bosonic quantum-error-correcting codes (bosonic codes for short) embed discrete quantum information into CV systems in a way that maps general CV noise into ef- 1 Note: These states cannot be prepared [24] using the Gaussian fective noise acting on the encoded qubits [14{17]. Such measurements that correspond to logical-Pauli measurements be- codes are promising for fault-tolerant computation [18{ cause these measurements are destructive. 2 state. Our result applies to both square- and hexagonal- error correction strives to avoid these unintended logical lattice GKP encodings [16]. operations (residual errors). In what follows, we apply Notation and conventions|Here we define nota- the machinery of error correction to a known Gaussian tion and conventions to be used throughout this Let- state, which means the outcome-dependent final state is ter. We define positionq ^ := p1 (^a +a ^y) and momen- known perfectly. 2 GKP error correction [15, 18] proceeds in two steps: tump ^ := p−i (^a − a^y) for any modea ^. This means 2 First, the one quadrature is corrected, then the conjugate [^q; p^] = i, with a vacuum variance of 1 in each quadrature 2 quadrature. We define the Kraus operator that corrects and ~ = 1. ^ q just the q quadrature KEC(t) via the circuit (read right The Weyl-Heisenberg displacement operators to left): X^(s) := e−isp^ and Z^(s) := eisq^ displace a state by +s in position and momentum, respectively. For brevity, in in Kq (t) = X(−t) • we also define a joint displacement V^ (s) := Z^(sp)X^(sq), EC where s = (s ; s )T. q p t htj • j0Li ~ t p The functions (s) := qhsj i and (s) := phsj i denote position- and momentum-space wave functions for a ^ iq^⊗q^ where the controlled operation is CZ = e , and t 2 R state j i, respectively (tilde indicates momentum space). is the measurement outcome. This circuit differs from Any function, including wave functions, can be evalu- R the original [15] in that the correction here is a negative ated with respect to position, '(^q) := ds '(s)jsiqqhsj, to displacement by t ratherp than by t rounded to the nearest produce an operator diagonal in the position basis|and integer multiple of π. The outputs may differ by a log- : p similarly for momentum. Finally, we define X∆(x) = ical operation X^(± π), but this is unimportant because P δ(x − n∆) as a Dirac comb with spacing ∆. n2Z the input state is known. The GKP encoding|In the original square-lattice ^ q ~ ^ Direct evaluation shows KEC(t) = 0;L(^q)X(−t). A GKP encoding [15], the wave functions for the logical ba- similar calculation shows that the Kraus operator for cor- sis states fj0 i ; j1 ig are Dirac combs in position space ^ p ~ ^ L L p recting the p quadrature is K (t) = 0;L(^p)Z(−t). Ap- p EC with state-dependent offset: j;L(s) = X2 π(s − j π) plying both corrections (in either order since they com- for j 2 f0; 1g. Their momentum-space wave functions mute up to a phase) performs full GKP error correction: are also Dirac combs but with no offset, different spac- ing, and a relative phase between the spikes: ~ (s) = ^ ^ p ^ q ^ ^ p j;L KEC(t) = KEC(tp)KEC(tq) = ΠGKPV (−t); (3) p1 (−1)js= π p (s). Note that the momentum-space 2 X π T spikes for j1Li alternate sign, and those for j0Li are uni- with measurement outcomes t = (tq; tp) . This Kraus form. operator (i) displaces the state by an outcome-dependent ^ GKP logical operators X^L and Z^L are implemented by amount, V (−t), and then (ii) projects it back into the p p ^ displacements X^( π) and Z^( π), respectively, while dis- GKP logical subspace with ΠGKP. p ^ placements by integer multiples of 2 π in either quadra- Applying KEC(t) to an input stateρ ^in produces ^ ^ ^ y ture leave the GKP logical subspace invariant. For later the unnormalized state ρ¯(t) = KEC(t)^ρinKEC(t), use, we define the four GKP-encoded logical Paulis where the bar indicates lack of normalization. The joint probability density function (pdf) for the out- µ X µ ^ σ^L := σjk jjLihkLj ; (1) comes, pdf(t) = Tr[ρ¯(t)], normalizes the output state: jk ρ^(t) = ρ¯^(t)= pdf(t). Bloch vector for the error-corrected state.| where σµ is the jk'th element of Pauli matrix σµ (with jk Using the logical basis in Eqs. (1) we represent the out- 0 µ µ σ = I). Note thatσ ^L have support only on the GKP put stateρ ^(t) = 1 P r (t)^σ by a 4-component Bloch ^ ^ 2 µ µ L logical subspace, while XL and ZL have full support. Fi- vector r(t) with outcome-dependent coefficients rµ(t) := nally, we denote the (rank-two) projector onto the GKP µ Tr[^ρ(t)^σL]. For the unnormalized state,r ¯0(t) = pdf(t), logical subspace [15, 25] and for the normalized state, r0(t) = 1. In what follows, we use the notation r = (r ; ~r), where ~r is the ordinary ^ 0 ~ ~ ~ ~ 0 ΠGKP :=σ ^L = 0;L(^q) 0;L(^p) = 0;L(^p) 0;L(^q): (2) (3-component) Bloch vector within r.