We present thermodynamic phase diagrams showing magnetic analog of ``three states of matter,'' namely, spin liquid, paramagnetic, and magnetically ordered phases, obtained by unbiased quantum Monte Carlo simulations. Our simulations are carried out for Kitaev's toric codes in two and three dimensions, i.e., the simplest realizations of gapped topological
Z2 spin liquids, extended by a nearest-neighbor ferromagnetic Ising coupling. We find that the ordered phase borders on the spin liquid by a discontinuous transition line in three dimensions, while it grows continuously from the quantum critical point in two dimensions. In both cases, our results elucidate peculiar proximity effects to the nearby spin liquids in the high-temperature paramagnetic phase, even when the ground state is magnetically ordered. The thorough study of magnetic three states of matter is achieved by introducing the ``fictitious vertex'' method into the directed loop algorithm. This provides a generic prescription to simulate models with off-diagonal multispin interactions, in which the conventional scheme may suffer from intrinsic ergodicity breakdown as in the present case.
We calculate the Wilson ratio of the one-dimensional Fermi gas with spin imbalance. The Wilson ratio of attractively interacting fermions is solely determined by the density stiffness and sound velocity of pairs and of excess fermions for the two-component Tomonaga-Luttinger liquid phase. The ratio exhibits anomalous enhancement at the two critical points due to the sudden change in the density of states. Despite a breakdown of the quasiparticle description in one dimension, two important features of the Fermi liquid are retained; namely, the specific heat is linearly proportional to temperature, whereas the susceptibility is independent of temperature. In contrast to the phenomenological Tomonaga-Luttinger liquid parameter, the Wilson ratio provides a powerful parameter for testing universal quantum liquids of interacting fermions in one, two, and three dimensions.