arXiv:2108.13143v1 [cond-mat.quant-gas] 30 Aug 2021 ri eleprmn)teprobability the experiment) real a in or state h sn oe ihfroantcnaetniho pnsi co spin-spin nearest-neighbor syste ferromagnetic underlying with the model of breaking Ising symmetr symmetry Moreover, the spontaneous systems. and these tries of behavior macroscopic the overlaps n rniinpoaiiis nmr oceetrs o ie tim given a for terms, a concrete to more [15], In Ref. in o probabilities. ensemble described transition random previously and was the which of approach, instead this model, employ statistical the of definition rniin nmr eea em,i sipratt nesadt understand to important is it terms, symmetry. sm this general a arbitrarily near near more an dynamics already In complex that a transition. fact indicates The T degeneracy macroscopic symmetry. states. the the quantum pe the breaks two small that of the arbitrarily conditions) of symmetry an boundary the add special we or of field when consequence breaking a symmetry is spontaneous degeneracy two-fold This hc ersnsrno ntr arcsadhsbe sda t [10–14]. R´enyi a entropy related as cir the used and Dyson’s been matrix on has density based and the is matrices of unitary It random measurements. transformation represents the random which of symmetry of application to Another description respect the transformations. Th with symplectic systems. o or ensemble quantum ensemble unitary of kn random random class of the a a way of Hamiltonian, of features atom no specific chara generic an some a is the have as describes guessing there though, such of that complex instead system, systems, is Thus, these quantum such RMT of many-body systems, the spectra small the physical for relatively many motivation a to successf The a even very applied data a been [1–7]. found experimental has and systems of new It sets not (RMT). large is Theory idea that q This fact averaged statistically. the using treated f by evolution, sugge behavior quantum a supported a the for also Such about Even walk. random transit information classical o dynamics. and generic a freedom complex to return similar of a the qua erratic degrees have quite are the or systems exceptions about quantum properties Typical information Moreover, all of experiment. not extraction the words, by other the In that idea limited. the on based nytetm famaueeti adm hra h nrylevels energy the whereas random, is measurement a of time the only ntefloigw ilcnie h ntr vlto fcoe quantu closed of evolution unitary the consider will we following the In ncnrs oteeRTapoce,w osdri h olwn a following the in consider we approaches, RMT these to contrast In ymtispa eta oei lsia swl si unu many- quantum in as well as classical in role central a play Symmetries | oeho ucin trvasfrtelte characterist latter conc the general for a localization. reveals as space It discussed is junction. matrix Josephson correlation qua time the smalle Probing the the degeneracies. of by spectral the determined pres at is diverges the rate rate in decay instance, slowest extract for the p to Moreover, decaying occur, us exponentially degeneracies enables an spectral and This when part times. recurrent random a evolution, at measured are that Ψ i h ignleeet ftetm orlto arxaeuse are matrix correlation time the of elements diagonal The h E r o adm hslast h iecreainmti TM sth as (TCM) matrix correlation time the to leads This random. not are j | Ψ rbn aybd ytm erseta degeneracies spectral near systems many-body Probing 0 i , h ntttfu hsk nvri¨tAgbr,D815Augs D-86135 Universit¨at Augsburg, f¨ur Physik, Institut E j | Ψ i fteeeg eigenstates energy the of .INTRODUCTION I. Dtd uut3,2021) 31, August (Dated: p k httesse si eti tt.Te eevaluate we Then state. certain a in is system the that .Ziegler K. {| E j i} ihagvniiilstate initial given a with cpoete ttetasto oHilbert- to transition the at properties ic tmeouinwt h ignlelements diagonal the with evolution ntum neo pnaeu ymtybreaking. symmetry spontaneous of ence r.Ti eaaini togyaffected strongly is separation This art. p n etdi h aeo bosonic a of case the in tested and ept teeg ee pcn,adti decay this and spacing, level energy st trsi etrs uha ee repulsion. level as such features, cteristic ls scaatrzdb h invariance the by characterized is class e e fmcocpcsae eetsymme- reflect states macroscopic of ies opoecoe unu systems quantum closed probe to d w itntprso h quantum the of parts distinct two i ed oalfigo h degeneracy the of lifting a to leads his .Frisac,tegonsae of groundstates the instance, For m. pigare upling eaalbewoepoete a be can properties whose available re o rbblte o unu states. quantum for probabilities ion M a enrcnl rpsdfor proposed recently been has RMT lraiaini h admMatrix Random the in realization ul o odtrietetaeo powers of trace the determine to ool aitnasi h M.W will We RMT. the in Hamiltonians f l ymtybekn emcnlift can term breaking symmetry all eeouino aybd systems many-body of evolution he h unu oe r accessible are model quantum the f ymtyadteascae phase associated the and symmetry .Teeaetpclyorthogonal, typically are These s. atte.Asaitclapoc is approach statistical A uantities. { ua arxesml 3 ,8 9]), 8, 6, [3, ensemble matrix cular aitna.Te eexperience we Then Hamiltonian. t ttsia praht extract to approach statistical a sts E tmsse na xeietis experiment an in system ntum train(.. ml magnetic small a (e.g., rturbation e wprilsteeouincnbe can evolution the particles ew ayeteeouino h return the of evolution the nalyze snce,aosadmesoscopic and atoms nuclei, as j ug Germany burg, cnces nteohrhand, other the On nucleus. ic t aitnasi hsn which chosen, is Hamiltonians f } k oysses hydetermine They systems. body yaia prahi which in approach dynamical fteHamiltonian the of eeaut i calculation a (in evaluate we aybd ytm.I is It systems. many-body m wn h aitna of Hamiltonian the owing ↑ | , . . . , eta olfrthe for tool central e | Ψ 0 ↑i i n measured a and and ↓ | H , . . . , n the and ↓i . 2 the probabilities {p1,p2, ...} at different discrete and randomly chosen times {t1,t2, ...}. This can be translated into practical observations, where it is assumed that each experiment is prepared in the same initial state and all measurements are performed for the same final state of the evolution at different times. These experiments provide an ensemble of probabilities {p1,p2, ...} with the corresponding times {t1,t2, ...}. For given overlaps hEj |Ψ0i, hEj |Ψi we can immediately predict some restrictions for the evolution of the probability pk in the N–dimensional Hilbert space. When the overlaps vanish for some of the eigenstates |Eni, the evolution cannot reach those states and the accessible Hilbert space is restricted to the states |Ej i with j 6= n. This reduction of the Hilbert space can be interpreted as Hilbert-space localization [16] or Hilbert-space fragmentation [17]. This effect can be associated with spontaneous symmetry breaking, induced by the choice of the initial and measured states. In the case that the overlap to some states is not strictly zero but very small, the access to those states may be negligible and can be ignored. This corresponds to a complex dynamical behavior and requires a careful analysis. It will be addressed briefly for the example of a bosonic Josephson junction in Sect. III, where the mirror symmetry of the junction is spontaneously broken. This paper is organized as follows. After the definition of the TCM in Sect. II we focus on the properties of its diagonal elements (Sect. II A). Then the effect of spectral degeneracies on the diagonal TCM elements are discussed in Sect. IIB. In Sect. III we analyze the diagonal TCM elements in the specific example of a bosonic Josephson junction.

II. TIME CORRELATION MATRIX

We consider the transition amplitude from |Ψ0i to |Ψi

−iHtk uk = hΨ|e |Ψ0i , (1) which is based on the unitary evolution with the Hamiltonian H from the initial state |Ψ0i. The prob- −iHtk 2 2 ability to measure the state |ψi at time tk is given by |hΨ|e |Ψ0i| . In other words, pk = |uk| is the probability to find the quantum system in the state |Ψi after the unitary evolution from the initial state |Ψ0i over the time tk. Since the evolution is defined by the Hamiltonian H, we consider its eigen- states {|Ej i}j=1,...,N and its corresponding eigenvalues {Ej}j=1,...,N and write the amplitude in spectral representation as

N N −iEj tk −iEj tk uk = hΨ|Ej ihEj |Ψ0ie ≡ qj e . (2) j=1 j=1 X X Although the phases are not directly accessible experimentally, their effect is observable through the interference of quantum states. This enables us, for instance, to relate the product of amplitudes at different times with probabilities of interfering amplitudes:

∗ ∗ 1 2 2 ∗ ∗ 1 2 2 u u ′ + u ′ u = |u + u ′ | − |u − u ′ | , i(u u ′ − u ′ u )= |u + iu ′ | − |u − iu ′ | , (3) k k k k 2 k k k k k k k k 2 k k k k

where the probabilities are observable in interferometric measurements. This relation suggests to consider ′ ∗ ′ the correlation of the amplitudes uk, uk at different times through the TCM hukuk iτ , where the average h...iτ is taken with respect to the distribution of times {tk} as a result of inaccurate clocks: The time is measured by a clock in each laboratory, which counts time steps {τn}. These clocks have a limited accuracy, such that the time steps vary randomly. This implies a sequence of measurements in each laboratory, where the clock indicates k time steps corresponding to the total evolution time tk = τ1 + ··· + τk for different values of k. Now we compare the measurements of different laboratories. This ∗ ′ provides a distribution of results for ukuk due to different inaccurate clocks, where we assume that the fluctuations of the time steps {τn} are independently and equally distributed. Then the TCM reads in spectral representation

′ ′ ∗ ′ ∗ ′ iEj (τ1+···+τk) −iEj (τ1+···+τk ) hukuk iτ = qj qj he e iτ (4) ′ Xj,j 3

i(Ej −E ′ )(τ1+···+τk) −iE ′ (τk+1+···+τ ′ ) ′ he j iτ he j k iτ k > k ∗ i(E −E ′ )(τ +···+τ ′ ) iE ′ (τ ′ +···+τ ) ′ = q qj′ he j j 1 k i he j k +1 k i k < k . (5) j  τ τ ′ i(Ej −Ej′ )(τ1+···+τk) ′ Xj,j  he iτ k = k

iEj τ i(Ej −E ′ )τ Defining λj = he iτ and λjj′ = he j iτ the TCM elements become ′ k k −k ′ ∗ ∗ λjj′ λj′ k ≥ k hu u ′ i = q q ′ ′ ′ . (6) k k τ j j k ∗k−k ′ ′ ( λjj′ λj k < k Xj,j ′ ′ The TCM decays exponentially with |k − k |, provided |λj | < 1. For fixed |k − k | the TCM is constant for the diagonal elements λjj = 1, though. This reflects the fact that a unitary evolution between the same energy eigenstates gives just a phase factor e−iEj τ (cf. Eq. (2)). On the other hand, these phase factors lead to a decay of different energy states due to interference effects after the time average.

A. Diagonal elements of the TCM

2 The diagonal TCM element h|uk| iτ is the probability to measure the state |Ψi at time tk. Before time 2 averaging, the expression |uk| is a diagonal element of the density matrix ρ(tk) with respect to the state 2 |Ψi. The trace of |uk| with respect to all states |Ψi of the underlying Hilbert space is the spectral form factor, often used for the characterization of many-body quantum chaos [18, 19]. We only mention this but will not study it here. 2 According to Eq. (6) the average transition probability h|uk| iτ reads N N

2 ∗ ′ k 2 h|uk| iτ = PN + qj qj λjj′ , PN = |qj | . (7) ′ ′6 j=1 j,j =1;Xj =j X The term PN describes the recurrent behavior, which does not depend on time. It is the asymptotic transition probability for k → ∞ 2 PN = lim h|uk| iτ , (8) k→∞ provided that the energy levels are not degenerated. The case of degenerate energy levels is discussed in the next section. The second term in Eq. (7) decays exponentially with time due to |λjj′ | < 1, and only this term describes a change of the transition probability during the evolution of the quantum system. This result provides a separation of the diagonal elements of the TCM into a static recurrent term PN and a dynamic term that decays quickly. The recurrent term PN stores important information regarding the properties of the quantum system. 2 2 2 Since |qj | = |hΨ|Ej i| |hΨ0|Ej i| is a product of the overlaps between the energy eigenstate |Ej i with the initial state and with the measured state, it provides a measure of how much this energy eigenstate con- tributes to the transition |Ψ0i → |Ψi during the unitary evolution. For instance, the asymptotic behavior of the return probability to the initial state |Ψ0i → |Ψ0i with the dimensionality N of the underlying Hilbert space describes Anderson localization when limN→∞ PN > 0 and the absence of Anderson local- ization when limN→∞ PN = 0 [16]. This can be understood by noting that the normalization of quantum N 2 states implies j=1 |hΨ0|Ej i| = 1 and that for a localized state only a few energy eigenstates have a nonzero overlap with |Ψ0i. For a delocalized state, on the other hand, the overlap is nonzero for a large number of energyP eigenstates which is of the order of N. An extreme case is given when these overlaps 2 are equal. Then we have |hΨ0|Ej i| =1/N due to the normalization, which implies PN =1/N. Anderson localization is associated with a random Hamiltonian [20]. According to the above described picture, we can also consider Hilbert-space localization for a deterministic Hamiltonian, which depends strongly on the initial state. For an energy eigenstate, the system will always remain in the latter under unitary evo- lution. More general, is the initial state a superposition of m energy eigenstates, the system will always remain inside the m–dimensional Hilbert space, spanned by these energy eigenstates. In the case, where the initial state is eigenstate of H0 of the Hamiltonian H = H0 + ηH1 and η is a small parameter, ηH1 is a small perturbation. In that case it is possible that this perturbation provides an exponentially decaying evolution away the initial state. This would be considered as exponential Hilbert-space localization. In the subsequent discussion we will focus on the diagonal elements of the TCM, since the off-diagonal TCM elements decay exponentially with |k − k′| according to Eq. (6). 4

B. Effect of spectral degeneracies

Assuming that there is a spectral degeneracy E1 = E2, we have λ12 = λ21 = 1, and the diagonal TCM elements in Eq. (7) read in this case

N N

2 ∗ ′ k ∗ ∗ ∗ ′ k h|uk| iτ = qj qj λjj′ = PN + q1 q2 + q2 q1 + qj qj λjj′ , (9) ′ ′ ′ ′ j,jX=1 j,j =1;j =6 j;(Xj,j )=(16 ,2),(2,1) such that the recurrent part of the transition probability becomes

N 2 ∗ ∗ 2 2 lim h|uk| iτ = PN + q1 q2 + q2 q1 = |q1 + q2| + |qj | . (10) k→∞ j=3 X Thus, the effect of a spectral degeneracy is a change of the recurrent and the decaying behavior, where 2 2 2 the recurrent term changes by |q1 + q2| − |q1| − |q2| . This means that the diagonal elements of the TCM are very sensitive in terms of spectral degeneracies. 2 After applying a discrete Fourier transformation to the decaying part of h|uk| iτ we obtain the function

N N λjj′ ˜ iω iωk ∗ ′ k ∗ ′ Ud(e )= e qj qj λjj′ = qj qj , (11) e−iω − λ ′ ≥ ′ ′ ′ ′ ′ ′ jj kX1 j,j =1;k =6 k;(Xj,j )=(16 ,2),(2,1) j,j =1;j =6 j;(Xj,j )=(16 ,2),(2,1)

which is a function of ω on the interval [0, 2π). In other words, U˜d(z) is a sum of poles inside the unit circle due to |λjj′ | < 1. The poles λnn′ and λn′n approach the unit circle when we are getting iω closer to a degeneracy of En and En′ This should be visible in U˜d(e ). The corresponding decay time Td = −1/ log |λnn′ | diverges due to |λnn′ |∼ 1. Therefore, the decay time Td is a measure for the distance from a spectral degeneracy; it diverges when we approach the degeneracy. In general, we can define 1 Tm =′ max − (12) j,j =1,...,N log |λjj′ | as the largest decay time as a measure of level degeneracy.

III. EXAMPLE: BOSONIC JOSEPHSON JUNCTION

In this section we will study the diagonal TCM elements of a bosonic Josephson junction (BJJ) with N bosons as a closed quantum system. The motivation for choosing this example is at least threefold: The model is (i) simple enough but not trivial with interesting features based on tunneling and boson- boson interaction, (ii) it can be solved exactly and (iii) it has been realized experimentally [21, 22] with applications to commercial quantum computers [23]. The BJJ consists of two identical wells filled with interacting bosons and a tunneling junction between them. More formally, it is defined by the Bose-Hubbard Hamiltonian [24]

J H = − (a†a + a†a )+ U(n2 + n2), n = a† a , (13) 2 l r r l l r l,r l,r l,r

† where al,r (al,r) are the creation (annihilation) operators in the left and right well, respectively. The first term of H describes tunneling of atoms between the wells, and for U > 0 the second term represents a repulsive particle-particle interaction that favors energetically a symmetric distribution of bosons in the double well. Without tunneling (i.e., for J = 0) there are two-fold degenerate energy levels Ek = U[(N − k)2 + k2]/2 with eigenstates that are superpositions of the product Fock state |k,N − ki (≡ |ki⊗ |N − ki) and its mirror image |N −k, ki. This two-fold degeneracy is similar to the two-fold degeneracy of the Ising model, mentioned in the Introduction. Thus, tunneling between the wells plays the role of the symmetry- breaking term. In contrast to the macroscopic Ising model though, an arbitrarily small tunneling term may not be sufficient to cause symmetry breaking. Hence the following study is more related to an 5

Ising model with a finite number of spins than to a macroscopic system. Nevertheless, the sensitive dynamics due to tunneling between the two wells with degenerate energy levels will give us an insight into the evolution near symmetry points. Another difference between the Ising model and the BJJ is that without tunneling in the latter all energy levels are two-fold degenerate. Therefore, the initial state can be prepared in any of these energy levels to follow the evolution due to tunneling in the vicinity of degenerate levels. This is important because the degenerate groundstate may not be reached due to the energy conservation in the closed quantum system. For the following we will use |Ψ0i = |0,Ni as the initial state and |Ψi = |N, 0i as the measured state. Then we define the return probability (RP) and the transition probability (TP) as

2 −iHtk 2 2 −iHtk 2 |ur,k| = |h0,N|e |0,Ni| , |ut,k| = |hN, 0|e |0,Ni| . (14)

Since both states |Ψ0i, |Ψi are eigenstates of H in the absence of tunneling (J = 0), we get

2 2 |ur,k| =1 , |ut,k| = 0 (15) for any k due to orthogonality. This reflects the fact that this pair of states breaks the mirror symmetry of the double well. The opposite extreme is the BJJ without boson-boson interaction (U = 0), which is more complex and will be discussed in the next subsection. As we will see, this case can be described by 2 2 simple functions for |ur,k| , |ut,k| . For the interplay of tunneling and boson-boson interaction (J, U 6= 0) 2 2 the behavior is more complex and we will rely on the time averaged expressions h|ur,k| iτ , h|ut,k| iτ with ′ ′ τk =τ ¯ + τk and an exponential distribution for τk:

∞ ′ −τn ′ h...iτ = ... e dτn , (16) 0 ≥ Z nY1 where the time is measured in units of ¯h/J.

A. Non-interacting bosons

For U = 0 the spectrum of H consists of equidistant energy levels Ej = −J(N/2 − j) (j =0, 1, ..., N) and eigenstates

−N/2 2 † † j † † N−j |Ej i = (a + ar) (a − ar) |0, 0i . j!(N − j)! l l

Then the RP and the TP at time tpk for N bosons read

2 −iHtk 2 N 2 2 −iHtk 2 N 2 |ur,k| = |h0,N|e |0,Ni| = | cos (Jtk/2)| , |ut,k| = |hN, 0|e |0,Ni| = | sin (Jtk/2)| . (17) ′ ′ For the average TP we get with tk = kτ¯ + τ1 + ··· + τk and the exponential distribution of Eq. (16)

N ′ k ′ −iJτ¯(l−l ) 2 −2N N N l+l e h|ut,k| iτ =2 ′ (−1) ′ ′ l l "1 − iJ(l − l )# l,lX=0   

′ k N 2 N −iJτ¯(l−l ) N N N ′ e =2−2N +2−2N (−1)l+l , (18) l l l′ 1 − iJ(l − l′) ′ ′6 " # Xl=0   l,l =0;Xl =l    ′ while for the average RP is the same without the factor (−1)l+l . From these results we get for the asymptotic value at k ∼ ∞ for both probabilities

N N 2 P ∼ h|u |2i ∼ h|u |2i ∼ 2−2N , (19) N r,k τ t,k τ l Xl=0   which reflects the mirror symmetry of the BJJ. 6

1.1 u=1 1 u=2 0.9 0.8 0.7 0.6 0.5 0.4 return probability 0.3 0.2 0.1 0 0 2 4 6 8 10 12 14 16 time 0.5 u=1 0.45 u=2

0.4

0.35

0.3

0.25

0.2

transition probability 0.15

0.1

0.05

0 0 2 4 6 8 10 12 14 16 time

FIG. 1. Signatures of a qualitative change of the evolution of a bosonic Josephson junction upon an increasing interaction strength u. The plots represent the dynamics of 20 bosons for u = 1, 2, where the top panel gives the 2 2 return probability |ur,k| and the bottom panel the transition probability |ut,k| .

B. Interacting bosons

Exact solutions exist for this model also for J, U 6= 0 but in contrast to the non-interacting case they are complex and difficult to present in general. For instance, the resolvent is a meromorphic function with polynomials of order N and N + 1 [16, 25]. Therefore, we only plot the results for the RP and the TP and their averaged counterparts here. For transparency, we choose for all subsequent plots N = 20 bosons. First, the evolution of the RP and the TP are presented in Fig. 1 for two values u = 1, 2 of the interaction parameter u = NU/J. This clearly reveals that the RP dominates over the TP for increasing u, as we expect from the results of the two limits J = 0 and U = 0. It is interesting to note that in a mean-field (classical) approximation of the BJJ there is a sharp phase transition in terms of the interaction parameter, where the mean-field TP is completely suppressed when u ≥ uc = 2 [26]. The strong interaction phase is also called the self-trapping phase. The analogue of the latter in the quantum BJJ is Hilbert-space localization, reflected by the scaling behavior of the inverse participation ratio [16]. This also indicates the existence of a critical uc. 2 2 In Fig. 2 the effect of time averaging on |ur,k| and |ut,k| for u = 1 is visualized. It reflects the smoothing of the strongly fluctuating dynamics with a recurrent and a decaying contribution according to Eqs. (7), (9). It is obvious that the separation of the recurrent and the decaying behavior is not accessible without time averaging. The existence of a critical interaction strength uc ≈ 1.89... is demonstrated in Fig. 3, where the 7

0.9 average return prob. return prob. 0.8

0.7

0.6

0.5

0.4

return probability 0.3

0.2

0.1

0 0 10 20 30 40 50 60 70 k 4.5×10−7 average transition prob. transition prob. 4×10−7

3.5×10−7

3×10−7

2.5×10−7

2×10−7

1.5×10−7 transition probability

1×10−7

5×10−8

0 0 10 20 30 40 50 60 70 k

2 2 FIG. 2. Comparison of the return probability |ur,k| and the average return probability h|ur,k| iτ (top panel) and of the corresponding transition probabilities (bottom panel). The average was taken with respect to the exponential distribution of Eq. (16). The interaction parameter is u = 1 andτ ¯ = 1/10.

2 2 h|ur,k| iτ jumps up upon increasing u at uc. Moreover, h|ut,k| iτ develops a characteristic peak at uc. This behavior reflects the appearance of nearly degenerate energy levels, as described in Sect. II B. 2 Finally, in Fig. 4 the change of the time scales for the decay of the average TP h|ut,k| iτ is visualized for u = 1.7, ..., 2.2. The decay is reduced by an increasing interaction strength u. This reflects the fact that the splitting of the energy levels is reduced by the interaction, as we would have expected.

IV. DISCUSSION AND CONCLUSIONS

Our analysis of the quantum unitary evolution is strictly focused on the result of a single measurement in each of many identical experiments, which are subject of a unitary evolution. Averaging with respect to the statistical outcome due to measurements at randomly distributed times leads to the TCM. We have focused on the diagonal TCM elements to study the evolution of the quantum system. The analysis of the off-diagonal TCM elements was the subject of a previous work [15]. Similar to the off-diagonal TCM elements, the diagonal TCM elements reveal a separation of the evolution into a static recurrent part and a dynamic decaying part. We have found that the decay rate of the latter is related to the spacing between energy levels, which diverges when the spacing vanishes near a degeneracy. Thus, the decay rate is a quantity that can be used to detect symmetry changes or the appearance of spontaneous symmetry breaking. This has been observed in the example of the BJJ: In this model the energy levels are two-fold degenerate in the limit u → ∞. This is reflected in Fig. 4, where the decay decreases with 8

0.1 average transition prob. 0.09 average return prob.

0.08

0.07

0.06

0.05

0.04

0.03

0.02 average transition/return probability 0.01

0 1.82 1.84 1.86 1.88 1.9 1.92 1.94 1.96 u

2 FIG. 3. The critical regime of the Hilbert-space localization with uc ≈ 1.89 is visualized with h|ur,k| iτ and 2 h|ut,k| iτ at k = 70.

0.0014 u=1.7 u=1.8 0.0012 u=2.0 u=2.2

0.001

0.0008

0.0006

0.0004 average transition probability 0.0002

0 15 20 25 30 35 40 45 50 k

2 FIG. 4. Decay of h|ut,k| iτ for different interaction parameters u = 1.7,..., 2.2.

increasing u. Another interesting aspect of the BJJ is the transition to Hilbert-space localization [16]. This transition was also detected with the helps of the average RP and average TP in Fig. 3, where the average RP experience a jump to a higher value for u>uc. On the other hand, the average TP has only a sharp peak near uc but has the same value away from uc. We can conclude that time averaging over an ensemble of measurements is crucial for extracting the (static) recurrent behavior and the (dynamic) decaying behavior. This can be formulated in terms of the TCM. The separation of static and dynamic behavior is generic for the unitary quantum evolution. Then the TCM provides a tool to analyze properties of the quantum system, which are associated with spectral degeneracy in the case of phase transitions. It can be applied to theoretical calculations as well as to experimental data, collected from many experiments at identical quantum systems. We have demonstrated in the case of the BJJ that the TCM approach delivers interesting generic information. The BJJ can be considered as a building block of the Bose-Hubbard model on a lattice. Therefore, the TCM approach should be applicable to more complex quantum models, including bosonic and fermionic Hubbard models as well as quantum spin systems. 9

ACKNOWLEDGMENTS

I am grateful to Eli Barkai for interesting discussions regarding the Hilbert-space fragmentation. This work was supported by the Julian Schwinger Foundation.

[1] E. P. Wigner, Characteristic vectors of bordered matrices with infinite dimensions, Annals of Mathematics 62, 548 (1955). [2] C. E. Porter and R. G. Thomas, Fluctuations of nuclear reaction widths, Phys. Rev. 104, 483 (1956). [3] F. J. Dyson, Statistical theory of the energy levels of complex systems. i, Journal of Mathematical Physics 3, 140 (1962), https://doi.org/10.1063/1.1703773. [4] C. E. Porter, Book Review: Statistical theories of spectra: fluctuations. C.E. PORTER (Academic Press, New York, 1965. xv-576 p. 5.95 paper, 9.50 cloth), Nuclear Physics 78, 696 (1966). [5] E. P. Wigner, Random matrices in physics, SIAM Review 9, 1 (1967). [6] M. Mehta, Random Matrices, Pure and applied mathematics No. v. 142 (Elsevier/Academic Press, 2004). [7] C. W. J. Beenakker, Random-matrix theory of quantum transport, Rev. Mod. Phys. 69, 731 (1997). [8] F. J. Dyson, Statistical theory of the energy levels of complex systems. ii, Journal of Mathematical Physics 3, 157 (1962), https://doi.org/10.1063/1.1703774. [9] F. J. Dyson, Statistical theory of the energy levels of complex systems. iii, Journal of Mathematical Physics 3, 166 (1962), https://doi.org/10.1063/1.1703775. [10] S. J. van Enk and C. W. J. Beenakker, Measuring Trρn on single copies of ρ using random measurements, Phys. Rev. Lett. 108, 110503 (2012). [11] A. Elben, B. Vermersch, M. Dalmonte, J. I. Cirac, and P. Zoller, R´enyi entropies from random quenches in atomic hubbard and spin models, Phys. Rev. Lett. 120, 050406 (2018). [12] B. Vermersch, A. Elben, M. Dalmonte, J. I. Cirac, and P. Zoller, Unitary n-designs via random quenches in atomic hubbard and spin models: Application to the measurement of r´enyi entropies, Phys. Rev. A 97, 023604 (2018). [13] Y. Li, X. Chen, and M. P. A. Fisher, Measurement-driven entanglement transition in hybrid quantum circuits, Phys. Rev. B 100, 134306 (2019). [14] B. Skinner, J. Ruhman, and A. Nahum, Measurement-induced phase transitions in the dynamics of entan- glement, Phys. Rev. X 9, 031009 (2019). [15] K. Ziegler, Quantized dynamics in closed quantum systems, Journal of Physics A: Mathematical and Theo- retical 54, 205303 (2021). [16] D. Cohen, V. I. Yukalov, and K. Ziegler, Hilbert-space localization in closed quantum systems, Phys. Rev. A 93, 042101 (2016). [17] S. Moudgalya and O. I. Motrunich, Hilbert space fragmentation and commutant algebras (2021), arXiv:2108.10324 [cond-mat.stat-mech]. [18] P. Kos, M. Ljubotina, and T. c. v. Prosen, Many-body quantum chaos: Analytic connection to random matrix theory, Phys. Rev. X 8, 021062 (2018). [19] A. Chan, A. De Luca, and J. T. Chalker, Solution of a minimal model for many-body quantum chaos, Phys. Rev. X 8, 041019 (2018). [20] P. W. Anderson, Absence of diffusion in certain random lattices, Phys. Rev. 109, 1492 (1958). [21] N. Bar-Gill, G. Kurizki, M. Oberthaler, and N. Davidson, Dynamic control and probing of many-body decoherence in double-well bose-einstein condensates, Phys. Rev. A 80, 053613 (2009). [22] B. Juli´a-D´ıaz, T. Zibold, M. K. Oberthaler, M. Mel´e-Messeguer, J. Martorell, and A. Polls, Dynamic gener- ation of spin-squeezed states in bosonic josephson junctions, Phys. Rev. A 86, 023615 (2012). [23] R. Barends, J. Kelly, A. Megrant, D. Sank, E. Jeffrey, Y. Chen, Y. Yin, B. Chiaro, J. Mutus, C. Neill, P. O’Malley, P. Roushan, J. Wenner, T. C. White, A. N. Cleland, and J. M. Martinis, Coherent josephson qubit suitable for scalable quantum integrated circuits, Phys. Rev. Lett. 111, 080502 (2013). [24] R. Gati and M. K. Oberthaler, A bosonic josephson junction, Journal of Physics B: Atomic, Molecular and Optical Physics 40, R61 (2007). [25] K. Ziegler, Dynamical creation of entangled bosonic states in a double well, Journal of Physics B: Atomic, Molecular and Optical Physics 44, 145302 (2011). [26] G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, Quantum dynamics of an atomic bose-einstein condensate in a double-well potential, Phys. Rev. A 55, 4318 (1997).