Quantum measurements with prescribed symmetry

Wojciech Bruzda Institute of Physics, Jagiellonian University, Krak´ow,Poland

Dardo Goyeneche Institute of Physics, Jagiellonian University, Krak´ow,Poland and Faculty of Applied Physics and Mathematics, Technical University of Gda´nsk,80-233 Gda´nsk,Poland

Karol Zyczkowski˙ Institute of Physics, Jagiellonian University, Krak´ow,Poland and Center for Theoretical Physics, Polish Academy of Sciences, Warsaw, Poland (Dated: July 4, 2017) We introduce a method to determine whether a given generalised quantum measurement is iso- lated or it belongs to a family of measurements having the same prescribed symmetry. The technique proposed reduces to solving a linear system of equations in some relevant cases. As consequence, we provide a simple derivation of the maximal family of Symmetric Informationally Complete measure- ment (SIC)-POVM in dimension 3. Furthermore, we show that the following remarkable geometrical structures are isolated, so that free parameters cannot be introduced: (a) maximal sets of mutually unbiased bases in prime power dimensions from 4 to 16, (b) SIC-POVM in dimensions from 4 to 16 and (c) contextuality Kochen-Specker sets in dimension 3, 4 and 6, composed of 13, 18 and 21 vectors, respectively.

Keywords: Mutually unbiased bases, SIC- for different convenient purposes. For example, POVM, defect of a unitary matrix. from the one-parametric family of SIC-POVM ex- isting in dimension three [1] only a single member maximizes the informational power, that is, the classical capacity of a quantum-classical channel I. INTRODUCTION generated by the SIC-POVM [14]. Furthermore, inequivalent sets of MUB provide different estima- Positive Operator Valued Measure (POVM) is tion of errors in [15]. the most general kind of measurement in quantum Highly symmetric quantum measurements, like mechanics, which generalizes projective measure- SIC-POVM or MUBs play an important role in ments. Some POVM having a prescribed symme- and foundations of quan- try play a crucial role in quantum mechanics: Sym- tum theory. On one hand, it is interesting itself metric Informationally Complete SIC-POVM [1, 2] to design a mathematical tool that allows one to and Mutually Unbiased Bases (MUB) [3, 4]. These construct a family of POVM having a prescribed geometrical structures have important applications symmetry from a given particular solution. On in quantum theory: SIC-POVM and MUB allow us the other hand, construction of such families of to unambiguously reconstruct any density matrix solutions provides flexibility when designing ex- of size d [1, 3] and define entropic uncertainty re- perimental implementations of these measurement arXiv:1704.04609v3 [quant-ph] 5 Jul 2017 lations [5, 6]. Even more, MUB are important to sets. For instance, a detailed description of a com- detect entanglement [7], bound entanglement [8], plete list of solutions of a set of k MUB in di- and to lock classical information in quantum states mension d may be helpful in tackling the problem [9]. whether an extended set of k +1 MUB exists. Fur- Finitely many SIC-POVM are known in dimen- thermore, it is also interesting to highlight those sion d ≤ 64 [10], including a 1-parametric family in quantum measurements having a prescribed sym- dimension three [1]. The existence of SIC-POVM metry that do not belong to a family, which makes with free parameters in dimension d > 3 is still an them special. A possible application of such iso- open problem. Moreover, maximal sets of MUB lated cases is the existence of a unique solution are known to be isolated in dimensions two to five optimizing a given function. For instance, isolated [11] and finitely many non-equivalent maximal sets solutions might define sets of measurement having of MUB are known for dimensions three to five [12] a unique maximal violation of a Bell inequality, and N qudit systems [13]. On the other hand, fam- which is a fundamental ingredient for self-testing ilies of symmetric POVM are useful for practical [16]. applications, as the parameters can be optimized In this work, we present a method to intro- 2 duce free parameters in generalized measurements jectors having the prescribed symmetry (1) given having a predefined geometrical structure. The by a real symmetric matrix S. For example, the N−d method proposed divides the entire non-linear case S = IN + d(N−1) (JN − IN ) corresponds to (1) (2) problem, called PNL, into a linear problem PL equiangular tight frames composed by N vectors (3) in dimension d. Here, N denotes the matrix of and a secondary non-linear problem PNL, which is J (1) size N having all entries equal to unity. We recall simpler than PNL. Remarkably, in some cases the (2) that a set of N vectors {|φii} defined in dimen- linear problem PL provides a definite answer to sion d forms an equiangular tight frame (ETF) if the full problem P . 2 NL |hφi|φji| = d(N − 1)/(N − d), for every i 6= j = The paper is organized as follows: in Section 0, . . . , d2 − 1. II we establish a connection between any given A remarkably important subclass of ETF is POVM and certain hermitian unitary matrices given by the so-called Symmetric Informationally having constant diagonal. Complete (SIC)-POVM [1], corresponding to the In Section III we apply the notion of a defect of 2 case N = d . Also, two orthonormal bases |φii a unitary matrix to identify isolated cases of gener- and |ψji in dimension d define a pair of mutually alized quantum measurements having a prescribed 2 unbiased bases (MUB) if |hφi|ψji| = 1/d, for ev- symmetry, for which no free parameters can be in- ery i, j = 0, . . . , d − 1. A set of m orthonormal troduced. Furthermore, in other cases we present bases is mutually unbiased if every pair of the set a constructive method to extend known solutions is mutually unbiased. Also, a set of m MUB in to an entire family by introducing free parameters. dimension d has associated the symmetric matrix In Section IV we show that known maximal sets 1 S = Idm + d (Jm − Im) ⊗ Jd. For a recent review of MUB in dimensions 4, 8, 9 and 16 and known on discrete structures in Hilbert spaces, including SIC-POVM in dimensions 4 to 16 are isolated. We MUB and SIC-POVM, see Ref. [17]. also study the robustness of our results for a given Let us recall a close connection existing between accuracy in specifying the POVM, which allows POVM and tight frames. A set of rank-one pro- us to derive conclusive results from approximate jectors {Πj} defines a tight frame if there exists solutions. In Section V we find an upper bound PN−1 a real number A > 0 such that j=0 Tr(ΩΠj) = for the maximal number of free parameters that 2 can be introduced in sets of 2 ≤ m ≤ d + 1 MUB ATr(Ω ) = A, for any rank-one projector Ω act- in dimension d and in some classes of equiangular ing on dimension d. Therefore, POVM are tight tight frames. Moreover, we show how the method frames for A = c. A crucial property for our work works to give the known one-parameter families of is the fact that the Gram matrix associated to a SIC-POVM in dimension three. In Section VI we tight frame, or POVM, is closely related to an her- prove that some existing Kochen-Specker sets from mitian unitary matrix, as we will see in Proposition quantum contextuality are isolated. In Section VII 1. We recall that the Gram matrix of a set of N we summarize our results and pose open questions. vectors |φji is given by Gij = hφi|φji, (2)

II. QUANTUM MEASUREMENTS AND where i, j = 0,...N − 1. For example, the Gram TIGHT FRAMES matrix of an equiangular tight frame composed by N vectors in dimension d has the form

 iα12 iα1N  A POVM {Γj} is a set of N positive semidefi- 1 re . . . re nite subnormalized operators defined in dimension  re−iα12 1 . . . reiα2N  N−1   d such that P Γ = , where is the iden- GETF = . . . . , j=0 j Id Id  . . .. .  tity matrix of size d. Along the work, we will re-   re−iα1N re−iα2N − 1 ... 1 strict our attention to rank-one POVM and con- (3) sider rank-one projectors Π , being proportional to j where r2 = d(N − 1)/(N − d). Furthermore, the elements of POVM. That is, Π = cΓ , where j j the Gram matrix of a set of m + 1 MUB c = N/d. For simplicity, we will refer to the set of {Id,H1,H2,...,Hm} in dimension d is given by projectors {Πj} as a POVM, understanding that   they are formally proportional to the elements of Id H1 H2 ··· Hm a POVM. The rank-one projectors Π satisfy the † † † j  H1 Id H1 H2 ··· H1 Hm   † † †  geometrical relation  H H H ··· H H  GMUB =  2 2 1 Id 2 m  , Tr(Π Π ) = S , (1)  . . . . .  i j ij  ......  † † † where S is a real symmetric matrix of size N. It Hm HmH1 HmH2 ··· Id is interesting to ask about the most general pro- (4) 3 where H1,H2,...,Hm are suitable unitary com- Eq.(1). That is, U is an hermitian unitary ma- † plex Hadamard matrices, so that HiHi = I and trix having constant diagonal Uii = 1 − 2d/N and 2 2d p |(Hi)jk| = 1/d for every j, k = 0,...,N − 1 and |Uij| = δij − N Sij. Note that we have introduced i = 1, . . . , m, see [19]. Let us establish a connection a parameter t in Eq.(6) for convenience, which can between Gram matrices of POVM and a special be set to t = 1 after applying our method. The kind of unitary hermitian matrices. full problem is given as follows: (1) PROPOSITION 1. Let Πj be a rank-one POVM Problem PNL : Find the most general matrix composed by N vectors in dimension d and G be Vτ (t) of size N of the form (6), initially depending the corresponding Gram matrix. Then, the matrix on τ parameters, such that 2d U = N − G is unitary. I N † Vτ (t)Vτ (t) = IN . (7) Proof. A Gram matrix G represents a POVM com- posed by N vectors in dimension d ≤ N if and only This problem implies to solve a system of non- 2 N if G = d G (cf. Prop. 1 in Ref.[18]). From this linear coupled equations, which depends on τ = property and taking into account that Tr(G) = N, [N(N − 1)/2 − z] − (N − 1) non-trivial variables, the spectrum of G satisfies where z is the number of zeros existing in the strictly upper triangular part of the matrix R. λ(G) = (N/d, . . . , N/d, 0,..., 0). (5) Note that τ is composed by the total number | {z } | {z } d N−d of parameters Rij, i.e. N(N − 1)/2 − z, minus the number of trivial variables (N − 1). These 2d Therefore, U = IN − N G is a unitary matrix. trivial parameters can be absorbed by applying the transformation V → EVE†, where E = From Proposition (1) we realize that the exis- Diag(1, eitR01 , . . . , eitR0(N−1) ). In order to simplify tence of a POVM with a symmetry prescribed by the resolution of problem PNL we define the fol- the matrix S from (1) is equivalent to prove the lowing linear problem: existence of a unitary hermitian matrix U having (2) positive constant diagonal Uii = 1−2d/N and sat- Problem PL : Find the most general matrix 2d p isfying |Uij| = N Sij for i 6= j. Unitary matrices Vτ (t) of size N, initially depending on τ param- 2d eters, such that U = IN − N G have been recently studied for the particular case of equiangular tight frames [20]. In the Bloch sphere associated to a one- system d † lim [Vτ (t)Vτ (t) ] = 0. (8) we have some relevant geometrical structures: or- t→0 dt thonormal basis (line), 3 MUB (3 orthogonal lines) Using (6), we can explicitly write Eq.(8) as and SIC-POVM (tetrahedron). All these struc- tures are unique, up to a global rotation. In higher X −2V V R + V V (R −R ) = 0, (9) dimensions, however, some geometrical structures k,k k,j j,k k,l l,j k,l j,l allow one to introduce free parameters, that can- l6=j,k not be absorbed in a global rotation. In Section for 1 j < k N and 1 ≤ l ≤ N, which is a linear III we introduce the method, which considerably 6 6 problem on variables R . Note that P(2) ⊂ P(1) , simplifies the study of this problem. ij L NL as Eq.(8) is a necessary condition to obtain Eq.(7). (2) The linear problem PL allows us to simplify III. RESTRICTED DEFECT AND FREE the full problem PNL by determining r out of τ PARAMETERS variables Rij, where r is the number of linearly (1) independent equations (9). After solving PL , the In this section, we derive the method to intro- remaining number of free parameters Rij lead us duce the maximal possible number of free param- to the definition of the restricted defect ∆ of the eters into a given POVM composed by N vectors hermitian unitary matrix U. It reads, in dimension d associated to a given symmetric matrix S, see Eq.(1). Using Proposition 1, this ∆ = τ − r, (10) problem is equivalent to finding the most general real antisymmetric matrix R of size N such that where τ = (N − 1)(N − 2)/2 − z and z is the num- ber of zeros existing in the strictly upper triangu- itRij Vij(t) = Uije , (6) lar part of the matrix U. Note that this quantity can be considered as a defect of a unitary matrix 2d is a unitary matrix, provided that U = IN − N G is [21], adopted to the case of matrices with a special associated to a given particular POVM satisfying structure. The standard defect was used to define 4 an upper bound on the number of free parameters E we find that allowed by complex Hadamard matrices [19] and  0 0 1 1  forms, by construction, an upper bound for the re- 1  0 0 1 −eitR24  stricted defect. In both cases, the defect equal to Vτ (t) =   . (14) zero implies that a given solution is isolated, so no 2  1 1 0 0  −itR24 free parameters can be introduced. 1 −e 0 0 In general, the restricted defect represents an up- Problem P(2) implies the following equation per bound for the maximal number of free param- L eters allowed by the full problem P . If ∆ = 0, NL R24 = 0, (15) then the full problem PNL is solved by the linear (1) where r = 1 and, therefore, ∆ = 0. Thus, we can- problem PL . In this case, it is not possible to in- troduce free parameters into the matrix V . On the not introduce free parameters in Eq.(14), which other hand, if ∆ > 0, it is necessary to solve an ad- implies that the considered pair of MUB is iso- ditional non-linear problem in order to determine lated. Indeed,√ |φii = |ii, i = 0, 1 and |ψ±i = the continuous family of solutions. (|0i ± |1i)/ 2 is the unique pair of MUB existing in dimension 2, up to a global rotation [11]. (3) Problem PNL : Find the most general matrix V∆(t) of size N, initially depending on ∆ parame- ters, such that IV. ISOLATED MUB AND SIC-POVM

† V∆(t)V∆(t) = IN . (11) In this section, we study the problem of intro- ducing free parameters in MUB and SIC-POVM. Note that Problem P(2) is simpler than Problem Our first main result consists in proving that max- NL imal sets of MUB existing in low prime power di- PNL as ∆ < τ. This is so because r > 0 in Eq.(10). (2) mensions. After solving Problem PNL we can assume that t = 1, without loss of generality. In Section IV we PROPOSITION 2. Maximal sets of d+1 MUBs will apply our results to SIC-POVM and maximal in dimensions d = 4, 8, 9 and 16 are isolated. sets of MUB. The results arises as follows. The upper triangu- Let us first illustrate the method in action by lar part of the Gram matrix associated to a set of m MUB in dimension d contains z = md(d − 1)/2 considering two MUB for a single qubit√ system: zero entries. Given that the size of the Gram ma- |φii = |ii, i = 0, 1 and |ψ±i = (|0i ± |1i)/ 2. The Gram matrix (4) associated to this set of m = 2 trix of the set is NMUB = md and m = d + 1 we 2 MUB is given by have z = d(d − 1)/2. Therefore, the matrix G contains τ = (NMUB − 1)(NMUB − 2)/2 − z pa-  2 0 1 1  rameters. The next step consists in determining 1 0 2 1 −1 how many of these parameters are remaining af- GMUB =   . (12) (1) 2  1 1 2 0  ter solving the linear problem PL . To this end, 1 −1 0 2 we calculated the number of linearly independent equations of the linear system defined in Eq.(9) A family of two MUB stemming from this fixed set for the cases d = 4, 8, 9 and 16, finding r4 = 141, would have associated a Gram matrix of the form r8 = 2233 and r9 = 3556 and r16 = 34545, respec- tively. By using these results and Eq.(10) we find  2 0 eitR13 eitR14  that the restricted defect ∆ vanished in all these 1  0 2 eitR23 −eitR24  cases. We based our calculation of the restricted GMUB =   , 2  e−itR13 e−itR23 2 0  defect ∆ on the maximal sets of MUB provided in e−itR14 −e−itR24 0 2 Refs. [4, 12]. In dimension d = 9, we considered (13) Ref. [22] to obtain simpler expressions of results and, from Prop.(1), the unitary matrix U = presented in [4]. 2d (1) IN − N GMUB. Note that full problem PNL Let us now consider the case of SIC-POVM. It initially depends on τ = 1 non-trivial param- is well-known that SIC-POVM for a qubit system eter, as R13, R14 and R23 can be absorbed is essentially unique, as it represents the regular by considering the diagonal unitary operator tetrahedron inscribed into the Bloch sphere, up E = diag[1, e−iR23 , eiR13 , eiR14 ] and the redefini- to a global rotation. Furthermore, single parame- tion V → EVE†. Therefore, according to Eq.(6), ter families of SIC-POVM for a system ex- and after considering the diagonal transformation ists [1], which represent the most general solution 5

[23]. For higher dimensions, the problem of in- for the restricted defect (10) we need to consider troducing free parameters in SIC-POVM is still the following two bounds: (i) upper bound for the open. As a preliminary result, exhaustive numeri- maximal perturbation of σ0 and (ii) lower bound cal simulations indicate that free parameters can- for the maximal perturbation of σ1. In Appendix not be introduced in SIC-POVM, at least in low A we show that dimensions higher than three. For this problem, (1) 0 |σ − σi| ≤ f(d, N) s, (18) we have solved the linear problem PL in dimen- i sions d = 4,..., 16, obtaining the following results. for i = 0, 1, where s is the inaccuracy quantificator PROPOSITION 3. SIC-POVM in dimensions defined in Eq.(16) and d = 4,..., 16 are isolated. s 26d2  2d2 N − d This result also includes the Hoggar lines [24], f(d, N) = 1 − (19) a special class of SIC-POVM defined in dimension N N N(N − 1) d = 8. We considered the total number of param- eters τ = (NSIC − 1)(NSIC − 2)/2, as there is no for N > 2d. Let us now find the smallest pos- pair of orthogonal vectors in SIC-POVM (z = 0), sible value of s such that the critical condition 2 0 0 and NSIC = d . In order to prove Proposition 3 σ0 = σ1 holds, which imposes an upper bound (1) for the confidence region of the restricted defect we solved the linear problem PL for both ana- lytical [24–28] and highly accurate numerical SIC- ∆. By considering Eq.(18) we find that ∆ does POVM [10]. In all the cases we have found ∆ = 0, not change its value for 0 ≤ s ≤ smax, where −1 which implies that free parameters cannot be intro- smax = σ1(2f(d, N)) . Note that σ1 depends on duced. Calculations of the restricted defect have the exact solution, which is not known if the exact been done in Matlab. solution is not available. By using Eq.(18) we find 0 Let us now study the robustness of the restricted that σ1 ≤ σ1 + f(d, N) s, which implies defect ∆ under the presence of inaccuracies in de- 0 scribing the POVM. Given the set of vectors {φ } σ1 + f(d, N) s j smax ≤ . (20) defined in Eq.(2) we quantify the inaccuracy in ap- 2f(d, N) 0 proximate vectors {φj} ≈ {φj} by introducing the inaccuracy factor: Note that this inequality provides a confidence re- gion only if f(d, N) s  1. 1 0 Confidence regions for SIC-POVM and maximal s = √ max kφj − φjk. (16) d j sets of MUB for 2 and 3 are depicted in Fig. 1 and Fig. 2, respectively. For the case of The factor sµ quantifies the maximal allowed in- 4 qubit systems we have solutions with precision accuracy in entries of vectors φj. For example, in s = 10−32, whereas the upper bound for the con- the case of approximate solutions {φ0 } having k −3 j fidence region is s4 ≈ 4 × 10 . For MUB, we digits of precision we have s ≈ 10−k. In our study calculated the restricted defect ∆ by considering of robustness of the defect ∆, we simulate the in- analytic solutions in all the cases [4, 12, 22]. troduction of inaccuracies by considering

0 (φj)i = (φj)i + sξi, (17) V. FREE PARAMETERS IN MUB AND SIC-POVM where (x)i denotes the i-th entry of the vector x and ξi are random numbers uniformly distributed in the interval [−1, 1]. In Section IV we have proven that some maxi- Let us assume that R is the real matrix associ- mal sets of MUB and SIC-POVM are isolated in ated to the linear system of equations (9). Note low dimensions. In this section, we first calculate that the number of linearly independent equations the restricted defect for sets of m = 2, . . . , d + 1 of such system is given by rank(R) = r. When con- MUB in dimensions d = 2,..., 8 (see Table I). In- sidering inaccuracies in the POVM the rank of the terestingly, the restricted defect for a pair of MUB 0 B1 and B2 coincides with the standard defect of perturbed matrix R and R may differ. Therefore, † we need to study how much the singular values of the complex Hadamard matrix H = B1B2. In- R0 are affected under the presence of inaccuracies. deed, the maximal number of free parameters that In particular, we are interested on the perturba- can be introduce in the Gram matrix tion of the two smallest singular values σ0 = 0 and ! H σ1 > 0, which are the responsible for the variation Id GMUB = † , (21) of the rank. In order to obtain a confidence region H Id 6

m \ d 2 3 4 5 6 7 8 9 2 0 0 3 0 4 0 21 16 3 0 0 3 0  0 27 20 4 - 0 0 0 ? 0 19 32 5 - - 0 0 ? 0 7 0 6 - - - 0 ? 0 0 0 7 - - - - ? 0 0 0 8 - - - - - 0 0 0 9 ------0 0 10 ------0

FIG. 1: Smallest singular values σ0 and σ1 of R as TABLE I: Upper bound on the maximal number of a function of the inaccuracy factor s for SIC-POVM. free parameters ∆ allowed by subsets of m MUB in Contour and filled symbols represent σ0 and σ1, re- dimension d. The results do not depend on selecting m subsets of MUB out of the full set of d + 1 MUB. spectively, for 2 (◦) and 3 () qubit systems. Each case is averaged over 8 samples randomly chosen and As a remarkable observation, maximal sets of MUB generated from the approximate SIC-POVM provided are isolated. Also, subsets of m ≥ 6 in dimension 8 in Ref. [10], which has accuracy s = 10−32. The con- and m ≥ 5 MUB in dimension 9 are isolated in all the fidence regions (blue and red rectangles) are given by cases. Question marks denote our lack of knowledge values of s existing between zero and the value de- about given number of MUB, while  indicates the case termined by the intersection of the lower and upper of three MUB in dimension six which is still considered bounds. Here, s2 and s3 stand for 2 and 3-qubits sys- as unresolved. tems, respectively. Outside the confidence regions it is not possible to discriminate between singular values σ0 and σ1. a pair of MUB in dimension 6 is four, coinciding with generic defect of complex Hadamard matri- ces of size 6 [21]. However, note that there is an exceptional pair of MUB for which ∆ = 0, as an isolated complex Hadamard matrix of size 6 exists [29]. Generic defect for triplets of MUB is not well understood (), see Table I and Ref. [30]. In dimensions 8 and 9 we restricted our attention to subsets of MUB arising from the maximal sets defined in Refs. [4, 12]. Note that subsets of m ≥ 6 MUB are isolated in dimension 8, whereas several families of m = 5 MUB exist [31]. Another obser- vation is that the restricted defect ∆ for maximal sets of d+1 MUB in dimension d coincides with the defect for d MUB. This is so because the (d+1)-th FIG. 2: Smallest singular values σ0 and σ1 of R as a MUB is univocally determined by the first d MUB. function of the inaccuracy factor s for maximal sets of For subsets of m < d, the restricted defect for m MUB. Contour and filled symbols represent σ and σ , 0 1 MUB may depend on the subset chosen. However, respectively, for 2 (◦) and 3 () qubit systems. Each case is averaged over 8 samples randomly chosen and results presented in Table I are consistent for every generated from analytic solutions. subset of MUB. We also studied the restricted defect for equian- gular tight frames composed by N = k2 vectors coincides with the maximal number of parameters in dimension d = k(k − 1)/2, typically denoted as that can be introduced in H. Let us explain some ETF(d, N) [32]. These ETF have associated the details concerning Table I. First, note that ∆ = 0 following hermitian unitary matrices [20]: for every subset of 2 ≤ m ≤ d + 1 MUB in prime i1j2−j1i2 dimensions d = 2, 3, 5, 7. This is so because ev- Ui1+ki2+1,j1+kj2+1 = ω , (22) ery complex Hadamard matrix involved in the set 2πi/k is equivalent to the Fourier matrix, which is iso- where i1, i2, j1, j2 ∈ {0, .., k − 1} and ω = e . lated in prime dimensions [19]. For triplets in di- Matrix (22) is equivalent to the tensor product of st mension 4, we have ∆ = 3, which coincides with Fourier matrices, Fk ⊗Fk, where (Fk)st = ω . Ta- the maximal number of free parameters that can ble II summarizes the restricted defect for matrix be introduced [11]. Generic restricted defect for U in low dimensions. Results are shown in Table 7

II. For prime values of k, the formula Therefore, we generate six 1-parametric families of SIC-POVM in dimension three: 1 ∆ = (k + 1)(k − 1)(k − 2), (23) 2 S1 : V∆(R23 = 0,R48), S2 : V∆(R23 = π, R48), S : V (R = 0,R ), S : V (R = π, R ), matches all solutions presented in Table II, so we 3 ∆ 26 48 4 ∆ 26 48 are tempted to believe that it holds for any prime S5 : V∆(R89 = 0,R48), S6 : V∆(R23 = π, R48)(25). k. Here we note that ∆ = 4 and six 1-parametric real solutions exist. These six solutions belong to k ∆ k ∆ k ∆ the 4-dimensional tangent plane defined by Eq.(8) 2 0 7 120 12 1237 and do not fit into a lower dimensional tangent 3 4 8 273 13 924 space, which explains why ∆ cannot take a lower 4 21 9 352 14 1632 value. Furthermore, solutions (25) are equivalent, 5 36 10 576 in the sense that we can transform one into the 6 112 11 540 other by applying permutation of rows or columns and multiplication of diagonal unitary operations TABLE II: Maximal number of free parameters that to the Gram matrix, which is equivalent to relabel can be introduced in ETF composed by N = k2 vec- and apply global phases to vectors. Solution (25) tors in dimension d = k(k − 1)/2. The case k = 2 represents the most general SIC-POVM existing corresponds to an ETF(3,4) which is isolated (regular in dimension three [23], up to equivalence. We simplex in dimension 3). Also, k = 3 has associated a remark that the generic hermitian defect for a SIC- SIC-POVM in dimension 3, where ∆ = 4 but only 2- POVM in dimension 3 is ∆ = 2, with the only parametric families exist [1]. For the case k = 4 there exception of the particular vector |φ i, where ∆ = exists several 6-parametric families of ETF(6,16) [20]. 00 4, however from this fact we cannot define a larger family. Let us now study the SIC-POVM problem in di- mension d = 3. By√ considering the fiducial state VI. ISOLATED KOCHEN-SPECKER SETS |φ00i = (1, −1, 0)/ 2 a SIC-POVM is given by [25]: In this section, we apply our method presented |φ i = XsZt|φ i, (24) in Section III to show that some sets of vectors st 00 used in a proof of the Kochen-Specker contextually where X|ji = |[j ⊕ 1]i, Z|ji = ωj|ji, ω = e2πi/3 theorem [33], typically called KS sets, are isolated. and ⊕ means addition modulo 3. The 9 × 9 Gram KS sets are collections of N vectors in dimension d, which contain m subsets of d vectors forming or- matrix GSIC of the SIC-POVM and its associated 2d thonormal basis. Some of these orthonormal bases unitary matrix U = IN − N GSIC depends on τ = 36−8 = 28 parameters, where the N −1 = 8 trivial have common vectors, so that N < md. These intersections are crucial to prove that a determin- parameters R1,j for j = 2,..., 9 have been set as zero. The linear system of equations (8) associate istic local hidden variable theory is not possible [33]. That is, for a system prepared in a quantum to problem P(1) has r = 24 linearly independent L state ρ and a set of KS vectors {φ , . . . , φ } equations, which provides a 4-dimensional complex 0 N−1 it is not possible to end up with N determinis- set of solutions R , depending on four parameters: ij tic probabilities P = Tr(ρ|φ ihφ |) ∈ {0, 1}, for R ,R ,R , and R . The additional restriction k k k 23 26 48 89 k = 0,...,N − 1. Therefore, the assumption of to have real parameters imply hidden determinism in quantum mechanics implies that predefined values of observables depend on the R − 3R = R − 3R = R − R = 0, 23 89 23 26 89 26 context in which measurements were implemented. The original proof given by Kochen and Specker in- which is equivalent to R = R = R /3. Af- 26 89 23 volves N = 117 vectors in dimension d = 3 [33]. ter setting t = 1, we obtain three solutions to (1) Subsequently, examples exhibiting a lower number problem PL : V∆(R23,R48), V∆(R26,R48) and of vectors were found. Some remarkable examples V∆(R89,R48). Now, we are in position to solve are KS sets composed by N = 13 vectors in dimen- (2) the non-linear problem PNL, which is much sim- sion d = 3 (Yo and Oh [34]), N = 18 vectors in pler than the full non-linear problem PNL. Indeed, dimension d = 4 (Cabello et al. [35]) and N = 21 (2) PNL implies to solve trivial trigonometric equa- vectors in dimension d = 6 (Lisonek et al. [36]). tions, which give us the solutions R26 ∈ {0, π}, Let us now apply our method to prove that these R89 ∈ {0, π} and R23 ∈ {0, π}, respectively. three inequivalent KS sets are isolated. The first 8 important observation is that the three KS sets maximal sets of MUB in dimension 4, 8, 9 and 16 form three POVM. This means that Proposition 1 and known SIC-POVM in dimensions 4−16 are iso- holds for these sets and, therefore, the method to lated. In particular, a special class of SIC-POVM introduce free parameters presented in Section III existing for 3-qubit systems, called Hoggar lines, can be applied. In order to do so we have to cal- is isolated (see Section IV). Moreover, we calcu- culate the restricted defect ∆ defined in Eq. (10), lated an upper bound for the maximal number of which is a function of the total number of param- free parameters that can be introduced in subsets eters τ and the number of linearly independent of 2 ≤ m ≤ d + 1 MUB in dimensions d = 2 − 9 (2) (see Section V). The same study has been done for equations associated to Problem PL (see Section III). The geometrical structure is determined by equiangular tight frames in low dimensions, which the orthogonality restrictions imposed by the KS define equiangular POVM (see Table II). sets. For the above mentioned three KS sets we As a further result, we studied the robust- have shown that they are isolated. The way to ness of our method under the presence of in- proceed is similar to the proof that maximal sets accuracies in defining the generalised measure- of MUB or SIC-POVM are isolated (see Proposi- ment, which allowed us to establish a confidence tion 2). However, there is a minor additional re- region for the maximal possible number of free mark: the sets are isolated despite the restricted parameters that can be introduced (see Section defect ∆ of the sets is non-zero. This is so be- IV). The importance of robustness relies on the cause the apparently remaining ∆ free parame- fact that some geometrical structures, like SIC- ters can be absorbed by considering a sequence of POVM, are established analytically in low dimen- non-trivial emphasing in the Gram matrices, which sion only, whereas accurate numerical solutions means that the free parameters can be absorbed as exist in every dimension d ≤ 121 and also in global phases of the KS vectors. Table III resumes d = 124, 143, 147, 168, 172, 195, 199, 228, 259 and the details of our calculations. 323 [10, 37]. A Matlab source code to support cal- culation of the restricted defect and additional fea- N d z τ r ∆ # free parameters tures is available on the GitHub platform: https: 13 3 24 78 66 12 0 //github.com/matrix-toolbox/defect. 18 4 63 90 83 7 0 Additionally, have proven that three Kochen- 21 6 105 105 103 2 0 Specker contextuality sets are isolated (see Section VI). Namely, 13 vectors in dimension 3 [34], 18 TABLE III: Isolated KS sets composed by 13 vectors in vectors in dimension 4 [35] and 21 vectors in di- dimension 3, 18 vectors in dimension 4 and 21 vectors mension 6 [36]. in dimension 6. The number of zeros (z) appearing Finally, we pose some intriguing open questions: into the upper triangular part of the Gram matrix of (i) Are maximal sets of MUB isolated in every KS set, total number of parameters (τ), rank of the prime power dimension? (ii) Are SIC-POVM iso- linear system defined in Eq. 8 (r), and restricted defect lated in every dimension d > 3? Furthermore, it (∆ = τ − r) are defined in Section III. For these three would be welcome to develop a more efficient soft- KS sets the free parameters produced by a positive (1) restricted defect ∆ can be absorbed as global phases ware to solve the linear problem PL for POVM of the vectors. having N > 300 elements, e.g. maximal sets of MUB or SIC-POVM in dimension d > 16.

VII. CONCLUSION Acknowledgements

We studied the problem to introduce free pa- We are grateful to Markus Grassl for his ad- rameters in a given POVM having prescribed sym- vice to study the robustness of the hermitian de- metry, where mutually unbiased bases (MUB) and fect under the presence of inaccuracies, to Woj- symmetric informationally complete (SIC)-POVM ciech Tadej for discussions and fruitful correspon- are relevant examples (see Section II). In particu- dence concerning the restricted defect of a uni- lar, our method allows us to determine whether a tary matrix and to Ingemar Bengtsson for discus- given quantum t-design having prescribed symme- sions about the restricted defect for SIC-POVM try [25] forms an isolated structure. We introduced in dimension 3. Financial support by Narodowe a powerful method that divides this full non-linear Centrum Nauki under the grant number DEC- problem into a linear problem and a simpler non- 2015/18/A/ST2/00274 and by the John Temple- linear problem (see Section III). ton Foundation under the project No. 56033 is Using our method, we have proven that known gratefully acknowledged. 9

Appendix A: Robustness of restricted defect Using this result, we have

† † † In this Appendix we derive the function f(d, N) |δ(R R)ij| = |δ(R )ijRij + (R )ijδ(R)ij| † which appears in Eq.(18) and allows us to show ≤ 2 max |δ(R )ijRij| that the restricted defect of a unitary matrix asso- Rij † ciated to a given generalised measurement is sta- ≤ 2 max |δ(R )ij| max |Rij| ble with respect to small perturbations. Consider Rij Rij s a set of vectors φj and the approximate vectors 26d5/2  2d2 N − d φ0 = φ + δφ . The perturbed Gram matrix is ≤ 1 − s. j j j N 2 N d(N − 1) given by

(G + δG)ij = hφi + δφi|φj + δφji Now we are in position to estimate the maximal perturbation on the eigenvalues of R†R ≈ hφi|φji + hδφi|φji + hφi|δφji, 0 † which implies that λ = λk + δλk ≈ λk + hδ(R R)i. (A3) √ |δGij| ≤ kφikkδφik + kφjkkδφjk ≤ 2 d s. (A1) From the Gerschgorin circle theorem [38] we have

Here, we used Eq.(16). Also, from U = − 2d G we † X † I N |hδ(R R)i| ≤ |δ(R R)ij|. (A4) 3/2 have |δUii| = 0 and |δUij| ≤ 4d s/N for i 6= j. ij Let us now calculate the perturbations on entries of the matrix R, which defines the system of equa- From combining Eqs.(A3), (A3) and (A4) we find tions (9). It is simple to show that if N > 2d the 0 that |λi − λi| ≤ f(d, N) s, where maximal perturbations are produced by the entries of R = −2U U , associated to the left term of s jk kk kj 26d5/2  2d2 N − d Eq.(9). Therefore f(d, N) = 1 − s. N 2 N d(N − 1)  2d 4d3/2 |δR | = 2|U ||δU | ≤ 2 1 − s ij kk ij N N Given that R†R is a positive operator, its eigenval- 8 d3/2  2d ues λi coincide with its singular values σi. There- ≤ 1 − s. (A2) fore |σ0 − σ | ≤ f(d, N) s, which proves Eq.(18). N N i i

[1] J. Renes, R. Blume-Kohout, A. Scott, C. Caves, [9] M. Ballester, S. Wehner, Entropic uncertainty re- Symmetric Informationally Complete Quantum lations and locking: tight bounds for mutually un- Measurements, J. Math. Phys. 45, 2171 (2004). biased bases, Phys. Rev. A 75, 022319 (2007). [2] C. A. Fuchs and M. C. Hoang and B. C. Stacey, [10] A. Scott, M. Grassl, SIC-POVMs: A new com- The SIC Question: History and State of Play, puter study, J. Math. Phys. 51, 042203 (2010). preprint arXiv:1703.07901 (2017) [11] S. Brierley, S. Weigert, I. Bengtsson, All Mutu- [3] I.D. Ivanovic, Geometrical description of quantal ally Unbiased Bases in Dimensions Two to Five, state determination, J. Phys. A, 14, 3241 (1981). Quantum Info. & Comp. 10, 0803-0820 (2010). [4] W. Wootters, B. Fields, Optimal State- [12] A. Sehrawat, A. Klimov, Unitarily inequivalent Determination by Mutually Unbiased Mea- mutually unbiased bases for n qubits, Phys. Rev. surements, Ann. Phys. 191, 363 (1989). A 90, 062308 (2014). [5] S. Wu, S. Yu, K. Mølmer, Entropic uncertainty [13] W. Kantor, MUBs inequivalence and affine planes, relation for mutually unbiased bases, Phys. Rev. J. Math. Phy. 53, 032204 (2012). A 79, 022104 (2009). [14] A. Szymusiak, Maximally informative ensembles [6] A. Rastegin, Uncertainty relations for MUBs and for SIC-POVMs in dimension 3, J. Phys. A: Math. SIC-POVMs in terms of generalized entropies, Theor. 47, 445301 (2014). Eur. Phys. J. D 67, 269 (2013). [15] J. Rehacek, Z. Hradil, A. Klimov, G. Leuchs, L. [7] C. Spengler, M. Huber, S. Brierley, T. Adakty- Sanchez-Soto, Sizing up entanglement in mutu- los, B. Hiesmayr, Entanglement detection via mu- ally unbiased bases with Fisher information, Phys. tually unbiased bases, Phys. Rev. A 86, 022311 Rev. A 88, 052110 (2013) (2012). [16] D. Mayers, A. Yao, Self testing quantum appara- [8] B. Hiesmayr, W. L¨offler,Mutually Unbiased Bases tus, Quant. Inf. & Comp. 4, 4, 273–286, (2004). and Bound Entanglement, Phys. Scr., 014017 [17] I. Bengtsson, K. Zyczkowski,˙ On discrete (2014). 10

structures in finite Hilbert spaces, preprint che Universit¨atBerlin, July 2006, available online: arXiv:1701.07902. http://magma.maths.usyd.edu.au/Magma2006/ [18] J. Tropp, I. Dhillon, R. Heath Jr., T. Strohmer, [29] I. Bengtsson, W. Bruzda, A. Ericsson, J. Lars- Designing structured tight frames via alternating son, W. Tadej, K. Zyczkowski.˙ Mutually unbiased projection, IEEE Trans. Info. Theory 51, 1, 188- bases and Hadamards of order six. J. Math. Phys. 209 (2005). 48, 052106 (2007). [19] W. Tadej, K. Zyczkowski,˙ A concise guide to com- [30] D. Goyeneche, Mutually unbiased triplets from plex Hadamard matrices, Open Syst. Inf. Dyn. 13, non-affine families of complex Hadamard matri- 133-177 (2006) ces in dimension six, J. Phys. A: Math. Theor. 46 [20] D. Goyeneche, O. Turek, Equiangular tight frames 105301 (2013). and unistochastic matrices, arXiv:1607.04528 [31] D. Goyeneche, S. Gomez, Mutually unbiased bases [math-ph] (2016). with free parameters, Phys. Rev. A 92, 062325 [21] W. Tadej, K. Zyczkowski,˙ Defect of a unitary ma- (2015). trix, Linear Alg. Appl., 429, 2-3, 447-481 (2008). [32] M. Fickus, D. Mixon, Tables of the existence [22] S. Chaturvedi, Aspects of mutually unbiased bases of equiangular tight frames, arXiv:1504.00253 in odd prime power dimensions, Phys. Rev. A 65, (2015). 044301 (2002). [33] S. Kochen and E. Specker, J. Math. Mech. 17, 59 [23] F. Sz¨oll˝osi,All complex equiangular tight frames (1967). in dimension 3, arXiv:1402.6429 [math.FA] (2014) [34] S. Yu, C. Oh, State-independent proof of Kochen- [24] S. Hoggar, 64 lines from a Quaternionic Polytope, Specker theorem with 13 rays, Phys. Rev. Lett. Geom. Dedicata 69 287 (1998). 108, 030402 (2012). [25] G. Zauner, Quantum Designs: Foundations of [35] A. Cabello, J. Estebaranz, G. Garcia Alcaine, a non-commutative Design Theory, PhD The- Bell-Kochen-Specker theorem: A proof with 18 sis, University of Vienna (1999). Available online vectors, Phys.Lett. A, 212, 183 (1996) at http://www.gerhardzauner.at/documents/gz- [36] P. Lisonek, P. Badziag, J. Portillo, A. Cabello, quantumdesigns.pdf Kochen-Specker set with seven contexts, Phys. [26] M. Grassl, On SIC-POVMs and MUBs in Dimen- Rev. A 89, 042101 (2014) sion 6, Proc. ERATO Conf. on Quantum Infor- [37] A. J. Scott, SICs: Extending the list of solutions, mation Science, Tokyo, 60-1, (2004) (arXiv:quant- arXiv:1703.03993 (quant-ph), 2017. ph/0406175) [38] S. Gerschgorin, Uber¨ die Abgrenzung der Eigen- [27] D. Appleby, SIC-POVMs and the Extended Clif- werte einer Matrix, Izv. Akad. Nauk. USSR Otd. ford Group, J. Math. Phys. 46 052107 (2005). Fiz.-Mat. Nauk, 6, 749-754 (1931). [28] M. Grassl, Finding equiangular lines in com- plex space, MAGMA 2006 Conference, Technis-