Quantum Zeno Effect in Noisy Integrable Quantum Circuits for Impurity Models
Yicheng Tang
[email protected]Department of Physics and Astronomy, Center for Materials Theory, Rutgers University,
Piscataway, New Jersey 08854, United States of America
Pradip Kattel
Department of Physics and Astronomy, Center for Materials Theory, Rutgers University,
Piscataway, New Jersey 08854, United States of America
J. H. Pixley
Department of Physics and Astronomy, Center for Materials Theory, Rutgers University,
Piscataway, New Jersey 08854, United States of America
Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, NY 10010
Natan Andrei
Department of Physics and Astronomy, Center for Materials Theory, Rutgers University,
Piscataway, New Jersey 08854, United States of America
Abstract
We theoretically study the open quantum system dynamics (in the Trotterized limit) of integrable quantum circuits in the presence of onsite dephasing noise with a spin-1/2 “impurity” interacting at the edge. Using a combination of Bethe Ansatz (BA) and exact diagonalization (ED), we study the dynamics of both the bulk and the impurity for the XXX (Heisenberg) and the XX qubit chains in the presence and absence of bulk noise. In the absence of noise, we show that the impurity exhibits two distinct phases, the bound mode phase where the impurity keeps oscillating in time, and the Kondo phase where it decays with Kondo time . Turning on the bulk dephasing noise, we find for the two models that in the long time limit in both regimes the quantum Zeno effect takes place where the dynamics of the impurity magnetization slows down as the noise strength increases. The impurity magnetization in the bound mode regime shows the opposite effect, decaying faster as the noise strength increases for short times (). We show that the bulk KPZ dynamics of the XXX model is converted to diffusive dynamics as in the XX case studied before by V. Alba, driving both systems to the Zeno effect slowing down the impurity dynamics in the long time limit.
Introduction:
The rise of noisy intermediate scale quantum (NISQ) [1, 2, 3] devices in a wide range of physical platforms including superconducting qubits [4, 5, 6], trapped ions [7, 8, 9, 10], neutral atoms [11, 12, 13], and circuit quantum electrodynamics [14, 15, 16, 17, 18], has made understanding the effects of noise on these devices and on quantum many-body systems (QMBS) in general a question of paramount importance [19, 20, 21, 22].
Simulating large generic QMBS on a classical computer is extremely difficult [23, 24, 25] as the size of the Hilbert space grows exponentially with the size of the system. To circumvent this challenge, quantum computers can be utilized to study QMBS [24, 25] (as they are built to hold such an exponentially large Hilbert space). This can be achieved by Trotterizing the continuous time evolution of a QMBS as a quantum circuit [26, 27].
While, generically solving a quantum circuit remains an equally challenging problem, integrable quantum circuits [28, 29, 30, 31, 32, 33], in particular, based on the -matrices that solve the Yang-Baxter equation [34, 35], can simulate the unitary evolution of initial states under integrable Hamiltonians [36, 37, 38, 39]. Integrability allows more control and insights into the behavior of these complex systems
garnering significant attention recently
[40, 41, 42].
All quantum devices currently accessible are under a significant influence of noise [43, 44, 45, 46, 47, 48, 49, 50, 51, 52].
Noise affects the dynamics and transport properties of quantum circuits resulting in dramatic changes [19, 20, 21]. Interesting effects arise such as dynamical phase transitions [53, 54], the quantum Zeno effect (where increasing the strength of noise slows down quantum dynamics [55, 56]) and the quantum Mpemba effect (where non-equilibrium quantum systems relax faster when they are further from equilibrium [57, 58, 59, 60]). Thus, studying the effect of noise on quantum circuits is of crucial importance to understand the experiments performed on NISQ devices [61].
Here, we study a circuit that can be solved exactly in the presence of dephasing noise.
Models: To begin with, we construct an integrable quantum circuit
that corresponds to a unitary operator
(1)
where is the number of qubits, given explicitly in [62] is the braid form of the trigonometric -matrix of the six-vertex model [63, 36] acting on the and qubits, is the boundary -matrix that solves the reflection algebra [64, 65, 66], and , and are free parameters.
In the Trotterziation limit [67, 68], choosing , this quantum circuit simulates the unitary evolution for an integrable Hamiltonian [25, 27, 69], where
(2)
with , and . This Hamiltonian describes the spin chain with a boundary impurity.
At the isotropic XXX ( and the XX () points, the model is integrable for arbitrary boundary coupling . For generic , the boundary coupling () and the boundary anisotropy parameter () are related and expressed in terms of a single free parameter to maintain the integrability [70, 71, 72, 73, 74, 75, 76]. In the ground state, the impurity is screened due to the Kondo effect when is purely imaginary or and it is screened by a bound mode exponenitially localized at the boundary when [70].
To study the effect of dephasing noise on the integrable quantum circuit, we solve the continuous-time dynamics of the noisy integrable Hamiltonian using the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation, [20, 19] that governs the dynamics of the density matrix with a Liouvillian [77, 78, 79, 80]
(3)
where the jump operators are , is the dephasing rate, is the Hamiltonian in Eq. (2), and
denote the anticommutator and the commutator, respectively. In the following, we will evolve the system from two initial states, Néel state and domain wall state expressed as:
where denote the set of sites with for the Néel (N) state , and for the domain wall (D) state .
Here, refers to the state with all spins down (qubits ) state. In both cases, we compute the time-dependent expectation value of magnetization at site , .
XX chain: Let us start from the XX point as it allows for a simpler treatment with similar phenomenology [71]. It is convenient to work in the fermionic language via a Jordan-Wigner transformation . The magnetization is then related to the two-point correlator as .
A closed set of equations gives the dynamics of the two-point correlators because the Hamiltonian and the jump operators are both quadratic in fermions [81]. Therefore, the dynamics of the two-point correlation function in the noisy XX model can be mapped to a two-particle quantum mechanical problem [62]. The correlator can be vectorized as a state (see [62]) where and forms an orthonormal basis in . Its evolution from an initial state reads with Liouvillian and
(4)
where .
The Liouvillian right eigenmodes
satisfy the stationary Schödinger equation with eigenvalue in terms of which we expand the time-dependent state . Here is the left eigenmode of the Liouvillian satisfying and it forms a bi-orthognal basis with as .
The time-dependent impurity magnetization (at site ) can be written as
(5)
where the amplitude depends on the initial state as follows .
Thus, requires the knowledge of the Liouvillian spectrum and of the amplitude which encodes the initial conditions.
Figure 1: (a) The parameter space is divided into two distinct regimes: the bound-mode regime, and the Kondo regime. The real part of the double-bound mode energy is computed for is given in the heatmap. For large , the boundary between two regimes is almost as shown in the dashed line in (a). The Liouvillian spectrum and amplitude for Néel state in the two-particle sector with dephasing rate and impurity coupling (b) in the Kondo regime, , (c) in the Kondo regime, and , (d) in the bound mode regime due to the existence of the double-bound mode. The data points labeled in the legend with different shapes and numbers represent four different kinds (respectively correspond to the four kinds of solutions solved analytically in the main text). The scattering color plots are numerical results solved with system size where the heatmap shows absolute values of amplitude for each eigenmode. The modes with larger dominate the impurity dynamics: the double-bound mode in the bound mode regime (d), and the modes within the diffusive band in the Kondo regime (b, c). With a small noise , the diffusive band is mixed with the other bands (b), and when the diffusive band is separated from the other bands (c,d).
Figure 2: The dynamics of the XX and XXX impurity models evolving from a Néel initial state under a Lindbladian evolution with and without dephasing noise. Different colors represent different dephasing rates , solid lines for in the Kondo regime, and dashed lines for in the bound-mode regime. The upper panel shows the data for the XXX chain for and the bottom panel for the XX chain for . For both XX and XXX chains, in the bound mode phase, the signal decays faster as the dephasing rate increases. However, in the Kondo phase, the impurity decay is slower as the dephasing rate increases.
The insets in each panel show the data on a log-linear scale displaying the same power law decay independent of in the long-time limit.
We now proceed to solve for the spectrum of the Liouvillian defined in Eq. (4).
In the absence of the impurity,
this Liouvillian can be mapped to an integrable non-Hermitian Hubbard model after a unitary transformation [81, 82, 83, 84, 85], here we generalize it to include the impurity as discussed in [62].
The Liouvillian in Eq. (4) can thus be solved via the Bethe Ansatz (BA) [86, 87, 88]. It has eigenvalue with the eigenmodes wavefunction taking the BA form
(6)
where , is the normalization factor such that and is the symmetrizer. The Liouvilian in Eq.(4) also has antisymmetrized eigenmodes [82], but only symmetrized eigenmodes participate in dynamics since the initial configuration is symmetric, i.e. . Here the amplitudes are related to the boundary scattering matrix and the bulk scattering matrix as: and [81, 82, 83, 84] with
(7)
The quantization condition for quasimomenta and are obtained from the Bethe Ansatz equations (BAE) [84]
(8)
If is a solution of the BAE (8), then is also a solution. This allows us to restrict the Brillouin zone to with () as the real (imaginary) part of . Furthermore, its complex conjugate also being a solution, indicates that if is a Liouvillian eigenvalue, then is also an eigenvalue [82].
The quasimomenta satisfying the BAEs fall into one of the following four cases in the limit . Here, we briefly summarize the resuls and details can be found in [62].
1.Free band: When both S-matrices are regular , the quasimomenta and are complex, with imaginary parts vanishing as , i.e and , for . These solutions exist for all parameters.
2. Bound mode band: When , quasimomenta are and for .
If , it is a bound mode; otherwise, it is not a valid solution. When both , another mode with exists.
3. Diffusive band: When , the eigenspectra are with . Gapless excitations with [81, 82], dominate the late time dynamics. These solutions with exist for all parameters.
4. Double-bound mode: There are two solutions where the particles are boundary-localized. When , the eigenvalues can be solved perturbtively , where and . When , the eigenvalues solved with an effective two-site Liouvillian are .
Fig. 1 shows the spectrum of the Liouvillian in the Kondo and bound regimes. We present the cross-over diagram that delineates these two regimes in Fig. 1(a) through the value of the real part of the Double bound mode . A non-zero signifies the bound mode regime, while a zero value indicates the Kondo regime.
We further analyze the spectrum’s real and imaginary parts in the Kondo regime in Fig. 1(a) and (b) as well as the bound mode regime in Fig. 1 (c). We see the free band (spectrum 1; circles), the bound mode band (spectrum 2; crosses), and the diffusive band (spectrum 3; star) in the Kondo regime, where the majority of the weight in lies in the diffusive band. In the bound mode regime, we also see the appearance of double bound modes (spectrum 4; diamonds) that have the largest weight in the coefficients and
hence they dominate the early-time impurity dynamics .
In [71], we showed that at the XX point without noise, there is a boundary phase transition at . When , the impurity is screened by many-body Kondo effect whereas when , the impurity is screened by single particle bound mode localized at the edge. Here, we turn to study the dynamical characterization of the boundary phases. In the bound mode phase (), the time-dependent magnetization shows oscillatory behavior when . However in the Kondo phase (), the impurity magnetization leaks to the bulk with the time scale proportional to the Kondo time . Detailed calculation showing this sharp phase transition in dynamical quantities is in [62].
When noise is introduced (), there is no longer a sharp phase transition. Instead, the impurity magnetization decays similarly in both regimes at late times, exhibiting the Zeno effect, where the dynamics slows down as noise strength increases. This behavior is driven by the diffusive band as it has largest overlap . The modes in the diffusive band have , so the dynamics slows down as increases and exhibits the quantum Zeno effect. In contrast, the early-time dynamics differ between the two regimes. In the bound mode regime, the impurity dynamics is dominated by the double-bound mode which has the largest overlap and . Thus, the impurity decays faster as increases. However, in the Kondo regime, the early-time dynamics are also influenced by the diffusive band—similar to the long-time limit where the decay rate decreases with increasing noise, resulting the quantum Zeno effect. The model parameters and corresponding to these two regimes are shown in Fig. 1(a).
Numerical results for starting from the Néel state at the XX point are shown in the bottom panel of Fig. 2 with two situations when the system is in the Kondo regime and in the bound-mode regime .
As seen from the data, for the early time in the bound-mode regime, the noise adds a decaying envelope that is accelerated as the dephasing rate is increased. This is an expected result because the boundary mode is essentially a two-qubit state, and qubits subjected to dephasing noise decay faster as the decay rate increases. In the Kondo regime or after the double bound mode decays i.e. , the noise slows the leakage of magnetization and leads to a quantum Zeno effect as the dephasing rate increases. In the long time limit, independent of the boundary coupling , the impurity magnetization falls off as with in both regimes. This behavior of the impurity in the Kondo regime is consistent with the heating of the bulk induced by Hermitian Lindblad operators [89, 90, 91], driving the impurity to (asymptotically free) weak coupling behavior [92, 93, 94, 95].
Figure 3: Data collapse for with in the XX case (filled markers, with system size ) and XXX case (hollow markers, with system size ) starting from a domain wall initial state (as shown in the color plot in the right inset for XX case with ) under various dephasing rates . Different marker shapes correspond to the spin configuration at different time slices, and different marker color refers to different dephasing rates. The left inset shows the dependence of the diffusion constant on (filled markers for XX case and hollow for XXX case) and the dashed lines are the analytical result for the noisy XX chain .
Although no noise is applied at the impurity site, the Zeno effect on impurity is present in the Kondo regime at all times and at late times () in the bound mode phase, because noise changes the bulk transport property from ballistic when the system is at the XX point to diffusive [96] with diffusion constant as shown in [82].
Heisenberg (XXX) chain: Next, we examine the Trotterized dynamics of the interacting circuit corresponding to the isotropic Heisenberg model () with a boundary impurity. The equilibrium boundary phases have been studied in [72, 76, 97]. We show here that the bound mode and Kondo phases also exist dynamically, manifesting themselves in a similar way as in the XX case: the bound mode phase exists when where the impurity magnetization oscillates in time, and the Kondo phase when where it decays as shown in Fig. 2 in the noiseless case. The bulk transport properties in the noiseless case belong to the KPZ universality class as shown in [98, 99].
Turning on the noise we find that it impacts both the bulk and boundary dynamics similarly as in the XX case discussed above.
To show the bulk dynamics is diffusive for XXX case, we prepare the system in a domain wall state where all sites on the left (right) half are up (down) and compute the magnetization as time evolves. When the system is diffusive, in the continuum limit as lattice spacing and system size , the magnetization satisfies the diffusion equation i.e. , with initial condition .
The numerical results for the dynamics of both XX and XXX cases are shown in Fig. 3, where all the data collapses on the curve with which is the solution of diffusion equation from the domain wall initial configuration. For large values of , we find clear numerical evidence that as in the XX case, but for smaller values of , it deviates (as shown in the left inset in Fig. 3) due to finite size and time effects [62] and we therefore expect that the relation holds for each value of considered.
At the same time, the bulk diffusive dynamics reflects itself in the impurity dynamics as shown in the impurity magnetization dynamics from the initial Néel state with and without noise in the upper panel of Fig. 2.
In the presence of noise, both Kondo and bound-mode regimes exist, and the impurity in each regime shows a similar behavior as in the free XX chain. In the bound mode regime, the impurity magnetization initially decays faster as the noise strength is increased, and after the bound modes vanish (), the impurity experiences the quantum Zeno effect whereas in the Kondo regime, the impurity magnetization exhibits the Zeno effect in all time scales.
Conclusion: We found that dephasing noise affects impurity behavior differently in the Kondo and bound mode regimes. In the Kondo regime, increasing slows the impurity magnetization decay, leading to a quantum Zeno effect. In the bound mode regime, magnetization initially decays faster but eventually exhibits the Zeno effect, similar to the Kondo regime. Additionally, the KPZ dynamics [98, 99] in the isotropic Heisenberg model become diffusive when the diffusion rate is introduced [100, 101].
Acknowledgment: We thank T. Giamarchi, V. Alba, F. Essler, and D. Huse for useful discussions. J.H.P. is partially supported by the Army Research Office Grant No. W911NF-23-1-0144 and the Office of Naval Research grant No. N00014-23-1-2357 (J.H.P.).
References
Preskill [2018]J. Preskill, Quantum computing in the
nisq era and beyond, Quantum 2, 79
(2018).
Lau et al. [2022]J. W. Z. Lau, K. H. Lim, H. Shrotriya, and L. C. Kwek, Nisq computing: where are we and where
do we go?, AAPPS
bulletin 32, 27
(2022).
Bharti et al. [2022]K. Bharti, A. Cervera-Lierta, T. H. Kyaw, T. Haug, S. Alperin-Lea, A. Anand, M. Degroote, H. Heimonen, J. S. Kottmann, T. Menke, et al., Noisy
intermediate-scale quantum algorithms, Reviews of Modern Physics 94, 015004 (2022).
Krantz et al. [2019]P. Krantz, M. Kjaergaard,
F. Yan, T. P. Orlando, S. Gustavsson, and W. D. Oliver, A quantum engineer’s guide to superconducting qubits, Applied physics
reviews 6 (2019).
Devoret et al. [2004]M. H. Devoret, A. Wallraff, and J. M. Martinis, Superconducting qubits: A short
review, arXiv
preprint cond-mat/0411174 (2004).
Kjaergaard et al. [2020]M. Kjaergaard, M. E. Schwartz, J. Braumüller, P. Krantz, J. I.-J. Wang,
S. Gustavsson, and W. D. Oliver, Superconducting qubits: Current state of play, Annual Review of
Condensed Matter Physics 11, 369 (2020).
Blatt and Roos [2012a]R. Blatt and C. F. Roos, Quantum simulations with
trapped ions, Nature Physics 8, 277
(2012a).
Bruzewicz et al. [2019]C. D. Bruzewicz, J. Chiaverini, R. McConnell, and J. M. Sage, Trapped-ion quantum
computing: Progress and challenges, Applied Physics Reviews 6 (2019).
Leibfried et al. [2003]D. Leibfried, R. Blatt,
C. Monroe, and D. Wineland, Quantum dynamics of single trapped ions, Reviews of Modern Physics 75, 281 (2003).
Singer et al. [2010]K. Singer, U. Poschinger,
M. Murphy, P. Ivanov, F. Ziesel, T. Calarco, and F. Schmidt-Kaler, Colloquium: Trapped ions as quantum bits: Essential numerical tools, Reviews of Modern
Physics 82, 2609
(2010).
Henriet et al. [2020]L. Henriet, L. Beguin,
A. Signoles, T. Lahaye, A. Browaeys, G.-O. Reymond, and C. Jurczak, Quantum computing with neutral atoms, Quantum 4, 327 (2020).
Jessen et al. [2004]P. Jessen, I. Deutsch, and R. Stock, Quantum information processing with trapped
neutral atoms, Quantum Information Processing 3, 91 (2004).
Wintersperger et al. [2023]K. Wintersperger, F. Dommert, T. Ehmer,
A. Hoursanov, J. Klepsch, W. Mauerer, G. Reuber, T. Strohm, M. Yin, and S. Luber, Neutral
atom quantum computing hardware: performance and end-user perspective, EPJ Quantum
Technology 10, 32
(2023).
Blais et al. [2020]A. Blais, S. M. Girvin, and W. D. Oliver, Quantum information processing and
quantum optics with circuit quantum electrodynamics, Nature Physics 16, 247 (2020).
Clerk et al. [2020]A. Clerk, K. Lehnert,
P. Bertet, J. Petta, and Y. Nakamura, Hybrid quantum systems with circuit quantum electrodynamics, Nature Physics 16, 257 (2020).
Blais et al. [2021]A. Blais, A. L. Grimsmo,
S. M. Girvin, and A. Wallraff, Circuit quantum electrodynamics, Reviews of Modern Physics 93, 025005 (2021).
Haroche et al. [2020]S. Haroche, M. Brune, and J. Raimond, From cavity to circuit quantum
electrodynamics, Nature Physics 16, 243
(2020).
Schmidt and Koch [2013]S. Schmidt and J. Koch, Circuit qed lattices: Towards quantum
simulation with superconducting circuits, Annalen der Physik 525, 395 (2013).
Breuer and Petruccione [2002]H.-P. Breuer and F. Petruccione, The theory of open
quantum systems (Oxford University Press, USA, 2002).
Banerjee et al. [2018]S. Banerjee, S. Banerjee, and Ahmad, Open quantum systems (Springer, 2018).
Schaller [2014]G. Schaller, Open quantum systems
far from equilibrium, Vol. 881 (Springer, 2014).
Mirrahimi and Rouchon [2015]M. Mirrahimi and P. Rouchon, Dynamics and control of
open quantum systems, Lecture notes (2015).
Williams [2010]C. P. Williams, Explorations in
quantum computing (Springer Science & Business
Media, 2010).
Nielsen and Chuang [2010]M. A. Nielsen and I. L. Chuang, Quantum computation and
quantum information (Cambridge university press, 2010).
Ljubotina et al. [2019a]M. Ljubotina, L. Zadnik, and T. Prosen, Ballistic spin transport in a
periodically driven integrable quantum system, Physical review letters 122, 150605 (2019a).
Paletta and Prosen [2024]C. Paletta and T. Prosen, Integrability of open
boundary driven quantum circuits, arXiv preprint arXiv:2406.12695 (2024).
Hudomal et al. [2024]A. Hudomal, R. Smith,
A. Hallam, and Z. Papić, Integrability breaking and bound states in google’s
decorated xxz circuits, PRX Quantum 5, 010316 (2024).
Hutsalyuk et al. [2024a]A. Hutsalyuk, Y. Jiang,
B. Pozsgay, H. Xu, and Y. Zhang, Exact spin correlators of integrable quantum circuits from algebraic
geometry, arXiv
preprint arXiv:2405.16070 (2024a).
Miao and Vernier [2023a]Y. Miao and E. Vernier, Integrable quantum circuits from the
star-triangle relation, Quantum 7, 1160 (2023a).
Gombor and Pozsgay [2024]T. Gombor and B. Pozsgay, Integrable deformations
of superintegrable quantum circuits, SciPost Physics 16, 114 (2024).
Perk and Au-Yang [2006]J. H. Perk and H. Au-Yang, Yang-baxter equations, arXiv preprint
math-ph/0606053 (2006).
Jimbo [1994]M. Jimbo, Introduction to the
Yang-Baxter equation, Vol. 17 (World Scientific River Edge, NJ, 1994).
Baxter [2007]R. J. Baxter, Exactly solved models in
statistical mechanics (Courier Corporation, 2007).
Reshetikhin [2010]N. Reshetikhin, Lectures on the
integrability of the six-vertex model, Exact methods in low-dimensional statistical physics
and quantum computing , 197 (2010).
Jimbo and Miwa [1994]M. Jimbo and T. Miwa, Algebraic analysis of solvable
lattice models, Vol. 85 (American Mathematical Soc., 1994).
[39]V. E. Korepin, V. E. Korepin, N. Bogoliubov, and A. Izergin, Quantum inverse
scattering method and correlation functions, .
Hutsalyuk et al. [2024b]A. Hutsalyuk, Y. Jiang,
B. Pozsgay, H. Xu, and Y. Zhang, Exact spin correlators of integrable quantum circuits from algebraic
geometry, arXiv
preprint arXiv:2405.16070 (2024b).
Miao and Vernier [2023b]Y. Miao and E. Vernier, Integrable quantum circuits from the
star-triangle relation, Quantum 7, 1160 (2023b).
Zhang et al. [2024]J. Zhang, G. Xia, C.-W. Wu, T. Chen, Q. Zhang, Y. Xie, W.-B. Su, W. Wu,
C.-W. Qiu, P.-x. Chen, et al., Observation of quantum strong mpemba
effect, arXiv
preprint arXiv:2401.15951 (2024).
Monroe et al. [2021]C. Monroe, W. C. Campbell, L.-M. Duan,
Z.-X. Gong, A. V. Gorshkov, P. W. Hess, R. Islam, K. Kim, N. M. Linke, G. Pagano, et al., Programmable quantum simulations of spin systems with trapped ions, Reviews of Modern
Physics 93, 025001
(2021).
Aspuru-Guzik and Walther [2012]A. Aspuru-Guzik and P. Walther, Photonic quantum
simulators, Nature physics 8, 285
(2012).
Bernien et al. [2017]H. Bernien, S. Schwartz,
A. Keesling, H. Levine, A. Omran, H. Pichler, S. Choi, A. S. Zibrov, M. Endres, M. Greiner,
et al., Probing many-body
dynamics on a 51-atom quantum simulator, Nature 551, 579 (2017).
Blatt and Roos [2012b]R. Blatt and C. F. Roos, Quantum simulations with
trapped ions, Nature Physics 8, 277
(2012b).
Bloch et al. [2012]I. Bloch, J. Dalibard, and S. Nascimbene, Quantum simulations with ultracold
quantum gases, Nature Physics 8, 267
(2012).
Gross and Bloch [2017]C. Gross and I. Bloch, Quantum simulations with ultracold
atoms in optical lattices, Science 357, 995 (2017).
Houck et al. [2012]A. A. Houck, H. E. Türeci, and J. Koch, On-chip quantum simulation
with superconducting circuits, Nature Physics 8, 292 (2012).
Joshi et al. [2022]L. K. Joshi, A. Elben,
A. Vikram, B. Vermersch, V. Galitski, and P. Zoller, Probing many-body quantum chaos with quantum simulators, Physical Review
X 12, 011018 (2022).
Semeghini et al. [2021]G. Semeghini, H. Levine,
A. Keesling, S. Ebadi, T. T. Wang, D. Bluvstein, R. Verresen, H. Pichler, M. Kalinowski, R. Samajdar, et al., Probing topological spin liquids on a programmable quantum
simulator, Science 374, 1242
(2021).
Zhang et al. [2017]J. Zhang, G. Pagano,
P. W. Hess, A. Kyprianidis, P. Becker, H. Kaplan, A. V. Gorshkov, Z.-X. Gong, and C. Monroe, Observation of a many-body dynamical phase transition with a 53-qubit
quantum simulator, Nature 551, 601 (2017).
Kawabata et al. [2023a]K. Kawabata, A. Kulkarni,
J. Li, T. Numasawa, and S. Ryu, Dynamical quantum phase transitions in sachdev-ye-kitaev
lindbladians, Physical Review B 108, 075110 (2023a).
Kawabata et al. [2023b]K. Kawabata, A. Kulkarni,
J. Li, T. Numasawa, and S. Ryu, Dynamical quantum phase transitions in sachdev-ye-kitaev
lindbladians, Phys. Rev. B 108, 075110 (2023b).
Lerner [2018]L. Lerner, Quantum zeno effect at
finite measurement strength and frequency, Physical Review A 98, 052132 (2018).
Vanhoecke and Schirò [2024]M. Vanhoecke and M. Schirò, Kondo-zeno crossover
in the dynamics of a monitored quantum dot, arXiv preprint arXiv:2405.17348 (2024).
Wang and Wang [2024]X. Wang and J. Wang, Mpemba effects in nonequilibrium open
quantum systems, arXiv preprint arXiv:2401.14259 (2024).
Chatterjee et al. [2023]A. K. Chatterjee, S. Takada, and H. Hayakawa, Quantum mpemba effect in a quantum dot
with reservoirs, Physical Review Letters 131, 080402 (2023).
Rylands et al. [2024]C. Rylands, K. Klobas,
F. Ares, P. Calabrese, S. Murciano, and B. Bertini, Microscopic origin of the quantum mpemba effect in integrable
systems, Physical Review Letters 133, 010401 (2024).
Joshi et al. [2024]L. K. Joshi, J. Franke,
A. Rath, F. Ares, S. Murciano, F. Kranzl, R. Blatt, P. Zoller, B. Vermersch, P. Calabrese, C. F. Roos, and M. K. Joshi, Observing the
quantum mpemba effect in quantum simulations, Phys. Rev. Lett. 133, 010402 (2024).
Schlimgen et al. [2022]A. W. Schlimgen, K. Head-Marsden, L. M. Sager, P. Narang, and D. A. Mazziotti, Quantum simulation of the lindblad
equation using a unitary decomposition of operators, Physical Review Research 4, 023216 (2022).
Tang et al. [2024]Y. Tang, P. Kattel,
J. H. Pixley, and N. Andrei, Supplementary information of ‘quantum zeno effect in
noisy integrable quantum circuits for impurity models’ (2024).
[63]B. Sutherland, Beautiful models: 70
years of exactly solved quantum many-body problems, .
Cherednik [1984]I. V. Cherednik, Factorizing particles
on a half-line and root systems, Teoreticheskaya i Matematicheskaya Fizika 61, 35 (1984).
Sklyanin [1988]E. K. Sklyanin, Boundary conditions for
integrable quantum systems, Journal of Physics A: Mathematical and General 21, 2375 (1988).
Wang et al. [2015]Y. Wang, W.-L. Yang,
J. Cao, and K. Shi, Off-diagonal Bethe ansatz for exactly solvable models (Springer, 2015).
Trotter [1959]H. F. Trotter, On the product of
semi-groups of operators, Proceedings of the American Mathematical Society 10, 545 (1959).
Suzuki [1991]M. Suzuki, General theory of fractal
path integrals with applications to many-body theories and statistical
physics, Journal
of mathematical physics 32, 400 (1991).
Whitfield et al. [2011]J. D. Whitfield, J. Biamonte, and A. Aspuru-Guzik, Simulation of
electronic structure hamiltonians using quantum computers, Molecular Physics 109, 735 (2011).
Kattel et al [2024]P. Kattel et al, Integrable and
non-integrable boundary defects in spin- anisotropic heisenberg
chain, in
preparation (2024).
Kattel et al. [2024a]P. Kattel, Y. Tang,
J. H. Pixley, and N. Andrei, The kondo effect in the quantum xx spin chain, Journal of Physics
A: Mathematical and Theoretical (2024a).
Wang [1997]Y. Wang, Exact solution of the open
heisenberg chain with two impurities, Physical Review B 56, 14045 (1997).
Furusaki and Hikihara [1998]A. Furusaki and T. Hikihara, Kondo effect in xxz spin
chains, Physical
Review B 58, 5529
(1998).
Hu and Pu [1998]Z.-N. Hu and F.-C. Pu, Two magnetic impurities in a spin
chain, Physical
Review B 58, R2925
(1998).
Laflorencie et al. [2008]N. Laflorencie, E. S. Sørensen, and I. Affleck, The kondo effect in spin
chains, Journal
of Statistical Mechanics: Theory and Experiment 2008, P02007 (2008).
Kattel et al. [2024b]P. Kattel, P. R. Pasnoori, J. H. Pixley, P. Azaria, and N. Andrei, Kondo effect in the isotropic
heisenberg spin chain, Phys. Rev. B 109, 174416 (2024b).
Lindblad [1976]G. Lindblad, On the generators of
quantum dynamical semigroups, Communications in Mathematical Physics 48, 119 (1976).
Gorini et al. [1976]V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, Completely
positive dynamical semigroups of n-level systems, Journal of Mathematical Physics 17, 821 (1976).
Sá et al. [2021]L. Sá, P. Ribeiro, and T. c. v. Prosen, Integrable nonunitary open quantum
circuits, Phys. Rev. B 103, 115132 (2021).
Su and Martin [2022]L. Su and I. Martin, Integrable nonunitary quantum
circuits, Physical Review B 106, 134312 (2022).
Medvedyeva et al. [2016a]M. V. Medvedyeva, F. H. L. Essler, and T. c. v. Prosen, Exact bethe ansatz
spectrum of a tight-binding chain with dephasing noise, Phys. Rev. Lett. 117, 137202 (2016a).
Alba [2023]V. Alba, Free fermions with dephasing
and boundary driving: Bethe ansatz results, arXiv preprint arXiv:2309.12978 (2023).
Ziolkowska and Essler [2020]A. A. Ziolkowska and F. Essler, Yang-baxter integrable
lindblad equations, SciPost Physics 8, 044 (2020).
Karnaukhov and Egorov [2005]I. Karnaukhov and B. Egorov, Exact solution of the
hubbard model with boundary hoppings and fields, Europhysics Letters 70, 218 (2005).
Prosen [2008]T. Prosen, Third quantization: a
general method to solve master equations for quadratic open fermi systems, New Journal of
Physics 10, 043026
(2008).
Faddeev [1996]L. Faddeev, How algebraic bethe
ansatz works for integrable model, arXiv preprint hep-th/9605187 (1996).
Li et al. [2014]Y.-Y. Li, J. Cao, W.-L. Yang, K. Shi, and Y. Wang, Exact solution of the one-dimensional hubbard model with arbitrary
boundary magnetic fields, Nuclear Physics B 879, 98 (2014).
Lieb and Wu [1968]E. H. Lieb and F.-Y. Wu, Absence of mott transition in an exact
solution of the short-range, one-band model in one dimension, Physical Review Letters 20, 1445 (1968).
Medvedyeva et al. [2016b]M. V. Medvedyeva, T. c. v. Prosen, and M. Žnidarič, Influence of dephasing on many-body
localization, Phys. Rev. B 93, 094205 (2016b).
Fischer et al. [2016]M. H. Fischer, M. Maksymenko, and E. Altman, Dynamics of a
many-body-localized system coupled to a bath, Phys. Rev. Lett. 116, 160401 (2016).
Levi et al. [2016]E. Levi, M. Heyl, I. Lesanovsky, and J. P. Garrahan, Robustness of many-body localization in the
presence of dissipation, Phys. Rev. Lett. 116, 237203 (2016).
Hewson [1997]A. C. Hewson, The Kondo problem to
heavy fermions, 2 (Cambridge
university press, 1997).
Kondo [2012]J. Kondo, The physics of dilute
magnetic alloys (Cambridge University Press, 2012).
Andrei et al. [1983]N. Andrei, K. Furuya, and J. Lowenstein, Solution of the kondo problem, Reviews of modern
physics 55, 331
(1983).
Tsvelick and Wiegmann [1983]A. Tsvelick and P. Wiegmann, Exact results in the
theory of magnetic alloys, Advances in Physics 32, 453 (1983).
Jin et al. [2020]T. Jin, M. Filippone, and T. Giamarchi, Generic transport formula for a system
driven by markovian reservoirs, Physical Review B 102, 205131 (2020).
Frahm and Zvyagin [1997]H. Frahm and A. A. Zvyagin, The open spin chain with
impurity: an exact solution, Journal of Physics: Condensed Matter 9, 9939 (1997).
Ljubotina et al. [2019b]M. Ljubotina, M. Žnidarič, and T. Prosen, Kardar-parisi-zhang physics in the quantum heisenberg magnet, Physical review letters 122, 210602 (2019b).
Rosenberg et al. [2024]E. Rosenberg, T. Andersen,
R. Samajdar, A. Petukhov, J. Hoke, D. Abanin, A. Bengtsson, I. Drozdov, C. Erickson, P. Klimov, et al., Dynamics of magnetization at infinite temperature in a heisenberg
spin chain, Science 384, 48
(2024).
Esposito and Gaspard [2005]M. Esposito and P. Gaspard, Emergence of diffusion in
finite quantum systems, Physical Review B—Condensed Matter and Materials Physics 71, 214302 (2005).
Eisler [2011]V. Eisler, Crossover between
ballistic and diffusive transport: the quantum exclusion process, Journal of
Statistical Mechanics: Theory and Experiment 2011, P06007 (2011).
Supplementary information of ‘Quantum Zeno Effect in Noisy Integrable Quantum Circuits for Impurity Models’
Yicheng Tang,1 Pradip Kattel,1 J. H. Pixley,1,2 and Natan Andrei1
1Department of Physics and Astronomy, Center for Material Theory,
Rutgers University, Piscataway, New Jersey, 08854, United States of America
2Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, NY 10010
In the supplementary material, we first map the circuit dynamics to the Lindbladian evolution. Then at the XX point, we show that the corresponding Lindbladian can be mapped to an integrable Hubbard Hamiltonian with imaginary onsite interaction in and boundary hopping defect in doubled Liouville-Fock space. We then solve the dynamics of the noiseless model with the free-fermion technique. After that, we solve the non-Hermition Hubbard model in the two-particle sector to discuss the dynamics in the presence of noise. We also discuss the derivation of the Bethe-Ansatz equations for the complex onsite Hubbard model in an arbitrary particle sector. In the final section, we show the details of numerical results.
I Mapping from Quantum circuits to Liouvillian dynamics
Consider the trigonometric six vertex R-matrix [36, 63]
(S.1)
which is a solution of the Yang-Baxter equation [86]
(S.2)
Notice that the Yang-Baxter equation remains satisfied if we shift where are arbitrary inhomogeneous parameters.
Let us consider the following two single-row transfer matrices
(S.3)
where are the labels for the physical space and is an auxiliary space.
Furthermore, denoting the physical space of impurity as , we define
(S.4)
such that the matrices satisfy the reflection equations
(S.5)
Now, we define the monodromy matrix
(S.6)
The trace of the monodromy matrix over the auxiliary space is defined as the double-row transfer matrix [66]
Such that when the transfer matrix Eq.(S.7) is evaluated at , we obtain
(S.9)
where
(S.10)
Here is the permutation matrix and explicit form of the matrix reads
(S.11)
These local unitary can be written as
(S.12)
(S.13)
upon identifying
(S.14)
Figure S.1: Illustrations of an integrable quantum circuit with spin-1/2 qubits in an initial classical Néel state (empty circles are qubits in the up states and black filled circles are in the down qubit state). Each rectangle represents a local unitary matrix; blue ones are the braid form of the trigonometric six-vertex R-matrix and red ones are the boundary reflection matrix .
The dashed rectangle shows one layer of the time evolution operator given by Eq.(1).
The circuit in Fig. S.1 corresponds to the unitary matrix of the following form
(S.15)
where in the Trotterization limit, we obtain an integrable Hamiltonian with boundary impurity of the form
(S.16)
where the parameters in two representations are related by and .
Notice that both and are expressed through a single free parameter hence the model is not integrable for arbitrary and . However, at the isotropic [76, 72] and XX [71] limit, the boundary and bulk anisotropy parameter becomes equal i.e. and becomes an arbitrary free parameter.
Taking into account the presence of noise while applying the gates, the density matrix of the qubits can be described by a discrete-time evolution as
(S.17)
with , for , and and are the Lindbladian jump operators. Here, we will concentrate solely on the case where . Taking , we obtain the Lindbladian-master equation [19]
(S.18)
where the Liouvillian operator can be written as
(S.19)
with and being the anticommutator and the commutator, respectively.
II XX chain point of the Hamiltonian
At the XX chain point i.e when , we end up with a Hamiltonian of a XX chain with an XX impurity [71]
(S.20)
After performing the Jordan-Wigner transformation, one can write the model in terms of Fermionic operators i.e. the Hamiltonian becomes
(S.21)
and the Liouvillian operator takes the form
(S.22)
The density matrix is an operator defined in the dimensional Hilbert space which can be written as . We can rewrite it as a vector in a dimensional Hilbert space, so-called the Liouville-Fock space, as without losing information [85]. Then the Liouvillian operator acts as an operator in the enlarged vector space in the following way
(S.23)
with the Liouvillian as superoperator
(S.24)
Here we introduce a new superoperator acting on the Liouville-Fock space in the following way:
(S.25)
If we rename and , and perform the unitary transformation
(S.26)
the Liouvillian operator becomes the Hubbard model with imaginary coupling [81]
(S.27)
where , and all other combinations anti-commute. In this formalism, the Lindbladian-master equation becomes an eigenvalue problem
(S.28)
One thing worthy to mention is that the dynamics of 2-point correlation function , under this mapping, can be shown to be governed by the 2-particle Schrodinger equation [81]. Particularly for the two-point correlation function we are interested in,
(S.29)
where we introduced a compact notation when , otherwise and .
III Noiseless XX-impurity Quench
Let us first discuss the dynamics of the model in the absence of the noise. We consider the quench dynamics of the free model governed by a quantum mechanical Hamiltonian
(S.30)
where when , other wise and is the discrete derivative. This first quantized Hamiltonian has eigenfunction
There exists a bound mode with , described by such that the wave function is normalizable in the thermodynamic limit, i.e. which gives with bound mode energy and gives . Likewise, gives .
After studying the spectrum and equilibrium physics, now we can study the quench dynamics of this problem. Since the model is non-interacting, we need to solve only the single-mode evolution and take the algebraic sum of each mode at the end of the calculation.
We prepare the initial state as , and we solve for . The dynamics of a wave function is given by the following equation
(S.34)
Let us look at the simplest case where we prepare the fermion (spin flip) at the impurity site initial and ask how the fermion number (magnetization) changes as time evolves, In the Kondo phase
(S.35)
with .
Changing the variable , we obtain
(S.36)
where the Kondo time is the only scale in this integral. We verify this result numerically as shown in Fig.S.2.
In the Bound Mode Phase, the same quantity can be expressed as the following
(S.37)
The last cosine term is from the bound mode with energy .
Figure S.2: Dynamics of the spin at the impurity site with initial Néel state. The system is in bound mode phase (left panel) where the analytic curve of the time-dependent spin profile is given by Eq. (S.40) and numerical cure is calculated for with exact diagonalization. The system is in the Kondo phase with (right panel) where the analytic curve of time-dependent spin profile is given by Eq. (S.40) and numerical cure is calculated for with exact diagonalization.
In the Kondo phase, the wavefunction at the spin, if we flip initially, is given by
(S.38)
In the bound mode phase,
(S.39)
The magnetization at site with spin initially flipped is given by
(S.40)
with or depending on which phase the system is in.
IV Liouvillian Spectrum
In the previous section we solved XX-impurity Model quenched from an initial state with only some spins flipped, i.e. . Now we want to ask the question of how noise affects this model. Using Eq. (S.29), we need to study the two-particle sector of this model to solve for the dynamics of the two-point correlation function.
In the Liouville-Fock space, our initial state (initial condition of the correlation matrix) with spin-flip at site corresponds to the initial density matrix . If we initially flip more spins, we just need to superpose them, i.e. . Due to the symmetry of the Liouvillian operator, the number of fermions and total magnetization are conserved separately. Thus, we only need to diagonalize the eigenvalue problem in Eq. (S.28) in two particles and one spin-flip sector.
The time-independent Schrödinger equation read from Eq.(S.29) is
(S.41)
with eigenvector
(S.42)
where is the normalization factor such that and is the symmetrizer. Note if the transformation in Eq.(S.26) is not performed, then there is a factor multiplying the wavefunction shown in the main text.
The amplitudes are related as follows:
(S.43)
where is given by the boundary condition
which gives
(S.44)
and the bulk interaction scattering matrix is that of Hubbard models given by
(S.45)
If we impose the boundary condition on a finite chain, then and are eigenstates if they satisfy the Bethe Ansatz equations
(S.46)
(S.47)
which when written explicitly takes the form
(S.48)
The solution of and can be written as and .
IV.1 1. The free band (Roots with vanishing Imaginary part)
By taking the norm of Eq. (S.48), it is easy to see that
(S.49)
and the same for . One can see that as long as the scattering matrix is not singular, the imaginary part of the Bethe roots vanishes in the thermodynamic limit, and results in a real band of quasimomenta. In the thermodynamic limit, the imaginary part of the roots is given by replacing and .
IV.2 Roots with Finite Imaginary Part
As mentioned earlier, when the scattering matrix is singular, the Bethe roots have a non-vanishing imaginary part in the thermodynamic limit. Now let us consider different situations with different parts of the S-matrix being singular.
Figure S.3: The upper two panels show the Liouvillian spectrum (top-left) and amplitude for the wavefunction in the double-bound mode (top-right) for , and . The perturbative calculation matches the exact diagonalization result better since is small. The lower two panels show the Liouvillian spectrum (bottom-left) and amplitude for the wavefunction in the double-bound mode (bottom-right) for , and . The two-site Hamiltonian calculation matches the exact diagonalization result better since is big. Also, wavefunctions for the both cases show that the particles are localized at the boundary more strongly and when is bigger.
IV.3 2. Boundary Mode Band
In the case where the boundary matrix is singular, the bound mode exists and the singular term can be written as
(S.50)
which leads to the constraint . Changing the variable to , the equation becomes
(S.51)
which has a solution
(S.52)
under the assumption .
Plugging in into the other equation gives the solution of as
(S.53)
In the thermodynamic limit, and .
IV.4 3. Diffusive Band
In the limit, without lack of generality, we choose . Then the bulk scattering matrix being singular indicates the solutions
satisfy the string hypothesis in the way . These solutions form a band of solutions whose corresponding eigenvalues have a vanishing real part. Now we define , then the S-matrix being singular also gives the relationship between and as
(S.54)
As a result, the eigenenergy is , and the Bethe Ansatz equation can be written as the form of a Bethe-Gaudin-Takahashi (BGT) equation [82] as following
(S.55)
To solve this equation, it is easier to change the variable and take the logarithm on both sides of the equation which gives
(S.56)
with . This equation which look extremely simple gives the spectrum , which match the numerical solution extremely well as shown in Fig. S.3.
This figure shows that the analytical calculation fits the exact diagonalization reasonably well when which the system is experiencing a true diffusive bulk. However, when , the diffusive band mixes with the free band and results in a worse but reasonable prediction because we are not taking care about the contribution from the impurity.
Noticing gives the constraint , and all solutions beyond this limit are not physical and should be dropped. Since , the largest imaginary energy is when , which gives and thus resulting the energy . This is the lower bound of the dissipation band, and when which can only happen when the dissipation band mixes with the free band.
IV.5 4. Double-Boundary Modes
As discussed in the main text, there is another solution when both the boundary S-matrix and the bulk S-matrix are singular. Such a mode physically describes two particles (in the imaginary Hubbard model) that are both exponentially localized at the of the chain and they experience a strong imaginary Hubbard interaction. The existence of these two states is discussed in the main text with a phase diagram in Fig.1(d). However, we are not able to separate these two singularities, so we treat them in two limits. When , the double-bound mode is perturbatively connected to the case when both particles form a noiseless bound mode, which means that
(S.57)
with or and being the normalization factor. Then the eigenvalue of this double-bound mode reads
(S.58)
and when , the bound mode is more localized on the boundary (decays faster) as one can see from the single boundary bound mode. Then we can treat it by solving a two-site Hamiltonian () as
(S.59)
which has two boundary bound mode with energy . It is a double-bound mode only when which gives the boundary between the Kondo phase and the bound mode phase at the large limit as .
V N-particle Bethe Ansatz
As mentioned earlier, the XX chain with boundary impurity maps to the open boundary condition Hubbard model with boundary impurity given by
(S.60)
As shown in [71], the boundary scattering matrix is
(S.61)
and the Bulk S-matrix, as first calculated by Lieb and Wu for real values of U, is given by
(S.62)
where is the permutation operator.
Now, the usual Bethe Ansatz routine can be employed. It is very important to notice that the boundary matrix does not depend on the spin variable but the bulk matrix depends on the spin variable (through the permutation operator). The Yang-Baxter equation is therefore trivially satisfied, thereby proving that the model is integrable.
Following [87, 66], we construct the N-particle Bethe Ansatz of the model using the functional Bethe Ansatz method. Let us consider the rational solution of the six-vertex model
(S.63)
where we will consider the crossing parameter .
Then for the open boundary condition, we consider the double-row transfer matrix
Now, we define the monodromy matrix
(S.64)
The trace of the monodromy matrix over the auxiliary space is defined as the double-row transfer matrix
(S.65)
It is quite easy to see that the transfer matrix forms a one-parameter family of commuting operators i.e.
(S.66)
First notice that using coordinate Bethe Ansatz, one can immediately construct the eigenvalue problem for the transfer matrix
(S.67)
Upon identification of the relation between the spectral parameter valued double-row transfer matrix and the transfer matrix obtained from the coordinate Bethe Ansatz via
(S.68)
Once again, the boundary matrix is independent of a spin variable; it is just a phase that trivially commutes with bulk matrices.
The eigenvalue of the transfer matrix satisfies the relation
(S.69)
where the function is
(S.70)
Imposing the regularity condition of the relation immediately gives the Bethe equation
(S.71)
To make the equations more symmetric, we change the variable , such that the Bethe equations becomes
If one is interested in the many-point correlation function, one needs to solve the many-particle BAE given by Eq. (S.75) and Eq. (S.76). In the main text, we only considered the two-point correlation function, and we leave the general problem of the 2n-point correlation function for future work. The Bethe equations for the real Hubbard interaction were first obtained in [84].
Figure S.4: The magnetization for time and site in the XXX case with system size (right panel) and XX case with system size (left panel) for dephasing rate .
VI Numerical details
Here we will show the detailed calculation to get the diffusion constant numerically for the Heisenberg (XXX) and XX case. We start from a Doman wall configuration initially with and . Then we evolve the state via the Lindbladian master equation (S.19) from the initial state . Then we calculated the magnetization profile for each time and site . The results for in the XX case and XXX case are shown in Fig. S.4. For the XXX case, the domain wall hit the boundary when and resulting in a finite size effect.
If the bulk is diffusive, then in the limit lattice spacing and system size should satisfy the diffusion equation . The initial condition is a domain wall meaning and the solution reads . This is obvious that for different times and spaces, the function collapses on the curve with a normalized variable in the way that . The system is not perfectly diffusive for early-time dynamics and due to the finite size effect, we cannot evolve the system for an extremely long time. Therefore, we
will fit the magnetization to for a finite time to get the diffusion constant . To get the long-time behavior, we expand and scale the diffusion constant when . The fitting and scaling are shown in Fig.S.5. We found for both XX and XXX case which is the diffusion constant for long-time dynamics is for big gamma. As for small , due to the finite size effect, we are not able to conclude for the XXX case. For the XX case, where we can go for larger system sizes, we can show .
Figure S.5: Scaling of the diffusion constant with time and for the XXX circuit. The left panel shows as a fit parameter obtained by fitting the spin profile to the . The right panel shows as a function for different times and we scale by taking . The dashed line is the fitting for and different colors for different .