Stabilizer entropy of quantum tetrahedra

Abstract

How complex is the structure of quantum geometry? In several approaches, the spacetime atoms are obtained by the $SU(2)$ intertwiner called quantum tetrahedron. The complexity of this construction has a concrete consequence in recent efforts to simulate such models and toward experimental demonstrations of quantum gravity effects. There are, therefore, both a computational and an experimental complexity inherent to this class of models. In this paper, we study this complexity under the lens of Stabilizer Entropy (SE). We calculate the SE of the gauge-invariant basis states and its average in the $SU(2)-$gauge invariant subspace. We find that the states of definite volume are singled out by the (near) maximal SE and give precise bounds to the verification protocols for experimental demonstrations on available quantum computers.

I Introduction

Understanding whether gravity admits a quantum formulation is one of the most intriguing challenges of modern physics. In the last decade, quantum information theory has provided new conceptual and mathematical tools to investigate the structure of spacetime at the quantum scale. Entanglement entropy has been used to probe the holographic architecture of spacetime, supporting the idea of entanglement as an essential resource to the emergence of classical spacetime geometry [1, 2, 3, 4, 5]. Recently, also fostered by new perspectives in quantum gravity phenomenology, the use of quantum information tools to design and investigate experimental evidence for quantum features of the gravitational field has attracted much attention [8, 9, 10].

Within the limits of current experimental technology, the first widely available quantum computers today allow to simulate quantum gravity states, providing suggestions, predictions and setup ideas for future experiments [11, 12, 13]. In this scenario, we expect non-stabilizer resources to play a double key role. Gently speaking, non-stabilizerness is a core property of quantum states describing the complexity of the expression of their density operator in a specific operator basis (in the case of qubit systems, the Pauli operator basis), and its interplay with entanglement is known to be the essential ingredient needed to unlock quantum advantage [14, 15]. Recently, this resource has been shown to be given an entropic meaning, as Stabilizer Entropy (SE), making it both computable [16] and measurable [17]. SE directly affects the cost (in terms of classical resources) of simulating a quantum state or process: a $n$-qubit state or circuit using a number $t$ of non-stabilizer resources can be simulated with a classical computer at a computational cost that scales as $\exp(t)\operatorname{poly}(n)$[18]. In particular, it provides bounds on the fidelity reachable in experimental realizations of quantum states [19]. Moreover, SE is involved in the onset of universal, complex patterns of entanglement [20, 21], quantum chaos [22, 23, 24, 25], complexity in the wave function of quantum many-body systems [26, 27], and decoding algorithms from black hole’s Hawking radiation [28, 29, 30, 31]. States and processes with high non-stabilizer resources are generically exponentially harder to simulate on a classical computer, and they are harder to certificate in experimental protocols. An analysis of such resources is therefore vital in order to assess the simulability of quantum gravity states. At the same time, we expect such property to provide a new tool, in addition to entanglement, to investigate the emergence of classical spacetime in quantum gravity.

In this work, we explore the novel direction of looking at the non-stabilizerness of quantum gravity states in the setting of non-perturbative theories. We consider a description of quantum geometry given in terms of spin network states, a general tool shared by lattice gauge theory [32] and several background-independent approaches to quantum gravity (like loop quantum gravity [33, 34], state sum models [35, 36], and group field theories [37, 38, 39]), where they provide a gauge-invariant basis for the field. A spin network is represented by a graph $\Gamma$, with edges and nodes colored respectively by $SU(2)$ spin halves and intertwiner operators [34].

Each node of the network graph is dual to a quantum polyhedron geometry, with a number of faces equal to the valence of the node. We focus on a single $4-$valent node, that is an $SU(2)$-gauge invariant state corresponding to a quantum tetrahedron [40, 41, 42, 43].

We study the non-stabilizerness of quantum tetrahedron states using Stabilizer Entropy (SE) [16]. The first result of the work is that the states that diagonalize the oriented volume operator are the ones with highest value of SE: this result provides a new lower bound for the number of preparations needed in future experimental setups of quantum gravity states. To our knowledge, at present, this bound, obtained from the calculated intertwiner states non-stabilizerness has not been fully reached yet (see e.g. the data from the experiments realized in [11]). As a second result, we show that the projection into the gauge invariant Hilbert space associated to the process of constructing the quantum tetrahedron out of a collection of four qubits inherently requires non-stabilizer resources. Such resources become an intrinsic feature of the quantum geometry state, reflecting in the computational complexity of a simulation of such processes, in a way that is ultimately dependent on the structure of the gauge invariant space itself.

The paper is organized as follows. Section II provides the basic notions of the stabilizer formalism necessary for our analysis. We recall the definition of Pauli operators, the Clifford group and construct the set of pure stabilizer states, highlighting the necessity to go beyond this set of states in order to unlock quantum advantage. Then, we introduce the definition of Stabilizer Entropy as an entropic measure of nonstabilizerness of a pure quantum state, as well as its properties. In Section III we introduce the setting of quantum gravity. We realize a quantum tetrahedron via projection into the SU$(2)$ gauge invariant (intertwiner) subspace of a spin network Hilbert space. We show the most general intertwiner state for any SU$(2)$ spin$-j$ irrep and then focus on the case of $j=1/2$. In Section IV and Section VI we compute the non-stabilizerness of the logical basis and of the volume eigenstates basis elements. Then, these numerical results provide an estimate of the upper bound of the fidelity of the experiment in [11]; we compare our estimations with the experimental fidelity obtained. In Section V we extend the analysis to non-stabilizerness of subspaces, in order to investigate the cost of projection in terms of non-stabilizer resources. To this end, we introduce the average SE gap onto a subspace, and we show that this quantity is directly dependent from the internal structure of said subspace in the form of its projector. Finally, we apply our obtained results to the intertwiner subspace and conclude that imposing the $SU(2)-$gauge invariance has an intrinsic cost in terms of non-stabilizer resources.

II Stabilizer formalism and Stabilizer Entropy

In this section, we review the stabilizer formalism and its role in the quantum computation framework.

Let $\mathcal{H}\simeq\mathbb{C}^{2\otimes n}$ a $n$-qubit system and $\mathbb{P}_{n}$ be the Pauli group acting on $\mathcal{H}$. Define the Clifford group $\mathcal{C}(n)\subset\mathcal{U}(n)$ as the normalizer of the Pauli group, namely, $\mathcal{C}(n):=\{C\in\mathcal{U}(n)\,,\,\text{s.t. }\,\forall P\in\mathbb{P}_{n}\,,CPC^{{\dagger}}=P^{\prime}\in\mathbb{P}_{n}\}$ [44]. Hence, given a computational basis $\{\ket{i}\}$ of $\mathcal{H}$ as the common eigenbasis of the operators belonging to $\mathcal{Z}_{2}:=\{\mathbbm{1},Z\}^{\otimes n}$, one can define the set of pure stabilizer states of $\mathcal{H}$ as the full Clifford orbit of $\{\ket{i}\}$[45], namely

$${\operatorname{STAB}}=\{C\ket{i}\,,C\in\mathcal{C}(n)\}\,.$$ (1)

Stabilizer states share some properties with regard to the computational complexity of simulating quantum processes using classical resources; these properties are summarized by the Gottesman-Knill theorem, which states that any quantum process that can be represented with initial stabilizer states upon which one performs (i) Clifford unitaries, (ii) measurements of Pauli operators, (iii) Clifford operations conditioned on classical randomness, can be perfectly simulated by a classical computer in polynomial time [44]. This means that stabilizer states and Clifford operators are not actually “quantum" from a computational perspective, since they do not provide any advantage over classical computers. Since the set of stabilizer states is by definition closed under Clifford operations, a certain amount of resources beyond the Clifford group is needed to prepare a generic state in the Hilbert space: this quantity is referred to as non-stabilizerness of this state, which has been proven to be a useful resource for universal quantum computation [46] and for which several measures have been proposed [47, 45]. For our analysis, we are going to use two entropic non-stabilizerness measures called 2-Stabilizer Rényi Entropy (SE) [16] and its linear counterpart. They are defined starting from the probability distribution $\Xi_{P}(\ket{\psi}):=d^{-1}\operatorname{tr}^{2}(P\ket{\psi}\bra{\psi})$, with $P\in\mathbb{P}_{n}$ and $d=\operatorname{dim}(\mathcal{H})=2^{n}$, associated to the tomography of the quantum state $\psi$. Then, the $2$-Stabilizer Rényi Entropy for pure states is defined as

$$\begin{split}M_{2}(\ket{\psi})&:=-\log{d\|\Xi_{P}(\ket{\psi})\|^{2}_{2}}\\ &=-\log{d^{-1}\sum_{P\in\mathbb{P}_{n}}\operatorname{tr}^{4}(P\ket{\psi}\bra{\psi})}\,,\end{split}$$ (2)

whereas the linear SE is defined as

$$\begin{split}M_{\operatorname{lin}}(\ket{\psi})&:=1-d\|\Xi_{P}(\ket{\psi})\|^{2}_{2}\\ &=1-d^{-1}\sum_{P\in\mathbb{P}_{n}}\operatorname{tr}^{4}(P\ket{\psi}\bra{\psi})\,.\end{split}$$ (3)

Both $M_{2}$ and $M_{\operatorname{lin}}$ are: (i) faithful, i.e. $M(\psi)=0\Leftrightarrow\psi\in{\operatorname{STAB}}$, otherwise $M(\psi)>0$; (ii) invariant under Clifford operators, namely $M(C\psi C^{\dagger})=M(\psi)$, whereas $M_{2}$ is also additive under tensor product of quantum states [16]. Both measures can also be written in a more compact form:

$$\begin{split}M_{2}(\ket{\psi})&=-\log{d\operatorname{tr}Q\psi^{\otimes 4}}\,,\\ M_{\operatorname{lin}}(\ket{\psi})&=1-d\operatorname{tr}Q\psi^{\otimes 4}\,,\end{split}$$ (4)

with $Q:=d^{-2}\sum_{P\in\mathbb{P}_{n}}P^{\otimes 4}$. This construction can be generalized in a straightforward way to the Pauli group (namely, the Heisenberg-Weyl group) for qudits, that is, $l-$level systems [48], and the SE are defined exactly as in Eq.(4), see Appendix C for details.

Stabilizer entropies quantify the computational complexity of qubit states by the entropy of the distribution over the Pauli basis: states with high values of $M_{2}$ or $M_{\operatorname{lin}}$ require an exponential amount of classical resources to be simulated, and hence are those which may exhibit quantum advantage.

Moreover, a result shown in [19] establishes a bound between SE and minimum number of copies needed to achieve a certain fidelity in certification protocols: therefore, from a computational perspective, the knowledge of SE of quantum gravity states is an essential tool to optimize time and resources involved in the preparation of the states on a quantum computer, once a desired value of fidelity has been established.

III Quantum tetrahedron

Spin networks are symmetric tensor network states defined by a graphs $\Gamma$, labeled by SU$(2)$ irreducible representations and intertwining operators. In loop quantum gravity, such states encode the quantum description of the 3D space manifold into purely combinatorial and algebraic variables [49, 50]. More generally, spin networks can be defined as abstract quantum many-body-like collections of fundamental quanta of space, connected by maximally entangled states to describe quantized discrete spatial geometries [39, 51, 52].

Consider a given graph $\Gamma$. To each edge $e$ of $\Gamma$ we associate a half-integer spin variable $j_{e}$ labeling a $(2j_{e}+1)$-dimensional $SU(2)$ irreducible representation space $\mathcal{H}_{{j_{e}}}$. At the same time, each $N$-valent node $n$ carries an intertwiner state $\ket{I_{n}}$ in the $SU(2)$-invariant Hilbert space $\mathcal{H}_{I}=\text{Inv}_{\text{SU}(2)}\left[\bigotimes_{i=1}^{N}\mathcal{H}_{{j_{i}}}\right]$, which is the degeneracy space associated to the recoupling of the $N$ spins meeting at the node into a singlet (gauge invariant) representation. A spin-network basis state is the triple $\ket{\Gamma;\{j_{e}\};\{I_{n}\}}$, defined by the direct sum over $j_{e}$ of the tensor product of the gauge invariant states $\ket{I_{n}}$ at all nodes:

$$\ket{\Gamma;\{j_{e}\};\{I_{n}\}}:=\bigoplus_{\{j_{e}\}}\bigotimes_{n}\ket{I_{n}}\,.$$ (5)

Spin network states can be enriched with a geometric interpretation: each $N$-valent node is dual to a $(N-1)$-simplex, represented at the quantum level by the intertwiner state. Accordingly, a 4-valent intertwiner state $\ket{I}$ describes the quantized geometry of a 3-simplex, namely a quantum tetrahedron [53, 54]. In the context of LQG, the presence of clearly defined geometric operators, such as the area operator and the volume operator [49, 43], further strengthens the geometric interpretation; in particular, the volume operator acting on a $4-$valent node is given by

$$\hat{V}=\left(\frac{\sqrt{2}}{3}\right)^{2}(8\pi\gamma)^{3}\left(-i\left[\vec{J}_{1}\cdot\vec{J}_{2},\vec{J}_{1}\cdot\vec{J}_{3}\right]\right)\,,$$ (6)

and has a diagonal representation on the intertwiner basis [55]. Note that this operator is Hermitian but not positive, since it also keeps information of the space orientation [36]. The two possible signs split the degeneracy between the eigenvalues.

For the sake of our work, we focus on the case of a single 4-valent intertwiner with all spins fixed to the same value $j$. The Hilbert space of the tensor product of the four spins $j$ recoupling into the total spin $J$ can be written as

$$\mathcal{H}_{{j}}^{\otimes 4}=\bigoplus_{J=0}^{4j}D^{J}_{j}\mathcal{H}_{{J}}\,,$$ (7)

where the multiplicity spaces (or degeneracy spaces) $D^{J}_{j}$ consist in the spaces of $SU(2)$-invariant intetrtwiner states in the tensor product of the total spin Hilbert space $\mathcal{H}_{{J}}$ with the individual spins $\mathcal{H}_{{j}}^{\otimes 4}$.

A state $\ket{I}\in\text{Inv}_{\text{SU}(2)}[\mathcal{H}_{{j}}^{\otimes 4}]$ can be written as the recoupling of four spins $\ket{\vec{m}}=\bigotimes_{i=1}^{4}\ket{m_{i}}\in\mathcal{H}_{{j}}^{\otimes 4}$ into the singlet ($J=0$):

$$\ket{I}=N\sum_{K=0}^{2j}\sum_{M=-K}^{K}\sum_{\{\vec{m}\}}C^{K,M}_{jm_{1}jm_{2}}C^{K,-M}_{jm_{3}jm_{4}}\ket{\vec{m}}\,,$$ (8)

where $N$ is a normalization factor and $C_{jm_{i}jm_{k}}^{KM}$ are the Clebsch-Gordan coefficients involved in the intermediate recoupling of the two pairs of spins $j$ into two states with spin $K$; to ensure the final recoupling into a singlet state, the magnetic indices $M$ have opposite signs. We use the shorthand notation $\sum_{\{\vec{m}\}}$ to indicate the sum running over the $(2j+1)^{4}$ basis element in $\mathcal{H}_{j}^{\otimes 4}$. In the following, we set all spins to the value $j=1/2$ (see Appendix C for a short discussion on the case $j\not=1/2$). In this case, the spin representation space reduces to $\mathcal{H}_{{\frac{1}{2}}}=\text{span}\left\{\ket{\uparrow}\equiv\ket{0},\ket{\downarrow}\equiv\ket{1}\right\}\simeq\mathbb{C}^{2}$, therefore, each spin, which is dual to one of the faces of the tetrahedron, is described by a qubit. The Hilbert space of the tensor product of the four $1/2$-spins decomposes as

$$\mathcal{H}_{{\frac{1}{2}}}^{\otimes 4}=\bigoplus_{J}D^{J}_{\frac{1}{2}}\mathcal{H}_{{J}}=2\mathcal{H}_{{0}}\oplus 3\mathcal{H}_{{1}}\oplus\mathcal{H}_{{2}}\,.$$ (9)

The $SU(2)$-invariant subspace $\mathcal{H}_{{0}}$ comes with degeneracy $D^{0}_{1/2}=2$, hence any gauge invariant state $\ket{I}$ is a one qubit state in the 2-dimensional intertwiner Hilbert space

$$\mathcal{H}_{I}:=\text{Inv}_{\text{SU}(2)}[\mathcal{H}_{{1/2}}^{\otimes 4}]=\mathcal{H}_{{0}}\oplus\mathcal{H}_{{0}}\simeq\mathbb{C}^{2}\,,$$ (10)

where we can choose a suitable basis $\{\ket{0_{s}},\,\ket{1_{s}}\}$. We can represent a 4-valent intertwiner state $\ket{I}\in\mathcal{H}_{I}$ both in the computational basis $\{\ket{0},\ket{1}\}^{\otimes 4}$ of the 4-qubits space $\mathcal{H}_{{1/2}}^{\otimes 4}$ and in the logical basis $\{\ket{0_{s}},\ket{1_{s}}\}\in\mathcal{H}_{{0}}\oplus\mathcal{H}_{{0}}$. In the following we will refer to a state $\ket{I}$ written in the logical basis as a logical intertwiner qubit (LIQ) state [56].

In terms of the computational basis, the expression of the elements of logical basis can be found using (8). Explicitly, one finds [57, 11]:

$$\ket{0_{s}} =\frac{1}{2}\left(\ket{0101}+\ket{1010}-\ket{0110}-\ket{1001}\right)\,,$$

$$\ket{1_{s}} =\frac{1}{\sqrt{3}}\big{[}\ket{0011}+\ket{1100}-\frac{1}{2}\big{(}\ket{0101}$$

$$+\ket{1010}+\ket{0110}+\ket{1001}\big{)}\big{]}\,.$$ (11)

A generic LIQ state is given by the following Bloch representation

$$\ket{I(\theta,\phi)}=\cos{\frac{\theta}{2}}\ket{0_{s}}+\sin{\frac{\theta}{2}}e^{i\phi}\ket{1_{s}}\,,$$ (12)

where $\theta\in[0,\pi]$ and $\phi\in[0,2\pi)$ are angles on the Bloch sphere. Finally, it can be shown [11] that the eigenstates of the volume operator take the following form:

$$\ket{V_{+}} =\frac{1}{\sqrt{2}}(\ket{0_{s}}-i\ket{1_{s}})\,,$$ (13)

$$\ket{V_{-}} =\frac{1}{\sqrt{2}}(\ket{0_{s}}+i\ket{1_{s}})\,,$$ (14)

that means $\hat{V}\ket{V_{\pm}}=\pm V_{0}\ket{V_{\pm}}$. These states represent remarkable points on the Bloch sphere, as they are placed at the equator and their angular coordinates are $(\theta=\frac{\pi}{2},\phi=\frac{\pi}{2})$ and $(\theta=\frac{\pi}{2},\phi=\frac{3\pi}{2})$.

IV Stabilizer entropy of a quantum tetrahedron

We shall investigate whether the basis states of the gauge invariant subspace $\mathcal{H}_{I}\subset\mathcal{H}_{{1/2}}^{\otimes 4}$ possess SE. In this sense, we investigate whether the gates of the quantum circuit associated with the construction of LIQ states $\ket{I}$ from the computational basis belong to the Clifford group. Without loss of generality, we can start from the reference state $\ket{0}^{\otimes 4}$, and then look for unitary transformations such that $\ket{0_{s}}=\operatorname{U}_{0_{s}}\ket{0}^{\otimes 4}$ and $\ket{1_{s}}=\operatorname{U}_{1_{s}}\ket{0}^{\otimes 4}$. Using the relations in Eq. 11, we can express the unitary operators $\operatorname{U}_{0_{s}}$ and $\operatorname{U}_{1_{s}}$ in terms of a set of unitary gates acting on a stabilizer reference state in the 4-qubit Hilbert space. Thereby, it is possible to realize the generic intertwiner state $\ket{{I(\theta,\phi)}}$ via the quantum circuit given in Fig. 1 [11] (see also [58] and [59, 60] for different descriptions of quantum spin network circuits). The only non-Clifford gates in the circuit are the two special unitary operators, $U$ and $V$, which depend on the parameters $\theta$ and $\phi$ of the output state $\ket{{I(\theta,\phi)}}$.

Figure 1: Circuit realization of the intertwiner state $\ket{I(\theta,\phi)}$. $H$ is the Hadamard gate, $U$ and $V$ are specific gates depending on $\theta$ and $\phi$, and the others are CNOT and anti-CNOT gates.

Such operators can be written as $2\times 2$ matrices [11]:

$$U=\begin{pmatrix}c_{0}&\sqrt{\left|{c_{+}}\right|^{2}+\left|{c_{-}}\right|^{2}}\\ -\sqrt{\left|{c_{+}}\right|^{2}+\left|{c_{-}}\right|^{2}}&c_{0}^{*}\end{pmatrix}\,,$$ (15)

$$V=\frac{1}{\sqrt{\left|{c_{+}}\right|^{2}+\left|{c_{-}}\right|^{2}}}\begin{pmatrix}-c_{+}&c_{-}^{*}\\ -c_{-}&-c_{+}^{*}\end{pmatrix},$$ (16)

where the coefficients $c_{0}$ and $c_{\pm}$ are given by the following functions of $\theta$ and $\phi$:

$$c_{0}=\sqrt{\frac{2}{3}}e^{i\phi}\sin{\frac{\theta}{2}}\,,$$ (17)

$$c_{\pm}=\frac{1}{\sqrt{2}}\left[-\frac{1}{\sqrt{3}}e^{i\phi}\sin{\frac{\theta}{2}}\pm\cos{\frac{\theta}{2}}\right]\,.$$ (18)

Our first remark is that the SE of an intertwiner state $\ket{I(\theta,\phi)}$ seen from the full space is entirely determined by the gates $U$ and $V$, which in general are not Clifford gates. Indeed, the calculation of the SE for the LIQ basis states (see Appendix A for full details) gives

$$M_{2}(\ket{0_{s}})=0\,,$$ (19)

$$M_{2}(\ket{1_{s}})=0,847997\,.$$

Hence, a generic superposition of these two basis states will not be in general a stabilizer state. We can plot the SE of all the intertwiner states in the basis of the Hilbert space $\mathcal{H}_{{1/2}}^{\otimes 4}$ as a function of the Bloch sphere representation of the gauge invariant space $\mathcal{H}_{I}$ (Fig. 2).

Figure 2: Contour plot of the $2-$Stabilizer Rényi entropy $M_{2}(\ket{I(\theta,\phi)})$ as a function of the Bloch sphere angles $\theta\in[0,\pi]$ and $\phi\in[0,2\pi)$. The two volume eigenstates, $\ket{V_{+}}$ and $\ket{V_{-}}$, are highlighted in the plot.

In particular, we remark the behavior of the volume eigenstates $\ket{V_{\pm}}$: from the SE standpoint, they are located at the center of the regions of maximal non-stabilizerness and they belong to the same Clifford orbit, having exactly the same value of SE:

$$M_{2}(\ket{V_{+/-}})=1,16993\,.$$ (20)

We recall that the SE of these states is calculated with respect to the full Pauli operator basis of $\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$, namely the Pauli group of $4$-qubits $\mathbb{P}_{4}$. This distinction is necessary since, as we have seen in Eq. (19), the basis states themselves are not stabilizer states in $\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$, whereas they are stabilizers in $\mathcal{H}_{I}$. Finally, notice that these states are also non separable in the tensor product structure associated to $\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$, as a demonstration that both entanglement and SE are necessary for the gauge structure[61, 62, 63, 64].

V Average SE of $SU(2)$-gauge invariant subspace

The reason for having different values of SE in the two bases states, as shown in Eq. (11), is rooted in the gauge structure of the intertwiner state. We shall explore and generalize this result further. The $SU(2)-$gauge invariant intertwiner space $\mathcal{H}_{I}$ is a subspace of the Hilbert space of four qubits $\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$. However, the SE in a particular basis state has hardly any physical meaning as any superposition of states in $\mathcal{H}_{I}$ is allowed and SE is not constant in any given subspace. To associate SE in a meaningful way to a subspace, one is also confronted with the choice of the Pauli basis with respect to the SE must be computed. To be concrete, if one has a state $\ket{\psi}\in\mathcal{H}_{I}$, do we want to know the SE of this state as a state expressed in the computational basis of $\mathcal{H}^{\otimes 4}$ or in the computational basis of $\mathcal{H}_{I}$ given by $\{\ket{0_{s}},\ket{1_{s}}\}$?

In the construction of a quantum gauge theory, one starts with an ambient Hilbert space $\mathcal{H}_{\operatorname{tot}}$. The gauge constraints are expressed as local projectors $\Pi^{(s)}$ and the gauge invariant subspace is the global projection over all the local gauge constraints,

$$\mathcal{H}_{G}=\Pi_{G}\mathcal{H}_{\operatorname{tot}}\,,$$ (21)

with $\Pi_{G}:=\prod_{s}\Pi^{(s)}$.

We associate to a subspace $\mathcal{H}_{G}$ its average SE with respect to the Heisenberg-Weyl basis $\{\mathcal{D}^{(i)}\}\subset\mathcal{L}(\mathcal{H}_{{i}})$. In this notation, $i=0$ refers to the Heisenberg-Weyl basis in $\mathcal{H}_{\operatorname{tot}}$, while $i=G$ refers to the gauge invariant Hilbert space $\mathcal{H}_{G}$. With this notation, the linear SE is $M^{(i)}_{\operatorname{lin}}(\psi):=1-d_{i}\operatorname{tr}[Q^{(i)}\psi^{\otimes 4}]\,,i\in\{0,G\}$ and the operator

$$Q^{(i)}:=d_{i}^{-2}\sum_{(p,q)}D_{(p,q)}^{(i)\otimes 4}$$ (22)

is defined accordingly with the corresponding Heisenberg-Weyl (see Appendix C for details) basis. The average SE in the subspace $\mathcal{H}_{G}$ is then defined as

$$M^{i} := \mathbb{E}_{U_{G}}M_{\operatorname{lin}}^{(i)}(\psi_{U_{G}})$$ (23)

$$= 1-d_{i}\operatorname{tr}[Q^{(i)}\mathbb{E}_{U}\psi_{U_{G}}^{\otimes 4}]\,,$$ (24)

with $\mathbb{E}_{U_{G}}$ denoting the unitary group average with respect to the Haar measure over $\mathcal{H}_{G}$, and $\psi_{U_{G}}:=U_{G}\psi U_{G}^{\dagger}$. In order to perform the Haar average over the subspace $\mathcal{H}_{G}$ we need the following:

Lemma 1.

Given any Hilbert space $\mathcal{H}=\mathcal{H}_{R}\oplus\mathcal{H}_{R}^{\perp}$, with $\operatorname{dim}(\mathcal{H})=d$ and $\operatorname{dim}(\mathcal{H}_{R})=d_{R}$, the unitary Haar average of $k$ copies of the state over the subspace $\mathcal{H}_{R}$ is given by

$$\mathbb{E}_{U_{R}}\psi_{U_{R}}^{\otimes k}=c_{R}(d,d_{R})\Pi_{R}^{\otimes k}\mathbb{E}_{U}\psi_{U}^{\otimes k}\,,$$ (25)

with $U_{R}\in\mathcal{U}(\mathcal{H}_{R})$, $U\in\mathcal{U}(\mathcal{H})$, $\Pi_{R}$ the projector onto $\mathcal{H}_{R}$, $\psi\in\mathcal{H}_{R}$ and $c_{R}(d,d_{R})=\binom{d+k-1}{k}\binom{d_{R}+k-1}{k}^{-1}$.

A proof of this lemma is given in Section B.2. Based on the above lemma, we obtain the general result

$$M^{i}=1-d_{i}c_{i}(d,d_{i})\operatorname{tr}[Q^{(i)}\Pi_{G}^{\otimes 4}\mathbb{E}_{U}\psi_{U}^{\otimes 4}]\,,$$ (26)

where now $\mathbb{E}_{U}$ denotes the Haar average over the full unitary group on $\mathcal{H}_{\operatorname{tot}}$.

Eq.(26) shows explicitly how the gauge structure enters the SE through the projector $\Pi_{G}$. When $i=G$, this projector becomes the identity map and the average SE has no recollection of the gauge structure. In order to quantify the amount of SE due to this structure, we define the SE-gap as

$$\Delta M(\mathcal{H}_{G}) := M^{0}-M^{G}$$

$$= \mathbb{E}_{U_{G}}M_{\operatorname{lin}}^{(0)}(\psi_{U_{G}})-\mathbb{E}_{U_{G}}M_{\operatorname{lin}}^{(G)}(\psi_{U_{G}})\,.$$

For a general $SU(2)-$gauge structure, namely the spin $j$ intertwiner states, the projector $\Pi_{I}$ reads (see Appendix C)

$$\Pi_{I} = \ket{0_{s}}\bra{0_{s}}+\ket{1_{s}}\bra{1_{s}}+\dots+\ket{2j_{s}}\bra{2j_{s}}$$ (28)

Let us now specialize these formulae to the case of $j=1/2$, that is, the quantum tetrahedron. Using Lemma 1, with $k=4$, $\mathcal{H}=\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$, $\mathcal{H}_{R}=\mathcal{H}_{I}$ and $\Pi_{I}=\ket{0_{s}}\bra{0_{s}}+\ket{1_{s}}\bra{1_{s}}$ being the projector onto the intertwiner space, (see Appendix B for details) we find

$$M^{0}=17/45.$$ (29)

Now, notice that the logical states $\ket{0_{s}},\ket{1_{s}}$ in the Pauli basis in which $Z_{I}=\ket{0_{s}}\bra{0_{s}}-\ket{1_{s}}\bra{1_{s}}$ are obviously stabilizer states with zero stabilizer entropy. Nevertheless even according to this Pauli basis there will be a non-zero value for the average stabilizer entropy. This space is just a generic qubit $\mathbb{C}^{2}$ from the point of view of the stabilizer entropy. The average SE in this space is thus just the average stabilizer entropy of a qubit, namely $M^{I}:=\mathbb{E}_{U_{I}}M_{\operatorname{lin}}^{(I)}(\psi_{U_{I}})=1/5$, as calculated in [16] (we also explicitly show this calculation in Appendix B).

Putting the pieces together, we are able to calculate the average SE gap, which reads

$$\Delta M(\mathcal{H}_{I})=8/45\,.$$ (30)

A value of $\Delta M(\mathcal{H}_{I})$ greater than zero tells us that projecting a generic $4$-qubit state onto this gauge invariant subspace has a cost in terms of non-stabilizer resources.

In particular, this means that the gauge structure bears a cost in terms of simulability, which is very important as one scales the system to many nodes.

VI Simulations of quantum gravity states

Very recently, quantum gravity states have been physically implemented on quantum computers [11, 12], and in particular, quantum tetrahedra states. The very first layer of difficulty that must be faced in a laboratory when attempting to conduct a quantum experiment (including one regarding simulations of quantum gravity states) is the preparation of an initial state that is faithful to the theoretical one $\ket{\psi}$. Typically, a large initial sample must be prepared, resulting in a mixed output $\tilde{\psi}$ from the processor. Ensuring the correct functioning of quantum devices in terms of the accuracy of the output requires a certification protocol [65]; one of the possible measures of the quality of the realization $\tilde{\psi}$ of a state $\ket{\psi}$ is the fidelity $\mathcal{F}(\ket{\psi},\tilde{\psi})=\operatorname{tr}(\psi\tilde{\psi})$ that measures the precision of preparation. It is known that SE can provide useful indications in an experimental setting.

In this last section, we argue that the numerical results found in (19) and (20) have a direct use in the recent results on quantum gravity states simulations. Indeed, we can use the SE of these states to estimate the maximum fidelity achievable with a given number of preparations or, conversely, the minimum number of preparations needed to achieve a desired value of fidelity within a desired error [19].

In [11], the authors present a realization of the intertwiner states $\ket{0_{s}},\ket{1_{s}}$ as of Eq. (11) as well as the volume eigenstates $\ket{V_{+}},\ket{V_{-}}$ on a $5$-qubit (Yorktown) and a $15$-qubit (Melbourne) IBM superconducting quantum computer. By performing 10 rounds of 1024 quantum measurements, they obtained fidelity values $\mathcal{F}_{exp}$ of the prepared states with respect to the theoretical ones. As anticipated in Sec. II, the explicit lower bound on the number $N_{\tilde{\psi}}^{\text{min}}$ of preparations needed to achieve an accuracy $\epsilon$, with a probability of failure $\delta$, is related to the SE of the theoretical state $\ket{\psi}$ as follows:

$$N_{\tilde{\psi}}\geq\frac{2}{\epsilon^{2}}\ln\left(\frac{2}{\delta}\right)\exp[M_{2}(\ket{\psi})]\,.$$ (31)

Inverting this relation, one can get an upper bound on the maximum achievable fidelity with a given number of preparations $N$ with failure probability $\delta$:

$$\mathcal{F}_{\text{max}}\leq 1-\sqrt{\frac{2}{N}\ln\left(\frac{2}{\delta}\right)}\exp\left[\frac{M_{2}(\psi)}{2}\right]\,.$$ (32)

$\ket{I}$	$\mathcal{F}_{\text{exp}}$	$\sigma$	$M_{2}$	$(1-\delta)$	$\mathcal{F}_{\text{max}}$	$N_{\tilde{\psi}}^{\text{min}}$
$\ket{0_{s}}$	0,906	0,005	0	0,05	0,973	835
$\ket{1_{s}}$	0,916	0,007	0,847997	0,05	0,959	2441
$\ket{V_{+}}$	0,918	0,009	1,16993	0,05	0,952	3535
$\ket{V_{-}}$	0,917	0,008	1,16993	0,05	0,952	3450

Table 1: Comparison between the experimentally obtained results of the fidelity in [11] and the bounds given by Eqs. (31) and (32). $\mathcal{F}_{\operatorname{max}}$ is calculated using $N=10240$ and the probability of failure $\delta$ is inferred using the Bienaymé-Čebyšëv inequality assuming the statistical error $\sigma$ as the variance of the probability distribution.

As one can see from the Table 1, the experimentally obtained data for the the fidelity is perfectly compatible with the bounds provided by the SE. However, the authors used many more copies than the minimum shown in the table. However, it is important to note that the constraints provided by Eqs. (31) and (32) are purely theoretical, hence independent of the inherent noise sources peculiar to the specific implementation protocols and hardware used in the experimental setting at hand. Nevertheless, by considering the SE of the states one wishes to prepare, the hardware resources and the number of preparations can be managed more efficiently in future experiments of quantum gravity states.

VII Discussion

In this paper, we have shown that the gauge invariant structure of quantum geometry states has a cost in terms of non-stabilizer resources. This implies that simulations of quantum geometry states can run more efficiently on a quantum computer and that preparations of future experiments can be more efficient if the non-stabilizer property of the state is taken into account. Moreover, we have seen that eigenstates of the oriented volume have near-maximal amount of SE: this begs the question of why such states possess greater quantum complexity, suggesting that a correspondence between entanglement and geometry in quantum gravity may extend at a deeper layer of quantumness. Concretely, the first step to answer this question will be to repeat this analysis for a generic spin-$j$ intertwiner and see if, also in that setting, the volume eigenstates are states with maximum SE. A further intriguing direction to explore is the role of non-stabilizer resources when taking into account the quantum state associated to an actual spin network, that is a collection of intertwiners describing a quantum simplicial complex: in that general setting, one should expect additional non-stabilizer resources coming from the graph structure, that is the adjacency matrix, describing the connectivity and the non-trivial additional geometrical degrees of freedom described by the holonomies dressing the links. In particular, in this sense, we expect that non-stabilizerness can be further used to characterise the transition amplitudes, that is the evolution, of quantum geometry states (see e.g. [12]).

In general terms, the dependence of the SE gap from the projector onto the intertwiner subspace opens to more wide-reaching questions: how does the gauge structure affect non-stabilizerness in a more general setting? Can we use this formalism to characterize the SE of other quantum gauge theories? Does abelianity (or lack thereof) of the gauge group play a role in the non-stabilizer resources of the gauge invariant subspace?

VIII acknowledgments

The authors thank L. Leone and L. Vacchiano for important discussions. AH acknowledges support from the PNRR MUR project PE0000023-NQSTI and PNRR MUR project CN $00000013$-ICSC.

Appendix A Calculation of the SE of the intertwiner basis states

In this section, we show the details of the calculations of the SE of the intertwiner basis state, namely the south and north poles of the Bloch sphere (i.e. $\theta=0$ and $\theta=\pi$). Recalling the circuit realization of the intertwiner states in (Fig. 1), we focus on the unitary operators $U_{0}$, $V_{0}$, $U_{1}$, $V_{1}$ involved in the preparation of basis states $\ket{0_{s}}$ and $\ket{1_{s}}$, which can be calculated by inserting in Eq. 17 $\phi=0$ and $\theta=0,\pi\,$, respectively:

$$U_{0}=\left(\begin{array}[]{cc}0&1\\ -1&0\end{array}\right)\,\,\, V_{0}=\frac{1}{\sqrt{2}}\left(\begin{array}[]{cc}-1&-1\\ 1&-1\end{array}\right)$$ (37)

$$U_{1}=\sqrt{\frac{2}{3}}\left(\begin{array}[]{cc}1&\frac{1}{\sqrt{2}}\\ -\frac{1}{\sqrt{2}}&1\end{array}\right)\,\,\, V_{1}=\frac{1}{\sqrt{2}}\left(\begin{array}[]{cc}1&-1\\ 1&1\end{array}\right)$$ (42)

Notice that, with the exception of the unitaries that we just calculated, all the gates involved in the circuit in (Fig. 1) are Clifford gates. Hence, to estimate the magic of an intertwiner state, we can focus our analysis on the reduced 2-qubit system given by the action of $CV_{0}\,(U_{0}\otimes\mathbbm{1})$ on $\ket{0}^{\otimes 2}$, where $CV_{0}$ is the controlled-$V_{0}$ gate. The reduced circuit is represented in Fig. 3.

Figure 3: 2-qubit circuit with the non-trivial contribution to the magic of an Intertwiner state.

This circuit returns the realization of $\ket{0_{s}}$ as a 2-qubit state:

$$\ket{0_{s}}_{2}=\frac{1}{\sqrt{2}}\left(\ket{10}-\ket{11}\right)$$ (43)

In this way, we rule out all trivial contribution to the non-stabilizerness, isolating only the significant ones.

For the particular case of $\ket{0_{s}}$ we can prove that the magic produced by the operators $U_{0}$ and $V_{0}$ is 0; i.e. $\ket{0_{s}}$ is a stabilizer state. Let us write the matrix form of the operators:

$$U_{0}\otimes\mathbbm{1}=\left(\begin{array}[]{cccc}0&0&1&0\\ 0&0&0&1\\ -1&0&0&0\\ 0&-1&0&0\end{array}\right)$$ (44)

$$CV_{0}=\left(\begin{array}[]{cccc}1&0&0&0\\ 0&1&0&0\\ 0&0&-\frac{1}{\sqrt{2}}&-\frac{1}{\sqrt{2}}\\ 0&0&\frac{1}{\sqrt{2}}&-\frac{1}{\sqrt{2}}\end{array}\right)$$ (45)

We consider the state $\psi_{0_{s}}=\ket{0_{s}}_{2}\bra{0_{s}}_{2}$, which in the Pauli basis reads

$$\psi_{0_{s}} =\frac{1}{d}\sum_{P\in\mathbbm{P}_{2}}\operatorname{Tr}(\psi_{0_{s}}P)P$$

$$=\frac{1}{4}\left(\mathbbm{1}\otimes\mathbbm{1}-\mathbbm{1}\otimes\operatorname{X}-\operatorname{Z}\otimes\mathbbm{1}+\operatorname{Z}\otimes\operatorname{X}\right)$$ (46)

The magic of the basis state is

$$M_{2}(\psi_{0_{s}})=-\log d^{-1}\sum_{P\in\mathbb{P}_{2}}\operatorname{tr}^{4}(\psi_{0_{s}}P).$$ (47)

Since the trace of the product of Pauli matrices is equal to $d$ only if the product returns $\mathbbm{1}\otimes\mathbbm{1}$ and $0$ otherwise, among the terms of the sum in (47) there are only four non vanishing contributions, which are the ones with $P$ equal to one of the terms in the Pauli decomposition of the state. Equation (47) returns

$$M_{2}(\psi_{0_{s}})=-\log\left(d^{-1}\left(\frac{d}{4}+\frac{d}{4}+\frac{d}{4}+\frac{d}{4}\right)\right)=0$$ (48)

Namely, the intertwiner state $\ket{0_{s}}$ is a stabilizer state.

We now repeat the same procedure for the state $\ket{1_{s}}$. First, we realize the operators as matrix:

$$U_{1}\otimes\mathbbm{1}=\left(\begin{array}[]{cccc}\sqrt{\frac{2}{3}}&0&\frac{1}{\sqrt{3}}&0\\ 0&\sqrt{\frac{2}{3}}&0&\frac{1}{\sqrt{3}}\\ -\frac{1}{\sqrt{3}}&0&\sqrt{\frac{2}{3}}&0\\ 0&-\frac{1}{\sqrt{3}}&0&\sqrt{\frac{2}{3}}\end{array}\right)$$ (49)

$$CV_{1}=\left(\begin{array}[]{cccc}1&0&0&0\\ 0&1&0&0\\ 0&0&\frac{1}{\sqrt{2}}&-\frac{1}{\sqrt{2}}\\ 0&0&\frac{1}{\sqrt{2}}&\frac{1}{\sqrt{2}}\end{array}\right)$$ (50)

The action of this operators on $\ket{0}^{\otimes 2}$ returns

$$\ket{1_{s}}_{2}=\sqrt{\frac{2}{3}}\ket{00}-\frac{1}{\sqrt{6}}\ket{10}-\frac{1}{\sqrt{6}}\ket{11}$$ (51)

We write the state $\psi_{1_{s}}$ as

$$\psi_{1_{s}} =\frac{1}{4}\Big{(}\mathbbm{1}\otimes\mathbbm{1}+\frac{1}{3}\mathbbm{1}\otimes\operatorname{X}+\frac{2}{3}\mathbbm{1}\otimes\operatorname{Z}-\frac{2}{3}\operatorname{X}\otimes\mathbbm{1}$$

$$-\frac{2}{3}\operatorname{X}\otimes\operatorname{X}-\frac{2}{3}\operatorname{X}\otimes\operatorname{Z}+\frac{2}{3}\operatorname{Y}\otimes\operatorname{Y}+\frac{1}{3}\operatorname{Z}\otimes\mathbbm{1}$$

$$-\frac{1}{3}\operatorname{Z}\otimes\operatorname{X}+\frac{2}{3}\operatorname{Z}\otimes\operatorname{Z}\Big{)}$$ (52)

There are ten non-vanishing contributions to the magic of this state, each of which is equal to the fourth power of one of the coefficients of (52). Direct calculation returns

$$M_{2}(\psi_{1_{s}})=0,847997$$ (53)

Appendix B Average SE gap

In this section, we show the details of the calculations of the average SE gap shown in Section V,

$$\Delta M(\mathcal{H}_{I}):=\mathbb{E}_{U}M_{\operatorname{lin}}^{(0)}(\psi_{U})-\mathbb{E}_{U}M_{\operatorname{lin}}^{(I)}(\psi_{U})$$ (54)

with $\mathbb{E}_{U}$ denoting the unitary group average with respect to the Haar measure, $M_{\operatorname{lin}}^{(i)}(\psi):=1-d_{i}\operatorname{tr}Q^{(i)}\psi^{\otimes 4}\,,i=\{0,I\}$ being the linear SE, $d_{i}$ being the dimensions of the Hilbert spaces involved (recall that $d_{I}$ is the dimension of the gauge invariant intertwiner space, hence $d_{I}=2$, whereas $d_{0}$ is the dimension of the ambient $4$-qubit space, thus $d_{0}=16$), and

$$\begin{split}Q^{(I)}&=\frac{1}{d_{I}^{2}}\sum_{P\in\mathbb{P}_{1}}P^{\otimes 4}\\ Q^{(0)}&=\frac{1}{d_{0}^{2}}\sum_{P\in\mathbb{P}_{4}}P^{\otimes 4}\end{split}$$ (55)

with $\mathbb{P}_{n}$ being the Pauli group on $n$ qubits modulo the phases.

B.1 Haar averages and calculation of $M^{I}$

We start by calculating $\mathbb{E}_{U}M_{\operatorname{lin}}^{(I)}(\psi_{U})$ since it is the simplest one: plugging the definition of $M_{\operatorname{lin}}$ in the Haar average reads

$$\mathbb{E}_{U}M_{\operatorname{lin}}^{(I)}(\psi_{U})=\mathbb{E}_{U}1-d_{I}\operatorname{tr}Q^{(I)}\psi_{U}^{\otimes 4}=1-d_{I}\operatorname{tr}Q^{(I)}\mathbb{E}_{U}U^{\otimes 4}\psi^{\otimes 4}U^{{\dagger}\otimes 4}$$ (56)

since the Haar average is linear. Our focus, then, is on evaluating

$$\operatorname{SP}^{(I)}:=\mathbb{E}_{U}\operatorname{tr}Q^{(I)}\psi_{U}^{\otimes 4}$$ (57)

In general, carrying out the Haar average of operators of the form $U^{\otimes 4}XU^{{\dagger}\otimes 4}\,,X\in L(\mathcal{H}^{\otimes k})$ requires the knowledge of the commutant of the $k$-tensored representation of the unitary group, according to Schur’s Lemma. By the Schur-Weyl duality, the basis of the full commutant of $U^{\otimes k}$ is constituted by the permutation operators acting over the $k$ copies of the Hilbert space of interest. In particular, the average of the $k$ copies of a state is carried out in detail in [66], and reads

$$\mathbb{E}_{U}U^{\otimes k}\psi^{\otimes k}U^{{\dagger}\otimes k}=\binom{d+k-1}{k}^{-1}\Pi_{\text{sym}}^{(k)}\,,$$ (58)

with

$$\Pi_{\text{sym}}^{(k)}:=\frac{1}{k!}\sum_{\pi\in S_{k}}\operatorname{T}_{\pi}$$ (59)

being the projector onto the subspace of $\mathcal{H}^{\otimes k}$ which is symmetric under permutations of $k$ objects. The result of this average is to be expected by the fact that operators belonging $L(\mathcal{H}^{\otimes k})$ of the form $\psi^{\otimes k}$ are actually symmetric under permutation operators, so the weight associated to each of this operators must be the same and is hence only determined by the normalization.

The permutation operators relative to $\pi\in S_{k}$ can be written in the computational basis of $\mathcal{H}^{\otimes k}$, namely $\{\ket{i_{1}\dots i_{k}}\}_{i_{h}=1}^{d}$ in this way:

$$\operatorname{T}_{\pi}=\sum_{i_{1},\dots,i_{k}}\ket{\pi(i_{1})\dots\pi(i_{k})}\bra{i_{1}\dots i_{k}}\,.$$ (60)

Notice that the permutation operators are invariant under $k$-copies of unitary operators $U^{\otimes k}$: this means that such operators will have the same expression in any basis of $\mathcal{H}^{\otimes k}$. Now, starting In our case, we are interested in the Haar average of four copies of the state, namely

$$\begin{split}\mathbb{E}_{U}\psi_{U}^{\otimes 4}&=\binom{d+3}{4}^{-1}\frac{1}{4!}\sum_{\pi\in S_{4}}\operatorname{T}_{\pi}\\ &=[(d+3)(d+2)(d+1)d]^{-1}\sum_{\pi\in S_{4}}\operatorname{T}_{\pi}\,.\end{split}$$ (61)

Substituting this expression in Eq. (57) we get

$$\begin{split}\mathbb{E}_{U}\operatorname{tr}Q^{(I)}\psi_{U}^{\otimes 4}&=[(d_{I}+3)(d_{I}+2)(d_{I}+1)d_{I}]^{-1}\operatorname{tr}Q^{(I)}\Pi_{\text{sym}}^{(I)}\\ &=[(d_{I}+3)(d_{I}+2)(d_{I}+1)d_{I}]^{-1}\frac{1}{d_{I}^{2}}\sum_{\pi\in S_{4}}\sum_{P\in\mathbb{P}_{2}}\operatorname{tr}P^{\otimes 4}\operatorname{T}_{\pi}^{I}\,.\end{split}$$ (62)

In order to tackle this calculation, one calculates the sums of the Pauli operators permutation by permutation: by means of example, we show the calculation for one permutation, but the treatment is similar for all of them. Let’s take, say, the $3$-cycle $(123)$: the permutation operator associated to this permutation reads

$$\operatorname{T}_{(123)}=\sum_{ijkl}\ket{kijl}\bra{ijkl}\,,$$ (63)

hence

$$\begin{split}\operatorname{tr}Q^{(I)}\operatorname{T}_{(123)}&=\frac{1}{d_{I}^{2}}\sum_{P\in\mathbb{P}_{2}}\sum_{ijkl}\operatorname{tr}P^{\otimes 4}\ket{kijl}\bra{ijkl}\\ &=\frac{1}{d_{I}^{2}}\sum_{P\in\mathbb{P}_{2}}\sum_{ijk}\operatorname{tr}(P\ket{k}\bra{i})\operatorname{tr}(P\ket{i}\bra{j})\operatorname{tr}(P\ket{j}\bra{k})\sum_{l}\operatorname{tr}(P\ket{l}\bra{l})\\ &=\frac{1}{d_{I}^{2}}\sum_{P\in\mathbb{P}_{2}}\sum_{ijk}\operatorname{tr}(P\ket{k}\bra{i})\operatorname{tr}(P\ket{i}\bra{j})\operatorname{tr}(P\ket{j}\bra{k})\operatorname{tr}P\\ &=\frac{1}{d_{I}}\sum_{ijk}\operatorname{tr}(\ket{k}\bra{i})\operatorname{tr}(\ket{i}\bra{j})\operatorname{tr}(\ket{j}\bra{k})\\ &=\frac{1}{d_{I}}\sum_{ijk}\delta_{ki}\delta_{ij}\delta_{jk}=\frac{1}{d_{I}}\sum_{jk}\delta_{jk}=1\,,\end{split}$$ (64)

where we used the fact that $\sum_{l}\ket{l}\bra{l}=\mathbbm{1}$ and $\operatorname{tr}P=d_{I}\delta_{P,\mathbbm{1}}$. Executing similar calculations for the other $23$ elements of $S_{4}$ and summing all the contributions, one gets

$$\mathbb{E}_{U}\operatorname{tr}Q^{(I)}\psi_{U}^{\otimes 4}=\frac{4d_{I}^{2}+12d_{I}+8}{d_{I}(d_{I}+1)(d_{I}+2)(d_{I}+3)}=\frac{2}{5}\,,$$ (65)

hence

$$\mathbb{E}_{U}M_{\operatorname{lin}}^{(I)}(\psi_{U})=1-d_{I}\mathbb{E}_{U}\operatorname{tr}Q^{(I)}\psi_{U}^{\otimes 4}=\frac{1}{5}\,.$$ (66)

B.2 Haar Averages on subspaces: proof of Lemma 1

Let generic $\mathcal{H}=\mathcal{H}_{R}\oplus\mathcal{H}_{R}^{\perp}$ Hilbert space decomposed in two orthogonal subspaces, with $\operatorname{dim}(\mathcal{H})=d$, $\mathcal{H}_{R}=\operatorname{span}\{\ket{1},\dots,\ket{d_{R}}\}$ and $\mathcal{H}_{R}^{\perp}=\operatorname{span}\{\ket{d_{R}+1},\dots,\ket{d}\}$, and the projector $\Pi_{R}$ onto $\mathcal{H}_{R}$ being

$$\Pi_{R}=\sum_{i=1}^{d_{R}}\ket{i}\bra{i}\,.$$ (67)

Applying the formula shown in Eq. (61), the average onto the subspace $\mathcal{H}_{R}$ reads

$$\mathbb{E}_{U_{R}}\psi_{U_{R}}^{\otimes k}=\binom{d_{R}+k-1}{k}^{-1}\frac{1}{k!}\sum_{\pi\in S_{k}}\operatorname{T}_{\pi}^{R}\,,$$ (68)

whereas if we calculate the average onto the full Hilbert space $\mathcal{H}$ we get

$$\mathbb{E}_{U}\psi_{U}^{\otimes k}=\binom{d+k-1}{k}^{-1}\frac{1}{k!}\sum_{\pi\in S_{k}}\operatorname{T}_{\pi}\,.$$ (69)

Now, it suffices to check that the action of $k$ copies of the projector onto the permutation operators representation on the full Hilbert space gives the permutation operators representation on the subspace $\mathcal{H}_{R}$. Applying the expression shown in Eq. (60) for the permutation operators to this case we get

$$\begin{split}\Pi_{R}^{\otimes k}\operatorname{T}_{\pi}&=\sum_{\{i_{1},\dots,i_{k}\}\in\{1,d\}^{\times k}}\Pi_{R}\ket{\pi(i_{1})}\dots\Pi_{R}\ket{\pi(i_{k})}\bra{i_{1}\dots i_{k}}\,.\end{split}$$ (70)

One then notices that

$$\Pi_{R}\ket{i}=\begin{cases}\ket{i}\quad\text{if}\;i\in\{1,d_{R}\}\\ 0\quad\text{otherwise}\end{cases}\,,$$ (71)

hence Eq. (70) reduces to

$$\Pi_{R}^{\otimes k}\operatorname{T}_{\pi}=\sum_{\{i_{1},\dots,i_{k}\}\in\{1,d\}^{\times k}}\ket{\pi(i_{1}\dots i_{k})}\bra{i_{1}\dots i_{k}}=\operatorname{T}_{\pi}^{R}\,.$$ (72)

Using this result, Eq. (68) reads

$$\begin{split}\mathbb{E}_{U_{R}}\psi_{U_{R}}^{\otimes k}&=\binom{d+k-1}{k}^{-1}\binom{d+k-1}{k}^{-1}\frac{1}{k!}\sum_{\pi\in S_{k}}\operatorname{T}_{\pi}^{R}\\ &=\binom{d_{R}+k-1}{k}^{-1}\frac{1}{k!}\Pi_{R}^{\otimes k}\sum_{\pi\in S_{k}}\operatorname{T}_{\pi}\\ &=\binom{d+k-1}{k}\binom{d_{R}+k-1}{k}^{-1}\mathbb{E}_{U}\psi_{U}^{\otimes k}\\ &\equiv c_{R}(d,d_{R})\mathbb{E}_{U}\psi_{U}^{\otimes k}\,,\end{split}$$ (73)

which completes the proof.∎

B.3 Calculation of $M^{0}$

In this subsection, we use the result of Eq. (73) to calculate $\mathbb{E}_{U}M_{\operatorname{lin}}^{(0)}(\psi_{U})$, which reads

$$\mathbb{E}_{U}M_{\operatorname{lin}}^{(0)}(\psi_{U})=1-d_{0}\mathbb{E}_{U_{I}\in\mathcal{U}(\mathcal{H}_{I})}\operatorname{tr}Q^{(0)}U_{I}^{\otimes 4}\psi_{I}^{\otimes 4}U_{I}^{{\dagger}\otimes 4}\,,$$ (74)

and in particular, we focus on the evaluation of

$$\operatorname{SP}^{(0)}:=\operatorname{tr}Q^{(0)}\mathbb{E}_{U_{I}\in\mathcal{U}(\mathcal{H}_{I})}U_{I}^{\otimes 4}\psi_{I}^{\otimes 4}U_{I}^{{\dagger}\otimes 4}\,,$$ (75)

Using the formula in Eq. (73), Eq. (75) reads

$$\operatorname{SP}^{(0)}=c_{I}\operatorname{tr}Q^{(0)}\Pi_{I}^{\otimes 4}\mathbb{E}_{U\in\mathcal{U}(\mathcal{H}_{\frac{1}{2}}^{\otimes 4})}U^{\otimes 4}\psi^{\otimes 4}U^{{\dagger}\otimes 4}\,,$$ (76)

with $\Pi_{I}=\ket{0_{s}}\bra{0_{s}}+\ket{1_{s}}\bra{1_{s}}$. We can plug the generic formula shown in Eq. (61) and getting

$$\begin{split}\operatorname{SP}^{(0)}&=c_{I}[d(d+1)(d+2)(d+3)]^{-1}\operatorname{tr}Q^{(0)}\Pi_{I}^{\otimes 4}\sum_{\pi\in S_{4}}\operatorname{T}_{\pi}\\ &=[(d_{I}+3)(d_{I}+2)(d_{I}+1)d_{I}]\operatorname{tr}Q^{(0)}\Pi_{I}^{\otimes 4}\sum_{\pi\in S_{4}}\operatorname{T}_{\pi}^{I}\end{split}$$ (77)

The relevant difference between this calculation and that carried out in Eq. (66) is the representation of the permutation operators $\operatorname{T}^{I}_{\pi}$: they now read

$$\operatorname{T}_{\pi}^{I}=\sum_{i_{s},j_{s},k_{s},l_{s}}\ket{\pi(i_{s}j_{s}k_{s}l_{s})}\bra{i_{s}j_{s}k_{s}l_{s}}\,.$$ (78)

with all indices belonging to $\{\ket{0_{s}},\ket{1_{s}}\}$ as written in Eq. (11). These are indeed operators acting on $\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$, but they are not the full permutation operators of $\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$, since the indices of the sums do not run on all the basis elements of $\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$, but only on the intertwiner basis elements. Substituting this expression (and the one in Eq. (73)) in Eq. (75) we get

$$\operatorname{SP}^{(0)}=[(d_{I}+3)(d_{I}+2)(d_{I}+1)d_{I}]^{-1}\frac{1}{d_{0}^{2}}\sum_{\pi\in S_{4}}\sum_{P\in\mathbb{P}_{4}}\operatorname{tr}P^{\otimes 4}\operatorname{T}^{I}_{\pi}$$ (79)

The previous observation renders the calculation of objects like $\operatorname{tr}Q^{(0)}\operatorname{T}_{\pi}^{I}$ slightly more difficult, since we cannot exploit the completeness relationship of basis elements like we did in Eq. (64), simply because $\{\ket{0_{s}},\ket{1_{s}}\}$ is not a complete basis for $\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$. However, by pure brute-force methods, we are able to evaluate this object. We proceed permutation by permutation, as before: the trace of $Q^{(0)}$ with a single permutation operator reads

$$\begin{split}\operatorname{tr}Q^{(0)}\operatorname{T}^{(s)}_{\pi}=\frac{1}{d_{0}^{2}}\sum_{P}\sum_{i_{s}j_{s}k_{s}l_{s}}\bra{\pi(i_{s})}P\ket{i_{s}}\bra{\pi(j_{s})}P\ket{j_{s}}\bra{\pi(k_{s})}P\ket{k_{s}}\bra{\pi(l_{s})}P\ket{l_{s}}\,.\end{split}$$ (80)

This sum is constituted by products of matrix elements of $16\times 16$ matrices between two intertwiner basis elements, which are $16$-component vectors in the original $\mathcal{H}_{\frac{1}{2}}^{\otimes 4}$ basis. By separately calculating the four possible matrix elements for each and every of the $256$ operators of $\mathbb{P}_{4}$ (namely, $\{\bra{0_{S}}P\ket{0_{s}},\bra{0_{S}}P\ket{1_{s}},\bra{1_{S}}P\ket{0_{s}},\bra{1_{S}}P\ket{1_{s}}\}_{P\in\mathbb{P}_{4}}$), and combining them according to the permutation $\pi$, and them summing over the Pauli operators $P$, one is able to compute objects of the form (80). Repeating this method for the $24$ permutation operators of $S_{4}$ and summing the results one gets

$$\mathbb{E}_{U}\operatorname{tr}Q^{(0)}\psi_{U}^{\otimes 4}=\frac{7}{180}\,,$$ (81)

hence

$$\mathbb{E}_{U}M^{(0)}_{\operatorname{lin}}(\psi_{U})=1-d_{0}\operatorname{tr}Q^{(0)}\psi_{U}^{\otimes 4}=\frac{17}{45}\,,$$ (82)

and finally, we can evaluate the average SE gap, which reads

$$\Delta M(\mathcal{H}_{I})=\mathbb{E}_{U}M^{(0)}_{\operatorname{lin}}(\psi_{U})-\mathbb{E}_{U}M^{(I)}_{\operatorname{lin}}(\psi_{U})=\frac{17}{45}-\frac{1}{5}=\frac{8}{45}\,.$$ (83)

Appendix C SE of 4-valent intertwiner with generic spin

In this section, we introduce a generalized version of the Stabilizer Entropy for qudit systems, following the lines of [48]. This generalization is needed since the dimension of the intertwiner Hilbert space associated with a quantum tetrahedron with all spins equal to $j\in\frac{\mathbb{N}}{2}$ is $2j+1$, according to Peter-Weyl’s theorem. Indeed, we can write explicitly the intertwiner state in (8) as

$$\ket{I} = N\sum_{K=0}^{2j}\sum_{M=-K}^{K}\sum_{\{\vec{m}\}}C^{K,M}_{jm_{1}jm_{2}}C^{K,-M}_{jm_{3}jm_{4}}\ket{\vec{m}}$$ (84)

$$= N\bigg{(}\sum_{\{\vec{m}\}}C^{0,0}_{jm_{1}jm_{2}}C^{00}_{jm_{3}jm_{4}}\ket{\vec{m}}+\sum_{M=-1}^{1}\sum_{\{\vec{m}\}}C^{1,M}_{jm_{1}jm_{2}}C^{1,-M}_{jm_{3}jm_{4}}\ket{\vec{m}}$$

$$+ \ldots+\sum_{M=-2j}^{2j}\sum_{\{\vec{m}\}}C^{2j,M}_{jm_{1}jm_{2}}C^{2j,-M}_{jm_{3}jm_{4}}\ket{m}\bigg{)}$$

$$= N\big{(}\ket{0_{s}}+\ket{1_{s}}+\ldots+\ket{2j_{s}}\big{)}$$

where the states $\{\ket{0_{s}},\ket{1_{s}},\ldots,\ket{2j_{s}}\}$ form a set of mutually orthogonal basis elements in the singlet subspace of the Peter-Weyl decomposition of the original Hilbert space. Hence, the projector on the gauge invariant subspace is

$$\Pi_{I} = \ket{0_{s}}\bra{0_{s}}+\ket{1_{s}}\bra{1_{s}}+\dots+\ket{2j_{s}}\bra{2j_{s}}$$ (85)

$$= \sum_{l=0}^{2j}\Pi_{l}$$

For $j\not=1/2$, this implies that the intertwiner Hilbert space is not a $2$-level system; thus, the usual formulation of stabilizer entropy for qubit systems cannot be applied to such a case.

In order to establish a definition of Stabilizer entropy for generic qudit systems, we introduce the generalization of Pauli group, namely the discrete Heisenberg-Weyl group. Consider an $l$-dimensional Hilbert space $\mathcal{H}\simeq\mathbb{C}^{l}\,,\,l\in\mathbb{N}$ and the space of linear operators acting on it, namely $\mathcal{L}(\mathcal{H})$: we define the boost and shift operators $X,Z\in\mathcal{L}(\mathcal{H})$ as

$$X\ket{j} =\ket{j\oplus_{l}1}$$ (86)

$$Z\ket{j} =\omega^{j}\ket{j}\,,$$ (87)

with $\omega=e^{\frac{2\pi i}{l}}$. The Heisenberg-Weyl operators are defined as

$$D_{(p,q)}=\tau^{-pq}Z^{p}X^{q}$$ (88)

with $\tau=e^{\frac{i}{l}\pi}$. They form a basis for $\mathcal{L}(\mathcal{H})$, given by the $l^{2}$ operators $\{D_{(0,0)}=Id_{l^{2}},\,D_{(p,q)}\lvert\,\,\,p,q\in\mathbb{Z}_{l}\}$, where the orthogonality relation reads

$$\operatorname{tr}(D_{(p,q)}D^{\dagger}_{(p^{\prime},q^{\prime})})=l\delta_{p,p^{\prime}}\delta_{q,q^{\prime}}.$$ (89)

Finally, the Heisenberg-Weyl group is defined as the group generated by such operators:

$$\mathcal{D}^{(l)}_{1}=\left<\{D_{(p,q)}\}\right>$$ (90)

Notice that for $l=2$ the Heisenberg-Weyl group reduces to the usual Pauli group for qubit systems.

The Heisenberg-Weyl group acting on a $n-$qudit system is simply given by the tensor product of $n$ copies of the Heisenberg-Weyl group for one qudit.

$$\mathcal{D}^{(l)}_{n}=\mathcal{D}^{(l)\otimes n}_{1}$$ (91)

Any $n-$qudit state $\psi\in\mathcal{H}^{\otimes n}$ can be written in the Heisenberg-Weyl basis as follows

$$\psi=\frac{1}{d}\sum_{\vec{p}\in\mathbb{Z}_{l}^{n}}\sum_{\vec{q}\in\mathbb{Z}_{l}^{n}}\operatorname{tr}(\psi D_{(\vec{p},\vec{q})})D_{(\vec{p},\vec{q})}\,,$$ (92)

with $d:=\operatorname{dim}(\mathcal{H}^{\otimes n})=l^{n}$. We can define the normalized expected value of a $n-$qudit pure state over the discrete Heisenberg–Weyl operators as $\Xi_{D_{(\vec{p},\vec{q})}}(\ket{\psi}):=d^{-1}\bra{\psi}D_{(\vec{p},\vec{q})}\ket{\psi}^{2}$ and the $\alpha$-Stabilizer Rényi Entropy on qudit systems reads

$$M_{\alpha}(\ket{\psi})=(1-\alpha)^{-1}\log\sum_{\vec{p}\in\mathbb{Z}_{l}^{n}}\sum_{\vec{q}\in\mathbb{Z}_{l}^{n}}\Xi_{D_{(\vec{p},\vec{q})}}^{\alpha}(\ket{\psi})-\log d\,.$$ (93)

In order to compute the SE of a quantum tetrahedron state with $j\not=1/2$, we can refer to (8) and write the intertwiner density matrix and its tomography in the Heisenberg-Weyl basis

$$I=\sum_{K,K^{\prime},M,M^{\prime}}\sum_{\{\vec{m}\}\{\vec{m}^{\prime}\}} C^{*\,K^{\prime},M^{\prime}}_{jm^{\prime}_{1}jm^{\prime}_{2}}C^{*\,K^{\prime},-M^{\prime}}_{jm^{\prime}_{3}jm^{\prime}_{4}}C^{K,M}_{jm_{1}jm_{2}}C^{K,-M}_{jm_{3}jm_{4}}\ket{\vec{m}}\bra{\vec{m}^{\prime}}\,,$$ (94)

$$I=\frac{1}{d}\sum_{D_{(\vec{p},\vec{q})}\in D_{4}}\sum_{K,K^{\prime},M,M^{\prime}}\sum_{\{\vec{m}\}\{\vec{m}^{\prime}\}} C^{*\,K^{\prime},M^{\prime}}_{jm^{\prime}_{1}jm^{\prime}_{2}}C^{*\,K^{\prime},-M^{\prime}}_{jm^{\prime}_{3}jm^{\prime}_{4}}C^{K,M}_{jm_{1}jm_{2}}C^{K,-M}_{jm_{3}jm_{4}}\operatorname{tr}(\ket{\vec{m}}\bra{\vec{m}^{\prime}}D_{(\vec{p},\vec{q})})D_{(\vec{p},\vec{q})}\,,$$ (95)

where the sums over $K$ and $K^{\prime}$ run from $0$ to $2j$ and the ones over $M$ and $M^{\prime}$ run respectively from $-K$ to $K$ and from $-K^{\prime}$ to $K$; $d=(2j+1)^{4}$ is the dimension of the Hilbert space over which the trace is performed.