The dynamics of quantum many-body systems in the chaotic regime are of particular interest due to the associated phenomena of information scrambling and entanglement generation within the system. While these systems are typically intractable using traditional numerical methods, an effective framework can be implemented based on dual-unitary circuits which have emerged as a minimal model for maximally chaotic dynamics. In this work, we investigate how individual two-body operators influence the global dynamics of circuits composed of dual-unitaries. We study their effect on entanglement generation while examining it from both bipartite and multipartite perspectives. Here we also highlight the significant role of local unitaries in the dynamics when paired with operators from the dual-unitary class, showing that systems with identical entangling power can exhibit a range of differing entanglement growth rates. Furthermore, we present calculations establishing time-step-dependent lower bounds, which depend on both the initial state and the entangling power of the constituent operators. Finally, we find that time-evolving an initial state composed of pair products generates a state with nearly maximal multipartite entanglement content, approaching the bounds established by Absolutely Maximally Entangled (AME) states.
The dynamics of quantum many-body systems in the chaotic regime are of particular interest due to the associated phenomena of information scrambling and entanglement generation within the system. While these systems are typically intractable using traditional numerical methods, an effective framework can be implemented based on dual-unitary circuits which have emerged as a minimal model for maximally chaotic dynamics. In this work, we investigate how individual two-body operators influence the global dynamics of circuits composed of dual-unitaries. We study their effect on entanglement generation while examining it from both bipartite and multipartite perspectives. Here we also highlight the significant role of local unitaries in the dynamics when paired with operators from the dual-unitary class, showing that systems with identical entangling power can exhibit a range of differing entanglement growth rates. Furthermore, we present calculations establishing time-step-dependent lower bounds, which depend on both the initial state and the entangling power of the constituent operators. Finally, we find that time-evolving an initial state composed of pair products generates a state with nearly maximal multipartite entanglement content, approaching the bounds established by Absolutely Maximally Entangled (AME) states.
Dual-unitary circuits are a class of locally interacting quantum many-body systems displaying unitary dynamics also when the roles of space and time are exchanged. These systems have recently emerged as a remarkable framework where certain features of many-body quantum chaos can be studied exactly. In particular, they admit a class of ``solvable'' initial states for which, in the thermodynamic limit, one can access the full nonequilibrium dynamics. This reveals a surprising property: when a dual-unitary circuit is prepared in a solvable state the quantum entanglement between two complementary spatial regions grows at the maximal speed allowed by the local structure of the evolution. Here we investigate the fate of this property when the system is prepared in a generic pair-product state. We show that in this case, the entanglement increment during a time step is submaximal for finite times, however, it approaches the maximal value in the infinite-time limit. This statement is proven rigorously for dual-unitary circuits generating high enough entanglement, while it is argued to hold for the entire class.
Operator entanglement is a well-established measure of operator complexity across a system bipartition. In this work, we introduce a measure of the ability of a unitary channel to generate operator entanglement, representing an operator-level generalization of the state-space entangling power. This operator space entangling power is demonstrated to be linked to the scrambling properties of the unitary channel via the recently introduced concept of mutual averaged noncommutativity of quantum operator algebras. An upper bound for the operator space entangling power is identified that corresponds to unitary channels with scrambling properties akin to those of typical unitaries. Additionally, for Hamiltonian dynamics, we find that the short-time growth rate of the operator space entangling power matches the Gaussian scrambling rate of the bipartite out-of-time-order correlator, establishing a direct link between information scrambling and operator entanglement generation for short timescales. Finally, we examine the average growth of local operator entanglement across a symmetric bipartition of a spin chain. For dual-unitary circuits, a combination of analytical and numerical investigations demonstrates that the average growth of local operator entanglement exhibits two distinct regimes in relation to the operator space entangling power of the building-block gate.
Entanglement is a defining feature of many-body quantum systems and is an essential requirement for quantum computing. It is therefore useful to study physical processes which generate entanglement within a large system, as they maybe replicated for applications involving the said requirements in quantum information processing. A possible avenue to maximize entanglement generation is to rely on the phenomena of information scrambling, i.e. transport of initially localized information throughout the system. Here the rationale is that the spread of information carries with it an inherent capacity of entanglement generation. Scrambling greatly depends upon the dynamical nature of the system Hamiltonian, and the interplay between entanglement generation and information scrambling maybe investigated taking a chain of interacting spins on a one dimensional lattice. This system is analogous to an array of qubits and this relative simplicity implies that the resulting unitary dynamics can be efficiently simulated using present-day cloud based NISQ devices. In our present work, we consider such a spin model which is made up of nearest and next nearest neighbor XXZ Model, along with an introduced coupling term λ. This coupling term serves as a tuning parameter which modifies the dynamical nature of the system from the integrable to the quantum chaotic regime. In order to quantify the entanglement generated within the system we use the more general multipartite metric which computes the average entanglement across all system bipartitions to obtain a global picture of the entanglement structure within the entire system. Further to make our study initial state independent, we consider an ensemble of random initial states and evaluate the multipartite entanglement generated for each of these states under time evolution. We find clear differences in the entanglement measures as the system undergoes through different regimes. We also analyze the fluctuations present within our entanglement metric and find a distinct pattern as we head from the integrable to the chaotic regime.
We investigate operator dynamics and entanglement growth in dual-unitary circuits, a class of locally scrambled quantum systems that enables efficient simulation beyond the exponential complexity of the Hilbert space. By mapping the operator evolution to a classical Markov process, we perform Monte Carlo simulations to access the time evolution of local operator density and entanglement with polynomial computational cost. Our results reveal that the operator density converges exponentially to a steady-state value, with analytical bounds that match our simulations. Additionally, we observe a volume-law scaling of operator entanglement across different subregions, and identify a critical transition from maximal to sub-maximal entanglement growth, governed by the circuit’s gate parameter. This transition, confirmed by both mean-field theory and Monte Carlo simulations, provides new insights into operator entanglement dynamics in quantum many-body systems. Our work offers a scalable computational framework for studying long-time operator evolution and entanglement, paving the way for deeper exploration of quantum information dynamics.
The entanglement in operator space is a well established measure for the complexity of the quantum many-body dynamics. In particular, that of local operators has recently been proposed as dynamical chaos indicator, i.e. as a quantity able to discriminate between quantum systems with integrable and dynamics. For systems the local-operator entanglement is expected to grow linearly in time, while it is expected to grow at most logarithmically in the integrable case. Here we study local-operator entanglement in dual-unitary quantum circuits, a class of statistically solvable quantum circuits that we recently introduced. We identify a class of completely chaotic dual-unitary circuits where the local-operator entanglement grows linearly and we provide a conjecture for its asymptotic behaviour which is in excellent agreement with the numerical results. Interestingly, our conjecture also predicts a phase transition in the slope of the local-operator entanglement when varying the parameters of the circuits.
After quantum quenches in many-body systems, finite subsystems evolve nontrivially in time, eventually approaching a stationary state. In typical situations, the reduced density matrix of a given subsystem begins and ends this endeavor as a low-entangled vector in the space of operators. This means that if its entanglement in the operator space initially grows (which is generically the case), it must eventually decrease, describing a barrier-shaped curve. Understanding the shape of this ``entanglement barrier'' is interesting for three main reasons: (i) it quantifies the dynamics of entanglement in the (open) subsystem; (ii) it gives information on the approximability of the reduced density matrix by means of matrix product operators; and (iii) it shows qualitative differences depending on the type of dynamics undergone by the system, signaling quantum chaos. Here we compute exactly the shape of the entanglement barriers described by different R\'enyi entropies after quantum quenches in dual-unitary circuits initialized in a class of solvable matrix product states (MPS)s. We show that, for free (SWAP-like) circuits, the entanglement entropy behaves as in rational conformal field theories (CFT)s. On the other hand, for completely chaotic dual-unitary circuits, it behaves as in holographic CFTs, exhibiting a longer entanglement barrier that drops rapidly when the subsystem thermalizes. Interestingly, the entanglement spectrum is nontrivial in the completely chaotic case. Higher R\'enyi entropies behave in an increasingly similar way to rational CFTs, such that the free and completely chaotic barriers are identical in the limit of infinite replicas (i.e., for the so called min-entropy). We also show that, upon increasing the bond dimension of the MPSs, the barrier maintains the same shape. It simply shifts to the left to accommodate for the larger initial entanglement.