Statistical properties of the localization measure of chaotic eigenstates in the Dicke model

Chaos and quantization of the three-particle generic Fermi-Pasta-Ulam-Tsingou model. II. Phenomenology of quantum eigenstates

We undertake a thorough investigation into the phenomenology of quantum eigenstates, in the three-particle Fermi-Pasta-Ulam-Tsingou model. Employing different Husimi functions, our study focuses on both the α-type, which is canonically equivalent to the celebrated Hénon-Heiles Hamiltonian, a nonintegrable and mixed-type system, and the general case at the saddle energy where the system is fully chaotic. Based on Husimi quantum surface of sections, we find that in the mixed-type system, the fraction of mixed eigenstates in an energy shell [E-δE/2,E+δE/2] with δE≪E shows a power-law decay with respect to the decreasing Planck constant ℏ. Defining the localization measures in terms of the Rényi-Wehrl entropy, in both the mixed-type and fully chaotic systems, we find a better fit with the β distribution and a lesser degree of localization, in the distribution of localization measures of chaotic eigenstates, as the controlling ratio α_{L}=t_{H}/t_{T} between the Heisenberg time t_{H} and the classical transport time t_{T} increases. This transition with respect to α_{L} and the power-law decay of the mixed states, together provide supporting evidence for the principle of uniform semiclassical condensation in the semiclassical limit. Moreover, we find that in the general case which is fully chaotic, the maximally localized state, is influenced by the stable and unstable manifold of the saddles (hyperbolic fixed points), while the maximally extended state notably avoids these points, extending across the remaining space, complementing each other.

Mixed eigenstates in the Dicke model: Statistics and power-law decay of the relative proportion in the semiclassical limit

How the mixed eigenstates vary when approaching the semiclassical limit in mixed-type many-body quantum systems is an interesting but still less known question. Here, we address this question in the Dicke model, a celebrated many-body model that has a well defined semiclassical limit and undergoes a transition to chaos in both quantum and classical cases. Using the Husimi function, we show that the eigenstates of the Dicke model with mixed-type classical phase space can be classified into different types. To quantitatively characterize the types of eigenstates, we study the phase space overlap index, which is defined in terms of the Husimi function. We look at the probability distribution of the phase space overlap index and investigate how it changes with increasing system size, that is, when approaching the semiclassical limit. We show that increasing the system size gives rise to a power-law decay in the behavior of the relative proportion of mixed eigenstates. Our findings shed more light on the properties of eigenstates in mixed-type many-body systems and suggest that the principle of uniform semiclassical condensation of Husimi functions should also be valid for many-body quantum systems.

Analysis of chaos and regularity in the open Dicke model

We present an analysis of chaos and regularity in the open Dicke model, when dissipation is due to cavity losses. Due to the infinite Liouville space of this model, we also introduce a criterion to numerically find a complex spectrum which approximately represents the system spectrum. The isolated Dicke model has a well-defined classical limit with two degrees of freedom. We select two case studies where the classical isolated system shows regularity and where chaos appears. To characterize the open system as regular or chaotic, we study regions of the complex spectrum taking windows over the absolute value of its eigenvalues. Our results for this infinite-dimensional system agree with the Grobe-Haake-Sommers (GHS) conjecture for Markovian dissipative open quantum systems, finding the expected 2D Poisson distribution for regular regimes, and the distribution of the Ginibre unitary ensemble (GinUE) for the chaotic ones, respectively.

Power-law decay of the fraction of the mixed eigenstates in kicked top model with mixed-type classical phase space

The properties of mixed eigenstates in a generic quantum system with a classical counterpart that has mixed-type phase space, although important to understand several fundamental questions that arise in both theoretical and experimental studies, are still not clear. Here, following a recent work [Č. Lozej, D. Lukman, and M. Robnik, Phys. Rev. E 106, 054203 (2022)2470-004510.1103/PhysRevE.106.054203], we perform an analysis of the features of mixed eigenstates in a time-dependent Hamiltonian system, the celebrated kicked top model. As a paradigmatic model for studying quantum chaos, the kicked top model is known to exhibit both classical and quantum chaos. The types of eigenstates are identified by means of the phase-space overlap index, which is defined as the overlap of the Husimi function with regular and chaotic regions in classical phase space. We show that the mixed eigenstates appear due to various tunneling precesses between different phase-space structures, while the regular and chaotic eigenstates are, respectively, associated with invariant tori and chaotic components in phase space. We examine how the probability distribution of the phase-space overlap index evolves with increasing system size for different kicking strengths. In particular, we find that the relative fraction of mixed states exhibits a power-law decay as the system size increases, indicating that only purely regular and chaotic eigenstates are left in the strict semiclassical limit. We thus provide further verification of the principle of uniform semiclassical condensation of Husimi functions and confirm the correctness of the Berry-Robnik picture.

Statistics of phase space localization measures and quantum chaos in the kicked top model

Quantum chaos plays a significant role in understanding several important questions of recent theoretical and experimental studies. Here, by focusing on the localization properties of eigenstates in phase space (by means of Husimi functions), we explore the characterizations of quantum chaos using the statistics of the localization measures, that is the inverse participation ratio and the Wehrl entropy. We consider the paradigmatic kicked top model, which shows a transition to chaos with increasing the kicking strength. We demonstrate that the distributions of the localization measures exhibit a drastic change as the system undergoes the crossover from integrability to chaos. We also show how to identify the signatures of quantum chaos from the central moments of the distributions of localization measures. Moreover, we find that the localization measures in the fully chaotic regime apparently universally exhibit the beta distribution, in agreement with previous studies in the billiard systems and the Dicke model. Our results contribute to a further understanding of quantum chaos and shed light on the usefulness of the statistics of phase space localization measures in diagnosing the presence of quantum chaos, as well as the localization properties of eigenstates in quantum chaotic systems.

Spectral Form Factor and Dynamical Localization

Quantum dynamical localization occurs when quantum interference stops the diffusion of wave packets in momentum space. The expectation is that dynamical localization will occur when the typical transport time of the momentum diffusion is greater than the Heisenberg time. The transport time is typically computed from the corresponding classical dynamics. In this paper, we present an alternative approach based purely on the study of spectral fluctuations of the quantum system. The information about the transport times is encoded in the spectral form factor, which is the Fourier transform of the two-point spectral autocorrelation function. We compute large samples of the energy spectra (of the order of 106 levels) and spectral form factors of 22 stadium billiards with parameter values across the transition between the localized and extended eigenstate regimes. The transport time is obtained from the point when the spectral form factor transitions from the non-universal to the universal regime predicted by random matrix theory. We study the dependence of the transport time on the parameter value and show the level repulsion exponents, which are known to be a good measure of dynamical localization, depend linearly on the transport times obtained in this way.

A Physical Measure for Characterizing Crossover from Integrable to Chaotic Quantum Systems

In this paper, a quantity that describes a response of a system's eigenstates to a very small perturbation of physical relevance is studied as a measure for characterizing crossover from integrable to chaotic quantum systems. It is computed from the distribution of very small, rescaled components of perturbed eigenfunctions on the unperturbed basis. Physically, it gives a relative measure to prohibition of level transitions induced by the perturbation. Making use of this measure, numerical simulations in the so-called Lipkin-Meshkov-Glick model show in a clear way that the whole integrability-chaos transition region is divided into three subregions: a nearly integrable regime, a nearly chaotic regime, and a crossover regime.

Chaos in the three-site Bose-Hubbard model: Classical versus quantum

We consider a quantum many-body system-the Bose-Hubbard system on three sites-which has a classical limit, and which is neither strongly chaotic nor integrable but rather shows a mixture of the two types of behavior. We compare quantum measures of chaos (eigenvalue statistics and eigenvector structure) in the quantum system, with classical measures of chaos (Lyapunov exponents) in the corresponding classical system. As a function of energy and interaction strength, we demonstrate a strong overall correspondence between the two cases. In contrast to both strongly chaotic and integrable systems, the largest Lyapunov exponent is shown to be a multivalued function of energy.

Chaos and Thermalization in the Spin-Boson Dicke Model

We present a detailed analysis of the connection between chaos and the onset of thermalization in the spin-boson Dicke model. This system has a well-defined classical limit with two degrees of freedom, and it presents both regular and chaotic regions. Our studies of the eigenstate expectation values and the distributions of the off-diagonal elements of the number of photons and the number of excited atoms validate the diagonal and off-diagonal eigenstate thermalization hypothesis (ETH) in the chaotic region, thus ensuring thermalization. The validity of the ETH reflects the chaotic structure of the eigenstates, which we corroborate using the von Neumann entanglement entropy and the Shannon entropy. Our results for the Shannon entropy also make evident the advantages of the so-called "efficient basis" over the widespread employed Fock basis when investigating the unbounded spectrum of the Dicke model. The efficient basis gives us access to a larger number of converged states than what can be reached with the Fock basis.

Canonical Density Matrices from Eigenstates of Mixed Systems

One key issue of the foundation of statistical mechanics is the emergence of equilibrium ensembles in isolated and closed quantum systems. Recently, it was predicted that in the thermodynamic (N→∞) limit of large quantum many-body systems, canonical density matrices emerge for small subsystems from almost all pure states. This notion of canonical typicality is assumed to originate from the entanglement between subsystem and environment and the resulting intrinsic quantum complexity of the many-body state. For individual eigenstates, it has been shown that local observables show thermal properties provided the eigenstate thermalization hypothesis holds, which requires the system to be quantum-chaotic. In the present paper, we study the emergence of thermal states in the regime of a quantum analog of a mixed phase space. Specifically, we study the emergence of the canonical density matrix of an impurity upon reduction from isolated energy eigenstates of a large but finite quantum system the impurity is embedded in. Our system can be tuned by means of a single parameter from quantum integrability to quantum chaos and corresponds in between to a system with mixed quantum phase space. We show that the probability for finding a canonical density matrix when reducing the ensemble of energy eigenstates of the finite many-body system can be quantitatively controlled and tuned by the degree of quantum chaos present. For the transition from quantum integrability to quantum chaos, we find a continuous and universal (i.e., size-independent) relation between the fraction of canonical eigenstates and the degree of chaoticity as measured by the Brody parameter or the Shannon entropy.

Phenomenology of quantum eigenstates in mixed-type systems: Lemon billiards with complex phase space structure

The boundary of the lemon billiards is defined by the intersection of two circles of equal unit radius with the distance 2B between their centers, as introduced by Heller and Tomsovic [E. J. Heller and S. Tomsovic, Phys. Today 46, 38 (1993)0031-922810.1063/1.881358]. This paper is a continuation of our recent papers on a classical and quantum ergodic lemon billiard (B=0.5) with strong stickiness effects [Č. Lozej et al., Phys. Rev. E 103, 012204 (2021)2470-004510.1103/PhysRevE.103.012204], as well as on the three billiards with a simple mixed-type phase space and no stickiness [Č. Lozej et al., Nonlin. Phenom. Complex Syst. 24, 1 (2021)1817-245810.33581/1561-4085-2021-24-1-1-18]. Here we study two classical and quantum lemon billiards, for the cases B=0.1953,0.083, which are mixed-type billiards with a complex structure of phase space, without significant stickiness regions. A preliminary study of their spectra was published recently [ Č. Lozej, D. Lukman, and M. Robnik, Physics 3, 888 (2021)10.3390/physics3040055]. We calculate a very large number (10^{6}) of consecutive eigenstates and their Poincaré-Husimi (PH) functions, and analyze their localization properties by studying the entropy localization measure and the normalized inverse participation ratio. We introduce an overlap index, which measures the degree of the overlap of PH functions with classically regular and chaotic regions. We observe the existence of regular states associated with invariant tori and chaotic states associated with the classically chaotic regions, and also the mixed-type states. We show that in accordance with the Berry-Robnik picture and the principle of uniform semiclassical condensation of PH functions, the relative fraction of mixed-type states decreases as a power law with increasing energy, thus, in the strict semiclassical limit, leaving only purely regular and chaotic states. Our approach offers a general phenomenological overview of the structural and localization properties of PH functions in quantum mixed-type Hamiltonian systems.

Quantum Chaos in the Extended Dicke Model

We systematically study the chaotic signatures in a quantum many-body system consisting of an ensemble of interacting two-level atoms coupled to a single-mode bosonic field, the so-called extended Dicke model. The presence of the atom-atom interaction also leads us to explore how the atomic interaction affects the chaotic characters of the model. By analyzing the energy spectral statistics and the structure of eigenstates, we reveal the quantum signatures of chaos in the model and discuss the effect of the atomic interaction. We also investigate the dependence of the boundary of chaos extracted from both eigenvalue-based and eigenstate-based indicators on the atomic interaction. We show that the impact of the atomic interaction on the spectral statistics is stronger than on the structure of eigenstates. Qualitatively, the integrablity-to-chaos transition found in the Dicke model is amplified when the interatomic interaction in the extended Dicke model is switched on.

Effective dimensions of infinite-dimensional Hilbert spaces: A phase-space approach

By employing Husimi quasiprobability distributions, we show that a bounded portion of an unbounded phase space induces a finite effective dimension in an infinite-dimensional Hilbert space. We compare our general expressions with numerical results for the spin-boson Dicke model in the chaotic energy regime, restricting its unbounded four-dimensional phase space to a classically chaotic energy shell. This effective dimension can be employed to characterize quantum phenomena in infinite-dimensional systems, such as localization and scarring.

Multifractality in Quasienergy Space of Coherent States as a Signature of Quantum Chaos

We present the multifractal analysis of coherent states in kicked top model by expanding them in the basis of Floquet operator eigenstates. We demonstrate the manifestation of phase space structures in the multifractal properties of coherent states. In the classical limit, the classical dynamical map can be constructed, allowing us to explore the corresponding phase space portraits and to calculate the Lyapunov exponent. By tuning the kicking strength, the system undergoes a transition from regularity to chaos. We show that the variation of multifractal dimensions of coherent states with kicking strength is able to capture the structural changes of the phase space. The onset of chaos is clearly identified by the phase-space-averaged multifractal dimensions, which are well described by random matrix theory in a strongly chaotic regime. We further investigate the probability distribution of expansion coefficients, and show that the deviation between the numerical results and the prediction of random matrix theory behaves as a reliable detector of quantum chaos.

Quantum localization measures in phase space

Measuring the degree of localization of quantum states in phase space is essential for the description of the dynamics and equilibration of quantum systems, but this topic is far from being understood. There is no unique way to measure localization, and individual measures can reflect different aspects of the same quantum state. Here we present a general scheme to define localization in measure spaces, which is based on what we call Rényi occupations, from which any measure of localization can be derived. We apply this scheme to the four-dimensional unbounded phase space of the interacting spin-boson Dicke model. In particular, we make a detailed comparison of two localization measures based on the Husimi function in the regime where the model is chaotic, namely, one that projects the Husimi function over the finite phase space of the spin and another that uses the Husimi function defined over classical energy shells. We elucidate the origin of their differences, showing that in unbounded spaces the definition of maximal delocalization requires a bounded reference subspace, with different selections leading to contextual answers.

Eigenstate thermalization scaling in approaching the classical limit

According to the eigenstate thermalization hypothesis (ETH), the eigenstate-to-eigenstate fluctuations of expectation values of local observables should decrease with increasing system size. In approaching the thermodynamic limit-the number of sites and the particle number increasing at the same rate-the fluctuations should scale as ∼D^{-1/2} with the Hilbert space dimension D. Here, we study a different limit-the classical or semiclassical limit-by increasing the particle number in fixed lattice topologies. We focus on the paradigmatic Bose-Hubbard system, which is quantum-chaotic for large lattices and shows mixed behavior for small lattices. We derive expressions for the expected scaling, assuming ideal eigenstates having Gaussian-distributed random components. We show numerically that, for larger lattices, ETH scaling of physical midspectrum eigenstates follows the ideal (Gaussian) expectation, but for smaller lattices, the scaling occurs via a different exponent. We examine several plausible mechanisms for this anomalous scaling.

Characterizing the Lipkin-Meshkov-Glick model excited-state quantum phase transition using dynamical and statistical properties of the diagonal entropy

Using the diagonal entropy, we analyze the dynamical signatures of the Lipkin-Meshkov-Glick model excited-state quantum phase transition (ESQPT). We first show that the time evolution of the diagonal entropy behaves as an efficient indicator of the presence of an ESQPT. We also compute the probability distribution of the diagonal entropy values over a certain time interval and we find that the resulting distribution provides a clear distinction between the different phases of ESQPT. Moreover, we observe that the probability distribution of the diagonal entropy at the ESQPT critical point has a universal form, well described by a beta distribution, and that a reliable detection of the ESQPT can be obtained from the diagonal entropy central moments.

Effects of stickiness in the classical and quantum ergodic lemon billiard

We study the classical and quantum ergodic lemon billiard introduced by Heller and Tomsovic in Phys. Today 46(7), 38 (1993)PHTOAD0031-922810.1063/1.881358, for the case B=1/2, which is a classically ergodic system (without a rigorous proof) exhibiting strong stickiness regions around a zero-measure bouncing ball modes. The structure of the classical stickiness regions is uncovered in the S-plots introduced by Lozej [Phys. Rev. E 101, 052204 (2020)10.1103/PhysRevE.101.052204]. A unique classical transport or diffusion time cannot be defined. As a consequence the quantum states are characterized by the following nonuniversal properties: (i) All eigenstates are chaotic but localized as exhibited in the Poincaré-Husimi (PH) functions. (ii) The entropy localization measure A (also the normalized inverse participation ratio) has a nonuniversal distribution, typically bimodal, thus deviating from the beta distribution, the latter one being characteristic of uniformly chaotic systems with no stickiness regions. (iii) The energy-level spacing distribution is Berry-Robnik-Brody (BRB), capturing two effects: the quantally divided phase space (because most of the PH functions are either the inner-ones or the outer-ones, dictated by the classical stickiness, with an effective parameter μ_{1} measuring the size of the inner region bordered by the sticky invariant object, namely, a cantorus), and the localization of PH functions characterized by the level repulsion (Brody) parameter β. (iv) In the energy range considered (between 20 000 states to 400 000 states above the ground state) the picture (the structure of the eigenstates and the statistics of the energy spectra) is not changing qualitatively, as β fluctuates around 0.8, while μ_{1} decreases almost monotonically, with increasing energy.

Ubiquitous quantum scarring does not prevent ergodicity

In a classically chaotic system that is ergodic, any trajectory will be arbitrarily close to any point of the available phase space after a long time, filling it uniformly. Using Born’s rules to connect quantum states with probabilities, one might then expect that all quantum states in the chaotic regime should be uniformly distributed in phase space. This simplified picture was shaken by the discovery of quantum scarring, where some eigenstates are concentrated along unstable periodic orbits. Despite that, it is widely accepted that most eigenstates of chaotic models are indeed ergodic. Our results show instead that all eigenstates of the chaotic Dicke model are actually scarred. They also show that even the most random states of this interacting atom-photon system never occupy more than half of the available phase space. Quantum ergodicity is achievable only as an ensemble property, after temporal averages are performed.

It is generally believed that most eigenstates of quantum chaotic models are ergodic. In this work, the authors disprove this by showing that all eigenstates of the Dicke model in the chaotic regime are scarred, and that ergodicity is an ensemble property, achievable only in the temporal average.