State Quantum Phase Transitions Research Articles

Spacing distribution for quantum Rabi models * *Dedicated to Mr. Hiromichi Shinohara, who was instrumental in establishing NTT-IFM, on his 70th birthday.

The asymmetric quantum Rabi model (AQRM) is a fundamental model in quantum optics describing the interaction of light and matter. Besides its immediate physical interest, the AQRM possesses an intriguing mathematical structure which is far from being completely understood. In this paper, we focus on the distribution of the level spacing, the difference between consecutive eigenvalues of the AQRM in the limit of high energies, i.e. large quantum numbers. In the symmetric case, that is the quantum Rabi model (QRM), the spacing distribution for each parity (given by the Z2 -symmetry) is fully clarified by an asymptotic expression derived by de Monvel and Zielinski, though some questions remain for the full spectrum spacing. However, in the general AQRM case, there is no parity decomposition for the eigenvalues. In connection with numerically exact studies and recent theoretical results, we describe the spacing distribution for the AQRM which is characterized by a new type of periodicity and symmetric behavior of the distribution with respect to the bias parameter, reflecting the hidden symmetry of the AQRM. In addition, we observe in the AQRM the excited state quantum phase transition for large values of the bias parameter, analogous to the QRM with large qubit energy, and an internal symmetry of the level spacing distribution for fixed bias. This novel symmetry is independent from the symmetry for half-integer bias and not explained by current theoretical knowledge.

Generalized quantum geometric tensor for excited states using the path integral approach

The quantum geometric tensor, composed of the quantum metric tensor and Berry curvature, fully encodes the parameter space geometry of a physical system. We first provide a formulation of the quantum geometrical tensor in the path integral formalism that can handle both the ground and excited states, making it useful to characterize excited state quantum phase transitions (ESQPT). In this setting, we also generalize the quantum geometric tensor to incorporate variations of the system parameters and the phase-space coordinates. This gives rise to an alternative approach to the quantum covariance matrix, from which we can get information about the quantum entanglement of Gaussian states through tools such as purity and von Neumann entropy. Second, we demonstrate the equivalence between the formulation of the quantum geometric tensor in the path integral formalism and other existing methods. Furthermore, we explore the geometric properties of the generalized quantum metric tensor in depth by calculating the Ricci tensor and scalar curvature for several quantum systems, providing insight into this geometric information.

Spectral kissing and its dynamical consequences in the squeeze-driven Kerr oscillator

Transmon qubits are the predominant element in circuit-based quantum information processing, such as existing quantum computers, due to their controllability and ease of engineering implementation. But more than qubits, transmons are multilevel nonlinear oscillators that can be used to investigate fundamental physics questions. Here, they are explored as simulators of excited state quantum phase transitions (ESQPTs), which are generalizations of quantum phase transitions to excited states. We show that the spectral kissing (coalescence of pairs of energy levels) experimentally observed in the effective Hamiltonian of a driven SNAIL-transmon is an ESQPT precursor. We explore the dynamical consequences of the ESQPT, which include the exponential growth of out-of-time-ordered correlators, followed by periodic revivals, and the slow evolution of the survival probability due to localization. These signatures of ESQPT are within reach for current superconducting circuits platforms and are of interest to experiments with cold atoms and ion traps.

Phase transitions of the anisotropic Dicke model

We systematically analyze the various phase transitions of the anisotropic Dicke model that is endowed with both rotating and counter-rotating light-matter couplings. In addition to the ground state quantum phase transition (QPT) from the normal to the super-radiant phase, the anisotropic Dicke model also exhibits other transitions namely the excited state quantum phase transition (ES- QPT), ergodic to non-ergodic transition (ENET) and the temperature dependent phase transition. We show that these phase transitions are profitably studied not only with the standard consecutive level spacing ratio, but also with the aid of various eigenvector quantities such as von Neumann entanglement entropy, the participation ratio, multifractal dimension and mutual information. For ENET, both the statics and dynamics of the participation ratio offer a consistent and useful picture. An exciting finding from our work is that the ESQPT and the ENET are closely related to each other. We show this with the aid of two characteristic energies in the spectrum corresponding to jumps in von Neumann entropy.

Spin-nematic squeezing for dynamical quantum phase transitions in a spinor Bose-Einstein condensate

We study spin-nematic squeezing in a spinor Bose-Einstein condensate and use it to probe dynamical quantum phase transitions. We show that the quench dynamics of the spin-nematic squeezing exhibits a nonanalytical change with respect to a final quadratic Zeeman energy, and demonstrate that the squeezing detects not only the ground state phase diagram, but also the excited state quantum phase transition. We further analyze the dynamic stabilization of the system and present its relation to the dynamical quantum phase transitions. Our results are applicable to the spinor condensates with both of the sodium and rubidium atoms.

Floquet time crystals in driven spin systems with all-to-all p -body interactions

We show the emergence of Floquet time crystal (FTC) phases in the Floquet dynamics of periodically driven $p$-spin models, which describe a collection of spin-1/2 particles with all-to-all $p$-body interactions. Given the mean-field nature of these models, we treat the problem exactly in the thermodynamic limit and show that, for a given $p$, these systems can host various robust time-crystalline responses with period $nT$, where $T$ is the period of the drive and $n$ an integer between 2 and $p$. In particular, the case of four-body interactions ($p=4$) gives rise to both a usual period-doubling crystal, and also a novel period-quadrupling phase. We develop a comprehensive framework to predict robust subharmonic response in classical area-preserving maps, and use this as a basis to predict the occurrence and characterize the stability of the resulting mean-field FTC phases in the quantum regime. Our analysis reveals that the robustness of the time-crystal behavior is reduced as their period increases, and establishes a connection between the emergence of time crystals, described by eigenstate ordering and robust subharmonic response, and the phenomenology of excited state and dynamical quantum phase transitions. Finally, for the models hosting two or more coexisting time crystal phases, we define protocols where the periodic subharmonic response of the system can be varied in time via the non-periodic modulation of an external control parameter.

Revisiting the phase transitions of the Dicke model

The Dicke model exhibits a variety of phase transitions. The quantum phase transition from the normal phase to the super-radiant phase is marked by a dramatic change in the scaling of the participation ratio. We find that the ground state in the super-radiant phase exhibits multifractality manifest in the participation ratio scaling as the square root of the full Hilbert space dimension. The thermal phase transition temperature, for which we obtain an exact analytical expression, is strikingly captured by the mutual information between two spins. In the excited state quantum phase transition within the super-radiant phase, we discover a new upper cut-off energy; the central energy band between the lower and upper cut-off energies shows distinctly different behaviour. This finding is corroborated with the aid of several eigenvector properties: von Neumann entanglement entropy between spins and bosons, the mean photon number, concurrence between two spins, and participation ratio. Thus we obtain a unified picture for the three different kinds of phase transitions.

Anharmonicity-induced excited-state quantum phase transition in the symmetric phase of the two-dimensional limit of the vibron model

In most cases, excited state quantum phase transitions can be associated with the existence of critical points (local extrema or saddle points) in a system's classical limit energy functional. However, an excited-state quantum phase transition might also stem from the lowering of the asymptotic energy of the corresponding energy functional. One such example occurs in the 2D limit of the vibron model, once an anharmonic term in the form of a quadratic bosonic number operator is added to the Hamiltonian. The study of this case in the broken-symmetry phase was presented in Phys. Rev. A. 81 050101 (2010). In the present work, we delve further into the nature of this excited-state quantum phase transition and we characterize it in the, previously overlooked, symmetric phase of the model making use of quantities such as the effective frequency, the expected value of the quantum number operator, the participation ratio, the density of states, and the quantum fidelity susceptibility. In addition to this, we extend the usage of the quasilinearity parameter, introduced in molecular physics, to characterize the phases in the spectrum of the anharmonic 2D limit of the vibron model and a down-to-earth analysis has been included with the characterization of the critical energies for the linear isomers HCN/HNC.

Chaos in a deformed Dicke model

The critical behavior in an important class of excited state quantum phase transitions is signaled by the presence of a new constant of motion only at one side of the critical energy. We study the impact of this phenomenon in the development of chaos in a modified version of the paradigmatic Dicke model of quantum optics, in which a perturbation is added that breaks the parity symmetry. Two asymmetric energy wells appear in the semiclassical limit of the model, whose consequences are studied both in the classical and in the quantum cases. Classically, Poincaré sections reveal that the degree of chaos not only depends on the energy of the initial condition chosen, but also on the particular energy well structure of the model. In the quantum case, Peres lattices of physical observables show that the appearance of chaos critically depends on the quantum conserved number provided by this constant of motion. The conservation law defined by this constant is shown to allow for the coexistence between chaos and regularity at the same energy. We further analyze the onset of chaos in relation with an additional conserved quantity that the model can exhibit.

Quantum fidelity susceptibility in excited state quantum phase transitions: application to the bending spectra of nonrigid molecules

We characterize excited state quantum phase transitions in the two dimensional limit of the vibron model with the quantum fidelity susceptibility, comparing the obtained results with the information provided by the participation ratio. As an application, we locate the eigenstate closest to the barrier to linearity and determine the linear or bent character of the different overtones for particular bending modes of six molecular species. We perform a fit and use the optimized eigenvalues and eigenstates in three cases and make use of recently published results for the other three cases.

Quantum geometric tensor and quantum phase transitions in the Lipkin-Meshkov-Glick model

We study the quantum metric tensor and its scalar curvature for a particular version of the Lipkin-Meshkov-Glick model. We build the classical Hamiltonian using Bloch coherent states and find its stationary points. They exhibit the presence of a ground state quantum phase transition, where a bifurcation occurs, showing a change of stability associated with an excited state quantum phase transition. Symmetrically, for a sign change in one Hamiltonian parameter, the same phenomenon is observed in the highest energy state. Employing the Holstein-Primakoff approximation, we derive analytic expressions for the quantum metric tensor and compute the scalar and Berry curvatures. We contrast the analytic results with their finite-size counterparts obtained through exact numerical diagonalization and find an excellent agreement between them for large sizes of the system in a wide region of the parameter space, except in points near the phase transition where the Holstein-Primakoff approximation ceases to be valid.

Work statistics and symmetry breaking in an excited-state quantum phase transition

We examine how the presence of an excited-state quantum phase transition manifests in the dynamics of a many-body system subject to a sudden quench. Focusing on the Lipkin-Meshkov-Glick model initialized in the ground state of the ferromagnetic phase, we demonstrate that the work probability distribution displays non-Gaussian behavior for quenches in the vicinity of the excited-state critical point. Furthermore, we show that the entropy of the diagonal ensemble is highly susceptible to critical regions, making it a robust and practical indicator of the associated spectral characteristics. We assess the role that symmetry breaking has on the ensuing dynamics, highlighting that its effect is only present for quenches beyond the critical point. Finally, we show that similar features persist when the system is initialized in an excited state and briefly explore the behavior for initial states in the paramagnetic phase.

Excited state quantum phase transitions in the bending spectra of molecules

We present an extension of the Hamiltonian of the two dimensional limit of the vibron model to encompass all possible interactions up to four-body operators. We apply this Hamiltonian to the modeling of the bending spectrum of four molecules: HNC, H2S, Si2C, and NCNCS. The selected molecular species include linear, bent, and nonrigid equilibrium structures, proving the versatility of the algebraic approach which allows for the consideration of utterly different physical cases within a single Hamiltonian and a general formalism. For each case we compute predicted bending energies and wave functions, that we use to depict the associated quantum monodromy diagram, Birge-Sponer plot, and participation ratio. In nonrigid cases, we also show the bending energy functional obtained using the coherent –or intrinsic– state formalism.

Chaos and quantum scars in a coupled top model

We consider a coupled top model describing two interacting large spins, which is studied semiclassically as well as quantum mechanically. This model exhibits a variety of interesting phenomena such as a quantum phase transition (QPT), a dynamical transition, and excited-state quantum phase transitions above a critical coupling strength. Both classical dynamics and entanglement entropy reveal ergodic behavior at the center of the energy density band for an intermediate range of coupling strength above QPT, where the level spacing distribution changes from Poissonian to Wigner-Dyson statistics. Interestingly, in this model we identify quantum scars as reminiscent of unstable collective dynamics even in the presence of an interaction. The statistical properties of such scarred states deviate from the ergodic limit corresponding to the random matrix theory and violate Berry's conjecture. In contrast to ergodic evolution, the oscillatory behavior in the dynamics of the unequal time commutator and survival probability is observed as the dynamical signature of a quantum scar, which can be relevant for its detection.

Measuring out-of-time-ordered correlation functions with a single impurity qubit in a bosonic Josephson junction

We calculate the out-of-time-ordered correlation function (OTOC) of a single impurity qubit coupled to a fully connected many-particle system, such as a bosonic Josephson junction, or spins with long-range interactions. In these systems, the qubit OTOC can be used to detect both ground state and excited state quantum phase transitions (QPTs), making it a robust order parameter that is considerably more sensitive than the standard one-body correlation function. Finite size scaling exponents for an N body system can also be accurately extracted from the long-time OTOC dynamics, however, for short times there is a discrepancy due to the fact that the qubit has not had enough time to couple to the larger system. Our results show that the OTOC of even the smallest probe is enough to diagnose a QPT in fully connected models but, like a continuous measurement, can still cause a backaction effect which leads to weakly chaotic dynamics and gradual information scrambling.

Collapse of the Kondo state and ferromagnetic quantum phase transition in YbFe2Zn20

We present the electrical resistivity data under application of pressures up to $\sim$ 26 GPa and down to 50 mK temperatures on YbFe$_2$Zn$_{20}$. We find a pressure induced magnetic phase transition with an onset at $p_c$=18.2$\pm$0.8 GPa. At ambient pressure, YbFe$_2$Zn$_{20}$ manifests a heavy fermion, nonmagnetic ground state and the Fermi liquid behavior at low temperatures. As pressure is increased, the power law exponent in resistivity, $n$, deviates significantly from Fermi liquid behavior and tends to saturate with $n$ = 1 near $p_c$. A pronounced resistivity maximum, $T_\text{max}$, which scales with Kondo temperature is observed. $T_\text{max}$ decreases with increasing pressure and flattened out near $p_c$ indicating the suppression of Kondo exchange interaction. For $p>p_c$, $T_\text{max}$ shows a sudden upward shift, most likely becoming associated with crystal electric field scattering. Application of magnetic field for $p>p_c$ broadens the transition and shifts it toward the higher temperature, which is a typical behavior of the ferromagnetic transition. The magnetic transition appears to abruptly develop above $p_c$, suggesting probable first-order (with changing pressure) nature of the transition; once stabilized, the ordering temperature does not depend on pressure up to $\sim$ 26 GPa. Taken as a whole, these data suggest that YbFe$_2$Zn$_{20}$ has a quantum phase transition at $p_c$ = 18.2 GPa associated with the avoided quantum criticality in metallic ferromagnets.

Excited state quantum phase transitions and monodromy

Quantum phase transitions (QPT) and excited state quantum phase transitions (ESQPT) in algebraic models are briefly discussed, in particular those occurring in the U(3) model of bending vibrations of molecules. Experimental evidence in the water molecule, H2O, is presented.

Nonadiabatic dynamics of the excited states for the Lipkin-Meshkov-Glick model.

We theoretically investigate the impact of the excited state quantum phase transition on the adiabatic dynamics for the Lipkin-Meshkov-Glick model. Using a time-dependent protocol, we continuously change a model parameter and then discuss the scaling properties of the system especially close to the excited state quantum phase transition where we find that these depend on the energy eigenstate. On top, we show that the mean-field dynamics with the time-dependent protocol gives the correct scaling and expectation values in the thermodynamic limit even for the excited states.

Excited-state quantum phase transitions in a two-fluid Lipkin model

Background: Composed systems have became of great interest in the framework of the ground state quantum phase transitions (QPTs) and many of their properties have been studied in detail. However, in these systems the study of the so called excited-state quantum phase transitions (ESQPTs) have not received so much attention. Purpose: A quantum analysis of the ESQPTs in the two-fluid Lipkin model is presented in this work. The study is performed through the Hamiltonian diagonalization for selected values of the control parameters in order to cover the most interesting regions of the system phase diagram. [Method:] A Hamiltonian that resembles the consistent-Q Hamiltonian of the interacting boson model (IBM) is diagonalized for selected values of the parameters and properties such as the density of states, the Peres lattices, the nearest-neighbor spacing distribution, and the participation ratio are analyzed. Results: An overview of the spectrum of the two-fluid Lipkin model for selected positions in the phase diagram has been obtained. The location of the excited-state quantum phase transition can be easily singled out with the Peres lattice, with the nearest-neighbor spacing distribution, with Poincar\'e sections or with the participation ratio. Conclusions: This study completes the analysis of QPTs for the two-fluid Lipkin model, extending the previous study to excited states. The ESQPT signatures in composed systems behave in the same way as in single ones, although the evidences of their presence can be sometimes blurred. The Peres lattice turns out to be a convenient tool to look into the position of the ESQPT and to define the concept of phase in the excited states realm.

Excited-state quantum phase transitions studied from a non-Hermitian perspective

A main distinguishing feature of non-Hermitian quantum mechanics is the presence of exceptional points (EPs). They correspond to the coalescence of two energy levels and their respective eigenvectors. Here, we use the Lipkin-Meshkov-Glick (LMG) model as a testbed to explore the strong connection between EPs and the onset of excited state quantum phase transitions (ESQPTs). We show that for finite systems, the exact degeneracies (EPs) obtained with the non-Hermitian LMG Hamiltonian continued into the complex plane are directly linked with the avoided crossings that characterize the ESQPTs for the real (physical) LMG Hamiltonian. The values of the complex control parameter $\alpha$ that lead to the EPs approach the real axis as the system size $N\rightarrow \infty$. This happens for both, the EPs that are close to the separatrix that marks the ESQPT and also for those that are far away, although in the latter case, the rate the imaginary part of $\alpha$ reduces to zero as $N$ increases is smaller. With the method of Pad\'e approximants, we can extract the critical value of $\alpha$.