We consider the time-dependent thermal pure quantum state method and introduce the efficient scheme to evaluate the change in physical quantities induced by the time-dependent perturbations, which has been proposed in our previous paper [H. Taguchi et al., Phys. Rev. B 105, 125137 (2022)]. Here, we treat the Kitaev model to consider the Majorana-mediated spin transport, as an example. We demonstrate how efficient our scheme is to evaluate spin oscillations induced by the magnetic field pulse.
We consider the time-dependent thermal pure quantum state method and introduce the efficient scheme to evaluate the change in physical quantities induced by the time-dependent perturbations, which has been proposed in our previous paper [H. Taguchi et al., Phys. Rev. B 105, 125137 (2022)]. Here, we treat the Kitaev model to consider the Majorana-mediated spin transport, as an example. We demonstrate how efficient our scheme is to evaluate spin oscillations induced by the magnetic field pulse.
We study how stable the Majorana-mediated spin transport in a quantum spin Kitaev model is against thermal fluctuations. Using the time-dependent thermal pure quantum state method, we examine finite-temperature spin dynamics in the Kitaev model. The model exhibits two characteristic temperatures, TL and TH, which correspond to energy scales of the local flux and the itinerant Majorana fermion, respectively. At low temperatures $(T\ensuremath{\ll}{T}_{L})$, an almost flux-free state is realized, and the spin excitation propagates in a similar way to that for the ground state. Namely, after the magnetic pulse is introduced at one of the edges, the itinerant Majorana fermions propagate the spin excitations even through the quantum spin liquid state region, and oscillations in the spin moment appear in the other edge with a tiny magnetic field. When $T\ensuremath{\sim}{T}_{L}$, larger oscillations in the spin moments are induced in the other edge, compared with the results at the ground state. At higher temperatures, excited Z2 fluxes disturb the coherent motion of the itinerant Majorana fermions, which suppresses the spin propagation. Our results demonstrate a crucial role of thermal fluctuations in the Majorana-mediated spin transport.
We study the spin transport through the quantum spin liquid (QSL) by investigating the real-time and real-space dynamics of the Kitaev spin system with zigzag edges using the time-dependent Majorana mean-field theory. After the magnetic-field pulse is introduced to one of the edges, spin moments are excited in the opposite edge region although spin moments are never induced in the Kitaev QSL region. This unusual spin transport originates from the fact that the S=1/2 spins are fractionalized into the itinerant and localized Majorana fermions in the Kitaev system. Although both Majorana fermions are excited by the magnetic pulse, only the itinerant ones flow through the bulk regime without spin excitations, resulting in the spin transport in the Kitaev system despite the presence of a nonzero spin gap. We also demonstrate that this phenomenon can be observed in the system with small Heisenberg interactions using the exact diagonalization.
The proliferation of quantum fluctuations and long-range entanglement presents an outstanding challenge for the numerical simulation of interacting spin systems with exotic ground states. Here, we present a toolset of Chebyshev polynomial-based iterative methods that provides a unified framework to study the thermodynamical properties, critical behavior and dynamics of frustrated quantum spin models with controlled accuracy. Similar to previous applications of the Chebyshev spectral methods to condensed matter systems, the algorithmic complexity scales linearly with the Hilbert space dimension and the Chebyshev truncation order. Using this approach, we study two paradigmatic quantum spin models on the honeycomb lattice: the Kitaev-Heisenberg (K-H) and the Kitaev-Ising (K-I) models. We start by applying the Chebyshev toolset to compute nearest-neighbor spin correlations, specific heat and entropy of the K-H model on a 24-spin cluster. Our results are benchmarked against exact diagonalization and a popular iterative method based on thermal pure quantum states. The transitions between a variety of magnetic phases, namely ferromagnetic, Néel, zigzag and stripy antiferromagnetic and quantum spin liquid phases are obtained accurately and efficiently. We also determine the temperature dependence of the spin correlations, over more than three decades in temperature, by means of a finite temperature Chebyshev polynomial method introduced here. Finally, we report novel dynamical signatures of the quantum phase transitions in the K-I model. Our findings suggest that the efficiency, versatility and low-temperature stability of the Chebyshev framework developed here could pave the way for previously unattainable studies of quantum spin models in two dimensions.
Kitaev quantum spin liquids host Majorana fermions via the fractionalization of spins. In a magnetic field, the Majorana fermions were predicted to comprise a topological state, which has attracted great attention by the discovery of the half-quantized thermal Hall conductivity. Nevertheless, a reliable theory remains elusive for the field effect, especially at finite temperature. Here we present unbiased large-scale numerical results for the Kitaev model in a wide range of magnetic field and temperature. We find that the unconventional paramagnetic region showing fractional spin dynamics extends at finite temperature, far beyond the field range where the topological state is expected at zero temperature. Our results show the confinement-deconfinement behavior between the fractional Majorana excitations and the conventional magnons.
The long-sought quantum spin liquid is a quantum-entangled magnetic state leading to the fractionalization of spin degrees of freedom. Quasiparticles emergent from the fractionalization affect not only the ground-state properties but also thermodynamic behavior in a peculiar manner. Here we investigate how the spin dynamics evolves from the high-temperature paramagnet to the quantum spin-liquid ground state, for the Kitaev spin model describing the fractionalization into itinerant matter fermions and localized Z2 fluxes. Beyond the previous studies [J. Yoshitake, J. Nasu, and Y. Motome, Phys. Rev. Lett. 117, 157203 (2016); J. Yoshitake, J. Nasu, Y. Kato, and Y. Motome, Phys. Rev. B 96, 024438 (2017)], in which the mean-field nature of the cluster dynamical mean-field theory prevented us from studying low-temperature properties, we develop a numerical technique by applying the continuous-time quantum Monte Carlo (CTQMC) method to statistical samples generated by the quantum Monte Carlo (QMC) method in a Majorana fermion representation. This QMC+CTQMC method is fully unbiased and enables us to investigate the low-temperature spin dynamics dominated by thermally excited Z2 fluxes, including the unconventional phase transition caused by Z2 flux loops in three dimensions, which was unreachable by the previous methods. We apply this technique to the Kitaev model in both two and three dimensions. Our results clearly distinguish two cases: while the dynamics changes smoothly through the crossover in the two-dimensional honeycomb case, it exhibits singular behaviors at the phase transition in the three-dimensional hyperhoneycomb case. We show that the low-temperature spin dynamics is a sensitive probe for thermally fluctuating Z2 fluxes that behave very differently between two and three dimensions.
The S=1/2 Kitaev model provides an exactly solvable setting in which spin fractionalization produces itinerant Majorana fermions coupled to an emergent Z2 gauge field. In realistic Kitaev magnets, however, non-Kitaev interactions break integrability and make their finite-temperature dynamics difficult to calculate. We develop a finite-temperature parton mean-field theory based on the Abrikosov-fermion representation and incorporate the effects of exchange interactions through the random phase approximation to calculate the dynamical spin structure factor. For the pure Kitaev model, the method reproduces a strong, nearly dispersionless peak at a finite low energy together with a broad continuum at higher energies. It also reproduces the momentum dependence of the dynamical spin structure factor and its temperature evolution, including the downward shift of the low-energy spectral weight, in qualitative agreement with continuous-time quantum Monte Carlo calculations. Within the present theory, the contraction of the continuum is traced to the temperature-induced reduction of the bond amplitudes, which narrows the parton bands and compresses the two-parton phase space. For ferromagnetic Kitaev coupling, the energy of the low-energy peak is minimized at the Γ point for a ferromagnetic Heisenberg interaction and maximized there for an antiferromagnetic interaction. In the ferromagnetic Kitaev-Γ model, the static response develops a finite-momentum softening, indicating an instability toward the period-three magnetic order reported in numerical studies. For uniform-ferromagnetic states selected by a [111] magnetic field, we find that a broad parton continuum coexists with magnonlike dispersive structures, producing spectral features consistent with a spinon-to-magnon crossover.