PhysRevE 104 054105
PhysRevE 104 054105
PhysRevE 104 054105
We analyze a one-dimensional XXZ spin chain in a disordered magnetic field. As the main probes of the
system’s behavior, we use the sensitivity of eigenstates to adiabatic transformations, as expressed through the
fidelity susceptibility, in conjunction with the low-frequency asymptotes of the spectral function. We identify a
region of maximal chaos—with exponentially enhanced susceptibility—which separates the many-body local-
ized phase from the diffusive ergodic phase. This regime is characterized by slow transport, and we argue that
the presence of such slow dynamics highly constrains any possible localization transition in the thermodynamic
limit. Rather, the results are more consistent with absence of the localized phase.
DOI: 10.1103/PhysRevE.104.054105
and ergodic regimes in finite systems and prevents the system the anisotropy parameter at = 1.1 such that the model is
from localization in the thermodynamic limit. close to the extensively studied Heisenberg spin chain [24,42]
In the last decade it was realized that fidelity susceptibility and at the same time has broken SU (2) symmetry, even in the
χ , or more generally, the quantum geometric tensor, is a very absence of disorder. It is expected that these minor modifi-
efficient measure for detecting zero-temperature quantum cations of the model do not affect any results related to the
phase transitions (see, for example, Refs. [34–37]). Fidelity MBL transition apart from a small shift of the critical disorder
susceptibility defines the sensitivity of the (ground) state to strength.
small perturbations. Near the phase transition this sensitivity
is usually amplified, leading to singular behavior of χ ; often a II. FIDELITY SUSCEPTIBILITY
divergence. Physically χ is determined by the low-frequency
tail of the spectral function, which is enhanced near phase The fidelity susceptibility, or equivalently, the diagonal
transition boundaries due to the critical slowing down. At component of the quantum geometric tensor, with respect to
finite energy densities the utility of the fidelity susceptibility some coupling λ of a given eigenstate n is defined as [34,37]
was less obvious because of its exponential divergence with ←− |n|∂λ H|m|2
the system size in ergodic systems [37]. We note that there are χn = n| ∂λ ∂λ |nc ≡ . (2)
finite-temperature generalizations of the fidelity susceptibility m=n
(En − Em )2
measuring distance between density matrices like the Bures
The Hilbert-Schmidt norm of the AGP is defined as an av-
metric, which in turn is related to the Fisher information (see,
erage over the eigenstates of χn . In order to avoid large
e.g., Refs. [38,39]). However, the latter are not sensitive to
eigenstate-to-eigenstate fluctuations but keeping the exponen-
chaos or integrability, but rather probing equilibrium thermo-
tial sensitivity of this probe, it was proposed in Ref. [32] to
dynamic properties of the system.
additionally regularize χn by introducing an energy cutoff μ
In recent work we have shown that the norm of the adia-
that is exponentially small in the system size.
batic gauge potential (AGP), which is equivalent to the fidelity
As a probe λ we will use the local longitudinal magnetic
susceptibility χn averaged over different eigenstates n, can
field acting on a single spin, i.e.,
serve as a very sensitive probe of quantum chaos [32]. The
AGP norm is able to pick up tiny (exponentially small in H → H + λSlz , ∂λ H = Slz . (3)
the system size) integrability breaking perturbations, which
are not necessarily visible to traditional measures of quantum All calculations will be done at λ = 0, i.e., we will analyze
chaos like the level repulsion [40], the spectral form factor the sensitivity of eigenstates to an infinitesimal increase of a Z
[24], or the survival probability [41]. This sensitivity comes magnetic field on a single site. By direct inspection it becomes
from the fact that the low-frequency tail of the spectral func- clear that the disorder average fidelity susceptibility exponen-
tion defining χn can detect changes at very (exponentially) tially diverges with the system size as long as W < ∞. At
long timescales, where small perturbations should leave the large W , this divergence comes from rare resonances where
most pronounced effect. In particular, dependence of the AGP |Em − En | 2−N , which dominate the average χn in the ab-
norm on the system size changes from polynomial for inte- sence of level repulsion. In contrast, since the resonances are
grable systems to exponential for chaotic systems, and this rare, the typical fidelity susceptibility does not diverge with
change is very easy to detect numerically. One of the key find- the system size at sufficiently large W . As such, it should
ings of that work was that the transition from the integrable become small, since eigenstates are almost polarized along
to ergodic phase across various models happens through an the Z axis and adding an extra magnetic field does not affect
intermediate phase, where the fidelity susceptibility diverges them. For this reason, it is much more convenient to analyze
even faster than in the ergodic regime, and this divergence is the scaling behavior of the typical χn obtained by averaging its
accompanied by exponentially long in the system size relax- logarithm. We will define the typical log-fidelity susceptibility
ation times. In another recent work [33], based on a similar as
analysis, it was further argued that in a particular disordered ζ = log(χn ) , (4)
central spin model this chaotic behavior coexists with the non-
thermalizing nature of the individual eigenstates, and hence where . . . stands for averaging over both different dis-
to nonergodic behavior of the system even in the absence of order realizations and different eigenstates.. We note that in
small parameters. the Eigenstate thermalization hypothesis (ETH) regime, ζ is
In this work we adopt the approach of Ref. [32] to analyze equivalent to the logarithm of the AGP norm analyzed in
properties of a disordered one-dimensional XXZ spin chain Ref. [32], since the susceptibility is concentrated around the
and specifically analyze the fate of the many-body localization mean [see, e.g., Fig. 6(A)].
transition and the nature of the nonergodic regime there. We It is a simple exercise to extract asymptotes of ζ at weak
use the standard Hamiltonian and strong disorder. At small W the system is expected to be
ergodic and fully obey ETH [43]. Then (see, e.g., Ref. [32])
1 x x one expects
H= S j S j+1 + S yj S yj+1 + S zj S zj+1 + h j S zj , (1)
W j j ζ = L log 2 + A. (5)
where S x,y,z
j are the spin-1/2 operators and h j are uncorre- The constant A generally depends on the couplings in the
lated random numbers uniformly distributed in the interval Hamiltonian but is insensitive to the system size (up to possi-
[−1, 1]. We assume periodic boundary conditions, and we fix ble logarithmic corrections). In the opposite regime of strong
054105-2
DYNAMICAL OBSTRUCTION TO LOCALIZATION IN A … PHYSICAL REVIEW E 104, 054105 (2021)
054105-3
DRIES SELS AND ANATOLI POLKOVNIKOV PHYSICAL REVIEW E 104, 054105 (2021)
(a)
054105-4
DYNAMICAL OBSTRUCTION TO LOCALIZATION IN A … PHYSICAL REVIEW E 104, 054105 (2021)
054105-5
DRIES SELS AND ANATOLI POLKOVNIKOV PHYSICAL REVIEW E 104, 054105 (2021)
been numerically extracted from the spectral function (up to where a > 0 is some constant and the constants C and C
L = 20) by Serbyn et al. [50]. Whereas the data is consistent are related by the continuity of the spectral function. For
with ours, it is analyzed differently and the authors had a example, such a spectral function with a = 3 describes the
different interpretation, supporting the Griffiths phase picture zero-temperature disordered model analyzed in Ref. [57] and
on the ergodic side, leading up to the critical fan where the with a = 3/2 describes an infinite temperature-disordered
Thouless energy decays exponentially with the system size spin model on a random regular graph [58]. Then as z di-
ωth ∝ exp(−κL). verges near the critical point, the sum rule is satisfied as
Let us end this section by noting that the commonly ac- long as a > 1. While this behavior of the spectral function is
cepted Griffiths picture, with a critical point at finite disorder, possible, it implies that the dynamical exponent near the tran-
is not internally consistent. The argument requires only three sition, z(ω) ≡ 1/(1 + d log | f (ω)|2 /d log |ω|) = log(1/ω)/a,
inputs: locality, conservation of energy, and bounded norm increases as ω becomes smaller. Physically this behavior of
of local operators. Consider a Griffiths scenario [52] where z(ω) implies that already on the ergodic side, the system’s
local correlations, of a conserved operator like magnetization dynamics is getting slower and slower with time instead of
or simply an operator coupled to the energy, decay in time accelerating as the system “realizes” that it is in the ergodic
like C(t ) ∼ t −1/z , where the dynamical exponent z depends on regime. Also, such scenarios are inconsistent with existing nu-
disorder W . In the Griffiths picture, z would be proportional merical results, both reported here in Fig. 2 and in other works
to the correlation length z ∼ ξ (W ). Since energy is conserved, [53,54]. The only remaining possibility to satisfy the spectral
decaying correlations imply spreading of this magnetization sum rule with the system-size-independent C(W ) at z → ∞ is
in space like x ∼ t 1/z , and the spreading stops at the Thou- when C(W ) → 0 at some critical disorder Wc . This scenario
less time when the dynamics detects the boundaries, i.e., the corresponds to opening a spectral gap in the system. Such
Thouless energy obeys ωth ∼ L −z . In order for ergodicity to be a spectral gap is present in normal Bethe-ansatz integrable
broken, one needs the Thouless energy to reach the level spac- models [32], suppressing the exponential divergence of the
ing, which implies z L/ log L. In the Griffiths picture, this fidelity susceptibility with the system size. However, this gap
condition is equivalent to saying that the correlation length usually exists only for special classes of perturbations, which
ξ (W ) > L has to exceed the system size, and as such the keep the system integrable. There is no numerical indication
system is critical. In frequency space the spectral function for that a spectral gap opens up in disordered systems.
general subdiffusion with dynamical exponent z simply reads
V. FIDELITY SUSCEPTIBILITY REVISITED
C
| f (ω)|2 = , For the sake of completeness we now revisit the fidelity
ω 1−1/z
susceptibility, as there are a few more interesting points to
where C = | f (ω = 1)|2 is defined by the UV spectral weight. make. From Fig. 1 we clearly discern three different scaling
Note that, thus far, this scenario would also yield exponential regimes of the typical fidelity susceptibility denoted by three
enhancement of the susceptibility in the critical region. How- different colors. These regimes can be loosely labeled as (i)
ever, it suffers from the same problem as before, namely, the completely ergodic or ETH type, (ii) glassy with slow relax-
total amount of spectral weight must be conserved. The latter ation, and (iii) localized. In this section we will analyze these
directly implies three regimes in more detail.
First we analyze the total conserved part of the local Z
Cz 1. magnetization, which is directly related to the integral of the
spectral function [see Eq. (9)]. We define the latter as
This inequality can be satisfied for any finite z, but for the 2
system to become critical, z has to diverge with L, implying Z = 4 exp log n Slz n . (12)
C has to decay to zero at least like log(L)/L. Not only is this We choose the normalization such that Z = 1 in the fully
inconsistent with the numerical data (see Fig. 3), it is simply localized regime and Z ∼ 2−L in the ergodic regime in the
impossible by locality. Recall that C is the UV spectral weight zero-magnetization sector. To avoid effects of the tail of the
at the onset of the asymptotic hydrodynamic scaling behavior. distribution of magnetization, we first take its logarithm, then
As such it cannot possibly know about the system size L. Of average it over different eigenstates and disorder realizations,
course, the spectral weight can be a function of disorder W , and then exponentiate. In Fig. 4 we show 1 − Z, measuring
and while we argued that C ∼ 1/W such that W ∗ ∼ L, faster the nonconserved part of the local magnetization. It is evident
decay will simply result in a slower drift of the critical point. that as disorder increases, the magnetization starts to localize.
This argument is consistent with the observed existence As we increase the system size, the transition to the localized
of the localization transition in the nonlocal models like the regime happens for larger and larger values of disorder, and
random regular graphs [55,56]. Note that in principle one can the observed drift is consistent with the drift of the peak in the
imagine other scenarios, where the exponent z for example fidelity susceptibility (see Fig. 1). At very large W all curves
would be frequency dependent, approaching unity at the crit- collapse into one, and thus the conserved magnetization be-
ical point only in the limit ω → 0. Such a scenario could comes insensitive to the system size. The black line is shown
correspond to logarithmic corrections to the 1/ω scaling: for reference, representing the leading perturbative correction
to the infinite disorder limit.
C
ω loga (1/ω)
ω > e−az Next in Fig. 5 we show the mean scaled fidelity suscep-
| f (ω)|2 = C
ω1−1/z
0 < ω < e−az , tibility χ̄ /2L ≡ χn /2L and the typical one exp[ζ ] versus
054105-6
DYNAMICAL OBSTRUCTION TO LOCALIZATION IN A … PHYSICAL REVIEW E 104, 054105 (2021)
(a)
(b)
disorder. The plot of the latter (solid lines black with circles)
reproduce the data from Fig. 1, except disorder is now shown
in the linear scale. In the ergodic regime at low disorder,
the mean and the typical susceptibilities are parallel to each
other with only a constant offset between them, which comes
from small Gaussian fluctuations of χn around the mean.
As disorder increases the typical susceptibility reaches the
FIG. 6. AGP distribution. The distribution of the fidelity suscep-
maximum and goes down, while χ̄ keeps growing until it
tibility is shown for a chain of L = 16. To make the structure more
saturates at large W . As we discussed earlier, the saturated
apparent we have divided disorder range in three intervals. (a) Disor-
value is determined by rare resonances in Eq. (2), where der increasing from W = 0.5 (blue) to W = 1.75 (red). (b) Disorder
|En − En±1 | 2−L . Such resonances inevitably occur due to increasing from W = 2.1 (blue) to W = 4.3 (red), and (c) shows the
absence of level repulsion (see also Fig. 7). As a result, the disorder increasing from W = 5 (blue) to W = 15 (red). The disorder
mean fidelity susceptibility fluctuates much more at strong is logarithmically equally spaced. Various exponentials are shown by
disorder than the typical one. dashed lines for comparison. The different panels are also labeled on
One can extract additional information about the system Fig. 7.
by analyzing the tail of the probability distribution of χn , as
it was first discussed in Ref. [32]. Both in ETH and localized
regimes, this tail is determined by anomalously small nearest c exp[−z/(2 + α)], we find
z z
energy differences |En − En+1 | and thus contains additional
P(z) ∝ exp − Q c exp − ,
information about the level statistics. For random matrix 2+α 2+α
ensembles, the distribution has been derived in [59] and a where Q(s) is the probability to have the level spacing s. At
numerical survey on disordered systems has been presented in small values of s, corresponding to large z, this probability
Ref. [60]. To avoid exponential factors, it is more convenient usually takes a universal power law form: Q(s) ∝ sβ with
to deal with the distribution of ζn ≡ log(χn ). The probability β = 0, indicating Poisson statistics, β > 0 corresponding to
that ζn = z at large values of z can be estimated as the level repulsion and β < 0 corresponding to the level at-
P(z) ≡ Pr(ζn = z) ≈ Pr(log(| f 2 (s)|/s2 ) = z), traction. Combining all these factors together one finds
where s is the smallest of |En − En+1 | and |En − En−1 |. 1+β
P(z) ∝ exp −z . (13)
Assuming that | f (s)|2 = c|s|−α at small s such that s = 2+α
054105-7
DRIES SELS AND ANATOLI POLKOVNIKOV PHYSICAL REVIEW E 104, 054105 (2021)
054105-8
DYNAMICAL OBSTRUCTION TO LOCALIZATION IN A … PHYSICAL REVIEW E 104, 054105 (2021)
054105-9
DRIES SELS AND ANATOLI POLKOVNIKOV PHYSICAL REVIEW E 104, 054105 (2021)
care in computing the typical susceptibility ζn = log χn , as size system L from our previous argument, as it becomes more
this would diverge. For those states we need to consider higher and more unlikely to have have an all polarized region since
order corrections in perturbation theory, since the first order the total magnetization ultimately needs to sum up to zero.
contribution vanishes. Fortunately that is straightforward to Consider the ith step, in which we are looking at the faith of
do, as long as there are no resonances such that we can neglect a polarized region consisting of 2i − 1 spins. The probability
renormalization of the denominators. To first order we had that both the left and right neighbors have the same polariza-
tion is
ζ = 3
4
log W −2 + 41 log O(W −4 ) + B, (A5)
[L/2 − (2i − 1)] (L/2 − 2i)
where B is some universal constant that is independent of pi = . (A7)
L − (2i − 1) L − 2i
L, W and comes from the disorder average over h. In order
for the polarized states to contribute in the next order, they One can go through exactly the same argument as before,
need to move an antialigned spin to the center by an extra resulting in a generalization of expression (A6):
application of the flip-flip Hamiltonian HX X . Consequently,
either the left or right must have opposite spin to the polarized
L/4 i−1
[1] G. P. Berman and F. M. Izrailev, Chaos: An Interdisciplinary [12] L. D’Alessio and A. Polkovnikov, Ann. Phys. 333, 19 (2013).
Journal of Nonlinear Science 15, 015104 (2005). [13] A. Haldar, R. Moessner, and A. Das, Phys. Rev. B 97, 245122
[2] M. A. Porter, N. J. Zabusky, B. Hu, and D. K. Campbell, Am. (2018).
Sci. 97, 214 (2009). [14] M. Heyl, P. Hauke, and P. Zoller, Sci. Adv. 5 (2019).
[3] V. Arnol’d, Mathematical Methods of Classical Mechanics [15] I. V. Gornyi, A. D. Mirlin, and D. G. Polyakov, Phys. Rev. Lett.
(Springer-Verlag, New York, 1989). 95, 206603 (2005).
[4] D. A. Abanin, W. De Roeck, and F. Huveneers, Phys. Rev. Lett. [16] D. Basko, I. Aleiner, and B. Altshuler, Ann. Phys. 321, 1126
115, 256803 (2015). (2006).
[5] T. Mori, T. Kuwahara, and K. Saito, Phys. Rev. Lett. 116, [17] V. Oganesyan and D. A. Huse, Phys. Rev. B 75, 155111 (2007).
120401 (2016). [18] J. Z. Imbrie, J. Stat. Phys. 163, 988 (2016).
[6] P. Weinberg, M. Bukov, L. D’Alessio, A. Polkovnikov, S. [19] R. Nandkishore and D. A. Huse, Annu. Rev. Condens. Matter
Vajna, and M. Kolodrubetz, Phys. Rep. 688, 1 (2017). Phys. 6, 15 (2015).
[7] A. Rajak, R. Citro, and E. G. D. Torre, J. Phys. A: Math. Theor. [20] D. A. Abanin, E. Altman, I. Bloch, and M. Serbyn, Rev. Mod.
51, 465001 (2018). Phys. 91, 021001 (2019).
[8] O. Howell, P. Weinberg, D. Sels, A. Polkovnikov, and M. [21] M. Žnidarič and M. Ljubotina, Proc. Natl. Acad. Sci. 115, 4595
Bukov, Phys. Rev. Lett. 122, 010602 (2019). (2018).
[9] S. Vajna, K. Klobas, T. Prosen, and A. Polkovnikov, Phys. Rev. [22] E. V. H. Doggen, I. V. Gornyi, A. D. Mirlin, and D. G. Polyakov,
Lett. 120, 200607 (2018). arXiv:2101.05651.
[10] W. D. Roeck and V. Verreet, arXiv:1911.01998. [23] P. Sierant, D. Delande, and J. Zakrzewski, Phys. Rev. Lett. 124,
[11] T. Prosen, Phys. Rev. Lett. 80, 1808 (1998). 186601 (2020).
054105-10
DYNAMICAL OBSTRUCTION TO LOCALIZATION IN A … PHYSICAL REVIEW E 104, 054105 (2021)
[24] J. Šuntajs, J. Bonča, T. Prosen, and L. Vidmar, Phys. Rev. E [48] R. Nandkishore, S. Gopalakrishnan, and D. A. Huse, Phys. Rev.
102, 062144 (2020). B 90, 064203 (2014).
[25] J. Šuntajs, J. Bonča, T. Prosen, and L. Vidmar, Phys. Rev. B [49] O. S. Barišić, J. Kokalj, I. Balog, and P. Prelovšek, Phys. Rev.
102, 064207 (2020). B 94, 045126 (2016).
[26] D. Abanin, J. H. Bardarson, G. D. Tomasi, S. Gopalakrishnan, [50] M. Serbyn, Z. Papić, and D. A. Abanin, Phys. Rev. B 96,
V. Khemani, S. A. Parameswaran, F. Pollmann, A. C. Potter, M. 104201 (2017).
Serbyn, and R. Vasseur, Ann. Phys. 427, 168415 (2021). [51] The Thouless energy is extracted numerically by identifying the
[27] R. K. Panda, A. Scardicchio, M. Schulz, S. R. Taylor, and M. point of maximal curvature of the data on log-log scale. This
Žnidarič, Europhys. Lett. 128, 67003 (2020). procedure is robust up to W ≈ 1.5. The coefficient C, for the
[28] F. H. L. Essler, R. van den Berg, and V. Gritsev, Phys. Rev. B asymptotic C/W ω behavior of the spectral function, is extracted
98, 024203 (2018). by multiplying the spectral function by W ω and extracting the
[29] M. Kiefer-Emmanouilidis, R. Unanyan, M. Fleischhauer, and J. value of the inflection point. This estimate smoothly goes from
Sirker, Phys. Rev. Lett. 124, 243601 (2020). 0.0245 to 0.0179 when disorder changes from W = 1 to W =
[30] M. Kiefer-Emmanouilidis, R. Unanyan, M. Fleischhauer, and 2.5.
J. Sirker, Absence of true localization in many-body localized [52] K. Agarwal, E. Altman, E. Demler, S. Gopalakrishnan, D. A.
phases, Phys. Rev. B 103, 024203 (2021). Huse, and M. Knap, Ann. Phys. 529, 1600326 (2017).
[31] D. J. Luitz and Y. B. Lev, Phys. Rev. B 102, 100202(R) (2020). [53] S. Bera, G. De Tomasi, F. Weiner, and F. Evers, Phys. Rev. Lett.
[32] M. Pandey, P. W. Claeys, D. K. Campbell, A. Polkovnikov, and 118, 196801 (2017).
D. Sels, Phys. Rev. X 10, 041017 (2020). [54] G. De Tomasi, D. Hetterich, P. Sala, and F. Pollmann, Phys. Rev.
[33] T. Villazon, P. W. Claeys, M. Pandey, A. Polkovnikov, and A. B 100, 214313 (2019).
Chandran, Sci. Rep. 10, 16080 (2020). [55] V. E. Kravtsov, I. M. Khaymovich, E. Cuevas, and M. Amini,
[34] L. Campos Venuti and P. Zanardi, Phys. Rev. Lett. 99, 095701 New J. Phys. 17, 122002 (2015).
(2007). [56] V. Kravtsov, B. Altshuler, and L. Ioffe, Ann. Phys. 389, 148
[35] A. F. Albuquerque, F. Alet, C. Sire, and S. Capponi, Phys. Rev. (2018).
B 81, 064418 (2010). [57] O. Motrunich, K. Damle, and D. A. Huse, Phys. Rev. B 63,
[36] M. Kolodrubetz, V. Gritsev, and A. Polkovnikov, Phys. Rev. B 134424 (2001).
88, 064304 (2013). [58] K. S. Tikhonov and A. D. Mirlin, Phys. Rev. B 99, 024202
[37] M. Kolodrubetz, D. Sels, P. Mehta, and A. Polkovnikov, Phys. (2019).
Rep. 697, 1 (2017). [59] P. Sierant, A. Maksymov, M. Kuś, and J. Zakrzewski, Phys.
[38] J. S. Sidhu and P. Kok, AVS Quantum Sci. 2, 014701 Rev. E 99, 050102(R) (2019).
(2020). [60] A. Maksymov, P. Sierant, and J. Zakrzewski, Phys. Rev. B 99,
[39] D. Liska and V. Gritsev, SciPost Phys. 10, 020 (2021). 224202 (2019).
[40] L. D’Alessio, Y. Kafri, A. Polkovnikov, and M. Rigol, Adv. [61] R. Vosk, D. A. Huse, and E. Altman, Phys. Rev. X 5, 031032
Phys. 65, 239 (2016). (2015).
[41] F. Borgonovi, F. M. Izrailev, L. F. Santos, and V. G. Zelevinsky, [62] A. C. Potter, R. Vasseur, and S. A. Parameswaran, Phys. Rev. X
Phys. Rep. 626, 1 (2016). 5, 031033 (2015).
[42] M. Serbyn, Z. Papić, and D. A. Abanin, Phys. Rev. X 5, 041047 [63] T. Thiery, M. Müller, and W. D. Roeck, arXiv:1711.09880.
(2015). [64] A. Rubio-Abadal, J.-y. Choi, J. Zeiher, S. Hollerith, J. Rui, I.
[43] We exclude very small values of disorder W 0.1, where the Bloch, and C. Gross, Phys. Rev. X 9, 041014 (2019).
system becomes very close to the integrable XXZ model. [65] J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, P. Lugan,
[44] K. Agarwal, S. Gopalakrishnan, M. Knap, M. Müller, and E. D. Clément, L. Sanchez-Palencia, P. Bouyer, and A. Aspect,
Demler, Phys. Rev. Lett. 114, 160401 (2015). Nature (London) 453, 891 (2008).
[45] I. Khait, S. Gazit, N. Y. Yao, and A. Auerbach, Phys. Rev. B 93, [66] D. Sels and E. Demler, Quantum generative model for sam-
224205 (2016). pling many-body spectral functions, Phys. Rev. B 103, 014301
[46] D. J. Luitz and Y. B. Lev, Ann. Phys. 529, 1600350 (2017). (2021).
[47] M. Žnidarič, A. Scardicchio, and V. K. Varma, Phys. Rev. Lett. [67] D. V. Vasilyev, A. Grankin, M. A. Baranov, L. M. Sieberer, and
117, 040601 (2016). P. Zoller, PRX Quantum 1, 020302 (2020).
054105-11