Search

When subject to a non-local unitary evolution, qubits in a quantum circuit become increasingly entangled. Conversely, measurements applied to individual qubits lead to their disentanglement from the collective system. The extent of entanglement reduction depends on the frequency of local projective measurements. A delicate balance emerges between unitary evolution, which enhances entanglement, and measurements which diminish it. In the thermodynamic limit, there is a phase transition from volume law entanglement to area law entanglement at a critical value of measurement frequency. This phenomenon, occurring in hybrid quantum circuits with both unitary gates and measurements, is termed as measurement-induced phase transition (MIPT). We study the behavior of MIPT in circuits comprising of two qubit unitary gates parameterized by Cartan decomposition. We show that the entangling power and gate typicality of the two-qubit local unitaries employed in the circuit can be used to explain the behavior of global bipartite entanglement the circuit can sustain. When the two qubit gate throughout the circuit is the identity and measurements are the sole driver of the entanglement behavior, we obtain analytical estimate for the entanglement entropy that shows remarkable agreement with numerical simulations. We also find that the entangling power and gate typicality enable the classification of the two-qubit unitaries by different universality classes of phase transitions that can occur in the hybrid circuit. For all unitaries in a particular universality class, the transition from volume to area law of entanglement occurs with same exponent that characterizes the phase transition., Comment: 14 page, 14 figures

Quantum channels describe subsystem or open system evolution. Using the classical Koopman operator that evolves functions on phase space, 4 classical Koopman channels are identified that are analogs of the 4 possible quantum channels in a bipartite setting. Thus when the complete evolution has a quantum-classical correspondence the correspondence at the level of the subunitary channels can be studied. The channels, both classical and quantum can be interpreted as noisy single particle systems. Having parallel classical and quantum operators gives us new access to study fine details of these major limiting theories. Using a coupled kicked rotor as a generic example, we contrast and compare spectra of the quantum and classical channel. The largest nontrivial mode of the quantum channel is seen to be mostly determined by the stable parts of the classical phase space, even those that are surprisingly small in relation to the scale of an effective $\hbar$. In the case when the dynamics has a significant fraction of chaos the spectrum has a prominent annular density that is approximately described by the single-ring theorem of random matrix theory, and the ring shrinks in size when the classical limit is approached. However, the eigenvalues and modes that survive the classical limit seem to be either scarred by unstable manifolds or, if they exist, stable periodic orbits., Comment: 47 pages, 12 figures

The effects of disorder and chaos on quantum many-body systems can be superficially similar, yet their interplay has not been sufficiently explored. This work finds a continuous phase transition when disorder breaks permutation symmetry, with details of the transition being controlled by the degree of chaos in the clean limit. The system changes from an area law entangled phase in the permutation symmetric subspace where collective variables exist to volume law entanglement in the full Hilbert space, beyond a critical strength of the disorder. The critical strength tends to zero when the original disorder free system is fully chaotic. We study this mainly via the scaling of the collective spin of non-equilibrium states which transit to have properties of what has been dubbed "deep Hilbert space". This has potential implications for general many body physics, as well as technologies such as transmon qubits.

Dual-unitary circuits are being vigorously studied as models of many-body quantum chaos that can be solved exactly for correlation functions and time evolution of states. Here we define their classical counterparts as dual-canonical transformations and associated dual-Koopman operators. Like their quantum counterparts, the correlations vanish everywhere except on the light cone, on which they decay with rates governed by a simple contractive map. Providing a large class of such dual-canonical transformations, we study in detail the example of a coupled standard map and show analytically that arbitrarily away from the integrable case, in the thermodynamic limit the system is mixing. We also define ``perfect" Koopman operators that lead to the correlation vanishing everywhere including on the light cone and provide an example of a cat-map lattice which would qualify to be a Bernoulli system at the apex of the ergodic hierarchy., Comment: Revised version, to be published in Phys. Rev. Lett

Non-KAM (Kolmogorov-Arnold-Moser) systems, when perturbed by weak time-dependent fields, offer a fast route to classical chaos through an abrupt breaking of invariant phase space tori. In this work, we employ out-of-time-order correlators (OTOCs) to study the dynamical sensitivity of a perturbed non-KAM system in the quantum limit as the parameter that characterizes the $\textit{resonance}$ condition is slowly varied. For this purpose, we consider a quantized kicked harmonic oscillator (KHO) model, which displays stochastic webs resembling Arnold's diffusion that facilitate large-scale diffusion in the phase space. Although the Lyapunov exponent of the KHO at resonances remains close to zero in the weak perturbative regime, making the system weakly chaotic in the conventional sense, the classical phase space undergoes significant structural changes. Motivated by this, we study the OTOCs when the system is in resonance and contrast the results with the non-resonant case. At resonances, we observe that the long-time dynamics of the OTOCs are sensitive to these structural changes, where they grow quadratically as opposed to linear or stagnant growth at non-resonances. On the other hand, our findings suggest that the short-time dynamics remain relatively more stable and show the exponential growth found in the literature for unstable fixed points. The numerical results are backed by analytical expressions derived for a few special cases. We will then extend our findings concerning the non-resonant cases to a broad class of near-integrable KAM systems., Comment: 17 pages, 6 figures. Close to the accepted version in Phys. Rev. E

Ordering and classifying multipartite quantum states by their entanglement content remains an open problem. One class of highly entangled states, useful in quantum information protocols, the absolutely maximally entangled (AME) ones, are specially hard to compare as all their subsystems are maximally random. While, it is well-known that there is no AME state of four qubits, many analytical examples and numerically generated ensembles of four qutrit AME states are known. However, we prove the surprising result that there is truly only {\em one} AME state of four qutrits up to local unitary equivalence. In contrast, for larger local dimensions, the number of local unitary classes of AME states is shown to be infinite. Of special interest is the case of local dimension 6 where it was established recently that a four-party AME state does exist, providing a quantum solution to the classically impossible Euler problem of 36 officers. Based on this, an infinity of quantum solutions are constructed and we prove that these are not equivalent. The methods developed can be usefully generalized to multipartite states of any number of particles., Comment: Rewritten as a regular article and few changes in the title from first version. Close to the published version

TThe problem of finding the resource free, closest local unitary, to any bipartite unitary gate $U$ is addressed. Previously discussed as a measure of nonlocality, the distance $K_D(U)$ to the nearest product unitary has implications for circuit complexity and related quantities. Dual unitaries, currently of great interest in models of complex quantum many-body systems, are shown to have a preferred role as these are maximally and equally away from the set of local unitaries. This is proved here for the case of qubits and we present strong numerical and analytical evidence that it is true in general. An analytical evaluation of $K_D(U)$ is presented for general two-qubit gates. For arbitrary local dimensions, that $K_D(U)$ is largest for dual unitaries, is substantiated by its analytical evaluations for an important family of dual-unitary and for certain non-dual gates. A closely allied result concerns, for any bipartite unitary, the existence of a pair of maximally entangled states that it connects. We give efficient numerical algorithms to find such states and to find $K_D(U)$ in general., Comment: 9+2 pages, 5 Figures. Many parts from previous version are rearranged and the current version is rewritten as a regular article

While quantum circuits built from two-particle dual-unitary (maximally entangled) operators serve as minimal models of typically nonintegrable many-body systems, the construction and characterization of dual-unitary operators themselves are only partially understood. A nonlinear map on the space of unitary operators was proposed in PRL.~125, 070501 (2020) that results in operators being arbitrarily close to dual unitaries. Here we study the map analytically for the two-qubit case describing the basins of attraction, fixed points, and rates of approach to dual unitaries. A subset of dual-unitary operators having maximum entangling power are 2-unitary operators or perfect tensors, and are equivalent to four-party absolutely maximally entangled states. It is known that they only exist if the local dimension is larger than $d=2$. We use the nonlinear map, and introduce stochastic variants of it, to construct explicit examples of new dual and 2-unitary operators. A necessary criterion for their local unitary equivalence to distinguish classes is also introduced and used to display various concrete results and a conjecture in $d=3$. It is known that orthogonal Latin squares provide a ``classical combinatorial design" for constructing permutations that are 2-unitary. We extend the underlying design from classical to genuine quantum ones for general dual-unitary operators and give an example of what might be the smallest sized genuinely quantum design of a 2-unitary in $d=4$., Comment: 26+4 pages, 13+2 Figures. Main changes: Proof of Theorem 1 modified, Theorem 2 concerning two-qubit gates is added, and summary of main results included. Version accepted for publication in PRX Quantum

A bipartite system whose subsystems are fully quantum chaotic and coupled by a perturbative interaction with a tunable strength is a paradigmatic model for investigating how isolated quantum systems relax towards an equilibrium. It is found that quantum coherence of the initial product states in the uncoupled eigenbasis can be viewed as a resource for equilibration and approach to thermalization as manifested by the entanglement. Results are given for four distinct perturbation strength regimes, the ultra-weak, weak, intermediate, and strong regimes. For each, three types of initially unentangled states are considered, coherent random-phase superpositions, random superpositions, and eigenstate products. A universal time scale is identified involving the interaction strength parameter. Maximally coherent initial states thermalize for any perturbation strength in spite of the fact that in the ultra-weak perturbative regime the underlying eigenstates of the system have a tensor product structure and are not at all thermal-like; though the time taken to thermalize tends to infinity as the interaction vanishes. In contrast to the widespread linear behavior, in this regime the entanglement initially grows quadratically in time., Comment: 22 Pages, 13 figures (24 including subfigures)

The famous combinatorial problem of Euler concerns an arrangement of $36$ officers from six different regiments in a $6 \times 6$ square array. Each regiment consists of six officers each belonging to one of six ranks. The problem, originating from Saint Petersburg, requires that each row and each column of the array contains only one officer of a given rank and given regiment. Euler observed that such a configuration does not exist. In recent work, we constructed a solution to a quantum version of this problem assuming that the officers correspond to quantum states and can be entangled. In this paper, we explain the solution which is based on a partition of 36 officers into nine groups, each with four elements. The corresponding quantum states are locally equivalent to maximally entangled two-qubit states, hence each officer is entangled with at most three out of his $35$ colleagues. The entire quantum combinatorial design involves $9$ Bell bases in nine complementary $4$-dimensional subspaces., Comment: 26 pages

Out-of-time-order correlators (OTOC) in the Ising Floquet system, that can be both integrable and nonintegrable is studied. Instead of localized spin observables, we study contiguous symmetric blocks of spins or random operators localized on these blocks as observables. We find only power-law growth of OTOC in both integrable and nonintegrable regimes. In the non-integrable regime, beyond the scrambling time, there is an exponential saturation of the OTOC to values consistent with random matrix theory. This motivates the use of "pre-scrambled" random block operators as observables. A pure exponential saturation of OTOC in both integrable and nonintegrable system is observed, without a scrambling phase. Averaging over random observables from the Gaussian unitary ensemble, the OTOC is found to be exactly same as the operator entanglement entropy, whose exponential saturation has been observed in previous studies of such spin-chains., Comment: 13 pages, 6 figures

The negative solution to the famous problem of $36$ officers of Euler implies that there are no two orthogonal Latin squares of order six. We show that the problem has a solution, provided the officers are entangled, and construct orthogonal quantum Latin squares of this size. As a consequence, we find an example of the long-elusive Absolutely Maximally Entangled state AME$(4,6)$ of four subsystems with six levels each, equivalently a $2$-unitary matrix of size $36$, which maximizes the entangling power among all bipartite unitary gates of this dimension, or a perfect tensor with four indices, each running from one to six. This special state deserves the appellation golden AME state as the golden ratio appears prominently in its elements. This result allows us to construct a pure nonadditive quhex quantum error detection code $(\!(3,6,2)\!)_6$, which saturates the Singleton bound and allows one to encode a $6$-level state into a triplet of such states., Comment: 14 pages, 12 figures

Deterministic classical dynamical systems have an ergodic hierarchy, from ergodic through mixing, to Bernoulli systems that are "as random as a coin-toss". Dual-unitary circuits have been recently introduced as solvable models of many-body nonintegrable quantum chaotic systems having a hierarchy of ergodic properties. We extend this to include the apex of a putative quantum ergodic hierarchy which is Bernoulli, in the sense that correlations of single and two-particle observables vanish at space-time separated points. We derive a condition based on the entangling power $e_p(U)$ of the basic two-particle unitary building block, $U$, of the circuit, that guarantees mixing, and when maximized, corresponds to Bernoulli circuits. Additionally we show, both analytically and numerically, how local-averaging over random realizations of the single-particle unitaries, $u_i$ and $v_i$ such that the building block is $U^\prime = (u_1 \otimes u_2 ) U (v_1 \otimes v_2 )$ leads to an identification of the average mixing rate as being determined predominantly by the entangling power $e_p(U)$. Finally we provide several, both analytical and numerical, ways to construct dual-unitary operators covering the entire possible range of entangling power. We construct a coupled quantum cat map which is dual-unitary for all local dimensions and a 2-unitary or perfect tensor for odd local dimensions, and can be used to build Bernoulli circuits., Comment: 48 pages (38 pages main text + 10 pages Appendix ), 16 (10+6) figures. Contains new results and substantially modified

The out-of-time-ordered correlators (OTOC) and the Loschmidt echo are two measures that are now widely being explored to characterize sensitivity to perturbations and information scrambling in complex quantum systems. Studying few qubits systems collectively modelled as a kicked top, we solve exactly the three- and four- qubit cases, giving analytical results for the OTOC and the Loschmidt echo. While we may not expect such few-body systems to display semiclassical features, we find that there are clear signatures of the exponential growth of OTOC even in systems with as low as 4 qubits in appropriate regimes, paving way for possible experimental measurements. We explain qualitatively how classical phase space structures like fixed points and periodic orbits have an influence on these quantities and how our results compare to the large-spin kicked top model. Finally we point to a peculiar case at the border of quantum-classical correspondence which is solvable for any number of qubits and yet has signatures of exponential sensitivity in a rudimentary form., Comment: 14 pages, 15 figures

Coined discrete-time quantum walks are studied using simple deterministic dynamical systems as coins whose classical limit can range from being integrable to chaotic. It is shown that a Loschmidt echo like fidelity plays a central role and when the coin is chaotic this is approximately the characteristic function of a classical random walker. Thus the classical binomial distribution arises as a limit of the quantum walk and the walker exhibits diffusive growth before eventually becoming ballistic. The coin-walker entanglement growth is shown to be logarithmic in time as in the case of many-body localization and coupled kicked rotors, and saturates to a value that depends on the relative coin and walker space dimensions. In a coin dominated scenario, the chaos can thermalize the quantum walk to typical random states such that the entanglement saturates at the Haar averaged Page value, unlike in a walker dominated case when atypical states seem to be produced.

The extreme-value statistics of the entanglement spectrum in disordered spin chains possessing a many-body localization transition is examined. It is expected that eigenstates in the metallic or ergodic phase, behave as random states and hence the eigenvalues of the reduced density matrix, commonly referred to as the entanglement spectrum, are expected to follow the eigenvalue statistics of a trace normalized Wishart ensemble. In particular, the density of eigenvalues is supposed to follow the universal Marchenko-Pastur distribution. We find deviations in the tails both for the disordered XXZ with total $S_z$ conserved in the half-filled sector as well as in a model that breaks this conservation. A sensitive measure of deviations is provided by the largest eigenvalue, which in the case of the Wishart ensemble after appropriate shift and scaling follows the universal Tracy-Widom distribution. We show that for the models considered, in the metallic phase, the largest eigenvalue of the reduced density matrix of eigenvector, instead follows the generalized extreme-value statistics bordering on the Fisher-Tipett-Gumbel distribution indicating that the correlations between eigenvalues are much weaker compared to the Wishart ensemble. We show by means of distributions conditional on the high entropy and normalized participation ratio of eigenstates that the conditional entanglement spectrum still follows generalized extreme value distribution. In the deeply localized phase we find heavy tailed distributions and L\'evy stable laws in an appropriately scaled function of the largest and second largest eigenvalues. The scaling is motivated by a recently developed perturbation theory of weakly coupled chaotic systems.

Maximally entangled bipartite unitary operators or gates find various applications from quantum information to being building blocks of minimal models of many-body quantum chaos, and have been referred to as "dual unitaries". Dual unitary operators that can create the maximum average entanglement when acting on product states have to satisfy additional constraints. These have been called "2-unitaries" and are examples of perfect tensors that can be used to construct absolutely maximally entangled states of four parties. Hitherto, no systematic method exists, in any local dimension, which result in the formation of such special classes of unitary operators. We outline an iterative protocol, a nonlinear map on the space of unitary operators, that creates ensembles whose members are arbitrarily close to being dual unitaries, while for qutrits and ququads we find that a slightly modified protocol yields a plethora of 2-unitaries. We further characterize the dual unitary operators via their entangling power and the 2-unitaries via the distribution of entanglement created from unentangled states., Comment: 6+3 pages. 4 figures in the main text + 5 figures in supplementary material. Version accepted for publication in Phys. Rev. Lett

We study the quantized kicked top's dynamics under projective measurements especially with an aim to connect macroscopic realism (macrorealism) and chaos. The kicked top is a classically chaotic system, and macrorealism is a set of assumptions regarding how systems should behave according to classical intuition which is at present being tested in experiments. Using the no-signalling in time condition, a derivative of macrorealism, we define measures of disturbance due to measurements in order to study the effects of chaos. Restricting to cases which allow a meaningful study of this question, our results strongly suggest that chaos is unfriendly to macrorealism., Comment: On macroscopic realism vs. chaos. 18 pages, 14 figures

Out-of-time-order correlators (OTOC) being explored as a measure of quantum chaos, is studied here in a coupled bipartite system. Each of the subsystems can be chaotic or regular and lead to very different OTOC growths both before and after the scrambling or Ehrenfest time. We present preliminary results on weakly coupled subsystems which have very different Lyapunov exponents. We also review the case when both the subsystems are strongly chaotic when a random matrix model can be pressed into service to derive an exponential relaxation to saturation., Comment: 18 pages, 6 figures, version accepted in Acta Physica Polonica A

Entanglement properties of bipartite unitary operators are studied via their local invariants, namely the entangling power $e_p$ and a complementary quantity, the gate typicality $g_t$. We characterize the boundaries of the set $K_2$ representing all two-qubit gates projected onto the plane $(e_p, g_t)$ showing that the fractional powers of the \textsc{swap} operator form a parabolic boundary of $K_2$, while the other bounds are formed by two straight lines. In this way a family of gates with extreme properties is identified and analyzed. We also show that the parabolic curve representing powers of \textsc{swap} persists in the set $K_N$, for gates of higher dimensions ($N>2$). Furthermore, we study entanglement of bipartite quantum gates applied sequentially $n$ times and analyze the influence of interlacing local unitary operations, which model generic Hamiltonian dynamics. An explicit formula for the entangling power a gate applied $n$ times averaged over random local unitary dynamics is derived for an arbitrary dimension of each subsystem. This quantity shows an exponential saturation to the value predicted by the random matrix theory (RMT), indicating "thermalization" in the entanglement properties of sequentially applied quantum gates that can have arbitrarily small, but nonzero, entanglement to begin with. The thermalization is further characterized by the spectral properties of the reshuffled and partially transposed unitary matrices., Comment: 21 pages, 8 figures

The entanglement production in bipartite quantum systems is studied for initially unentangled product eigenstates of the subsystems, which are assumed to be quantum chaotic. Based on a perturbative computation of the Schmidt eigenvalues of the reduced density matrix, explicit expressions for the time-dependence of entanglement entropies, including the von Neumann entropy, are given. An appropriate re-scaling of time and the entropies by their saturation values leads a universal curve, independent of the interaction. The extension to the non-perturbative regime is performed using a recursively embedded perturbation theory to produce the full transition and the saturation values. The analytical results are found to be in good agreement with numerical results for random matrix computations and a dynamical system given by a pair of coupled kicked rotors., Comment: 18 pages, 16 figures , custom bibstyle file

Out-of-time-order correlators (OTOC), vigorously being explored as a measure of quantum chaos and information scrambling, is studied here in the natural and simplest multi-particle context of bipartite systems. We show that two strongly chaotic and weakly interacting subsystems display two distinct phases in the growth of OTOC.The first is dominated by intra-subsystem scrambling, when an exponential growth with a positive Lyapunov exponent is observed till the Ehrenfest time. This phase is essentially independent of the interaction, while the second phase is an interaction dominated exponential approach to saturation that is universal and described by a random matrix model. This simple random matrix model of weakly interacting strongly chaotic bipartite systems, previously employed for studying entanglement and spectral transitions, is approximately analytically solvable for its OTOC. The example of two coupled kicked rotors is used to demonstrate the two phases, and the extent to which the random matrix model is applicable. That the two phases correspond to delocalization in the subsystems followed by inter-subsystem mixing is seen via the participation ratio in phase-space. We also point out that the second, universal, phase alone exists when the observables are in a sense already scrambled. Thus, while the post-Ehrenfest time OTOC growth is in general not well-understood, the case of strongly chaotic and weakly coupled systems presents an, perhaps important, exception., Comment: 11 pages, 6 figures

Out-of time-ordered correlators (OTOC) have recently attracted significant attention from the physics of many-body systems, to quantum black-holes, with an exponential growth of the OTOC indicating quantum chaos. Here we consider OTOC in the context of coined discrete quantum walks, a very well studied model of quantization of classical random walks with applications to quantum algorithms. Three separate cases of operators, variously localized in the coin and walker spaces, are discussed in this context and it is found that the approximated behavior of the OTOCs is well described by simple algebraic functions in all these three cases with different time scale of growth. The quadratic increase of OTOC signals the absence of quantum chaos in these simplest forms of quantum walks., Comment: 19 pages, 3 figures

Exactly solvable models that exhibit quantum signatures of classical chaos are both rare as well as important - more so in view of the fact that the mechanisms for ergodic behavior and thermalization in isolated quantum systems and its connections to non-integrability are under active investigation. In this work, we study quantum systems of few qubits collectively modeled as a kicked top, a textbook example of quantum chaos. In particular, we show that the 3 and 4 qubit cases are exactly solvable and yet, interestingly, can display signatures of ergodicity and thermalization. Deriving analytical expressions for entanglement entropy and concurrence, we see agreement in certain parameter regimes between long-time average values and ensemble averages of random states with permutation symmetry. Comparing with results using the data of a recent transmons based experiment realizing the 3-qubit case, we find agreement for short times, including a peculiar step-like behaviour in correlations of some states. In the case of 4-qubits we point to a precursor of dynamical tunneling between what in the classical limit would be two stable islands. Numerical results for larger number of qubits show the emergence of the classical limit including signatures of a bifurcation., Comment: supersedes arXiv:1808.07741. More detailed analysis, new results, more figures

The out-of-time-ordered correlator (OTOC) is a measure of quantum chaos that is being vigorously investigated. Analytically accessible simple models that have long been studied in other contexts could provide insights into such measures. This paper investigates the OTOC in the quantum bakers map which is the quantum version of a simple and exactly solvable model of deterministic chaos that caricatures the action of kneading dough. Exact solutions based on the semiquantum approximation are derived that tracks very well the correlators till the Ehrenfest time. The growth occurs at the exponential rate of the classical Lyapunov exponent, but modulated by slowly changing coefficients. Beyond this time saturation occurs as a value close to that of random matrices. Using projectors for observables naturally leads to truncations of the unitary time-$t$ propagator and the growth of their singular values is shown to be intimately related to the growth of the out-of-time-ordered correlators., Comment: 20 pages, 6 figures

Spectral statistics of quantum systems have been studied in detail using the nearest neighbour level spacings, which for generic chaotic systems follows random matrix theory predictions. In this work, the probability density of the closest neighbour and farther neighbour spacings from a given level are introduced. Analytical predictions are derived using a $3 \times 3$ matrix model. The closest neighbour density is generalized to the $k-$th closest neighbour spacing density, which allows for investigating long-range correlations. For larger $k$ the probability density of $k-$th closest neighbour spacings is well described by a Gaussian. Using these $k-$th closest neighbour spacings we propose the ratio of the closest neighbour to the second closest neighbour as an alternative to the ratio of successive spacings. For a Poissonian spectrum the density of the ratio is flat, whereas for the three Gaussian ensembles repulsion at small values is found. The ordered spacing statistics and their ratio are numerically studied for the integrable circle billiard, the chaotic cardioid billiard, the standard map and the zeroes of the Riemann zeta function. Very good agreement with the predictions is found.

Exactly solvable models that exhibit quantum signatures of classical chaos are both rare as well as important - more so in view of the fact that the mechanisms for ergodic behavior and thermalization in isolated quantum systems and its connections to non-integrability are under active investigation. In this work, we study quantum systems of few qubits collectively modeled as a kicked top, a textbook example of quantum chaos. In particular, we show that the 3 and 4 qubit cases are exactly solvable and yet, interestingly, can display signatures of ergodicity and thermalization. Deriving analytical expressions for entanglement entropy and concurrence, we see agreement in certain parameter regimes between long-time average values and ensemble averages of random states with permutation symmetry. Comparing with results using the data of a recent transmons based experiment realizing the 3-qubit case, we find agreement for short times, including a peculiar step-like behaviour in correlations of some states. In the case of 4-qubits we point to a precursor of dynamical tunneling between what in the classical limit would be two stable islands. Numerical results for larger number of qubits show the emergence of the classical limit including signatures of a bifurcation., Comment: 6+5 pages, 7 figures

We derive universal entanglement entropy and Schmidt eigenvalue behaviors for the eigenstates of two quantum chaotic systems coupled with a weak interaction. The progression from a lack of entanglement in the noninteracting limit to the entanglement expected of fully randomized states in the opposite limit is governed by the single scaling transition parameter, $\Lambda$. The behaviors apply equally well to few- and many-body systems, e.g.\ interacting particles in quantum dots, spin chains, coupled quantum maps, and Floquet systems as long as their subsystems are quantum chaotic, and not localized in some manner. To calculate the generalized moments of the Schmidt eigenvalues in the perturbative regime, a regularized theory is applied, whose leading order behaviors depend on $\sqrt{\Lambda}$. The marginal case of the $1/2$ moment, which is related to the distance to closest maximally entangled state, is an exception having a $\sqrt{\Lambda}\ln \Lambda$ leading order and a logarithmic dependence on subsystem size. A recursive embedding of the regularized perturbation theory gives a simple exponential behavior for the von Neumann entropy and the Havrda-Charv{\' a}t-Tsallis entropies for increasing interaction strength, demonstrating a universal transition to nearly maximal entanglement. Moreover, the full probability densities of the Schmidt eigenvalues, i.e.\ the entanglement spectrum, show a transition from power laws and L\'evy distribution in the weakly interacting regime to random matrix results for the strongly interacting regime. The predicted behaviors are tested on a pair of weakly interacting kicked rotors, which follow the universal behaviors extremely well.

Many-body states that are invariant under particle relabelling, the permutation symmetric states, occur naturally when the system dynamics is described by symmetric processes or collective spin operators. We derive expressions for the reduced density matrix for arbitrary subsystem decomposition for these states and study properties of permutation symmetric states and their subsystems when the joint system is picked randomly and uniformly. Thus defining a new random matrix ensemble, we find the average linear entropy and von Neumann entropy which implies that random permutation symmetric states are marginally entangled and as a consequence the tripartite mutual information (TMI) is typically positive, preventing information from being shared globally. Applying these results to the quantum kicked top viewed as a multi-qubit system we find that entanglement, mutual information and TMI all increase for large subsystems across the Ehrenfest or log-time and saturate at the random state values if there is global chaos. During this time the out-of-time order correlators (OTOC) evolve exponentially implying scrambling in phase space. We discuss how positive TMI may coexist with such scrambling., Comment: Minor modifications to the text, references added; published in Phys. Rev. E 98, 052205, doi: 10.1103/PhysRevE.98.052205

The entangling power and operator entanglement entropy are state independent measures of entanglement. Their growth and saturation is examined in the time-evolution operator of quantum many-body systems that can range from the integrable to the fully chaotic. An analytically solvable integrable model of the kicked transverse field Ising chain is shown to have ballistic growth of operator von Neumann entanglement entropy and exponentially fast saturation of the linear entropy with time. Surprisingly a fully chaotic model with longitudinal fields turned on shares the same growth phase, and is consistent with a random matrix model that is also exactly solvable for the linear entropy entanglements. However an examination of the entangling power shows that its largest value is significantly less than the nearly maximal value attained by the nonintegrable one. The importance of long-range spectral correlations, and not just the nearest neighbor spacing, is pointed out in determing the growth of entanglement in nonintegrable systems. Finally an interesting case that displays some features peculiar to both integrable and nonintegrable systems is briefly discussed., Comment: Comments are welcome

A quantum walk on a toral phase space involving translations in position and its conjugate momentum is studied in the simple context of a coined walker in discrete time. The resultant walk, with a family of coins parametrized by an angle is such that its spectrum is exactly solvable with eigenangles for odd parity lattices being equally spaced, a feature that is remarkably independent of the coin. The eigenvectors are naturally specified in terms the $q-$Pochhammer symbol, but can also be written in terms of elementary functions, and their entanglement can be analytically found. While the phase space walker shares many features in common with the well-studied case of a coined walker in discrete time and space, such as ballistic growth of the walker position, it also presents novel features such as exact periodicity, and formation of cat-states in phase-space. Participation ratio (PR) a measure of delocalization in walker space is studied in the context of both kinds of quantum walks; while the classical PR increases as $\sqrt{t}$ there is a time interval during which the quantum walks display a power-law growth $\sim t^{0.825}$. Studying the evolution of coherent states in phase space under the walk enables us to identify an Ehrenfest time after which the coin-walker entanglement saturates.

Long-time average behavior of quantum correlations in a multi-qubit system, collectively modeled as a kicked top, is addressed here. The behavior of dynamical generation of quantum correlations such as entanglement, discord, concurrence, as previously studied, and Bell correlation function and tangle, as identified in this study, is a function of initially localized coherent states. Their long-time average reproduces coarse-grained classical phase space structures of the kicked top which contrast, often starkly, chaotic and regular regions. Apart from providing numerical evidence of such correspondence in the semiclassical regime of a large number of qubits, we use data from a recent transmons based experiment to explore this in the deep quantum regime of a 3-qubit kicked top. The degree to which quantum correlations can be regarded as a quantum signature of chaos, and in what ways the various correlation measures are similar or distinct are discussed., Comment: 6 pages, 4 figures

It is demonstrated here that local dynamics have the ability to strongly modify the entangling power of unitary quantum gates acting on a composite system. The scenario is common to numerous physical systems, in which the time evolution involves local operators and nonlocal interactions. To distinguish between distinct classes of gates with zero entangling power we introduce a complementary quantity called gate-typicality and study its properties. Analyzing multiple applications of any entangling operator interlaced with random local gates, we prove that both investigated quantities approach their asymptotic values in a simple exponential form. This rapid convergence to equilibrium, valid for subsystems of arbitrary size, is illustrated by studying multiple actions of diagonal unitary gates and controlled unitary gates., Comment: 7 pages, 3 figures

The entanglement and localization in eigenstates of strongly chaotic subsystems are studied as a function of their interaction strength. Excellent measures for this purpose are the von-Neumann entropy, Havrda-Charv{\' a}t-Tsallis entropies, and the averaged inverse participation ratio. All the entropies are shown to follow a remarkably simple exponential form, which describes a universal and rapid transition to nearly maximal entanglement for increasing interaction strength. An unexpectedly exact relationship between the subsystem averaged inverse participation ratio and purity is derived that infers the transition in the localization as well.

Local entanglement between pairs of spins, as measured by concurrence, is investigated in a disordered spin model that displays a transition from an ergodic to a many-body localized phase in excited states. It is shown that the concurrence vanishes in the ergodic phase and becomes nonzero and increases in the many-body localized phase. This happen to be correlated with the transition in the spectral statistics from Wigner to Poissonian distribution. A scaling form is found to exist in the second derivative of the concurrence with the disorder strength. It also displays a critical value for the localization transition that is close to what is known in the literature from other measures. An exponential decay of concurrence with distance between spins is observed in the localized phase. Nearest neighbor spin concurrence in this phase is also found to be strongly correlated with the disorder configuration of onsite fields: nearly similar fields implying larger entanglement., Comment: 8 pages, 7 figures

Quantities invariant under local unitary transformations are of natural interest in the study of entanglement. This paper deduces and studies a particularly simple quantity that is constructed from a combination of two standard permutations of the density matrix, namely realignment and partial transpose. This bipartite quantity, denoted here as $R_{12}$, vanishes on large classes of separable states including classical-quantum correlated states, while being maximum for only maximally entangled states. It is shown to be naturally related to the 3-tangle in three qubit states via their two-qubit reduced density matrices. Upper and lower bounds on concurrence and negativity of two-qubit density matrices for all ranks are given in terms of $R_{12}$. Ansatz states satisfying these bounds are given and verified using various numerical methods. In rank-2 case it is shown that the states satisfying the lower bound on $R_{12}$ {\it vs} concurrence define a class of three qubit states that maximizes the tripartite entanglement (the 3-tangle) given an amount of entanglement between a pair of them. The measure $R_{12}$ is conjectured, via numerical sampling, to be always larger than the concurrence and negativity. In particular this is shown to be true for the physically interesting case of $X$ states., Comment: Substantially improved version (now 16 pages, 11 figures) that is accepted for publication in Phys. Rev. A

The spectral properties of interacting strongly chaotic systems are investigated for growing interaction strength. A very sensitive transition from Poisson statistics to that of random matrix theory is found. We introduce a new random matrix ensemble modeling this dynamical symmetry breaking transition which turns out to be universal and depends on a single scaling parameter only. Coupled kicked rotors, a dynamical systems paradigm for such transitions, are compared with this ensemble and excellent agreement is found for the nearest-neighbor-spacing distribution. It turns out that this transition is described quite accurately using perturbation theory., Comment: 5 pages, 3 figures

We study the properties of the eigenvalues of real random matrices and their products. It is known that when the matrix elements are Gaussian-distributed independent random variables, the fraction of real eigenvalues tends to unity as the number of matrices in the product increases. Here we present numerical evidence that this phenomenon is robust with respect to the probability distribution of matrix elements, and is therefore a general property that merits detailed investigation. Since the elements of the product matrix are no longer distributed as those of the single matrix nor they remain independent random variables, we study the role of these two factors in detail. We study numerically the properties of the Hadamard (or Schur) product of matrices and also the product of matrices whose entries are independent but have the same marginal distribution as that of normal products of matrices, and find that under repeated multiplication, the probability of all eigenvalues to be real increases in both cases, but saturates to a constant below unity showing that the correlations amongst the matrix elements are responsible for the approach to one. To investigate the role of the non-normal nature of the probability distributions, we present a thorough analytical treatment of the $2 \times 2$ single matrix for several standard distributions. Within the class of smooth distributions with zero mean and finite variance, our results indicate that the Gaussian distribution has the maximum probability of real eigenvalues, but the Cauchy distribution characterised by infinite variance is found to have a larger probability of real eigenvalues than the normal. We also find that for the two-dimensional single matrices, the probability of real eigenvalues lies in the range [5/8,7/8]., Comment: To appear in J. Phys. A: Math, Theor

Record statistics is the study of how new highs or lows are created and sustained in any dynamical process. The study of the highest or lowest records constitute the study of extreme values. This paper represents an exploration of record statistics for certain aspects of the classical and quantum standard map. For instance the momentum square or energy records is shown to behave like that of records in random walks when the classical standard map is in a regime of hard chaos. However different power laws is observed for the mixed phase space regimes. The presence of accelerator modes are well-known to create anomalous diffusion and we notice here that the record statistics is very sensitive to their presence. We also discuss records in random vectors and use it to analyze the {\it quantum} standard map via records in their eigenfunction intensities, reviewing some recent results along the way.

We present a solvable model of iterating cluster state protocols that lead to entanglement production, between contiguous blocks, of 1 ebit per iteration. This continues till the blocks are maximally entangled at which stage an unravelling begins at the same rate till the blocks are unentangled. The model is a variant of the transverse field Ising model and can be implemented with CNOT and single qubit gates. The inter qubit entanglement as measured by the concurrence is shown to be zero for periodic chain realizations while for open boundaries there are very specific instances at which these can develop. Thus we introduce a class of simply produced states with very large multipartite entanglement content of potential use in measurement based quantum computing., Comment: 7 pages, 3 figures

The eigenstates of a quantum spin glass Hamiltonian with long-range interaction are examined from the point of view of localisation and entanglement. In particular, low particle sectors are examined and an anomalous family of eigenstates is found that is more delocalised but also has larger inter-spin entanglement. These are then identified as particle-added eigenstates from the one-particle sector. This motivates the introduction and the study of random promoted two-particle states, and it is shown that they may have large delocalisation such as generic ran- dom states and scale exactly like them. However, the entanglement as measured by two-spin concurrence displays different scaling with the total number of spins. This shows how for different classes of complex quantum states entanglement can be qualitatively different even if localisation measures such as participation ratio are not., Comment: 7 pages, 3 figures, 1 table

Nonlocal properties of an ensemble of diagonal random unitary matrices of order $N^2$ are investigated. The average Schmidt strength of such a bipartite diagonal quantum gate is shown to scale as $\log N$, in contrast to the $\log N^2$ behavior characteristic to random unitary gates. Entangling power of a diagonal gate $U$ is related to the von Neumann entropy of an auxiliary quantum state $\rho=AA^{\dagger}/N^2$, where the square matrix $A$ is obtained by reshaping the vector of diagonal elements of $U$ of length $N^2$ into a square matrix of order $N$. This fact provides a motivation to study the ensemble of non-hermitian unimodular matrices $A$, with all entries of the same modulus and random phases and the ensemble of quantum states $\rho$, such that all their diagonal entries are equal to $1/N$. Such a state is contradiagonal with respect to the computational basis, in sense that among all unitary equivalent states it maximizes the entropy copied to the environment due to the coarse graining process. The first four moments of the squared singular values of the unimodular ensemble are derived, based on which we conjecture a connection to a recently studied combinatorial object called the "Borel triangle". This allows us to find exactly the mean von Neumann entropy for random phase density matrices and the average entanglement for the corresponding ensemble of bipartite pure states., Comment: 14 pages, 6 figures

The kicked Ising model with both a pulsed transverse and a continuous longitudinal field is studied numerically. Starting from a large transverse field and a state that is nearly an eigenstate, the pulsed transverse field is quenched with a simultaneous enhancement of the longitudinal field. The generation of multipartite entanglement is observed along with a phenomenon akin to quantum resonance when the entanglement does not evolve for certain values of the pulse duration. Away from the resonance, the longitudinal field can drive the entanglement to near maximum values that is shown to agree well with those of random states. Further evidence is presented that the time evolved states obtained do have some statistical properties of such random states. For contrast the case when the fields have a steady value is also discussed., Comment: 7 pages, 7 figures

The probability that there are $k$ real eigenvalues for an $n$ dimensional real random matrix is known. Here we study this for the case of products of independent random matrices. Relating the problem of the probability that the product of two real 2 dimensional random matrices has real eigenvalues to an issue of optimal quantum entanglement, this is fully analytically solved. It is shown that in $\pi/4$ fraction of such products the eigenvalues are real. Being greater than the corresponding known probability ($1/\sqrt{2}$) for a single matrix, it is shown numerically that the probability that {\it all} eigenvalues of a product of random matrices are real tends to unity as the number of matrices in the product increases indefinitely. Some other numerical explorations, including the expected number of real eigenvalues is also presented, where an exponential approach of the expected number to the dimension of the matrix seems to hold., Comment: Revised and published version, 5 pages, 4 figures

Motivated by link transformations of lattice gauge theory, a method for generating local unitary invariants, especially for a system of qubits, has been pointed out in an earlier work [M. S. Williamson {\it et. al.}, Phys. Rev. A {\bf 83}, 062308 (2011)]. This paper first points the equivalence of the so constructed transformations to the combined operations of partial transpose and realignment. This allows construction of local unitary invariants of any system, with subsystems of arbitrary dimensions. Some properties of the resulting operators and consequences for pure tripartite higher dimensional states are briefly discussed., Comment: 5 pages, 2 figures

The transition from a gapless liquid to a gapped dimerized ground state that occurs in the frustrated antiferromagnetic Majumdar-Ghosh (or $J_1-J_2$ Heisenberg) model is revisited from the point of view of entanglement. We study the evolution of entanglement spectra, a "projected subspace" block entropy, and concurrence in the Schmidt vectors through the transition. The standard tool of Schmidt decomposition along with the existence of the unique MG point where the ground states are degenerate and known exactly, suggests the projection into two orthogonal subspaces that is useful even away from this point. Of these, one is a dominant five dimensional subspace containing the complete state at the MG point and the other contributes marginally, albeit with increasing weight as the number of spins is increased. We find that the marginally contributing subspace has a minimum von Neumann entropy in the vicinity of the dimerization transition. Entanglement content between pairs of spins in the Schmidt vectors, studied via concurrence, shows that those belonging to the dominant five dimensional subspace display a clear progress towards dimerization, with the concurrence vanishing on odd/even sublattices, again in the vicinity of the dimerization, and maximizing in the even/odd sublattices at the MG point. In contrast, study of the Schmidt vectors in the marginally contributing subspace, as well as in the projection of the ground state in this space, display pair concurrence which decrease on both the sublattices as the MG point is approached. The robustness of these observations indicate their possible usefulness in the study of models that have similar transitions, and have hitherto been difficult to study using standard entanglement signatures., Comment: 25 pages, 8 figures

The record statistics of complex random states are analytically calculated, and shown that the probability of a record intensity is a Bernoulli process. The correlation due to normalization leads to a probability distribution of the records that is non-universal but tends to the Gumbel distribution asymptotically. The quantum standard map is used to study these statistics for the effect of correlations apart from normalization. It is seen that in the mixed phase space regime the number of intensity records is a power law in the dimensionality of the state as opposed to the logarithmic growth for random states., Comment: figures redrawn, discussion added