## A Search for High-Threshold Qutrit Magic State Distillation Routines

Shiroman Prakash and Rishabh Singhal

Department of Physics and Computer Science, Dayalbagh Educational Institute, Agra, India

## Abstract

Determining the best attainable threshold for qudit magic state distillation is directly related to the question of whether or not contextuality is sufficient for universal quantum computation. We carry out a search for high-threshold magic state distillation routines for a highly-symmetric qutrit magic state known as the strange state. Our search covers a large class of [[ n, 1]] 3 qutrit stabilizer codes with up to 23 qutrits, and is facilitated by a theorem that relates the distillation performance of a qudit stabilizer code to its weight-enumerators. We could not find any code with n &lt; 23 qutrits that distills the strange state with better than linear noise suppression, other than the 11-qutrit Golay code. However, for n = 23, we find over 600 CSS codes that can distill the qutrit strange state with cubic noise suppression. While none of these codes surpass the threshold of the 11-qutrit Golay code, their existence suggests that, for large codes, the ability to distill the qutrit strange state is somewhat generic.

## Contents

| 1 Introduction   | 1 Introduction                                          | 1 Introduction                                                        |   2 |
|------------------|---------------------------------------------------------|-----------------------------------------------------------------------|-----|
| 2                | Preliminaries                                           | Preliminaries                                                         |   4 |
|                  | 2.1                                                     | The Heisenberg-Weyl displacement group . . . . . . . . . . . .        |   4 |
|                  | 2.2                                                     | Qudit stabilizer codes . . . . . . . . . . . . . . . . . . . . . . .  |   4 |
|                  | 2.3                                                     | Weight enumerators . . . . . . . . . . . . . . . . . . . . . . . .    |   5 |
|                  | 2.4                                                     | Discrete Wigner functions . . . . . . . . . . . . . . . . . . . . .   |   6 |
| 3                | Distillation and weight enumerators                     | Distillation and weight enumerators                                   |   7 |
|                  | 3.1                                                     | A general formulation in terms of complete weight enumerators         |   7 |
|                  | 3.2                                                     | The strange state and simple-weight enumerators . . . . . . . .       |   9 |
|                  |                                                         | 3.2.1 Example: the [[11 , 1 , 5]] 3 Golay code . . . . . . . . . . .  |  10 |
|                  | 3.3                                                     | Conditions for magic state distillation . . . . . . . . . . . . . .   |  11 |
| 4                | Search for distillation routines                        | Search for distillation routines                                      |  12 |
|                  | 4.1                                                     | Equivalence of distillation routines . . . . . . . . . . . . . . . .  |  12 |
|                  | 4.2                                                     | Small quantum error-correcting codes . . . . . . . . . . . . . .      |  14 |
|                  | 4.3                                                     | CSS codes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |  14 |
| 5                | Discussion                                              | Discussion                                                            |  16 |
| A                | Projecting phase point operators onto a stabilizer code | Projecting phase point operators onto a stabilizer code               |  19 |
| B                | Some codes that do not distill the strange state        | Some codes that do not distill the strange state                      |  20 |

## 1 Introduction

Contextuality was identified as a necessary and possibly-sufficient condition for universal quantum computing in [1]. The argument of [1] is based on magic state distillation [2] for qudits of odd-prime dimension, and later extended to qudits of arbitrary odd dimensions in [3], and continuous variable systems in [4]. 1 The authors of these works showed that qudit states that do not exhibit contextuality with respect to stabilizer measurements have a non-negative discrete Wigner function [5,6]. The set of such states is known as the Wigner polytope. Because Clifford unitaries and stabilizer measurements are efficiently simulable for states in the Wigner polytope [7,8], they thus cannot be distilled into pure magic states.

In order to demonstrate that contextuality is not only necessary but also sufficient for qutrit quantum computation, one must demonstrate that a supply of qudits that do exhibit contextuality with respect to stabilizer measurements may be used to achieve universal quantum computation. In the language of magic state distillation, this translates into the question, can any qudit mixed state outside the Wigner polytope be distilled into a pure magic state?

For qudits of odd prime dimension p , the Wigner polytope is a convex polytope with p 2 facets that lives in the p 2 -1 dimensional space of qudit density matrices. In [9], it was shown that no finite magic state distillation routine can distill all states that lie outside one of the faces of the Wigner polytope , generalizing the analogous result for qubit states that 2 lie outside the stabilizer polytope [11]. However, the possibility remains that a sequence of magic state distillation routines, based on stabilizer codes of increasing length n may distill states arbitrarily close to a face of the Wigner polytope. Is there any evidence that such a sequence of magic state distillation routines exist?

The problem simplifies if one focuses on qutrits. There exists a qutrit magic state, first identified by Howard and van Dam [12], sometimes known in the literature in the qutrit strange state | S ⟩ [8], that lies directly above one of the facets of the Wigner polytope . 3 As discussed in [14], noisy | S ⟩ states can be twirled via Clifford unitaries to lie on a line connecting a pure | S ⟩ state to the maximally mixed state:

$$\hat { \rho } ( \epsilon ) = ( 1 - \epsilon ) \left | S \right \rangle \left \langle S \right | + \epsilon \frac { 1 } { 3 },$$

with all noise converted to depolarizing noise, parameterized by ϵ . Any state ˆ( ) ρ ϵ for ϵ &lt; 3 4 lies / outside the Wigner polytope and exhibits contextuality with respect to stabilizer measurements. We then ask whether or not there exists a family of n -to-1 magic state distillation routines that distill the strange state with a threshold approaching ϵ = 3 4 as /

1 There are certain subtleties associated with state-independent contextuality for qubits and qudits of even dimension - for simplicity, we focus exclusively on qudits of odd-prime dimension in this paper.

2 Some distillation routines that distill qutrit states up to one of the hyperedges of the Wigner polytope were found in [10].

3 See [13] for another application of the strange state.

## n →∞ ? 4

At the time [1] was published, and for several years thereafter, while some qutrit and qudit magic state distillation routines had been proposed [15-18, 10], no magic state distillation routine that distilled the Howard van Dam strange state was known. It was later discovered that an 11-qutrit CSS code based on the ternary Golay code can distill the | S ⟩ state, with a threshold of ϵ ∗ = 0 38 [19]. . Do there exist any other qutrit stabilizer codes that distill the strange state? If so, how do their thresholds compare to that of the 11-qutrit Golay code? 5 In this paper, we carry out a computational search over reasonably small qutrit error-correcting codes to help answer these questions.

To facilitate our search, we prove a simple theorem connecting the performance of a stabilizer code for qudit magic state distillation its complete weight enumerator. For distillation routines for the qutrit strange state, this formula simplifies and depends only on the stabilizer code's simple weight enumerator.

Our search is limited by practical considerations, as well as the existing classifications of qutrit error-correcting codes in the literature - namely, the classification of qutrit stabilizer states in [21,22] and a classification of self-orthogonal classical ternary codes available on [23]. We carried out a search over all [[ n, 1]] 3 stabilizer codes with n ≤ 9, and a search over all [[11 , 1]] 3 stabilizer codes that can be obtained from a [[12 0 6]] , , 3 stabilizer state. For such codes, we demand transversality of a particular single-qutrit gate (the square of the qutrit Hadamard gate), which allows us to restrict our search to projection onto the trivial syndrome of each stabilizer code. We also searched over all [[ n, 1]] 3 CSS codes for odd n ≤ 23 that possess a complete set of transversal single-qudit Clifford gates.

We found that none of the codes we searched with n &lt; 23 could distill the | S ⟩ state with better-than-linear 6 noise suppression, other that the 11-qutrit Golay code of [19]. However, for n = 23, we found over 600 CSS codes that could distill the strange state with cubic noise suppression - which is approximately 1 / 3 of all the codes we could construct from the ternary self-orthogonal codes listed in [23] - suggesting that for large codes, magic state distillation is somewhat generic. None of these 23-qutrit codes, however, had a threshold that exceeds that of the 11-qutrit Golay code.

Let us point out that the new distillation routines we find here are mainly of theoretical interest. The success probabilities are quite low, and far better yields are obtained via triorthogonal codes [24-27] (see [16,28-30] for constructions of qutrit and qudit triorthogonal codes). Nevertheless, the CSS codes we study have a complete set of transversal Clifford gates, and may turn out to be useful for fault-tolerant quantum computation in other settings.

Our paper is organized as follows. In section 2 we briefly review some background material.

4 An additional consideration is the increase in overhead cost of distillation as n →∞ .

5 Recently, [20] reported the existence of 13-qutrit and 29-qutrit codes that could distill the qutrit strange state with high thresholds. However, we attempted to reproduce these results and found that neither of the codes reported in [20] are able to distill the qutrit strange state. See Appendix B.

6 We also found a few 9 and 11 qutrit codes that could distill the state with linear noise suppression.

In section 3 we derive a relation between weight-enumerators and the performance of magic state distillation routines. In section 4 we describe our search space and the results. In section 5 we conclude with some brief discussion.In addition, in appendix A we present a useful lemma, and in appendix B we study a two qutrit codes, conjectured in [20] to distill the | S ⟩ state, that we find do not distill the | S ⟩ state.

## 2 Preliminaries

Here we review some well-done definitions and results concerning stabilizer codes for qudits of dimension p where p is any odd prime.

## 2.1 The Heisenberg-Weyl displacement group

The qudit version of the Pauli group is generated by

$$\hat { X } = \sum _ { k } \left | k + 1 \right \rangle \left \langle k \right |, \ \ \hat { Z } = \sum _ { k } \omega ^ { k } \left | k \right \rangle \left \langle k \right |,$$

and multiplication by ω = e 2 πi/p , and is known as the Heisenberg-Weyl displacement group. The conventions for Heisenberg Weyl displacement operators we will follow are those in [1],

$$\hat { D } ( u, v ) = \omega ^ { - 2 ^ { - 1 } u v } \hat { X } ^ { u } \hat { Z } ^ { v }.$$

n -qudit Heisenberg Weyl-displacement operators are denoted as

$$\hat { D } ( \vec { u }, \vec { v } ) = \hat { D } ( u _ { 1 }, v _ { 1 } ) \otimes \hat { D } ( u _ { 2 }, v _ { 2 } ) \otimes \dots \otimes \hat { D } ( u _ { n }, v _ { n } ),$$

where ⃗ u = ( u , u 1 2 , . . . , u n ) and ⃗ v = ( v , v 1 2 , . . . , v n ).

It is convenient to combine ⃗ u and ⃗ v into a symplectic vector χ = ( ⃗ u, ⃗ v ), and write ˆ ( D χ ) = ˆ ( D ⃗ u, ⃗ v ). Then

$$\hat { D } ( \chi ) \hat { D } ( \chi ^ { \prime } ) = \omega ^ { 2 ^ { - 1 } [ \chi, \chi ^ { \prime } ] } D _ { \chi + \chi ^ { \prime } },$$

where

$$[ \chi, \chi ^ { \prime } ] = \vec { u } \cdot \vec { v } ^ { \prime } - \vec { u } ^ { \prime } \cdot \vec { v }$$

is the symplectic inner product (and · is the usual dot product). For commuting HeisenbergWeyl operators [ χ, χ ′ ] vanish.

## 2.2 Qudit stabilizer codes

An [[ n, k ]] p stabilizer code is defined by the simultaneous eigenspace of n -k independent commuting n -qudit Heisenberg-Weyl displacement operators, and has dimension p k . We will assume that each of these operators are of the form given in (2.1), with no additional overall phase. As such, to completely specify the stabilizer code, we also have to specify

the eigenvalues of each stabilizer, which are of the form ω a , where a can be taken to be an element of Z p . In this paper we will restrict our attention to codes defined by the 1 = ω 0 eigenspace of each operator, for reasons we explain below. The set of commuting operators { ˆ ( D ⃗ u 1 , ⃗ v 1 ) , . . . , D ˆ ( ⃗ u n -k , ⃗ v n -k ) } can be specified by a symplectic matrix

$$H = \left ( \begin{array} { c | c } \vec { u } _ { 1 } & \vec { v } _ { 1 } \\ \vec { u } _ { 2 } & \vec { v } _ { 2 } \\ \vdots & \vdots \\ \vec { u } _ { n - k } & \vec { v } _ { n - k } \end{array} \right ),$$

such that the symplectic inner product of any two rows vanishes. These operators generate an abelian group of order p n -k , which we refer to as the stabilizer group S . Because S is abelian, there are no overall phases arising from multiplication and S consists of all Heisenberg-Weyl Displacement operators of the form ˆ ( D χ ) where χ is in the row-span of H .

A stabilizer code can equivalently be thought of as an additive classical code over GF p ( 2 ), that is self-orthogonal under the Hermitian inner product [31,32].

We use the notation N ( S ) to denote all n -qudit Heisenberg-Weyl displacement operators that commute with S . S ⊆ N ( S ). In general, N ( S ) will not form an abelian group, and multiplication of two different elements may give rise to Heisenberg-Weyl operators with overall phases. It is nevertheless useful to define N ( S ) ∗ as N ( S ) modulo overall phases (i.e., identify Heisenberg-Weyl operators that differ by an overall phase), and specify N ( S ) ∗ as the row-span of a symplectic matrix. In this way N ( S ) ∗ can be thought of as a classical additive code in GF p ( 2 ), denoted as S ⊥ , that is the set of all codewords orthogonal to S under the Hermitian inner product. We use the notation S ⊥ and N ( S ) ∗ interchangeably.

Clearly, S ⊂ S ⊥ . Any L ∈ S ⊥ defines a coset of S in S ⊥ which is of the form 7 S + L = { L + M M | ∈ S} . All stabilizers in the same coset will commute with each other.

## 2.3 Weight enumerators

Let us briefly review the definition of the weight enumerator of a stabilizer code.

Let S denote the stabilizer group of an [[ n, k, d ]] 3 stabilizer code. The projector onto the codespace of S with trivial syndrome can be written explicitly as,

$$\hat { \Pi } _ { \mathcal { S } } ^ { 0 } = \frac { 1 } { p ^ { n - k } } \sum _ { M \in \mathcal { S } } M.$$

Let ˜ be a either a stabilizer code S S , or one of its cosets in S ⊥ . The complete weight enumerator of ˜ is a function of S p 2 formal variables { y ij } , for i, j = 0 , . . . , p -1, defined as

7 In this expression, addition means addition as symplectic vectors (or equivelently, vectors in GF p ( 2 ) n ), or multiplication as Heisenberg-Weyl displacement operators.

follows:

$$w ( \tilde { \mathcal { S } } ; \omega _ { i j } ) = \sum _ { ( \vec { x } | \vec { z } ) \in \tilde { \mathcal { S } } } \prod _ { i = 1 } ^ { n } y _ { x _ { i }, x _ { j } }$$

The simple weight enumerator of ˜ to be a function of two formal variables S x and y , defined as follows:

$$w ( \dot { \mathcal { S } } ; x, y ) = w ( \dot { \mathcal { S } } ; y _ { i j } ( x, y ) ).$$

where

̸

$$y _ { i j } ( x, y ) = \begin{cases} x & ( i, j ) = ( 0, 0 ) \\ y & ( i, j ) \neq ( 0, 0 ). \end{cases}$$

If we further set x = y = 1, then the simple weight enumerator becomes equal to |S| ˜ . These definitions coincide with the standard definitions of the simple and complete weight enumerators of a classical error-correcting code over the finite field GF p ( 2 ).

A MacWilliams identity relates w ( S ; 1 , z ) to w ( S ⊥ ; 1 , z ) [33, 32]. This is:

$$w ( S ^ { \perp } ; 1, z ) = \frac { ( 1 + ( p ^ { 2 } - 1 ) z ) ^ { n } } { p ^ { n - k } } w \left ( S ; 1, \frac { 1 - z } { 1 + ( p ^ { 2 } - 1 ) z } \right ).$$

It is conventional to define A S ( z ) = w ( S ; 1 , z ) and B S ( z ) = w ( S ⊥ ; 1 , z ). A S ( z ) -B S ( z ) is a polynomial in z with non-negative coefficients and is the simple weight enumerator of the set of all logical operators for the stabilizer code S . The lowest power of z that appears in A z ( ) -B z ( ) is z d , where d is the distance of the code.

## 2.4 Discrete Wigner functions

A discrete phase space formalism for qudits was formulated in [6,34,5], that we briefly review here.

Singe-qudit phase point operators are defined as,

$$\hat { A } ( 0, 0 ) = \frac { 1 } { p } \sum _ { x, z } \hat { D } ( u, v ), \quad \hat { A } ( u, v ) = \hat { D } ( u, v ) \hat { A } ( 0, 0 ) \hat { D } ( u, v ) ^ { \dagger }.$$

These are normalized so that tr ˆ ( A u, v ) = 1. Multi-qudit phase point operators are defined as

$$\hat { A } ( \vec { u }, \vec { v } ) = \bigotimes _ { i } \hat { A } ( u _ { i }, v _ { i } ),$$

and we will use the notation ˆ (0 A , 0) ⊗ n = ˆ (0 0) A , ⊗ ˆ (0 A , 0) ⊗ . . . ⊗ ˆ (0 A , 0).

Any qudit density matrix ˆ can be written as a linear combination of the phase point ρ operators with real, but possibly negative, coefficients. These coefficients define the discrete Wigner function of the qudit state:

$$\hat { \rho } = \sum _ { u, v } W ( \hat { \rho } ; u, v ) \hat { A } ( u, v ), \ \ W ( \hat { \rho } ; u, v ) = \frac { 1 } { p } \text{tr} \ \hat { \rho } \hat { A } ( u, v )$$

Note that tr ˆ = ρ ∑ u,v W ρ u,v (ˆ; ) = 1, and the discrete Wigner function defines a quasiprobability distribution.

The Clifford group acts covariantly on the discrete Wigner function, for qudits of odd prime dimension. Heisenberg-Weyl displacement operators act as translations in discrete phase, and a general Clifford unitary can be represented as a symplectic rotation followed by a translation. See [5,35,7,8], and references therein, for more details. This allows one to use the discrete Wigner function to define an efficient classical simulation for Clifford unitaries and stabilizer measurements acting on qudit states with non-negative Wigner functions [7], and develop resource theories of magic for qudits of odd prime dimension - see [8,36,37].

## 3 Distillation and weight enumerators

In the magic state model of fault-tolerant quantum computing [2, 38], one begins with a quantum computer that can initialize qudits in the computational basis, perform Clifford unitaries, and carry out stabilizer measurements. These operations are conventionally assumed to be noise-free for simplicity. To obtain universal quantum computing we supplement this quantum computer with the ability to initialize ancilla qudits in certain non-stabilizer states known as magic states, which are noisy. A magic state distillation protocol is a way to distill an arbitrarily pure magic state from many noisy magic states using only Clifford unitaries and stabilizer measurements. For the purposes of determining the best attainable threshold [39], any magic state distillation protocol for the qutrit strange state can be thought of as a procedure that projects n noisy qudits onto the codespace of an [[ n, k ]] p stabilizer code - if the projection is successful, one decodes the resulting qudits to obtain k hopefully-less-noisy magic states. 8

## 3.1 A general formulation in terms of complete weight enumerators

Here we provide a general formulation for qudit magic state distillation in the language complete weight enumerators. These results follow in part from the formulation in [9]. Let us also mention that a similar formulation in terms of signed-weight enumerators for qubit magic state distillation was given in [40].

Let S be a group of commuting Heisenberg-Weyl operators acting on n qudits. (Recall that, each element of S is assumed to be of the form given in equation (2.1), with no additional overall phase.) The projector onto the joint eigenspace of the code, with eigenvalues 1 = ω 0 , is given by,

$$\frac { } { } \hat { \Pi } _ { \mathcal { S } } ^ { 0 } = \frac { 1 } { | \mathcal { S } | } \sum _ { M \in \mathcal { S } } M$$

8 Sometimes, i.e., in the case of triorthogonal codes, error-correction is also possible prior to decoding, which can result in a substantial reduction in overhead. However, this is not possible for distillation routines for the qutrit strange state.

Theorem 1. Let ˆ ρ be a single qudit mixed state described by the Wigner function W ρ x,z (ˆ; ) . Then the probability for successful projection onto the eigenspace of S with trivial syndrome is given by the complete weight enumerator of S ⊥ , w ( S ⊥ ; { W ρ i,j (ˆ; ) } ) , with the formal variables y ij in the complete Weight enumerator replaced by the entries of the Wigner function of ρ :

$$\mu. \quad \text{tr} \, \left ( \hat { \Pi } _ { \mathcal { S } } ^ { 0 } \rho ^ { n } \right ) = w ( \mathcal { S } ^ { \$$

Proof. Let the single-qudit state ˆ be described by Wigner function ρ W in (ˆ; ρ u, v ). The n -qudit Wigner function for ρ ⊗ n is given by,

$$W _ { i n } ^ { ( n ) } ( \hat { \rho } ^ { \otimes n } ; \vec { u }, \vec { v } ) = \prod _ { i } W _ { i n } ( \$$

The probability for successful projection onto the codespace is,

$$\nu = \text{tr} \left ( \hat { \Pi } _ { \mathcal { S } } ^ { 0 } \rho _ { i n } ^ { \otimes n } \right ) \quad \quad \quad$$

$$= \sum _ { \vec { u }, \vec { v } } ^ { \theta ^ { \prime } \, u ^ { \prime } } \left ( \prod _ { i = 1 } ^ { n } W _ { i$$

In Appendix A, Lemma 2, we show that

$$\text{tr} \, \left ( \hat { \Pi } _ { \mathcal { S } } ^ { 0 } \hat { A } ( \vec { u }, \vec { v } ) \right ) = \begin{cases} 0 & ( \vec { u } | \vec { v } ) \notin \mathcal { S } ^ { \perp } \\ 1 & ( \vec { u } | \vec { v } ) \in \mathcal { S } ^ { \perp } \end{cases}.$$

We therefore find that,

$$\mu = \sum _ { ( \vec { u } | \vec { v } ) \in S ^ { \perp } } \prod _ { i } W _ { i n } ( \hat { \rho } ; u _ { i }, v _ { i } ).$$

In magic state distillation, we first project onto the codespace of S , which we assume for simplicity to be an [[ n, 1]] p code, then decode, to obtain a new state ρ ′ = f MSD ( ρ ). We can use the above lemma to compute the Wigner function of ρ ′ in terms of the complete weight enumerators of S ⊥ and its cosets.

Choose any representatives of the logical Pauli operators ¯ X and ¯ , Z to define logical Heisenberg-Weyl displacement operators ¯ ( D x,z ) and logical phase point operators ¯ ( A x, z ). The Wigner function of ˆ ρ ′ is given by,

$$W ( \hat { \rho } ^ { \prime } ; x, z ) = \frac { 1 } { \nu } \text{tr} \, \left ( p ^ { - 1 } \bar { A } ( x, z ) \hat { \Pi } _ { \mathcal { S } } ^ { 0 } \hat { \rho } ^ { \otimes n } \right ).$$

Let,

$$\iota \varkappa, \\ p ^ { - 1 } \bar { A } ( 0, 0 ) \hat { \Pi } _ { \mathcal { S } } ^ { 0 } = p ^ { - 2 } \sum _ { u, v } \bar { D } ( u, v ) \hat { \Pi } _ { \mathcal { S } } ^ { 0 } = \frac { 1 } { p ^ { n + 1 } } \sum _ { ( \vec { u } | \vec { v } ) \in \mathcal { S } ^ { \perp } } \hat { D } ( \vec { u }, \vec { v } ) \equiv \hat { \Pi } _ { \mathcal { S } ^ { \perp } } ^ { 0 }.$$

Although Π ˆ 0 S ⊥ is not a projector onto the codespace of a stabilizer code, in Appendix A, we also show that,

$$W ( \hat { \rho } ^ { \prime } ; 0, 0 ) = \frac { 1 } { \nu } \text{tr} \, \left ( \hat { \Pi } _ { \mathcal { S } ^ { \perp } } ^ { 0 } \hat { \rho } ^ { \otimes n } \right ) = \frac { 1 } { \nu } \text{tr} \, \left ( \hat { \Pi } _ { \mathcal { S } ^ { \perp } } ^ { 0 } \hat { \rho } ^ { \otimes n } \right ) = \frac { 1 } { \nu } w ( \mathcal { S } ; \{ W ( \hat { \rho }, i, j ) \} ).$$

In general,

$$W ( \hat { \rho } ^ { \prime } ; u, v ) = \frac { 1 } { \nu } w ( \mathcal { S } + \bar { D } ( u, v ) ; \{ W ( \hat { \rho } ; i, j ) \} ).$$

As a cross-check, note that ∑ u,v w ( S + ¯ ( D u,v ); { y ij } ) = w ( S ⊥ ; { y ij } ) , so the resulting Wigner function is correctly normalized.

The problem of computing the output state of a general qudit magic state distillation routine defined by a stabilizer S is thus reduced to computing the complete weight enumerators of S , and its cosets in S ⊥ .

## 3.2 The strange state and simple-weight enumerators

We now restrict our attention to the special case where ˆ a twirled qutrit strange state ρ | S ⟩ . For this particular magic state, we will show that the complete weight enumerators in the previous subsection can be replaced by simple weight enumerators.

The qutrit strange state is given by

$$| S \rangle = \frac { 1 } { \sqrt { 2 } } \left ( | 1 \rangle - | 2 \rangle \right ).$$

After twirling by symplectic rotations, as described in [14], all noise is converted to depolarizing noise, and noisy | S ⟩ states are described by the density matrix:

$$\hat { \rho } _ { S } ( \epsilon ) = \left ( 1 - \epsilon \right ) | S \rangle \left \langle S \right | + \epsilon \frac { 1 } { 3 }.$$

This corresponds to the discrete Wigner function:

̸

$$W ( \hat { \rho } _ { \epsilon } ; u, v ) = \begin{cases} x & ( u, v ) = ( 0, 0 ) \\ y & ( u, v ) \neq ( 0, 0 ) \end{cases} = \begin{pmatrix} y \ y \ y \\ y \ y \ y \\ x \ y \ y \end{pmatrix},$$

where the parameters x and y satisfy x +8 y = 1, and are related to ϵ via:

$$3 x = - 1 + 4 \frac { \epsilon } { 3 }$$

$$6 y = 1 - \frac { \epsilon } { 3 }. \tag* { \circ }$$

Applying Theorem 1 to ˆ ( ), we find the following lemma. ρ S ϵ

Lemma 1. The probability for successful projection of ˆ( ) ρ ϵ ⊗ n , ν , onto S , is given by the simple weight enumerator of S ⊥ ,

$$\text{tr} \, \left ( \hat { \Pi } _ { \mathcal { S } } \hat { \rho } _ { S } ( \epsilon ) ^ { n } \right ) = w ( \mathcal { S } ^ { \perp } ; x ( \epsilon ), y ( \epsilon ) ),$$

where x , and y are given by eqs. (3.15) and (3.16) .

̸

Let us assume we distilling using a stabilizer code that has a complete set of transversal Clifford gates. The output state ˆ ρ ′ = f MSD (ˆ) of such a distillation protocol will also be ρ of the form in equation (3.14), with parameter ϵ ′ . (Alternatively, if the stabilizer code does not have a complete set of Clifford gates, one could still twirl the output state by symplectic rotations as well to bring it into this form.) Let W ρ (ˆ ; 0 ′ , 0) = x ′ and W ρ (ˆ ; ′ i, j ) = y ′ for ( i, j ) = (0 , 0). Then using equation 3.10, we find,

$$x ^ { \prime } = \frac { w ( \mathcal { S } ; x, y ) } { w ( \mathcal { S } ^ { \perp } ; x, y ) }$$

$$y ^ { \prime } = \frac { 1 } { 8 } \frac { \stackrel { \cdot } { w } ( \mathcal { S } ^ { \perp } ; x, y ) - w ( \mathcal { S } ; x, y ) } { w ( \mathcal { S } ^ { \perp } ; x, y ) }.$$

The formal variables x and y used to define the simple weight enumerator in equation 2.6 are now reinterpreted as entries of the discrete Wigner function in equation 3.14.

Let us rewrite this expression in terms of the noise parameter ϵ . We find,

$$\epsilon ^ { \prime } = 3 \frac { 3 A ( z ( \epsilon ) ) + B ( z ( \epsilon ) ) } { 4 B ( z ( \epsilon ) }.$$

where,

$$z ( \epsilon ) = \frac { 3 - \epsilon } { 8 \epsilon - 6 },$$

A z ( ) = w ( S ; 1 , z ) and B Z ( ) = w ( S ⊥ ; 1 , z ).

The error suppression of the magic state distillation routine is determined by the largest power of ϵ that divides 33 A z ϵ ( ( )) + B z ϵ ( ( )). By contrast, the distance of the error-correcting code is given by the largest power of z that divides B z ( ) -A z ( ).

## 3.2.1 Example: the [[11 , 1 5]] , 3 Golay code

To illustrate the above formalism, let us apply it to the 11-qutrit Golay code of [19]. The 11-qutrit Golay code is an [[11 1 5]] , , 3 CSS-code formed using two copies of the (self-dual) classical ternary Golay code.

The weight enumerator of the 11-qutrit Golay code is computed (e.g., via Magma [41]) to be:

$$A ( z ) = 1 + 5 2 8 z ^ { 6 } + 7 9 2 0 z ^ { 8 } + 1 1 0 0 0 z ^ { 9 } + 2 3 7 6 0 z ^ { 1 0 } + 1 5 8 4 0 z ^ { 1 1 }.$$

Using the MacWilliams identity, we find,

$$B ( z ) = 1 + 5 2 8 z ^ { 5 } + 5 2 8 z ^ { 6 } + 1 5 8 4 0 z ^ { 7 } + 4 0 9 2 0 z ^ { 8 } + 1 2 9 8 0 0 z ^ { 9 } + 1 9 8 0 0 0 z ^ { 1 0 } + 1 4 5 8 2 4 z ^ { 1 1 }. \ ( 3. 2 3 )$$

Substitute these results into equation (3.20), to obtain

$$\sigma ^ { \prime } & = \frac { 3 \left ( 4 8 3 3 6 z ^ { 1 1 } + 6 7 3 2 0 z ^ { 1 0 } + 4 0 7 0 0 z ^ { 9 } + 1 6 1 7 0 z ^ { 8 } + 3 9 6 0 z ^ { 7 } + 5 2 8 z ^ { 6 } + 1 3 2 z ^ { 5 } + 1 \right ) } { 1 4 5 8 2 4 z ^ { 1 1 } + 1 9 8 0 0 0 z ^ { 1 0 } + 1 2 9 8 0 0 z ^ { 9 } + 4 0 9 2 0 z ^ { 8 } + 1 5 8 4 0 z ^ { 7 } + 5 2 8 z ^ { 6 } + 5 2 8 z ^ { 5 } + 1 } \\ & = \left ( 3 0 2 1 e ^ { 1 1 } - 2 4 8 1 6 e ^ { 1 0 } + 9 2 1 8 0 e ^ { 9 } - 2 0 3 2 8 0 e ^ { 8 } + 2 9 2 7 1 0 e ^ { 7 } - 2 8 3 5 3 6 e ^ { 6 } + 1 8 1 7 6 4 e ^ { 5 } - 7 1 2 8 0 e ^ { 4 } + 1 3 3 6 5 e ^ { 3 } \right ) / \\ & \left ( 9 9 0 e ^ { 1 1 } - 7 9 2 0 e ^ { 1 0 } + 2 7 5 0 0 e ^ { 9 } - 5 0 4 9 0 e ^ { 8 } + 3 7 6 2 0 e ^ { 7 } + 4 7 2 5 6 e ^ { 6 } - 1 7 2 6 5 6 e ^ { 5 } + 2 4 3 5 4 0 e ^ { 4 } - 2 0 4 9 3 0 e ^ { 3 } \\ & + 1 0 6 9 2 0 e ^ { 2 } - 3 2 0 7 6 e + 4 3 7 4 \right ) \\ & \approx 5 5 e ^ { 3 } / 1 8 + O ( \epsilon ^ { 4 } ).$$

(3.24)

The threshold of the code, ϵ ∗ , is the critical value of ϵ such that, ϵ &lt; ϵ ∗ implies ϵ ′ &lt; ϵ . Using (3.24), we find

$$\epsilon _ { * } & = \frac { 1 } { 4 5 } \left ( - 2 6 2 \, \sqrt { \frac { 2 } { 4 0 5 \sqrt { 1 0 9 } - 2 9 8 1 } } + 2 ^ { 2 / 3 } \, \sqrt { 4 0 5 \sqrt { 1 0 9 } - 2 9 8 1 } + 3 1 \right ) \\ & \approx 0. 3 8 7.$$

Interestingly, z ∗ = z ϵ ( ∗ ) satisfies a simple cubic equation, 11 z 3 ∗ +12 z 2 ∗ +3 z ∗ +1 = 0.

## 3.3 Conditions for magic state distillation

There are two conditions that a stabilizer code must satisfy for it to qualify as a magic state distillation routine for the strange state. We first require that the probability of successful projection to be nonzero in the limit ϵ → 0, which translates into the requirement

̸

$$B \left ( z ( \epsilon = 0 ) \right ) = B ( - 1 / 2 ) \neq 0.$$

We also require that the noise suppression be better than linear.

Assuming equation (3.26) is satisfied, the noise-suppression exponent, δ , of the magic state distillation routine, ϵ ′ = θ ϵ ( δ ), is determined by the smallest power of ϵ that divides 3 A z ϵ ( ( )) + B z ϵ ( ( )). Let us write

$$3 A ( z ( \epsilon ) ) + B ( z ( \epsilon ) ) = C _ { 0 } + C _ { 1 } \epsilon + C _ { 2 } \epsilon ^ { 2 } + \dots$$

Generically, we expect C 0 and C 1 will be non-zero. The necessary and sufficient condition for C 0 and C 1 to vanish are,

$$( 3 A ( z ( \epsilon ) ) + B ( z ( \epsilon ) ) ) \, | _ { \epsilon = 0 } & = 0 \\ \frac { d } { d \epsilon } | _ { \epsilon = 0 } ( 3 A ( z ( \epsilon ) ) + B ( z ( \epsilon ) ) ) & = 0.$$

Translated into z , these conditions become

$$3 A ( - 1 / 2 ) + B ( - 1 / 2 ) = 0,$$

$$3 A ^ { \prime } ( - 1 / 2 ) + B ^ { \prime } ( - 1 / 2 ) = 0.$$

As a check, observe that, for the weight enumerators of the 11-qutrit Golay code, A ( -1 2) = / 2187 / 64, B ( -1 2) / = 6561 64, / A ′ ( -1 2) / = -8019 / 16 and B ′ ( -1 2) / = 24057 16, / these conditions are satisfied.

The MacWilliams identity simplifies at z = -1 2 to: /

$$B ( - 1 / 2 ) = 3 ( - 1 ) ^ { n } A ( - 1 / 2 ).$$

For n odd, (3.29) is therefore automatically satisfied. We can also ask, when do we get cubic noise suppression? If condition in equation (3.30) is satisfied, the condition for C 2 to vanish is that

$$3 A ^ { \prime \prime } ( - 1 / 2 ) + B ^ { \prime \prime } ( - 1 / 2 ) = 0.$$

Taking derivatives of the MacWilliams identity, we find that,

$$B$$

$$B ^ { \prime } ( - 1 / 2 ) & = ( - 1 ) ^ { n + 1 } \left ( 3 A ^ { \prime } ( - 1 / 2 ) + 8 n A ( - 1 / 2 ) \right ), \\ B ^ { \prime \prime } ( - 1 / 2 ) & = \frac { 1 } { 3 } ( - 1 ) ^ { n } \left ( 9 A ^ { \prime \prime } ( - 1 / 2 ) + 4 8 ( n - 1 ) A ^ { \prime } ( - 1 / 2 ) + 6 4 ( n ^ { 2 } - n ) A ( - 1 / 2 ) \right ). \, ( 3. 3 4 )$$

When n is odd, then the MacWilliams identities above along with (3.30) imply that condition (3.32) is satisfied. Thus, cubic noise suppression is guaranteed for any odd-length stabilizer code whose weight enumerator satisfies the two conditions, (3.26) and (3.30).

## 4 Search for distillation routines

With the above results in place, a computational search for distillation routines for the qutrit strange state is straightforward. For each stabilizer code S in our search space, we compute the simple weight enumerator A z ( ), and then B z ( ) using the MacWilliams identity. We then check if conditions (3.26) and (3.30) are satisfied.

## 4.1 Equivalence of distillation routines

When are two stabilizer codes equivalent for the purposes of magic state distillation?

Conventionally, one considers two quantum error-correcting codes to be equivalent if they differ by permutations of coordinates and by local Clifford operations. However, two stabilizer codes which differ by local Clifford transformations will in general give rise to different magic state distillation protocols.

This observation increases the size of our search space substantially - to search for all magic state distillation routines associated with a given stabilizer code S , we need to search over all orbits of S under local Clifford transformations, and all possible eigenvalues of the stabilizers. For each [[ n, 1]] 3 code, this could increase the search space 9 by a factor as large as, | SL (2 , Z 3 ) | n 3 n -1 = 24 n 3 n -1 .

9 Any local Clifford operator can be written as a product of V ˆ F ˆ ( D χ ) where V ˆ F is a symplectic rotation. Acting on a stabilizer code with a transversal Heisenberg-Weyl displacement operator can only change the

It is possible, however, to substantially reduce this search space size, by demanding that the both the stabilizer codes we study and the magic states we wish to distill possess certain symmetries.

Let the magic state we wish to distill be invariant under a subgroup G of the single-qudit Clifford group. It is natural to require 10 our noisy input states first undergo a twirling procedure (as described in detail for qutrits in [14]) so that they are also invariant under G . Let C be a single-qudit Clifford-unitary, if a twirled noisy magic state ˆ ρ M , satisfies Cρ ˆ M ( C ) † = ˆ ρ M , then two stabilizer codes which differ via local Clifford transformations belonging to the subgroup of the Clifford group generated by by C are equivalent for distillation of the ˆ ρ M magic state. The qutrit magic state | S ⟩ is a simultaneous eigenvector of all symplectic rotations, as emphasized in [14]. Assuming noisy input magic states are twirled by applying a random symplectic rotation, two stabilizer codes are equivalent for magic state distillation if and only if they are related to each other by local symplectic transformations. This reduction in the size of search-space is only possible for the qutrit strange state, by virtue of its exceptional symmetry - and is not expected to be true for any other qutrit or qudit magic state.

It is also natural to require that some, or all, of the unitaries in G be transversal gates for the stabilizer code used for distillation. This is not a sufficient condition for distillation, but it is true for all magic state distillation routines studied in the literature that we are aware of, other than triorthogonal codes (which possess a transversal non-Clifford gate) [24].

In this paper, we are interested in distilling the qutrit magic state | S ⟩ . The symmetry group of | S ⟩ , is the set of all symplectic rotations. Therefore, it is natural to demand that our code possess a complete set of transversal Clifford gates. However, we can instead impose a slightly less restrictive requirement, and demand only that ˆ H 2 be a transversal gate, where ˆ H denotes the single-qutrit Hadamard gate. Notice that, acting with ( ˆ H 2 ) ⊗ n on a stabilizer code corresponds to multiplying each row of the symplectic matrix defining the code by -1, leaving the eigenvalues unchanged. Alternatively, its action can be thought of as changing the eigenvalues ω a of any stabilizer to ω -a . Demanding that ( H 2 ) ⊗ n commute with the codespace is therefore equivalent to demanding that the eigenvalues of all stabilizers in the code must be +1 = ω 0 . This condition allows us to restrict our search to projection onto the trivial syndrome.

Our search space thus consists of two classes of codes: a search over all small quantum error correcting codes, and a search over CSS codes with a complete set of transversal Clifford gates. For both classes, we restrict attention to projection onto the trivial syndrome, due to our requirement that ( ˆ H 2 ) ⊗ n be transversal.

eigenvalues of the stabilizers. Symplectic rotations leave the eigenvalues unchanged but transform the stabilizers by a product of local SL (2 , Z p ) transformations. Therefore, in principle, a search over magic state distillation routines seems to require a search over all local SL (2 , Z 3 ) orbits of each [[ n, 1]] 3 stabilizer codes, with all possible allowed eigenvalues for each of the stabilizers.

10 One, in principle, can also consider magic state distillation without twirling the input states - see [42].

## 4.2 Small quantum error-correcting codes

Any n -qudit stabilizer state can be thought of as a [[ n, 0 , d ]] 3 stabilizer code, and any [ n -1 1] , 3 stabilizer code can be obtained by shortening an [[ n, 0]] 3 stabilizer code. Using the graph-state formalism, [21,22] classified all [[ n, 0 , d ]] 3 stabilizer states up to local Clifford transformations for n ≤ 10, and all 11 12-qutrit stabilizer states of maximum distance 12 6.

For each [[ n, 0]] 3 code, we computed all inequivalent ways of shortening it to obtain an [[ n -1 1]] , 3 stabilizer code. For each resulting [[ n -1 1]] , 3 code, we computed the weight enumerator A z ( ) and the performance for magic state distillation, as described in the previous section. We thus searched over all [[ n, 1 , d ]] 3 codes for m ≤ 9, and all [[11 , 1 , d ]] codes that can be obtained from shortening a [[12 0 6]] , , 3 stabilizer state.

Remarkably the only code we found that could distill the qutrit strange state with betterthan-linear noise-suppression was the 11-qutrit Golay code. There were also a few [[9 1]] , 3 and [11 1]] , 3 codes that distilled the strange state with linear noise suppression, (i.e., ϵ ′ = αϵ with α &lt; 1), all with lower thresholds than the 11-qutrit Golay code.

## 4.3 CSS codes

A maximal self-orthogonal classical ternary code has parameters [ n, ⌊ n 2 ⌋ ] 3 . For n odd, therefore, using the CSS construction [43,44], two copies of a self-orthogonal ternary code of odd length can be combined to form a [[ n, 1 , d ]] 3 qutrit code, whose transversal gate set is the entire single-qutrit Clifford group. Explicitly, if G c is the generator matrix of the self-orthogonal classical ternary code, the quantum CSS code is given by symplectic matrix:

$$H = \begin{pmatrix} G _ { c } & 0 \\ 0 & G _ { c } \end{pmatrix}.$$

Any such CSS code is a natural candidate for a magic state distillation routine for the strange state.

Classical maximal self-orthogonal ternary codes up to size n = 23 have been classified in [45-50], and are conveniently available on a website maintained by Harada and Munemasa [23].

We computed the weight enumerators of all CSS-codes constructed this way from the codes given in [23]. We found that no indecomposable 13 13, 15, 17 or 19-qutrit CSS codes were able to distill the strange state.

Somewhat surprisingly, however, we found a total of 646 inequivalent indecomposable 23-qutrit CSS codes 14 that were able to distill the strange state with cubic noise suppression.

11 Danielsen has also classified all optimal 11-qutrit stabilizer states, but we could not obtain them.

12 The distance of a stabilizer state is the minimum weight of any non-trivial element of its stabilizer group.

13 A code is said to be indecomposable if its generator matrix cannot be written as the direct sum of two smaller generator matrices.

14 The 646 inequivalent CSS codes gave rise to 263 different simple weight enumerators.

Figure 1: A histogram of all the thresholds that arise from 23-qubit CSS codes that are able to distill the strange state. We found a total of 646 codes, and the highest threshold was ϵ ∗ = 0 318. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00436v1-figure-30b4139a-5783-43ab-895b-8ff6eddf63b1.jpg)

A complete list of the 646 classical ternary codes that gave rise to these codes is included in MAGMA format [41] as supplemental data along with this arXiv submission.

There are a total of 1928 indecomposable maximally-self-dual [23 , 11] 3 codes, so the probability that a randomly chosen code will give rise to a quantum CSS code that distills the strange state is 0 335, which is very close to 1 . / 3. The probability that a randomly chosen maximally self-orthogonal [11 , 5] 3 code gives rise to a CSS code that distills the strange state is also 1 / 3. This seems to suggest that quantum error correcting codes that distill the qutrit strange state are actually quite common.

The thresholds that arise from these codes range from 0 063 to 0 318, and are plotted in . . a histogram in Figure 1. None of these thresholds exceed that of the ternary Golay code.

The code with the highest threshold was a [[23 1 5]] , , 3 code. It is formed from two copies of the classical [23 , 11 6] , 3 code whose generator matrix is:

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00436v1-figure-4876c5ab-69b8-4b07-9908-462fa6200e9d.jpg)

The quantum CSS code has weight enumerator

A z ( ) = 2079023616 z 23 +6035662080 z 22 +8258226816 z 21 +7208904960 z 20 +4504066560 z 19 +2182781824 z 18 +790797312 z 17 +252077184 z 16 +52015680 z 15 +14590080 z 14 +2083104 z 13 +628800 z 12 +121824 z 11 +58320 z 10 +16120 z 9 +4608 z 8 +720 z 6 +1 .

(4.3)

It gives rise to a distillation performance

$$\epsilon ^ { \prime } \approx \frac { 7 3 \epsilon ^ { 3 } } { 1 8 } + O \left ( \epsilon ^ { 4 } \right ),$$

plotted in Figure 2.

Figure 2: The distillation performance of the [[23 , 1 5]] , 3 code defined via equation (4.2).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00436v1-figure-d5cc7cd5-ccce-4514-9cb3-95edb23830fa.jpg)

The probability of successful projection onto the codespace for these codes is very low. We refer to the probability for 23 pure copies of the strange state to successfully project onto the codespace as the success probability of the code. The success probabilities range from 1 / 5159780352 ≈ 1 9 . × 10 -10 to 1 / 35831808 ≈ 2 8 . × 10 -8 . The highest probability of successful distillation is attained for the 23-qutrit code in equation (4.2) is 1 / 35831808 - three other codes have the same success probability and very similar thresholds. It appears that the success probability is correlated with the threshold. A plot of success probability versus threshold for the 646 codes that distill the strange state is shown in Figure 3.

## 5 Discussion

Magic state distillation [2,38] is a somewhat mysterious approach to fault-tolerant quantum computing. [2] identified two magic states for qubits, which they denote as | H ⟩ and | T ⟩ . For the qubit | H ⟩ state of [2] and its qudit generalizations defined in [16, 28] - magic state distillation is better understood, and systematic approaches to construct distillation routines

Figure 3: A scatter plot of the success probability (logarithmic scale) and threshold for the 643 CSS codes that distill the strange state.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00436v1-figure-0b6744be-a97e-41ff-a9a2-bec09aaedc3f.jpg)

with low overhead have been proposed in the literature via the formalism of triothogonal codes. However, for the | T ⟩ state of [2] few distillation routines are known [51], and the mechanism behind distillation remains unclear. The qutrit analogue of the | T ⟩ state, which lies directly above a face of the stabilizer polytope, is the strange state | S ⟩ [12, 14] that we study here. The study of distillation of the strange state may in fact be more tractable than that of the qubit | T ⟩ state, thanks to the simple relation between distillation performance and simple enumerators valid only for the qutrit strange state. We hope that, by enlarging the landscape of codes known to distill this state, will lead, in the future, to a better understanding of magic state distillation.

In this work, we carried out a computational search and found that over 600 [[ n, 1 , d ]] 3 codes that distill the qutrit strange state, with cubic noise suppression. This seems to suggest that the ability to distill the qutrit strange state is a generic feature, at least of codes constructed to have a complete set of transversal Clifford gates. While none of these codes have thresholds that exceed that of the 11-qutrit Golay code, their existence does provide some hope that there may exist larger codes with thresholds better than the 11-qutrit Golay code. These results also do provide some support for the belief that a sequence of distillation routines of increasing size whose threshold approaches the limit set by contextuality exists, as required to prove the conjecture that contextuality is sufficient for qutrit quantum computation [1].

## Acknowledgements

SP thanks Prof. P.S. Satsangi for inspiration and guidance. We also thank Mark Howard and Amolak Kalra for comments on a draft of this manuscript. We acknowledge the support of MeitY QCAL, Amazon Braket and DST-SERB grant (CRG/2021/009137).

## Appendices

## A Projecting phase point operators onto a stabilizer code

In this appendix we prove a lemma which is needed for the proof of Theorem 1 in the main text.

Lemma 2. Let S be a group of (conventionally-normalized) commuting n -qudit HeisenbergWeyl operators.

$$\text{tr} \, \left ( \hat { \Pi } _ { \mathcal { S } } ^ { 0 } \hat { A } ( \vec { u }, \vec { v } ) \right ) = \begin{cases} 0 & ( \vec { u } | \vec { v } ) \notin \mathcal { S } ^ { \perp } \\ 1 & ( \vec { u } | \vec { v } ) \in \mathcal { S } ^ { \perp } \end{cases}.$$

Proof. S can be generated by n -k independent Heisenberg-Weyl displacement operators, which we assume have been chosen.

Let

$$\hat { \Pi } _ { \mathcal { S } } ^ { 0 } = \frac { 1 } { p ^ { n - k } } \sum _ { M \in \mathcal { S } } M$$

denote the projector onto the ω 0 eigenspace of each M ∈ S . Let Π ˆ ⃗ s S denote the projector onto a non-trivial subspace, labeled by the vector of syndromes ⃗ s of the n -k generators of S .

Explicitly, we can write

$$\hat { \Pi } _ { \mathcal { S } } ^ { \vec { s } } = \frac { 1 } { p ^ { n - k } } \sum _ { M \in \mathcal { S } } \omega ^ { f ( \vec { s }, M ) } M,$$

̸

for some phases f ⃗ s, M ( ) that depend on the syndrome ⃗ s . Clearly f (0 , M ) = 0. For ⃗ s = 0, then it is easy to see that f ⃗ s, M ( ) takes on all the values j ∈ Z p , i.e., { 0 , 1 , . . . , p - } 1 , and

$$| \{ M \ | f ( \vec { s }, M ) = j \} | = p ^ { n - k - 1 }.$$

Note also that

$$\text{tr} \, \hat { \Pi } _ { \mathcal { S } } ^ { \vec { s } } M = \begin{cases} 0 & M \notin \mathcal { S } \\ p ^ { k } \omega ^ { f ( \vec { s }, M ^ { - 1 } ) } & M \in \mathcal { S } \end{cases}.$$

Therefore

$$\text{tr} \, \left ( \hat { \Pi } _ { \mathcal { S } } ^ { \vec { s } } \hat { A } ( 0, 0 ) ^ { \otimes n } \right ) = \frac { 1 } { p ^ { n } } \sum _ { \vec { u }, \vec { v } } \text{tr} \, \left ( \hat { \Pi } _ { \mathcal { S } } ^ { \vec { s } } \hat { D } ( \vec { u }, \vec { v } ) \right )$$

$$= \frac { 1 } { p ^ { n } } \sum _ { M \in \mathcal { S } } ^ { \cdot } \omega ^ { f ( \vec { s }, M ^ { - 1 } ) } p ^ { k }$$

̸

$$= \begin{cases} \frac { p ^ { k } } { p ^ { n } } \cdot p ^ { n - k - 1 } \left ( 1 + \omega + \omega ^ { 2 } + \dots \omega ^ { p - 1 } \right ) = 0 & \vec { s } \neq 0 \\ 1 & \vec { s } = 0 \end{cases}$$

Note that

$$\hat { D } ( \vec { u }, \vec { v } ) ^ { - 1 } \hat { \Pi } _ { \mathcal { S } } ^ { 0 } \hat { D } ( \vec { u }, \vec { v } ) = \hat { \Pi } _ { \mathcal { S } } ^ { \vec { s } },$$

̸

where the syndrome ⃗ s = 0 if ˆ ( D ⃗ u, ⃗ v ) commutes with all M ∈ S , and ⃗ s = 0 otherwise. The condition that ˆ ( D ⃗ u, ⃗ v ) commutes with all M ∈ S is equivalent to the condition ( ⃗ u, ⃗ v ) ∈ S ⊥ . Therefore,

$$\text{tr } \hat { \Pi } _ { \mathcal { S } } \hat { A } ( \vec { u }, \vec { v } ) = \text{tr } \left ( \hat { D } ( \vec { u }, \vec { v } ) \hat { \Pi } _ { S } \hat { D } ( - \vec { u }, - \vec { v } ) \hat { A } ( 0, 0 ) ^ { \otimes n } \right )$$

̸

$$( A. 1 0 ) \\ ( A. 1 1 )$$

$$= \begin{cases} \text{tr} \, ( \hat { \Pi } _ { S } ^ { \vec { s } \neq 0 } \hat { A } ( 0, 0 ) ^ { \otimes n } ) & ( \vec { u } | \vec { v } ) \notin \mathcal { S } ^ { \perp } \\ \text{tr} \, ( \hat { \Pi } _ { S } ^ { 0 } \hat { A } ( 0, 0 ) ^ { \otimes n } ) & ( \vec { u } | \vec { v } ) \in \mathcal { S } ^ { \perp } \end{cases},$$

and the result follows using equation (A.8).

Let S ⊥ be the dual of S . |S ⊥ | = p n + k . Let us also point out that, if we define,

$$\hat { \Pi } _ { \mathcal { S } ^ { \perp } } ^ { 0 } \equiv \sum _ { u, v } \bar { D } ( u, v ) \hat { \Pi } _ { \mathcal { S } } ^ { 0 } = \frac { 1 } { p ^ { n + k } } \sum _ { ( \vec { u } | \vec { v } ) \in \mathcal { S } ^ { \perp } } \hat { D } ( \vec { u }, \vec { v } ),$$

then one can also show

$$\text{tr} \, \left ( \hat { \Pi } _ { \mathcal { S } ^ { \perp } } ^ { 0 } \hat { A } ( \vec { u }, \vec { v } ) \right ) = \begin{cases} 0 & ( \vec { u } | \vec { v } ) \notin \mathcal { S } \\ 1 & ( \vec { u } | \vec { v } ) \in \mathcal { S } \end{cases}.$$

## B Some codes that do not distill the strange state

It was recently [20] proposed that two qutrit stabilizer codes - a [[13 , 1]] 3 and a [[29 , 1]] 3 CSS code - distill the strange state, and estimated their thresholds to be of 0 425 and . ≈ 0 7, . respectively. We computed the distillation performance of these codes exactly, via their simple weight enumerators, and found that, contrary to the estimates in [20], these two codes do not, in fact, distill the strange state.

The [[13 1 4]] , , 3 CSS-code of [20] generated from two copies of the [13 , 6 3] , 3 maximal self-orthogonal ternary code with generator matrix is, in row-reduced form,

$$M _ { 1 3 } = \begin{pmatrix} 1 & 0 & 0 & 0 & 0 & 0 & 2 & 0 & 2 & 0 \\ 0 & 1 & 0 & 0 & 0 & 0 & 2 & 2 & 1 & 1 & 2 & 0 \\ 0 & 0 & 1 & 0 & 0 & 0 & 1 & 2 & 0 & 1 & 0 & 2 & 2 \\ 0 & 0 & 0 & 1 & 0 & 0 & 1 & 0 & 2 & 0 & 1 & 2 & 2 \\ 0 & 0 & 0 & 0 & 1 & 0 & 0 & 0 & 2 & 2 & 1 & 1 & 0 \\ 0 & 0 & 0 & 0 & 1 & 2 & 1 & 2 & 0 & 1 & 1 & 0 \end{pmatrix}. \quad \begin{pmatrix} & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & && & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & - & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & } & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & + & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & _ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & 0 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ; & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & = & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &, & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ^ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & b & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ( & & & & & & & & & & & & & & & &$$

The weight enumerator of the [[13 1 4]] , , 3 code is

A z ( ) = 1+8 z 3 +600 z 6 +720 z 7 +4320 z 8 +18320 z 9 +61200 z 10 +151200 z 11 +178144 z 12 +116928 z 13 . (B.2)

Figure 4: The distillation performance for the [[13 , 1 4]] , 3 code of [20].

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00436v1-figure-a8576187-877b-49d9-9220-77cb8f39d5b6.jpg)

## This gives rise to a relation

ϵ ′ = 19683 -157464 ϵ +612360 ϵ 2 -1516320 ϵ 3 +2708235 ϵ 4 -3744144 ϵ 5 +4143960 ϵ 6 -3674592 ϵ 7 +2543865 ϵ 8 -1312920 ϵ 9 +471264 ϵ 10 -104384 ϵ 11 +10713 ϵ 12 2 ( 6561 -52488 ϵ +204120 ϵ 2 -495720 ϵ 3 +844425 ϵ 4 -1078704 ϵ 5 +1074888 ϵ 6 -849744 ϵ 7 +529395 ϵ 8 -251400 ϵ 9 +85152 ϵ 10 -18184 ϵ 11 +1827 ϵ 12 ) , (B.3)

plotted in Figure 4. As ϵ → 0, we see that ϵ ′ → 3 2. / The reason for this is that the pure state | S ⟩ ⊗ 13 is orthogonal to the codespace of the stabilizer code, so the code is completely unsuitable for magic state distillation. This can be seen from the fact that B ( -1 2) = 0. This / feature was shared by 6 out of the 7 possible 13-qutrit codes that arise from our construction using the classical ternary codes listed in [23].

The [[29 1 7]] , , 3 CSS code has weight enumerator

A z ( ) =1 + 40 z 6 +4280 z 9 +96 z 10 +2832 z 11 +196584 z 12 +198768 z 13 +1773408 z 14 +15542368 z 15

+91797024

z

16

+565547232

z

17

+3037545272

z

18

+13979050848

z

19

+55970778960

z

20

+192507694176

- +559711606992 z 22 +1361197350960 z 23 +2723501140720 z 24 +4358977591776 z 25
- +5363568387600 z 26 +4767481212256 z 27 +2724627154368 z 28 +751557878400 z 29 .

(B.4)

(Computing the above weight enumerator on the machines we had access to took approximately one day of computation time.)

The resulting relation between ϵ ′ and ϵ is ϵ ′ = 1937 ϵ 224 + O ϵ ( 2 ) , and is shown in Figure 5 below. The threshold for distillation is zero.

Therefore, at present, the 11-qutrit Golay code has the highest threshold for distillation of the strange state, which is still quite far from the threshold set by contextuality [1]. We

z

21

Figure 5: The distillation performance for the [[29 1 7]] , , 3 code of [20]. The threshold for distillation is zero.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00436v1-figure-68150ecb-b7aa-4062-b7de-5b04ad1e1b0a.jpg)

have not yet carried out a search over 29-qutrit codes for distillation, as, to our knowledge a complete list of [29 , 14] 3 self-orthogonal ternary codes is not yet available.

## References

- [1] M. Howard, J. Wallman, V. Veitch, and J. Emerson, 'Contextuality supplies the 'magic' for quantum computation,' Nature 510 (2014), no. 7505 351-355.
- [2] S. Bravyi and A. Kitaev, 'Universal quantum computation with ideal Clifford gates and noisy ancillas,' Phys. Rev. A 71 (2005) 022316.
- [3] N. Delfosse, C. Okay, J. Bermejo-Vega, D. E. Browne, and R. Raussendorf, 'Equivalence between contextuality and negativity of the Wigner function for qudits,' New J. Phys. 19 (2017) 123024.
- [4] J. Haferkamp and J. Bermejo-Vega, 'Equivalence of contextuality and wigner function negativity in continuous-variable quantum optics,' arXiv preprint arXiv:2112.14788 (2021) 2112.14788 .
- [5] D. Gross, 'Hudson's theorem for finite-dimensional quantum systems,' Journal of Mathematical Physics 47 (2006) 122107-122107, quant-ph/0602001 .
- [6] W. K. Wootters, 'A Wigner-function formulation of finite-state quantum mechanics,' Annals of Physics 176 (1987) 1-21.

| [7]   | V. Veitch, C. Ferrie, D. Gross, and J. Emerson, 'Negative quasi-probability as a resource for quantum computation,' New J. Phys. 14 (2012) 113011.                                        |
|-------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [8]   | V. Veitch, S. A. H. Mousavian, D. Gottesman, and J. Emerson, 'The resource theory of stabilizer quantum computation,' New J. Phys. 16 (2014) 013009.                                      |
| [9]   | S. Prakash and A. Gupta, 'Contextual bound states for qudit magic state distillation,' Phys. Rev. A 101 (2020) 010303.                                                                    |
| [10]  | H. Dawkins and M. Howard, 'Qutrit Magic State Distillation Tight in Some Directions,' Phys. Rev. Lett. 115 (2015) 030501.                                                                 |
| [11]  | E. T. Campbell and D. E. Browne, 'Bound states for magic state distillation in fault-tolerant quantum computation,' Phys. Rev. Lett. 104 (2010), no. 3 030503.                            |
| [12]  | W. van Dam and M. Howard, 'Noise thresholds for higher-dimensional systems using the discrete Wigner function,' Phys. Rev. A 83 (2011) 032310.                                            |
| [13]  | P.-E. Emeriau, M. Howard, and S. Mansfield, 'Quantum Advantage in Information Retrieval,' PRX Quantum 3 (2022) 020307.                                                                    |
| [14]  | A. Jain and S. Prakash, 'Qutrit and ququint magic states,' Phys. Rev. A 102 (2020) 042409.                                                                                                |
| [15]  | H. Anwar, E. T. Campbell, and D. E. Browne, 'Qutrit magic state distillation,' New J. Phys. 14 (2012) 063006.                                                                             |
| [16]  | E. T. Campbell, H. Anwar, and D. E. Browne, 'Magic-State Distillation in All Prime Dimensions Using Quantum Reed-Muller Codes,' Phys. Rev. X 2 (2012) 041021.                             |
| [17]  | E. T. Campbell, 'Enhanced Fault-Tolerant Quantum Computing in d -Level Systems,' Phys. Rev. Lett. 113 (2014) 230501.                                                                      |
| [18]  | M. Howard and H. Dawkins, 'Small codes for magic state distillation,' The European Physical Journal D 70 (2016), no. 3 55.                                                                |
| [19]  | S. Prakash, 'Magic state distillation with the ternary Golay code,' Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 476 (2020) 20200187, 2003.02717 . |
| [20]  | A. K. Sharma and S. S. Garani, 'Near-threshold qudit stabilizer codes with efficient encoding circuits for magic-state distillation,' Phys. Rev. A 109 (2024) 062426.                     |
| [21]  | L. Eirik Danielsen, 'Graph-based classification of self-dual additive codes over finite fields,' Advances in Mathematics of Communications 3 (2009), no. 4 329-348.                       |
| [22]  | L. E. Danielsen, 'On the Classification of Hermitian Self-Dual Additive Codes Over GF(9),' IEEE Transactions on Information Theory 58 (2012) 5500-5511.                                   |

- [23] A. Munemasa, 'Database of ternary maximal self-orthogonal codes.' http://www.math.is.tohoku.ac.jp/ ~ munemasa/research/codes/mso3.htm , 2024. Accessed: 2024-07-06.
- [24] S. Bravyi and J. Haah, 'Magic-state distillation with low overhead,' Phys. Rev. A 86 (2012) 052329.
- [25] J. Haah, M. B. Hastings, D. Poulin, and D. Wecker, 'Magic state distillation with low space overhead and optimal asymptotic input count,' Quantum 1 (2017) 31.
- [26] J. Haah and M. B. Hastings, 'Codes and Protocols for Distilling T , controlledS , and Toffoli Gates,' Quantum 2 (2018) 71.
- [27] M. B. Hastings and J. Haah, 'Distillation with Sublogarithmic Overhead,' Phys. Rev. Lett. 120 (2018) 050504.
- [28] E. T. Campbell, 'Enhanced Fault-Tolerant Quantum Computing in d -Level Systems,' Phys. Rev. Lett. 113 (2014) 230501.
- [29] A. Krishna and J.-P. Tillich, 'Towards Low Overhead Magic State Distillation,' Phys. Rev. Lett. 123 (2019) 070507.
- [30] S. Prakash and T. Saha, 'Low Overhead Qutrit Magic State Distillation,' arXiv preprint arXiv:2403.06228 (2024) 2403.06228 .
- [31] A. R. Calderbank, E. M. Rains, P. M. Shor, and N. J. Sloane, 'Quantum error correction via codes over GF (4),' IEEE Transactions on Information Theory 44 (1998), no. 4 1369-1387.
- [32] A. Ketkar, A. Klappenecker, S. Kumar, and P. K. Sarvepalli, 'Nonbinary stabilizer codes over finite fields,' IEEE transactions on information theory 52 (2006), no. 11 4892-4914.
- [33] P. Shor and R. Laflamme, 'Quantum analog of the MacWilliams identities for classical coding theory,' Phys. Rev. Lett. 78 (1997), no. 8 1600.
- [34] K. S. Gibbons, M. J. Hoffman, and W. K. Wootters, 'Discrete phase space based on finite fields,' Phys. Rev. A 70 (2004), no. 6 062101.
- [35] D. M. Appleby, 'Symmetric informationally complete-positive operator valued measures and the extended Clifford group,' Journal of Mathematical Physics 46 (2005) 052107, quant-ph/0412001 .
- [36] X. Wang, M. M. Wilde, and Y. Su, 'Efficiently computable bounds for magic state distillation,' arXiv preprint arXiv:1812.10145 (2018).

| [37]   | N. Koukoulekidis and D. Jennings, 'Constraints on magic state protocols from the statistical mechanics of Wigner negativity,' npj Quantum Information 8 (2022), no. 1 42.                                               |
|--------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [38]   | E. Knill, 'Quantum computing with realistically noisy devices,' Nature 434 (2004) 39-44.                                                                                                                                |
| [39]   | E. T. Campbell and D. E. Browne, 'On the structure of protocols for magic state distillation,' in Workshop on Quantum Computation, Communication, and Cryptography , pp. 20-32, Springer, 2009. 0908.0838 .             |
| [40]   | P. Rall, 'Signed quantum weight enumerators characterize qubit magic state distillation,' arXiv preprint arXiv:1702.06990 (2017).                                                                                       |
| [41]   | W. Bosma, J. Cannon, and C. Playoust, 'The Magma Algebra System I: The User Language,' Journal of Symbolic Computation 24 (1997), no. 3 235-265.                                                                        |
| [42]   | P. Rall, 'Fractal Properties of Magic State Distillation,' arXiv preprint arXiv:1708.09256 (2017) 1708.09256 .                                                                                                          |
| [43]   | A. R. Calderbank and P. W. Shor, 'Good quantum error-correcting codes exist,' Physical Review A 54 (1996), no. 2 1098.                                                                                                  |
| [44]   | A. Steane, 'Multiple-particle interference and quantum error correction,' Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 452 (1996), no. 1954 2551-2577.         |
| [45]   | C. L. Mallows, V. Pless, and N. J. A. Sloane, 'Self-Dual Codes over $ GF( 3 ) $ ,' SIAM Journal on Applied Mathematics 31 (1976), no. 4 649-666, https://doi.org/10.1137/0131058 .                                      |
| [46]   | J. Conway, V. Pless, and N. Sloane, 'Self-dual codes overGF(3)andGF(4)of length not exceeding 16,' IEEE Transactions on Information Theory 25 (1979), no. 3 312-322.                                                    |
| [47]   | V. Pless, N. Sloane, and H. Ward, 'Ternary codes of minimum weight 6 and the classification of the self-dual codes of length 20,' IEEE Transactions on Information Theory 26 (1980), no. 3 305-316.                     |
| [48]   | M. Harada and A. Munemasa, 'A complete classification of ternary self-dual codes of length 24,' Journal of Combinatorial Theory, Series A 116 (2009), no. 5 1063-1072.                                                  |
| [49]   | M. Harada and A. Munemasa, 'On the classification of weighing matrices and self-orthogonal codes,' Journal of Combinatorial Designs 20 (2012), no. 1 40-57, https://onlinelibrary.wiley.com/doi/pdf/10.1002/jcd.20295 . |
| [50]   | M. Araya, M. Harada, and Y. Suzuki, 'Ternary maximal self-orthogonal codes of lengths 21, 22 and 23,' Journal of Algebra Combinatorics Discrete Structures and Applications 5 (2018), no. 1 1-4.                        |

[51] E. Knill, 'Fault-Tolerant Postselected Quantum Computation: Schemes,' 2004.
Shiroman Prakash Rishabh Singhal
A Search for High-Threshold Qutrit Magic State Distillation Routines
Determining the best attainable threshold for qudit magic state distillation
is directly related to the question of whether or not contextuality is
sufficient for universal quantum computation. We carry out a search for
high-threshold magic state distillation routines for a highly-symmetric qutrit
magic state known as the strange state. Our search covers a large class of
$[[n,1]]_3$ qutrit stabilizer codes with up to 23 qutrits, and is facilitated
by a theorem that relates the distillation performance of a qudit stabilizer
code to its weight-enumerators. We could not find any code with $n<23$ qutrits
that distills the strange state with better than linear noise suppression,
other than the 11-qutrit Golay code. However, for $n=23$, we find over 600 CSS
codes that can distill the qutrit strange state with cubic noise suppression.
While none of these codes surpass the threshold of the 11-qutrit Golay code,
their existence suggests that, for large codes, the ability to distill the
qutrit strange state is somewhat generic.