© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2601.13872 · Journal of High Energy Physics

A phase space approach to the wavefunction spreading and operator growth in the Krylov basis

Kunal Pal, Kuntal Pal, Keun-Young Kim

quant-phSubmitted Jan 20, 2026• Announced Jan 21, 2026• 1 citation• Revised Mar 17, 2026 (v2)

v1 · Jan 20, 2026 v2 · Mar 17, 2026

Abstract

In the Wigner-Weyl phase space formulation of quantum mechanics, we analyse the problem of the spreading of an initial state or an initial operator under time evolution when described in terms of the Krylov basis. After constructing the phase space functions corresponding to the Krylov basis states generated by a Hamiltonian from a given initial state by using the Weyl transformation, we subsequently use them to cast the Krylov state complexity as an integral over the phase space in terms of the Wigner function of the time-evolved initial state, so that the contribution of the classical Liouville equation and higher-order quantum corrections to the Wigner function time evolution equation towards the Krylov state complexity can be identified. Next, we construct the double phase space functions associated with the Krylov basis for operators by using a suitable generalisation of the Weyl transformation applicable for superoperators, and use them to rewrite the Krylov operator complexity as an integral over the double phase space in terms of a generalisation of the usual Wigner function. These results, in particular, show that the complexity measures based on the expansion of a time-evolved state (or an operator) in the Krylov basis can be thought to belong to a general class of complexity measures constructed from the expansion coefficients of the time-dependent Wigner function in an orthonormal basis in the phase space, and help us to connect these complexity measures with measures of complexity of time evolved state based on harmonic expansion of the time-dependent Wigner function.

Comments: V2. Added section 4.1

Related papers

Statistics and complexity of wavefunction spreading in quantum dynamical systems

Y. W. Fu, Keun-young Kim A, Kunal Pal +1 more

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.

arXiv:2411.09390•Journal of High Energy Physics•Jun 13, 2025•10 citations

Complexity growth and the Krylov-Wigner function

Ritam Basu, Anirban Ganguly, Souparna Nath +1 more

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\left(\sqrt{D}\right) $$ O D , then hits a sharp ramp and finally saturates close to its upper bound of $$ \sqrt{D} $$ D .

Journal of High Energy Physics•May 23, 2024•18 citations

On Minimizing Krylov Complexity Using Higher-Order Generators

Saud Čindrak, Kathy Lüdge

Krylov complexity provides a powerful framework for characterizing the dynamical evolution of quantum systems through the spreading of states in Krylov space. The motivation for this is rooted in the optimality of the Krylov basis for the analyzed cost function. In this work, we reinterpret the motivation for the Krylov basis from a dynamical perspective and show that it corresponds to a first-order approximation of the time-evolution operator. We extend this framework to higher-order generators and analytically disprove the optimality assumption by showing that an infinite-order generator can be constructed to exhibit smaller spread for arbitrary times. We propose a natural time scale for the construction of these higher-order generators and discuss results for matrices sampled from Gaussian Unitary Ensembles, demonstrating smaller Krylov complexity at all higher orders. These results extend the framework of Krylov complexity beyond the conventional Krylov basis by disproving the widely held assumption of optimality, extending the construction to higher-order generators, and introducing a physically motivated method for their construction. Our findings therefore suggest that previous statements and results concerning Krylov complexity may need to be reconsidered.

quant-pharXiv:2603.07183•Mar 10, 2026•0 citations

Complexity Growth and the Krylov-Wigner function

Ritam Basu, Anirban Ganguly, Nath, Souparna +1 more

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 $\sqrt{D}$.

arXiv:2402.13694•Feb 21, 2024•0 citations

Krylov complexity in the Schrödinger field theory

Peng-Zhang He, Hai-Qing Zhang

A bstract We investigate the Krylov complexity of Schrödinger field theories, focusing on both bosonic and fermionic systems within the grand canonical ensemble which includes a chemical potential. Krylov complexity measures operator growth in quantum systems by analyzing how operators spread within the Krylov space, a subspace of the Hilbert space spanned by successive applications of the superoperator [ H , ·] on an initial operator. Using the Lanczos algorithm, we construct an orthonormal Krylov basis and derive the Lanczos coefficients, which govern the operator connectivity and thus characterize the complexity. Our study reveals that the Lanczos coefficients { b n } are almost independent of the chemical potential, while { a n } are dependent on the chemical potential. Both { a n } and { b n } show linear relationships with respect to n . For both bosonic and fermionic systems, the Krylov complexities behave similarly over time, especially at late times, due to the analogous profiles of the squared absolute values of their autocorrelation functions | φ 0 ( t )| 2 . The Krylov complexity grows exponentially with time, but its asymptotic scaling factor λ K is significantly smaller than the twice of the slope of the { b n } coefficients, contrasting to the relativistic field theories where the scaling aligns more closely with the twice of the slope of { b n }.

arXiv:2411.16302•Journal of High Energy Physics•Mar 19, 2025•4 citations

Wigner negativity in Krylov space and emergent semiclassicality

Vijay Balasubramanian, Paweł Caputa, Onkar Parrikar +1 more

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 $t^{1/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 ($t^{1/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.

arXiv:2607.01351•Jul 1, 2026•0 citations

In the Wigner-Weyl phase space formulation of quantum mechanics, we analyse the problem of the spreading of an initial state or an initial operator under time evolution when described in terms of the Krylov basis. After constructing the phase space functions corresponding to the Krylov basis states generated by a Hamiltonian from a given initial state by using the Weyl transformation, we subsequently use them to cast the Krylov state complexity as an integral over the phase space in terms of the Wigner function of the time-evolved initial state, so that the contribution of the classical Liouville equation and higher-order quantum corrections to the Wigner function time evolution equation towards the Krylov state complexity can be identified. Next, we construct the double phase space functions associated with the Krylov basis for operators by using a suitable generalisation of the Weyl transformation applicable for superoperators, and use them to rewrite the Krylov operator complexity as an integral over the double phase space in terms of a generalisation of the usual Wigner function. These results, in particular, show that the complexity measures based on the expansion of a time-evolved state (or an operator) in the Krylov basis can be thought to belong to a general class of complexity measures constructed from the expansion coefficients of the time-dependent Wigner function in an orthonormal basis in the phase space, and help us to connect these complexity measures with measures of complexity of time evolved state based on harmonic expansion of the time-dependent Wigner function.

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.

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\left(\sqrt{D}\right)

O D , then hits a sharp ramp and finally saturates close to its upper bound of

\sqrt{D}

D .

Krylov complexity provides a powerful framework for characterizing the dynamical evolution of quantum systems through the spreading of states in Krylov space. The motivation for this is rooted in the optimality of the Krylov basis for the analyzed cost function. In this work, we reinterpret the motivation for the Krylov basis from a dynamical perspective and show that it corresponds to a first-order approximation of the time-evolution operator. We extend this framework to higher-order generators and analytically disprove the optimality assumption by showing that an infinite-order generator can be constructed to exhibit smaller spread for arbitrary times. We propose a natural time scale for the construction of these higher-order generators and discuss results for matrices sampled from Gaussian Unitary Ensembles, demonstrating smaller Krylov complexity at all higher orders. These results extend the framework of Krylov complexity beyond the conventional Krylov basis by disproving the widely held assumption of optimality, extending the construction to higher-order generators, and introducing a physically motivated method for their construction. Our findings therefore suggest that previous statements and results concerning Krylov complexity may need to be reconsidered.

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

\sqrt{D}

.

A bstract We investigate the Krylov complexity of Schrödinger field theories, focusing on both bosonic and fermionic systems within the grand canonical ensemble which includes a chemical potential. Krylov complexity measures operator growth in quantum systems by analyzing how operators spread within the Krylov space, a subspace of the Hilbert space spanned by successive applications of the superoperator [ H , ·] on an initial operator. Using the Lanczos algorithm, we construct an orthonormal Krylov basis and derive the Lanczos coefficients, which govern the operator connectivity and thus characterize the complexity. Our study reveals that the Lanczos coefficients { b n } are almost independent of the chemical potential, while { a n } are dependent on the chemical potential. Both { a n } and { b n } show linear relationships with respect to n . For both bosonic and fermionic systems, the Krylov complexities behave similarly over time, especially at late times, due to the analogous profiles of the squared absolute values of their autocorrelation functions | φ 0 ( t )| 2 . The Krylov complexity grows exponentially with time, but its asymptotic scaling factor λ K is significantly smaller than the twice of the slope of the { b n } coefficients, contrasting to the relativistic field theories where the scaling aligns more closely with the twice of the slope of { b n }.

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

t^{1/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 (

t^{1/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.