Quantum localization of chaotic eigenstates and the level spacing distribution

Correspondence Principle, Ergodicity, and Finite-Time Dynamics

The quantum-classical correspondence stands out as one of the most intriguing challenges in the realm of quantum mechanics. Within the context of quantum chaos, the conventional approach, which relies on spectral statistics to infer distinct classical dynamical properties, fails in typical scenarios, such as slow relaxation processes, dynamical localization, or other deviations from ergodicity. To address this difficulty, we propose a novel approach, the essence of which lies in associating a quantum energy shell with classical finite-time trajectories. Our results extend the knowledge of quantum-classical correspondence, offer explanations for previously conflicting findings, and also bear significance for applications in quantum metrology, sensing, and computation when the determination of the ergodicity degree of the quantum states is in need.

Chaos and quantization of the three-particle generic Fermi-Pasta-Ulam-Tsingou model. I. Density of states and spectral statistics

We study the mixed-type classical dynamics of the three-particle Fermi-Pasta-Ulam-Tsingou (FPUT) model in relationship with its quantum counterpart and present new results on aspects of quantum chaos in this system. First we derive for the general N-particle FPUT system the transformation to the normal mode representation. Then we specialize to the three-particle FPUT case and derive analytically the semiclassical energy density of states, and its derivatives in which different singularies are determined, using the Thomas-Fermi rule. The result agrees with the numerical energy density from the Krylov subspace method, as well as with the energy density obtained by the method of quantum typicality. Here, in paper I, we concentrate on the energy level statistics (level spacing and spacing ratios), in all classical dynamical regimes of interest: the almost entirely regular, the entirely chaotic, and the mixed-type regimes. We clearly confirm, correspondingly, the Poissonian statistics, the Gaussian orthogonal ensemble statistics, and the Berry-Robnik-Brody (BRB) statistics in the mixed-type regime. It is found that the BRB level spacing distribution perfectly fits the numerical data. The extracted quantum Berry-Robnik parameter is found to agree with the classical value within better than one percent. We discuss the role of localization of chaotic eigenstates, and its appearances, in relation to the classical phase space structure (Poincaré and smaller alignment index plots), whose details will be presented in paper II, where the structure and the statistical properties of the Husimi functions in the quantum phase space will be studied.

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.

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.

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.

Entropic comparison of Landau-Zener and Demkov interactions in the phase space of a quadrupole billiard

We investigate two types of avoided crossings in a chaotic billiard within the framework of information theory. The Shannon entropy in the phase space for the Landau-Zener interaction increases as the center of the avoided crossing is approached, whereas for the Demkov interaction, the Shannon entropy decreases as the center of avoided crossing is passed by with an increase in the deformation parameter. This feature can provide a new indicator for scar formation. In addition, it is found that the Fisher information of the Landau-Zener interaction is significantly larger than that of the Demkov interaction.

Time-reversal-invariant hexagonal billiards with a point symmetry

A biparametric family of hexagonal billiards enjoying the C_{3} point symmetry is introduced and numerically investigated. First, the relative measure r(ℓ,θ;t) in a reduced phase space was mapped onto the parameter plane ℓ×θ for discrete time t up to 10^{8} and averaged in tens of randomly chosen initial conditions in each billiard. The resulting phase diagram allowed us to identify fully ergodic systems in the set. It is then shown that the absolute value of the position autocorrelation function decays like |C_{q}(t)|∼t^{-σ}, with 0<σ⩽1 in the hexagons. Following previous examples of irrational triangles, we were able to find billiards for which σ∼1. This is further evidence that, although not chaotic (all Lyapunov exponents are zero), billiards in polygons might exhibit a near strongly mixing dynamics in the ergodic hierarchy. Quantized counterparts with distinct classical properties were also characterized. Spectral properties of singlets and doublets of the quantum billiards were investigated separately well beyond the ground state. As a rule of thumb, for both singlet and doublet sequences, we calculate the first 120 000 energy eigenvalues in a given billiard and compute the nearest neighbor spacing distribution p(s), as well as the cumulative spacing function I(s)=∫_{0}^{s}p(s^{'})ds^{'}, by considering the last 20 000 eigenvalues only. For billiards with σ∼1, we observe the results predicted for chaotic geometries by Leyvraz, Schmit, and Seligman, namely, a Gaussian unitary ensemble behavior in the degenerate subspectrum, in spite of the presence of time-reversal invariance, and a Gaussian orthogonal ensemble behavior in the singlets subset. For 0<σ<1, formulas for intermediate quantum statistics have been derived for the doublets following previous works by Brody, Berry and Robnik, and Bastistić and Robnik. Different regimes in a given energy spectrum have been identified through the so-called ergodic parameter α=t_{H}/t_{C}, the ratio between the Heisenberg time and the classical diffusive-like transport time, which signals the possibility of quantum dynamical localization when α<1. A good quantitative agreement is found between the appropriate formulas with parameters extracted from the classical phase space and the data from the calculated quantum spectra. A rich variety of standing wave patterns and corresponding Poincaré-Husimi representations in a reduced phase space are reported, including those associated with lattice modes, scarring, and high-frequency localization phenomena.

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.

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.

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

The quantum localization is one of the remarkable phenomena in the studies of quantum chaos and plays an important role in various contexts. Thus, an understanding of the properties of quantum localization is essential. In spite of much effort dedicated to investigating the manifestations of localization in the time-dependent systems, the features of localization in time-independent systems are still less explored, particularly in quantum systems which correspond to the classical systems with smooth Hamiltonian. In this work, we present such a study for a quantum many-body system, namely, the Dicke model. The classical counterpart of the Dicke model is given by a smooth Hamiltonian with two degrees of freedom. We examine the signatures of localization in its chaotic eigenstates. We show that the entropy localization measure, which is defined in terms of the information entropy of Husimi distribution, behaves linearly with the participation number, a measure of the degree of localization of a quantum state. We further demonstrate that the localization measure probability distribution is well described by the β distribution. We also find that the averaged localization measure is linearly related to the level repulsion exponent, a widely used quantity to characterize the localization in chaotic eigenstates. Our findings extend the previous results in billiards to the quantum many-body system with classical counterpart described by a smooth Hamiltonian, and they indicate that the properties of localized chaotic eigenstates are universal.

Statistical properties of the localization measure of chaotic eigenstates and the spectral statistics in a mixed-type billiard

We study the quantum localization in the chaotic eigenstates of a billiard with mixed-type phase space [J. Phys. A: Math. Gen. 16, 3971 (1983)JPHAC50305-447010.1088/0305-4470/16/17/014; J. Phys. A: Math. Gen. 17, 1049 (1984)JPHAC50305-447010.1088/0305-4470/17/5/027], after separating the regular and chaotic eigenstates, in the regime of slightly distorted circle billiard where the classical transport time in the momentum space is still large enough, although the diffusion is not normal. This is a continuation of our recent papers [Phys. Rev. E 88, 052913 (2013)PLEEE81539-375510.1103/PhysRevE.88.052913; Phys. Rev. E 98, 022220 (2018)2470-004510.1103/PhysRevE.98.022220]. In quantum systems with discrete energy spectrum the Heisenberg time t_{H}=2πℏ/ΔE, where ΔE is the mean level spacing (inverse energy level density), is an important timescale. The classical transport timescale t_{T} (transport time) in relation to the Heisenberg timescale t_{H} (their ratio is the parameter α=t_{H}/t_{T}) determines the degree of localization of the chaotic eigenstates, whose measure A is based on the information entropy. We show that A is linearly related to normalized inverse participation ratio. The localization of chaotic eigenstates is reflected also in the fractional power-law repulsion between the nearest energy levels in the sense that the probability density (level spacing distribution) to find successive levels on a distance S goes like ∝S^{β} for small S, where 0≤β≤1, and β=1 corresponds to completely extended states. We show that the level repulsion exponent β is empirically a rational function of α, and the mean 〈A〉 (averaged over more than 1000 eigenstates) as a function of α is also well approximated by a rational function. In both cases there is some scattering of the empirical data around the mean curve, which is due to the fact that A actually has a distribution, typically with quite complex structure, but in the limit α→∞ well described by the beta distribution. The scattering is significantly stronger than (but similar as) in the stadium billiard [Nonlin. Phenom. Complex Syst. (Minsk) 21, 225 (2018)] and the kicked rotator [Phys. Rev. E 91, 042904 (2015)PLEEE81539-375510.1103/PhysRevE.91.042904]. Like in other systems, β goes from 0 to 1 when α goes from 0 to ∞. β is a function of 〈A〉, similar to the quantum kicked rotator and the stadium billiard.

Structure, size, and statistical properties of chaotic components in a mixed-type Hamiltonian system

We perform a detailed study of the chaotic component in mixed-type Hamiltonian systems on the example of a family of billiards [introduced by Robnik in J. Phys. A: Math. Gen. 16, 3971 (1983)JPHAC50305-447010.1088/0305-4470/16/17/014]. The phase space is divided into a grid of cells and a chaotic orbit is iterated a large number of times. The structure of the chaotic component is discerned from the cells visited by the chaotic orbit. The fractal dimension of the border of the chaotic component for various values of the billiard shape parameter is determined with the box-counting method. The cell-filling dynamics is compared to a model of uncorrelated motion, the so-called random model [Robnik et al. J. Phys. A: Math. Gen. 30, L803 (1997)JPHAC50305-447010.1088/0305-4470/30/23/003], and deviations attributed to sticky objects in the phase space are found. The statistics of the number of orbit visits to the cells is analyzed and found to be in agreement with the random model in the long run. The stickiness of the various structures in the phase space is quantified in terms of the cell recurrence times. The recurrence time distributions in a few selected cells as well as the mean and standard deviation of recurrence times for all cells are analyzed. The standard deviation of cell recurrence time is found to be a good quantifier of stickiness on a global scale. Three methods for determining the measure of the chaotic component are compared and the measure is calculated for various values of the billiard shape parameter. Lastly, the decay of correlations and the diffusion of momenta is analyzed.

Exact distribution of spacing ratios for random and localized states in quantum chaotic systems

Typical eigenstates of quantum systems, whose classical limit is chaotic, are well approximated as random states. Corresponding eigenvalue spectra are modeled through an appropriate ensemble described by random matrix theory. However, a small subset of states violates this principle and displays eigenstate localization, a counterintuitive feature known to arise due to purely quantum or semiclassical effects. In the spectrum of chaotic systems, the localized and random states interact with one another and modify the spectral statistics. In this work, a 3×3 random matrix model is used to obtain exact results for the ratio of spacing between a generic and localized state. We consider time-reversal-invariant as well as noninvariant scenarios. These results agree with the spectra computed from realistic physical systems that display localized eigenmodes.

Aspects of diffusion in the stadium billiard

We perform a detailed numerical study of diffusion in the ɛ stadium of Bunimovich, and propose an empirical model of the local and global diffusion for various values of ɛ with the following conclusions: (i) the diffusion is normal for all values of ɛ (≤0.3) and all initial conditions, (ii) the diffusion constant is a parabolic function of the momentum (i.e., we have inhomogeneous diffusion), (iii) the model describes the diffusion very well including the boundary effects, (iv) the approach to the asymptotic equilibrium steady state is exponential, (v) the so-called random model (Robnik et al., 1997) is confirmed to apply very well, (vi) the diffusion constant extracted from the distribution function in momentum space and the one derived from the second moment agree very well. The classical transport time, an important parameter in quantum chaos, is thus determined.

Statistical properties of the localization measure in a finite-dimensional model of the quantum kicked rotator

We study the quantum kicked rotator in the classically fully chaotic regime K=10 and for various values of the quantum parameter k using Izrailev's N-dimensional model for various N≤3000, which in the limit N→∞ tends to the exact quantized kicked rotator. By numerically calculating the eigenfunctions in the basis of the angular momentum we find that the localization length L for fixed parameter values has a certain distribution; in fact, its inverse is Gaussian distributed, in analogy and in connection with the distribution of finite time Lyapunov exponents of Hamilton systems. However, unlike the case of the finite time Lyapunov exponents, this distribution is found to be independent of N and thus survives the limit N=∞. This is different from the tight-binding model of Anderson localization. The reason is that the finite bandwidth approximation of the underlying Hamilton dynamical system in the Shepelyansky picture [Phys. Rev. Lett. 56, 677 (1986)] does not apply rigorously. This observation explains the strong fluctuations in the scaling laws of the kicked rotator, such as the entropy localization measure as a function of the scaling parameter Λ=L/N, where L is the theoretical value of the localization length in the semiclassical approximation. These results call for a more refined theory of the localization length in the quantum kicked rotator and in similar Floquet systems, where we must predict not only the mean value of the inverse of the localization length L but also its (Gaussian) distribution, in particular the variance. In order to complete our studies we numerically analyze the related behavior of finite time Lyapunov exponents in the standard map and of the 2×2 transfer matrix formalism. This paper extends our recent work [Phys. Rev. E 87, 062905 (2013)].

Exponential Fermi acceleration in general time-dependent billiards

We show, that under very general conditions, a generic time-dependent billiard, for which a phase space of corresponding static (frozen) billiards is of the mixed type, exhibits the exponential Fermi acceleration in the adiabatic limit. The velocity dynamics in the adiabatic regime is represented as an integral over a path through the abstract space of invariant components of corresponding static billiards, where the paths are generated probabilistically in terms of transition-probability matrices. We study the statistical properties of possible paths and deduce the conditions for the exponential Fermi acceleration. The exponential Fermi acceleration and theoretical concepts presented in the paper are demonstrated numerically in four different time-dependent billiards.