Next Article in Journal
Analytic Formulae for T Violation in Neutrino Oscillations
Previous Article in Journal
View-Driven Multi-View Clustering via Contrastive Double-Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hilbert Space Delocalization under Random Unitary Circuits

1
Institut für Theoretische Physik, Universität zu Köln, Zülpicher Strasse 77a, 50937 Cologne, Germany
2
ICFO—Institut de Ciències Fotòniques, The Barcelona Institute of Science and Technology, Av. Carl Friedrich Gauss 3, 08860 Castelldefels, Barcelona, Spain
*
Author to whom correspondence should be addressed.
Entropy 2024, 26(6), 471; https://doi.org/10.3390/e26060471
Submission received: 22 April 2024 / Revised: 21 May 2024 / Accepted: 24 May 2024 / Published: 29 May 2024
(This article belongs to the Section Quantum Information)

Abstract

:
The unitary dynamics of a quantum system initialized in a selected basis state yield, generically, a state that is a superposition of all the basis states. This process, associated with the quantum information scrambling and intimately tied to the resource theory of coherence, may be viewed as a gradual delocalization of the system’s state in the Hilbert space. This work analyzes the Hilbert space delocalization under the dynamics of random quantum circuits, which serve as a minimal model of the chaotic dynamics of quantum many-body systems. We employ analytical methods based on the replica trick and Weingarten calculus to investigate the time evolution of the participation entropies which quantify the Hilbert space delocalization. We demonstrate that the participation entropies approach, up to a fixed accuracy, their long-time saturation value in times that scale logarithmically with the system size. Exact numerical simulations and tensor network techniques corroborate our findings.

1. Introduction

Starting from a vector of a generically chosen basis of the Hilbert space, unitary quantum dynamics generate a superposition spanning the entire basis of the Hilbert space. This phenomenon, referred to as Hilbert space delocalization, can be viewed as the spreading of the many-body wave function over the Hilbert space under quantum dynamics. Hilbert space delocalization is tied to non-equilibrium processes in quantum mechanics and, hence, pivotal in our understanding of quantum foundations [1], quantum technologies [2,3], and condensed matter theory [4,5,6]. Being intimately connected to the resource theories of quantum coherence and entanglement [7,8], Hilbert space delocalization provides a valuable tool for the characterization of quantum phases of matter [9,10,11,12,13,14], quantum chaos and thermalization [15,16,17,18,19,20,21,22], and their violations in non-ergodic systems—including monitored systems [23,24,25,26,27,28] and disorder-induced localization [29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46].
Natural quantifiers for the localization and delocalization properties of a state ρ = | Ψ Ψ | over the many-body Hilbert space are the inverse participation ratios (IPRs) and the participation entropies, which measure the spreading of the state distribution p n n | ρ | n over the basis B = { | n } . The IPR and participation entropy are given, respectively, by
I q n B n | ρ | n q = n B p n q , S q = 1 1 q ln I q .
(We note the connection with the relative entropies of coherence [47] C q ( ρ , B ) S q ( 1 q ) 1 ln tr ( ρ q ) , a scalable quantifier in quantum coherence resource theory. In particular, for pure states, these two quantities coincide, C q ( ρ , B ) = S q ).
Despite the explicit basis dependence, extensive studies have demonstrated that the IPR captures the structural properties of quantum matter, including universal behavior at phase transitions [26,48]. However, the investigations of the IPR and participation entropies have focused so far on the equilibrium and stationary state features, leaving the question of the time evolution of quantum coherence, and, hence, of Hilbert space delocalization, unresolved.
In this work, we investigate Hilbert space delocalization under the dynamics of ( 1 + 1 ) D quantum circuits comprising 2-qudit Haar random gates arranged in a brick wall pattern [23]. The locality and unitarity of this setup constitute the minimal requirements for chaotic evolution in many-body systems [49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70], hence allowing us to gain a phenomenological understanding of the Hilbert space delocalization under generic quantum dynamics. Our analysis combines rigorous analytical methods, based on the mapping between the average of the IPR over the random circuits and a statistical mechanics model, with exact numerical simulations, including tensor network techniques [71,72,73,74] (here, our implementation is based on the open-source library ITensor [75,76] and available at [77]). We find that the saturation of participation entropy S q to the long-time stationary value, equal to the random Haar state average S ˜ q H , occurs exponentially quickly in time (i.e., in circuit depth):
S q = S ˜ q H α q N ( 2 K d ) t 1 ,
where K d > 0 is a constant characterizing the properties of 2-qudit gates, and α q > 0 is a constant. We find an exact expression for the constant K d for Haar random 2-qudit gates and provide analytic arguments for the value of α q in the limit of large on-site Hilbert space dimension. In other words, Equation (2) implies that the participation entropy S q approaches its stationary value, up to a fixed accuracy, at a time scaling logarithmically with the system size, τ HSD ln ( N ) .
The paper is structured as follows. In Section 2, we outline the toolbox employed in the rest of the paper, and, in Section 3, we discuss the stationary value of participation entropy that the deep quantum circuit reaches. Section 4 presents the core of our analytical approach based on the mapping of the IPR calculation to a statistical mechanics problem. We resolve the Rényi-2 participation entropy evolution in Section 5, where we present the analytic prediction. These findings are further corroborated by the numerical analysis in Section 6 for generic Rényi index q. We present the outlook of our work and a discussion of its further implications in Section 7. The more technical parts of the paper are detailed in the appendix. In Appendix A, we present a self-contained discussion about Haar averages. In Appendix B, we review the interface problem related to the entanglement propagation and, in Appendix C, briefly discuss the numerical implementations.

2. Methods

We begin by presenting the analytical toolbox we will employ in the rest of the paper. The latter is based on the replica trick, the superoperator formalism [78], and Weingarten calculus [79,80,81,82] (see also Refs. [54,61]). Since n | | m = tr ( | m n | ( ) ) and tr ( A B ) q = tr ( A q B q ) , the IPR definition (1) implies that I q = tr ( Λ B ρ q ) , where Λ B = n B | n n | q is a replica operator acting on H N q . Throughout this manuscript, we will study qudits with local Hilbert space dimension d; hence, H N C d N . For concreteness, we will fix B to be the computational basis, where | n = | b 1 , , b N with b j = 0 , , d 1 for each j. Nevertheless, due to the invariance of the considered circuits under local basis rotations, our results remain valid for any basis of H N obtained by a unitary transformation U = U 1 U N of the computational basis, where U k belongs to the unitary group U ( d ) for a single qudit. For the computational basis, Λ B = j = 1 N Λ q ( j ) , which defines the “book” boundary condition for each qudit [83,84]:
Λ q ( j ) b j = 0 , , d 1 ( | b j b j | ) q ,
where the superscript in Λ q ( j ) bookkeeps the number of the affected qudit. For later convenience, we recast the problem in the superoperator formalism. For any operator U and A, we have A = n , m = 0 d N 1 A n , m | n m | | A = n , m = 0 d N 1 A n , m | n , m , and U A U ( U U * ) | A (we denote with ( ) * and ( ) , respectively, the complex and Hermitian conjugation of •). In this representation, B | A = tr ( B A ) , so the IPR can be written as
I q = Λ q ( 1 ) , Λ q ( 2 ) , , Λ q ( N ) | ρ q .
For convenience, we introduce a graphical notation. We implicitly define it for the replica boundary Λ q and the replica state ρ q via the IPR as follows:
Entropy 26 00471 i001
Lastly, since we will consider the unitary evolution | ρ = ( U U * ) | ρ 0 for some | ρ 0 values selected as the initial state and unitary operation U acting on the Hilbert space H N , we introduce the following graphical notation for U and U * , respectively:
Entropy 26 00471 i002
We will omit the legs indices, which can be inferred from the context. Moreover, we will typically reserve thick (thin) lines for multi-replica (single-replica) objects.

3. Delocalization Properties of Random Haar States

We begin our discussion by reviewing the Hilbert space delocalization of random states [21,83] which correspond to the stationary ensemble of states obtained under the action of sufficiently deep random circuits. Uniformly distributed random states in the Hilbert space | ρ Haar = ( U U * ) | ρ 0 are obtained from a reference state | ρ 0 = | Ψ 0 , Ψ 0 via a unitary operation U acting globally on the system of N qudits, drawn with the Haar measure from the N qudit unitary group U ( d N ) , where d is the qudit local Hilbert space dimension. Let us compute the average IPR over the Haar ensemble, which reads, using the graphical notation introduced in Section 2, as
Entropy 26 00471 i003
where E Haar denotes the expectation value over the unitary group U ( d N ) taken with the Haar measure. Here, we have defined, via the linearity of the average and of the expectation value, the q-replica transfer matrix on N qudits:
Entropy 26 00471 i004
The explicit formula for W N ( q ) can be obtained by a direct evaluation of the Haar average, as detailed in Appendix A. The final result is expressed in terms of the permutation operators | σ j acting on each qudit j, with the matrix elements in the superoperator computational basis given as
b 1 j , b ¯ 1 j , b 2 j , b ¯ 2 j , , b q j , b ¯ q j | σ j = k = 1 q δ b k j , b ¯ k σ ( j ) ,
fixed by the permutation σ S q . We employ the Weingarten function Wg ( D ; σ ) allows us to define the dual states | σ ˜ 1 , 2 , , N = τ S q Wg ( d N ; σ τ 1 ) | τ 1 | τ N acting on the whole system of N qudits. With these states, the transfer matrix reads
Entropy 26 00471 i005
We are now in position to compute the IPR. Since the initial state ρ 0 is pure, σ , , σ | ρ 0 q = 1 . Furthermore, using σ Wg ( d N ; σ τ 1 ) = ( d N 1 ) ! / ( d N + q 1 ) [81], we have
Entropy 26 00471 i006
where we use | S q | = q ! and Λ q | τ = d for all τ S q . From this calculation, it follows that the annealed averaged Rényi participation entropy for the Haar states is given by
S ˜ q H 1 1 q ln I ¯ q H = ln [ ( d N + 1 ) ( d N + q 1 ) ] ln [ q ! ] q 1 N N ln [ d ] 1 q 1 ln [ q ! ] .
As expected and already discussed in [21], the Haar random states are (almost) fully delocalized over the many-body basis B . Indeed, (12) shows that S ˜ q H differs only by a sub-leading constant term ( q 1 ) 1 ln [ q ! ] from the maximal value of participation entropy N ln ( d ) achieved for a uniformly distributed state with n | ρ | n = d N . Since the logarithm is non-linear, in principle, one should expect corrections when considering the quenched average of the participation entropy S ¯ q H E Haar [ S q ] . We briefly estimate these fluctuations via the variance (for the explicit expression of I q r ¯ , see Ref. [83]).
std ( I q ) I q 2 ¯ I q ¯ 2 .
For this expression, we need to compute I q 2 ¯ , which requires the use of 2 q replicas of H N and the boundary conditions Λ q Λ q | acting on each site. Algebraic manipulations analogous to the ones described above give
Entropy 26 00471 i007
It follows that std ( I q ) O d ( 1 2 q ) N / 2 , and, in the thermodynamic limit, these fluctuations are irrelevant. In particular, in the scaling limit N , the annealed and the quenched averages coincide, S ˜ q H S ¯ q H N 0 , and, with probability equal to unity, approximate the value of the participation entropy calculated for a single Haar random state with an arbitrarily small fixed accuracy.

4. Hilbert Space Delocalization in Brick Wall Quantum Circuits

After this preliminary discussion, relevant to the deep circuit limit, we now discuss the Hilbert space delocalization under random quantum circuits. As the initial state, we fix | Ψ 0 = | 0 N and study the approach of the annealed and quenched average participation entropy to the asymptotic value S ˜ q H = ( 1 q ) 1 ln [ I q H ] with I q H q ! ( d N ) ! / ( d N + q 1 ) ! , cf. Equations (11) and (12). (For the circuit case, the quenched and annealed averages are given, respectively, by
S ¯ q = E Haar [ ( 1 q ) 1 ln I q ( t ) ] , S ˜ q = ( 1 q ) 1 ln E Haar [ I q ( t ) ] ,
with I q ( t ) = b | b | U t | Ψ 0 | 2 q ). The evolution operator corresponding to the brick wall quantum circuit of depth t (referred to also as time) is given by
U t = τ = 0 t 1 U ( τ ) , U ( τ ) = i = 0 N / 2 ( τ % 2 ) U 2 i + ( τ % 2 ) , 2 i + 1 + ( τ % 2 ) ,
where m % k denotes the integer m modulo and the integer k. Here, each U m , n is an independent identically distributed (i.i.d.) Haar random unitary from U ( d 2 ) acting on the m-th and n-th qudit. Our goal is to compute the
I ¯ q E Haar Λ q , , Λ q | ( U t U t * ) q | ψ 0 , ψ 0 , , ψ 0 ,
where the N copies of | ψ 0 = | 0 2 q implement the initial condition | Ψ 0 on the replica space. We note that the E Haar denotes here the average over the realizations of the two-body gates U m , n which comprise the considered brick wall circuit. We adapt the graphical notation of Section 2 for a generic two-body gate U with indices i , j , k , l { 0 , d 1 } numbering the states of the qudits:
Entropy 26 00471 i008
Recalling that the two-body gates U m , n are i.i.d. variables, we observe that their averages factorize, implying that the transfer matrix corresponding to the circuit (16) of depth t reads
Entropy 26 00471 i009
Contracting with the replica boundaries, and recalling the decomposition for W 2 ( q ) in Equation (18), we obtain the final expression
Entropy 26 00471 i010
where, denoting | σ ˜ i , i + 1 = τ S q Wg ( d 2 ; σ τ 1 ) | τ i | τ i + 1 , we use σ | ψ 0 = 1 for each site, ( Λ q | i Λ q | i + 1 ) | σ ˜ i , i + 1 = d 2 , and we define | + i , i + 1 = σ S q | σ i | σ i + 1 and | + ˜ i , i + 1 = σ S q | σ ˜ i , i + 1 . (We note that the contraction between the W 2 ( q ) requires τ | i · | σ i = d # ( σ τ 1 ) , with # ( σ ) the number of cycles of the permutation σ . These weights can be reabsorbed in W 2 ( q ) via a local basis change, see Appendix A).
Equation (20) is one of the main results of this work as it already allows for efficient tensor network implementations [75]. At the same time, Equation (20) can be viewed as the partition function of a statistical mechanics model once the W 2 ( q ) values are replaced with decomposition Equation (18). To highlight this interpretation, we can unravel the upper boundary condition as
| Q in = { σ j S q } | Q in ( { σ j } ) | Q in ( { σ j } ) i = 1 N / 2 | σ j i | σ j i + 1 .
It follows that W t is a sum of backward paths from each initial state | Q in ( { σ j } ) to each final | Q fin ( { σ j } ) i = 1 N / 2 | σ ˜ j i | σ j i + 1 .
We conclude this section by highlighting the similarities and the differences between the IRP calculation and the computations for propagation of entanglement entropy [54,55], defined as S q ent = ln [ tr ( ρ A q ) ] , where ρ A = tr A c ( ρ ) is the reduced density matrix of the pure state ρ = | Ψ Ψ | and A A c a bipartition of the system. As for the participation entropy, we focus on the annealed averages, determined by the purities, i.e., the reduced density matrix moments P q E Haar [ tr ( ρ A q ) ] . Denoting the identity permutation as I S q ( ) and the cyclic permutation as S S q ( 12 q ) , in the superoperator formalism, we have
P q ( t ; A ) = Q in ( { S j : j A , I l : l A c } ) | W t ( q ) | ψ 0 N ,
where | Q in ( { S j : j A , I l : l A c } ) is a product of identity (cyclic) operators on the qudit j, for j A c ( j A ). The statistical mechanics model corresponding to the calculation of P q ( t ; A ) is exactly the same as Equation (20) but with different boundary conditions; for the IPR, we have free boundaries, cf. Equation (21), while, for the P q ( t ; A ) , the boundary conditions are of the domain wall type. This fact leads to substantial differences between the Hilbert space delocalization and the entanglement propagation in the considered random circuits, which we will discuss in the next section focusing on the two-replica limit. We will further corroborate these findings with a numerical analysis for various replica numbers (Rényi index) q.

5. Two-Replica Computations

In the two-replica limit, the parallelism between the inverse participation ratios and the reduced density matrix moments is more direct. Since S 2 = { I , S } , Equation (21) translates to the sum over all domain wall boundary conditions (22). Therefore, the IPR, up to the overall constant, is given by the sum of P 2 over all possible choices of subsets A of contiguous pairs of sites Q { { 1 , 2 } , { 3 , 4 } , , { N 1 , N } } , cf. Equations (20) and (21), resulting in
I ¯ 2 ( t + 1 ) = 1 ( d 2 + 1 ) N / 2 A Q P 2 ( t ; A ) ,
where the shift to time t + 1 comes from the simplification of one layer in Equation (20). The time evolution of the purities [54,55] can be split into two terms:
P 2 ( t ; A ) = d 2 N A + d 2 ( N / 2 N A ) d N + 1 + f A ( t ) ,
where the first term describes the long-time saturation value of the purity, the second term parametrizes the approach to the saturation value, and the factor 2 in the exponents arises because the set Q contains N / 2 pairs of neighboring lattice sites. Combining Equations (23) and (24) leads to
I ¯ 2 ( t + 1 ) = 1 ( d 2 + 1 ) N / 2 A Q d 2 N A + d 2 ( N / 2 N A ) d N + 1 + f A ( t ) = I 2 H + 1 ( d 2 + 1 ) N / 2 A Q f A ( t ) ,
where the first term in Equation (24) is independent of the partition. This fact implies that
A Q d 2 N A + d 2 ( N / 2 N A ) = 2 N A = 0 N / 2 N / 2 N A d 2 N A = 2 ( d 2 + 1 ) N / 2 ,
and, as a result, we have
S ˜ 2 H S ˜ 2 ( t ) = ln 1 + d N + 1 2 ( d 2 + 1 ) N / 2 A Q f A ( t 1 ) .
We are interested in identifying the leading terms governing the decay of S ˜ 2 H S ˜ 2 ( t ) at long times. Fixing the system size N and considering the limit of d N , the factor ( d N + 1 ) / [ 2 ( d 2 + 1 ) N / 2 ] in Equation (26) simplifies to unity (we note that the d limit has been extensively considered in the literature, see Ref. [49]). We focus first on the single domain wall configurations, corresponding to bipartitions of the system into A = { 1 , 2 , , 2 N A } and A c = { 2 N A + 1 , , N } . The exact solution for the single domain wall case, detailed in Appendix B, shows that f A ( t ) ( 2 K d ) t for domains localized in the bulk of the system, where
K d = d d 2 + 1 .
The time evolution of the purities P 2 ( t ; A ) , where A corresponds to a configuration with n d domains, is suppressed by the factor ( 2 K d ) n d t , which vanishes at any t > 0 in the large d limit. Therefore, the leading contributions to the S ˜ 2 H S ˜ 2 ( t ) arise due to the single domain configurations. The number of the relevant single domain wall configurations scales proportionally to the system size N, translating Equation (26) to
S ˜ 2 H S ˜ 2 ( t ) d 1 α d N ( 2 K d ) t 1
where the coefficient α d d 1 1 / 2 at sufficiently large system size N. While we have used the large on-site Hilbert space dimension limit to derive Equation (28), we anticipate that this equation applies to any d 2 with a properly chosen factor α d .
We are unable to analytically demonstrate the validity of Equation (28) at a fixed finite value of d and for N 1 . Nevertheless, we outline the underlying calculations, which allow us to pin-point the leading factors in the time dependence of S ˜ 2 H S ˜ 2 ( t ) . First, we observe that, at fixed d and for N 1 , the contributions from the exponentially many configurations A Q are exponentially in N suppressed by the term ( d N + 1 ) / [ 2 ( d 2 + 1 ) N / 2 ] . Performing tensor network contractions, we find for a generic initial configuration that f A ( t ) = a A ( 2 K d ) t in the long-time limit. Interpreting the tensor network contraction as the wandering domain walls problem (cf. Appendix B), the factor ( 2 K d ) t corresponds to the single domain wall configuration which dominates f A ( t ) at long times, while the factor a A depends on the processes required to reach the single domain wall configuration from the initial condition. The summation over all initial conditions yields an overall factor ( 2 K d ) t which describes the behavior of S ˜ 2 H S ˜ 2 ( t ) at any d 2 . We conjecture that the competition of the renormalizing factors a A with the exponential terms appearing in the summation over A Q yields a factor proportional to the system size N at any d 2 . This conjecture is equivalent to the validity of Equation (28) at any d 2 .
To corroborate our analytical considerations, we compare our predictions with the numerical simulations of the brick wall random circuits varying the on-site Hilbert space dimension d. A brief summary of our numerical implementations is detailed in Appendix C. First, we numerically calculate the exact time evolution of the system’s state | Ψ up to the system size N = 24 . Our findings are reported in Figure 1 (left). We observe that, already at the time scale O ( 1 ) , the difference between the annealed average S ¯ 2 and the quenched average S ˜ 2 of the participation entropy is negligible. Thus, the self-averaging properties of the IPR and participation entropies derived in the long-time limit in Section 3 apply also for the relatively shallow circuits, allowing us to use the S ˜ 2 ( t ) = ln [ I ¯ 2 ( t ) ] as an accurate proxy for the circuit-averaged participation entropy S ¯ 2 ( t ) .
To test our analytical prediction (28), we focus on the difference between the stationary Haar value S ˜ 2 H and the annealed average S ˜ 2 ( t ) . We implement Equation (20) as a tensor network contraction, which allows us to reach system sizes N 1024 for any on-site Hilbert space dimension d. In Figure 1 (right) we can see that the difference S ˜ 2 H S ˜ 2 ( t ) decreases exponentially in time (circuit depth) t proportionally to ( 2 K d ) t 1 or at a slightly quicker pace at the smallest considered system sizes N = 8 , 16 . The prefactor of the exponential decay increases monotonically with system size N. The growth of the prefactor becomes clearly linear beyond N = 64 . At N > 64 , the behavior of S ˜ 2 H S ˜ 2 ( t ) at longer times is captured with good accuracy by Equation (28) with α 2 = 0.291 ( 5 ) .
Performing tensor network calculations for the on-site Hilbert space dimension d = 3 , , 13 and for system sizes up to N = 256 , we verify that S ˜ 2 H S ˜ 2 ( t ) decreases exponentially in time, matching the prediction in Equation (28) and reproducing the analytically found value of the constant K d . We also find that the prefactor α d increases with d, as shown in the inset in Figure 1 (right). In particular, the numerical results indicate that α d d 1 / 2 , in accordance with the analytical prediction for the large on-site Hilbert space dimension limit.
The results discussed in this section demonstrate that the Hilbert space delocalization under the brick wall random quantum circuits is an abrupt process. Indeed, for a given tolerance error ε 1 , the participation entropy reaches S ˜ 2 ( t ) = S ˜ 2 H ε at time t HSD ln ( N ) , which scales logarithmically with the system size N.

6. Numerical Results for Any Replica

For generic values of the Rényi index q, the analysis presented above cannot be trivially extended. First, the calculation of I ¯ q requires all domain wall configurations | Q 0 ( { σ j } ) with σ j S q . Thus, the calculation of the IPR involves a sum over many more initial configurations than just the configurations appearing in the computation of the purities P q , which includes only the identity permutation I and the cyclic permutation S . Moreover, W 2 ( q ) involves contributions with negative signs and additional weights, which affect the value of K d , resulting in more complex analysis.
Nonetheless, the main features of the participation entropy growth under brick wall quantum circuits discussed in Section 5 may be expected to hold for a generic replica number q. This leads us to the following conjecture:
S ˜ q H S ˜ q ( t ) α d , q N β d , q t S ˜ q H S ˜ q ( t ) = ε at t HSD ln ( N ) ,
where α d , q is a coefficient such that α d , 2 = α d in Equation (28), β d , q < 1 is a constant with β d , 2 = 2 K d , and ε 1 is a fixed tolerance. The heuristic idea behind the above conjecture is that the limit d 1 leads to positive weights in W 2 ( q ) , and, up to corrections in O ( 1 / d ) , the main argument leading to Equation (28) should apply for Equation (29). We corroborate this hypothesis by analyzing the dynamics for various values of q. The results are summarized in Figure 2.
For qubits ( d = 2 ), we first benchmark the self-averaging properties of the circuit, as shown in Figure 2a for q = 3 . The self-averaging holds similarly to q = 2 , i.e., the annealed S ˜ 3 and the quenched average S ¯ 3 of participation entropy rapidly approach each other with the increase in the circuit’s depth. We find analogous results for q = 1 , 4 and expect that a similar phenomenology will arise for any Rényi index q. We therefore limit our discussion to the difference between the saturation value and the annealed average of the participation entropy, i.e., S ˜ q H S ˜ q ( t ) .
In panels (b–d) of Figure 2, we demonstrate that the exponential decay of S ˜ q H S ˜ q ( t ) with the circuit depth t is observed for all considered values of the Rényi index q. For d = 2 and q = 3 (see Figure 2b), our tensor network implementation allows us to reach N 512 . Surprisingly, we find that the exponential decay of S ˜ q H S ˜ q ( t ) is very accurately fitted with α 2 , 3 N ( 2 K 2 ) t , even though we expected a decay β 2 , 3 t with a different coefficient β 2 , 3 fixed by the weights in | σ ˜ i , i + 1 , cf. Equation (29) and Section 4.
For larger values of q 4 , as well as for the von Neumann entropy limit q = 1 , we are limited to the relatively small system sizes N 24 accessible with the exact numerical simulation of the brick wall random circuits. The time evolution of S ˜ q H S ˜ q ( t ) for q = 1 , 4 approaches, with increasing system size N, the exponential decay ( 2 K 2 ) t . The separation of S ˜ q H S ˜ q ( t ) at a fixed time for different system sizes N suggests that the prefactor in front of the exponential decay ( 2 K 2 ) t may be scaling proportionally to the system size N. However, the range of available system sizes does not allow us to fully confirm this observation.
Finally, performing exact numerical calculations for d = 3 (up to N = 14 ) and d = 4 (up to N = 10 ), we verify that the decay of S ˜ q H S ˜ q ( t ) is accurately fitted with ( 2 K d ) t for q = 1 and q = 4 . This finding is analogous to the results for qubits ( d = 2 ) and corroborates further the observation that the base 2 K d characterizing the exponential decay of S ˜ q H S ˜ q ( t ) is independent of the Rényi index q.
Overall, the numerical results for the Rényi index q 2 presented in this section support our analysis, confirming that the approach of the participation entropy S ˜ q ( t ) to its long-time saturation value S ˜ q H occurs exponentially in time and with a prefactor that scales extensively with the number of qubits. Hence, also for q 2 , the approach, up to a fixed tolerance, of the participation entropy to its saturation value ε occurs at time t HSD ln ( N ) , which scales logarithmically with the system size N.

7. Discussion and Conclusions

We have investigated how the dynamics of random unitary circuits delocalize an initially localized state over a basis of the many-body Hilbert space. We combined analytical and numerical methods to calculate the time evolution of the participation entropies (1) that characterize the spread of state ρ = | Ψ Ψ | over a basis B of many-qudit systems under the dynamics of brick wall random circuits. Our main finding is that the process of Hilbert space delocalization occurs abruptly so that the long-time saturation value of participation entropies is approached up to a fixed tolerance at times t HSD ln ( N ) , scaling logarithmically with the system size N. These results may appear surprising from the point of view of the close relations between the participation entropy and the entanglement entropy [26,36,43,85]. While our analytical considerations show close parallels between entanglement entropy and participation entropy calculation, the boundary conditions of the relevant statistical mechanics model are different. This difference amounts to the uncovered contrasts between the time evolution of participation entropy and entanglement entropy. Indeed, given a connected subsystem, the entanglement entropy under local circuits grows linearly in time [49], saturating at times scaling proportionally to the system size N. However, when the subsystem is the unions of multiple disjoint intervals, the saturation time is much faster, of the order of the biggest domain configuration. These typical configurations are the leading contribution for the two-replica IPR computation and explain the logarithmic behavior of the saturation time in N.
We focused on times that were short compared to the system size. A complementary path considering random circuits in the limit of t and N and demonstrating the onset of universal features in the distributions of the overlaps between the output states generalizing the Porter–Thomas distribution was considered in Ref. [86]. We note that the IPR I 2 coincides with a collision probability which characterizes anticoncentration properties of many-body wave functions. The anticoncentration constitutes a necessary ingredient in formal arguments of the classical hardness of the sampling problems [87,88,89,90], and analytical bounds for the collision probability under random quantum circuit evolution were previously obtained in Refs. [91,92].
Random quantum circuits may be perceived as minimal models of local unitary dynamics. Hence, similar to the entanglement entropy case [93], we expect that our results about Hilbert space delocalization extend to generic chaotic non-integrable many-body systems. Moreover, our results may provide a relevant reference point for understanding the Hilbert space delocalization in non-ergodic many-body systems, cf. [94,95,96]. The statistical invariance of the considered circuits under a product of arbitrary on-site unitary transformations shows that our results about the growth of participation entropies hold for any basis B U obtained from the computational basis B by a product of independent on-site rotations.
The participation entropies considered here are not only directly available in numerical simulations of many-body systems but are also of experimental relevance. Indeed, recent progress in the stochastic sampling of many-body wave functions in ultracold atomic and solid-state experiments [97,98] allows, in principle, for their direct experimental evaluation. Similarly, the quantum processors realizing quantum circuits enable high-frequency sampling of the output state [99,100], which is a basis of the cross-entropy benchmarking [101], and could enable direct measuring of the participation entropies. Finally, the process of the estimation of the IPR, which, in general, requires resources scaling exponentially with the system size N, may be simplified by use of appropriate quantum algorithms [102].
This work has identified the leading terms relevant to the participation of entropy growth in random unitary dynamics. There are several interesting directions for further research. Analogously to the case of entanglement entropies [49,66], one may ask a question about fluctuations of the participation entropies around the identified mean values. Properties of the dynamics of participation entropies under higher dimensional circuits are another unresolved problem. Moreover, stabilizer Rényi entropies [103,104,105,106,107], which determine the number of beyond-classical (non-Clifford) operations needed to perform a quantum task, are related to participation entropies [83]. The long-time limit of stabilizer Rényi entropies under local unitary dynamics has been recently understood [108]. However, understanding the properties of the growth of stabilizer Rényi entropies is an exciting direction for further research facilitated by the results of this work.

Author Contributions

X.T. commenced the project and performed analytical calculations and tensor network contraction computations. P.S. performed exact numerical simulations. All authors have read and agreed to the published version of the manuscript.

Funding

X.T. acknowledges DFG under Germany’s Excellence Strategy–Cluster of Excellence Matter and Light for Quantum Computing (ML4Q) EXC 2004/1–390534769 and DFG Collaborative Research Center (CRC) 183 project no. 277101999—project B01. P.S. acknowledges support from European Research Council AdG NOQIA; MCIN/AEI (PGC2018-0910.13039/501100011033, CEX2019-000910-S/10.13039/501100011033, Plan National FIDEUA PID2019-106901GB-I00, Plan National STAMEENA PID2022-139099NB, I00, project funded by MCIN/AEI/10.13039/501100011033 and by the “European Union NextGenerationEU/PRTR" (PRTR-C17.I1), FPI); QUANTERA MAQS PCI2019-111828-2; QUANTERA DYNAMITE PCI2022-132919, QuantERA II Programme co-funded by European Union’s Horizon 2020 program under grant agreement no. 101017733); Ministry for Digital Transformation and of Civil Service of the Spanish Government through the QUANTUM ENIA project call—Quantum Spain project; and the European Union through the Recovery, Transformation and Resilience Plan—NextGenerationEU within the framework of the Digital Spain 2026 Agenda; Fundació Cellex; Fundació Mir-Puig; Generalitat de Catalunya (European Social Fund FEDER and CERCA program, AGAUR Grant No. 2021 SGR 01452, QuantumCAT U16-011424, co-funded by ERDF Operational Program of Catalonia 2014-2020); the computing resources at Urederra and technical support provided by NASERTIC (RES-FI-2024-1-0043); funded by the European Union. Views and opinions expressed are, however, those of the author(s) only and do not necessarily reflect those of the European Union, European Commission, European Climate, Infrastructure and Environment Executive Agency (CINEA), or any other granting authority. Neither the European Union nor any granting authority can be held responsible for them (EU Quantum Flagship PASQuanS2.1, 101113690, EU Horizon 2020 FET-OPEN OPTOlogic, grant no. 899794), EU Horizon Europe Program (this project has received funding from the European Union’s Horizon Europe research and innovation program under grant agreement no. 101080086 NeQSTGrant Agreement 101080086—NeQST); ICFO Internal “QuantumGaudi” project; European Union’s Horizon 2020 program under the Marie Sklodowska-Curie grant agreement no. 847648; “La Caixa” Junior Leaders fellowships, La Caixa” Foundation (ID 100010434): CF/BQ/PR23/11980043.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The code for the numerical simulations and the data are available at [77].

Acknowledgments

We thank A. De Luca and L. Piroli for enlightening discussions and collaborations on related topics. We are grateful to M. Ippoliti for pointing out the connection of our results to the anticoncentration measures in quantum information theory.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Schur-Weyl Duality and Integration over Haar Unitaries

In this Appendix, we briefly review the average over the q-tensor product of Haar distributed unitary U U ( d N ) . Here, U ( d N ) is the unitary group acting on N qudits with local Hilbert space dimension d. At the heart of our discussion is the averaged channel
Φ Haar ( q ) ( ) = E Haar U q ( ) ( U ) q d μ ( U ) U q ( ) ( U ) q ,
where μ ( U ) is the Haar measure over U ( d N ) , whose defining feature is the invariance under left group multiplication: μ ( U ) = μ ( V U ) for any V U ( d N ) . It is noteworthy that Φ Haar ( q ) ( ) commutes with all unitaries V q with V U ( d N ) , a fact following from
V q Φ Haar ( q ) ( ) ( V ) q = Haar d μ ( U ) V q U q ( ) ( U ) q ( V ) q = Haar d μ ( U ˜ ) U ˜ q ( ) ( U ˜ ) q = Φ Haar ( q ) ( ) .
By Schur–Weyl duality, it then can be expressed in terms of permutation operators
Φ Haar ( q ) ( ) = σ , τ S q Wg ( d N ; σ τ 1 ) tr T σ 1 ( ) T τ ,
where S q is the permutation group over q elements, and each permutation σ S k is represented as a replica operator T σ | n 1 , , n q = | n σ ( 1 ) , , n σ ( q ) for any basis { n j = 0 , , d N 1 } over the many-body Hilbert space H N = C d N . Without loss of generality, guaranteed by the Haar invariance of Φ Haar ( q ) ( ) , we can focus on the computational basis. Here, each n j is represented by the d-modular string | n j = | b j 1 , , b j N , with b j k = 0 , , d 1 . It is then easy to see that T σ = t σ N is reducible to permutation operators acting on each qudit digit, namely, t σ ( j ) | b j 1 , b j 2 , , b j q = | b j σ ( 1 ) , b j σ ( 2 ) , , b j σ ( q ) . (We will occasionally use the superscript ( j ) to clarify which qudit the action is on, stressing that the operator representation is the same.)
The expression in Equation (A3) was originally determined in [79], and Wg ( d N ; σ τ 1 ) are known as Weingarten coefficients. Denoting the projectors of the irreducible representations of the permutation group S k onto the D dimensional Hilbert space by Π λ , their dimensions by d λ , and their character functions by χ λ , the Weingarten function can be written as sums over the integer partitions of k), defined as follows (see Refs. [79,82] for further details):
Wg ( D ; σ ) = λ k d λ 2 ( k ! ) 2 χ λ ( σ ) tr ( Π λ ) .
Equation (A3) is then implicitly determined from Φ Haar ( q ) ( T τ ) = T τ , leading to linear equation
δ σ , τ = ω S q Wg ( d N ; σ ω 1 ) d N # ( ω 1 τ ) ,
where tr ( t ω t τ ) = d # ( ω τ ) , and # ( σ ) is the number of cycles of σ . Passing to the superoperator formalism, we have t τ ( j ) | τ j . Furthermore, since U q ( ) ( U ) q ( U U * ) q ( ) , it follows that Φ Haar ( q ) ( ) W N ( q ) , cf. Equation (10). There, the decomposition in terms of the dual operators is recast using
σ ˜ | 1 , , N = τ S q Wg ( d N ; σ τ 1 ) τ | 1 τ | N .
Figure A1. Annealed averaged Rényi-2 entanglement entropy for N A = N / 2 . Markers are numerical simulations obtained with tensor networks, and orange lines are the exact expression Equation (A11) in quantitative agreement.
Figure A1. Annealed averaged Rényi-2 entanglement entropy for N A = N / 2 . Markers are numerical simulations obtained with tensor networks, and orange lines are the exact expression Equation (A11) in quantitative agreement.
Entropy 26 00471 g0a1

Appendix B. Entanglement Growth for Finite Circuits

We detail the computation of the annealed averaged entanglement entropy. We consider the bipartition A = { 1 , , N A } and A c = { N A + 1 , , N } . In this case, the average purity of the reduced density matrix evolves as
P 2 ( t ) = S N A I ( N N A ) | W t | ψ 0 , , ψ 0 .
Recalling the structure in Equation (20), and that, upon absorption of σ | τ   r a n g l e = 2 # ( σ τ ) the W 2 ( 2 ) is
W 2 ( 2 ) = 1 K d K d 0 0 0 0 0 0 0 0 0 0 K d K d 1 ,
it follows that P 2 ( t ) is mapped to a random walk (of the entanglement cut), with absorbing boundaries at zero and N. Let us define the probability u z , t , where the random walk is absorbed exactly at time t when starting from the position z = N A . It follows that
u z , t + 1 = 1 2 u z 1 , t + 1 2 u z + 1 , t , u z , 0 = δ z , 0 + δ z , N , u 0 , t = u N , t = 1 .
Solving this stochastic process, we have
u z , t = ν = 0 N / 2 1 2 N sin π ( 2 ν + 1 ) N cos t 1 π ( 2 ν + 1 ) N sin π ( 2 ν + 1 ) z N .
Since each drift has an energy cost of 2 K , we have the final result
P 2 ( t ) = ( 2 K d ) t s = t + 1 u N A , s + s = 0 t ( 2 K d ) s u N A , s .
In particular, for N , we recast the findings in Ref. [55], and, for a sufficiently short time t N / 2 , we find that S 2 ent = v d t with the entanglement velocity v d = ln ( 2 K ) .

Appendix C. Details of the Numerical Implementation

This appendix discusses briefly the numerical methods employed in this manuscript.

Appendix C.1. Exact Computational Methods

For generic values of q and qubit ( d = 2 ) circuits, we evolve the exact representation of the state on the full Hilbert space. Each gate action translates to a sparse matrix multiplication, resulting in N qubits and M gates to O ( 2 N ) memory and O ( M 2 N ) runtime requirements. The knowledge of the state at each time step (circuit depth) allows us to compute I q for any value of q. We store the values of I q ( t ) and S q ( t ) for each realization of the circuit and average over N H = 10 4 realizations. This allows us to obtain both annealed and quenched averages of the participation entropies.

Appendix C.2. Tensor Network Simulations

For integer q, we can represent the average circuit in the tensor network formalism. As discussed in Section 4, the local degree of freedom is a “spin” of dimension d eff = q ! . Similarly, the W 2 ( q ) transfer matrix is a d eff × d eff gate. We store the initial state | Q in , cf. Section 5, in a matrix product state (MPS) [71]:
| Q in = { s j = 0 , , d eff 1 } A 1 , α 1 s 1 A α 1 , α 2 s 2 A α N 1 , 1 s N | s 1 , , s N .
Each layer in W t is an MPO, which we apply to the state, and compress with tolerance ε = 10 15 . The final contraction is with the state | Q fin . Our implementation is based on the open-source library ITensor [75].

References

  1. Zurek, W.H. Decoherence, einselection, and the quantum origins of the classical. Rev. Mod. Phys. 2003, 75, 715–775. [Google Scholar] [CrossRef]
  2. Nielsen, M.A.; Chuang, I.L. Quantum Computation and Quantum Information; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  3. Preskill, J. Quantum Computing in the NISQ era and beyond. Quantum 2018, 2, 79. [Google Scholar] [CrossRef]
  4. Amico, L.; Fazio, R.; Osterloh, A.; Vedral, V. Entanglement in many-body systems. Rev. Mod. Phys. 2008, 80, 517–576. [Google Scholar] [CrossRef]
  5. Laflorencie, N. Quantum entanglement in condensed matter systems. Phys. Rep. 2016, 646, 1–59. [Google Scholar] [CrossRef]
  6. D’Alessio, L.; Kafri, Y.; Polkovnikov, A.; Rigol, M. From quantum chaos and eigenstate thermalization to statistical mechanics and thermodynamics. Adv. Phys. 2016, 65, 239–362. [Google Scholar] [CrossRef]
  7. Chitambar, E.; Gour, G. Quantum resource theories. Rev. Mod. Phys. 2019, 91, 025001. [Google Scholar] [CrossRef]
  8. Kelly, S.P.; Poschinger, U.; Schmidt-Kaler, F.; Fisher, M.P.A.; Marino, J. Coherence requirements for quantum communication from hybrid circuit dynamics. SciPost Phys. 2023, 15, 250. [Google Scholar] [CrossRef]
  9. Stéphan, J.M.; Furukawa, S.; Misguich, G.; Pasquier, V. Shannon and entanglement entropies of one- and two-dimensional critical wave functions. Phys. Rev. B 2009, 80, 184421. [Google Scholar] [CrossRef]
  10. Stéphan, J.M.; Misguich, G.; Pasquier, V. Rényi entropy of a line in two-dimensional Ising models. Phys. Rev. B 2010, 82, 125455. [Google Scholar] [CrossRef]
  11. Alcaraz, F.C.; Rajabpour, M.A. Universal Behavior of the Shannon Mutual Information of Critical Quantum Chains. Phys. Rev. Lett. 2013, 111, 017201. [Google Scholar] [CrossRef]
  12. Luitz, D.J.; Plat, X.; Laflorencie, N.; Alet, F. Improving entanglement and thermodynamic Rényi entropy measurements in quantum Monte Carlo. Phys. Rev. B 2014, 90, 125105. [Google Scholar] [CrossRef]
  13. Laflorencie, N.; Rachel, S. Spin-resolved entanglement spectroscopy of critical spin chains and Luttinger liquids. J. Stat. Mech. Theory Exp. 2014, 2014, P11013. [Google Scholar] [CrossRef]
  14. Luitz, D.J.; Alet, F.; Laflorencie, N. Shannon-Rényi entropies and participation spectra across three-dimensional O(3) criticality. Phys. Rev. B 2014, 89, 165106. [Google Scholar] [CrossRef]
  15. Srednicki, M. Chaos and quantum thermalization. Phys. Rev. E 1994, 50, 888–901. [Google Scholar] [CrossRef] [PubMed]
  16. Deutsch, J.M. Quantum statistical mechanics in a closed system. Phys. Rev. A 1991, 43, 2046–2049. [Google Scholar] [CrossRef] [PubMed]
  17. Pappalardi, S.; Foini, L.; Kurchan, J. Eigenstate Thermalization Hypothesis and Free Probability. Phys. Rev. Lett. 2022, 129, 170603. [Google Scholar] [CrossRef] [PubMed]
  18. Fava, M.; Kurchan, J.; Pappalardi, S. Designs via Free Probability. arXiv 2023, arXiv:2308.06200. [Google Scholar]
  19. Pappalardi, S.; Fritzsch, F.; Prosen, T. General Eigenstate Thermalization via Free Cumulants in Quantum Lattice Systems. arXiv 2023, arXiv:2303.00713. [Google Scholar]
  20. Foini, L.; Dymarsky, A.; Pappalardi, S. Out-of-equilibrium Eigenstate Thermalization Hypothesis. arXiv 2024, arXiv:2406.04684. [Google Scholar]
  21. Bäcker, A.; Haque, M.; Khaymovich, I.M. Multifractal dimensions for random matrices, chaotic quantum maps, and many-body systems. Phys. Rev. E 2019, 100, 032117. [Google Scholar] [CrossRef]
  22. Pausch, L.; Carnio, E.G.; Rodríguez, A.; Buchleitner, A. Chaos and Ergodicity across the Energy Spectrum of Interacting Bosons. Phys. Rev. Lett. 2021, 126, 150601. [Google Scholar] [CrossRef]
  23. Fisher, M.P.; Khemani, V.; Nahum, A.; Vijay, S. Random Quantum Circuits. Annu. Rev. Condens. Matter Phys. 2023, 14, 335–379. [Google Scholar] [CrossRef]
  24. Potter, A.C.; Vasseur, R. Entanglement Dynamics in Hybrid Quantum Circuits. In Entanglement in Spin Chains; Quantum Science and Technology; Bayat, A., Bose, S., Johannesson, H., Eds.; Springer: Cham, Switzerland, 2022. [Google Scholar] [CrossRef]
  25. Lunt, O.; Richter, J.; Pal, A. Quantum Simulation Using Noisy Unitary Circuits and Measurements. In Entanglement in Spin Chains; Quantum Science and Technology; Bayat, A., Bose, S., Johannesson, H., Eds.; Springer: Cham, Switzerland, 2022. [Google Scholar] [CrossRef]
  26. Sierant, P.; Turkeshi, X. Universal Behavior beyond Multifractality of Wave Functions at Measurement-Induced Phase Transitions. Phys. Rev. Lett. 2022, 128, 130605. [Google Scholar] [CrossRef]
  27. Sierant, P.; Schirò, M.; Lewenstein, M.; Turkeshi, X. Measurement-induced phase transitions in (d+1)-dimensional stabilizer circuits. Phys. Rev. B 2022, 106, 214316. [Google Scholar] [CrossRef]
  28. Chahine, K.; Buchhold, M. Entanglement phases, localization and multifractality of monitored free fermions in two dimensions. arXiv 2023, arXiv:2309.12391. [Google Scholar]
  29. Abanin, D.A.; Altman, E.; Bloch, I.; Serbyn, M. Colloquium: Many-body localization, thermalization, and entanglement. Rev. Mod. Phys. 2019, 91, 021001. [Google Scholar] [CrossRef]
  30. Sierant, P.; Lewenstein, M.; Scardicchio, A.; Vidmar, L.; Zakrzewski, J. Many-Body Localization in the Age of Classical Computing. arXiv 2024, arXiv:2403.07111. [Google Scholar]
  31. Evers, F.; Mirlin, A.D. Anderson transitions. Phys. Rev. Lett. 2000, 84, 3690. [Google Scholar] [CrossRef] [PubMed]
  32. Evers, F.; Mirlin, A.D. Anderson transitions. Rev. Mod. Phys. 2008, 80, 1355. [Google Scholar] [CrossRef]
  33. Rodriguez, A.; Vasquez, L.J.; Römer, R.A. Multifractal Analysis with the Probability Density Function at the Three-Dimensional Anderson Transition. Phys. Rev. Lett. 2009, 102, 106406. [Google Scholar] [CrossRef]
  34. Rodriguez, A.; Vasquez, L.J.; Slevin, K.; Romer, R.A. Multifractal finite-size scaling and universality at the Anderson transition. Phys. Rev. Lett. 2010, 105, 046403. [Google Scholar] [CrossRef] [PubMed]
  35. Luca, A.D.; Scardicchio, A. Ergodicity breaking in a model showing many-body localization. Europhys. Lett. 2013, 101, 37003. [Google Scholar] [CrossRef]
  36. Macé, N.; Alet, F.; Laflorencie, N. Multifractal Scalings Across the Many-Body Localization Transition. Phys. Rev. Lett. 2019, 123, 180601. [Google Scholar] [CrossRef]
  37. Pietracaprina, F.; Laflorencie, N. Hilbert-space fragmentation, multifractality, and many-body localization. Ann. Phys. 2021, 435, 168502. [Google Scholar] [CrossRef]
  38. Altshuler, B.L.; Kravtsov, V.E.; Scardicchio, A.; Sierant, P.; Vanoni, C. Renormalization group for Anderson localization on high-dimensional lattices. arXiv 2024, arXiv:2403.01974. [Google Scholar]
  39. Scoquart, T.; Gornyi, I.V.; Mirlin, A.D. The role of Fock-space correlations in many-body localization. arXiv 2024, arXiv:2402.10123. [Google Scholar] [CrossRef]
  40. Vanoni, C.; Altshuler, B.L.; Kravtsov, V.E.; Scardicchio, A. Renormalization Group Analysis of the Anderson Model on Random Regular Graphs. arXiv 2023, arXiv:2306.14965. [Google Scholar]
  41. Sierant, P.; Lewenstein, M.; Scardicchio, A. Universality in Anderson localization on random graphs with varying connectivity. SciPost Phys. 2023, 15, 45. [Google Scholar] [CrossRef]
  42. Vanoni, C.; Vitale, V. An analysis of localization transitions using non-parametric unsupervised learning. arXiv 2024, arXiv:2311.16050. [Google Scholar]
  43. De Tomasi, G.; Khaymovich, I.M. Multifractality Meets Entanglement: Relation for Nonergodic Extended States. Phys. Rev. Lett. 2020, 124, 200602. [Google Scholar] [CrossRef]
  44. Styliaris, G.; Anand, N.; Campos Venuti, L.; Zanardi, P. Quantum coherence and the localization transition. Phys. Rev. B 2019, 100, 224204. [Google Scholar] [CrossRef]
  45. Anand, N.; Styliaris, G.; Kumari, M.; Zanardi, P. Quantum coherence as a signature of chaos. Phys. Rev. Res. 2021, 3, 023214. [Google Scholar] [CrossRef]
  46. Andreadakis, F.; Anand, N.; Zanardi, P. Scrambling of algebras in open quantum systems. Phys. Rev. A 2023, 107, 042217. [Google Scholar] [CrossRef]
  47. Baumgratz, T.; Cramer, M.; Plenio, M.B. Quantifying Coherence. Phys. Rev. Lett. 2014, 113, 140401. [Google Scholar] [CrossRef]
  48. Luitz, D.J.; Alet, F.; Laflorencie, N. Universal Behavior beyond Multifractality in Quantum Many-Body Systems. Phys. Rev. Lett. 2014, 112, 057203. [Google Scholar] [CrossRef] [PubMed]
  49. Nahum, A.; Ruhman, J.; Vijay, S.; Haah, J. Quantum Entanglement Growth under Random Unitary Dynamics. Phys. Rev. X 2017, 7, 031016. [Google Scholar] [CrossRef]
  50. Nahum, A.; Vijay, S.; Haah, J. Operator Spreading in Random Unitary Circuits. Phys. Rev. X 2018, 8, 021014. [Google Scholar] [CrossRef]
  51. Chan, A.; De Luca, A.; Chalker, J.T. Spectral Statistics in Spatially Extended Chaotic Quantum Many-Body Systems. Phys. Rev. Lett. 2018, 121, 060601. [Google Scholar] [CrossRef]
  52. Chan, A.; De Luca, A.; Chalker, J.T. Solution of a Minimal Model for Many-Body Quantum Chaos. Phys. Rev. X 2018, 8, 041019. [Google Scholar] [CrossRef]
  53. Shivam, S.; De Luca, A.; Huse, D.A.; Chan, A. Many-Body Quantum Chaos and Emergence of Ginibre Ensemble. Phys. Rev. Lett. 2023, 130, 140403. [Google Scholar] [CrossRef]
  54. Zhou, T.; Nahum, A. Emergent statistical mechanics of entanglement in random unitary circuits. Phys. Rev. B 2019, 99, 174205. [Google Scholar] [CrossRef]
  55. Zhou, T.; Nahum, A. Entanglement Membrane in Chaotic Many-Body Systems. Phys. Rev. X 2020, 10, 031066. [Google Scholar] [CrossRef]
  56. Bertini, B.; Piroli, L. Scrambling in random unitary circuits: Exact results. Phys. Rev. B 2020, 102, 064305. [Google Scholar] [CrossRef]
  57. von Keyserlingk, C.W.; Rakovszky, T.; Pollmann, F.; Sondhi, S.L. Operator Hydrodynamics, OTOCs, and Entanglement Growth in Systems without Conservation Laws. Phys. Rev. X 2018, 8, 021013. [Google Scholar] [CrossRef]
  58. Rakovszky, T.; Pollmann, F.; von Keyserlingk, C.W. Diffusive Hydrodynamics of Out-of-Time-Ordered Correlators with Charge Conservation. Phys. Rev. X 2018, 8, 031058. [Google Scholar] [CrossRef]
  59. Khemani, V.; Vishwanath, A.; Huse, D.A. Operator Spreading and the Emergence of Dissipative Hydrodynamics under Unitary Evolution with Conservation Laws. Phys. Rev. X 2018, 8, 031057. [Google Scholar] [CrossRef]
  60. Sünderhauf, C.; Piroli, L.; Qi, X.L.; Schuch, N.; Cirac, J.I. Quantum chaos in the Brownian SYK model with large finite N: OTOCs and tripartite information. J. High Energy Phys. 2019, 2019, 38. [Google Scholar] [CrossRef]
  61. Piroli, L.; Sünderhauf, C.; Qi, X.L. A random unitary circuit model for black hole evaporation. J. High Energy Phys. 2020, 2020. [Google Scholar] [CrossRef]
  62. Foligno, A.; Zhou, T.; Bertini, B. Temporal Entanglement in Chaotic Quantum Circuits. Phys. Rev. X 2023, 13, 041008. [Google Scholar] [CrossRef]
  63. Dowling, N.; Kos, P.; Modi, K. Scrambling Is Necessary but Not Sufficient for Chaos. Phys. Rev. Lett. 2023, 131, 180403. [Google Scholar] [CrossRef]
  64. Kos, P.; Bertini, B.; Prosen, T.c.v. Correlations in Perturbed Dual-Unitary Circuits: Efficient Path-Integral Formula. Phys. Rev. X 2021, 11, 011022. [Google Scholar] [CrossRef]
  65. Claeys, P.W.; Lamacraft, A. Ergodic and Nonergodic Dual-Unitary Quantum Circuits with Arbitrary Local Hilbert Space Dimension. Phys. Rev. Lett. 2021, 126, 100603. [Google Scholar] [CrossRef] [PubMed]
  66. Sierant, P.; Schirò, M.; Lewenstein, M.; Turkeshi, X. Entanglement Growth and Minimal Membranes in (d+1) Random Unitary Circuits. Phys. Rev. Lett. 2023, 131, 230403. [Google Scholar] [CrossRef]
  67. Rampp, M.A.; Claeys, P.W. Hayden-Preskill recovery in chaotic and integrable unitary circuit dynamics. arXiv 2023, arXiv:2312.03838. [Google Scholar]
  68. Fritzsch, F.; Ghosh, R.; Prosen, T. Boundary chaos: Exact entanglement dynamics. SciPost Phys. 2023, 15, 92. [Google Scholar] [CrossRef]
  69. Bulchandani, V.B.; Sondhi, S.L.; Chalker, J.T. Random-matrix models of monitored quantum circuits. J. Stat. Phys. 2024, 191, 55. [Google Scholar] [CrossRef]
  70. Luca, A.D.; Liu, C.; Nahum, A.; Zhou, T. Universality classes for purification in nonunitary quantum processes. arXiv 2023, arXiv:2312.17744. [Google Scholar]
  71. Schollwöck, U. The density-matrix renormalization group in the age of matrix product states. Ann. Phys. 2011, 326, 96–192. [Google Scholar] [CrossRef]
  72. Orús, R. Tensor networks for complex quantum systems. Nat. Rev. Phys. 2019, 1, 538–550. [Google Scholar] [CrossRef]
  73. Orús, R. A practical introduction to tensor networks: Matrix product states and projected entangled pair states. Ann. Phys. 2014, 349, 117–158. [Google Scholar] [CrossRef]
  74. Bañuls, M.C. Tensor Network Algorithms: A Route Map. Annu. Rev. Condens. Matter Phys. 2023, 14, 173–191. [Google Scholar] [CrossRef]
  75. Fishman, M.; White, S.R.; Stoudenmire, E.M. The ITensor Software Library for Tensor Network Calculations. arXiv 2020, arXiv:2007.14822. [Google Scholar] [CrossRef]
  76. Fishman, M.; White, S.R.; Stoudenmire, E.M. Codebase release 0.3 for ITensor. SciPost Phys. Codebases 2022, 4–r0.3. [Google Scholar] [CrossRef]
  77. Turkeshi, X.; Sierant, P. Code for “Hilbert space delocalisation under random unitary circuits”. arXiv 2024, arXiv:2404.10725. [Google Scholar] [CrossRef]
  78. Breuer, H.P.; Petruccione, F. The Theory of Open Quantum Systems; Oxford University Press on Demand: Oxford, UK, 2002. [Google Scholar]
  79. Weingarten, D. Asymptotic behavior of group integrals in the limit of infinite rank. J. Math. Phys. 1978, 19, 999–1001. [Google Scholar] [CrossRef]
  80. Collins, B.; Sniady, P. Integration with respect to the Haar measure on unitary, orthogonal and symplectic group. Commun. Math. Phys. 2006, 264, 773–795. [Google Scholar] [CrossRef]
  81. Köstenberger, G. Weingarten Calculus. arXiv 2021, arXiv:2101.00921. [Google Scholar]
  82. Collins, B.; Matsumoto, S.; Novak, J. The Weingarten Calculus. arXiv 2021, arXiv:2109.14890. [Google Scholar] [CrossRef]
  83. Turkeshi, X.; Schirò, M.; Sierant, P. Measuring nonstabilizerness via multifractal flatness. Phys. Rev. A 2023, 108, 042408. [Google Scholar] [CrossRef]
  84. Turkeshi, X.; Sierant, P. Error-Resilience Phase Transitions in Encoding-Decoding Quantum Circuits. Phys. Rev. Lett. 2024, 132, 140401. [Google Scholar] [CrossRef]
  85. Roy, S. Hilbert-space correlations beyond multifractality and bipartite entanglement in many-body localized systems. Phys. Rev. B 2022, 106, L140204. [Google Scholar] [CrossRef]
  86. Christopoulos, A.; Chan, A.; Luca, A.D. Universal distributions of overlaps from unitary dynamics in generic quantum many-body systems. arXiv 2024, arXiv:2404.10057. [Google Scholar]
  87. Aaronson, S.; Arkhipov, A. The computational complexity of linear optics. In Proceedings of the Forty-Third Annual ACM Symposium on Theory of Computing, San Jose, CA, USA, 6–8 June 2011; pp. 333–342. [Google Scholar]
  88. Bremner, M.J.; Montanaro, A.; Shepherd, D.J. Average-Case Complexity Versus Approximate Simulation of Commuting Quantum Computations. Phys. Rev. Lett. 2016, 117, 080501. [Google Scholar] [CrossRef] [PubMed]
  89. Bouland, A.; Fefferman, B.; Nirkhe, C.; Vazirani, U. On the complexity and verification of quantum random circuit sampling. Nat. Phys. 2019, 15, 159–163. [Google Scholar] [CrossRef]
  90. Oszmaniec, M.; Dangniam, N.; Morales, M.E.; Zimborás, Z. Fermion Sampling: A Robust Quantum Computational Advantage Scheme Using Fermionic Linear Optics and Magic Input States. PRX Quantum 2022, 3, 020328. [Google Scholar] [CrossRef]
  91. Dalzell, A.M.; Hunter-Jones, N.; Brandão, F.G.S.L. Random Quantum Circuits Anticoncentrate in Log Depth. PRX Quantum 2022, 3, 010333. [Google Scholar] [CrossRef]
  92. Bertoni, C.; Haferkamp, J.; Hinsche, M.; Ioannou, M.; Eisert, J.; Pashayan, H. Shallow shadows: Expectation estimation using low-depth random Clifford circuits. arXiv 2023, arXiv:2209.12924. [Google Scholar]
  93. Kim, H.; Huse, D.A. Ballistic Spreading of Entanglement in a Diffusive Nonintegrable System. Phys. Rev. Lett. 2013, 111, 127205. [Google Scholar] [CrossRef] [PubMed]
  94. Torres-Herrera, E.J.; García-García, A.M.; Santos, L.F. Generic dynamical features of quenched interacting quantum systems: Survival probability, density imbalance, and out-of-time-ordered correlator. Phys. Rev. B 2018, 97, 060303. [Google Scholar] [CrossRef]
  95. Creed, I.; Logan, D.E.; Roy, S. Probability transport on the Fock space of a disordered quantum spin chain. Phys. Rev. B 2023, 107, 094206. [Google Scholar] [CrossRef]
  96. Hopjan, M.; Vidmar, L. Scale-Invariant Survival Probability at Eigenstate Transitions. Phys. Rev. Lett. 2023, 131, 060404. [Google Scholar] [CrossRef]
  97. Brydges, T.; Elben, A.; Jurcevic, P.; Vermersch, B.; Maier, C.; Lanyon, B.P.; Zoller, P.; Blatt, R.; Roos, C.F. Probing Rényi entanglement entropy via randomized measurements. Science 2019, 364, 260–263. [Google Scholar] [CrossRef] [PubMed]
  98. Ebadi, S.; Wang, T.T.; Levine, H.; Keesling, A.; Semeghini, G.; Omran, A.; Bluvstein, D.; Samajdar, R.; Pichler, H.; Ho, W.W.; et al. Quantum phases of matter on a 256-atom programmable quantum simulator. Nature 2021, 595, 227–232. [Google Scholar] [CrossRef] [PubMed]
  99. Arute, F.; Arya, K.; Babbush, R.; Bacon, D.; Bardin, J.C.; Barends, R.; Biswas, R.; Boixo, S.; Brandao, F.G.S.L.; Buell, D.A.; et al. Quantum supremacy using a programmable superconducting processor. Nature 2019, 574, 505–510. [Google Scholar] [CrossRef] [PubMed]
  100. Wu, Y.; Bao, W.S.; Cao, S.; Chen, F.; Chen, M.C.; Chen, X.; Chung, T.H.; Deng, H.; Du, Y.; Fan, D.; et al. Strong Quantum Computational Advantage Using a Superconducting Quantum Processor. Phys. Rev. Lett. 2021, 127, 180501. [Google Scholar] [CrossRef] [PubMed]
  101. Boixo, S.; Isakov, S.V.; Smelyanskiy, V.N.; Babbush, R.; Ding, N.; Jiang, Z.; Bremner, M.J.; Martinis, J.M.; Neven, H. Characterizing quantum supremacy in near-term devices. Nat. Phys. 2018, 14, 595–600. [Google Scholar] [CrossRef]
  102. Liu, Y.; Sierant, P.; Stornati, P.; Lewenstein, M.; Płodzień, M. Quantum Algorithms for Inverse Participation Ratio Estimation in multi-qubit and multi-qudit systems. arXiv 2024, arXiv:2405.03338. [Google Scholar]
  103. Leone, L.; Oliviero, S.F.E.; Hamma, A. Stabilizer Rényi Entropy. Phys. Rev. Lett. 2022, 128, 050402. [Google Scholar] [CrossRef] [PubMed]
  104. Haug, T.; Lee, S.; Kim, M.S. Efficient quantum algorithms for stabilizer entropies. arXiv 2023, arXiv:2305.19152. [Google Scholar] [CrossRef]
  105. Oliviero, S.F.E.; Leone, L.; Hamma, A.; Lloyd, S. Measuring magic on a quantum processor. NPJ Quantum Inf. 2022, 8, 148. [Google Scholar] [CrossRef]
  106. Gu, A.; Leone, L.; Ghosh, S.; Eisert, J.; Yelin, S.; Quek, Y. A little magic means a lot. arXiv 2023, arXiv:2308.16228. [Google Scholar]
  107. Tirrito, E.; Tarabunga, P.S.; Lami, G.; Chanda, T.; Leone, L.; Oliviero, S.F.E.; Dalmonte, M.; Collura, M.; Hamma, A. Quantifying non-stabilizerness through entanglement spectrum flatness. Phys. Rev. Lett. 2024, 109, L040401. [Google Scholar] [CrossRef]
  108. Turkeshi, X.; Dymarsky, A.; Sierant, P. Pauli Spectrum and Magic of Typical Quantum Many-Body States. arXiv 2023, arXiv:2312.11631. [Google Scholar]
Figure 1. Participation entropy for the Rényi index q = 2 . The left panel shows the quenched S ¯ 2 (in red) and the annealed S 2 ˜ (in blue) average of the participation entropy rescaled by the system size N and plotted as a function of circuit depth t. The inset shows the exponential in time decay of the difference between the average participation entropy S ˜ 2 H for random Haar states and the averages S ¯ 2 , S ˜ 2 as a function of time. The right panel shows the approach to S ˜ 2 H of the annealed average S 2 ˜ (blue symbols) and the exact solution obtained with the tensor network approach (solid purple lines) compared with the asymptotic formula S ˜ 2 H S ˜ 2 ( t ) = α d N ( 2 K d ) t 1 (red dashed line), where, for qubits, d = 2 , α 2 0.291 ( 5 ) and K 2 = 2 5 . The inset shows the prefactor α d as function of the on-site Hilbert space dimension d and the α d extrapolated from tensor network simulations, while α d 1 2 is shown by the black dashed line.
Figure 1. Participation entropy for the Rényi index q = 2 . The left panel shows the quenched S ¯ 2 (in red) and the annealed S 2 ˜ (in blue) average of the participation entropy rescaled by the system size N and plotted as a function of circuit depth t. The inset shows the exponential in time decay of the difference between the average participation entropy S ˜ 2 H for random Haar states and the averages S ¯ 2 , S ˜ 2 as a function of time. The right panel shows the approach to S ˜ 2 H of the annealed average S 2 ˜ (blue symbols) and the exact solution obtained with the tensor network approach (solid purple lines) compared with the asymptotic formula S ˜ 2 H S ˜ 2 ( t ) = α d N ( 2 K d ) t 1 (red dashed line), where, for qubits, d = 2 , α 2 0.291 ( 5 ) and K 2 = 2 5 . The inset shows the prefactor α d as function of the on-site Hilbert space dimension d and the α d extrapolated from tensor network simulations, while α d 1 2 is shown by the black dashed line.
Entropy 26 00471 g001
Figure 2. Participation entropy growth for various Rényi indices q. Panel (a) shows the approach of the quenched S ¯ 3 (in red) and annealed S 3 ˜ (in blue) average of the participation entropy to the random Haar state value S ˜ 3 H as a function of the circuit depth t. Panel (b) presents the approach to S ˜ 3 H of the annealed average S 3 ˜ (blue symbols) calculated with the exact numerical simulation compared with the tensor network contraction results (purple lines) and the fitted formula S ˜ 3 H S ˜ 3 ( t ) 0.76 N ( 2 K d ) t 1 (red dashed line) consistent with Equation (29). Panels (c,d) show the approach of the circuit-averaged participation entropy S ¯ q to the random Haar state value S ˜ q H , respectively, for q = 1 and q = 4 .
Figure 2. Participation entropy growth for various Rényi indices q. Panel (a) shows the approach of the quenched S ¯ 3 (in red) and annealed S 3 ˜ (in blue) average of the participation entropy to the random Haar state value S ˜ 3 H as a function of the circuit depth t. Panel (b) presents the approach to S ˜ 3 H of the annealed average S 3 ˜ (blue symbols) calculated with the exact numerical simulation compared with the tensor network contraction results (purple lines) and the fitted formula S ˜ 3 H S ˜ 3 ( t ) 0.76 N ( 2 K d ) t 1 (red dashed line) consistent with Equation (29). Panels (c,d) show the approach of the circuit-averaged participation entropy S ¯ q to the random Haar state value S ˜ q H , respectively, for q = 1 and q = 4 .
Entropy 26 00471 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Turkeshi, X.; Sierant, P. Hilbert Space Delocalization under Random Unitary Circuits. Entropy 2024, 26, 471. https://doi.org/10.3390/e26060471

AMA Style

Turkeshi X, Sierant P. Hilbert Space Delocalization under Random Unitary Circuits. Entropy. 2024; 26(6):471. https://doi.org/10.3390/e26060471

Chicago/Turabian Style

Turkeshi, Xhek, and Piotr Sierant. 2024. "Hilbert Space Delocalization under Random Unitary Circuits" Entropy 26, no. 6: 471. https://doi.org/10.3390/e26060471

APA Style

Turkeshi, X., & Sierant, P. (2024). Hilbert Space Delocalization under Random Unitary Circuits. Entropy, 26(6), 471. https://doi.org/10.3390/e26060471

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop