Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2602.22336

Basis-independent stabilizerness and maximally noisy magic states

Michael Zurel, Jack Davis

quant-phSubmitted Feb 25, 2026• Announced Feb 27, 2026• 0 citations• Revised Sep 30, 2026 (v2)
v1 · Feb 25, 2026v2 · Sep 30, 2026
arXivHTMLPDF

Abstract

Absolutely stabilizer states are those that remain convex mixtures of stabilizer states after conjugation by any unitary. Here we give a characterization of such states for multiple qudits of all prime dimensions by introducing a polytope of their allowed spectra. We illustrate this through the examples of one qubit, two qubits, and one qutrit. In particular, the set of absolutely stabilizer states for a single qubit is a ball inscribed in the stabilizer octahedron, but for higher dimensions the geometry is more complicated. For odd-prime-dimensional qudits, we also give a complete characterization of absolutely Wigner-positive states, i.e., states whose Wigner function remains nonnegative after conjugation by any unitary. In so doing, we show there are absolutely Wigner-positive states that are not absolutely stabilizer, which can be seen as a unitarily invariant version of bound magic. We then study the radii of the largest balls contained in the sets of absolutely stabilizer states and absolutely Wigner-positive states. These radii respectively tell us the lowest possible purity of nonstabilizer and Wigner-negative states. Conversely, we also find the radius of the smallest ball containing the set of absolutely Wigner-positive states, giving a tight purity-based necessary condition thereof. Finally, we apply our methods to magic witnessing and find simple sufficient conditions for determining if an $n$-qubit state is reducible to an $m$-qubit magic state for any $m<n$.

Comments
28 pages, 11 figures

Related papers

A Phase-Space Geometric Measure of Magic in Qubit Systems

Soumyojyoti Dutta, Tushar

Magic -- the resource enabling quantum computational advantage beyond stabilizer circuits -- has a clean phase-space characterization in odd prime dimensions that qubits notoriously lack. We study C(rho), the l_1 distance from a state's discrete Wigner function to the stabilizer polytope, and determine its exact geometry. For a single qubit with Bloch vector r we prove the closed form C = max(0, ||r||_1 - 1)/2, so the maximum is (sqrt(3)-1)/2, attained at the eight face states. For products of single-qubit states, C(rho_1 x ... x rho_n) <= prod_i (1 + C(rho_i)) - 1, with equality whenever every factor has Bloch sign product r_x r_y r_z <= 0; this includes the exact tensor powers ((1+sqrt(3))/2)^n - 1 and the exact maximum of C over fully separable n-qubit states. Because no discrete Wigner function for qubits is Clifford covariant, any such measure is frame dependent, and we identify which of its features are not. For a single qubit there are exactly two valid frames up to phase-space translation and C is identical on both. Pauli-covariant two-qubit frames number 2048 up to translation, and across them the ratio C(rho_Rx)/C(rho_Ry) takes the values 1/2, 1 and 2. In every one of these frames some stabilizer state has negative Wigner entries, a consequence of the Mermin-Peres square, so the two-qubit bound Gamma >= 1 + C/2 against the robustness of magic Gamma is frame independent. Within the Wootters frame we prove integer values 1, 2, 1 of the tightness ratio kappa := (Gamma-1)/C for three families in the two-qubit repetition codespace, with explicit certificates. For n-qubit GHZ phase states we prove Gamma = |sin phi| + |cos phi| for all n and C >= 2^(floor(n/2)-1) (|sin phi| + |cos phi| - 1), tight for n = 2 and, by exact linear programming, for n = 3, 4, so kappa decays exponentially. Finally, C is not a magic monotone, so it cannot replace Gamma in distillation bounds.

quant-pharXiv:2603.20792•Mar 24, 2026•0 citations

Geometric optimization problems in quantum computation and discrete mathematics: Stabilizer states and lattices

Arne Heimendahl

This thesis consists of two parts: Part I deals with properties of stabilizer states and their convex hull, the stabilizer polytope. Stabilizer states, Pauli measurements and Clifford unitaries are the three building blocks of the stabilizer formalism whose computational power is limited by the Gottesman- Knill theorem. This model is usually enriched by a magic state to get a universal model for quantum computation, referred to as quantum computation with magic states (QCM). The first part of this thesis will investigate the role of stabilizer states within QCM from three different angles. The first considered quantity is the stabilizer extent, which provides a tool to measure the non-stabilizerness or magic of a quantum state. It assigns a quantity to each state roughly measuring how many stabilizer states are required to approximate the state. It has been shown that the extent is multiplicative under taking tensor products when the considered state is a product state whose components are composed of maximally three qubits. In Chapter 2, we will prove that this property does not hold in general, more precisely, that the stabilizer extent is strictly submultiplicative. We obtain this result as a consequence of rather general properties of stabilizer states. Informally our result implies that one should not expect a dictionary to be multiplicative under taking tensor products whenever the dictionary size grows subexponentially in the dimension. In Chapter 3, we consider QCM from a resource theoretic perspective. The resource theory of magic is based on two types of quantum channels, completely stabilizer preserving maps and stabilizer operations. Both classes have the property that they cannot generate additional magic resources. We will show that these two classes of quantum channels do not coincide, specifically, that stabilizer operations are a strict subset of the set of completely stabilizer preserving channels. This might have the consequence that certain tasks which are usually realized by stabilizer operations could in principle be performed better by completely stabilizer preserving maps. In Chapter 4, the last one of Part I, we consider QCM via the polar dual stabilizer polytope (also called the Lambda-polytope). This polytope is a superset of the quantum state space and every quantum state can be written as a convex combination of its vertices. A way to classically simulate quantum computing with magic states is based on simulating Pauli measurements and Clifford unitaries on the vertices of the Lambda-polytope. The complexity of classical simulation with respect to the polytope is determined by classically simulating the updates of vertices under Clifford unitaries and Pauli measurements. However, a complete description of this polytope as a convex hull of its vertices is only known in low dimensions (for up to two qubits or one qudit when odd dimensional systems are considered). We make progress on this question by characterizing a certain class of operators that live on the boundary of the Lambda-polytope when the underlying dimension is an odd prime. This class encompasses for instance Wigner operators, which have been shown to be vertices of Lambda. We conjecture that this class contains even more vertices of Lambda. Eventually, we will shortly sketch why applying Clifford unitaries and Pauli measurements to this class of operators can be efficiently classically simulated. Part II of this thesis deals with lattices. Lattices are discrete subgroups of the Euclidean space. They occur in various different areas of mathematics, physics and computer science. We will investigate two types of optimization problems related to lattices. In Chapter 6 we are concerned with optimization within the space of lattices. That is, we want to compare the Gaussian potential energy of different lattices. To make the energy of lattices comparable we focus on lattices with point density one. In particular, we focus on even unimodular lattices and show that, up to dimension 24, they are all critical for the Gaussian potential energy. Furthermore, we find that all n-dimensional even unimodular lattices with n 24 are local minima or saddle points. In contrast in dimension 32, there are even unimodular lattices which are local maxima and others which are not even critical. In Chapter 7 we consider flat tori R^n/L, where L is an n-dimensional lattice. A flat torus comes with a metric and our goal is to approximate this metric with a Hilbert space metric. To achieve this, we derive an infinite-dimensional semidefinite optimization program that computes the least distortion embedding of the metric space R^n/L into a Hilbert space. This program allows us to make several interesting statements about the nature of least distortion embeddings of flat tori. In particular, we give a simple proof for a lower bound which gives a constant factor improvement over the previously best lower bound on the minimal distortion of an embedding of an n-dimensional flat torus. Furthermore, we show that there is always an optimal embedding into a finite-dimensional Hilbert space. Finally, we construct optimal least distortion embeddings for the standard torus R^n/Z^n and all 2-dimensional flat tori.

•Jan 1, 2023•0 citations

On the geometry and typicality of quantum magic

Zhenhuan Liu, Zi-Wen Liu

We prove that, for an $n$-qubit system of dimension $d=2^n$, every state satisfying $\operatorname{Tr}(ρ^2)\le 1/(d-a_\ast)$, with $a_\ast=0.458327\cdots$, lies inside the stabilizer polytope and is therefore magic-free. Combining this result with general geometric properties of high-dimensional polytopes, we establish quantitative estimates for the Hilbert--Schmidt inradius and volume radius of the stabilizer polytope, and use them to characterize the typicality of magic in random induced states obtained by tracing out a $k$-dimensional subsystem from a $d\times k$-dimensional Haar-random pure state. We prove a sharp phase transition in the probability of such states having magic, whose transition dimension $k_\star$ is bounded between $Ω(d^2/\log^2d)$ and $\mathcal{O}(d^2)$. We further prove that the number of facets of the stabilizer polytope lies between $\exp[Ω(d^2/\log^2 d)]$ and $\exp[\mathcal{O}(d^2\log^2 d)]$ employing a result of Bourgain and Milman in convex geometry, substantially improving upon the previous quasipolynomial lower bound and implying that doubly-exponentially many linear inequalities in the number of qubits are required for an exact description of the magic-free region. Overall, our results reveal the near-extremal geometry of the stabilizer polytope and provide a quantitative foundation for understanding the typicality, robustness, and detectability of magic.

quant-pharXiv:2609.03944•Sep 4, 2026•0 citations

Polytopes of Absolutely Wigner Bounded Spin States

Jérôme Denis, J. Davis, Robert B. Mann +1 more

Quasiprobability has become an increasingly popular notion for characterising non-classicality in quantum information, thermodynamics, and metrology. Two important distributions with non-positive quasiprobability are the Wigner function and the Glauber-Sudarshan function. Here we study properties of the spin Wigner function for finite-dimensional quantum systems and draw comparisons with its infinite-dimensional analog, focusing in particular on the relation to the Glauber-Sudarshan function and the existence of absolutely Wigner-bounded states. More precisely, we investigate unitary orbits of mixed spin states that are characterized by Wigner functions lower-bounded by a specified value. To this end, we extend a characterization of the set of absolutely Wigner positive states as a set of linear eigenvalue constraints, which together define a polytope centred on the maximally mixed state in the simplex of spin- j states. The lower bound determines the relative size of such absolutely Wigner bounded (AWB) polytopes and we study their geometric characteristics. In each dimension a Hilbert-Schmidt ball representing a tight purity-based sufficient condition to be AWB is exactly determined, while another ball representing a necessary condition to be AWB is conjectured. Special attention is given to the case where the polytope separates orbits containing only positive Wigner functions from other orbits because of the use of Wigner negativity as a witness of non-classicality. Comparisons are made to absolute symmetric state separability and spin Glauber-Sudarshan positivity, with additional details given for low spin quantum numbers.

arXiv:2304.09006•Quantum•Dec 4, 2024•3 citations

The geometry of absolute separability and other convex matrix properties from spectrum

Jennifer Ahiable, Naga Bhavya Teja Kothakonda, Andreas Winter

We investigate the geometric structure of the set of spectra of bipartite absolute separable states ($\mathrm{ASEP}_{m,n}$) and absolute positive partial transpose states ($\mathrm{APPT}_{m,n}$), \emph{i.e.}, bipartite quantum states that remain separable or PPT, respectively, under all global unitary transformations. First, we establish general geometric properties of absolute convex sets of matrices, their spectra and extreme points. Regarding absolute separability, we present a permutation-symmetric reformulation of the absolute PPT criterion and use it to demonstrate that $\mathrm{APPT}_{m,n}$ is a spectrahedron for all $m\leq n$: in particular, all its faces are exposed. While $\mathrm{ASEP}_{2,n}$ is a spectrahedron, too, for $\mathrm{ASEP}_{m,n}$ we can in general only show that it is a semialgebraic set for all $m,n\geq 3$. Furthermore, we provide a complete characterization of the faces and extreme points of $\mathrm{APPT}_{m,n}$ and demonstrate that the dimension of a face is determined by the rank of a certain matrix, with maximal proper faces having dimension $mn-m-1$. In the quantitative setting, we provide a rigorous lower bound on the maximal attainable purity of $\mathrm{APPT}_{m,n}$ via an inscribed polytope $\mathcal{P}_{m,n}$ and conjecture that the maximal purity of $\mathrm{APPT}_{m,n}$ (along with its spectra) coincides with the polytope for arbitrary dimensions except when $m=n=2$. Additionally, we also provide a rigorous upper bound on the minimal von Neumann entropy of $\mathrm{APPT}_{m,n}$ and demonstrate numerically that the minimum entropy eventually coincides with the polytope $\mathcal{P}_{m,n}$ as the local system dimension $n$ increases. Finally, we show that the relative spectral volume of $\mathrm{APPT}_{m,n}$ decays exponentially in $n$ by a constant multiplicative factor of the relative volume of the inscribed polytope $\mathcal{P}_{m,n}$.

quant-pharXiv:2608.03390•Aug 5, 2026•0 citations

Wigner’s theorem for stabilizer states and quantum designs

Valentin Obst, Arne Heimendahl, Tanmay Singal +1 more

We describe the symmetry group of the stabilizer polytope for any number n of systems and any prime local dimension d. In the qubit case, the symmetry group coincides with the linear and anti-linear Clifford operations. In the case of qudits, the structure is somewhat richer: For n = 1, it is a wreath product of permutations of bases and permutations of the elements within each basis. For n > 1, the symmetries are given by affine symplectic similitudes. These are the affine maps that preserve the symplectic form of the underlying discrete phase space up to a non-zero multiplier. We phrase these results with respect to a number of a priori different notions of “symmetry,” including Kadison symmetries (bijections that are compatible with convex combinations), Wigner symmetries (bijections that preserve inner products), and symmetries realized by an action on Hilbert space. Going beyond stabilizer states, we extend an observation of Heinrich and Gross [Quantum 3, 132 (2019)] and show that the symmetries of fairly general sets of Hermitian operators are constrained by certain moments. In particular: the symmetries of a set that behaves like a 3-design preserve Jordan products and are therefore realized by conjugation with unitaries or anti-unitaries. (The structure constants of the Jordan algebra are encoded in an order-three tensor, which we connect to the third moments of a design). This generalizes Kadison’s formulation of the classic Wigner theorem on quantum mechanical symmetries.

arXiv:2405.17565•Journal of Mathematical Physics•Nov 1, 2024•3 citations
Absolutely stabilizer states are those that remain convex mixtures of stabilizer states after conjugation by any unitary. Here we give a characterization of such states for multiple qudits of all prime dimensions by introducing a polytope of their allowed spectra. We illustrate this through the examples of one qubit, two qubits, and one qutrit. In particular, the set of absolutely stabilizer states for a single qubit is a ball inscribed in the stabilizer octahedron, but for higher dimensions the geometry is more complicated. For odd-prime-dimensional qudits, we also give a complete characterization of absolutely Wigner-positive states, i.e., states whose Wigner function remains nonnegative after conjugation by any unitary. In so doing, we show there are absolutely Wigner-positive states that are not absolutely stabilizer, which can be seen as a unitarily invariant version of bound magic. We then study the radii of the largest balls contained in the sets of absolutely stabilizer states and absolutely Wigner-positive states. These radii respectively tell us the lowest possible purity of nonstabilizer and Wigner-negative states. Conversely, we also find the radius of the smallest ball containing the set of absolutely Wigner-positive states, giving a tight purity-based necessary condition thereof. Finally, we apply our methods to magic witnessing and find simple sufficient conditions for determining if an nnn-qubit state is reducible to an mmm-qubit magic state for any m<nm<nm<n.
Magic -- the resource enabling quantum computational advantage beyond stabilizer circuits -- has a clean phase-space characterization in odd prime dimensions that qubits notoriously lack. We study C(rho), the l_1 distance from a state's discrete Wigner function to the stabilizer polytope, and determine its exact geometry. For a single qubit with Bloch vector r we prove the closed form C = max(0, ||r||_1 - 1)/2, so the maximum is (sqrt(3)-1)/2, attained at the eight face states. For products of single-qubit states, C(rho_1 x ... x rho_n) <= prod_i (1 + C(rho_i)) - 1, with equality whenever every factor has Bloch sign product r_x r_y r_z <= 0; this includes the exact tensor powers ((1+sqrt(3))/2)^n - 1 and the exact maximum of C over fully separable n-qubit states. Because no discrete Wigner function for qubits is Clifford covariant, any such measure is frame dependent, and we identify which of its features are not. For a single qubit there are exactly two valid frames up to phase-space translation and C is identical on both. Pauli-covariant two-qubit frames number 2048 up to translation, and across them the ratio C(rho_Rx)/C(rho_Ry) takes the values 1/2, 1 and 2. In every one of these frames some stabilizer state has negative Wigner entries, a consequence of the Mermin-Peres square, so the two-qubit bound Gamma >= 1 + C/2 against the robustness of magic Gamma is frame independent. Within the Wootters frame we prove integer values 1, 2, 1 of the tightness ratio kappa := (Gamma-1)/C for three families in the two-qubit repetition codespace, with explicit certificates. For n-qubit GHZ phase states we prove Gamma = |sin phi| + |cos phi| for all n and C >= 2^(floor(n/2)-1) (|sin phi| + |cos phi| - 1), tight for n = 2 and, by exact linear programming, for n = 3, 4, so kappa decays exponentially. Finally, C is not a magic monotone, so it cannot replace Gamma in distillation bounds.
This thesis consists of two parts: Part I deals with properties of stabilizer states and their convex hull, the stabilizer polytope. Stabilizer states, Pauli measurements and Clifford unitaries are the three building blocks of the stabilizer formalism whose computational power is limited by the Gottesman- Knill theorem. This model is usually enriched by a magic state to get a universal model for quantum computation, referred to as quantum computation with magic states (QCM). The first part of this thesis will investigate the role of stabilizer states within QCM from three different angles. The first considered quantity is the stabilizer extent, which provides a tool to measure the non-stabilizerness or magic of a quantum state. It assigns a quantity to each state roughly measuring how many stabilizer states are required to approximate the state. It has been shown that the extent is multiplicative under taking tensor products when the considered state is a product state whose components are composed of maximally three qubits. In Chapter 2, we will prove that this property does not hold in general, more precisely, that the stabilizer extent is strictly submultiplicative. We obtain this result as a consequence of rather general properties of stabilizer states. Informally our result implies that one should not expect a dictionary to be multiplicative under taking tensor products whenever the dictionary size grows subexponentially in the dimension. In Chapter 3, we consider QCM from a resource theoretic perspective. The resource theory of magic is based on two types of quantum channels, completely stabilizer preserving maps and stabilizer operations. Both classes have the property that they cannot generate additional magic resources. We will show that these two classes of quantum channels do not coincide, specifically, that stabilizer operations are a strict subset of the set of completely stabilizer preserving channels. This might have the consequence that certain tasks which are usually realized by stabilizer operations could in principle be performed better by completely stabilizer preserving maps. In Chapter 4, the last one of Part I, we consider QCM via the polar dual stabilizer polytope (also called the Lambda-polytope). This polytope is a superset of the quantum state space and every quantum state can be written as a convex combination of its vertices. A way to classically simulate quantum computing with magic states is based on simulating Pauli measurements and Clifford unitaries on the vertices of the Lambda-polytope. The complexity of classical simulation with respect to the polytope is determined by classically simulating the updates of vertices under Clifford unitaries and Pauli measurements. However, a complete description of this polytope as a convex hull of its vertices is only known in low dimensions (for up to two qubits or one qudit when odd dimensional systems are considered). We make progress on this question by characterizing a certain class of operators that live on the boundary of the Lambda-polytope when the underlying dimension is an odd prime. This class encompasses for instance Wigner operators, which have been shown to be vertices of Lambda. We conjecture that this class contains even more vertices of Lambda. Eventually, we will shortly sketch why applying Clifford unitaries and Pauli measurements to this class of operators can be efficiently classically simulated. Part II of this thesis deals with lattices. Lattices are discrete subgroups of the Euclidean space. They occur in various different areas of mathematics, physics and computer science. We will investigate two types of optimization problems related to lattices. In Chapter 6 we are concerned with optimization within the space of lattices. That is, we want to compare the Gaussian potential energy of different lattices. To make the energy of lattices comparable we focus on lattices with point density one. In particular, we focus on even unimodular lattices and show that, up to dimension 24, they are all critical for the Gaussian potential energy. Furthermore, we find that all n-dimensional even unimodular lattices with n 24 are local minima or saddle points. In contrast in dimension 32, there are even unimodular lattices which are local maxima and others which are not even critical. In Chapter 7 we consider flat tori R^n/L, where L is an n-dimensional lattice. A flat torus comes with a metric and our goal is to approximate this metric with a Hilbert space metric. To achieve this, we derive an infinite-dimensional semidefinite optimization program that computes the least distortion embedding of the metric space R^n/L into a Hilbert space. This program allows us to make several interesting statements about the nature of least distortion embeddings of flat tori. In particular, we give a simple proof for a lower bound which gives a constant factor improvement over the previously best lower bound on the minimal distortion of an embedding of an n-dimensional flat torus. Furthermore, we show that there is always an optimal embedding into a finite-dimensional Hilbert space. Finally, we construct optimal least distortion embeddings for the standard torus R^n/Z^n and all 2-dimensional flat tori.
We prove that, for an nnn-qubit system of dimension d=2nd=2^nd=2n, every state satisfying Tr⁡(ρ2)≤1/(d−a∗)\operatorname{Tr}(ρ^2)\le 1/(d-a_\ast)Tr(ρ2)≤1/(d−a∗​), with a∗=0.458327⋯a_\ast=0.458327\cdotsa∗​=0.458327⋯, lies inside the stabilizer polytope and is therefore magic-free. Combining this result with general geometric properties of high-dimensional polytopes, we establish quantitative estimates for the Hilbert--Schmidt inradius and volume radius of the stabilizer polytope, and use them to characterize the typicality of magic in random induced states obtained by tracing out a kkk-dimensional subsystem from a d×kd\times kd×k-dimensional Haar-random pure state. We prove a sharp phase transition in the probability of such states having magic, whose transition dimension k⋆k_\stark⋆​ is bounded between Ω(d2/log⁡2d)Ω(d^2/\log^2d)Ω(d2/log2d) and O(d2)\mathcal{O}(d^2)O(d2). We further prove that the number of facets of the stabilizer polytope lies between exp⁡[Ω(d2/log⁡2d)]\exp[Ω(d^2/\log^2 d)]exp[Ω(d2/log2d)] and exp⁡[O(d2log⁡2d)]\exp[\mathcal{O}(d^2\log^2 d)]exp[O(d2log2d)] employing a result of Bourgain and Milman in convex geometry, substantially improving upon the previous quasipolynomial lower bound and implying that doubly-exponentially many linear inequalities in the number of qubits are required for an exact description of the magic-free region. Overall, our results reveal the near-extremal geometry of the stabilizer polytope and provide a quantitative foundation for understanding the typicality, robustness, and detectability of magic.
Quasiprobability has become an increasingly popular notion for characterising non-classicality in quantum information, thermodynamics, and metrology. Two important distributions with non-positive quasiprobability are the Wigner function and the Glauber-Sudarshan function. Here we study properties of the spin Wigner function for finite-dimensional quantum systems and draw comparisons with its infinite-dimensional analog, focusing in particular on the relation to the Glauber-Sudarshan function and the existence of absolutely Wigner-bounded states. More precisely, we investigate unitary orbits of mixed spin states that are characterized by Wigner functions lower-bounded by a specified value. To this end, we extend a characterization of the set of absolutely Wigner positive states as a set of linear eigenvalue constraints, which together define a polytope centred on the maximally mixed state in the simplex of spin- j states. The lower bound determines the relative size of such absolutely Wigner bounded (AWB) polytopes and we study their geometric characteristics. In each dimension a Hilbert-Schmidt ball representing a tight purity-based sufficient condition to be AWB is exactly determined, while another ball representing a necessary condition to be AWB is conjectured. Special attention is given to the case where the polytope separates orbits containing only positive Wigner functions from other orbits because of the use of Wigner negativity as a witness of non-classicality. Comparisons are made to absolute symmetric state separability and spin Glauber-Sudarshan positivity, with additional details given for low spin quantum numbers.
We investigate the geometric structure of the set of spectra of bipartite absolute separable states (ASEPm,n\mathrm{ASEP}_{m,n}ASEPm,n​) and absolute positive partial transpose states (APPTm,n\mathrm{APPT}_{m,n}APPTm,n​), \emph{i.e.}, bipartite quantum states that remain separable or PPT, respectively, under all global unitary transformations. First, we establish general geometric properties of absolute convex sets of matrices, their spectra and extreme points. Regarding absolute separability, we present a permutation-symmetric reformulation of the absolute PPT criterion and use it to demonstrate that APPTm,n\mathrm{APPT}_{m,n}APPTm,n​ is a spectrahedron for all m≤nm\leq nm≤n: in particular, all its faces are exposed. While ASEP2,n\mathrm{ASEP}_{2,n}ASEP2,n​ is a spectrahedron, too, for ASEPm,n\mathrm{ASEP}_{m,n}ASEPm,n​ we can in general only show that it is a semialgebraic set for all m,n≥3m,n\geq 3m,n≥3. Furthermore, we provide a complete characterization of the faces and extreme points of APPTm,n\mathrm{APPT}_{m,n}APPTm,n​ and demonstrate that the dimension of a face is determined by the rank of a certain matrix, with maximal proper faces having dimension mn−m−1mn-m-1mn−m−1. In the quantitative setting, we provide a rigorous lower bound on the maximal attainable purity of APPTm,n\mathrm{APPT}_{m,n}APPTm,n​ via an inscribed polytope Pm,n\mathcal{P}_{m,n}Pm,n​ and conjecture that the maximal purity of APPTm,n\mathrm{APPT}_{m,n}APPTm,n​ (along with its spectra) coincides with the polytope for arbitrary dimensions except when m=n=2m=n=2m=n=2. Additionally, we also provide a rigorous upper bound on the minimal von Neumann entropy of APPTm,n\mathrm{APPT}_{m,n}APPTm,n​ and demonstrate numerically that the minimum entropy eventually coincides with the polytope Pm,n\mathcal{P}_{m,n}Pm,n​ as the local system dimension nnn increases. Finally, we show that the relative spectral volume of APPTm,n\mathrm{APPT}_{m,n}APPTm,n​ decays exponentially in nnn by a constant multiplicative factor of the relative volume of the inscribed polytope Pm,n\mathcal{P}_{m,n}Pm,n​.
We describe the symmetry group of the stabilizer polytope for any number n of systems and any prime local dimension d. In the qubit case, the symmetry group coincides with the linear and anti-linear Clifford operations. In the case of qudits, the structure is somewhat richer: For n = 1, it is a wreath product of permutations of bases and permutations of the elements within each basis. For n > 1, the symmetries are given by affine symplectic similitudes. These are the affine maps that preserve the symplectic form of the underlying discrete phase space up to a non-zero multiplier. We phrase these results with respect to a number of a priori different notions of “symmetry,” including Kadison symmetries (bijections that are compatible with convex combinations), Wigner symmetries (bijections that preserve inner products), and symmetries realized by an action on Hilbert space. Going beyond stabilizer states, we extend an observation of Heinrich and Gross [Quantum 3, 132 (2019)] and show that the symmetries of fairly general sets of Hermitian operators are constrained by certain moments. In particular: the symmetries of a set that behaves like a 3-design preserve Jordan products and are therefore realized by conjugation with unitaries or anti-unitaries. (The structure constants of the Jordan algebra are encoded in an order-three tensor, which we connect to the third moments of a design). This generalizes Kadison’s formulation of the classic Wigner theorem on quantum mechanical symmetries.