---
title: Magic–Entanglement Complementarity
url: https://www.emergentmind.com/topics/magic-entanglement-complementarity
type: topic
---

# Magic–Entanglement Complementarity

Searching arXiv for the core papers to ground the article in current literature.
Magic–entanglement complementarity denotes a family of quantitative relations between nonstabilizerness and entanglement, rather than a single universal trade-off. In recent work, complementarity has appeared in at least four distinct forms: entanglement can facilitate the nonlocal spreading of locally injected magic, entanglement and magic can undergo distinct dynamical transitions, local dissipation can irreversibly destroy entanglement while later reviving magic, and operator entanglement can be upper-bounded by magic monotones in the Heisenberg picture [2503.20873; 2312.02039; 2605.22603; 2501.18679]. Taken together, these results place magic and entanglement in a joint resource-theoretic landscape in which they are neither reducible to one another nor generically monotone with each other.

## 1. Resource-theoretic setting and measures

A pure stabilizer state is stabilized by an Abelian subgroup of the \(N\)-qubit Pauli group with \(2^N\) elements; such states are prepared by Clifford circuits and admit efficient classical simulation. Magic quantifies deviation from stabilizer structure, and several recent works use stabilizer-entropic measures as the principal diagnostics. In the stabilizer-state setting, the linear stabilizer entropy is
\[
Y(\rho) := 2^{-N} \sum_{P\in P_N} \operatorname{tr}(P\rho)^4,
\]
with \(Y(\rho)=1\) iff \(\rho\) is a stabilizer state, \(Y(\rho)<1\) otherwise, and \(M_2(\rho)=-\log Y(\rho)\) the second stabilizer Rényi entropy. The more general family
\[
Y_\alpha(\rho) := 2^{-N} \sum_{P\in P_N} \operatorname{tr}(P\rho)^{2\alpha}, \qquad
M_\alpha(\rho) := \frac{1}{1-\alpha}\log Y_\alpha(\rho)
\]
is used to track the leading dependence of magic under local injections [2503.20873].

Entanglement is quantified differently across settings. For bipartite stabilizer states on \(AB\), the von Neumann entropy \(E=-\operatorname{tr}(\rho_A\log_2\rho_A)\) is an integer and counts Bell pairs across the cut after local Clifford reduction. In monitored circuits, the same entropy across a half cut is used as the entanglement order parameter. In two-qubit studies, entanglement is quantified by the concurrence \(\Delta\), while dissipative GHZ analyses use bipartite negativity and separability thresholds. In operator-space formulations, local operator entanglement (LOE) is the entanglement entropy of the vectorized Heisenberg operator across a doubled bipartition \((H_A\otimes H_A)\otimes(H_{\bar A}\otimes H_{\bar A})\) [2312.02039; 2603.24902; 2605.22603; 2501.18679].

The same plurality holds on the magic side. Hybrid-circuit work uses the stabilizer \(\alpha\)-Rényi entropy defined from the Pauli-string expectation distribution \(\Xi_P(|\psi\rangle)=2^{-N}\langle\psi|P|\psi\rangle^2\). Dissipative studies use the robustness of magic \(\mathcal R(\rho)\), defined by the minimal signed-stabilizer decomposition cost. Two-qubit Pareto analyses use \(M_2\), and operator-space work uses operator stabilizer Rényi entropy, unitary nullity, and \(T\)-count. This suggests that “complementarity” is measure-dependent and protocol-dependent, not a single invariant law.

## 2. Entanglement as a conduit for magic spreading

The clearest “highway” formulation appears for an initial pure stabilizer state \(|\psi\rangle\) on \(AB\) with bipartite entanglement \(E\), followed by a Haar random unitary \(U_A\) on a local subregion \(A\). To leading order for \(|A|\gg 1\),
\[
\overline{Y}
=
\mathbb E_{U_A}\,Y(U_A|\psi\rangle\langle\psi|U_A^\dagger)
=
4\cdot 2^{-|A|-E}\,[1+O(2^{-|A|-E})].
\]
Since lower \(Y\) means higher magic, increasing \(E\) exponentially suppresses \(\overline{Y}\) and therefore enhances global magic generation. The same leading factor persists for the stabilizer Rényi family,
\[
\overline{Y_\alpha}=(1+c_\alpha)2^{-|A|-E},\qquad c_2=3,\qquad \overline{Y_\infty}=2^{-|A|-E},
\]
so \(M_\alpha\) grows linearly in \(|A|+E\) at leading order [2503.20873].

The mechanism is stabilizer counting. Stabilizers fully supported on \(B\) commute through \(U_A\) and contribute coherently to \(Y\). For a bipartite stabilizer with entanglement \(E\), the number of stabilizers supported on \(B\) is \(2^{|B|-E}\). After normalization by \(2^N\), this yields the factor \(2^{-|A|-E}\). In the operator-spreading picture, a Pauli on \(A\) conjugated by \(U_A\) proliferates into many Pauli strings, while the \(E\) logical pairs between \(A\) and \(B\) make that mixing extend into \(B\). The result is delocalization of locally injected magic across the full system.

A second complementarity appears when independent Haar random unitaries act on both sides. For \(U_A\otimes U_B\) and \(|A|,|B|\gg 1\),
\[
\overline{Y}
=
4\cdot 2^{-N}
\Bigl[
1+3\cdot 2^{-2E}
+
O(2^{-|A|-E})
+
O(2^{-|B|-E})
\Bigr].
\]
Because the Haar-random pure-state value is \(4\cdot 2^{-N}\), the deviation from global Haar magic is controlled solely by \(3\cdot 2^{-2E}\). Numerically, \(E\approx 3\) already saturates to the Haar value for all \(N\). In this regime, entanglement suppresses the residual gap left by locality, so a product local unitary can act “as if” it were global Haar on the metric \(Y\) [2503.20873].

The same qualitative structure extends beyond the bipartite stabilizer case. For tripartite stabilizer states, the relevant data are the GHZ count \(g\), pairwise Bell counts \(b_{AB},b_{AC},b_{BC}\), and local singles. Under \(U_A\otimes U_B\), the leading exponent is controlled by the acted region \(|A|+|B|\) together with the entanglement across \(AB|C\), while the product-vs-global gap is suppressed by \(2^{-2b_{AB}-g}\). Non-stabilizer entanglement built from imperfect Bell pairs, and shallow brickwork circuits replacing Haar \(U_A\), preserve the same qualitative conclusion: increasing total entanglement enhances global magic spreading, although finite-depth architectures exhibit saturation effects at very small depths [2503.20873].

## 3. Distinct dynamical regimes and phase structure

In hybrid Clifford+\(T\) circuits with measurements, complementarity takes the form of separated critical behavior. The architecture is a 1D chain of \(N\) qubits with a brickwork pattern of random two-qubit Clifford gates on nearest neighbors, projective measurements in the computational basis with probability \(p\) per qubit per time step, and \(T\)-gate injections with probability \(q(N)=\eta/N^\beta\). For the principal case \(\beta=1\), entanglement undergoes the usual monitored-induced transition at
\[
p_c^{\mathrm{entgl}}=0.15995(10),
\]
while magic collapses only at a higher critical rate,
\[
p_c^{\mathrm{magic}}\simeq 0.22
\]
for \(\eta=2.0\). The regime \(p=0.18\) is diagnostic: entanglement already shows area-law scaling, while magic remains (sub)-extensive. The resulting intermediate phase demonstrates that entanglement alone does not diagnose nonstabilizerness or computational hardness in monitored dynamics [2312.02039].

A related decoupling appears in gauge theory. In the \((1+1)\)-dimensional SU(2) lattice gauge theory formulated in a dressed-site basis, gauge-invariant entanglement entropy \(S(A)\) and stabilizer Rényi entropy \(\mathcal M_2\) were computed for systems up to \(L=100\) (\(300\) qubits). At the CFT point \((m,g)=(0,0)\), the entanglement scaling yields \(c=0.990(5)\), while magic obeys
\[
\mathcal M_2(L)=\alpha L+\beta\ln L+\gamma
\]
with \(\alpha=0.29(1)\), \(\beta=-0.69(22)\), and \(\gamma=0.88(59)\). At \(m=0.2\), the entanglement entropy decreases monotonically with \(g\), its sharpest change occurs near \(g_\star\approx 1.9\), and \(\mathcal M_2/L\) shows a broad plateau before dropping rapidly only after the same crossover. The strong-coupling limit \(g\to\infty\) is a product stabilizer state with both \(S(A)\to 0\) and \(\mathcal M_2\to 0\). The intermediate region is therefore “magic-rich but low-entanglement” rather than maximally nonclassical in both senses [2606.09971].

The light-front formulation of the \((1+1)\)-dimensional transverse-field Ising model provides a sharper basis-dependent contrast. In instant-form momentum space, the ground state is a product over \(k>0\) of two-mode BCS-like states,
\[
|\mathrm{GS}_{IF}\rangle
=
\bigotimes_{k>0}
\left(
\cos\phi_k\,|00\rangle_{k,-k}
-
i\sin\phi_k\,|11\rangle_{k,-k}
\right),
\]
which carry pairwise entanglement between \(+k\) and \(-k\). The light-front Hamiltonian is diagonal in \(k^+\), so the light-front ground state is the separable Fock vacuum \(\bigotimes_{k^+}|0\rangle_{k^+}\). Away from criticality, the instant-form ground state has nonzero momentum-space magic,
\[
M^{IF}_2
=
-\sum_{k>0}
\ln\!\left(
1-\left(\frac{k\,m}{k^2+m^2}\right)^2
\right),
\qquad
M^{LF}_2=0.
\]
At the quantum critical point, however, both ground states are stabilizers: the instant-form state is a product of maximally entangled Bell-like pairs in momentum space, whereas the light-front ground state remains separable. High entanglement with zero magic is therefore possible, and the quantization scheme itself determines which resource is used [2507.10777].

## 4. Dissipative complementarity and reborn magic

Under local amplitude damping, magic and entanglement respond in qualitatively different ways because separability is preserved by local CPTP maps, whereas stabilizer membership is not. For the \(n\)-qubit GHZ family
\[
|\psi_n\rangle=\alpha|0^n\rangle+\beta|1^n\rangle,\qquad 0<\alpha<\beta,
\]
the evolved state under \(\mathcal E_\gamma^{\otimes n}\) has only one surviving off-diagonal GHZ coherence \(c=\alpha\beta(1-\gamma)^{n/2}\). On the real GHZ–X slice, stabilizer membership is exact:
\[
\rho\in\mathcal S
\quad\Longleftrightarrow\quad
|c|\le \min(p_{0^n},p_{1^n}).
\]
This produces two stabilizer thresholds: the lower entry \(\gamma_-^{(n)}\), defined by \(P_0=c\), and the upper exit
\[
\gamma_+^{(n)}=1-r^{2/n},\qquad r=\alpha/\beta.
\]
Magic is absent exactly on \([\gamma_-^{(n)},\gamma_+^{(n)}]\) and present on the two open branches \([0,\gamma_-^{(n)})\cup(\gamma_+^{(n)},1)\) [2605.22603].

Entanglement, by contrast, exhibits sudden death at a single threshold independent of the bipartition:
\[
\gamma_e^{(n)}=r^{2/n}.
\]
At this same point the state is fully separable. The central complementarity identity is therefore
\[
\gamma_e^{(n)}+\gamma_+^{(n)}=1,
\]
valid for every \(n\ge 2\) in the re-entrant regime \(r<1\). Its origin is the system–environment duality of amplitude damping,
\[
\rho_E(\gamma)=\rho_S(1-\gamma),
\]
which mirrors the system-side magic-rebirth condition against the environment-side entanglement-death condition. In Regimes I and II of the ordering analysis, the entire reborn branch lies in fully separable states, and all proper marginals are stabilizer. Reborn magic is then nonlocal but not entanglement-bearing [2605.22603].

This nonlocal reborn magic can nevertheless be concentrated. Measuring the \(n-1\) commuting parity stabilizers \(Z_iZ_{i+1}\) and postselecting the trivial syndrome projects onto the GHZ subspace with success probability
\[
p_{\rm succ}=P_0+P_n\ge \alpha^2.
\]
A subsequent CNOT cascade and discarding of spectators yields a single-qubit state \(\tilde\rho(\gamma)\). The extraction is lossless in expectation for the robustness of magic:
\[
(P_0+P_n)\bigl[\mathcal R(\tilde\rho)-1\bigr]
=
\mathcal R\bigl(\rho_n(\gamma)\bigr)-1.
\]
For sufficiently large \(n\), the decoded branch enters the standard distillable regions for \(|H\rangle\)-type and \(|T\rangle\)-type protocols [2605.22603].

The same dissipative framework also splits pure stabilizer inputs into magic-generators and magic-insulators. Under homogeneous local amplitude damping, a pure stabilizer state is an insulator iff its computational-basis support has constant Hamming weight; otherwise it is a generator. At two qubits, the Bell states
\[
\mathcal R_{\Phi^+}(\gamma)=1+\gamma(1-\gamma),\qquad
\mathcal R_{\Psi^+}(\gamma)=1
\]
exemplify the split: \(|\Phi^+\rangle\) exits the stabilizer polytope immediately for any \(0<\gamma<1\), whereas \(|\Psi^+\rangle\) remains stabilizer for all \(\gamma\) [2605.22603].

## 5. Operator-space, complexity, and algorithmic consequences

In the Heisenberg picture, complementarity becomes a hierarchy of rigorous inequalities. For an \(N\)-qubit unitary \(U\), a nontrivial Pauli \(O\), and a bipartition \(A:\bar A\), the local operator entanglement of the Heisenberg-evolved operator \(O_U=U^\dagger O U\) satisfies
\[
E_A^{(\alpha)}(O_U)\le M^{(\alpha)}(O_U)\le \nu(U)\le \tau(U),
\]
where \(M^{(\alpha)}\) is the operator stabilizer Rényi entropy, \(\nu(U)\) is unitary nullity, and \(\tau(U)\) is the \(T\)-count. Thus large LOE is impossible without large magic in operator space. For random ensembles, the same quantities nearly coincide on average. For the \(T\)-doped Clifford ensemble \(\mu_\tau\),
\[
\int_{U\sim\mu_\tau} E_A^{(\mathrm{pur})}(O_U)
=
\exp\!\left(
-\log(4/3)\,\tau + O\!\left((4/3)^\tau/D\right)
\right),
\]
and for the \(\nu\)-compressible ensemble \(\mu_\nu\),
\[
\int_{U\sim\mu_\nu} E_A^{(\mathrm{pur})}(O_U)
=
\exp\!\left(
-\nu + 2 + O(4^{-\nu}+2^\nu/D)
\right).
\]
The implication is operational: any dynamics that are hard for tensor-network simulation because LOE is large must also be hard for stabilizer and Pauli-truncation methods [2501.18679].

A closely related result states that a unitary map generates nonlocal magic if and only if it generates operator entanglement on Pauli strings. On that basis, an average measure of a unitary’s Pauli-entangling power is introduced as a proxy for nonlocal magic generation, with analytical formulas, a typical value, and upper bounds in terms of the nonstabilizerness properties of the evolution [2504.09360].

The computational consequences in the Schrödinger picture are more asymmetric. Using stabilizer nullity \(\nu\) as the magic order parameter, Hilbert space can be partitioned into an entanglement-dominated (ED) phase, where the entanglement \(S_1(\psi_A)\) significantly exceeds the state’s magic, and a magic-dominated (MD) phase, where magic dominates entanglement. For pure states with nullity \(\nu\), one has
\[
n_A-|S_A|-\nu\le S_1(\psi_A)\le n_A-|S_A|,
\]
and more generally
\[
n_A-|S_A|-2\nu\le S_\alpha(\psi_A)\le n_A-|S_A|.
\]
In the ED phase there are sample- and time-efficient, input-agnostic quantum algorithms for entanglement estimation, distillation, and dilution, with near-reversible resource conversion. In the MD phase, constant-relative-error estimation and constant-fraction distillation are provably intractable in general. This establishes a computational phase separation rather than a simple resource inequality [2403.19610].

## 6. Two-qubit frontier geometry and general lessons

For pure two-qubit states, the interplay between concurrence \(\Delta\) and stabilizer Rényi-2 magic \(M_2\) can be solved exactly. The minimal-magic frontier is a single continuous curve,
\[
M_2^{(\min)}(\Delta)
=
-\ln(\Delta^4-\Delta^2+1),
\]
which vanishes at \(\Delta=0\) and \(\Delta=1\) and is strictly positive for \(0<\Delta<1\). A Schmidt-family representative,
\[
|\psi_{\min}(\theta)\rangle=\cos\theta\,|00\rangle+\sin\theta\,|11\rangle,
\qquad
\Delta=\sin 2\theta,
\]
saturates this lower boundary. Partially entangled two-qubit pure states therefore cannot have zero magic when \(M_2\) is the chosen nonstabilizerness measure [2603.24902].

The maximal-magic frontier is piecewise and has three segments:
\[
M_2^{(\max)}(\Delta)
=
\begin{cases}
\ln\!\bigl( 9/(3\Delta^4-2\Delta^3+4) \bigr), & 0\le \Delta \le \Delta_G,\\[4pt]
\ln\!\bigl( 16/(8\Delta^4-8\Delta^2+9) \bigr), & \Delta_G\le \Delta \le \sqrt3/2,\\[4pt]
\ln\!\bigl( 18/(7\Delta^4-6\Delta^2+9) \bigr), & \sqrt3/2\le \Delta \le 1,
\end{cases}
\]
where \(\Delta_G\approx 0.63726445\). The global maximum \(M_2=\ln(16/7)\) occurs at two distinct entanglement values, \(\Delta=1/2\) and \(\Delta=1/\sqrt2\). At \(\Delta=1\), the upper boundary gives \(M_2=\ln(9/5)\), while Bell stabilizers still realize \(M_2=0\). Entanglement is therefore neither a monotone lower bound nor a monotone upper bound for magic, even in the smallest nontrivial bipartite system [2603.24902].

Taken together, the contemporary literature supports a plural notion of magic–entanglement complementarity. In some settings, entanglement is a conduit that exponentially enhances the global spread of local magic. In others, entanglement and magic separate into distinct phases, distinct critical points, or even opposite responses to noise. In operator space, magic bounds entanglement from above; in two-qubit geometry, the feasible region is bounded by exact Pareto frontiers rather than a single curve. A plausible implication is that any general theory of quantum resources must track not only how much entanglement and magic are present, but also where they reside, how they are generated, and which operational task probes them.

Source: https://www.emergentmind.com/topics/magic-entanglement-complementarity