The SIC Question: History and State of Play

TL;DR

This paper reviews recent progress in solving the SIC quantum measurement problem, presenting new solutions up to dimension 151, revealing additional symmetries, and confirming a longstanding conjecture in quantum information theory.

Contribution

It provides new solutions to SIC problems in higher dimensions, demonstrates additional symmetries, and verifies Zauner's conjecture in many cases, expanding the known landscape of SIC solutions.

Findings

Solutions in all dimensions up to 151

Verification of Zauner's conjecture in many cases

Identification of new symmetries in SIC solutions

Abstract

Recent years have seen significant advances in the study of symmetric informationally complete (SIC) quantum measurements, also known as maximal sets of complex equiangular lines. Previously, the published record contained solutions up to dimension 67, and was with high confidence complete up through dimension 50. Computer calculations have now furnished solutions in all dimensions up to 151, and in several cases beyond that, as large as dimension 844. These new solutions exhibit an additional type of symmetry beyond the basic definition of a SIC, and so verify a conjecture of Zauner in many new cases. The solutions in dimensions 68 through 121 were obtained by Andrew Scott, and his catalogue of distinct solutions is, with high confidence, complete up to dimension 90. Additional results in dimensions 122 through 151 were calculated by the authors using Scott's code. We recap the history of the problem, outline how the numerical searches were done, and pose some conjectures on how the search technique could be improved. In order to facilitate communication across disciplinary boundaries, we also present a comprehensive bibliography of SIC research.

Figures

Equations

$${\left|⟨{\psi }_j|{\psi }_k⟩\right|}^2=\frac{1}{d+1}$$

$$d=2–24,28,30,31,35,37,39,43,48,124.$$

$$A||u|{|}^2\le \sum _j{\left|⟨v_j,u⟩\right|}^2\le B||u|{|}^2,$$

$$⟨A,B⟩=tr(AB).$$

$${\Pi }_j=\left|{\psi }_j⟩⟨{\psi }_j\right|,$$

$$tr({\Pi }_j{\Pi }_k)=\frac{d{\delta }_{jk}+1}{d+1},$$

$$Q_j^{\pm }=\pm \sqrt{d+1}{\Pi }_j+\frac{1\mp \sqrt{d+1}}{d}I.$$

$$X|j⟩=|j+1⟩,Z|j⟩={\omega }^j|j⟩,$$

$$X^l{Z}^{\alpha }={\omega }^{-l\alpha }{Z}^{\alpha }{X}^l.$$

$$D_{l\alpha }=(-e^{i\pi /d}{)}^{l\alpha }{X}^l{Z}^{\alpha }.$$

$$D_{l\alpha }{D}_{m\beta }=(-e^{i\pi /d}{)}^{\alpha m-\beta l}{D}_{l+m,\alpha +\beta }.$$

$${\left|\psi_{0}^{\rm(qubit)}\right\rangle}=\frac{1}{\sqrt{6}}\left(\begin{array}[]{c}\sqrt{3+\sqrt{3}}\\ e^{i\pi/4}\sqrt{3-\sqrt{3}}\end{array}\right).$$

$${\left|\psi_{0}^{(\rm Hesse)}\right\rangle}=\frac{1}{\sqrt{2}}\left(\begin{array}[]{c}0\\ 1\\ -1\end{array}\right)$$

$$|⟨{\hat{v}}_i,{\hat{v}}_j⟩|=\alpha \forall i\ne j,$$

$$|{\psi }_0⟩=\sum _j{a}_j|j⟩,$$

$$[\mathbf{F}{]}_{\beta l}={\left|⟨{\psi }_0|X^l{Z}^{\beta }|{\psi }_0⟩\right|}^2={\left|\sum _j{a}_j^{\ast }{a}_{j+l}{\omega }^{-j\beta }\right|}^2.$$

$$f_j^{(l)}=a_j^{\ast }{a}_{j+l}.$$

$$[\mathbf{G}{]}_{kl}=(f^{(l)}\star f^{(l)}{)}_k=\sum _j{a}_j{a}_{j+k}^{\ast }{a}_{j+l}^{\ast }{a}_{j+k+l}.$$

$$[\mathbf{G}{]}_{kl}=[\mathbf{G}{]}_{lk}.$$

$$[\mathbf{F}{]}_{lm}=\frac{1}{d}\sum _{\alpha ,\beta }{\omega }^{-\alpha l+\beta m}[\mathbf{F}{]}_{\beta \alpha }.$$

$$\mathcal{F}{\left(\{[\mathbf{G}{]}_{kl}\}\right)}_{\beta }=\sum _k{\omega }^{-k\beta }[\mathbf{G}{]}_{kl},$$

$$[\mathbf{G}{]}_{kl}=\sum _j{a}_j{a}_{j+k}^{\ast }{a}_{j+l}^{\ast }{a}_{j+k+l}=\frac{{\delta }_{k0}+{\delta }_{l0}}{d+1}.$$

$$\sum _{k,l}|[\mathbf{G}{]}_{kl}{|}^2=\frac{2}{d+1}.$$

$$\sum _{k,l}|[\mathbf{G}{]}_{kl}{|}^2\ge \frac{2}{d+1},$$

$$\left|[\mathbf{G}{]}_{k,l}\right|=\left|[\mathbf{G}{]}_{-k,l}\right|=\left|[\mathbf{G}{]}_{k,-l}\right|=\left|[\mathbf{G}{]}_{-k,-l}\right|=\left|[\mathbf{G}{]}_{l,k}\right|=\left|[\mathbf{G}{]}_{-l,k}\right|=\left|[\mathbf{G}{]}_{l,-k}\right|=\left|[\mathbf{G}{]}_{-l,-k}\right|.$$

$$\left|[\mathbf{G}{]}_{kl}\right|\le [\mathbf{G}{]}_{0l}.$$

$$[\mathbf{G}{]}_{-k,l}=[\mathbf{G}{]}_{k,l}^{\ast }.$$

$$[\mathbf{G}{]}_{0k}$$

$$[\mathbf{G}{]}_{1k}$$

$$[\mathbf{G}{]}_{2k}$$

$$X^m{Z}^n\mapsto X^{-n}{Z}^{m-n},$$

Full text

I Introduction

The problem of symmetric, informationally complete quantum measurements Zauner:1999 ; Caves:2002 ; Renes:2004 ; Scott:2010a stands at the confluence of multiple areas of physics and mathematics. SICs, as they are known for short, tie into algebraic number theory Appleby:2016 ; Bengtsson:2016 ; FOOP:2017 ; Appleby:2017 , higher-dimensional sphere packing Stacey:2016 , Lie and Jordan algebras Appleby:2011a ; Appleby:2015 , finite groups Zhu:2015a ; Stacey:2016b and quantum information theory Fuchs:2010 ; Tabia:2012 ; Tabia:2013 ; Fuchs:2013 ; Zhu:2012 ; Graydon:2016 ; Zhu:2016a ; Szymusiak:2015 ; Appleby:2016b ; Fuchs:2016a . Without the study of SICs, one might think that the intersection of all these subjects would have to be the empty set. And yet, for all that, a SIC is a remarkably simple mathematical structure, as structures go. Consider the complex vector space $\mathbb{C}^{d}$. To a physicist, this is the Hilbert space associated with a $d$-level quantum system. Let $\{{\left|\psi_{j}\right\rangle}\}$ be a set of exactly $d^{2}$ unit vectors in $\mathbb{C}^{d}$ such that

[TABLE]

whenever $j\neq k$. The set $\{{\left|\psi_{j}\right\rangle}\}$, which can be associated with a set of pairwise equiangular lines through the origin, is a SIC.

One can prove that no more than $d^{2}$ vectors in a $d$-dimensional Hilbert space can be equiangular. That is, if $\{{\left|\psi_{j}\right\rangle}\}$ is a set of vectors, and $\left|{\left\langle\psi_{j}\middle|\psi_{k}\right\rangle}\right|^{2}=\alpha$ for every $j\neq k$, then that set can have at most $d^{2}$ elements. In addition, for a maximal set the value of $\alpha$ is fixed by the dimension; it must be $1/(d+1)$. So, a SIC is a maximal equiangular set in $\mathbb{C}^{d}$; the question is whether they can be constructed for all values of the dimension. Despite a substantial number of exact solutions, as well as a longer list of high-precision numerical solutions Scott:2010a ; Scott:2017 ; Appleby:2017 , the problem remains open.

Exact solutions, found by hand in a few cases and by computer algebra software in the others, are known in the following dimensions:

[TABLE]

The historical record of exact solutions has been spread over several publications Scott:2010a ; Chien:2015b ; Appleby:2017 ; Grassl:2017 . For several years, the most extensive published set of numerical results went as high as dimension $d=67$ Scott:2010a . Now, numerical solutions are known in all dimensions up to and including $d=151$, as well as a handful of other dimensions up to $d=844$. These numerical solutions were found using code designed and written by Scott, who extended the results of Scott:2010a through $d=121$ using his personal computer over several years of dedicated effort. In addition, Scott found solutions in a set of dimensions ($d=124$, 143, 147, 168, 172, 195, 199, 228, 259, 323) by taking advantage of particular simplifying assumptions that are applicable in those dimensions Scott:2017 . Further close study of these properties led to a solution for $d=844$ Grassl:2017 . Because dimension $d=121$ was pushing the limits of what was computationally feasible without those simplifying assumptions, the authors calculated solutions in dimensions 122 through 151 by running Scott’s code on the Chimera supercomputer at UMass Boston. In turn, Scott was able to employ another algorithm (outlined below) to refine the numerical precision of these results.

The solutions from all of these search efforts are available together at the following website:

http://www.physics.umb.edu/Research/QBism/

An intriguing feature of the SIC problem is that some numerical solutions, if extracted to sufficiently high precision, can be converted to exact ones Chien:2015b ; Appleby:2017 . Most recently, this technique was used to derive an exact solution in dimension $d=48$. Another interesting aspect is that the number of distinct SIC constructions varies from one dimension to another. (The sense in which two SICs can be equivalent will be discussed in detail below.) One reason computational research is valuable, beyond extending the list of dimensions in which SICs are known, is that it provides what is likely a complete picture for many values of the dimension. This is important for understanding the subtle connection between SICs and algebraic number theory Appleby:2016 , a connection that brings a new angle of illumination to Hilbert’s twelfth problem Bengtsson:2016 ; FOOP:2017 .

SICs are so called because, thanks to the rules of quantum theory, a SIC in $\mathbb{C}^{d}$ specifies a measurement procedure that can, in principle, be applied to a $d$-level quantum system. For example, a SIC in $\mathbb{C}^{2}$ is a set of four equiangular lines, and it is a mathematical model of a measurement that a physicist can perform on a single qubit. The term “informationally complete”—the “IC” in “SIC”—means that if one has a probability distribution for the possible outcomes of a SIC experiment, one can compute the probabilities for the possible outcomes of any other experiment carried out on the target system Fuchs:2013 . So, while one can pose the question of their existence using pure geometry, SICs are relevant to applied physics. Indeed, SIC measurements have recently been performed or approximated in the laboratory Du:2006 ; Durt:2008 ; Medendorp:2011 ; Pimenta:2013 ; Bian:2014 ; Zhao:2015 ; Bent:2015 ; SosaMartinez:2017 , and they are known to be optimal measurements for quantum-state tomography Scott:2006 .

A SIC provides a frame—more specifically, an equiangular tight frame—for the vector space $\mathbb{C}^{d}$. Given a finite-dimensional Hilbert space $\mathcal{H}$ with an inner product $\left\langle\cdot,\cdot\right\rangle$, a frame for $\mathcal{H}$ is a set of vectors $\{v_{j}\}\subset\mathcal{H}$ such that for any vector $u\in\mathcal{H}$,

[TABLE]

for some positive constants $A$ and $B$. The frame is equal-norm if all the vectors $\{v_{j}\}$ have the same norm, and the frame is tight if the “frame bounds” $A$ and $B$ are equal. The ratio of the number of vectors to the dimension of the space is known as the redundancy of the frame Et-Taoui:2016 . For more on this terminology and its history, we refer to Kovačević and Chebira Kovacevic:2007 ; Kovacevic:2007b . In our experience, the language of frames is more common among those who come to SICs from pure mathematics or from signal processing than among those motivated by quantum physics.

Any vector in $\mathbb{C}^{d}$ can be represented by its inner products with all the SIC vectors. In quantum physics, one also considers the set of Hermitian operators on $\mathbb{C}^{d}$. This set in fact forms a Hilbert space itself, with a dimension of $d^{2}$, and the inner product given by the Hilbert–Schmidt formula

[TABLE]

Rewriting the SIC vectors $\{{\left|\psi_{j}\right\rangle}\}$ as rank-1 projection operators,

[TABLE]

we construct a nonorthogonal basis for the Hilbert space of Hermitian operators. Because the inner products of these projectors are uniform, given by

[TABLE]

then it is straightforward to find a shifting and rescaling that orthogonalizes the basis $\{\Pi_{j}\}$, at the cost of making the operators non-positive-semidefinite. In fact, there are two choices:

[TABLE]

The bases $\{Q_{j}^{\pm}\}$ have interesting properties with regard to Lie algebra theory Appleby:2015 and the study of quantum probability Zhu:2016a ; Fuchs:2017 .

II Generating SICs with Groups

All known SICs have an additional kind of symmetry, above and beyond their definition: They are group covariant. Each SIC can be constructed by starting with a single vector, known as a fiducial vector, and acting upon it with the elements of some group. It is not known in general whether a SIC must be group covariant. Because such an assumption greatly reduces the search space Scott:2010a ; Appleby:2016 , it has been the only method used so far: The fact that we only know of group-covariant SICs could potentially be an artifact of this. (However, we do have a proof that all SICs in $d=2$ and $d=3$ are group covariant Hughston:2016 .)

In all cases but one, the group that generates a SIC from a fiducial is an instance of a Weyl–Heisenberg group. We can define this group as follows. First, fix a value of $d$, and let $\omega=e^{2\pi i/d}$. Let $\{{\left|0\right\rangle},{\left|1\right\rangle},\ldots,{\left|d-1\right\rangle}\}$ be an orthonormal basis for the Hilbert space $\mathcal{H}_{d}=\mathbb{C}^{d}$. Then, construct the shift and phase operators

[TABLE]

where the shift is modulo $d$. These operators satisfy the Weyl commutation relation,

[TABLE]

In a sense, the operators $X$ and $Z$ come as close as possible to commuting, without actually doing so: The only cost to exchanging their order is a phase factor determined by the dimension.

The Weyl–Heisenberg displacement operators in dimension $d$ are defined by

[TABLE]

The product of two displacement operators is, up to a phase factor, a third:

[TABLE]

Therefore, by allowing the generators to be multiplied by phase factors, we can define a group, known as the Weyl–Heisenberg group in dimension $d$. This group dates back to the early days of quantum physics. Weyl introduced the generators $X$ and $Z$ as long ago as 1925 in order to define what one might mean by the quantum theory of discrete degrees of freedom Weyl:1931 ; Scholz:2006 ; Scholz:2007 (see also (Fuchs:2014, , pp. 2055–56)). This group, and structures derived from it, are critically important in quantum information and computation; for example, this is the basic prerequisite for the Gottesman–Knill theorem, which indicates when a quantum computation can be efficiently simulated classically Gottesman:1999 . The close relationship between SICs and the Weyl–Heisenberg group suggests that SICs are a kind of structure that quantum physics should have been studying all along.

Zhu has proved that in prime dimensions, group covariance implies Weyl–Heisenberg covariance Zhu:2010 . The one known exception to the rule of Weyl–Heisenberg covariance is the Hoggar SIC Hoggar:1981 ; Hoggar:1998 , which lives in a prime-power dimension, $d=8$. As in all other dimensions, there is a Weyl–Heisenberg SIC, but there is also the Hoggar SIC. Like many other exceptions to mathematical classifications, it is related to the octonions Stacey:2016 ; Stacey:2016b .

One example of a Weyl–Heisenberg SIC can be constructed by taking the orbit of the following two-dimensional vector under the Weyl–Heisenberg displacement operators:

[TABLE]

This orbit is a set of four vectors. In the Bloch sphere representation, they form the vertices of a regular tetrahedron inscribed within the sphere.

An example in dimension $d=3$, one which is remarkable for the further subtle symmetries it possesses beyond even group covariance, is the orbit of

[TABLE]

under the Weyl–Heisenberg displacements. This set of vectors is known as the Hesse SIC Hughston:2007 ; Dang:2013 ; Hughston:2016 , thanks to its relation with the Hesse configuration familiar from design theory and the study of cubic curves, specifically nonsingular cubic curves in complex projective two-space Coxeter:1991 ; Tabia:2013 ; Zhu:2015a ; Stacey:2016c .

III Historical Overview

In order to understand the current state of SIC research, one must grasp how people came to the SIC question, what other structures they think are related, what tools they suspect are applicable, and so forth. A physicist, motivated by quantum information theory, is apt to have a different mental context than a pure mathematician driven by the abstract appeal of geometry. To attempt to foster an interdisciplinary discussion, we provide in this section a brief historical overview.

The Hesse SIC can be extracted from a 1940 article by H. S. M. Coxeter Coxeter:1940 , which discusses what he later called the “Hessian polyhedron” Coxeter:1991 . This polyhedron lives in three-dimensional complex vector space and has 27 vertices (which correspond to the 27 lines on a cubic surface Baez:2016 ). The vectors in the Hesse SIC are a subset of those vertices. Mathematicians were studying SICs in $d=2$ and $d=3$ more explicitly as early as the 1970s, in the terminology of complex equiangular lines Delsarte:1975 . Stuart Hoggar found a SIC in dimension $d=8$ a few years later Hoggar:1981 , by considering the diagonals of a quaternionic polytope and converting their coordinates to complex numbers. The SICs in $d=2$ and $d=3$, together with the Hoggar lines in $d=8$, still stand out among the known SICs; various unusual attributes they possess have led them to be designated the sporadic SICs Stacey:2016 ; FOOP:2017 .

In a 1987 article, Richard Feynman used a construction that is in retrospect a $d=2$ SIC to study the probability theory of a qubit. SICs entered quantum theory more generally starting with the work of Gerhard Zauner, who began to consider the problem in the 1990s (Fuchs:2014, , p. 1941). By 1999, Zauner had found the connection with the Weyl–Heisenberg group and proven SIC existence up to $d=5$ Zauner:1999 . He also posed a conjecture that the search for Weyl–Heisenberg SICs could be simplified by considering a particular unitary operator Zauner:1999 , a conjecture we will describe in detail below.

Independently of Zauner, Carlton Caves developed the idea of a SIC as representing a quantum measurement Caves:2002 , motivated originally by attempts to prove the “quantum de Finetti theorem” Caves:2002c . SICs turned out to have more symmetry than was necessary for that proof, but they soon took on a life of their own in quantum information theory Fuchs:2016a . The term “SIC” (first pronounced as “sick,” and later like “seek”) dates to this period. A 2004 paper, “Symmetric Informationally Complete Quantum Measurements” Renes:2004 , introduced them to the mainstream of the quantum information community, and reported numerical solutions for dimensions up to $d=45$. In a later survey, Andrew Scott and Markus Grassl extended these numerical results up to dimension $d=67$ Scott:2010a . Further work by Scott established Weyl–Heisenberg SICs in all cases up to $d=121$ without exception. Moreover, these results together are probably a complete catalogue of all distinct Weyl–Heisenberg SICs up to dimension $d=90$ Scott:2017 . In later sections, we will give an overview of how these searches were done.

Exact solutions were found more slowly. Having exact expressions for SIC fiducial vectors allowed Appleby et al. to discover a connection with Galois theory Appleby:2013 . In turn, this led to further relations with algebraic number theory, a frankly mysterious development that is still under active investigation Appleby:2016 ; FOOP:2017 ; Bengtsson:2016 ; Appleby:2017 .

An exhaustive treatment of SIC research would require transgressing many disciplinary boundaries, and would doubtless grow to an intimidating length. In order to provide at least a helpful ordering on the literature, we have made our bibliography as comprehensive as possible. We note in particular a selection of papers that address SICs in the context of abstract algebra Appleby:2011a ; Zhu:2014 ; Appleby:2015 ; Appleby:2017 , algebraic number theory Appleby:2013 ; Appleby:2016 ; FOOP:2017 ; Bengtsson:2016 , category theory Wetering:2017 , finite group theory Zhu:2010 ; Appleby:2012b ; Bengtsson:2012b ; Zhu:2015a ; Szymusiak:2015 ; Stacey:2016 ; Stacey:2016b ; Chien:2017 , quantum computing and contextuality Bengtsson:2012 ; Cabello:2012 ; Zhu:2016a ; Stacey:2016b ; Fuchs:2016b ; Howard:2016 ; Fuchs:2017 ; Andersson:2017 , and quantum entanglement Chen:2015 ; Shen:2015 ; Graydon:2016 ; Graydon:2016a . We also note the significant number of student theses written in whole or in part on the SIC problem Zauner:1999 ; Renes:2004a ; MKF:2008 ; YadsanAppleby:2012 ; Zhu:2012 ; Tabia:2013b ; Andersson:2014 ; Blanchfield:2014 ; Chien:2015b ; Dang:2015 ; Stacey:2015 .

Before moving on, we note that the real analogue of the SIC problem, i.e., finding maximal sets of equiangular lines in real vector spaces, has also been of considerable interest to mathematicians Sloane:2015 ; Greaves:2016 ; Balla:2016 . The maximal number of equiangular lines in a $d$-dimensional vector space is not $d^{2}$, but only $d(d+1)/2$. That is, if we have a set of $N$ unit vectors $\{\hat{v}_{i}\}$ in a $d$-dimensional vector space, such that

[TABLE]

then the size $N$ of the set cannot exceed $d(d+1)/2$. Moreover, while the complex bound of $d^{2}$ has been saturated in every dimension that we have been able to check, it is known that the real bound of $d(d+1)/2$ is not even attained for all values of $d$. For example, in $d=7$, one can construct a set of $7\cdot 8/2=28$ equiangular lines, but this is also the best that can be done in $d=8$. In fact, the only known instances where the bound of $d(d+1)/2$ can be attained are dimensions 2, 3, 7 and 23 Balla:2016 .

There is a sign freedom in this definition of the angle, since Eq. (14) is satisfied if the inner product $\langle\hat{v}_{i},\hat{v}_{j}\rangle$ is either $+\alpha$ or $-\alpha$. The presence of this discrete choice means that investigations of real equiangular lines often have a rather combinatorial flavor. In contrast, when we take the magnitude of a complex inner product, we discard a continuous quantity, a phase that in principle can be anywhere from [math] to $2\pi$. Generally speaking, the “feel” of the real and complex problems differ, as is evidenced by the different areas of mathematical expertise brought to bear upon them. However, subtle and unanticipated points of contact between the real and complex cases do exist Stacey:2016b .

IV How to Search for SICs Numerically

As before, let $\{{\left|j\right\rangle}\}$ be an orthonormal basis for the Hilbert space $\mathcal{H}_{d}=\mathbb{C}^{d}$. In this basis, the fiducial vector can be written

[TABLE]

for some set of coefficients $\{a_{j}\}$.

Acting with the Weyl–Heisenberg operator $D_{l\beta}$ on the fiducial vector ${\left|\psi_{0}\right\rangle}$ produces a new vector, whose squared innner product with the fiducial vector is

[TABLE]

The right-hand side has the form of the magnitude squared of a Fourier coefficient, i.e., of a power spectrum. Specifically, the set of squared inner products between ${\left|\psi_{(l,\beta)}\right\rangle}$ and ${\left|\psi_{0}\right\rangle}$ for any given value of $l$ is the power spectrum of the sequence

[TABLE]

By the Wiener–Khinchin theorem, we know that the power spectrum of a sequence is the Fourier transform of the autocorrelation of that sequence Wiener:1964 . Therefore, the autocorrelation of the sequence $f^{(l)}_{j}$, when put through the Fourier transform, will yield the sequence $[{\bf F}]_{\beta l}$. The set of autocorrelation sequences for all values of $l$ forms a matrix. Using $\star$ to denote the correlation of two sequences, we can write the elements of this matrix as

[TABLE]

The matrix ${\bf G}$ is in many situations more convenient to work with than the original matrix ${\bf F}$, because ${\bf G}$ lacks phase factors and treats both of its indices on equal footing. For example, it is apparent from the definition of ${\bf G}$ that

[TABLE]

This is equivalent to a property of the matrix ${\bf F}$ that is less obvious to the eye:

[TABLE]

If we take the Fourier transform of the columns of the matrix ${\bf G}$,

[TABLE]

we recover the squared inner products between the candidate SIC vectors and the fiducial. This means that if the vectors $\{{\left|\psi_{(l,\beta)}\right\rangle}\}$ really do comprise a SIC, then the matrix ${\bf G}$ must take a very specific form. Every entry in $\mathcal{F}\left(\{[{\bf G}]_{kl}\}\right)_{\beta}$ must equal $1/(d+1)$, except for the element at $l=\beta=0$, which equals 1. Recalling that a constant sequence is the discrete Fourier transform of a Kronecker delta function, we can deduce the desired values of $[{\bf G}]_{kl}$.

The result is that if ${\left|\psi\right\rangle}$ is a Weyl–Heisenberg fiducial vector, then

[TABLE]

This implication also works in reverse, thanks to the transitivity of the group action.

The basic idea of finding SICs numerically is to use standard optimization methods to find a fiducial vector that makes $[{\bf G}]_{kl}$ as close to the desired form as possible. Note that, when ${\bf G}$ is constructed from a SIC fiducial,

[TABLE]

One can prove Appleby:2014b that this is a lower bound. In general,

[TABLE]

and the inequality is saturated if and only if the input vector is truly a SIC fiducial.

This naturally suggests a way to find SIC fiducial vectors: Minimize the LHS of the inequality in Eq. (24), aiming for the lower bound given on the RHS. During our time investigating the SIC question, we have at various points implemented this idea in Mathematica, in Python and in C++ using the GNU Scientific Library. We find in general that numerical optimization finds a local minimum quickly, but a local minimum might only imply inner products between the vectors that are correct to a few decimal digits. A way around this problem is to repeat the optimization many times, starting from different points in the search space. Since these trials can run concurrently, the problem is amenable to parallelization. This is the approach we followed when using the Chimera supercomputer to obtain solutions in dimensions 122 through 151. Scott’s implementation, which we employed on Chimera, uses a C++ code for a limited-memory quasi-Newton optimization algorithm, L-BFGS, due to Liu Liu:2012 .

As is evident from Figure 1, the time required to obtain solutions did not increase steadily with the dimension. For example, $d=146$ took eleven days of computer time and $d=148$ required twelve days, but $d=147$ took only 18 hours. Likewise, Chimera spent 28 days trying to find a $d=151$ solution before succeeding, but it found a SIC in $d=150$ in only two hours. (These figures are all for “wall clock” elapsed time. The number of processor-hours devoted was greater, since we ran Scott’s code in parallel on 96 of Chimera’s cores.) We suspect the variation is due to different numbers of inequivalent solutions existing in different dimensions: The more solutions, the easier it is to hit upon one of them.

Once we have a numerical result in hand, we can refine its precision. This requires a code that uses multi-precision arithmetic, which will run more slowly than the optimization in the first step Scott:2017 . The fiducial vectors available at Scott:2017 and at the website referenced above were obtained in this way and are accurate to 150 digits.

Before moving on, we note a conjecture, based on numerical evidence, that hints at additional hidden structure in the SIC problem. Note that the definition of ${\bf G}$ implies

[TABLE]

Here, indices are to be interpreted modulo $d$. Because any autocorrelation attains its maximum at zero offset, we also know immediately that the elements of ${\bf G}$ cannot be larger in the middle of the matrix than they are on the edges:

[TABLE]

The Fourier transform of $[{\bf G}]_{kl}$ over the index $k$ is, by the Wiener–Khinchin theorem, the power spectrum of $f^{(l)}$. Because power spectra are nonnegative, we can say that

[TABLE]

Are there additional symmetries or redundancies, not so apparent from the definition? By happenstance, one of the authors (CAF) observed that imposing a subset of the constraints in Eq. (22) was sufficient to find a SIC fiducial vector (Fuchs:2014, , pp. 1252–59). Specifically, by finding a solution to the simultaneous equations

[TABLE]

one finds a solution to all the equations in (22). The redundancies in Eq. (25) are sufficient to imply that this holds up to $d=5$. We call the idea that it remains true in all dimensions the “$3d$ conjecture.” It has been verified numerically up to dimension $d=28$ Appleby:2014b . If the $3d$ conjecture is indeed true, it would reduce the complexity of the problem, as measured by the number of simultaneous equations to solve, from quadratic in the dimension to linear.

V Zauner Symmetry

Is there any way to narrow the search space for SIC fiducial vectors? To see how to answer this in the affirmative, we must elaborate upon the group theory we discussed in the previous sections. The Clifford group for dimension $d$ is the “normalizer” of the Weyl–Heisenberg group: It is the set of all unitary operators that, acting by conjugation, map the set of Weyl–Heisenberg operators in dimension $d$ to itself. We saw earlier how the orbit of a vector under the Weyl–Heisenberg group can be a SIC; likewise, we can study the orbit of a vector under the entire Clifford group.

For our purposes in this note, the important point is that if we conjugate a Weyl–Heisenberg operator $D_{kl}$ by a Clifford unitary $U$, we obtain a Weyl–Heisenberg operator $D_{k^{\prime}l^{\prime}}$, possibly with an additional phase. Details on the construction and representation of the Clifford group in any finite dimension $d$ can be found in Appleby Appleby:2005b .

It was conjectured by Zauner Zauner:1999 , and independently by Appleby, that in every dimension $d>2$, a Weyl–Heisenberg SIC fiducial exists that is an eigenvector of a certain order-3 Clifford unitary, which is now known as the Zauner unitary. Acting on the Weyl–Heisenberg generators, the Zauner unitary effects the change

[TABLE]

up to an overall phase factor. (See equation (3.10b) of Zauner Zauner:1999 , or equation (127) of Appleby Appleby:2005b .) Applying this again yields

[TABLE]

and a third iteration gives

[TABLE]

confirming that this operation has order 3.

How might assuming the Zauner conjecture simplify the search for SICs? First, we will make some remarks on this from an algebraic perspective, and then we will address the point in a way suited to numerical optimization. Let ${\left|\psi\right\rangle}$ be a candidate fiducial vector, and suppose that it is an eigenvector of the Zauner unitary $U$ with unit eigenvalue:

[TABLE]

Consequently,

[TABLE]

As $U$ is a Clifford unitary, requiring that ${\left|\psi\right\rangle}$ is an eigenvector of $U$ implies degeneracies among the elements of the matrix ${\bf F}$.

Because

[TABLE]

the Zauner unitary sends the left edge of ${\bf F}$ to the top edge and then to the main diagonal. More generally, specifying a column of ${\bf F}$ (which is equivalent to fixing a column of ${\bf G}$) and imposing the Zauner condition means that the same constraint also simultaneously fixes a row and a diagonal.

Earlier, in Eq. (20), we saw that a symmetry of ${\bf G}$ implied a Fourier-type relation among the elements of ${\bf F}$. We have expressed the Zauner condition as a degeneracy within ${\bf F}$, but what does it imply for ${\bf G}$? The result can be found straightforwardly, and it closely resembles Eq. (20):

[TABLE]

This the expression of Zauner symmetry in the ${\bf G}$ matrix. A special case of note: If we set $k=l=0$, then

[TABLE]

We note that the assumption that the fiducial is a Zauner eigenvector is enough to prove some additional cases of the $3d$ conjecture, up to dimension $d=9$, in a straightforward way. In dimension $d=5$, the basic symmetries of ${\bf G}$ imply that the $3d$ constraints automatically specify the entire matrix ${\bf G}$, and thus also fix ${\bf F}$ to have the desired form for a SIC. By imposing the condition that our initial vector is a Zauner eigenvector, we can extend this up to dimension $d=8$. This can be seen directly by drawing an $8\times 8$ grid and shading in the appropriate squares.

In fact, we can carry this argument a little further. We obtain ${\bf G}$ by Fourier transforming the columns of ${\bf F}$. Therefore, if we specify $\{[{\bf G}]_{k0}\}$, we automatically fix $\{[{\bf F}]_{k0}\}$, and by Eq. (25), we set the values of $\{[{\bf G}]_{0k}\}$ as well. Imposing the Zauner condition fixes $\{[{\bf F}]_{0k}\}$ in terms of $\{[{\bf F}]_{k0}\}$. Specifically, for $k\neq 0$, we have

[TABLE]

Recalling that

[TABLE]

we therefore find that

[TABLE]

In other words, the Zauner condition implies that if we add up the entries in a column, leaving out the entry on the top row, they must all cancel out and leave zero. We knew already, thanks to the Wiener–Khinchin theorem and Eq. (27), that the imaginary parts will sum to zero. Now we can establish this for the real parts as well.

The left-most column of ${\bf G}$ is also its top row, which tells us the averages of each column of ${\bf F}$. So, if the Zauner orbits leave only one element in a column unspecified, then we can fill in that element, because we know the average over the whole column vector. This proves the $3d$ conjecture in dimension $d=9$.

We now turn to the simplification that the Zauner conjecture provides for numerical search efforts. By postulating that the SIC fiducial we are looking for is a Zauner eigenvector, we can significantly reduce the effective size of the search space. First, suppose that $U$ is a unitary of order $n$, so that

[TABLE]

and the eigenvalues of $U$ can all be written

[TABLE]

with $m$ an integer. The projector onto the eigenspace with this eigenvalue is

[TABLE]

We can restrict our numerical search to this subspace by projecting our vectors into it,

[TABLE]

at each iteration of the optimization algorithm.

Most of the known solutions were found by postulating Zauner symmetry. Scott has also found several solutions by assuming that the fiducial was an eigenvector of another Clifford unitary. For an in-depth exposition of these variations, see Scott:2017 .

VI Exhaustive Searches

Suppose that ${\left|\psi_{0}\right\rangle}\in\mathcal{H}_{d}$ is a vector in a Weyl–Heisenberg SIC, and let $U$ be a Clifford unitary. Applying $U$ to the vector ${\left|\psi_{0}\right\rangle}$ will yield some vector,

[TABLE]

The Weyl–Heisenberg orbit of ${\left|\psi_{0}\right\rangle}$ is a SIC, so what about the orbit of ${\left|\chi_{0}\right\rangle}$ under the same group? We define

[TABLE]

and we consider the squared magnitudes of the inner products

[TABLE]

Because $U$ is a Clifford unitary, it maps $D_{kl}$ to some Weyl–Heisenberg operator $D_{k^{\prime}l^{\prime}}$, with any phase factor dropping out when we take the magnitude of the inner product. So,

[TABLE]

meaning that the image of our original SIC under the mapping $U$ is also a SIC. One way in which the Hesse SIC is remarkable is that it is invariant under the entire Clifford group. For contrast, we can take the vector

[TABLE]

which differs from the Hesse SIC fiducial in Eq. (13) by a sign. The orbit of this vector under the Clifford group is a set of four separate SICs, comprising 36 vectors in all—the so-called Norrell states, which are significant in the theory of quantum computation Veitch:2014 ; Stacey:2016c .

We consider two SICs equivalent if they can be mapped into each other by a Clifford unitary. In fact, it is convenient to extend the Clifford group by including the anti-unitary operation of complex conjugation. The extended Clifford group for dimension $d$, $\hbox{EC}(d)$, is the set of all unitary and anti-unitary operators that send the Weyl–Heisenberg group to itself. For (extensive) details, we again refer to Appleby Appleby:2005b ; Appleby:2009 .

In order to search the space as exhaustively as possible and create a catalogue of all essentially unique SICs, Scott’s code chooses initial vectors at random under the unitarily invariant Haar measure on the complex projective space $\mathbb{C}P^{d-1}$. Once enough solutions are found—generally, this means hundreds of them—the code then refines their precision, as described above. Then, we must identify unique orbits under the extended Clifford group. This last step is computationally demanding, because we must translate each solution vector ${\left|\psi\right\rangle}$ by each element in the extended Clifford group $\hbox{EC}(d)$. However, in the process, Scott’s algorithm also finds the stabilizer group of each fiducial, which is important information. The task of determining when two SICs are equivalent up to a unitary or anti-unitary transformation has been discussed in depth by Zhu Zhu:2012 , and we expect that additional theoretical insights may lead to an improved algorithm for this step.

Following this procedure, Scott has carried out exhaustive searches in dimensions up to $d=90$. We strongly expect his catalogue of solutions to be complete up to that point: All Weyl–Heisenberg SICs in those dimensions are equivalent to the ones tabulated, up to equivalence under the extended Clifford group.

VII Discussion

In the preceding sections, we have described the process of finding SIC fiducial vectors numerically. However, some patterns among SICs have only become apparent when exact solutions were studied carefully. Suppose we refrain from taking the magnitude-squared in our definition of a SIC, Eq. (1). Then

[TABLE]

for some set of phases $\{e^{i\theta_{jk}}\}$. (In fact, one can reconstruct the SIC from knowing the phases Appleby:2011a .) It was recently discovered that when $d>3$, for all the known Weyl–Heisenberg SICs, these phases have a remarkable meaning in algebraic number theory: They are units in ray class fields and extensions thereof Appleby:2016 . This is a topic to which we can hardly do justice here, and indeed, treatments accessible to anyone who is not already an algebraic number theorist have only recently been attempted Bengtsson:2016 ; FOOP:2017 . For now, we content ourselves with the observation that this area of number theory is the territory of Hilbert’s twelfth problem, one of the still outstanding questions on history’s most influential list of mathematical challenges Gray:2001 . (Specialists may recall that according to the Kronecker–Weber theorem, any abelian extension of the rationals is contained in a cyclotomic field. When we instead consider abelian extensions of real quadratic fields, the analogue of the cyclotomic fields are the ray class fields. The phases of Weyl–Heisenberg SICs appear to be playing a role regarding ray class fields much like the role that roots of unity play with cyclotomic fields. Moreover, recalling Eq. (14), it is intriguing that in the real-vector-space version of equiangular lines, we discard a phase factor that is a unit among the ordinary integers, while in the complex Weyl–Heisenberg case, the phases turn out to be units among algebraic integers.) From Hilbert space to Hilbert’s twelfth problem! What physicist would ever have anticipated that? And who could turn down the opportunity to intermingle two subjects that had seemed so widely separated?

SICs have found relevance, not just in quantum computation, but in signal-processing tasks like high-precision radar Howard:2006 and speech recognition Balan:2009 . In February 2016, our colleague Marcus Appleby attended a conference in Bonn, Germany on uses of the Weyl–Heisenberg group. Many participants were engineers, including representatives from the automotive and cell-phone industries. Appleby was told that if he managed to construct a SIC in dimension 2048, he should patent it Marcus-email . At the moment, dimension 2048 is beyond our abilities for algebraic or numerical solutions, but this may not always be the case.

VIII Acknowledgments

We are deeply indebted to Andrew J. Scott, the guru of SIC numerical solutions, for code and for discussions. We also thank Marcus Appleby for many conversations, and Gary McConnell for email feedback on the original arXiv version of this article. This research was supported in part by MCH’s Oracle Undergraduate Research Fellowship at UMass Boston.

IX Author Contributions

MCH ran the calculations on Chimera to find SICs in dimensions 122 through 151. BCS wrote the paper. CAF directed the research, contributed to the bibliography and worked with BCS in revising the paper.