Entangled Frozen States in Quantum Dynamics

• Entangled Frozen States represent regimes where entanglement remains nearly constant over time, achieved via various mechanisms in closed and non-Hermitian quantum systems.
• Various formulations—including interaction-free evolution, eigenstate stabilization, and Hilbert-space fragmentation—explain EFS behavior in spin chains and optical lattices.
• Research on EFS utilizes diagnostics such as Schmidt weight, quantum Fisher information, and logarithmic negativity to verify robust entanglement storage and control.

Searching arXiv for recent and foundational papers on Entangled Frozen States and closely related formulations. arxiv_search(query="\"Entangled Frozen States\" OR entanglement freezing spin chain IFE fragmentation", max_results=10) arxiv_search(query="Entangled Frozen States", max_results=10) Entangled Frozen States (EFS) denote entangled quantum states, or dynamical entanglement regimes, in which entanglement remains exactly constant or approximately time-independent under the relevant evolution. In the literature, the notion appears in several technically distinct but overlapping forms: freezing–thawing plateaus in closed lossless spin chains and optical lattices (Qian et al., 2021), entangled eigenmodes engineered and immobilized by local phase flips in branched XY networks (0708.2794), interaction-free evolving entangled states of bipartite Hamiltonians (Napoli et al., 2014), gap-protected entanglement-storage units prepared by optimal control (Caneva et al., 2011), high-quantum-Fisher-information many-body states stored near stable fixed points of the two-axis countertwisting model (Kajtoch et al., 2015), entangled zero-energy sectors generated by rank-deficient local constraints in quantum Hilbert-space fragmentation (Zhou et al., 6 Apr 2026), and frozen long-distance entangled links in PT-symmetric non-Hermitian spin systems (Ahuja et al., 12 Jun 2026).

1. Definition and conceptual scope

A common operational definition of freezing is dynamical rather than purely kinematic. In the lossless multiparty formulation, a “freezing” process is a “stable” behavior within a time interval where the physical quantity such as entanglement $E(t)$ is approximately independent of time with $\partial E(t)/\partial t \approx 0$, or where the fluctuation is negligible relative to the value of the quantity, $\sigma(E)/\Delta E \ll 1$ (Qian et al., 2021). In that sense, an EFS need not be a strictly stationary vector of the Hamiltonian; it may instead be a regime in which the chosen entanglement monotone stays on a plateau for finite or infinite times.

A stricter formulation arises in the theory of interaction-free evolving states. For a bipartite Hamiltonian $H=H_A+H_B+H_I=H_0+H_I$, a pure state $|\psi\rangle$ is interaction-free evolving if

$$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$

for some real $\alpha$. Equivalently, $|\psi\rangle$ must satisfy $H_I|\psi\rangle=\alpha|\psi\rangle$ together with the invariance conditions encoded by the subspaces

$$\mathcal N_\alpha=\bigcap_{n=0}^{\infty}\ker\!\big((H_I-\alpha\mathbb I)H_0^n\big) .$$

If such a state is entangled across $\partial E(t)/\partial t \approx 0$0, its entanglement is frozen with respect to the interaction, because the evolution is generated only by local free terms up to a global phase (Napoli et al., 2014).

A third formulation is algebraic and many-body: in classically fragmented systems with rank-deficient local Hamiltonian terms, the entangled-frozen subspace inside a mobile classical Krylov sector $\partial E(t)/\partial t \approx 0$1 is

$\partial E(t)/\partial t \approx 0$2

Its normalized vectors are exact zero-energy eigenstates embedded inside mobile sectors; they do not evolve under the Hamiltonian, but they are not product-frozen states and are therefore intrinsically entangled (Zhou et al., 6 Apr 2026).

These formulations are not identical. Some EFS are exact eigenstates, some are exact interaction-free states, and some are approximate long-lived plateaus. The shared feature is persistence of entanglement under continued dynamics.

2. Closed-system and spin-network realizations

In a lossless optical-lattice or spin-chain setting, EFS arise without decoherence. The model of an $\partial E(t)/\partial t \approx 0$3-site optical lattice in the single-occupation insulating phase is mapped, via a Schrieffer–Wolff transformation, to an effective XY spin model in which a distinguished qubit $\partial E(t)/\partial t \approx 0$4 couples to the end of a one-dimensional chain $\partial E(t)/\partial t \approx 0$5. The dynamics remain unitary and closed, and in the one-excitation sector the excitation amplitude at $\partial E(t)/\partial t \approx 0$6 is

$\partial E(t)/\partial t \approx 0$7

Entanglement is quantified by the Schmidt weight $\partial E(t)/\partial t \approx 0$8, leading to explicit expressions $\partial E(t)/\partial t \approx 0$9 and $\sigma(E)/\Delta E \ll 1$0 for the bipartitions $\sigma(E)/\Delta E \ll 1$1 and $\sigma(E)/\Delta E \ll 1$2 (Qian et al., 2021).

In the thermodynamic limit $\sigma(E)/\Delta E \ll 1$3, the amplitude becomes

$\sigma(E)/\Delta E \ll 1$4

so the local dynamics acquire a Bessel-function envelope. The derivatives of $\sigma(E)/\Delta E \ll 1$5 and $\sigma(E)/\Delta E \ll 1$6 become negligible at large finite times, yielding permanent entanglement freezing in a fully coherent, lossless many-body system. For finite $\sigma(E)/\Delta E \ll 1$7, the same model exhibits alternating freezing and thawing. The first freezing onset and first thawing are identified as

$\sigma(E)/\Delta E \ll 1$8

with approximate period

$\sigma(E)/\Delta E \ll 1$9

The first frozen interval therefore lengthens linearly with $H=H_A+H_B+H_I=H_0+H_I$0, and the finite chain displays intermittent EFS rather than permanent freezing (Qian et al., 2021).

A different closed-system route freezes entanglement by converting a dynamically generated state into an exact eigenstate. In branched Y-shaped XY spin chains, an initially localized single excitation evolves into a maximally entangled symmetric endpoint state. For the minimal Y network, imposing $H=H_A+H_B+H_I=H_0+H_I$1 maps the system to a perfect-transfer chain, and at the transfer time $H=H_A+H_B+H_I=H_0+H_I$2 the endpoint state becomes

$H=H_A+H_B+H_I=H_0+H_I$3

a Bell state of the two outputs. In extended Y structures with bifurcated ends, local Pauli-$H=H_A+H_B+H_I=H_0+H_I$4 phase flips applied at the appropriate arrival time transform the symmetric four-end state

$H=H_A+H_B+H_I=H_0+H_I$5

into

$H=H_A+H_B+H_I=H_0+H_I$6

which is an eigenstate of the full Hamiltonian. After that transformation, the entanglement is exactly frozen up to a global phase (0708.2794).

These two realizations already distinguish two major EFS mechanisms in closed systems: asymptotic freezing through dispersive many-body dynamics, and exact freezing through conversion into a symmetry-protected eigenmode.

3. Eigenstate, interaction-free, and phase-space protected formulations

The interaction-free evolving framework formalizes EFS as states invisible to the interaction Hamiltonian. The existence condition is

$H=H_A+H_B+H_I=H_0+H_I$7

and pure or mixed IFE states are precisely those supported on the direct sum $H=H_A+H_B+H_I=H_0+H_I$8. For such states, there is no energy exchange between the subsystems, since

$H=H_A+H_B+H_I=H_0+H_I$9

and entanglement cannot change except under local free evolution. The non-homogeneous spin-star model provides an explicit construction in which all IFE states satisfy $|\psi\rangle$0 and may be entangled across the “central spin vs bath” bipartition (Napoli et al., 2014).

A related but control-oriented formulation appears in entanglement-storage units. There, the objective is to prepare highly entangled eigenstates of a time-independent many-body Hamiltonian, so that after control is turned off the state acquires only a global phase and its entanglement is frozen. The protocol uses the Chopped RAndom Basis method with a cost function

$|\psi\rangle$1

for the Lipkin–Meshkov–Glick model, and

$|\psi\rangle$2

for the ordered Ising chain, where $|\psi\rangle$3 is an entanglement measure and $|\psi\rangle$4 is the energy fluctuation with respect to the storage Hamiltonian $|\psi\rangle$5. ESU states are described as “gap-protected entangled eigenstates of the system Hamiltonian in the absence of the control,” and in the LMG model their bipartite entropy scales as $|\psi\rangle$6, while the lifetime measure $|\psi\rangle$7 is effectively independent of system size for the optimized storage states (Caneva et al., 2011).

A semiclassical yet still many-body realization is provided by the two-axis countertwisting model

$|\psi\rangle$8

Starting from a spin coherent state at an unstable saddle, the dynamics generate Heisenberg-scaling entanglement on the timescale

$|\psi\rangle$9

At $$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$0, a collective rotation $$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$1 moves the state to the vicinity of stable center fixed points, after which the quantum Fisher information remains on a broad plateau with

$$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$2

The plateau persists under continued nonlinear evolution, and for $$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$3-axis noise the storage is particularly robust because parity in the $$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$4 basis is conserved (Kajtoch et al., 2015).

Taken together, these constructions show that exact eigenstate protection, interaction-free evolution, gap protection, and phase-space trapping are distinct routes to the same phenomenology: entanglement that remains constant or nearly constant while the system continues to evolve.

4. Entangled frozen sectors in quantum Hilbert-space fragmentation

Recent work places EFS at the center of quantum Hilbert-space fragmentation. The starting point is a classically fragmented model defined by a semigroup word problem with local rewriting rules and a local Hamiltonian

$$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$5

If the coupling matrix $$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$6 on an equivalence class $$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$7 is rank-deficient, local null directions exist. These null directions generate exact zero-energy entangled states inside otherwise mobile classical Krylov sectors. The global frozen space is

$$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$8

where $$e^{-iHt}|\psi\rangle = e^{-i\alpha t} e^{-iH_0 t}|\psi\rangle \quad \forall t ,$$9 is the product-frozen subspace and $\alpha$0 is the direct sum of entangled-frozen sectors (Zhou et al., 6 Apr 2026).

This mechanism is demonstrated in four models of increasing algebraic structure. In the asymmetric qubit projector model with relation $\alpha$1, the local projector onto $\alpha$2 has a null vector

$\alpha$3

and each mobile classical sector contains exactly one EFS. In the $\alpha$4-symmetric GHZ projector model, the same construction survives with the null vector $\alpha$5, but the mobile quantum subspace further splits into parity sectors. In the $\alpha$6-symmetric cyclic qutrit projector, local kernel conditions such as

$\alpha$7

produce entangled-frozen subspaces whose dimensions grow with system size. In the Temperley–Lieb model with

$\alpha$8

the bond projectors satisfy the Jones relations, and after removing EFS the mobile quantum space decomposes into many irreducible standard modules $\alpha$9 (Zhou et al., 6 Apr 2026).

This leads to a weak/strong distinction. After removing EFS, if the mobile quantum Krylov subspace splits into only $|\psi\rangle$0 irreducible blocks, each showing Gaussian Orthogonal Ensemble level statistics, the model is weakly quantum fragmented; the unresolved spectrum then follows an $|\psi\rangle$1GOE distribution. If the number of irreducible blocks grows with system size and the gap-ratio distribution tends toward Poisson, the fragmentation is strong. The key point is that EFS are not peripheral anomalies in this framework; they are the mechanism by which classical fragmentation becomes quantum fragmentation (Zhou et al., 6 Apr 2026).

5. Non-Hermitian PT-symmetric frozen long-distance entanglement

A non-Hermitian realization of EFS arises in spin-$|\psi\rangle$2 SSH or XX bulks weakly coupled to two distant links. The effective Hamiltonian in the one-excitation sector is

$|\psi\rangle$3

where $|\psi\rangle$4 is the Hermitian bulk-plus-link Hamiltonian and $|\psi\rangle$5 is an alternating $|\psi\rangle$6 field generated by the no-jump evolution of a continuously monitored GKLS dynamics. The model is PT-symmetric, with exceptional points separating the PT-unbroken and PT-broken regimes (Ahuja et al., 12 Jun 2026).

For the XX bulk ($|\psi\rangle$7), a Schrieffer–Wolff treatment in the weak-coupling limit yields an effective two-link Hamiltonian

$|\psi\rangle$8

with eigenvalues $|\psi\rangle$9. The exceptional point is therefore

$H_I|\psi\rangle=\alpha|\psi\rangle$0

In the PT-broken regime near this exceptional point, an initially fully separable state evolves into a link-dominated eigenstate with nearly equal weight on the two end spins. Because non-Hermitian evolution amplifies the eigenvector with largest imaginary part, the normalized long-time state converges to a near-Bell state of the two links, producing almost frozen long-distance entanglement (Ahuja et al., 12 Jun 2026).

The link entanglement is quantified by logarithmic negativity,

$H_I|\psi\rangle=\alpha|\psi\rangle$1

and its long-time average is

$H_I|\psi\rangle=\alpha|\psi\rangle$2

Near the exceptional line and in the broken phase, $H_I|\psi\rangle=\alpha|\psi\rangle$3 approaches unity, while the temporal fluctuations are negligible; in the unbroken phase, entanglement oscillates strongly and is not frozen. The non-Hermitian system systematically outperforms the corresponding Hermitian staggered-field model, and the XX point $H_I|\psi\rangle=\alpha|\psi\rangle$4 is identified as optimal for maximizing the frozen long-distance entanglement. The same framework also supports stationary EFS: in the broken regime the stationary state is the eigenstate with largest imaginary part, and in the unbroken regime it is the ground state; both can exhibit substantial long-distance entanglement, with the non-Hermitian maxima exceeding the Hermitian ones (Ahuja et al., 12 Jun 2026).

This non-Hermitian construction extends the EFS concept beyond unitary Hamiltonian settings. Freezing is no longer produced by a stationary eigenstate of a Hermitian operator alone, but by PT-symmetric spectral selection near exceptional points.

6. Relations, misconceptions, and general significance

Several misconceptions are corrected by the literature. First, freezing is not restricted to open-system decoherence scenarios. The lossless optical-lattice study explicitly states that entanglement freezing and thawing need not arise from decoherence, but can emerge in strictly unitary multiparty dynamics (Qian et al., 2021). Second, an EFS need not be a product dark state. Interaction-free evolving states may be entangled pure or non-separable mixed states (Napoli et al., 2014), and fragmentation-based EFS are entangled zero modes embedded inside mobile sectors rather than isolated product configurations (Zhou et al., 6 Apr 2026). Third, “frozen” does not always mean “completely static state”: in some settings the state continues to evolve under local or free dynamics while the relevant entanglement measure stays invariant, whereas in others the entire state is an eigenstate and only a global phase changes (Napoli et al., 2014).

The literature also uses different entanglement diagnostics. Lossless spin-chain freezing is expressed through Schmidt weight $H_I|\psi\rangle=\alpha|\psi\rangle$5 (Qian et al., 2021); engineered Y networks use entanglement of formation and concurrence (0708.2794); TACT storage uses quantum Fisher information as the operative entanglement and metrological monotone (Kajtoch et al., 2015); non-Hermitian long-distance EFS use logarithmic negativity (Ahuja et al., 12 Jun 2026). This suggests that EFS are best understood as a dynamical or structural property, not as a feature tied to a single entanglement measure.

A plausible synthesis is that EFS typically require one of a small number of organizing mechanisms: spectral isolation by large gaps, destructive interference or local null directions, confinement near stable classical structures, or spectral selection near exceptional points. The exact implementation differs sharply across models, but in every case the entanglement remains stable because the relevant state or subspace is dynamically decoupled from the channels that would otherwise redistribute amplitude. The result is a family of phenomena linking many-body dynamics, constrained quantum systems, non-Hermitian physics, and quantum information storage.

In that broader sense, Entangled Frozen States designate not one unique state class but a research program: the identification of entangled states or entanglement regimes whose persistence is enforced by Hamiltonian structure, symmetry, locality, algebraic fragmentation, optimal control, or PT-symmetric spectral physics.