Operator-Space Entanglement Entropy in Random Fermionic Gaussian Dynamics

arXiv:2610.21723 · quant-ph, cond-mat.stat-mech, math-ph, math.MP · Submitted 2026-10-15 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: I'm Kai.

Mira: And I'm Mira, here with Lev.

Kai: Today's paper: "Operator-Space Entanglement Entropy in Random Fermionic Gaussian Dynamics".

Mira: The gist The study derives an exact expression for time-dependent bipartite operator-space entanglement entropy in random Gaussian dynamics modeled by a Brownian SYK2 model,

Kai: First, who's behind it and why it matters.

Title and authors: Kai: So we're looking at this paper, "Operator-Space Entanglement Entropy in Random Fermionic Gaussian Dynamics." It’s about how entanglement behaves when the system itself is evolving randomly under certain rules.

Mira: Yeah, and the title says they are studying something called operator-space entanglement entropy, which is a way to measure entanglement not just for states, but for operators themselves.

Lev: That sounds like it might be very abstract stuff if you're thinking about running simulations on actual hardware right now.

Kai: It’s actually trying to connect that concept to real dynamics, specifically using a model called the Brownian SYK2 model which has random evolution described by Ito calculus and a covariance matrix formalism.

Mira: And they are aiming for something more precise than just saying entanglement grows linearly for chaotic systems; they want an exact expression for the time-dependent bipartite OSEE in the thermodynamic limit.

Lev: An exact expression in the thermodynamic limit is always a big ask, especially when you're dealing with stochastic processes where things get messy quickly.

Kai: Exactly, and they use a mapping to a domain-wall melting problem to get there, which they say helps them derive an analytic characterization of the corresponding gaussian Page curve emerging at late times.

The paper's summary: Mira: They are essentially taking this complex random Gaussian dynamics setup and using that domain-wall melting mapping to solve it exactly for the bipartite OSEE at any time, in the limit where you have an infinite system.

Kai: So what they’re doing is finding a way to track how much entanglement operators have between two parts of the system as this random process unfolds over time.

Lev: Tracking that evolution analytically sounds like it needs a really compact way to handle all those stochastic noise terms and the covariance matrix updates.

Mira: They also get this analytic characterization of the gaussian Page curve that appears at late times, which is important because it gives us a specific prediction for how entanglement scales in that chaotic regime.

Kai: And they note that these curves they find differ from those predicted by fully Haar-random many-body states, which is a key distinction.

Lev: So the result isn't just another curve you plug into a simulation; it’s something fundamentally different from what you get when you assume perfect randomness across the whole system.

Mira: Right, and they are also exploring how this OSEE relates to non-gaussianity measures in the Heisenberg picture, suggesting an upper bound based on operator-space nonstabilizerness.

The paper's improvements: Kai: One of the main things they push is that they found a way to relate the operator’s entanglement entropy directly to a measure of non-gaussianity in the Heisenberg picture, giving an upper bound on OSEE based on operator-space nonstabilizerness.

Mira: That seems like a way to constrain the complexity without having to calculate the full OSEE, which they say is often limited by how hard it is to compute beyond integrable dynamics.

Lev: Constraining it by nonstabilizerness sounds useful for error correction because if you know the measure of nonstabilizerness, you know how much entanglement you're dealing with before things get too complicated.

Kai: They also map the computation onto studying the spectrum of the Choi covariance matrix, which simplifies finding that exact solution for the time evolution.

Mira: That mapping is helpful because it reduces the problem from tracking complex operator evolution to just analyzing a spectrum of a matrix, which they say is more tractable.

Lev: If you can compress that into a nonlinear partial differential equation for the resolvent G, like they mention, that makes sense for finding an exact solution.

Conclusion: Kai: So to wrap up the "Operator-Space Entanglement Entropy in Random Fermionic Gaussian Dynamics" paper, they successfully derived an exact expression for the time-dependent bipartite OSEE in the thermodynamic limit using a mapping to domain-wall melting.

Mira: And they provided an analytic characterization of that gaussian Page curve at late times, which is different from what you see in fully Haar-random states.

Lev: For someone building real hardware, it means that understanding these scaling laws helps us predict the complexity overhead when running these kinds of dynamics on actual systems.

Kai: And they also suggested a way to relate OSEE to non-gaussianity measures, offering an upper bound based on operator-space nonstabilizerness.

Mira: It’s a solid result because it tackles the difficulty of calculating OSEE beyond integrable systems by finding an analytic path through the dynamics.

Lev: It shows that even in complex random environments, there are structured ways to extract information about entanglement that aren't just bounded by simple logarithmic growth or linear growth.

Kai: So they set up a framework where we can study operator-space non-Gaussianity as a proxy for how entangled the system really is.

Mira: It’s motivating because it suggests that these concepts might be related in a way that helps us understand the nature of chaotic dynamics in quantum systems.

Pablo Bayona-Penaa, –Elena Gramellinia, –Lorenzo Pirolia

Dipartimento di Fisica e Astronomia, Universita di Bologna · INFN, Sezione di Bologna

quant-ph, cond-mat.stat-mech, math-ph, math.MP

Submitted: 2026-10-15

Updated: 2026-10-15

Comments: 36 pages, 2 figures

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 79/100

The gist: The gist The study derives an exact expression for time-dependent bipartite operator-space entanglement entropy in random Gaussian dynamics modeled by a Brownian SYK2 model, providing an analytic

Key concepts

Operator-Space Entanglement Entropy (OSEE): OSEE extends standard entanglement entropy to operators by considering their Choi state in a doubled Hilbert space. It serves as a measure to distinguish between integrable and chaotic dynamics, allowing researchers to quantify how 'entangled' an operator is.
Brownian SYK2 Model: This model describes the evolution of fermionic systems under random Gaussian dynamics. It uses a time-dependent two-local Hamiltonian with random parameters, where the evolution follows stochastic processes best handled by Ito calculus and covariance matrix formalism.
Domain-Wall Melting Mapping: The authors map the operator dynamics to a domain-wall (DW) melting problem. This mapping is crucial because it allows them to use existing tools from studying entanglement in quantum simple exclusion processes, leading to an exact solution for the time evolution of bipartite entanglement.
Gaussian Page Curve: This curve describes how the average OSEE behaves at late times ($t o ext{infinity}$). The study finds an analytic characterization of this curve for the Brownian SYK2 model, showing it differs from predictions made for fully Haar-random many-body states.

Terminology

Entanglement — 1× in this paper · explained in 4 episodes
Entanglement Entropy — 1× in this paper · explained in 3 episodes

Summary

The gist The study derives an exact expression for time-dependent bipartite operator-space entanglement entropy in random Gaussian dynamics modeled by a Brownian SYK2 model, providing an analytic characterization of the corresponding gaussian Page curve emerging at late times.

Operator-Space Entanglement Entropy (OSEE) and its Importance

The operator-space entanglement entropy (OSEE) extends the notion of bipartite entanglement entropy to operators by mapping them onto their corresponding Choi state in a doubled Hilbert space, where it is defined as the entanglement entropy of that state. This metric is important because it allows for an extension to operators of a notion of classical simulability, distinguishing between integrable and chaotic dynamics. For local free-fermionic dynamics, the bipartite entanglement entropy is bounded by a function growing at most logarithmically in time, while for chaotic local Hamiltonians it is expected to grow linearly.

Modeling Random Gaussian Dynamics

The work considers fermionic systems and initially unentangled Majorana strings, studying their evolution under random gaussian evolution modeled by a Brownian SYK2 model. The dynamics is defined by a time-dependent two-local Hamiltonian of the form H(t) = i/2 X i<j J ij (t) ψ iψ j. This model involves random parameters J ij (t), where E[J ij (tn)] = 0, and the evolution is described by a stochastic process best described by means of Ito calculus. Along each trajectory the dynamics remains gaussian, and it can be described using the covariance matrix formalism, where G(t) = R t G(0) R T t.

Mapping to Domain-Wall Melting and Exact Solutions

After a mapping to a domainwall (DW) melting problem, the authors extend recent tools developed in the study of the bipartite entanglement in the quantum simple exclusion process. This mapping leads to an analytic expression for the time evolution of the bipartite entanglement entropy at all times in the thermodynamic limit. Furthermore, they compute gaussian Page curves parameterized by both the subsystem fraction and the length of the initial Majorana string, which differ qualitatively from those of fully Haar-random many-body states.

Thermodynamic Limit and Page Curves

When t → ∞, the Brownian motion converges to the Haar random distribution over SO(L), simplifying the analysis of the average OSEE. Using Coulomb-gas approach, they obtain a stationary gaussian Page curve E[Sop(l)] L = q Z lambda+ / lambda− rho⋆(λ) h2(λ) dλ. This curve differs from the predictions obtained for fully Haar-random many-body states.

Dynamics and Finite Size Corrections

The dynamics are analyzed by exploiting the stochastic differential equation (SDE) for the eigenvalues dλ i, which involves terms related to pairwise eigenvalue interaction. The resulting moment hierarchy is compressed into a nonlinear PDE for the resolvent G, which is equivalent to a compact nonlinear partial differential equation. The solution for the balanced case q = θ = 1/2 yields the exact parametric solution of the thermodynamic resolvent equation.

Implications for Operator-Space Non-Gaussianity

The work discusses the possibility of relating OSEE with non-gaussianity measures in the Heisenberg picture, suggesting an upper bound on OSEE in terms of operator-space nonstabilizerness. They conclude that a natural candidate for this measure is the total Majorana weight W(O) = X S S a S squared, although this measure is not a faithful monotone.

Conclusions and Future Directions

The study successfully identified a mapping between the Choi-state representation of the operator and that of DW states, reducing OSEE computation to studying the spectrum of the Choi covariance matrix. The difference between random-state and random-operator entanglement dynamics lies on finite size corrections stemming from different underlying group structures. Future work is suggested for obtaining an alternative derivation using replica calculation techniques or extending the approach to quantities like out-of-time order correlators (OTOCs).

Appendix Details

The analysis in Appendix A summarizes the steps leading from the finite-size Jacobi ensemble to Eq. (38), where the stationary density is obtained by minimizing an energy functional. Appendix B provides details on the real matrix Jacobi process, showing how Eq. (43) leads to a closed matrix Jacobi process. Appendix D extends the characteristic construction to arbitrary subsystem fraction q = l/L and filling fraction θ = M/L, yielding an exact parametric solution for the thermodynamic resolvent equation. The final entropy formula is given by Eq. (D.61), which is an exact one-dimensional quadrature for the thermodynamic OSEE at arbitrary subsystem and filling fractions due to the invariance of OSEE under q ↔ θ.

References Cited

[1] Zanardi P 2001 Phys. Rev. A 63(4) 040304(R) URL https://link.aps.org/doi/10.1103/PhysRevA.63.040304

[2] Prosen T and Pizorn I 2007 Phys. Rev. A 76(3) 032316 URL https://link.aps.org/doi/10.1103/PhysRevA.76.032316

[4] Dubail J 2017 J. Phys. A: Math. Theor. 50 234001 URL https://iopscience.iopscience/article/

[5] Nielsen M A, Chuang I L et al. 2000 Quantum computation and quantum information vol 1 (Cambridge university press Cambridge) URL https://www.cambridgeuniversitypress.com/

[6] Verstraete F and Cirac J I 2006 Phys. Rev. B 73(9) 094423 URL https://link.aps.org/doi/10.1103/PhysRevB.73.9

[8] Schollwock U 2011 Ann. Phys. 326 96–192 URL https://www.sciencedirect.com/science/article/abs/pii/S0003491610001752

[9] Cirac J I, Perez-Garcia D, Schuch N and Verstraete F 2017 Ann. Phys. 378 100–149 URL https://www.sciencedirect.com/science/article/abs/pii/S0003491616303013

[58] Demni N 2008 Journal of Theoretical Probability 21 118–143 URL https://www.jstor.org/stable/24904084

[59] Demni N, Hamdi T and Hmidi T 2012 Indiana Univ. Math. J. 61 1351–1368 URL https://www.jstor.org/stable/24904084

[67] Demni N and Hamdi T 2025 arXiv:2501.07552 URL https://arxiv.org/abs/2501.07552

[68] Demni N, Gilliers N and Hamdi T 2025 arXiv preprint arXiv:2511.02758 URL https://arxiv.org/abs/arXiv:2511.02758

[69] White G A L, Eisert J and Dowling N 2026 arXiv:2609.38300 URL https://arxiv.org/abs/2609.

Improvements for AI systems

Bold header for OSEE dynamics modeling: The improved AI system can model operator-space entanglement entropy (OSEE) evolution under random gaussian dynamics, specifically using the Brownian SYK2 model to derive an exact expression for the time-dependent bipartite OSEE in the thermodynamic limit.
Bold header for asymptotic behavior prediction: The system can analytically characterize the gaussian Page curve emerging at late times, providing predictions that are distinct from those of fully Haar-random many-body states, as shown in Figure 1.
Bold header for non-Gaussianity relation: The AI system can discuss the possibility of relating the OSEE with non-gaussianity measures in the Heisenberg picture, potentially constraining an upper bound using a natural extension of non-gaussianity to operator space.