A bstract Given a choice of an ordered, orthonormal basis for a D -dimensional Hilbert space, one can define a discrete version of the Wigner function — a quasi-probability distribution which represents any quantum state as a real, normalized function on a discrete phase space. The Wigner function, in general, takes on negative values, and the amount of negativity in the Wigner function gives an operationally meaningful measure of the complexity of simulating the quantum state on a classical computer. Further, Wigner negativity also gives a lower bound on an entropic measure of spread complexity. In this paper, we study the growth of Wigner negativity for a generic initial state under time evolution with chaotic Hamiltonians. In [1], a perturbative argument was given to show that the Krylov basis minimizes the early time growth of Wigner negativity in the large- D limit. Using tools from random matrix theory, here we show that for a generic choice of basis, the Wigner negativity for a classical initial state becomes exponentially large in an O (1) amount of time evolution. On the other hand, we show that in the Krylov basis the negativity grows at most as a power law, and becomes exponentially large only at exponential times. We take this as evidence that the Krylov basis is ideally suited for a dual, semi-classical effective description of chaotic quantum dynamics for large- D at sub-exponential times. For the Gaussian unitary ensemble, this effective description is the q → 0 limit of q -deformed JT gravity.
A bstract Given a choice of an ordered, orthonormal basis for a D -dimensional Hilbert space, one can define a discrete version of the Wigner function — a quasi-probability distribution which represents any quantum state as a real, normalized function on a discrete phase space. The Wigner function, in general, takes on negative values, and the amount of negativity in the Wigner function gives an operationally meaningful measure of the complexity of simulating the quantum state on a classical computer. Further, Wigner negativity also gives a lower bound on an entropic measure of spread complexity. In this paper, we study the growth of Wigner negativity for a generic initial state under time evolution with chaotic Hamiltonians. In [1], a perturbative argument was given to show that the Krylov basis minimizes the early time growth of Wigner negativity in the large- D limit. Using tools from random matrix theory, here we show that for a generic choice of basis, the Wigner negativity for a classical initial state becomes exponentially large in an O (1) amount of time evolution. On the other hand, we show that in the Krylov basis the negativity grows at most as a power law, and becomes exponentially large only at exponential times. We take this as evidence that the Krylov basis is ideally suited for a dual, semi-classical effective description of chaotic quantum dynamics for large- D at sub-exponential times. For the Gaussian unitary ensemble, this effective description is the q → 0 limit of q -deformed JT gravity.
A bstract For any state in a D -dimensional Hilbert space with a choice of basis, one can define a discrete version of the Wigner function — a quasi-probability distribution which represents the state on a discrete phase space. The Wigner function can, in general, take on negative values, and the amount of negativity in the Wigner function has an operational meaning as a resource for quantum computation. In this note, we study the growth of Wigner negativity for a generic initial state under time evolution with chaotic Hamiltonians. We introduce the Krylov-Wigner function, i.e., the Wigner function defined with respect to the Krylov basis (with appropriate phases), and show that this choice of basis minimizes the early time growth of Wigner negativity in the large D limit. We take this as evidence that the Krylov basis (with appropriate phases) is ideally suited for a dual, semi-classical description of chaotic quantum dynamics at large D . We also numerically study the time evolution of the Krylov-Wigner function and its negativity in random matrix theory for an initial pure state. We observe that the negativity broadly shows three phases: it rises gradually for a time of O(D) O D , then hits a sharp ramp and finally saturates close to its upper bound of D D .
For any state in a D-dimensional Hilbert space with a choice of basis, one can define a discrete version of the Wigner function -- a quasi-probability distribution which represents the state on a discrete phase space. The Wigner function can, in general, take on negative values, and the amount of negativity in the Wigner function has an operational meaning as a resource for quantum computation. In this note, we study the growth of Wigner negativity for a generic initial state under time evolution with chaotic Hamiltonians. We introduce the Krylov-Wigner function, i.e., the Wigner function defined with respect to the Krylov basis (with appropriate phases), and show that this choice of basis minimizes the early time growth of Wigner negativity in the large D limit. We take this as evidence that the Krylov basis (with appropriate phases) is ideally suited for a dual, semi-classical description of chaotic quantum dynamics at large D. We also numerically study the time evolution of the Krylov-Wigner function and its negativity in random matrix theory for an initial pure state. We observe that the negativity rises gradually for a time of O(D) and then saturates close to its upper bound of D.
The Wigner function, in general, takes on negative values, and the amount of negativity in the Wigner function gives an operationally meaningful measure of the complexity of simulating the quantum state on a classical computer. In this paper, we study the growth of Wigner negativity of the Hartle-Hawking state of Jackiw-Teitelboim gravity under time evolution. We work in the gravitational length basis, at genus zero, and at high temperature, β≪1. Our main analytic result is simple: to leading order in β the state is a Gaussian wavepacket of minimum uncertainty, sitting at rest at the turning point of the Liouville wall. Its Wigner function is therefore positive, and the negativity is 1+δ(β,t), where the correction δ is smaller than any power of β. We show that δ is an even function of time, so there is no linear growth at t=0, and that once the packet has reflected off the wall its leading-order evolution is a Clifford shear, which cannot change the negativity. We also numerically study the time evolution of the Wigner function and its negativity. We observe that δ stays below 10−11 at t=0 and through the reflection, then rises slowly over a few tens of β/π, and then approaches a late-time value consistent with δ∞(β)≃0.08e−4.4/β for 0.25≤β≤1: the lower the temperature, the larger the late-time negativity. The exact survival amplitude Z(β+it)/Z(β) fixes the spread complexity to order t4, and it also fixes the seed-normalised Wigner diagnostic of our earlier work (arXiv:2607.04065, arXiv:2607.17346). We take this as evidence that the length basis is ideally suited for a dual, semi-classical description of the dynamics at genus zero. Beyond genus zero the length basis is overcomplete, and we make no statement about finite eS0.
We propose that the Krylov basis gives a semiclassical representation of dynamics in general large-N, complex, many-body systems. As a probe of this semiclassicality, we study the growth of Wigner negativity -- a measure of the complexity of classical simulation -- under time evolution in Krylov space in several solvable models. We begin with 2d CFTs, initially in either the vacuum or the thermofield double state on a line excited by a primary operator. In both cases, Wigner negativity remains an O(1) constant and does not grow at late times, indicating approximately classical dynamics in the Krylov basis. We then study random matrix theory with the maximally entangled state between two copies as the initial state. For general one-cut matrix models, we argue that Wigner negativity in the Krylov basis grows as t1/2 at large O(1) times but does not scale with the Hilbert space dimension, thus indicating semiclassical dynamics in Krylov space. Finally, in the double-scaled SYK model, we find an approximately classical phase (constant negativity) at early times and a semiclassical phase (t1/2 growth) at late times. In all these examples, Wigner negativity either remains constant or grows slowly, demonstrating emergent semiclassicality of dynamics in Krylov space.
A bstract We consider the statistics of the results of a measurement of the spreading operator in the Krylov basis generated by the Hamiltonian of a quantum system starting from a specified initial pure state. We first obtain the probability distribution of the results of measurements of this spreading operator at a certain instant of time, and compute the characteristic function of this distribution. We show that the moments of this characteristic function are related to the so-called generalised spread complexities, and obtain expressions for them in several cases when the Hamiltonian is an element of a Lie algebra. Furthermore, by considering a continuum limit of the Krylov basis, we show that the generalised spread complexities of higher orders have a peak in the time evolution for a random matrix Hamiltonian belonging to the Gaussian unitary ensemble. We also obtain an upper bound in the change in generalised spread complexity at an arbitrary time in terms of the operator norm of the Hamiltonian and discuss the significance of these results.
The Normalised Wigner Negativity Rate as a Second-Moment Probe of Infall in AdS3
In spread complexity, the growth rate of the average position of a time-evolved state along its Krylov chain recovers the proper radial momentum of an infalling particle in AdS. It is, however, a measure of the first moment only, irrespective of the spread of the wavepacket away from its classical trajectory. The rate of a normalized Krylov-Wigner negativity can be proposed as a diagnostic of the second moment of the boundary state that captures this spreading. Starting with the seed-normalized Krylov-Wigner distribution, that is, the Wigner transform of the descendant cloud with the decaying return amplitude divided out, we obtain an analytic Bessel form in the macroscopic limit and compute its total negativity explicitly. Retaining the Bessel variable all the way through, we find that the negativity goes as sinh4Δ(πt/β). But the raw normalized-state negativity saturates to an O(1) constant, well below the O(D) bound. Using the exact negative binomial statistics of the Krylov chain, the normalized negativity is at late times a fixed power of the second moment of the Krylov wavepacket, N(t)∝[Var(N^Krylov)]Δ, for every Δ>1/2. The relation linearizes precisely at Δ=1. In this regime via the momentum dictionary of Caputa et al. (arXiv:2410.23334), the rate becomes the product of the proper radial position and momentum, R∝CPρ, which we interpret as the rate of the tidal stretch of nearby geodesics falling into the horizon. We also comment briefly on a speculative direction for future research, in particular the interpretation of the transverse string size operator in terms of the Krylov number operator through the common SU(1,1) discrete series.