1 Introduction

The study of nonequilibrium quantum dynamics in many-body systems has been a longstanding pursuit in contemporary condensed matter physics. The unitary time evolution of a generic nonintegrable quantum many-body system would approach the thermal equilibrium, a phenomenon closely associated with the eigenstate thermalization hypothesis (ETH) [1,2,3,4]. In recent years, extensive studies have been exploring quantum systems that violate the ETH. Notably, the many-body localization provides an interesting possibility of ETH violation by introducing disorders [5,6,7,8,9,10,11]. More recently, the ETH violation due to quantum many-body scar states and Hilbert space fragmentation has also greatly attracted both theoretical and experimental investigations [12,13,14,15,16,17,18,19,20,21].

Lattice gauge theory provides an alternative approach to ETH violation systems, leading to a wide class of dynamical phenomena associated with the configurations of gauge fields. In particular, the gauge degrees of freedom may induce disorder-free localization in quantum systems [22,23,24,25,26,27,28,29,30]. Besides, the lattice gauge theory can also hold quantum many-body scar states embedded in thermal eigenstates [31,32,33,34,35,36,37].

Recently, an intriguing quantum many-body system, called the quantum breakdown model, was proposed to describe the dielectric breakdown process from a microscopic perspective [38,39,40,41]. The one-dimensional (1D) fermionic quantum breakdown model features a spatially asymmetric breakdown interaction that annihilates a fermion at one site and simultaneously creates more fermions at the neighboring site on the right [38]. With an increasing number of fermion orbitals (flavors) at each site, the quantum breakdown model undergoes a crossover from the many-body localization phase to the quantum chaotic phase with scar states.

In this paper, we investigate the minimal quantum breakdown model, which necessitates the smallest number of fermion orbitals that are compatible with the spatially asymmetric breakdown interactions. Under this circumstance, we can identify an extensive number of local conserved quantities that contribute to the fragmentation of the total Hilbert space. Notably, we find that such a minimal quantum breakdown model in certain charge sectors can be mapped to a model of lattice gauge theory, known as the quantum link model [42,43,44,45]. The latter can be experimentally simulated in various quantum devices [46,47,48,49,50,51,52,53,54,55,56,57,58,59]. Through this mapping, the local conserved quantities in the quantum breakdown model play the role of the gauge symmetry generators in the lattice gauge theory. As a result, Krylov subspaces with distinct gauge configurations give rise to various subspace dynamics, ranging from free fermions on hypercubic lattices with boundary defects to strongly interacting sectors with quantum many-body scars. Our results reveal that the lattice gauge theory not only offers a theoretical perspective to understand the dynamics in the quantum breakdown model but also provides a practical approach to simulating this model in advanced quantum experiments.

The rest of this paper is organized as follows. In Sect. 2, we introduce the Hamiltonian of the quantum breakdown model and identify its symmetry and conserved quantities. Then in Sect. 3, we map the minimal quantum breakdown model to U(1) lattice gauge theory. Based on this mapping, we discuss various quantum dynamical behaviors in certain representative gauge sectors in Sect. 4. Our work is then concluded in Sect. 5.

2 Quantum breakdown model

2.1 Model Hamiltonian

The breakdown process of a dielectric gas subjected to a sufficiently strong electric field can be phenomenologically described as follows. Because of the strong electric field, the neutral atom can be ionized into one electron and one ion. Then, the free electron is immediately accelerated by the strong electric field. On the contrary, the produced ion is accelerated in the opposite direction, but experiences much slower dynamics because of its much heavier mass. Therefore, we ignore the ion dynamics and focus only on the fast motions of electrons. Subsequently, the fast electrons collide with other atoms, triggering the progressive generation of additional electrons and ions. As a result, more and more electrons are generated and accelerated by the electric field, leading to a Townsend particle avalanche of electrons [60]

By ignoring the ions, this breakdown process can be effectively described by a microscopic Hamiltonian called the quantum breakdown model [38]. We consider a 1D system with M sites, each site having N fermion orbitals (flavors). The generic Hamiltonian is given by

\(\begin{aligned} \hat{H}=\hat{H}_I+\hat{H}_\mu . \end{aligned}\)
(1)

With \(\hat{c}_{m,i}^\dagger\) and \(\hat{c}_{m,i}\) being the creation and annihilation operators of the ith fermionic mode at the mth site, the interacting part \(\hat{H}_I\) represents the spatially asymmetric breakdown interaction

\(\begin{aligned} \hat{H}_I=\sum \limits _{m=1}^{M-1}\sum \limits _{l=1}^{N}\sum \limits _{i_1<\cdots
(2)

Here, \(\text {h.c.}\) is the Hermitian conjugate, and q is a nonnegative integer. The values of the interaction strength \(J_{m, l}^{i_1 i_2\cdots i_{2q+1}}\) are complex numbers that are arranged antisymmetrically with respect to the indices \(i_1,\cdots ,i_{2q+1}\). This asymmetric interaction indicates that the annihilation of one fermion leads to the creation of \(2q+1\) fermions at the adjacent site. Therefore, Eq. (2) defines a class of quantum breakdown models with different q values. The asymmetrical interaction is defined to maintain the fermion parity. For the nontrivial \(\hat{H}_I\) to be valid, the number of fermion orbitals per site must satisfy

\(\begin{aligned} N\ge 2q+1\ . \end{aligned}\)
(3)

The second part \(\hat{H}_\mu\) is the on-site potential, which is given by

\(\begin{aligned} \hat{H}_\mu =\sum \limits _{m=1}^M\mu _m \hat{n}_{m},\quad \hat{n}_m=\sum \limits _{i=1}^N\hat{c}_{m,i}^\dagger \hat{c}_{m,i}. \end{aligned}\)
(4)

Here, \(\mu _m\) represents the potential at the mth site. Also, \(\hat{n}_m\) is the fermion number operator at the mth site.

The quantum breakdown model displays a wide variety of dynamical phases, including many-body localization, Hilbert space fragmentation, and quantum chaos [38,39,40,41]. As shown in Ref. [38], the \(q=1\) quantum breakdown model has almost all eigenstates solvable when \(N=3\), while it exhibits quantum chaos with many-body scar states when N is large. Interestingly, a dynamical breakdown transition is controlled by the ratio between the interaction strength and the on-site potential. A considerably large interaction is required to overcome the energy barrier induced by the on-site potential, leading to the proliferation of electrons.

In this paper, we study the minimal quantum breakdown models which are defined by the requirement

\(\begin{aligned} N=2q+1\ . \end{aligned}\)
(5)

In contrast to the quantum chaotic phase established in the large N regime, we will show that the \(N=2q+1\) case has a large number of conserved quantities that make this model almost exactly solvable.

2.2 Symmetry and conserved quantity

To obtain the conserved quantities of the minimal quantum breakdown model, we need to analyze its symmetry. In this section, we focus on the global symmetries of this model. We first consider the following spatially dependent unitary transformation:

\(\begin{aligned} \hat{V}_m\hat{c}_{m,i} \hat{V}_m^\dagger =e^{i\varphi _m}\hat{c}_{m,i},\quad \hat{V}_m=e^{-i\varphi _m\hat{n}_m}. \end{aligned}\)
(6)

Then, the invariance of the breakdown Hamiltonian requires that

\(\begin{aligned} \varphi _m=(2q+1)\varphi _{m+1} \mod 2\pi . \end{aligned}\)
(7)

Moreover, as we will show below, the boundary conditions have a significant effect on this equation, as the periodic boundary condition (PBC) necessitates an extra restriction between the first site and the last site, which is not present under the open boundary condition (OBC).

We first discuss the symmetry of the PBC. In this case, the phase relations of \(\hat{V}_m\) that keep the Hamiltonian invariant are given by Eq. (7) for \(m=1,2,\dots ,(M-1)\). Assuming \(\varphi _M=\varphi\) where \(\varphi\) is a site-independent constant angle, these relations immediately lead to a close solution \(\varphi _m=(2q+1)^{M-m}\varphi\). Furthermore, the PBC further imposes a constraint \(\varphi _M=(2q+1)\varphi _1 \mod 2\pi\), which requires \(\varphi\) to satisfy the following condition: \(\varphi =(2q+1)^M\varphi \mod 2\pi\). This condition implies that the angle \(\varphi\) can only take discrete values

\(\begin{aligned} \varphi =\frac{2\pi p}{(2q+1)^M-1}\ ,\ \ \ (p=0,1,\dots ,(2q+1)^M-2) \end{aligned}\)
(8)

Therefore, the quantum breakdown model under the PBC has a discrete \(\mathbb {Z}_{(2q+1)^M-1}\) symmetry, a global symmetry depending on the system size [61,62,63].

For the OBC, the phase relations in Eq. (7) give rise to a solution \(\varphi _m=(2q+1)^{M-m}\varphi\) with \(\varphi \in [0,2\pi )\) taking continuous values. Therefore, the quantum breakdown model with OBC has a spatially modulated global U(1) symmetry called the exponential symmetry [38, 40, 61,62,63,64,65,66]. This exponential symmetry is generated by an exponential U(1) charge

\(\begin{aligned} \hat{Q}=\sum \limits _{m=1}^M(2q+1)^{M-m}\hat{n}_m. \end{aligned}\)
(9)

The conserved charge \(\hat{Q}\) implies that the fermions at the mth site have an effective charge \(q_m=(2q+1)^{M-m}\). Intuitively, the asymmetric breakdown interactions annihilate one fermion at the mth site and create \(2q+1\) fermions at the adjoint site, splitting the effective charge into \(2q+1\) pieces. This conserved charge makes it possible to use exact diagonalization to study the energy spectrum and quantum dynamics within each charge sector [38].

2.3 Extensive conserved quantities in the minimal quantum breakdown model

While the symmetry analysis in Sect. 2.2 applies to the quantum breakdown model with \(N\ge 2q+1\), the minimal quantum breakdown model with \(N=2q+1\) (Eq. (5)) has a richer and more interesting structure which we will focus on in the rest of the paper. As we shall show below, an extensive number of local conserved quantities exist in the \(N=2q+1\) quantum breakdown model. These conserved quantities result in exponentially many disconnected Krylov subspaces, a hallmark of Hilbert space fragmentation. For simplicity, we impose OBC here, for which case we do not need to worry about the relation between the first site and the last site. The analysis for symmetries under PBC and the corresponding dynamical properties are however similar.

To extract the conserved quantities, we express the minimal quantum breakdown model in a simpler form. Specifically, we employ a local \(U(2q+1)\) unitary transformation among the \(2q+1\) fermion flavors on the mth site as

\(\begin{aligned} \hat{c}_{m,i}=\sum \limits _{j}U_{i,j}^{(m)}\hat{f}_{m,j}, \end{aligned}\)
(10)

where \(U^{(m)}\) is an element in the \(U(2q+1)\) group. On the one hand, the uniform on-site potential \(\hat{H}_\mu\) in Eq. (4) is invariant under this transformation, since \(\hat{n}_m=\sum \nolimits _{i=1}^{2q+1}\hat{c}_{m,i}^\dagger \hat{c}_{m,i}=\sum \nolimits _{i=1}^{2q+1} \hat{f}_{m,i}^\dagger \hat{f}_{m,i}\). On the other hand, since there are only \(2q+1\) fermion modes per site, we have \(\hat{c}_{m,1}\hat{c}_{m,2}\cdots \hat{c}_{m,2q+1}= \det [U^{(m)}]\hat{f}_{m,1}\hat{f}_{m,2}\cdots \hat{f}_{m,2q+1}\) and \(\hat{c}_{m,1}^\dagger \hat{c}_{m,2}^\dagger \cdots \hat{c}_{m,2q+1}^\dagger = \det [U^{(m)}]^*\hat{f}_{m,1}^\dagger \hat{f}_{m,2}^\dagger \cdots \hat{f}_{m,2q+1}^\dagger\), namely, \(c_{m,1}\hat{c}_{m,2}\cdots \hat{c}_{m,2q+1}\) and its Hermitian conjugate transform as a singlet in the on-site flavor space. Then, the breakdown interaction \(\hat{H}_I\) transforms as

\(\begin{aligned} \hat{H}_I=\sum \limits _{m=1}^{M-1}\sum \limits _{l,l'=1}^{2q+1}J_{m,l}^{12\cdots (2q+1)}[\det U^{(m+1)}]^*U_{l,l'}^{(m)}\left( \prod \limits _{k=1}^{2q+1}\hat{f}_{m+1, k}^{\dagger }\right) \hat{f}_{m,l'}+\text {h.c.} \end{aligned}\)
(11)

Since the coefficients \(J_{m,l}^{12\cdots (2q+1)}\) is a vector with component l under the \(U(2q+1)\) rotation, we can always choose the unitary transformation \(U^{(m)}\) such that [38]

\(\begin{aligned} \sum \limits _{l=1}^{2q+1} J_{m,l}^{12\cdots (2q+1)}[\det U^{(m+1)}]^*U_{l,l'}^{(m)}=J_m\delta _{l',1}. \end{aligned}\)
(12)

The real number \(J_m\) is the norm of the vector \((J_{m,1}^{12\cdots (2q+1)},J_{m,2}^{12\cdots (2q+1)},\cdots ,J_{m,2q+1}^{12\cdots (2q+1)})^T\), which is given by

\(\begin{aligned} J_m=\sqrt{\sum \limits _{l=1}^{2q+1}\left| J_{m,l}^{12\cdots (2q+1)}\right| ^2}\ . \end{aligned}\)
(13)

Then, the minimal quantum breakdown model takes a much simpler form:

\(\begin{aligned} \hat{H} & =\hat{H}_I+\hat{H}_\mu ,\nonumber \\ \hat{H}_I & = \sum \limits _{m=1}^{M-1}\left[ J_m\left( \prod \limits _{i=1}^{2q+1}\hat{f}_{m+1, i}^{\dagger }\right) \hat{f}_{m,1}+ \text{ h.c. } \right] ,\nonumber \\ \hat{H}_{\mu } & = \sum \limits _{m=1}^M \mu _m \hat{n}_m. \end{aligned}\)
(14)

After the basis transformation, only fermions on the first orbital can move between different sites. A fermion moving in the right direction generates additional 2q fermions at the right adjacent site (Fig. 1b). Once these 2q fermions are created, they become immobile, incapable of further moving rightward. The only possible dynamics for these 2q fermions is their simultaneous annihilation when a fermion on the first orbital moves to the left site, namely, the Hermitian conjugate of their creation process. We note that such a simplification of Hamiltonian similar to Eq. (14) via local unitary transformations is not applicable to generic quantum breakdown models with \(N>2q+1\).

Fig. 1
figure 1

a The \(q=2\) quantum breakdown model. b The quantum breakdown model after the basis rotation. A solid circle (\(\bullet\)) denotes an occupied orbital, while a hollow circle (\(\circ\)) indicates an unoccupied orbital


The simplified form of the minimal quantum breakdown model in Eq. (14) allows us to reveal many more hidden conserved charges. To see this, we attach an effective charge, denoted as \(q_{m,i}\), to the ith orbital f-fermion on the mth site. Consequently, it is straightforward to see that the following modulated charge \(\hat{Q}(\{q_{m,i}\})\) commutes with the Hamiltonian \(\hat{H}\) in Eq. (14) provided that \(q_{m,i}\) satisfy the following condition:

\(\begin{aligned} \hat{Q}(\{q_{m,i}\})=\sum \limits _{m=1}^M\sum \limits _{i=1}^{2q+1}q_{m,i}f_{m,i}^\dagger f_{m,i}\ ,\quad q_{m,1} = \sum \limits _{i=1}^{2q+1} q_{m+1,i}\ . \end{aligned}\)
(15)

The above charge \(\hat{Q}(\{q_{m,i}\})\) in Eq. (15) reduces to the conserved charge \(\hat{Q}\) for c-fermions in Eq. (9) if one chooses \(q_{m,i} = (2q+1)^{M-m}\). Clearly, the arbitrariness of \(q_{m,i}\) in Eq. (15) gives rise to many more conserved quantities.

To further extract the conserved quantities encoded in Eq. (15), we can reformulate the charge constraint as \(\sum \nolimits _{i=1}^{2q+1} q_{m+1,i}/q_{m,1} = 1\). In particular, we choose the ratios \(q_{m+1,i}/q_{m,1}\) to be given by the following parameters:

\(\begin{aligned} q_{m+1,1}/q_{m,1}&=\gamma ,\nonumber \\ q_{m+1,2}/q_{m,1} & =(1-\gamma )/(2q)+\beta _2,\nonumber \\ q_{m+1,3}/q_{m,1} & =(1-\gamma )/(2q)+\beta _3-\beta _2,\nonumber \\ \cdots & \nonumber \\ q_{m+1,i}/q_{m,1} & =(1-\gamma )/(2q)+\beta _{i}-\beta _{i-1},\nonumber \\ \cdots & \nonumber \\ q_{m+1,2q}/q_{m,1} & =(1-\gamma )/(2q)+\beta _{2q}-\beta _{2q-1},\nonumber \\ q_{m+1,2q+1}/q_{m,1} & =(1-\gamma )/(2q)-\beta _{2q}. \end{aligned}\)
(16)

Here, \(\gamma\) and \(\beta _i\) with \(i=2,3,\dots ,2q\) represent the specific 2q free parameters (which can be any complex numbers). Consequently, the conserved quantity in Eq. (15) transforms into:

\(\begin{aligned} \hat{Q}(\{q_{m,i}\}) & =\sum \limits _{i=2}^{2q+1}q_{1,i}\hat{f}_{1,i}^\dagger \hat{f}_{1,i}+q_{1,1}\sum \limits _{m=1}^M\gamma ^{m-1} \hat{f}_{m,1}^\dagger \hat{f}_{m,1}\nonumber \\ & \quad +q_{1,1}\sum \limits _{m=2}^M\gamma ^{m-2}\frac{1-\gamma }{2q}\sum \limits _{i=2}^{2q+1}\hat{f}_{m,i}^\dagger \hat{f}_{m,i}\nonumber \\ & \quad +q_{1,1}\sum \limits _{m=2}^M\gamma ^{m-2} \sum \limits _{i=2}^{2q}\beta _i(\hat{f}_{m,i}^\dagger \hat{f}_{m,i}-\hat{f}_{m,i+1}^\dagger \hat{f}_{m,i+1}). \end{aligned}\)
(17)

\(q_{1,i}\) for \(i=1,\cdots ,2q+1\) represent the effective charges for fermions at the first site. The conservation of the total charge is independent of the choices of the free parameters \(\gamma\) and \(\beta _i\); thus, the coefficient of each power of these free parameters can be identified as an independent conserved quantity. All of these coefficients give a large number of local conserved quantities. Notably, these local conserved quantities can be categorized into three distinct sets as follows.

The first set of conserved quantities is localized on the first site and is denoted as:

\(\begin{aligned} \hat{Q}_{a,i} = \hat{f}_{1,i}^\dagger~\hat{f}_{1,i}, \quad i = 2, \cdots , 2q+1. \end{aligned}\)
(18)

This set signifies that the fermions on 2q orbitals of the first site are entirely decoupled from the rest of the system. These 2q conserved quantities at the first site directly stem from taking the OBC.

The second set comprises on-site conserved quantities as:

\(\begin{aligned} \hat{Q}_{b,m,i} = \hat{f}_{m,i}^\dagger \hat{f}_{m,i} - \hat{f}_{m,i+1}^\dagger \hat{f}_{m,i+1}, \end{aligned}\)
(19)

where the site index \(m = 2,\dots ,M\) and orbital index \(i = 2,\dots ,2q\). These on-site conserved quantities indicate that the \(\hat{f}_{m,i>1}\) fermions at the mth site are subject to simultaneous annihilation and creation, revealing the conservation of the population imbalance between f-fermion orbitals with orbital indices \(i>1\). This behavior aligns directly with the breakdown interaction in Eq. (14).

The third set of conserved quantities is also local, but intriguingly not on-site, which is given by:

\(\begin{aligned} \hat{Q}_{c,m} = \hat{f}_{m,1}^\dagger \hat{f}_{m,1} + \frac{1}{2q}\sum \limits _{i=2}^{2q+1}\left( \eta _m\hat{f}_{m+1,i}^\dagger \hat{f}_{m+1,i} - \tilde{\eta }_m\hat{f}_{m,i}^\dagger \hat{f}_{m,i}\right) , \end{aligned}\)
(20)

where \(m = 1,\dots ,M\). These conserved quantities describe the interactions between the first fermion \(\hat{f}_{m,1}\) at each site and the other fermions. We impose the coefficients \(\eta _m = 1-\delta _{m,M}\) and \(\tilde{\eta }_m = 1-\delta _{m,1}\) such that these conserved quantities are compatible with the OBC.

As a result, we obtain an extensive number of local conserved operators for the minimal quantum breakdown Hamiltonian in Eq. (14). These operators in Eqs. (18) to (20) can further add or multiply to generate additional conserved quantities, serving as generators of the underlying commutant algebra [21]. For example, the total fermion number \(\hat{\mathcal {N}}_1 = \sum \nolimits _{m=1}^M \hat{f}_{m,1}^\dagger \hat{f}_{m,1}\) on the first orbital of all sites is conserved, which is equivalent to \(\hat{\mathcal {N}}_1 = \sum \nolimits _{m=1}^M \hat{Q}_{c,m}\).

As implied by the conserved quantities \(\hat{Q}_{a,i}\), the fermions created by \(\hat{f}_{1,i}^\dagger\) with \(i>1\) remain frozen in their dynamics. However, the remaining two groups of conserved quantities, \(\hat{Q}_{b,m,i}\) and \(\hat{Q}_{c,m}\), play a crucial role in shaping the connected Hilbert subspaces (i.e., Krylov subspaces) and, thereby, constraining the quantum dynamics of the minimal quantum breakdown model with \(N=2q+1\).

The breakdown interactions in Eq. (14) simultaneously annihilate or create all fermions with orbital index \(i>1\). Consequently, the subspace that exhibits nontrivial dynamics at the mth site must exclusively comprise states with an identical number of fermions with orbital index \(i>1\). Conversely, a site with an uneven distribution of these fermion modes becomes dynamically frozen. These two distinct types of states can be distinguished by the eigenvalues of the second set of conserved charges \(\hat{Q}_{b,m,i}\). In a subspace where any \(\hat{Q}_{b, m, i}\) has a nonzero eigenvalue, the fermions at the mth site with orbital index \(i>1\) will remain static. As a result, this site serves as a blocking site [38], effectively dividing the system into two dynamically isolated regions. This dynamical constraint precisely exemplifies the Hilbert space fragmentation in the \(N=2q+1\) quantum breakdown model.

Assume that there are two blocking sites at \(m_1\) and \(m_2\) and no other blocking site in between. The nontrivial dynamics then exist in the region \(m_1, in which any \(\hat{Q}_{b, m, i}\) with \(m_1+1\le m\le m_2-1\) and \(2\le i\le 2q+1\) would have a zero eigenvalue for the states in this subspace. Within such a subspace, the quantum dynamics is further influenced by the conserved charges \(\hat{Q}_{c,m}\). In the subsequent section, we will demonstrate that the minimal quantum breakdown model in certain subspaces is equivalent to a U(1) lattice gauge model. In the language of lattice gauge theory, the conserved quantities \(\hat{Q}_{c,m}\) play the role of gauge symmetry generators.

3 Mapping to the lattice gauge theory

Without loss of generality, we make a slight change of notation in the following discussion. Here, we assume \(m_1=1\) and \(m_2=M+1\) are two blocking sites, between which the region is connected without any other blocking sites. As explained in the above section, the states \(|{\psi }\rangle\) in this dynamical subspace are restricted by

\(\begin{aligned} \hat{Q}_{b,m,i}|{\psi }\rangle =0, \quad 2\le m\le M,\ i=2,\cdots ,2q+1. \end{aligned}\)
(21)

To further simplify the minimal quantum breakdown model in Eq. (14), we define a set of operators as

\(\begin{aligned} \hat{F}_m^\dagger & = \hat{f}_{m,2}^\dagger \hat{f}_{m,3}^\dagger \cdots \hat{f}_{m,2q}^\dagger \hat{f}_{m,2q+1}^\dagger ,\nonumber \\ \hat{F}_m & = \hat{f}_{m,2q+1}\hat{f}_{m,2q}\cdots \hat{f}_{m,3}\hat{f}_{m,2},\nonumber \\ \hat{N}_m & = \frac{1}{2q}\sum \limits _{i=2}^{2q+1}\hat{f}_{m,i}^\dagger \hat{f}_{m,i}. \end{aligned}\)
(22)

\(\hat{F}_m\) and \(\hat{F}_m^\dagger\) are the operators for collectively annihilating and creating 2q fermions at the mth site. Then, the quantum breakdown model in Eq.(14) can be re-expressed as

\(\begin{aligned} \hat{H}=\sum \limits _{m=1}^{M-1}( J_m~\hat{f}_{m+1,1}^\dagger \hat{F}_{m+1}^\dagger \hat{f}~_{m,1}+ \text{ h.c.})+\sum \limits _{m=1}^M \mu _m~(\hat{f}_{m,1}^\dagger \hat{f}~_{m,1}+2q\hat{N}_m). \end{aligned}\)
(23)

Since the breakdown interaction preserves the fermion parity, the two operators, \(\hat{F}_m\) and \(\hat{F}_m^\dagger\), behave like bosonic operators. They satisfy the following commutation relations:

\(\begin{aligned} \left[ \hat{F}_m,\hat{f}_{n,1}\right] & = \left[ \hat{F}_m^\dagger ,\hat{f}_{n,1}^\dagger \right] =0, \nonumber \\ \left[ \hat{F}_m,\hat{F}_n\right] & = \left[ \hat{F}_m^\dagger ,\hat{F}_n^\dagger \right] =0, \nonumber \\ \left[ \hat{F}_n^\dagger ,\hat{F}_m\right] & = \delta _{mn}\left( \prod \limits _{i=2}^{2q+1}\hat{f}_{m,i}^\dagger \hat{f}_{m,i}-\prod \limits _{i=2}^{2q+1}(1-\hat{f}_{m,i}^\dagger \hat{f}_{m,i})\right) . \end{aligned}\)
(24)

Upon imposing the constraint in Eq. (21) on the Hilbert space, only two configurations for the 2q immobile f-fermion modes (with orbital indices \(2\le i\le 2q+1\)) at each site are accessible: the vacuum state \(|{\downarrow }\rangle _m\equiv |{0}\rangle\) and the fully occupied state \(|{\uparrow }\rangle _m\equiv \hat{F}_m^\dagger |{0}\rangle\) (Fig. 2a). Consequently, it is straightforward to show that \([\hat{F}_m^\dagger ,\hat{F}_m]|{\downarrow }\rangle _m=-|{\downarrow }\rangle _m\) and \([\hat{F}_m^\dagger ,\hat{F}_m]|{\uparrow }\rangle _m=|{\uparrow }\rangle _m\). Therefore, the commutator \([\hat{F}_m^\dagger ,\hat{F}_m]\) can be expressed as \([\hat{F}_m^\dagger ,\hat{F}_m]_c= 2\hat{N}_m-1\), where the subscript c denotes the subspace spanned by \(|{\uparrow }\rangle _m\) and \(|{\downarrow }\rangle _m\). Furthermore, \([\hat{F}_m,2\hat{N}_m-1]_c=2\hat{F}_m\) and \([\hat{F}_m^\dagger ,2\hat{N}_m-1]_c=-2\hat{F}_m^\dagger\). As a result, the operators \(\hat{F}_m\), \(\hat{F}_m^\dagger\), and \(2\hat{N}_m-1\) constrained in this subspace follow the same algebra as the Pauli matrices \(\hat{\sigma }_{m,m+1}^{\pm ,z}\), with the correspondence

\(\begin{aligned} \hat{F}_{m+1}\rightarrow \hat{\sigma }_{m,m+1}^-\ ,\ \hat{F}_{m+1}^\dagger \rightarrow \hat{\sigma }_{m,m+1}^+\ ,\ 2\hat{N}_{m+1}-1\rightarrow \hat{\sigma }_{m,m+1}^z\ . \end{aligned}\)
(25)

Fig. 2
figure 2

a The correspondence between fermions and link spins established in Eq. (25). b The sketch of interactions in the minimal quantum breakdown model in Eq. (26). c A typical configuration in the minimal quantum breakdown model in Eq. (26)


We emphasize that such an identification between fermionic operators and spin operators is exclusively applicable within the constrained Hilbert space where Eq. (21) is satisfied.

After mapping to spin operators, we can express the minimal quantum breakdown model in the subspace of Eq. (21) into the following form:

\(\begin{aligned} \hat{H} & = \sum \limits _{m=1}^{M-1}( J_m\hat{f}_{m+1,1}^\dagger \hat{\sigma }_{m,m+1~}^+ \hat{f}_{m,1}+\text{ h.c.}) \nonumber \\ & \quad +\sum \limits _{m=1}^M \mu _m~\hat{f}_{m,1}^\dagger \hat{f}_{m,1}+\sum \limits _{m=1}^{M-1}q\mu _m~\hat{\sigma }_{m,m+1}^z. \end{aligned}\)
(26)

Here, we ignore the decoupled fermions on the first site with orbital indices \(2\le i\le 2q+1\) and omit some constant terms. The Hamiltonian, as shown in Fig. 2, resembles the quantum link model, which characterizes the motion of fermions along the lattice sites while simultaneously interacting with the spins on the lattice links [42,43,44,45]. This is a lattice version of 1+1D quantum electrodynamics. Correspondingly, the local conserved quantities \(\hat{Q}_{c,m}\) in Eq. (20) are transformed into

\(\begin{aligned} \hat{G}_m=\hat{f}_{m,1}^\dagger \hat{f}_{m,1}+\frac{1}{2}[\eta _m(\hat{\sigma }_{m,m+1}^z+1)-\tilde{\eta }_m(\hat{\sigma }_{m-1,m}^z+1)]. \end{aligned}\)
(27)

The quantum link model bridges the minimal quantum breakdown model with the 1D lattice gauge theory. To see this, we introduce the Hamiltonian of the U(1) lattice gauge theory [45, 46]. Here, we focus on the open boundary conditions, and the quantum link formulation of the U(1) lattice gauge theory is

\(\begin{aligned} \hat{H}_{\text {LGT}} & = -t\sum \limits _{m=1}^{M-1}(\hat{\Psi }_{m}^\dagger \hat{U}_{m,m+1} \hat{\Psi }_{m+1}+\text {h.c.})\nonumber \\ & \quad +\mu \sum \limits _{m=1}^M(-1)^m \hat{\Psi }_m^\dagger \hat{\Psi }_m+g\sum \limits _{m=1}^{M-1}\hat{E}_{m,m+1}^2 . \end{aligned}\)
(28)

In this context, \(\hat{\Psi }_m\) and \(\hat{\Psi }_m^\dagger\) denote the fermionic annihilation and creation operators on the lattice sites. Meanwhile, \(\hat{E}_{m,m+1}\) represents the electric field operator on the lattice links, and \(\hat{U}_{m,m+1}=e^{i\hat{A}_{m,m+1}}\) corresponds to a parallel transport operator induced by the U(1) link gauge field \(\hat{A}_{m,m+1}\). They adhere to the relation \([\hat{E}_{m,m+1}, \hat{U}_{m,m+1}]=\hat{U}_{m,m+1}\). This relationship arises from the fact that the electric field operator \(\hat{E}_{m,m+1}\) on the link is the canonical momentum of the link gauge field \(\hat{A}_{m,m+1}\).

The first term in Eq. (28) signifies the couplings between the fermions (matter fields) at the lattice sites and the gauge fields at the lattice links. The second term represents the fermion mass, while the final term characterizes the energy of electric fields. Notably, the U(1) gauge symmetry of this lattice model under open boundary conditions is generated by [45, 46]:

\(\begin{aligned} \hat{\mathcal {G}}_m=\hat{\Psi }_m^\dagger \hat{\Psi }_m+\tilde{\eta }_m\hat{E}_{m-1,m}- \eta _m\hat{E}_{m,m+1}. \end{aligned}\)
(29)

These gauge generators satisfy \([\hat{H}_{\text {LGT}},\hat{\mathcal {G}}_m]=0\). In the lattice gauge theory for staggered fermions, the physical states in the gauge-invariant subspace satisfy \((\hat{\mathcal {G}}_m-\frac{1-(-1)^m}{2})|{\psi }\rangle =0\), which is a lattice manifestation of Gauss law \(\triangledown \cdot E=\rho\) [45,46,47].

The infinite local Hilbert space dimension associated with the gauge fields on the links poses a challenge in simulating the U(1) lattice gauge theory in experiments. A truncated local Hilbert space with small \(\hat{E}_{m,m+1}^2\) eigenvalues is often adopted, which is more experimentally accessible and is justified for large coupling strength g in \(\hat{H}_{\text {LGT}}\). This truncation leads to the quantum link model, achieved by the substitutions: \(\hat{U}_{m,m+1}\rightarrow \hat{S}_{m,m+1}^+\), \(\hat{U}_{m,m+1}^\dagger \rightarrow \hat{S}_{m,m+1}^-\), and \(\hat{E}_{m,m+1}\rightarrow \hat{S}_{m,m+1}^z\), where \(S_{m,m+1}^\mu\) are spin-S operators on the link. Now, the gauge fields on the links are effectively substituted by a spin S. Consequently, the spin-S quantum link model can be expressed as [45, 46] :

\(\begin{aligned} \hat{H}_{S} & = -t^\prime \sum \limits _{m=1}^{M-1}(\hat{\Psi }_{m}^\dagger \hat{S}_{m,m+1} ^+ \hat{\Psi }_{m+1}+\text {h.c.})\nonumber \\ & \quad +\mu \sum \limits _{m=1}^M(-1)^m \hat{\Psi }_m^\dagger \hat{\Psi }_m+g^\prime \sum \limits _{m=1}^{M-1}(\hat{S}_{m,m+1}^z)^2. \end{aligned}\)
(30)

The precise form of the renormalized constants \(t^\prime\) and \(g^\prime\) is irrelevant in our discussion. In the case of the spin \(S=\frac{1}{2}\) quantum link model, the energy term of electric fields is a constant energy offset that can hence be dropped. When we map the \(\Psi\)-fermions to the f-fermions, the model in Eq. (30) shares the same interaction terms as in Eq. (26) by the following substitutions: \(\hat{S}_{m,m+1}^{\pm } \rightarrow \hat{\sigma }_{m,m+1}^{\mp }\) and \(2\hat{S}_{m,m+1}^z \rightarrow -\hat{\sigma }_{m,m+1}^z\). It is worth noting that we have chosen a different basis for the link spins in Eq. (26) to establish a connection between the spin-up state in the quantum link model and the occupied state \(\hat{F}_m^\dagger |{0}\rangle\) in the minimal quantum breakdown model.

By employing the above substitutions and setting the parameters \(J_m=-t^\prime\) and \(\mu _m=(-1)^m\mu\), the minimal quantum breakdown model in Eq. (26) can be mapped into the \(S=\frac{1}{2}\) version of the quantum link model in Eq. (30). However, this mapping is not exact. The distinction arises from the staggered chemical potential in the minimal quantum breakdown model, which results in a staggered magnetic field for link spins in Eq. (26). The staggered magnetic field for link spins is different from the energy term of electric fields \(g^\prime (\hat{S}_{m,m+1}^z)^2\) presented in Eq.(30). An interesting future direction is to investigate the effect of this effective magnetic field on the quantum dynamics of the underlying many-body system. Nevertheless, we emphasize that the minimal quantum breakdown model shares the same gauge symmetry structure as the U(1) lattice gauge theory, which can be seen in Eqs. (20), (27), and (29).

In general, quantum simulations of the lattice gauge theory make great efforts to enforce the gauge condition \((\hat{\mathcal {G}}_m-\frac{1-(-1)^m}{2})|{\psi }\rangle =0\) in experiments. However, the connection between lattice gauge theory and the quantum breakdown model we revealed here suggests that the other gauge sectors may also unveil intriguing physical phenomena. In the next section, we focus on the Hamiltonian Eq. (26) and delve into several typical gauge sectors, to illustrate the rich dynamical phenomena exhibited by the minimal quantum breakdown model.

4 Dynamically connected subspaces

In the \(N=2q+1\) minimal quantum breakdown model, the existence of local conserved quantities effectively partitions the entire Hilbert space into an exponential number of disconnected Krylov subspaces. Even after resolving the conserved quantities \(\hat{Q}_{a,i}\) and \(\hat{Q}_{b,m,i}\), it is still possible to make further fragmentations within the Hilbert subspace constrained by \(\hat{Q}_{b,m,i}|{\psi }\rangle =0\). As a result, the conserved quantities \(\hat{Q}_{c,m}\) play an indispensable role in determining the dynamically connected subspaces. To make it clear, we employ the minimal quantum breakdown model in the form of Eq. (26) to illustrate the dynamical structures in the original quantum breakdown model.

4.1 Blocking gauge configuration

In addition to the aforementioned blocking sites with at least one \(\hat{Q}_{b,m,i}\) satisfying \(\hat{Q}_{b,m,i}|{\psi }\rangle \ne 0\), there also exist two types of blocking configurations in the subspaces with all \(\hat{Q}_{b,m,i}|{\psi }\rangle = 0\). For simplicity, we take the form of the minimal quantum breakdown model shown in Eq. (26). Such blocking configurations are determined by the local gauge generators \(\hat{G}_m\) (i.e., \(\hat{Q}_{c,m}\)) in Eq. (27), which can have eigenvalues \(-1,0,1,2\) (if \(m\ne 1\) or M). As shown in Fig. 3, when the eigenvalue of \(\hat{G}_m\) at the mth site equals 2, the mth fermionic site and two surrounding spin sites admit the following configuration \(\downarrow \bullet \uparrow\); similarly, when the eigenvalue of \(\hat{G}_m\) at the mth site equals \(-1\), we obtain a configuration \(\uparrow \circ \downarrow\).Footnote 1 These two configurations are dynamically frozen. Other fermions cannot jump into this site and change the spin configuration on the nearby links. Therefore, the gauge configurations with the eigenvalue of bulk \(\hat{G}_m\) being 2 or \(-1\) further subdivide the lattice into spatially disconnected parts.

Fig. 3
figure 3

Different blocking configurations in the minimal quantum breakdown model in Eq. (26)


Furthermore, these two blocking configurations can generate more complicated blocking structures. As shown in Fig. 3c and d, the blocking configurations \(\downarrow \bullet \uparrow\) and \(\uparrow \circ \downarrow\) can be further extended by adding more parallel spins on the two sides. These extended blocking configurations form a domain-wall structure of the link spins, while the fermionic sites within the blocking region are either fully occupied or empty.

The existence of many blocking configurations further reduces the dimension of connected Hilbert subspaces. In the following, we only need to focus on the situations where the eigenvalues of all the bulk \(\hat{G}_m\) are 0 or 1. In other words, we are free to remove the blocking region because of their frozen dynamics. Then, the dynamically connected region consists of the following configurationsFootnote 2:

\(\begin{aligned} & \langle {\hat{G}_m}\rangle =0: \quad \downarrow \circ \downarrow ,\ \uparrow \circ \uparrow ,\ \uparrow \bullet \downarrow ;\nonumber \\ & \langle {\hat{G}_m}\rangle =1: \quad \downarrow \bullet \downarrow ,\ \uparrow \bullet \uparrow ,\ \downarrow \circ \uparrow . \end{aligned}\)
(31)

Here, \(\langle {\hat{G}_m}\rangle\) represents the expectation value of \(\hat{G}_m\) over a product state, which is also an eigenstate due to the diagonal structure of \(\hat{G}_m\). Moreover, the combined conserved quantity \(\hat{G}=\sum \nolimits _m \hat{G}_m\) corresponds to the total number of fermions in the quantum link model. This quantity equals to the number of fermions \(\hat{\mathcal {N}}_1\) on the first orbitals in the quantum breakdown model. Consequently, we can investigate the dynamically connected Krylov subspaces labeled with different fermion numbers, or fermion filling factors, on the first orbitals. In the following discussion about the subspace dynamics, we restrict ourselves into the Krylov subspaces below the half-filling, since the Krylov subspaces above half-filling can be readily generated via a particle-hole transformation.

4.2 Free-fermion sector

The simplest situation is the case with only one fermion, namely, the eigenvalue of \(\hat{G}=\sum \nolimits _m \hat{G}_m\) is equal to 1. Under this circumstance, the connected Krylov subspace is generated by a reference state like Fig. 4a. If we label the state by the position of the fermion as \(|{m}\rangle\), such a state in the original quantum breakdown model can be expressed as

\(\begin{aligned} |{m}\rangle =\hat{f}_{m,1}^\dagger \prod \limits _{j=2}^{m}\left( \prod \limits _{i=2}^{2q+1}\hat{f}_{j,i}^\dagger \right) |{0}\rangle =\hat{f}_{m,1}^\dagger \prod \limits _{j=2}^{m}\hat{F}_j^\dagger |{0}\rangle . \end{aligned}\)
(32)

Fig. 4
figure 4

Typical configurations in a Krylov subspace with one fermion. The integers above show the eigenvalues of local conserved \(\hat{G}_m\)


The reduced Hamiltonian in this Krylov subspace admits the following form:

\(\begin{aligned} H_{\text {free}}=\sum \limits _{m=1}^{M-1}(J_m|{m+1}\rangle \langle {m}|+\text {h.c.})+\sum \limits _{m=1}^MV_m|{m}\rangle \langle {m}|, \end{aligned}\)
(33)

where the onsite potential is given by \(V_m=\mu _m+\sum \nolimits _{i=2}^{m}2q\mu _i\). This is thus effectively a single-particle tight-binding model [38].

Notably, by tuning the parameters \(J_m\) and \(\mu _m\), this model can exhibit different behaviors. For example, if we take a uniform \(J_m=J\) and set \(\mu _m=0\), this subspace describes a tight-binding model whose eigenstates are Bloch waves:

\(\begin{aligned} |{\psi _k}\rangle =\sqrt{\frac{2}{M+1}}\sum \limits _{m=1}^{M}\sin \left( \frac{km\pi }{M+1}\right) |{m}\rangle ,\quad k=1,\dots ,M. \end{aligned}\)
(34)

Consequently, the eigenvalues are given by \(E_k=2J\cos (\frac{k\pi }{M+1})\).

Another interesting setup is to take a constant on-site potential \(\mu _m=\mu\). In this case, \(V_m=(2q(m-1)+1)\mu\) becomes a linear potential, resulting in the Wannier-Stark localization of the free fermion. Furthermore, if \(J_m\) and \(V_m\) are taken from some random distributions, the system will display the Anderson localization. The localization phenomenon in this subspace indicates that the original quantum breakdown model displays many-body localization in certain charge sectors.

4.3 Boundary interaction

We now analyze Krylov subspaces that involve more than one fermion. In the subspaces with two fermions, two sites are specified with their eigenvalue of \(\hat{G}_m\) being 1. A representative configuration is shown in Fig. 5. In this scenario, both particles move freely within their individual dynamical zones, with the exception of a contact interaction occurring near the boundary that separates them.

Fig. 5
figure 5

Typical configurations in a Krylov subspace with two fermions. The integers above show the eigenvalues of local conserved \(\hat{G}_m\)


For example, the dynamical region of the first fermion in Fig. 5a contains the leftmost six sites, while that of the second fermion comprises the rightmost six sites. Starting from the reference state in Fig. 5a, if the second particle moves to the right, the first particle can migrate to the initial site of the second particle (Fig. 5c) . However, due to the constraint imposed by the spin configurations, the first particle cannot move further in the right direction. On the other hand, if the second particle remains at its initial location, the first particle cannot occupy the same site due to the Pauli principle. Consequently, these two particles exhibit an effective interaction near the boundary between their dynamical regions [38]. This boundary interaction results in challenges for analytical solutions in this Krylov subspace.

The connectivity graph of this Krylov subspace is shown in Fig. 6. This graph takes the form of a square lattice, with the absence of a corner site. In the case of \(J_m=J\) and \(\mu _m=0\), if the distance between two fermions in the reference state (Fig. 5a) is sufficiently large, we can ignore the boundary defect and treat the dynamics of this Krylov subspace as a free particle on a 2D square lattice. This is reasonable because the probability of finding a particle at the boundary decreases as the length of its dynamical zone increases.

Fig. 6
figure 6

The connectivity of the Krylov subspace shown in Fig. 5. The nodes are product states in the Krylov subspace and the edges are the nontrivial interactions induced by the Hamiltonian. The locations of configurations in Fig. 5 are labeled by the corresponding letters


The picture based on the subspace connectivity can be generalized to situations involving generically n fermions. A representative connectivity graph for the subspace with \(n=3\) fermions is shown in Fig. 7. When the fermion density in the system is sufficiently small, the average distance between two adjacent particles in the reference state like Fig. 5a becomes significantly large. As a result, the particles have a neglected probability of simultaneously appearing at the boundary between their dynamical regions. In this context, such a boundary interaction may be considered weak and, therefore, can be neglected. Consequently, the effective dynamics can be viewed as a free particle moving on a n-dimensional hypercubic lattice.

Fig. 7
figure 7

The connectivity of the Krylov subspace with three fermions. The reference state \(|{\bullet \downarrow \circ \downarrow \circ \downarrow \bullet \downarrow \circ \downarrow \circ \downarrow \bullet \downarrow \circ \downarrow \circ \downarrow \circ }\rangle\) is shown below the graph


4.4 The intermediate cases

Consider increasing the fermion density in the reference state. While ensuring that the density remains below half-filling, an increase in fermion density leads to a decrease in the average distance between two adjacent particles. Roughly speaking, the decrease in particle separation leads to strong contact interactions between neighboring particles. Consequently, the Krylov subspace becomes strongly interacting, and the dynamics dramatically deviate from those of a nearly free particle moving on a hypercubic lattice.

It is very difficult to develop an analytical description for Krylov subspaces with high densities. Therefore, we perform numerical investigations on these Krylov subspaces here. For simplicity, we focus on the situations where \(J_m=J=1\) and \(\mu _m=0\). Figure 8 shows the subsystem entanglement entropy of the eigenstates. With an eigenstate \(|{\psi }\rangle\), we define \(\rho _A = \text {Tr}_{\bar{A}}[|{\psi }\rangle \langle {\psi }|]\) as the reduced density matrix of the subsystem A, where \(\bar{A}\) is the complementary set of A. Then, the subsystem entanglement entropy is obtained by \(S_A = -\text {Tr}[\rho _A\ln \rho _A]\). The numerical findings indicate that Krylov subspaces with many fermions exhibit a nonchaotic feature, as evidenced by many low-entangled eigenstates located in the central region of the spectrum. This property is linked to the boundary interactions between adjacent fermions in this kinetically constrained model.

Fig. 8
figure 8

The subsystem entanglement entropy of eigenstates with many fermions. With the Hamiltonian in Eq.(26), we set \(J_m=J=1\) and \(\mu _m=0\). In each plot, \(\nu =N_1/M\) denotes the filling factor of fermions, where \(N_1\) is the fermion number and M is the number of fermion sites. a \(N_1=5\) and \(M=25\). The gauge configuration for \(\hat{G}_m\) in Eq. (27) is set to 1000010000100001000010000. The entanglement entropy is evaluated for subsystem A which includes the leftmost 13 fermion sites and 12 link spins. b \(N_1=6\) and \(M=24\). The gauge configuration for \(\hat{G}_m\) is set to 100010001000100010001000. The entanglement entropy is evaluated for subsystem A containing the leftmost 12 fermion sites and 12 link spins. c \(N_1=8\) and \(M=24\). The gauge configuration for \(\hat{G}_m\) is set to 100100100100100100100100. The entanglement entropy is evaluated for subsystem A which includes the leftmost 12 fermion sites and 12 link spins. d \(N_1=11\) and \(M=22\). The gauge configuration for \(\hat{G}_m\) is set to 1010101010101010101010, similar to Fig. 9. The entanglement entropy is evaluated for subsystem A that includes the leftmost 11 fermion sites and 11 link spins


4.5 PXP model and quantum many-body scar

An interesting Krylov subspace with a high density of fermions is illustrated in Fig. 9. This corresponds to the half-filling case, where the reference state has fermions at odd sites. Within the Krylov subspace generated by this reference state, we can introduce a new gauge generator:

\(\begin{aligned} \widetilde{G}_m=\hat{G}_m-\frac{1-(-1)^m}{2}. \end{aligned}\)
(35)

Fig. 9
figure 9

Typical configurations of the minimal quantum breakdown model in Eq. (26) in the half-filling Krylov subspace given by Eq. (36), which maps to the PXP model. The integers above show the eigenvalues of local conserved \(\hat{G}_m\)


Then, all the states in this subspace satisfy

\(\begin{aligned} \widetilde{G}_m|{\psi }\rangle =0. \end{aligned}\)
(36)

This precisely corresponds to the physical gauge sector discussed in lattice gauge theory for staggered fermions [45,46,47].

Remarkably, despite the strong interactions within this subspace, the dynamics starting from the reference state reveal persistent revivals, a distinctive feature of quantum many-body scars [34,35,36,37, 56]. Notably, the Hamiltonian in this subspace of the quantum link model can be precisely mapped to the PXP model [56, 67], a system simulated with Rydberg atoms [12,13,14]. As shown in Fig. 9, if we map our spins on the links in Eq. (26) into Pauli matrices \(\hat{\tau }^{x,y,z}_m\) through [56]

\(\begin{aligned} & \hat{\tau }^z_m\leftrightarrow (-1)^{m}\hat{\sigma }^z_{m,m+1},\nonumber \\ & \hat{\tau }^x_m\leftrightarrow \hat{f}_{m+1,1}^\dagger \hat{\sigma }^+_{m,m+1}\hat{f}_{m,1}+\hat{f}_{m,1}^\dagger \hat{\sigma }^-_{m,m+1}\hat{f}_{m+1,1},\nonumber \\ & \hat{\tau }^y_m\leftrightarrow -i(-1)^{m}(\hat{f}_{m+1,1}^\dagger \hat{\sigma }^+_{m,m+1}\hat{f}_{m,1}-\hat{f}_{m,1}^\dagger \hat{\sigma }^-_{m,m+1}\hat{f}_{m+1,1}). \end{aligned}\)
(37)

With \(J_m=J\) and \(\mu _m=0\), the model Eq. (26) in the gauge sector of Eq. (36) maps to the PXP Hamiltonian for Rydberg atoms:

\(\begin{aligned}\hat{ H}_{\text{PXP}}=\frac{J}{4}\sum \limits _{m}(1-\hat{\tau }^z_{m-1})\hat{\tau }^x_{m}(1-\hat{\tau }^z_{m+1})\ . \end{aligned}\)
(38)

This mapping can be understood intuitively. With the identification \(\hat{\sigma }_{m,m+1}^z\leftrightarrow (-1)^m\hat{\tau }_m^z\), the bulk gauge symmetry generators Eq. (35) become \(\widetilde{G}_m=\hat{f}_{m,1}^\dagger \hat{f}_{m,1}+\frac{1}{2}[(-1)^m(\hat{\tau }_m^z+\hat{\tau }_{m-1}^z+1)-1]\), where \(1. With \(|{\psi }\rangle\) in the constrained Hilbert space being a basis vector, the gauge constraint Eq. (36) leads to \(\langle {\psi |\hat{\tau }_m^z+\hat{\tau }_{m-1}^z|\psi }\rangle =(-1)^m[1-2\langle {\psi |\hat{f}_{m,1}^\dagger \hat{f}_{m,1}|\psi }\rangle ]-1<2\), which corresponds to the Rydberg blockade that two neighboring atoms cannot be excited simultaneously to Rydberg states. Consequently, we anticipate observing many-body scar dynamics when evolving a \(Z_2\) configuration in Fig. 9a.

In Fig. 10a, we show the overlap between the \({Z}_2\) configuration and the energy eigenstates in the subspace under the OBC. Our numerical findings distinctly reveal the presence of scar tower structures. Additionally, the fidelity \(F(t) = |\langle {\psi _{Z_2}}|e^{-iHt}|{\psi _{Z_2}}\rangle |^2\) during the time evolution exhibits clear and notable revivals, as shown in Fig. 10b. These outcomes collectively support the existence of quantum many-body scar states within the quantum link model, i.e., the minimal quantum breakdown model.

Fig. 10
figure 10

Quantum many-body scars in the half-filling Krylov subspace. a Overlap of the \(Z_2\) configuration \(|{\psi _{Z_2}}\rangle\) like Fig. 9a with the energy eigenstates in this subspace. b The fidelity dynamics starting from \(|{\psi _{Z_2}}\rangle\). The parameters are \(J_m=J=1,\mu _m=0\), and the number of fermionic sites is \(L=22\). We take the open boundary condition in the numerical simulation


5 Discussions

In this paper, we study the minimal quantum breakdown model with \(N=2q+1\) and investigate its understanding from the perspective of lattice gauge theory. An extensive number (proportional to the system size) of locally conserved quantities in the model leads to the emergence of Hilbert space fragmentation, which is closely tied to numerous dynamically blocking sites. The mapping between the minimal quantum breakdown model and the U(1) lattice gauge theory offers a powerful framework for employing gauge configurations to describe the various dynamics within the connected Krylov subspaces of the original quantum breakdown model. These intrinsic connections among the minimal quantum breakdown model, quantum link model, and lattice gauge theory may motivate proposals for experimentally implementing the generic quantum breakdown model and more quantum models with generic modulated symmetries.