Chaos and quantum scars in coupled top model Debabrata Mondal, Sudip Sinha and S. Sinha Indian Institute of Science Education and Research-Kolkata, Mohanpur, Nadia-741246, India (Dated: September 2, 2020) We consider a coupled top model describing two interacting large spins, which is studied semiclas- sically as well as quantum mechanically. This model exhibits variety of interesting phenomena such as quantum phase transition (QPT), dynamical transition and excited state quantum phase transi- tions above a critical coupling strength. Both classical dynamics and entanglement entropy reveals ergodic behavior at the center of energy density band for an intermediate range of coupling strength above QPT, where the level spacing distribution changes from Poissonian to Wigner-Dyson statis- tics. Interestingly, in this model we identify quantum scars as reminiscence of unstable collective dynamics even in presence of interaction. Statistical properties of such scarred states deviate from ergodic limit corresponding to random matrix theory and violate Berry’s conjecture. In contrast to ergodic evolution, oscillatory behavior in dynamics of unequal time commutator and survival prob- ability is observed as dynamical signature of quantum scar, which can be relevant for its detection. PACS numbers: Introduction: After the recent experiment on a chain of field Ising model [40, 41], which is governed by the Hamil- Rydberg atoms [1], quantum scar (QS) in many body sys- tonian, tems has drawn significant interest due to its connection µ with ergodicity and non equilibrium dynamics. Observed ˆ = Sˆ1x Sˆ2x Sˆ1zSˆ2z, (1) athermal behavior and periodic revival phenomena in dy- H − − − S namics of a specific initial state has been attributed to the where Sˆia represents components (a = x, y, z) of two many body quantum scarring phenomena [2–8]. Many large spins of equal magnitude S denoted by index i = body quantum scar (MBQS) gives rise to the deviation 1, 2 and µ is the ferromagnetic coupling strength. from ergodicity [2–22] leading to the violation of eigen- Large magnitude of spin S 1 allows us to analyze the state thermalization hypothesis (ETH) [23] in quantum model semiclassically, where the spin vectors are repre- systems [5, 8–12]. Recent theoretical studies on interact- sented by S~i = (S sin θi cos φi,S sin θi sin φi,S cos θi). In ing spin systems reveal that emergent symmetries and terms of these dynamical variables, corresponding classi- symmetry protected many body states are the key ingre- cal Hamiltonian can be written as, dient for MBQS [11, 13–21]. Originally QS in non in- teracting quantum chaotic system has been identified as = 1 z 2 cos φ 1 z 2 cos φ µ z z (2) reminiscence of unstable classical orbits [24, 25]. However Hcl − − 1 1 − − 2 2 − 1 2 p p such connection between MBQS and underlying dynam- where zi = cos θi is conjugate momentum corresponding ics remains unclear for interacting many body systems. to the variable φi. Here cl and classical energy E are Similar problem also arises in understanding ergodicity scaled by S. The classicalH spin dynamics is described by in quantum many particle systems from the viewpoint of following equations of motion (EOM), chaos and phase space mixing. However some studies on quantum collective models reveal a connection between zi cos φi z˙ = 1 z 2 sin φ ; φ˙ = µz¯ (3) i i i i 2 i ergodicity and underlying chaotic dynamics [26–29], as − − √1 zi − well elucidate the role of underlying unstable collective p − dynamics in the formation of MBQS [30]. where ¯i = i. To understand the overall behavior of dy- namics, we6 first analyze the fixed points (FP) and their In this work we consider coupled top model which al- stability for varying coupling strength µ, as charted in lows us to study the collective dynamics of two large Fig.1(a). Due to the interaction between the spins, at spins semiclassically, revealing interesting phenomena a critical coupling strength µc = 1, the model under- like quantum phase transition (QPT), dynamical tran- goes a QPT to ferromagnetically ordered ground state sition (DT) [31–34], and chaotic behavior at an interme- arXiv:2005.07048v2 [cond-mat.quant-gas] 1 Sep 2020 (GS) as well as dynamical transition corresponding to diate regime of coupling strength. More importantly, in the highest excited state (ES) with anti-ferromagnetic the chaotic regime, we identify QS as a signature of un- ordering. For µ < µ , two symmetry unbroken sta- stable orbits of collective spin dynamics within certain c ble FPs are represented by: (I) z1 = 0, φ1 = 0, z2 = symmetry classes. Statistical and dynamical properties 0, φ = 0 with energy density E{ = 2 (GS) and (II) of such scarred states exhibit clear deviation from ergod- 2 } − z1 = 0, φ1 = π, z2 = 0, φ2 = π with E = 2 (ES). For icity. symmetry{ unbroken phase FP-I} (FP-II), both spins are Model and semiclassical analysis: The coupled top aligned to positive (negative) x-axis without any mag- model [35–39] describes the dynamics of ferromagneti- netization along z-axis. Both FP-I and II undergoes cally coupled two large spins similar to the transverse a pitchfork bifurcation at µc and become unstable for 2 µ > µc. As a result of bifurcation, two pairs of symme- shown in Fig.1(c). Whereas above the critical point, the try broken stable steady states appear above µc, namely: symmetry broken FPs with z = 0 in the phase portrait 2 6 (III) z1 = z2 = 1 1/µ , φ1 = φ2 = 0 with E = correspond to QPT and dynamical transition, however (µ +{ 1/µ) which± corresponds− to the ferromagnetic} GS intermediate region of phase space is filled with irreg- − p 2 and (IV) z1 = z2 = 1 1/µ , φ1 = φ2 = π rep- ular (chaotic) trajectories (see Fig.1(d)). Such mixed resenting{ a dynamically− ± stable− anti-ferromagnetic} state phase space structure indicates interesting ergodic be- with E = (µ+1/µ) correspondingp to the ES. Apart from havior [4, 43, 44], which leads us to further study the these states, there exists another pair of unstable FPs, model quantum mechanically. which we denote as ‘π-mode’, describing the relative an- Quantum chaos and ergodicity: To analyze the system gle between the spins: (V) z = z = 0, φ = π, φ = 0 quantum mechanically, we diagonalize the Hamiltonian { 1 2 1 2 } and z1 = z2 = 0, φ1 = 0, φ2 = π with E = 0. Eq.(1) for a fixed magnitude of spin S so that the dimen- { } sionality of Hilbert space becomes = (2S + 1)2. The N (a) (c)1.5 (e) eigenvalues n and corresponding eigenfunctions ψn of the HamiltonianE are obtained using the basis states| i FP-II FP-IV (AFM) 1.0 ˆ DT FP-II m1z, m2z , where miz is eigenvalue of Siz. The eigen- FP-V | i 0.5 functions and eigenvalues can be classified according to QPT FP-I the symmetry of the Hamiltonian (Eq.(1)) which remains 0.0 FP-I FP-III (FM) ıπ ˆ invariant under parity Πˆ = e− P where ˆ = Sˆ1x + Sˆ2x -0.5 and under exchange symmetry between twoP spins (S -1.0 -0.5 0.0 0.5 1.0 1 ↔ S2) [45]. From the quantum mechanical analysis of the (b) (d)1.5 (f) model, we notice for µ > µc, both QPT and dynami- 1.0 cal transition lead to excited state quantum phase tran- 0.5 sitions (ESQPT) [46–50] corresponding to the unstable symmetry unbroken FPs I and II with energy densities 0.0 E = 2 and 2, respectively, where the derivative of semi- − -0.5 classical density of state becomes singular [42, 46–50]. -1.0 -0.5 0.0 0.5 1.0 Next, we study the entanglement entropy (EE) of the eigenstates which contains useful information re- FIG. 1: (a) Schematic diagram of different steady states at lated to ergodicity, as it is expected that chaoticity different energies with increasing µ. Solid (dotted) lines repre- and phase space mixing can lead to enhancement of EE sent stable (unstable) branches and vertical pink dashed line [42, 43, 51, 52]. The reduced density matrix of -th spin S marks critical coupling µc. (b) Small oscillation frequency as ρˆ = Tr ¯ ψn ψn is obtained by tracing out other spin S S | i h | a function of µ for QPT (black solid line) and DT (red circles). sector ( ¯ = ) which yields the EE, Sen = Tr(ˆρ logˆρ ) S 6 S − S S Phase portraits ((c) and (d)) and trajectories on Bloch sphere of the corresponding state ψn . We quantify the de- ((e) and (f)) exhibiting QPT and DT. For (c),(e) µ = 0.8 < µc gree of ergodicity of a state by| comparingi it with the EE and (d),(f) µ = 1.3 > µc. of maximally random state partitioned into subsystems A(B) with dimension A( B). Maximal EE correspond- ing to subsystem A withD DD D [53] is given by, Interestingly, the dynamical transition corresponding A ≤ B to ES gives rise to the anti-ferromagnetic spin configura- S log( ) /2 (4) tion, which survives as a dynamically stable steady state max ' DA − DA DB (FP-IV) in spite of ferromagnetic interaction between where A = B = 2S + 1. For any value of coupling them. This can be understood from the transformation strengthDµ, theD eigenstates with maximum EE are located ˆ ˆ Sx Sx and µ µ in Eq.(1), which converts the ES at the band center with energy density E 0, as shown in → − → − to GS of anti-ferromagnetic coupled top model. As shown Fig.2(b). Moreover, EE corresponding to≈ the band center in Fig.1(b), small amplitude oscillation frequencies (ex- (E 0), attains its maximum possible limit Smax (given citation frequencies) corresponding to both GS and ES in Eq.(4))≈ around µ 2, indicating ergodic behavior.