Universal Hypothesis of Autocorrelation Function from Krylov Complexity

In a quantum many-body system, autocorrelation functions can determine linear responses nearby equilibrium and quantum dynamics far from equilibrium. In this letter, we bring out the connection between the operator complexity and the autocorrelation function. In particular, we focus on a particular kind of operator complexity called the Krylov complexity. We find that a set of Lanczos coefficients $\{b_n\}$ computed for determining the Krylov complexity can reveal the universal behaviors of autocorrelations, which are otherwise impossible. When the time axis is scaled by $b_1$, different autocorrelation functions obey a universal function form at short time. We further propose a characteristic parameter deduced from $\{b_n\}$ that can largely determine the behavior of autocorrelations at the intermediate time. This parameter can also largely determine whether the autocorrelation function oscillates or monotonically decays in time. We present numerical evidences and physical intuitions to support these universal hypotheses of autocorrelations. We emphasize that these universal behaviors are held across different operators and different physical systems.

The autocorrelation function refers to the temporal correlation of the same operators Ô at two different times, denoted by C(t) = Ô(t) Ô(0) .It plays an important role in studying quantum matters, including condensed matter materials, ultracold atomic gases, and NMR systems.The autocorrelation functions can be measured by spectroscopy methods, which reveals linear response nearby equilibrium.It can also be measured in far-fromequilibrium dynamics, for example, through quench experiments.In non-equilibrium situations, autocorrelation functions can determine quantum dynamics involving highly excited states.
The behavior of the autocorrelation function C(t) crucially depends on the time evolution of Ô(t), which follows the Heisenberg evolution Ô(t) = e i Ĥt Ôe −i Ĥt , with Ĥ being the Hamiltonian of the physical system.Usually, Ô(t) becomes more and more complicated as time evolves, consequently reducing the autocorrelation C(t).This observation brings out the connection between operator complexity and its autocorrelation function.In the past years, various measures have been proposed to describe operator complexity quantitatively and to study how the complexity of an operator grows under the Heisenberg evolution [1][2][3][4][5][6][7][8][9][10][11].In particular, the Krylov complexity has been proposed to quantify operator complexity recently [11].The advantage of the Krylov complexity is that it exhibits universal behavior for generic chaotic quantum many-body systems, and such universal behaviors are shared by a large class of different quantum many-body Hamiltonians .
Therefore, it is natural to ask whether these recent developments in measuring operator complexity can help us better understand the autocorrelation function's universal behavior.Especially, since the Krylov complexity exhibits universality, the question is, by utilizing the notation introduced for studying the Krylov complexity, whether we can reveal hidden universal behaviors of the autocorrelation function, which are otherwise impossible.
Let us first introduce the notation for describing the Krylov complexity [11].Given a Hilbert space spanned by {|i }, an operator Ô = ij O ij |i j| can be mapped to a state Ô in the double space as Ô = ij O ij |i ⊗ |j , and vice versa.By using the Baker-Campbell-Hausdorff formula, the Heisenberg evolution can be expressed as where L Ô ≡ [ Ĥ, Ô].Eq. ( 1) can be viewed as expanding a state Ô under a set of states { Ln Ô }.However, this set of states are neither orthogonal nor normalized.Applying the Gram-Schmidt procedure to this set of states yield a set of orthogonal basis, denoted by the Krylov basis Ŵn [11].Here the inner product is introduced as In general, the larger n, the more complicated the operator Ŵn .Expanding Ô(t) under the Krylov basis gives It has been shown that ϕ n (t) obeys the following differential equation The Krylov complexity is defined as Here {b n } are called the Lanczos coefficients introduced in the Gram-Schmidt procedure.They depend on both the choice of operator Ô and the system Hamiltonian Ĥ.The units of {b n } are energy.Note that Ŵ0 = 1 b0 Ô with a normalization factor b 2 0 = Ô Ô , arXiv:2305.02356v1 [cond-mat.stat-mech]3 May 2023 we have ϕ 0 = Ŵ0 Ô = 1 b0 C(t).That is to say, the time evolution of ϕ 0 gives rise to the autocorrelation function, normalized by its value at t = 0.
Eq. ( 3) is reminiscent of the Schrödinger equation for a single particle hopping along a half-infinite chain.At t = 0, only ϕ 0 = 1 and all other ϕ n =0 = 0.As time evolves, this particle hops away from n = 0, consequently reducing the autocorrelation C(t) and increasing the Krylov complexity.Inspired by this connection, it is realized that the dynamics of C(t) is governed by the Lanczos coefficients {b n }.Hence, if there exist certain universal behaviors of C(t), it is conceivable that their information is hidden inside {b n }.
Summary of the Main Results.Before presenting the details, we first summary the main findings of this work.First of all, we note that originally the Lanczos coefficients b n are only defined for all non-negative integers.However, we can use these data, especially the data with small n, to interpolate a smooth function b[x] defined for all x > 0. We require that the interpolated function has to be smooth enough and cannot strongly vary between two neighboring integers.1. Short-Time Universality: When the time t is scaled by b [1], the autocorrelation function C as a function of tb [1] exhibits universal behavior for tb [1] < ∼ 1 across all different choices of operators and Hamiltonians.
We have run extensive numerical tests which support these hypothesis.Below, we first present a set of representative numerical evidences, and the complete code to verify our hypothesis are also available [41].
Here we first present a Gedanken numerical simulation.This simulation does not involve physical Hamiltonians and only solves Eq. in terms of tα, C(t) deviates from each other once tα is non-zero.However, when C(t) is plotted in terms of tb 1 , all C(tb 1 ) collapse when tb 1 < ∼ 1.This difference can be seen clearly by comparing Fig. 1(b) with its inset.This demonstrates the first hypothesis.
However, these curves shown in Fig. 1 [1] for different choices of {b n }.We find that when b [1]/b[1] <∼ 0.5, the corresponding C(t) all exhibit oscillatory behavior.This can also be seen from all cases in Fig. 1(c) and (d).When b [1]/b[1] >∼ 0.5, the corresponding C(t) all monotonically decay, as one can also seen from all cases in Fig. 1(e) and (f).This demonstrates the third hypothesis.All together, the Gedanken simulation shows that the autocorrelation function exhibits universal behavior in terms of tb 1 and a single parameter b [1]/b [1] can largely determine the behavior of C(t).
Physical Model.The Gedanken numerical simulation is inspiring and generic because it does not depend on the concrete physical operator and physical Hamiltonian.However, it also has limitations because the Lanczos coefficients {b n } obtained from physical models are usually not smooth enough.To show how much the discussion above can hold for realistic models, we show results calculated with two typical physical models below.
The first model is the one-dimensional quantum Ising model with both transversal and longitudinal fields, whose Hamiltonian is written as h/J and g/J are two tunable parameters in this model.The second model is the one-dimensional spinless fermion Hubbard model, whose Hamiltonian is written as Here J and J denote the nearest and the next nearest neighbor hopping strengths.V 1 and V 2 are the nearest and the next nearest interaction strengths.In both models, we use J as the natural energy unit although they represent different energy scales in two models.In Fig. 2, we first compare C(t) of two different operators in the quantum Ising model, but with different model parameters, as shown in Fig. 2 [1], with two from the quantum Ising model and two from the spinless Hubbard model.Their corresponding C(t) all exhibit the oscillatory behavior.On contrast, Fig. 3(c) show four cases from these two models with larger b [1]/b [1].Their corresponding C(t) all monotonically decay.Moreover, we note that because all these cases have different values of b [1]/b [1], therefore their C(t) do not collapse when tb 1 > ∼ 1 but they all collapse when tb 1 < ∼ 1.Hence, the information presented in Fig. 2 and Fig. 3 demonstrate the three hypothesis in realistic models.
Intuitions.Finally we discuss the intuitions that lead to these hypothesis.First, since all ϕ n = 0 at t = 0 except for ϕ 0 , and ϕ 0 only couples to ϕ 1 , the short-time dynamics of ϕ 0 is dominated by its coupling to ϕ 1 .Hence, 2) for red circles and = (1, 0.5) blue squares, and (J /J, V1/J, V2/J) = (0.2, 3, 0.5) for black diamonds and = (0.21, 1, 1) for green triangles.The corresponding C(t) is plotted in terms of tb1 in (b).(c) bn for the operators σz 0 (red circles and blue squares) in the quantum Ising model, and the operator 2n0 − 1 for the spinless Hubbard model (black diamonds and green triangles).The model parameters are (h/J, g/J) = (1, 0) for red circles and = (1, 0.5) blue squares, and (J /J, V1/J, V2/J) = (0.2, 5, 0.5) for black diamonds and = (0.21, 1, 3.5) for green triangles.The corresponding C(t) is plotted in terms of tb1 in (d).Other information about the system size, the solid fitting lines and the numbers in the legend are the same as described in the caption of Fig. 2.
by ignoring all ϕ n>1 , we can obtain This gives rise to a solution ϕ 0 = cos(tb 1 ).This shows that b 1 is a natural unit to scale t, resulting in a universal function form for autocorrelation functions at short-time.Secondly, we take the continuum limit of Eq. ( 3), which results in the following differential equation Now if we make a frame transformation to redefine t = tb[x], we arrive at the following equation where ∂ is defined as  3) can be solved exactly for δ = 0, 1/2 and 1, which respectively give B 1 (2αt)/(αt), δ = 0; exp −α 2 t 2 /2 , δ = 1/2; sech(αt), δ = 1, (10) where B 1 denotes the Bessel function of the first kind.For generic δ, there is no analytical solution but the equation can be easily solved numerically.It is found that nonmonotonic oscillation exists when δ < 0.5 but disappears when δ > 0.5.
Summary and Outlook.In summary, we have found three universal properties of the autocorrelation functions with the help of the Krylov complexity.To reveal these universal properties, two key findings are the characteristic parameter b [1]/b [1] and scaling time with b [1].We emphasize that these universal properties are shared by different operators in different systems.On the theory side, our results bring out the generic connection between complexity and correlation in quantum many-body systems that deserves further theoretical investigations.On the experimental side, the quench experiments recently performed in NMR and cold atom systems directly measure the autocorrelation function [42,43], and our results can be straightforwardly verified in these experiments.
Note Added.Oscillatory versus non-oscillatory behavior of auto-correlation function in random spin model has been discussed in Ref. [44].
Such an interpolation allows us to take derivatives of b[x] and the first-order derivative is denoted by b [x].With the help of b[x] and b [x], we present the following three universal hypothesis of C(t).

( 3 ) 1 FIG. 1 :
FIG. 1: (a),(c) and (e) show several different choices of bn generated by b[n] = αf (n), with different functions f (n) shown in the legend.In (a), b [1]/b[1] take five different values as shown by the numbers in the legend.b [1]/b[1] = 1/10 in (b) and b [1]/b[1] = 2/3 in (c) for all choices of different functions of b[x].(b), (d) and (f) show the autocorrelation function C(t) as a function of tb1, corresponding to bn shown in (a), (c) and (e).The insets in (b), (d) and (f) show C(t) as a function of tα.
(b) do not collapse when tb 1 > ∼ 1.To exhibit universality in the intermediate time, we require that different set of {b n } share a common characteristic parameter.A main finding of this work is this characteristic parameter, which turns out to be b [1]/b[1].In Fig. 1(c) and (e), we show two sets of {b n }.Within each set, the values of {b n } are very different.However, they share the same value of b [1]/b[1].It is easy to check that the five different sets of {b n } shown in Fig. 1(c) and Fig. 1(e) respectively share b [1]/b[1] = 1/10 and b [1]/b[1] = 2/3.We then calculate C(t) using {b n } shown in Fig. 1(c) and (e), which respectively lead to results shown in Fig. 1(d) and (f).The inset of Fig. 1(d) and (f) also show that these autocorrelations C(t) behave very different in term of tα.However, when they are plotted in terms of tb 1 , Fig. 1(d) and (f) show that they perfectly collapse into a single curve up to tb 1 ∼ 5.This demonstrates the second hypothesis.The numbers shown in the legend of Fig. 1(a) are b [1]/b (a) and (b).Then, we compare C(t) for two different operators in two different physical systems, one in the quantum Ising model and the other in the spinless Hubbard model.The results are shown in Fig. 2(c) and (d).The Lanczos coefficients {b n } are now calculated with the chosen operators and the physical Hamiltonian.Then we use polynomial function to fit {b n } with respect to a few smallest integer n, up to n ∼ 10.With the fitted function, we can obtain b [1]/b[1] with an error bar from the fitting error.The fitting error also includes fluctuation due to varying the range of polynomial function and varying the number of integer points included in the fitting.

1 (
a-b) The Lanczos coefficients bn and the autocorrelation functions C(t) for two different operators σz 0 σx blue squares) and σy 0 σy 0 (red circles) in the quantum Ising model.(h/J, g/J) = (1, 0) for blue squares and = (1, 1) for red circles.(c-d) bn and C(t) for two different operators in two different systems.Red circles are results for σx 0 in the quantum Ising model with (h/J, g/J) = (1, 1.2), and blue squares are results for 2n0 − 1 in the spinlees Hubbard model with (J /J, V1/J, V2/J) = (0.5, 2, 0.5).Here we calculate 13 sites with periodic boundary condition for the quantum Ising model and 14 sites for the spinless Hubbard model, and the lower index 0 and 1 represent site indices.In (b) and (d), C(t) is plotted in terms of tb1 and is plotted in terms of tJ in its inset.The sold lines in (a) and (c) are fitting curve of b[x], and the numbers in the legend of (a) and (c) show b [1]/b[1] for different cases, with error bars from fitting error.Each figure of Fig. 2(a) and (c) shows two cases with same value of b [1]/b[1] within the error bars.Their corresponding C(t) in terms of tb 1 are respectively shown in Fig. 2(b) and (d), compared with C(t) plotted in terms of tJ in the insets.It is clear that by changing tJ to tb 1 , the horizontal axes is stretched such that the two cases shown in each figure are in good agreement with each other up to tb 1 ∼ 2 − 5.

Fig. 3 (
Fig.3(a)shows four cases with small b [1]/b[1], with two from the quantum Ising model and two from the spinless Hubbard model.Their corresponding C(t) all exhibit the oscillatory behavior.On contrast, Fig.3(c)show four cases from these two models with larger b [1]/b[1].Their corresponding C(t) all monotonically decay.Moreover, we note that because all these cases have different values of b [1]/b[1], therefore their C(t) do not collapse when tb 1 > ∼ 1 but they all collapse when tb 1 < ∼ 1.Hence, the information presented in Fig.2and Fig.3demonstrate the three hypothesis in realistic models.

1 FIG. 3 :(red circles) and σy 0 σy 1 (
FIG. 3: (a) The Lanczos coefficients bn for two operators σx 0 ) Hence, in the new frame, the differential equation is solely controlled by b [x]/b[x].Since the autocorrelation function only concerns the time dynamics of ϕ at x = 0, we conjecture that b [1]/b[1] largely determines the behavior of the autocorrelation.Thirdly, as for why b [1]/b[1] ∼ 0.5 separates oscillating decay and monotonically decay, the intuition comes from solving simple situation with b n = αn δ .In these situations, b[1] = α and b [1]/b[1] = δ.It turns out that Eq. (