---
title: Random Phased Subspace States
url: https://www.emergentmind.com/topics/random-phased-subspace-states
type: topic
---

# Random Phased Subspace States

Searching arXiv for the cited papers to ground the article in current arXiv records.
Random phased subspace states are most naturally understood as pure-state ensembles obtained by fixing an orthonormal basis and an amplitude profile on a chosen support subspace, and randomizing only the phases. In the formalism of phase-random states, one fixes \(\{ |u_n\rangle \}_{n=1}^{2^N}\), amplitudes \(r_n \ge 0\) with \(\sum_n r_n^2=1\), and draws phases \(\varphi_n\) independently and uniformly from \([0,2\pi]\), giving
\[
|\phi\rangle=\sum_{n=1}^{2^N} r_n e^{i\varphi_n}|u_n\rangle.
\]
If \(r_n=0\) for \(n\notin S\), then every state is supported on \(\mathcal H_S:=\mathrm{span}\{|u_n\rangle:n\in S\}\), yielding the fixed-support specialization most closely aligned with the phrase “random phased subspace states” [1111.2747]. Closely related, but technically distinct, are full-support binary-phase ensembles, randomized subspace-covering ensembles for degenerate eigenspaces, and random subset-phase states generated by pseudochaotic dynamics [1906.10611], [2604.09483], [2410.21268].

## 1. Definition and ensemble structure

The foundational ensemble is the phase-random state ensemble
\[
\Upsilon_{\rm phase}=\{\,|\phi\rangle\,\}_{\mathrm d\varphi},
\qquad
|\phi\rangle=\sum_{n=1}^{2^N} r_n e^{i\varphi_n}|u_n\rangle,
\]
with phase measure
\[
\mathrm d\varphi=\frac{\mathrm d\varphi_1}{2\pi}\cdots \frac{\mathrm d\varphi_{2^N}}{2\pi}.
\]
The defining structural constraint is that amplitudes \(r_n\) and basis vectors \(|u_n\rangle\) are fixed, whereas the phase vector is uniformly distributed on the torus. The ensemble-average density matrix is
\[
\hat\Phi:=\int \mathrm d\varphi\,|\phi\rangle\langle\phi|
=\sum_n r_n^2 |u_n\rangle\langle u_n|,
\]
so phase averaging removes off-diagonal terms in the chosen basis [1111.2747].

This formalism already contains the subspace-supported case. If \(r_n=0\) outside a prescribed support set \(S\), then the states lie in
\[
\mathcal H_S=\mathrm{span}\{|u_n\rangle:n\in S\}.
\]
When amplitudes are equal on the support, \(r_n=|S|^{-1/2}\) for \(n\in S\), the ensemble becomes the natural equal-amplitude random-phase model on a fixed subspace. The construction is basis-dependent: the notion of “subspace support” is defined relative to the basis in which phases are randomized [1111.2747].

The surrounding literature uses nearby constructions that should not be conflated.

| Ensemble | State form | Structural distinction |
|---|---|---|
| Phase-random state | \(\sum_n r_n e^{i\varphi_n}|u_n\rangle\) | Fixed amplitudes, continuous random phases in a fixed basis |
| Binary-phase full-support state | \(2^{-n/2}\sum_x(-1)^{f(x)}|x\rangle\) | Full support, real amplitudes, binary phases \(\{\pm1\}\) |
| QRSI branch state | \(R_i\,\mathcal{P}^{(q)}_{\omega_i}(R_i^\dagger H R_i;\mathcal R)\) | Random ambient-space rotations induce random directions in a degenerate eigenspace |
| Random subset-phase state | \(\sum_{b'}u_{b',b}(-1)^{f(ba)+f(b'a)}|p(b'a)\rangle\) | Support is a random subset of size \(2^k\), not generally a linear subspace |

A recurrent misconception is that all such ensembles are interchangeable. They are not. Fixed-support phase-random states randomize phases while preserving a predetermined amplitude profile; binary-phase constructions are full-support analogues; QRSI randomizes directions in a target eigenspace through branchwise unitaries rather than through phase-only perturbations; and random subset-phase states replace subspace support by a hidden computational-basis subset.

## 2. Subspace support, reduced states, and canonical typicality

An explicit subspace-restricted instance arises when the fixed basis is the Hamiltonian eigenbasis and the support is an energy shell
\[
\mathcal H_R=\mathrm{span}\left\{|e_\alpha\rangle \,\middle|\, e-\delta e<e_\alpha<e+\delta e\right\}\subset \mathcal H_S\otimes\mathcal H_E.
\]
If the initial state lies in \(\mathcal H_R\), then the corresponding phase-random ensemble has support only on that shell. In this setting, long-time dynamics under a time-independent Hamiltonian
\[
H=\sum_{n=1}^d e_n |e_n\rangle\langle e_n|
\]
is related to phase randomization through the evolution
\[
|\phi(t)\rangle=\sum_n r_n e^{-i e_n t/\hbar+i\varphi_n}|e_n\rangle,
\]
under the assumption of phase ergodicity [1111.2747].

For a bipartite split \(\mathcal H_S\otimes\mathcal H_E\), the reduced subsystem state associated with shell-supported phase-random states is
\[
\hat\rho_S=\sum_{\alpha=1}^{d_R} r_\alpha^2\,\mathrm{Tr}_E\bigl(|e_\alpha\rangle\langle e_\alpha|\bigr),
\]
where \(d_R=\dim\mathcal H_R\). The canonical comparison state is the reduced microcanonical state
\[
\mathrm{Tr}_E\!\left(\frac{\mathbb I_R}{d_R}\right).
\]
The condition derived for \(\hat\rho_S\) to coincide with the canonical state is
\[
\sum_{\alpha,\beta=1}^{d_R}
\left(\tilde r_\alpha^2-\frac{1}{d_R}\right)
\left(\tilde r_\beta^2-\frac{1}{d_R}\right)
\mathrm{Tr}\!\left[\hat e_S^\alpha \hat e_S^\beta\right]=0,
\]
with
\[
\hat e_S^\alpha:=\mathrm{Tr}_E\bigl(|e_\alpha\rangle\langle e_\alpha|\bigr).
\]
This is presented as a trade-off relation between the amplitude profile \(\tilde r_\alpha\) on the shell and the structure of the eigenstates in the restricted subspace [1111.2747].

Two limiting mechanisms satisfy the condition directly. One is equal amplitudes in the shell, \(\tilde r_\alpha^2=1/d_R\). The other is identical subsystem reductions of all eigenstates in the shell. This places random phased subspace states in direct contact with canonical typicality: for energy-shell support, phase randomization alone can make subsystem states typically thermal when the support weights or eigenstate structure satisfy the stated condition.

The same logic extends to any fixed support subspace specified in the phase-randomizing basis. A plausible implication is that “random phased subspace states” are best regarded not as a separate formal class, but as a specialization of phase-random states in which the support restriction carries the physically relevant geometry.

## 3. Entanglement structure, concentration, and circuit generation

For a bipartition \(A|\bar A\), the reduced state is
\[
\hat\phi_A=\mathrm{Tr}_{\bar A}|\phi\rangle\langle\phi|,
\]
and entanglement is measured by the linear entropy
\[
E_L^{(A)}(|\phi\rangle)=1-\mathrm{Tr}\,\hat\phi_A^2.
\]
The von Neumann entanglement is
\[
E^{(A)}(|\phi\rangle)=S(\hat\phi_A),\qquad
S(\rho)=-\mathrm{Tr}(\rho\log\rho),
\]
with lower bound
\[
-\log\bigl(1-E_L^{(A)}(|\phi\rangle)\bigr)\le E^{(A)}(|\phi\rangle).
\]
The exact phase-average formula is
\[
\left\langle E_L^{(A)}\right\rangle_{\Upsilon_{\rm phase}}
=
I_L^{(A)}(\hat\Phi)-\sum_{n=1}^{2^N} r_n^4\,E_L^{(A)}(|u_n\rangle),
\]
where
\[
I_L^{(A)}(\hat\Phi)
=
S_L(\mathrm{Tr}_{\bar A}\hat\Phi)+S_L(\mathrm{Tr}_{A}\hat\Phi)-S_L(\hat\Phi).
\]
This decomposition separates the contribution of the dephased mixed state \(\hat\Phi\) from a correction weighted by the entanglement of the basis vectors themselves [1111.2747].

Support size enters sharply. If the basis is separable and only \(R\) amplitudes are nonzero, then
\[
0\le \left\langle E_L^{(A)}\right\rangle_{\Upsilon_{\rm phase}^{\rm sep}}
\le S_L(\hat\Phi^{\rm sep}),
\qquad
\hat\Phi^{\rm sep}=\sum_n r_n^2 |u_n^{\rm sep}\rangle\langle u_n^{\rm sep}|.
\]
If the number of nonzero coefficients is \(R\), then
\[
S_L(\hat\Phi^{\rm sep})\le 1-\frac{1}{R}.
\]
This makes the support cardinality itself an upper-bound parameter for average linear entropy in the separable-basis case. A small support restricts the attainable phase-averaged entanglement unless the support basis states are already entangled [1111.2747].

For equal amplitudes in a separable basis, the ensemble exhibits both high average entanglement and concentration. The paper proves
\[
\mathrm{Prob}\!\left[
\Delta E_L^{(A)}(|\phi^{\rm eq,sep}\rangle)>
\frac{2}{2^N}+\epsilon
\right]
\le e^{-c\epsilon^4 2^N},
\qquad
c=\frac{1}{2^{11}\pi^2},
\]
where
\[
\Delta E_L^{(A)}(|\phi\rangle)=
\left|
E_L^{(A)}(|\phi\rangle)
-
\left\langle E_L^{(A)}\right\rangle_{\Upsilon_{\rm phase}^{\rm eq,sep}}
\right|.
\]
Combined with the lower bound
\[
\langle E^{(A)}\rangle_{T;\infty}
\ge
N_A-\log\!\left(1+2^{2N_A-N}-2^{N_A-N}\right),
\]
this links phase-random subspace ensembles to volume-law entanglement and matrix-product-state hardness in the large-support regime [1111.2747].

The same work also gives an explicit diagonal-circuit mechanism. Starting from
\[
|\phi_0\rangle=\sum_a r_a e^{i\varphi_a}|\bar a\rangle,
\]
one applies random two-qubit diagonal gates
\[
W_t=CZ_{i_t j_t}P_{i_t}(\alpha_t)P_{j_t}(\beta_t),
\]
with \(CZ_{ij}=\mathrm{diag}(1,1,1,-1)\) and \(P_k(\theta)=\mathrm{diag}(1,e^{i\theta})\), producing
\[
|\phi_T\rangle=W_T\cdots W_1|\phi_0\rangle.
\]
The asymptotic average linear entropy converges to that of the computational-basis phase-random ensemble:
\[
\lim_{T\to\infty}\mathbf E\!\left[E_L^{(A)}(|\phi_T\rangle)\right]
=
\left\langle E_L^{(A)}\right\rangle_{\Upsilon_{\rm phase}^{\rm comp}}.
\]
This establishes that phase-random subspace behavior can be generated by diagonal random circuits, not only by abstract sampling over phase variables [1111.2747].

## 4. Binary-phase full-support ensembles as a limiting contrast

A distinct but closely related ensemble replaces continuous phases on a fixed support by random binary signs on the entire computational basis:
\[
|\psi_f\rangle
=
2^{-n/2}\sum_{x\in\{0,1\}^n}(-1)^{f(x)}|x\rangle,
\qquad
f:\{0,1\}^n\to\{0,1\}.
\]
Equivalently, the amplitudes are real and take values \(\pm 2^{-n/2}\). This is a full-support, equal-magnitude, binary-phase ensemble. It is therefore a special full-space analogue of phased subspace states, not a subspace-state result [1906.10611].

The main theorem states that the distribution
\[
2^{-n/2}\sum_{x\in\{0,1\}^n}(-1)^{f(x)}|x\rangle
\]
is a
\[
\frac{4t^2}{2^n}\text{-approximate state }t\text{-design for all }t.
\]
Writing
\[
\rho_{(t,n,2)}
=
\mathbb{E}_{f:\{0,1\}^n\to\{0,1\}}
\big[(|f\rangle_{(2)}\langle f|_{(2)})^{\otimes t}\big]
\]
and \(\rho_{(t,n,H)}\) for the \(t\)-copy Haar moment state, the exact comparison given is
\[
TD\big(\rho_{(t,n,2)},\rho_{(t,n,H)}\big)\le \frac{4t^2}{2^n}.
\]
Thus for any polynomial \(k=t(n)\),
\[
TD\big(\rho_{(k,n,2)},\rho_{(k,n,H)}\big)\le 4k^2\,2^{-n}.
\]
The result shows that binary phases \(\{\pm1\}\) already suffice for near-Haar \(t\)-copy behavior in the full-support setting [1906.10611].

The contrast with random phased subspace states is structural rather than superficial. The support is always the entire computational basis, and the proof relies on that full-support structure. The entrywise moment formulas are expressed on tuples from \((\{0,1\}^n)^t\), and the combinatorics distinguishes exact multiplicity matching for complex phases from parity-of-multiplicity matching for binary phases. A plausible implication is that analogous parity-based moment mechanisms may exist for a fixed subspace support, but that extension is not proved in the paper [1906.10611].

Implementation is unusually weak. If \(F_k\) is a keyed Boolean function, the state
\[
|\phi_k\rangle=2^{-n/2}\sum_x (-1)^{F_k(x)}|x\rangle
\]
is generated by a Hadamard layer followed by reversible evaluation of \(F_k\) and phase kickback. If \(F\) is computable by a classical circuit of size \(s(n)\) and depth \(d(n)\), then the generator is computable by an HT circuit—one Hadamard layer followed by Toffoli gates—of size \(O(s(n))\) and depth \(d(n)+1\). The same construction yields pseudorandom states from post-quantum pseudorandom functions and \(\frac{4t^2}{2^n}\)-approximate state \(t\)-designs from \((2t)\)-wise independent Boolean functions [1906.10611].

## 5. Randomized subspace-covering ensembles and the phase-only boundary

The sharpest distinction between random phased subspace states and other randomized subspace ensembles is provided by quantum randomized subspace iteration. Here the target object is a degenerate eigenspace
\[
\mathcal G=\operatorname{eigenspace}(H,E_\star)=\operatorname{span}\{\ket{v_1},\dots,\ket{v_g}\},
\]
with ground-state specialization \(E_\star=E_0\) and spectral gap
\[
\gamma=E_{g+1}-E_0>0.
\]
The problem is to generate states \(\{\ket{\psi_i}\}_{i=1}^M\) such that each has substantial overlap with \(\mathcal G\) and the projected vectors span \(\mathcal G\) [2604.09483].

QRSI does not randomize phases in a fixed basis. Instead it samples independent branchwise random unitaries \(R_i\in U(N)\), defines rotated Hamiltonians
\[
H_i=R_i^\dagger H R_i,
\]
runs an eigenstate-preparation primitive on each branch, and maps the outputs back:
\[
\ket{\psi_i}
=
R_i\,\mathcal{P}^{(q)}_{\omega_i}\!\left(R_i^\dagger H R_i;\mathcal R\right).
\]
The resulting target manifold on branch \(i\) is
\[
\mathcal G_i=\operatorname{eigenspace}(H_i,E_0)=R_i^\dagger\mathcal G.
\]
The paper emphasizes that this is more general than random phases: each branch presents “a randomly oriented copy of \(\mathcal G\)” to the preparation primitive [2604.09483].

The key subspace object is the foot-point matrix
\[
F\in\mathbb C^{g\times M},
\qquad
F_{k,i}=\braket{v_k}{R_i|\tilde\phi_i},
\]
so that
\[
\Pi_{\mathcal G}R_i\ket{\tilde\phi_i}
=
\sum_{k=1}^g F_{k,i}\ket{v_k}.
\]
This gives each branch a random coefficient vector in \(\mathbb C^g\), but the coefficients do not arise from random phases over a fixed deterministic amplitude profile. They arise from the random orientation induced by \(R_i\) together with the branchwise preparation dynamics [2604.09483].

Under independent Haar-random branch rotations and a preparation primitive with nonzero projection onto the rotated target manifold, Proposition 1a states that for \(M\ge g\),
\[
\operatorname{rank}(F)=g \quad \text{a.s.}
\]
Equivalently, the projected ensemble spans \(\mathcal G\) almost surely. The proof uses a \(U(g)\)-equivariance argument inside the degenerate eigenspace, not a random-phase ansatz. The normalized foot-point is \(U(g)\)-invariant and hence uniform on the unit sphere \(S^{2g-1}\subset\mathcal G\), so the construction yields independent continuous random directions in \(\mathbb C^g\) [2604.09483].

The paper further weakens Haar randomness to an anti-concentration condition on the output foot-point distribution. If for every hyperplane \(V\subset\mathcal G\),
\[
\Pr_{R\sim\nu}[d_{FS}(\hat w,V)\ge \delta]\ge \eta,
\]
then \(M=g\) independent branches give
\[
\Pr[\operatorname{rank}(F)=g]\ge \eta^g,
\]
and \(1-\varepsilon\) success follows with
\[
M\le g\lceil \eta^{-g}\log(1/\varepsilon)\rceil.
\]
This establishes manifold coverage rather than exact Haar sampling on the subspace in every implementation [2604.09483].

A second misconception addressed by QRSI concerns spectral deformation. Randomization does not split the degeneracy or alter the branch difficulty: for any \(R_i\in U(N)\),
\[
\operatorname{spec}(H_i)=\operatorname{spec}(H),
\qquad
\gamma_i=\gamma.
\]
The recovered manifold is then identified through the corrected-state coefficient matrix
\[
C=[\bm c_1|\cdots|\bm c_M]
\]
or the Gram matrix
\[
G_{ij}=\braket{\psi_i}{\psi_j}.
\]
With branch leakage
\[
\varepsilon_q:=\max_i \|\Pi_{\mathcal G^\perp}\ket{\psi_i}\|^2,
\]
Proposition 3 gives
\[
\sigma_{g+1}(C)\le \sqrt{M\,\varepsilon_q},
\qquad
\sigma_g(C)\ge \sigma_g(F)-\sqrt{M\,\varepsilon_q}.
\]
The conceptual outcome is that QRSI is best understood as a randomized subspace-covering ensemble. It is relevant to random phased subspace states only in a broad sense: both produce random coefficient vectors in a target subspace, but the QRSI mechanism is random basis rotation in the ambient space, not phase-only randomization [2604.09483].

## 6. Random subset-phase states and the broader pseudorandom landscape

A further nearby notion is the random subset-phase state generated by random subsystem-embedded dynamics. The construction uses a subsystem size \(k\le n\), a random subset isometry
\[
O_a=\sum_{b\in\{0,1\}^k}(-1)^{f(ba)}\ketbra{p(ba)}{ba},
\]
and a full unitary
\[
U=\sum_{a\in\{0,1\}^{n-k}} O_a u O_a^\dagger.
\]
For an initial computational basis state \(\ket{p(ba)}\), the evolved state is
\[
U\ket{p(ba)}
=
\sum_{b'\in\{0,1\}^k}
u_{b',b}(-1)^{f(ba)+f(b'a)}\ket{p(b'a)}.
\]
Its support is the set \(\{p(b'a):b'\in\{0,1\}^k\}\), of size \(2^k\). This is a random computational-basis subset, not generally a linear subspace [2410.21268].

The especially important example is
\[
u=H^{\otimes k}P,
\]
where \(P\) is a random sign operator and the elements of \(u\) are “exactly \(2^{-k/2}\) up to signs.” In that case the evolved state is a uniform-magnitude random signed superposition over a subset of size \(2^k\). The phases are discrete signs rather than continuous \(e^{i\theta}\), and the support is subset-structured rather than linear-subspace-structured. Accordingly, the closest description is “random subset-phase state,” not random phased linear-subspace state [2410.21268].

The paper’s formal guarantees are pseudorandom rather than exact-Haar. One theorem states that if \(\mathcal E_k\) is an ensemble of embedded unitaries satisfying, for \(K=2^k\),
\[
\operatorname{Pr}\left[\abs{u_{b,b'}^2}\ge K^{-\epsilon}\right]\le \operatorname{negl}(n)
\]
for some \(\epsilon>0\), together with
\[
\mathbb E_{u\sim\mathcal E_k}\left[\abs{u_{b,b'}^2}\right]=K^{-1}
\]
for all \(b,b'\), then the corresponding RSED ensemble generates pseudorandom states. A second theorem states that “RSEDs with an ensemble of embedded unitary operators that generate a pseudorandom state ensemble in the subspace generate a pseudorandom state ensemble in the entire space” [2410.21268].

The same framework gives a pseudochaotic dynamical interpretation through the out-of-time-ordered correlator
\[
O_{VW}(U)=\frac{1}{2^n}\operatorname{tr}\left(VUWU^\dagger VUWU^\dagger\right).
\]
For \(k=\omega(\log n)\), local-operator OTOCs are negligible with probability higher than \(1-\operatorname{negl}(n)\) if the maximum averaged magnitude of the embedded operator elements is \(O(2^{-k/2})\). Circuit complexity is also explicit: “a \(\omega(\log n)O(t[\log t]^2)\) depth local circuit can implement a pseudochaotic RSED with an error negligible in \(n\) up to time \(\tau=O(n)\),” and the overview states that a circuit implementation “requires \(\operatorname{polylog}(n)\) depth with all-to-all connectivity” [2410.21268].

Within the broader taxonomy, random phased subspace states occupy the phase-only, fixed-support corner. Full-support binary-phase states show that very sparse phase alphabets can already yield near-Haar \(t\)-copy behavior, but only in the full-support regime. QRSI shows that randomized coverage of a target manifold can be achieved by branchwise ambient-space rotations rather than by phase-only perturbation. Random subset-phase states show that equal-magnitude random-sign superpositions over hidden support sets can be generated efficiently and can be computationally indistinguishable from Haar-random states under restricted access. The central invariant across these models is randomized coherence on a structured support; the central difference is whether the structure is a fixed subspace, the full space, a degenerate eigenspace reached by random rotation, or a random subset.

Source: https://www.emergentmind.com/topics/random-phased-subspace-states