Identification of quantum scars via phase-space localization measures Sa´ulPilatowsky-Cameo,1 David Villase~nor,1 Miguel A. Bastarrachea-Magnani,2 Sergio Lerma-Hern´andez,3 Lea F. Santos,4 and Jorge G. Hirsch1 1Instituto de Ciencias Nucleares, Universidad Nacional Aut´onomade M´exico, Apdo. Postal 70-543, C.P. 04510 CDMX, Mexico 2Departamento de F´ısica, Universidad Aut´onomaMetropolitana-Iztapalapa, San Rafael Atlixco 186, C.P. 09340 CDMX, Mexico 3Facultad de F´ısica, Universidad Veracruzana, Circuito Aguirre Beltr´ans/n, C.P. 91000 Xalapa, Veracruz, Mexico 4Department of Physics, Yeshiva University, New York, New York 10016, USA There is no unique way to quantify the degree of delocalization of quantum states in unbounded continuous spaces. In this work, we explore a recently introduced localization measure that quantifies the portion of the classical phase space occupied by a quantum state. The measure is based on the α-moments of the Husimi function and is known as the R´enyi occupation of order α. With this quantity and random pure states, we find a general expression to identify states that are maximally delocalized in phase space. Using this expression and the Dicke model, which is an interacting spin- boson model with an unbounded four-dimensional phase space, we show that the R´enyi occupations with large α are highly effective at revealing quantum scars. Furthermore, by analyzing the large moments of the Husimi function, we are able to find the unstable periodic orbits that scar some of the eigenstates of the model. I. INTRODUCTION in the chaotic regime and distinguish the eigenstates that are maximally delocalized from those that are highly lo- The notion of localization in quantum systems presup- calized in phase space. The Dicke model has a rich and poses a basis representation. Anderson localization [1, 2], unbounded phase space. This spin-boson model [34] was for example, refers to the suppression of the classical dif- initially introduced to explain the phenomenon of super- fusion of particles in real space due to quantum inter- radiance [35{38] and has since fostered a wide variety of ferences. This phenomenon has a dynamical counterpart theoretical studies, including the behavior of out-of-time- originally studied in the context of the kicked rotor [3{6]. ordered correlators [39{41], manifestations of quantum Localization is an extremely broad subject that extends scarring [42{46], non-equilibrium dynamics [38, 47{51], to the hydrogen atom in a monochromatic field [7, 8], and measures of quantum localization with respect to Rydberg atoms [9], quantum billiards [10{16], banded phase space [45, 52]. The model is also of great inter- random matrices [17], interacting systems [18{22], and est to experiments with trapped ions [53, 54], supercon- kicked interacting systems [23{26], among others. ducting circuits [55], and cavity assisted Raman transi- In finite spaces, the degree of delocalization of a quan- tions [56, 57]. tum state is based on how much it spreads on a chosen By comparing the R´enyi occupations of the high- finite basis representation, as quantified by participation energy eigenstates of the Dicke model with the R´enyi oc- ratios [5, 27] and R´enyi entropies [28, 29]. In unbounded cupations of random states, as a function of α, we range continuous spaces, such measures can reach arbitrarily the eigenstates according to their degree of delocaliza- large values. If one wishes to investigate maximally de- tion. While all eigenstates are scarred [45], the most localized states in these systems, one must define a finite localized states are the ones scarred by unstable periodic region as a reference volume. There are multiple ways orbits with short periods. The analysis of the fourth- to select this bounded region, and each choice produces moment of the Husimi functions of these highly localized a different localization measure [30]. A natural volume states allow us to identify the classical unstable periodic orbits that cause their scarring. These orbits are different arXiv:2107.06894v1 [quant-ph] 14 Jul 2021 of reference is obtained by measuring the degree of lo- calization of a quantum state over the classical phase- from the ones found in [46] and their uncovering repre- space energy shell, as proposed in [30]. This measure of sents a step forward in the difficult task of identifying localization has been named R´enyi occupation, because the unstable periodic orbits underlying quantum scars. it is related to the exponential of the R´enyi-Wehrl en- In addition, we study the dynamical properties of the tropies [31, 32]. Quantum states can be represented in unveiled orbits and use them to explain the structures of the phase space through quasi-probability distributions, the eigenstates and their degree of localization. such as the Husimi function [33]. The R´enyi occupations The paper is organized as follows. The Dicke model of order α refer to averages of normalized α-moments of and its phase space are presented in Sec. II. In Sec. III, the Husimi functions. we introduce a measure of localization of quantum states In this work, we use maximally delocalized random based on the R´enyi occupations of order α and provide pure states and obtain a general analytical expression an analytical expression for maximally delocalized states. for the R´enyi occupations of order α. We employ this In Sec. IV we discuss quantum scarring in the light of the expression to analyze the eigenstates of the Dicke model R´enyi occupations and identify the unstable periodic or- 2 bits that scar the most localized eigenstates. We also system size, interpret our localization measure in terms of the geo- hxjH^Djxi metric and dynamical properties of the classical periodic hcl(x) = = (2) orbits. Our conclusions are presented in Sec. V. j r ! ! P 2 + Q2 p2 + q2 + 0 P 2 + Q2 + 2γqQ 1 − − ! : 2 2 4 0 II. DICKE MODEL The bosonic Glauber and atomic Bloch coherent states are, respectively, hp i 2 2 y The Dicke Hamiltonian is a fundamental model in jq; pi = e−(j=4)(q +p )e j=2(q+ip) a^ j0i; (3) quantum optics that describes a single mode of the elec- 2 2 h p i P + Q j 2 2 ^ tromagnetic field interacting with a set of N two-level jQ; P i = 1 − e (Q+iP )= 4−P −Q J+ jj; −ji; atoms [34]. Setting ~ = 1, the Hamiltonian is given by 4 with j0i being the photon vacuum and jj; −ji the state y γ y with all the atoms in their ground state. H^D = !a^ a^ + !0J^z + p (^a +a ^)(J^+ + J^−); (1) N B. Phase space wherea ^y (^a) is the bosonic creation (annihilation) oper- ator of the field mode and J^+ (J^−) is the raising (lower- The classical Hamiltonian hcl(x) has a four- ing) operator defined by J^± = J^x ± iJ^y, where J^x;y;z = N dimensional phase space M in the coordinates x = P k 2 2 (1=2) k=1 σ^x;y;z are the collective pseudo-spin opera- (q; p; Q; P ), where P + Q ≤ 4, because the pseudo- tors satisfying the SU(2) algebra andσ ^x;y;z are the Pauli spin degree of freedom is bounded. This space can be matrices. The Hamiltonian parameters are the radiation partitioned into a family of classical energy shells frequency of the electromagnetic field !, the transition M() = fx 2 M j h (x) = g frequency of the two-level atoms !0, and the atom-field cl coupling strength γ. Since the Hamiltonian commutes in terms of the rescaled classical energy = E=j. Due ^2 ^2 ^2 ^2 with the total pseudo-spin operator J = Jx + Jy + Jz , to energy conservation, for a given initial condition x 2 the Hilbert space is separated into different invariant sub- M(), the classical evolution given by hcl remains in 2 spaces for each eigenvalue j(j +1) of J^ . We work within M(), that is, x(t) 2 M() for all times t. the totally symmetric subspace that includes the ground- The classical energy shells are bounded with respect to state of the system and is defined by the maximum value the three-dimensional surface measure dx δ(hcl(x) − ). of the pseudo-spin length j = N =2. The model also The finite volume of M() is given by y ^ ^ iπ(^a a^+Jz +j) Z possesses a discrete parity symmetry, Π = e , 2 which leads to an additional subspace separation accord- dx δ(hcl(x) − ) = (2π~eff) ν(); (4) M ing to the eigenvalues of the parity operator Π± = ±1. where ~eff = 1=j is the effective Planck constant [65], The model has been used in studies of the quantum 2 which determines the volume of the Plank cell (2π~eff) , phase transition from a normal (γ < γc) to a superradiant and ν() is the semiclassical density of states obtained by (γ > γ ) phase, which arises when the coupling strength c p taking only the first term in the Gutzwiller trace formula reaches the critical value γc = !!0=2 [35, 36, 58, 59]. [62, 66, 67] (see App. A for more details on this quantity). The model also exhibits regular or chaotic behavior de- To calculate the average value, hfi, of a phase-space pending on the Hamiltonian parameters and excitation function f(x) over the classical energy shell M(), we in- energies [60]. In this paper, we choose the resonant fre- tegrate with respect to the surface measure dx δ(hcl(x)− quency case ! = !0 = 1, system size j = N =2 = 100, ) and divide the result by the total volume (2π )2ν(), and γ = 2γ to work in the superradiant phase.