Fig. 2: Dynamically emerging Nagaoka ferromagnetism in a doped Kitaev spin liquid with FM spin couplings.

a Typical examples for snapshots of the Fock configurations occurring with large probability in the wave function; see Eq. (6); obtained from matrix product states (MPS) at times τ = 0 [1/J] and τ = 10 [1/J] are shown. We identify the FM cluster size NC by connecting all spins pointing up (purple dots) or down (orange dots) starting from the hole site (green dots). b Extracted from 10,000 snapshots of the state at each time, a large ferromagnetic (FM) cluster of spins forms around the hole. Of the N = 160 sites in the system, the cluster has a size of Nc ≈ 35 (40) sites for hopping constants t/J = 3 (6). The speed of the FM cluster expansions is set by the hopping constant t; see dashed lines fitted to the snapshot data. c While the FM cluster size increases, the flux operator average 〈W〉 decays with time. All data shown is for ferromagnetic couplings J > 0 and anisotropy Jz/J = 2.5.