Source author record

Jacobus J. M. Verbaarschot

Jacobus J. M. Verbaarschot appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2024arXiv

Quantum Chaos on Edge

In recent years, the physics of many-body quantum chaotic systems close to their ground states has come under intensified scrutiny. Such studies are motivated by the emergence of model systems exhibiting chaotic fluctuations throughout the entire spectrum (the Sachdev-Ye-Kitaev (SYK) model being a renowned representative) as well as by the physics of holographic principles, which likewise unfold close to ground states. Interpreting the edge of the spectrum as a quantum critical point, here we combine a wide range of analytical and numerical methods to the identification and comprehensive description of two different universality classes: the near edge physics of ``sparse'' and the near edge of ``dense'' chaotic systems. The distinction lies in the ratio between the number of a system's random parameters and its Hilbert space dimension, which is exponentially small or algebraically small in the sparse and dense case, respectively. Notable representatives of the two classes are generic chaotic many-body models (sparse) and single particle systems, invariant random matrix ensembles, or chaotic gravitational systems (dense). While the two families share identical spectral correlations at energy scales comparable to the level spacing, the density of states and its fluctuations near the edge are different. Considering the SYK model as a representative of the sparse class, we apply a combination of field theory and exact diagonalization to a detailed discussion of its edge spectrum. Conversely, Jackiw-Teitelboim gravity is our reference model for the dense class, where an analysis of the gravitational path integral and random matrix theory reveal universal differences to the sparse class, whose implications for the construction of holographic principles we discuss.

preprint2022arXiv

Replica Symmetry Breaking in Random Non-Hermitian Systems

Recent studies have revealed intriguing similarities between the contribution of wormholes to the gravitational path integral and the phenomenon of replica symmetry breaking observed in spin glasses and other disordered systems. Interestingly, these configurations may also be important for the explanation of the information paradox of quantum black holes. Motivated by these developments, we investigate the thermodynamic properties of a $PT$-symmetric system composed of two random non-Hermitian Hamiltonians with no explicit coupling between them. After performing ensemble averaging, we identify numerically and analytically a robust first-order phase transition in the free energy of two models with quantum chaotic dynamics: the elliptic Ginibre ensemble of random matrices and a non-Hermitian Sachdev-Ye-Kitaev (SYK) model. The free energy of the Ginibre model is temperature-independent in the low-temperature phase. The SYK model has a similar behavior for sufficiently low temperature, then it experiences a possible continuous phase transition to a phase with a temperature-dependent free energy before the first-order transition takes place at a higher temperature. We identify the order parameter of the first-order phase transition and obtain analytical expressions for the critical temperature. The mechanism behind the transition is the existence of replica symmetry breaking configurations coupling Left and Right replicas that control the low-temperature limit of the partition function. We speculate that quantum chaos may be necessary for the observed dominance of off-diagonal replica symmetry breaking configurations in the low-temperature limit.

preprint2022arXiv

Symmetry Classification and Universality in Non-Hermitian Many-Body Quantum Chaos by the Sachdev-Ye-Kitaev Model

Spectral correlations are a powerful tool to study the dynamics of quantum many-body systems. For Hermitian Hamiltonians, quantum chaotic motion is related to random matrix theory spectral correlations. Based on recent progress in the application of spectral analysis to non-Hermitian quantum systems, we show that local level statistics, which probe the dynamics around the Heisenberg time, of a non-Hermitian $q$-body Sachdev-Ye-Kitev (nHSYK) model with $N$ Majorana fermions, and its chiral and complex-fermion extensions, are also well described by random matrix theory for $q > 2$, while for $q = 2$, they are given by the equivalent of Poisson statistics. For that comparison, we combine exact diagonalization numerical techniques with analytical results obtained for some of the random matrix spectral observables. Moreover, depending on $q$ and $N$, we identify $19$ out of the $38$ non-Hermitian universality classes in the nHSYK model, including those corresponding to the tenfold way. In particular, we realize explicitly $14$ out of the $15$ universality classes corresponding to non-pseudo-Hermitian Hamiltonians that involve universal bulk correlations of classes ${\rm AI}^\dagger$ and ${\rm AII}^\dagger$, beyond the Ginibre ensembles. These results provide strong evidence of striking universal features in non-unitary many-body quantum chaos, which in all cases can be captured by nHSYK models with $q > 2$.

preprint2021arXiv

Replica Symmetry Breaking and Phase Transitions in a PT Symmetric Sachdev-Ye-Kitaev Model

We show that the low temperature phase of a conjugate pair of uncoupled, quantum chaotic, nonhermitian systems such as the Sachdev-Ye-Kitaev (SYK) model or the Ginibre ensemble of random matrices are dominated by replica symmetry breaking configurations with a nearly flat free energy that terminates in a first order phase transition. In the case of the SYK model, we show explicitly that the spectrum of the effective replica theory has a gap. These features are strikingly similar to those induced by wormholes in the gravity path integral which suggests a close relation between both configurations. For a non-chaotic SYK, the results are qualitatively different: the spectrum is gapless in the low temperature phase and there is an infinite number of second order phase transitions unrelated to the restoration of replica symmetry.

preprint2020arXiv

Spectral Fluctuations in the Sachdev-Ye-Kitaev Model

We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being the scale fluctuations for which we give a simple estimate. Our numerical results for $N=32$ show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion is obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations using a nonlinear $σ$-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlator to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the $σ$-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for $q=4$ but also for the $q=2$ SYK model while the correct result in this case should have been Poisson statistics. We analyze the number variance and spectral form factor for $N=32$ and $q=4$ numerically. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings or exponentially short times, respectively.

preprint2016arXiv

Spectral and thermodynamic properties of the Sachdev-Ye-Kitaev model

We study spectral and thermodynamic properties of the Sachdev-Ye-Kitaev model, a variant of the $k$-body embedded random ensembles studied for several decades in the context of nuclear physics and quantum chaos. We show analytically that the fourth and sixth order energy cumulants vanish in the limit of large number of particles $N \to \infty$ which is consistent with a Gaussian spectral density. However, for finite $N$, the tail of the average spectral density is well approximated by a semi-circle law. The specific heat coefficient, determined numerically from the low temperature behavior of the partition function, is consistent with the value obtained by previous analytical calculations. For energy scales of the order of the mean level spacing we show that level statistics are well described by random matrix theory. Due to the underlying Clifford algebra of the model, the universality class of the spectral correlations depends on $N$. For larger energy separations we identify an energy scale that grows with $N$, reminiscent of the Thouless energy in mesoscopic physics, where deviations from random matrix theory are observed. Our results are a further confirmation that the Sachdev-Ye-Kitaev model is quantum chaotic for all time scales. According to recent claims in the literature, this is an expected feature in field theories with a gravity-dual.

preprint2015arXiv

Dirac Spectrum of the Wilson Dirac Operator for QCD with Two Colors

We study the lattice artefacts of the Wilson Dirac operator for QCD with two colors and fermions in the fundamental representation from the viewpoint of chiral perturbation theory. These effects are studied with the help of the following spectral observables: the level density of the Hermitian Wilson Dirac operator, the distribution of chirality over the real eigenvalues, and the chiral condensate for the quenched as well as for the unquenched theory. We provide analytical expressions for all these quantities. Moreover we derive constraints for the level density of the real eigenvalues of the non-Hermitian Wilson Dirac operator and the number of additional real modes. The latter is a good measure for the strength of lattice artefacts. All computations are confirmed by Monte Carlo simulations of the corresponding random matrix theory which agrees with chiral perturbation theory of two color QCD with Wilson fermions.

preprint2015arXiv

Phase Diagram of Dynamical Twisted Mass Wilson Fermions at Finite Isospin Chemical Potential

We consider the phase diagram of twisted mass Wilson fermions of two-flavor QCD in the parameter space of the quark mass, the isospin chemical potential, the twist angle and the lattice spacing. This work extends earlier studies in the continuum and those at zero chemical potential. We evaluate the phase diagram as well as the spectrum of the (pseudo-)Goldstone bosons using the chiral Lagrangian for twisted mass Wilson fermions at non-zero isospin chemical potential. The phases are obtained from a mean field analysis. At zero twist angle we find that already an infinitesimal isospin chemical potential destroys the Aoki phase. The reason is that in this phase we have massless Goldstone bosons with a non-zero isospin charge. At finite twist angle only two different phases are present, one phase which is continuously connected to the Bose condensed phase at non-zero chemical potential and another phase which is continuously connected to the normal phase. For either zero or maximal twist the phase diagram is more complicated as the saddle point equations allow for more solutions.

preprint2014arXiv

Dirac Spectra of 2-dimensional QCD-like theories

We analyze Dirac spectra of two-dimensional QCD like theories both in the continuum and on the lattice and classify them according to random matrix theories sharing the same global symmetries. The classification is different from QCD in four dimensions because the anti-unitary symmetries do not commute with $γ_5$. Therefore in a chiral basis, the number of degrees of freedom per matrix element are not given by the Dyson index. Our predictions are confirmed by Dirac spectra from quenched lattice simulations for QCD with two or three colors with quarks in the fundamental representation as well as in the adjoint representation. The universality class of the spectra depends on the parity of the number of lattice points in each direction. Our results show an agreement with random matrix theory that is qualitatively similar to the agreement found for QCD in four dimensions. We discuss the implications for the Mermin-Wagner-Coleman theorem and put our results in the context of two-dimensional disordered systems.

preprint2013arXiv

A classification of 2-dim Lattice Theory

A unified classification and analysis is presented of two dimensional Dirac operators of QCD-like theories in the continuum as well as in a naive lattice discretization. Thereby we consider the quenched theory in the strong coupling limit. We do not only consider the case of a lattice which has an even number of lattice sites in both directions and is thus equivalent to the case of staggered fermions. We also study lattices with one or both directions with an odd parity to understand the general mechanism of changing the universality class via a discretization. Furthermore we identify the corresponding random matrix ensembles sharing the global symmetries of these QCD-like theories. Despite the Mermin-Wagner-Coleman theorem we find good agreement of lattice data with our random matrix predictions.

preprint2013arXiv

Discretization Effects in the $ε$ Domain of QCD

At nonzero lattice spacing the QCD partition function with Wilson quarks undergoes either a second order phase transition to the Aoki phase for decreasing quark mass or shows a first order jump when the quark mass changes sign. We discuss these phase transitions in terms of Wilson Dirac spectra and show that the first order scenario can only occur in the presence of dynamical quarks while in the quenched case we can only have a transition to the Aoki phase. The exact microscopic spectral density of the non-Hermitian Wilson Dirac operator with dynamical quarks is discussed as well. We conclude with some remarks on discretization effects for the overlap Dirac operator.

preprint2013arXiv

Random Matrix Models for the Hermitian Wilson-Dirac operator of QCD-like theories

We introduce Random Matrix Models for the Hermitian Wilson-Dirac operator of QCD-like theories. We show that they are equivalent to the $ε$-limit of the chiral Lagrangian for Wilson chiral perturbation theory. Results are obtained for two-color QCD with quarks in the fundamental representation of the color group as well as any-color QCD with quarks in the adjoint representation. For $N_c=2$ we also have obtained the lattice spacing dependence of the quenched average spectral density for a fixed value of the index of the Dirac operator. Comparisons with direct numerical simulations of the random matrix ensemble are shown.

preprint2013arXiv

Spectral Properties of the Wilson Dirac Operator and random matrix theory

Random Matrix Theory has been successfully applied to lattice Quantum Chromodynamics. In particular, a great deal of progress has been made on the understanding, numerically as well as analytically, of the spectral properties of the Wilson Dirac operator. In this paper, we study the infra-red spectrum of the Wilson Dirac operator via Random Matrix Theory including the three leading order $a^2$ correction terms that appear in the corresponding chiral Lagrangian. A derivation of the joint probability density of the eigenvalues is presented. This result is used to calculate the density of the complex eigenvalues, the density of the real eigenvalues and the distribution of the chiralities over the real eigenvalues. A detailed discussion of these quantities shows how each low energy constant affects the spectrum. Especially we consider the limit of small and large (which is almost the mean field limit) lattice spacing. Comparisons with Monte Carlo simulations of the Random Matrix Theory show a perfect agreement with the analytical predictions. Furthermore we present some quantities which can be easily used for comparison of lattice data and the analytical results.

preprint2013arXiv

The Effect of the Low Energy Constants on the Spectral Properties of the Wilson Dirac Operator

We summarize recent analytical results obtained for lattice artifacts of the non-Hermitian Wilson Dirac operator. Hereby we discuss the effect of all three low energy constants. In particular we study the limit of small lattice spacing and also consider the regime of large lattice spacing which is closely related to the mean field limit. Thereby we extract simple relations between measurable quantities like the average number of additional real modes and the low energy constants. These relations may improve the fitting for the low energy constants.

preprint2012arXiv

Eigenvalue Density of the non-Hermitian Wilson Dirac Operator

We find the lattice spacing dependence of the eigenvalue density of the non-Hermitian Wilson Dirac operator in the $ε$-domain. The starting point is the joint probability density of the corresponding random matrix theory. In addition to the density of the complex eigenvalues we also obtain the density of the real eigenvalues separately for positive and negative chiralities as well as an explicit analytical expression for the number of additional real modes.

preprint2012arXiv

Random Matrix Models for Dirac Operators at finite Lattice Spacing

We study discretization effects of the Wilson and staggered Dirac operator with $N_{\rm c}>2$ using chiral random matrix theory (chRMT). We obtain analytical results for the joint probability density of Wilson-chRMT in terms of a determinantal expression over complex pairs of eigenvalues, and real eigenvalues corresponding to eigenvectors of positive or negative chirality as well as for the eigenvalue densities. The explicit dependence on the lattice spacing can be readily read off from our results which are compared to numerical simulations of Wilson-chRMT. For the staggered Dirac operator we have studied random matrices modeling the transition from non-degenerate eigenvalues at non-zero lattice spacing to degenerate ones in the continuum limit.