Threefold Way for Typical Entanglement

A typical quantum state with no symmetry can be realized by letting a random unitary act on a fixed state, and the subsystem entanglement spectrum follows the Laguerre unitary ensemble (LUE). For integer-spin time reversal symmetry, we have an analogous scenario where we prepare a time-reversal symmetric state and let random orthogonal matrices act on it, leading to the Laguerre orthogonal ensemble (LOE). However, for half-integer-spin time reversal symmetry, a straightforward analog leading to the Laguerre symplectic ensemble (LSE) is no longer valid because that time-reversal symmetric state is forbidden by the Kramers' theorem. We devise a system in which the global time reversal operator is fractionalized on the subsystems, and show that LSE arises in the system. Extending this idea, we incorporate general symmetry fractionalization into the system, and show that the statistics of the entanglement spectrum is decomposed into a direct sum of LOE, LUE, and/or LSE. Here, various degeneracies in the entanglement spectrum may appear, depending on the non-Abelian nature of the symmetry group and the cohomology class of the nontrivial projective representation on the subsystem. Our work establishes the entanglement counterpart of Dyson's threefold way for Hamiltonians with symmetries.