“…Tensor-Network methods.-Tensor-Network (TN) methods [40][41][42][43][44][45] provide efficient representations for ground states of local Hamiltonians based on their entanglement structure [40]. In particular, the projected entangled-pair state (PEPS) method and its infinite version in the thermodynamic limit (iPEPS) [40,41,45,46], have played a major role in the characterization and discovery of many exotic phases, ranging from magnetically ordered states [47][48][49] to QSLs [19,21,50,51] and valence-bond crystals [52,53]. In particular, the modified version of the iPEPS algorithm, designed for the honeycomb structures [54], has been shown to be very successful for simulating and characterizing the Kitaev model and its variants, such as the Kitaev-Heisenberg model, in the thermodynamic limit [55][56][57].…”