Source author record

Yukitoshi Motome

Yukitoshi Motome appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.str-el cond-mat.stat-mech cond-mat.mtrl-sci quant-ph cond-mat.dis-nn cond-mat.mes-hall cond-mat.other cond-mat.quant-gas cond-mat.supr-con math-ph math.MP physics.comp-ph

Catalog footprint

What is connected

111works

12topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Asymmetric modulation of Majorana excitation spectra and nonreciprocal thermal transport in the Kitaev spin liquid under a staggered magnetic field

We present the theoretical results for the Majorana band structure and thermal transport properties in the Kitaev honeycomb model under a staggered magnetic field in addition to a uniform one. The Majorana band structure is asymmetrically modulated in momentum space, and the valley degree of freedom is activated and controlled by the magnitudes and directions of the uniform and staggered magnetic fields; as a result, the Majorana excitation is either gapped or gapless, the latter of which in general has Majorana Fermi surfaces, in contrast to the point nodes in the absence of the fields. We show that the asymmetric deformation of the Majorana band induces a nonlinear (nonreciprocal) thermal transport, in addition to the linear one. The field dependence of the linear thermal current is correlated with the magnitude of the Majorana gap, while that of the nonlinear thermal current is associated with the asymmetry of the Majorana excitation spectra. We show the estimates of the thermal currents and discuss that the linear response is experimentally measurable, whereas future improvement is required for the observation of the nonlinear one. We also discuss how the thermal currents are modulated by other additional effects of the magnetic field and contributions from conventional magnons, latter of which are expected in heterostructures to realize an internal staggered magnetic field.

preprint2022arXiv

Feasibility of Kitaev quantum spin liquids in ultracold polar molecules

Ultracold atoms and molecules trapped in optical lattices are expected to serve as simulators of strongly correlated systems and topological states of matter. A fascinating example is to realize the Kitaev quantum spin liquid by using ultracold polar molecules. However, although experimental implementation of the Kitaev-type interaction was proposed, the stability of the Kitaev quantum spin liquid has not been fully investigated thus far. Here we study a quantum spin model with long-range angle-dependent Kitaev-type interactions proposed for the polar molecules, by the pseudofermion functional renormalization group method. We reveal that the ground state is magnetically ordered in both ferromagnetic and antiferromagnetic models regardless of the spatial anisotropy of the interactions, while the isotropic case is most frustrated and closest to the realization of the Kitaev quantum spin liquid. Furthermore, by introducing a cutoff in the interaction range, we clarify how the Kitaev quantum spin liquid is destroyed by the long-range interactions. The results urge us to reconsider the feasibility of the Kitaev quantum spin liquid in ultracold polar molecules.

preprint2022arXiv

Magnetic field-temperature phase diagrams for multiple-$Q$ magnetic orderings: Exact steepest descent approach to long-range interacting spin systems

Multiple-$Q$ magnetic orderings represent magnetic textures composed of superpositions of multiple spin density waves or spin spirals, as represented by skyrmion crystals and hedgehog lattices. Such magnetic orderings have been observed in various magnetic materials in recent years, and attracted enormous attention, especially from the viewpoint of topology and emergent electromagnetic fields originating from noncoplanar magnetic structures. Although they often exhibit successive phase transitions among different multiple-$Q$ states while changing temperature and an external magnetic field, it is not straightforward to elucidate the phase diagrams, mainly due to the lack of concise theoretical tools as well as appropriate microscopic models. Here, we provide a theoretical framework for a class of effective spin models with long-range magnetic interactions mediated by conduction electrons in magnetic metals. Our framework is based on the steepest descent method with a set of self-consistent equations that leads to exact solutions in the thermodynamic limit, and has many advantages over existing methods such as biased variational calculations and numerical Monte Carlo simulations. Applying the framework to the models with instabilities toward triple- and hextuple-$Q$ magnetic orderings, we find that interesting reentrant phase transitions where the multiple-$Q$ phases appear only at finite temperature and/or nonzero magnetic field. Furthermore, we show that the multiple-$Q$ states can be topologically-nontrivial stacked skyrmion crystals or hedgehog lattices, which exhibit large net spin scalar chirality associated with nonzero skyrmion number. The results demonstrate that our framework could be a versatile tool for studying magnetic and topological phase transitions and related quantum phenomena in actual magnetic metals hosting multiple-$Q$ magnetic orderings.

preprint2022arXiv

Magnetic Hedgehog Lattice in a Centrosymmetric Cubic Metal

The hedgehog lattice (HL) is a three-dimensional topological spin texture hosting a periodic array of magnetic monopoles and antimonopoles. It has been studied theoretically for noncentrosymmetric systems with the Dzyaloshinskii-Moriya interaction, but the stability, as well as the magnetic and topological properties, remains elusive in the centrosymmetric case. We here investigate the ground state of an effective spin model with long-range bilinear and biquadratic interactions for a centrosymmetric cubic metal by simulated annealing. We show that our model stabilizes a HL composed of two pairs of left- and right-handed helices, resulting in no net scalar spin chirality, in stark contrast to the noncentrosymmetric case. We find that the HL turns into topologically-trivial conical states in an applied magnetic field. From the detailed analyses of the constituent spin helices, we clarify that the ellipticity and angles of the helical planes change gradually while increasing the magnetic field. We discuss the results in comparison with the experiments for a centrosymmetric cubic metal SrFeO$_3$.

preprint2022arXiv

Phase degree of freedom and topology in multiple-$Q$ spin textures

A periodic array of topological spin textures, such as skyrmions and hedgehogs, is called the multiple-$Q$ spin texture, as it is represented by a superposition of multiple spin density waves. Depending on the way of superposition, not only the magnetic but also the topological properties are modified, leading to a variety of quantum transport and optical phenomena caused by the emergent electromagnetic fields through the Berry phase. Among others, the phase degree of freedom of the superposed waves is potentially important for such modifications, but its effect has not been fully elucidated thus far. Here we perform systematic theoretical analyses of magnetic and topological properties of the multiple-$Q$ spin textures with the phase degree of freedom. By introducing a hyperspace with an additional dimension corresponding to the phase degree of freedom, we establish a generic framework to deal with the phase shift in the multiple-$Q$ spin textures. Applying the framework to the two-dimensional 3$Q$ spin textures, we clarify the complete topological phase diagram while changing the phase and magnetization, which depends on the types of the superposed waves. We also study the three-dimensional 4$Q$ spin textures and clarify even richer topological phase diagrams. In particular, we find novel topological phase transitions associated with the previously unidentified Dirac strings on which the hedgehogs and antihedgehogs cause pair creation and fusion. Moreover, we demonstrate that phase shifts are caused by an external magnetic field in both 3$Q$ and 4$Q$ cases by analyzing the numerical data in the previous studies. Our results illuminate the topological aspects of the skyrmion and hedgehog lattices with the phase degree of freedom, which would be extended to other multiple-$Q$ textures and useful for the exploration of topologically nontrivial magnetic phases and exotic quantum phenomena.

preprint2022arXiv

Spin dynamics in the square-lattice cupola system Ba(TiO)Cu$_4$(PO$_4$)$_4$

We report high-resolution single-crystal inelastic neutron scattering measurements on the spin-1/2 antiferromagnet Ba(TiO)Cu$_4$(PO$_4$)$_4$. This material is formed from layers of four-site \cupola" structures, oriented alternately upwards and downwards, which constitute a rather special realization of two-dimensional (2D) square-lattice magnetism. The strong Dzyaloshinskii-Moriya (DM) interaction within each cupola, or plaquette, unit has a geometry largely unexplored among the numerous studies of magnetic properties in 2D Heisenberg models with spin and spatial anisotropies. We have measured the magnetic excitations at zero field and in fields up to 5 T, finding a complex mode structure with multiple characteristic features that allow us to extract all the relevant magnetic interactions by modelling within the linear spin-wave approximation. We demonstrate that Ba(TiO)Cu$_4$(PO$_4$)$_4$ is a checkerboard system with almost equal intra- and inter-plaquette couplings, in which the intra-plaquette DM interaction is instrumental both in enforcing robust magnetic order and in opening a large gap at the Brillouin-zone center. We place our observations in the perspective of generalized phase diagrams for spin-1/2 square-lattice models and materials, where exploring anisotropies and frustration as routes to quantum disorder remains a frontier research problem.

preprint2021arXiv

Electrical conductivity in helical and conical magnetic states

We theoretically study the electrical conductivity in a one-dimensional helimagnet whose spin texture changes from helimagnetic to conical magnetic, and to forced ferromagnetic state while increasing the magnetic field along the helical axis. We find that the conductivity in the helimagnetic state at zero field depends on the electron filling and the coefficient of the spin-charge coupling. We also find that the conductivity in the conical magnetic state changes nonlinearly to the applied field, and the magnetoresistance becomes negative and positive depending on the model parameters.

preprint2021arXiv

Noncoplanar multiple-$Q$ spin textures by itinerant frustration: Effects of single-ion anisotropy and bond-dependent anisotropy

We theoretically investigate multiple-$Q$ spin textures, which are composed of superpositions of spin density waves with different wave numbers, for an effective spin model of centrosymmetric itinerant magnets. Our focus is on the interplay between biquadratic interactions arising from the spin-charge coupling and magnetic anisotropy caused by the spin-orbit coupling. Taking into account two types of the magnetic anisotropy, single-ion anisotropy and bond-dependent anisotropy, we elucidate magnetic phase diagrams for an archetypal triangular-lattice system in the absence and presence of an external magnetic field. In the case of the single-ion anisotropy, we find a plethora of multiple-$Q$ instabilities depending on the strength and the sign of the anisotropy (easy plane or easy axis), including a skyrmion crystal with topological number of two. In an external magnetic field, we find that a skyrmion crystal with topological number of one is stabilized by the in-plane (out-of-plane) magnetic field under the easy-plane (easy-axis) anisotropy. We also examine the stability of the field-induced skyrmion crystal by rotating the field direction. As a biproduct, we show that a triple-$Q$ state with nonzero chirality appears in the presence of the biquadratic interaction and the easy-axis anisotropy. Meanwhile, we find that the bond-dependent anisotropy also stabilizes both types of skyrmion crystals. We show that, however, for the skyrmion crystal with topological number of one, Bloch- and Néel-type skyrmion crystals are selectively realized depending on the sign of the bond-dependent anisotropy. Moreover, we find yet another multiple-$Q$ states, including two types of meron crystals with the skyrmion numbers of one and two. The systematic investigation will provide a reference to complex magnetic textures in centrosymmetric magnetic metals.

preprint2021arXiv

Phase shift, ellipticity, angle, and topological number in skyrmion lattices

We theoretically study skyrmion lattices realized in a Kondo lattice model on a triangular lattice, focusing on the phase, ellipticity, and angle of the constituent multiple-$Q$ waves. Analyzing the numerical data obtained in the previous study [Ozawa R, Hayami S and Motome Y 2017 Phys. Rev. Lett. 118 147205], we extract these parameters for the two types of skyrmion lattices with the skyrmion number of $1$ and $2$. We show that the topological transition between the two skyrmion lattices driven by an external magnetic field is accompanied by significant modulations of all three parameters.

preprint2020arXiv

Anomalous Hall effect in $κ$-type organic antiferromagnets

We theoretically propose a mechanism for the anomalous Hall effect (AHE) in an antiferrromagnetic (AFM) state of $κ$-type organic conductors. We incorporate the spin-orbit coupling in the effective Hubbard model on the $κ$-type lattice structure taking into account the orientation of the molecules and their arrangement with dimerization. Treating this model by means of the Hartree-Fock approximation and the linear response theory, we find that an intrinsic contribution to the Hall conductivity becomes nonzero in the electron-doped AFM metallic phase with a small canted ferromagnetic moment. We show that, contrary to the conventional wisdom, the spin canting is irrelevant to the Hall response; the nonzero Hall conductivity originates from the collinear component of the AFM order in the presence of the spin-orbit coupling. These features are well explained analytically in the limit of strong dimerization on the anisotropic triangular lattice. Furthermore, we present an intuitive picture for the present AHE by considering the real-space configuration of emergent magnetic fluxes. We also find that the Hall response appears even in the undoped AFM insulating phase at nonzero frequency as the magneto-optical Kerr effect, which is enhanced around the charge transfer excitations. We discuss possible detections of the AHE in ET based compounds.

preprint2020arXiv

Antiferromagnet-semiconductor van der Waals heterostructures: interlayer interplay of exciton with magnetic ordering

Van der Waals (vdW) heterostructures have attracted great interest because of their rich material combinations.The discovery of two-dimensional magnets has provided a new platform for magnetic vdW heterointerfaces; however, research on magnetic vdW heterointerfaces has been limited to those with ferromagnetic surfaces. Here we report a magnetic vdW heterointerface using layered intralayer-antiferromagnetic MPSe3 (M=Mn, Fe) and monolayer transition metal dichalcogenides (TMDs). We found an anomalous upshift of the excitonic peak in monolayer TMDs below the antiferromagnetic transition temperature in the MPSe3, capturing a signature of the interlayer exciton-magnon coupling. This is a concept extended from single materials to heterointerfaces. Moreover, this coupling strongly depends on the in-plane magnetic structure and stacking direction, showing its sensitivity to their magnetic interfaces. Our finding offers an opportunity to investigate interactions between elementary excitations in different materials across interfaces and to search for new functions of magnetic vdW heterointerfaces.

preprint2020arXiv

Magnetic hedgehog lattices in noncentrosymmetric metals

The magnetic hedgehog lattice (HL) is a noncoplanar magnetic texture with a periodic array of magnetic monopoles and anti-monopoles. Despite phenomenological and numerical studies thus far, there remain open issues on the microscopic origin, especially with respect to the recent experimental findings of two different types of HLs even at zero magnetic field. Here, we study the stability of the HLs for an effective spin model with long-range interactions arising from itinerant nature of electrons. By variational calculations and simulated annealing, we find that the HLs are stabilized in the ground state at zero magnetic field by the synergetic effect of the anti-symmetric exchange interactions generated by the spin-orbit coupling and the multiple-spin interactions generated by the spin-charge coupling. We also clarify the phase diagram in the magnetic fields, which includes topological phase transitions with pair annihilation of the monopoles and anti-monopoles depending on the field directions.

preprint2020arXiv

Majorana-Magnon Crossover by a Magnetic Field in the Kitaev Model: Continuous-Time Quantum Monte Carlo Study

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.

preprint2020arXiv

Optical Hall response in spin-orbit coupled metals: Comparative study of magnetic cluster monopole, quadrupole, and toroidal orders

The optical Hall response is theoretically studied for spin-orbit coupled metals with ferroic orders of cluster-type magnetic multipoles. We find that different magnetic multipoles give rise to distinct spectra in the optical Hall conductivity. In the cases of monopole and quadrupole orders, the optical Hall response appears predominantly in high- and low-energy regions, which correspond to the energy scales of electron correlation and kinetic energy, respectively, while the response is dispersed and rather weak in the case of toroidal order. By decomposing the spectra into different interband contributions, we reveal selection rules stemming from the interplay between the antisymmetric spin-orbit coupling and the underlying multipoles. Our results suggest that the optical Hall measurement is useful to detect and distinguish the cluster-type magnetic multipole orders.

preprint2020arXiv

Perovskite as a spin current generator

We theoretically show that materials with perovskite-type crystal structures provide a platform for spin current generation, taking advantage of a mechanism requiring neither the spin-orbit coupling nor a ferromagnetic moment, but is based on spin-split band structures in certain kinds of collinear antiferromagnetic states. By investigating a multiband Hubbard model for transition metal compounds by means of the Hartree-Fock approximation and the Boltzmann transport theory, we find that a pure spin current is induced by an electric field applied to a C-type antiferromagnetic metallic phase. The spin current generation originates from a cooperative effect of spatially anisotropic electron transfer integrals owing to the GdFeO$_3$-type lattice distortion, which is ubiquitous in many perovskites, and the collinear spin configuration. We discuss our finding from the symmetry point of view, in comparison with other spin current generator candidates with collinear antiferomagnetism. We also propose several ways to detect the phenomenon in candidate perovskite materials.

preprint2020arXiv

Phase transitions between helices, vortices, and hedgehogs driven by spatial anisotropy in chiral magnets

Superpositions of spin helices can yield topological spin textures, such as two-dimensional vortices and skyrmions, and three-dimensional hedgehogs. Their topological nature and spatial dimensionality depend on the number and relative directions of the constituent helices. This allows mutual transformation between the topological spin textures by controlling the spatial anisotropy. Here we theoretically study the effect of anisotropy in the magnetic interactions for an effective spin model for chiral magnetic metals. By variational calculations for both cases with triple and quadruple superpositions, we find that the hedgehog lattices, which are stable in the isotropic case, are deformed by the anisotropy, and eventually changed into other spin textures with reduced dimension, such as helices and vortices. We also clarify the changes of topological properties by tracing the real-space positions of magnetic monopoles and antimonopoles as well as the emergent magnetic field generated by the noncoplanar spin textures. Our results suggest possible control of the topological spin textures, e.g., by uniaxial pressure and chemical substitution in chiral materials.

preprint2020arXiv

Square skyrmion crystal in centrosymmetric itinerant magnets

We theoretically investigate the origin of the square-type skyrmion crystal in centrosymmetric itinerant magnets, motivated from the recent experimental finding in GdRu$_2$Si$_2$ [N. D. Khanh $\mathit{et}$ $\mathit{al.}$, Nat. Nanotech. $\mathbf{15}$, 444 (2020)]. By simulated annealing for an effective spin model derived from the Kondo lattice model on a square lattice, we find that a square skyrmion crystal composed of a superposition of two spin helices is stabilized in a magnetic field by synergy between the positive biquadratic, bond-dependent anisotropic, and easy-axis anisotropic interactions. This is in stark contrast to triangular skyrmion crystals which are stabilized by only one of the three, suggesting that the square skyrmion crystal is characteristic of itinerant magnets with magnetic anisotropy. We also show that a variety of noncollinear and noncoplanar spin textures appear depending on the model parameters as well as the applied magnetic field. The present systematic study will be useful not only for identifying the key ingredients in GdRu$_2$Si$_2$ but also for exploring further skyrmion-hosting materials in centrosymmetric itinerant magnets.

preprint2020arXiv

The range of non-Kitaev terms and fractional particles in RuCl$_3$

Significant efforts have focused on the magnetic excitations of relativistic Mott insulators, predicted to realize the Kitaev quantum spin liquid (QSL). This exactly solvable model involves a highly entangled state resulting from bond-dependent Ising interactions that produce excitations which are non-local in terms of spin flips. A key challenge in real materials is identifying the relative size of the non-Kitaev terms and their role in the emergence or suppression of fractional excitations. Here, we identify the energy and temperature boundaries of non-Kitaev interactions by direct comparison of the Raman susceptibility of RuCl3 with quantum Monte Carlo (QMC) results for the Kitaev QSLs. Moreover, we further confirm the fractional nature of the magnetic excitations, which is given by creating a pair of fermionic quasiparticles. Interestingly, this fermionic response remains valid in the non-Kitaev range. Our results and focus on the use of the Raman susceptibility provide a stringent new test for future theoretical and experimental studies of QSLs.

preprint2020arXiv

Thermodynamic and Transport Properties in Disordered Kitaev Models

Effects of bond randomness and site dilution are systematically investigated for the Kitaev model describing a quantum spin liquid with fractional excitations of itinerant Majorana fermions and localized fluxes. We find that, in the high-temperature region where the itinerant Majorana fermions release their entropy, both types of disorders suppress the longitudinal thermal conductivity while keeping the specific heat almost unchanged. This suggests that both disorders reduce the mean-free path of the Majorana fermions. On the other hand, in the low-temperature region, the other specific heat peak associated with the entropy release from the localized fluxes is suppressed for both cases, but it is broadened and shifted to the lower-temperature side by the bond randomness, while the position and the width are almost unchanged against the site dilution. Contrasting behavior is also found in the thermal Hall effect under a magnetic field; the half quantization of the thermal Hall conductivity is fragile against the site dilution, while it remains for the bond randomness despite the reduced onset temperature. We discuss the contrasting behavior from the stability of the topological nature by calculating flux condensation and Majorana excitation gap.

preprint2020arXiv

Thermodynamic classification of three-dimensional Kitaev spin liquids

In the field of frustrated magnetism, Kitaev models provide a unique framework to study the phenomena of spin fractionalization and emergent lattice gauge theories in two and three spatial dimensions. Their ground states are quantum spin liquids, which can typically be described in terms of a Majorana band structure and an ordering of the underlying $\mathbb{Z}_2$ gauge structure. Here we provide a comprehensive classification of the "gauge physics" of a family of elementary three-dimensional Kitaev models, discussing how their thermodynamics and ground state order depends on the underlying lattice geometry. Using large-scale, sign-free quantum Monte Carlo simulations we show that the ground-state gauge order can generally be understood in terms of the length of elementary plaquettes -- a result which extends the applicability of Lieb's theorem to lattice geometries beyond its original scope. At finite temperatures, the proliferation of (gapped) vison excitations destroys the gauge order at a critical temperature scale, which we show to correlate with the size of vison gap for the family of three-dimensional Kitaev models. We also discuss two notable exceptions where the lattice structure gives rise to "gauge frustration" or intertwines the gauge ordering with time-reversal symmetry breaking. In a more general context, the thermodynamic gauge transitions in such 3D Kitaev models are one of the most natural settings for phase transitions beyond the standard Landau-Ginzburg-Wilson paradigm.

preprint2020arXiv

Tracing Monopoles and Anti-monopoles in a Magnetic Hedgehog Lattice

The magnetic hedgehog lattice (HL), which was recently discovered in the $B$20-type chiral magnet MnSi$_{1-x}$Ge$_x$, is a topological spin texture with a periodic array of magnetic monopoles and anti-monopoles. Within the continuum approximation, the monopoles and anti-monopoles are predicted to move, collide, and pair annihilate in an applied magnetic field, but it remains unclear how the lattice discretization affects their motions. Here, we study the trajectories of monopoles and anti-monopoles in a lattice system by simulated annealing with field sweep. We show that the monopoles and anti-monopoles move and repel before pair annihilations. We also clarify that their motions are closely related with the field dependence of the scalar spin chirality.

preprint2019arXiv

Antiferromagnetic Kitaev Interactions in Polar Spin-Orbit Mott Insulators

A bond-directional anisotropic exchange interaction, called the Kitaev interaction, is a promising route to realize quantum spin liquids. The Kitaev interactions were found in Mott insulators with the strong spin-orbit coupling, in the presence of quantum interference between indirect electron transfers. Here we theoretically propose a different scenario by introducing a polar structural asymmetry that unbalances the quantum interference. We show that the imbalance activates additional exchange processes and gives rise to a dominant antiferromagnetic Kitaev interaction, in stark contrast to the conventional ferromagnetic ones. We demonstrate by ab initio calculations that polar Ru trihalides with multiple anions, $α$-RuH$_{3/2}X_{3/2}$ ($X$=Cl and Br), exhibit the antiferromagnetic Kitaev interaction whose magnitude is several times larger compared to existing candidates. Our proposal opens the way for materializing the Kitaev spin liquids in unexplored parameter regions.

preprint2019arXiv

Chiral spin liquids with crystalline $\mathbb{Z}_2$ gauge order in a three-dimensional Kitaev model

Chiral spin liquids (CSLs) are time-reversal symmetry breaking ground states of frustrated quantum magnets that show no long-range magnetic ordering, but instead exhibit topological order and fractional excitations. Their realization in simple and tractable microscopic models has, however, remained an open challenge for almost two decades until it was realized that Kitaev models on lattices with odd-length loops are natural hosts for such states, even in the absence of a time-reversal symmetry breaking magnetic field. Here we report on the formation of CSLs in a three-dimensional Kitaev model, which differ from their widely studied two-dimensional counterparts, namely, they exhibit a crystalline ordering of the $\mathbb{Z}_2$ gauge fluxes and thereby break some of the underlying lattice symmetries. We study the formation of these unconventional CSLs via extensive quantum Monte Carlo simulations and demonstrate that they are separated from the featureless paramagnet at high temperatures by a single first-order transition at which both time-reversal and lattice symmetries are simultaneously broken. Using variational approaches for the ground state, we explore the effect of varying the Kitaev couplings and find at least five distinct CSL phases, all of which possess crystalline ordering of the $\mathbb{Z}_2$ gauge fluxes. For some of these phases, the complementary itinerant Majorana fermions exhibit gapless band structures with topological features such as Weyl nodes or nodal lines in the bulk and Fermi arc or drumhead surface states.

preprint2016arXiv

Asymmetric Magnon Excitation by Spontaneous Toroidal Ordering

Effects of spontaneous toroidal ordering on magnetic excitation are theoretically investigated for a localized spin model that includes a staggered Dzyaloshinsky-Moriya interaction and anisotropic exchange interactions, which arise from the antisymmetric spin-orbit coupling and the multi-orbital correlation effect. We show that the model exhibits a Neel-type antiferromagnetic order, which simultaneously accompanies a ferroic toroidal order. We find that the occurrence of toroidal order modulates the magnon dispersion in an asymmetric way with respect to the wave number: a toroidal dipole order on the zigzag chain leads to a band-bottom shift, while a toroidal octupole order on the honeycomb lattice gives rise to a valley splitting. These asymmetric magnon excitations could be a source of unusual magnetic responses, such as nonreciprocal magnon transport. A variety of modulations are discussed while changing the lattice and magnetic symmetries. Implications of candidate materials for asymmetric magnon excitations are presented.

preprint2016arXiv

Emergent spin-valley-orbital physics by spontaneous parity breaking

The spin-orbit coupling in the absence of spatial inversion symmetry plays an important role in realizing intriguing electronic states in solids, such as topological insulators and unconventional superconductivity. Usually, the inversion symmetry breaking is inherent in the lattice structures, and hence, it is not easy to control these interesting properties by external parameters. We here theoretically investigate the possibility of generating the spin-orbital entanglement by spontaneous electronic ordering caused by electron correlations. In particular, we focus on the centrosymmetric lattices with local asymmetry at the lattice sites, e.g., zig-zag, honeycomb, and diamond structures. In such systems, conventional staggered orders, such as charge order and antiferromagnetic order, break the inversion symmetry and activate the antisymmetric spin-orbit coupling, which is hidden in a sublattice-dependent form in the paramagnetic state. Considering a minimal two-orbital model on a honeycomb lattice, we scrutinize the explicit form of the antisymmetric spin-orbit coupling for all the possible staggered charge, spin, orbital, and spin-orbital orders. We show that the complete table is useful for understanding of spin-valley-orbital physics, such as spin and valley splitting in the electronic band structure and generalized magnetoelectric responses in not only spin but also orbital and spin-orbital channels, reflecting in peculiar magnetic, elastic, and optical properties in solids.

preprint2016arXiv

Engineering chiral density waves and topological band structures by multiple-$Q$ superpositions of collinear up-up-down-down orders

Magnetic orders characterized by multiple ordering vectors harbor noncollinear and noncoplanar spin textures and can be a source of unusual electronic properties through the spin Berry phase mechanism. We theoretically show that such multiple-$Q$ states are stabilized in itinerant magnets in the form of superpositions of collinear up-up-down-down (UUDD) spin states, which accompany the density waves of vector and scalar chirality. The result is drawn by examining the ground state of the Kondo lattice model with classical localized moments, especially when the Fermi surface is tuned to be partially nested by the symmetry-related commensurate vectors. We unveil the instability toward the multiple-$Q$ UUDD states with chirality density waves, using the perturbative theory, variational calculations, and large-scale Langevin dynamics simulations. We also show that the chirality density waves can induce rich nontrivial topology of electronic structures, such as the massless Dirac semimetal, Chern insulator with quantized topological Hall response, and peculiar edge states which depend on the phase of chirality density waves at the edges.

preprint2016arXiv

Fractional Spin Fluctuation as a Precursor of Quantum Spin Liquids: Majorana Dynamical Mean-Field Study for the Kitaev Model

Experimental identification of quantum spin liquids remains a challenge, as the pristine nature is to be seen in asymptotically low temperatures. We here theoretically show that the precursor of quantum spin liquids appears in the spin dynamics in the paramagnetic state over a wide temperature range. Using the cluster dynamical mean-field theory and the continuous-time quantum Monte Carlo method, which are newly developed in the Majorana fermion representation, we calculate the dynamical spin structure factor, relaxation rate in nuclear magnetic resonance, and magnetic susceptibility for the honeycomb Kitaev model whose ground state is a canonical example of the quantum spin liquid. We find that dynamical spin correlations show peculiar temperature and frequency dependence even below the temperature where static correlations saturate. The results provide the experimentally-accessible symptoms of the fluctuating fractionalized spins evincing the quantum spin liquids.

preprint2016arXiv

Topological Insulators from Electronic Superstructures

The possibility of realizing topological insulators by spontaneous formation of electronic superstructure is theoretically investigated in a minimal two-orbital model including both the spin-orbit coupling and electron correlations on a triangular lattice. Using the mean-field approximation, we show that the model exhibits several different types of charge ordered insulators, where the charge disproportionation forms a honeycomb or kagome superstructure. We find that the charge ordered insulators in the presence of strong spin-orbit coupling can be topological insulators showing quantized spin Hall conductivity. Their band gap is dependent on electron correlations as well as the spin-orbit coupling, and even vanishes with showing the massless Dirac dispersion at the transition to a trivial charge ordered insulator. Our results suggest a new route to realize and control topological states of quantum matter by the interplay between the spin-orbit coupling and electron correlations.

preprint2015arXiv

A strong coupling theory for electron-mediated interactions in double-exchange models

We present a theoretical framework for evaluating effective interactions between localized spins mediated by itinerant electrons in double-exchange models. Performing the expansion with respect to the spin-dependent part of the electron hopping terms, we show a systematic way of constructing the effective spin model in the large Hund's coupling limit. As a benchmark, we examine the accuracy of this method by comparing the results with the numerical solutions for the spin-ice type model on a pyrochlore lattice. We also discuss an extension of the method to the double-exchange models with Heisenberg and $XY$ localized spins.

preprint2015arXiv

Analysis of Charge-spin-orbital Fluctuations by Ab Initio Calculation and Random Phase Approximation: Application to Non-coplanar Antiferromagnet Cd2Os2O7

We present a systematic analysis on the basis of ab initio calculations and many-body perturbation theory for clarifying the dominant fluctuation in complex charge-spin-orbital coupled systems. For a tight-binding multiband model obtained from the maximally-localized Wannier function analysis of the band structure by the local density approximation, we take into account electron correlations at the level of random phase approximation. To identify the dominant fluctuation, we carry out the eigenmode analysis of the generalized susceptibility that includes all the multiple degrees of freedom: charge, spin, and orbital. We apply this method to the paramagnetic metallic phase of a pyrochlore oxide Cd$_2$Os$_2$O$_7$, which shows a metal-insulator transition accompanied by a peculiar noncoplanar antiferromagnetic order of all-in all-out type. We find that the corresponding spin fluctuation is dominantly enhanced by the on-site Coulomb repulsions in the presence of strong spin-orbit coupling and trigonal crystal field splitting. Our results indicate that the combined method offers an effective tool for the systematic analysis of potential instabilities in strongly correlated electron materials.

preprint2015arXiv

Charge-spin-orbital Fluctuations in Mixed Valence Spinels: Comparative Study of AlV2O4 and LiV2O4

Mixed valence spinels provide a fertile playground for the interplay between charge, spin, and orbital degrees of freedom in strongly correlated electrons on a geometrically frustrated lattice. Among them, AlV$_2$O$_4$ and LiV$_2$O$_4$ exhibit contrasting and puzzling behavior: self-organization of seven-site clusters and heavy fermion behavior. We theoretically perform a comparative study of charge-spin-orbital fluctuations in these two compounds, on the basis of the multiband Hubbard models constructed by using the maximally-localized Wannier functions obtained from the ab initio band calculations. Performing the eigenmode analysis of the generalized susceptibility, we find that, in AlV$_2$O$_4$, the relevant fluctuation appears in the charge sector in $σ$-bonding type orbitals. In contrast, in LiV$_2$O$_4$, optical-type spin fluctuations in the $a_{\rm 1g}$ orbital are enhanced at an incommensurate wave number at low temperature. Implications from the comparative study are discussed for the contrasting behavior, including the metal-insulator transition under pressure in LiV$_2$O$_4$.

preprint2015arXiv

Magnetic three states of matter: A quantum Monte Carlo study of spin liquids

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 $Z_2$ 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.

preprint2015arXiv

Phase Diagram of the Kitaev-type Model on a Decorated Honeycomb Lattice in the Isolated Dimer Limit

An effective model in the isolated dimer limit of the Kitaev-type model on a decorated honeycomb lattice is investigated at finite temperature. The ground state of this model is exactly shown to be a chiral spin liquid with spontaneous breaking of time reversal symmetry. We elaborate the finite-temperature phase diagram by using the mean-field approximation and Monte Carlo simulation. We find that the phase transition between the high-temperature paramagnetic phase and the low-temperature chiral spin liquid phase is always of second order in the Monte Carlo results, although a tricritical point appears in the mean-field phase diagram. The finite-size scaling analysis of the Monte Carlo data indicates that the phase transition belongs to the two-dimensional Ising universality class.

preprint2015arXiv

Spin correlation and Majorana spectrum in chiral spin liquids in a decorated-honeycomb Kitaev model

Temperature evolution of the spin correlation and excitation spectrum is investigated for a Kitaev model defined on a decorated honeycomb lattice by using the quantum Monte Carlo simulation in the Majorana fermion representation. The ground state of this quantum spin model is given by two kinds of chiral spin liquids: one is topologically trivial with Abelian anyon excitations, and the other is topologically nontrivial with non-Abelian anyon excitations. While lowering temperature, the model exhibits several crossovers in the paramagnetic state, which originate from the fractionalization of quantum spins into Majorana fermions, in addition to a phase transition associated with time reversal symmetry breaking. We show that the spin correlation develops around the crossover temperatures, whereas it shows a slight change at the critical temperature, as in other Kitaev-type models. We also calculate the excitation spectrum in terms of Majorana fermions, and find that the excitation gap in the non-Abelian phase is fragile against thermal fluctuations of the $Z_2$ fluxes, while that in the Abelian phase remains open.

preprint2015arXiv

Spontaneous formation of kagome network and Dirac half-semimetal on a triangular lattice

In spin-charge coupled systems, geometrical frustration of underlying lattice structures can give rise to nontrivial magnetic orders and electronic states. Here we explore such a possibility in the Kondo lattice model with classical localized spins on a triangular lattice by using a variational calculation and simulated annealing. We find that the system exhibits a four-sublattice collinear ferrimagnetic phase at 5/8 filling for a large Hund's-rule coupling. In this state, the system spontaneously differentiates into the up-spin kagome network and the isolated down-spin sites, which we call the kagome network formation. In the kagome network state, the system becomes Dirac half-semimetallic: The electronic structure shows a massless Dirac node at the Fermi level, and the Dirac electrons are almost fully spin polarized due to the large Hund's-rule coupling. We also study the effect of off-site Coulomb repulsion in the kagome network phase where the system is effectively regarded as a 1/3-filling spinless fermion system on the kagome lattice. We find that, at the level of the mean-field approximation, a $\sqrt{3} \times \sqrt{3}$-type charge order occurs in the kagome network state, implying the possibility of fractional charge excitations in this triangular lattice system. Moreover, we demonstrate that the kagome network formation with fully-polarized Dirac electrons are controllable by an external magnetic field.

preprint2015arXiv

Spontaneous multipole ordering by local parity mixing

Broken spatial inversion symmetry in spin-orbital coupled systems leads to a mixing between orbitals with different parity, which results in unusual electronic structures and transport properties. We theoretically investigate the possibility of multipole ordering induced by a parity mixing. In particular, we focus on the system in which the parity mixing appears in a sublattice-dependent form. Starting from the periodic Anderson model with such a local parity mixing, we derive an extended Kondo lattice model with sublattice-dependent antisymmetric exchange couplings between itinerant electrons and localized spins. By the variational calculation, simulated annealing, and Monte Carlo simulation, we show that the model on a quasi-one-dimensional zig-zag lattice exhibits an odd-parity multipole order composed of magnetic toroidal and quadrupole components at and near half filling. The multipole order causes a band deformation with the band bottom shift and a magnetoelectric response. The results suggest that unusual odd-parity multipole orders will be widely observed in multi-orbital systems with local parity mixing.

preprint2015arXiv

Theory of Valence Transition in BiNiO$_3$

Motivated by the colossal negative thermal expansion recently found in BiNiO$_3$, the valence transition accompanied by the charge transfer between the Bi and Ni sites is theoretically studied. We introduce an effective model for Bi-$6s$ and Ni-$3d$ orbitals with taking into account the valence skipping of Bi cations, and investigate the ground-state and finite-temperature phase diagrams within the mean-field approximation. We find that the valence transition is caused by commensurate locking of the electron filling in each orbital associated with charge and magnetic orderings, and the critical temperature and the nature of the transitions are strongly affected by the relative energy between the Bi and Ni levels and the effective electron-electron interaction in the Bi sites. The obtained phase diagram well explains the temperature- and pressure-driven valence transitions in BiNiO$_3$ and the systematic variation of valence states for a series of Bi and Pb perovskite oxides.

preprint2015arXiv

Thermal fractionalization of quantum spins in a Kitaev model: $T$-linear specific heat and coherent transport of Majorana fermions

Finite-temperature ($T$) properties of a Kitaev model defined on a honeycomb lattice are investigated by a quantum Monte Carlo simulation, from the viewpoint of fractionalization of quantum $S=1/2$ spins into two types of Majorana fermions, itinerant and localized. In this system, the entropy is released successively at two well-separated $T$ scales, as a clear indication of the thermal fractionalization. We show that the high-$T$ crossover, which is driven by itinerant Majorana fermions, is closely related with the development of nearest-neighbor spin correlations. On the other hand, the low-$T$ crossover originates in thermal fluctuations of fluxes composed of localized Majorana fermions, by which the spectrum of itinerant Majorana fermions is significantly disturbed. As a consequence, in the intermediate-$T$ range between the two crossovers, the system exhibits $T$-linear behavior in the specific heat and coherent transport of Majorana fermions, which are unexpected for the Dirac semimetallic spectrum in the low-$T$ limit. We also show that the flux fluctuations tend to open an energy gap in the Majorana spectrum near the gapless-gapped phase boundary. Our results indicate that the fractionalization is experimentally observable in the specific heat, spin correlations, and transport properties.

preprint2015arXiv

Thermodynamics of Chiral Spin Liquids with Abelian and Non-Abelian Anyons

Thermodynamic properties of chiral spin liquids are investigated for a variant of the Kitaev model defined on a decorated honeycomb lattice. Using the quantum Monte Carlo simulation, we find that the model exhibits a finite-temperature phase transition associated with the time reversal symmetry breaking, in both topologically trivial and nontrivial regions. While changing the exchange constants, the phase transition changes from continuous to discontinuous one, apparently correlated with the change in the excitations from Abelian to non-Abelian anyons. We show this coincidence by computing the topological quantities: the Chern number and the thermal Hall conductivity. In addition, we find, as a diagnostic of the chiral spin liquids, successive crossovers with multi-stage entropy release above the critical temperature, which indicates that the hierarchical fractionalization of a quantum spin occurs differently between the two regions.

preprint2015arXiv

Vortex Crystals with Chiral Stripes in Itinerant Magnets

We study noncoplanar magnetic ordering in frustrated itinerant magnets. For a family of Kondo square lattice models with classical local moments, we find that a double-$Q$ noncoplanar vortex crystal has lower energy than the single-$Q$ helical order expected from the Ruderman-Kittel-Kasuya-Yosida interaction whenever the lattice symmetry dictates four global maxima in the bare magnetic susceptibility. By expanding in the small Kondo exchange and the degree of noncoplanarity, we demonstrate that this noncoplanar state arises from a Fermi surface instability, and it is generic for a wide range of electron filling fractions whenever the two ordering wave vectors connect independent sections of the Fermi surface.

preprint2014arXiv

3D Dirac Electrons on a Cubic Lattice with Noncoplanar Multiple-Q Order

Noncollinear and noncoplanar spin textures in solids manifest themselves not only in their peculiar magnetism but also in unusual electronic and transport properties. We here report our theoretical studies of a noncoplanar order on a simple cubic lattice and its influence on the electronic structure. We show that a four-sublattice triple-Q order induces three-dimensional massless Dirac electrons at commensurate electron fillings. The Dirac state is doubly degenerate, while it splits into a pair of Weyl nodes by lifting the degeneracy by an external magnetic field; the system is turned into a Weyl semimetal in applied field. In addition, we point out the triple-Q Hamiltonian in the strong coupling limit is equivalent to the 3D π-flux model relevant to an AIII topological insulator. We examine the stability of such a triple-Q order in two fundamental models for correlated electron systems: a Kondo lattice model with classical localized spins and a periodic Anderson model. For the Kondo lattice model, performing a variational calculation and Monte Carlo simulation, we show that the triple-Q order is widely stabilized around 1/4 filling. For the periodic Anderson model, we also show the stability of the same triple-Q state by using the mean-field approximation. For both models, the triple-Q order is widely stabilized via the couplings between conduction electrons and localized electrons even without any explicit competing magnetic interactions and geometrical frustration. We also show that the Dirac electrons induce peculiar surface states: Fermi "arcs" connecting the projected Dirac points, similarly to Weyl semimetals.

preprint2014arXiv

Effect of magnetoelastic coupling on spin-glass behavior in Heisenberg pyrochlore antiferromagnets with bond disorder

Motivated by puzzling aspects of spin-glass behavior reported in frustrated magnetic materials, we theoretically investigate effects of magnetoelastic coupling in geometrically frustrated classical spin models. In particular, we consider bond-disordered Heisenberg antiferromagnets on a pyrochlore lattice coupled to local lattice distortions. By integrating out the lattice degree of freedom, we derive an effective spin-only model, the bilinear-biquadratic model with bond disorder, which is analyzed by classical Monte Carlo simulations. First, we discuss the phase diagrams as well as thermodynamic and magnetic properties. We show that the spin-glass transition temperature is largely enhanced by the spin-lattice coupling $b$ in the weakly disordered regime. This enhancement is ascribed to the suppression of thermal fluctuations in semidiscrete degenerate manifold formed in the presence of the spin-lattice coupling. We also find that, as increasing the strength of disorder, the system shows a concomitant transition of the nematic order and spin glass at a temperature determined by $b$, being almost independent of bond disorder. Although further-neighbor exchange interactions originating in the cooperative lattice distortions result in the spin-lattice order in the weakly disordered regime, the concomitant transition remains robust. We investigate the nature of the concomitant transition by analyzing the hysteresis of the magnetic susceptibility, the nonlinear susceptibility, and the specific heat. Furthermore, we discuss linear and nonlinear magnetic susceptibilities in the high-temperature paramagnetic phase as well as single-spin flip dynamics in the nematic phase. All these results are discussed in comparison with experiments for typical pyrochlore magnets, such as Y2Mo2O7 and ZnCr2O4.

preprint2014arXiv

Entanglement Spectrum in Cluster Dynamical Mean-Field Theory

We study the entanglement spectrum of the Hubbard model at half filling on a kagome lattice. The entanglement spectrum is defined by the set of eigenvalues of reduced thermal density matrix, which is naturally obtained in the framework of the dynamical mean-field theory. Adopting the cluster dynamical mean-field theory combined with continuous-time auxiliary-field Monte Carlo method, we calculate the entanglement spectrum for a three-site triangular cluster in the kagome Hubbard model. We find that the results at the three-particle sector well captures the qualitative nature of the system. In particular, the eigenvalue of the reduced density matrix, corresponding to the chiral degrees of freedom, exhibits characteristic temperature scale T_{\rm chiral}, below which a metallic state with large quasiparticle mass is stabilized. The entanglement spectra at different particle number sectors also exhibit characteristic changes around T_{\rm chiral}, implying the development of inter-triangular ferromagnetic correlations in the correlated metallic regime.

preprint2014arXiv

Low-energy Majorana states in spin-liquid transitions in a three-dimensional Kitaev model

A three-dimensional Kitaev model on a hyperhoneycomb lattice is investigated numerically at finite temperature. The Kitaev model is one of the solvable quantum spin models, where the ground state is given by gapped and gapless spin liquids, depending on the anisotropy of the interactions. This model can be rewritten as a free Majorana fermion system coupled with $Z_2$ variables. The density of states of Majorana fermions shows an excitation gap in the gapped region, while it is semimetallic in the gapless region reflecting the Dirac node. Performing the Monte Carlo simulation, we calculate the temperature dependence of the Majorana spectra. We find that the semimetallic dip is filled up as temperature increases in the gapless region, but surprisingly, the spectrum develops an excitation gap in the region near the gapless-gapped boundary. Such changes of the low-energy spectrum appear sharply at the transition temperature from the spin liquid to the paramagnetic state. The results indicate that thermal fluctuations of the $Z_2$ fields significantly influence the low-energy state of Majorana fermions, especially in the spin liquid formation.

preprint2014arXiv

Magnetic Field Effect in One-Dimensional Charge Ordering Systems

We study effects of an external magnetic field on charge ordering in the one-dimensional extended Hubbard model at quarter filling by means of the quantum Monte Carlo method. We find that the Zeeman coupling enhances the charge order correlation, which is prominent when the system is located in the critical regime near the charge ordering transition at zero magnetic field. This behavior is interpreted in terms of the crossover to the fully spin-polarized limit where the model is exactly solvable. Furthermore, by incorporating the interchain Coulomb repulsion, we show that the charge-ordering transition temperature is raised by the magnetic field. We also discuss the relevance of our results to magnetoresistance effects observed in molecular conductors.

preprint2014arXiv

Multiple-Q Instability by (d-2)-dimensional Connections of Fermi Surfaces

We propose a general mechanism of multiple-Q ordering with noncollinear and noncoplanar spin textures in itinerant magnets. By analyzing the fourth-order perturbation with respect to the spin-charge coupling, we find that the instability toward multiple-Q ordering is caused by (d-2)-dimensional connections of the Fermi surfaces in the extended Brillouin zone by the multiple-Q wave vectors: zero(one)-dimensional point (line) connections in the two(three)-dimensional systems. The instabilities are obtained as the "fixed points" in the two-parameter flow diagram, irrespective of lattice structures. The hidden Fermi surface instability provides a universal origin of noncollinear and noncoplanar instabilities common to frustrated and unfrustrated lattice systems.

preprint2014arXiv

Quantum spin Hall effect in a two-orbital model on a honeycomb lattice

The spin Hall effect is investigated in a two-orbital tight-binding model on a honeycomb lattice. We show that the model exhibits three topologically-different insulating phases at half filling, which are distinguished by different quantized values of the spin Hall conductivity. We analytically determine the phase boundaries, where the valence and conduction bands touch with each other with forming the Dirac nodes at the Fermi level. The results are discussed in terms of the effective antisymmetric spin-orbit coupling. The relation to the Kane-Mele model and implications for a magnetoelectric effect are also discussed.

preprint2014arXiv

Reconstruction of chiral edge states in a magnetic Chern insulator

We study reconstruction of the chiral edge states in a magnetic Chern insulator in the Kondo lattice model on a triangular lattice at 1/4 filling. In this state, the spin scalar chirality associated with a four-sublattice noncoplanar magnetic order assures the topological nature characterized by a nonzero Chern number. Performing a Langevin-based numerical simulation for finite-size clusters with the open boundary condition in one direction, we clarify how the magnetic and electronic properties are modulated near the edges of the system. As a result, we find that the magnetic state near the edges is reconstructed to develop ferromagnetic spin correlations. At the same time, the electronic state is also modified and the chiral edge current is enhanced. We discuss this enhancement from the viewpoint of the electronic band structure for the gapless edge states.

preprint2014arXiv

Spin-orbit coupling in octamers in a spinel sulfide CuIr$_2$S$_4$: Competition between spin-singlet and quadrupolar states, and its relevance to remnant paramagnetism

We theoretically investigate magnetic properties in the low-temperature phase with the formation of eight-site clusters, octamers, in the spinel compound CuIr$_2$S$_4$. The octamer state was considered to be a spin-singlet state induced by a Peieirls instability through the strong anisotropy of $d$ orbitals, the so-called orbital Peierls state. We reexamine this picture by taking into account the spin-orbit coupling which was ignored in the previous study. We derive a low-energy effective model between $j_{\rm eff}=1/2$ quasispins on Ir$^{4+}$ cations in an octamer from the multiorbital Hubbard model with the strong spin-orbit coupling by performing the perturbation expansion from the strong correlation limit. The effective Hamiltonian is in the form of the Kitaev-Heisenberg model but with an additional interaction, a symmetric off-diagonal exchange interaction originating from the perturbation process including both d-d and d-p-d hopping. Analyzing the effective Hamiltonian on two sites and the octamer by the exact diagonalization, we find that there is a competition between a spin-singlet state and a quadrupolar state. The former singlet state is a conventional one, adiabatically connected to the orbital-Peierls state. On the other hand, the latter quadrupolar state is stabilized by the additional interaction, which consists of a linear combination of different total spin momenta along the spin quantization axis. In the competing region, the model exhibits paramagnetic behavior with a renormalized small effective moment at low temperature. This peculiar remnant paramagnetism is not obtained in the Kitaev-Heisenberg model without the additional interaction. Our results renew the picture of the octamer state and provide a scenario for the intrinsic paramagnetic behavior recently observed in a muon spin rotation experiment [K. M. Kojima et al., Phys. Rev. Lett. 112, 087203 (2014)].

preprint2014arXiv

Spontaneous Parity Breaking in Spin-Orbital Coupled Systems

Effects of spontaneous parity breaking by charge, spin, and orbital orders are investigated in a two-band Hubbard model on a honeycomb lattice. This is a minimal model in which the inter-orbital hopping, atomic spin-orbit coupling, and strong electron correlation give rise to fascinating properties, such as the magnetoelectric effects, quantum spin Hall effect, and spin or valley splitting in the band structure. We perform the symmetry analysis of possible broken-parity states and the mean-field analysis of their competition. We find that the model at 1/4 filling exhibits a spin-orbital composite ordered state and a charge ordered state, in addition to a paramagnetic quantum spin-Hall insulator. We show that the composite ordered phase exhibits two types of magnetoelectric responses. The charge ordered state shows spin splitting in the band structure, while the topological nature varies depending on electron correlations.

preprint2014arXiv

Surface and interface effects on a magnetic Chern insulator

Effects of an open surface on a magnetic Chern insulator are investigated in comparison with those of an interface to a capping magnetic layer. In magnets, an open surface often perturbs the magnetic order by a reconstruction of the magnetic moment directions near the surface. On the other hand, in topological insulators, it leads to the formation of topologically protected surface states. These two contrasting effects may coexist in magnetic Chern insulators, which give rise to nontrivial surface reconstruction. For instance, the chiral edge current is largely enhanced by the edge reconstruction in a two-dimensional magnetic Chern insulator realized in a quarter-filled Kondo lattice model on a triangular lattice. We here show that the edge reconstruction can be described semiquantitatively by a simple junction model between the bulk topological magnetic state and a ferromagnetic capping layer. We further clarify how the chiral edge current is affected by the magnetic structure in the capping layer. Our results indicate that the topological edge state can be controlled magnetically through the junctions.

preprint2014arXiv

Topological semimetal-to-insulator phase transition between noncollinear and noncoplanar multiple-Q states on a square-to-triangular lattice

Noncollinear and noncoplanar magnetic orders lead to unusual electronic structures and transport properties. We here investigate two types of multiple-Q magnetically ordered states and a topological phase transition between them in two dimensions. One is a coplanar but noncollinear double-Q state on a square lattice, which is a semimetal accommodating massless Dirac electrons. The other is a noncoplanar triple-Q state on a triangular lattice, which is a Chern insulator showing the quantum anomalous Hall effect. We discuss the peculiar electronic structures in these two multiple-Q states in a unified way on the basis of the Kondo lattice model, which suggests a quantum phase transition between the two states in a continuous change of lattice geometry between the square and triangular lattices. We systematically examine the possibility of such a transition by using the mean-field approximation for the ground state of the periodic Anderson model. After clarifying the parameter region in which the double-Q (triple-Q) state is stabilized on the square (triangular) lattice, we show that a continuous topological phase transition indeed takes place between the double-Q Dirac semimetal and the triple-Q Chern insulator on the square-to-triangular lattice. The nature of the transition is discussed by the topologically-protected edge states as well as the bulk magnetic and electronic properties. The results indicate that unusual critical phenomena may occur at finite temperature related with multiple-Q chiral spin-liquid states.

preprint2014arXiv

Toroidal order in a partially disordered state on a layered triangular lattice: implication to UNi4B

A partial disorder on a layered triangular lattice is theoretically investigated from a viewpoint of toroidal ordering and magnetoelectric effects. We consider an extended periodic Anderson model including a site-dependent antisymmetric spin-orbit coupling between conduction and localized electrons. We show that, by the mean-field approximation, the model exhibits a coplanar vortex-lattice-type magnetic order as observed in a hexagonal uranium compound UNi4B, in the parameter region with intermediate hybridization and electron correlation. This peculiar state accommodates a toroidal order, which leads to the linear magnetoelectric effect. We discuss the implications of our results to UNi4B, focusing on the possible source of the site-dependent antisymmetric spin-orbit coupling.

preprint2014arXiv

Toroidal order in metals without local inversion symmetry

Toroidal order, given by a composite of electric and magnetic orders, manifests itself not only in the peculiar magnetism but also in anomalous transport and magnetoelectric effect. We report our theoretical results on the influence and stability of a toroidal order in metals on the basis of a microscopic model. We consider an effective single-band Hubbard-type model with a site-dependent antisymmetric spin-orbit coupling, which is derived from a four-band tight-binding model including the atomic spin-orbit coupling, off-site hybridizations between orbitals with different parities, and odd-parity crystalline electric field. For this single-band model on a layered honeycomb lattice, we investigate the electronic structure, magnetotransport, and magnetoelectric effect in the toroidal ordered state with a vortex-like magnetic structure. The ferroic order of the microscopic toroidal moments acts as an effective gauge field for electrons, which modulates the electronic band structure with a shift of the band bottom in the momentum space. In addition, the site-dependent antisymmetric spin-orbit coupling gives rise to highly anisotropic Hall responses. The most salient feature is two different types of magnetoelectric responses: one is a magnetic order with net toroidal magnetization induced by an electric current perpendicular to the planes, and the other is a uniform transverse magnetization induced by an electric current within the planes. We examine the ground state of the effective model by the mean-field approximation, and show that the toroidal order is stabilized by strong electron correlations at low electron density. We also discuss the temperature dependence of the magnetoelectric effects associated with spontaneous toroidal ordering. Implications to experiments are also presented.

preprint2014arXiv

Vaporization of Kitaev spin liquids

Quantum spin liquid is an exotic quantum state of matter in magnets. This state is a spin analogue of the liquid helium which does not solidify down to the lowest temperature due to strong quantum fluctuations. In conventional fluids, liquid and gas possess the same symmetry and adiabatically connect to each other by bypassing the critical end point. We find that the situation is qualitatively different in quantum spin liquids realizing in a three-dimensional Kitaev model; both gapless and gapped quantum spin liquid phases at low temperatures are always distinguished from the high-temperature paramagnet (spin gas) by a phase transition. The results challenge common belief that the absence of thermodynamic singularity down to the lowest temperature is a symptom of a quantum spin liquid.

preprint2013arXiv

Charge Order in a Two-Dimensional Kondo Lattice Model

The possibility of charge order is theoretically examined for the Kondo lattice model in two dimensions, which does not include bare repulsive interactions. Using two complementary numerical methods, we find that charge order appears at quarter filling in an intermediate Kondo coupling region. The charge ordered ground state is an insulator exhibiting an antiferromagnetic order at charge-poor sites, while the paramagnetic charge-ordered state at finite temperatures is metallic with pseudogap behavior. We confirm that the stability of charge order is closely related with the local Kondo-singlet formation at charge-rich sites. Our results settle the controversy on charge order in the Kondo lattice model in realistic spatial dimensions.

preprint2013arXiv

Charge Order with a Noncoplanar Triple-Q Magnetic Order on a Cubic Lattice

Ground state of the periodic Anderson model on a cubic lattice is investigated by mean-field calculations, with focusing on the possibility of charge ordering. We show that three different charge-ordered phases appear at 3/2 filling, which are accompanied by ferromagnetic, noncoplanar triple-Q, and ferrimagnetic order, respectively. The triple-Q magnetic order is likely stabilized by emergent geometrical frustration in the charge-ordered phase with the ordering wave vector (pi,pi,pi). We show that the nearest-neighbor hopping between localized electrons stabilizes the triple-Q phase relative to the other two phases. We also discuss the electronic structure in the triple-Q charge-ordered phase.

preprint2013arXiv

Effect of Quantum Spin Fluctuation on Scalar Chiral Ordering in the Kondo Lattice Model on a Triangular Lattice

Spin scalar chiral ordering gives rise to nontrivial topological characters and peculiar transport properties. We here examine how quantum spin fluctuations affect the spin scalar chiral ordering in itinerant electron systems. We take the Kondo lattice model on a triangular lattice, and perform the linear spin wave analysis in the Chern insulator phases with spin scalar chiral ordering obtained in the case that the localized spins are classical. We find that, although the quantum fluctuation destabilizes the spin scalar chiral phase at 3/4 filling that originates from the perfect nesting of Fermi surface, it retains the phase at 1/4 filling that is induced by an effective positive biquadratic interaction. The reduction of the ordered magnetic moment by the zero-point quantum fluctuation is considerably small, compared with those in spin-only systems. The results suggest that the Chern insulator at 1/4 filling remains robust under quantum fluctuations.

preprint2013arXiv

Invariant Energy Levels and Flat Band Engineering in a Kondo Lattice Model on Geometrically Frustrated Lattices

We show the existence of invariant energy levels in a Kondo lattice model on an isolated complete graph, such as a triangle and a tetrahedron. These energy levels always have fixed eigenenergies $t \pm J/2$, irrespective of the configuration of localized moments ($t$ is the transfer integral of conduction electrons and $J$ is the spin-charge coupling constant). We also extend the analysis to geometrically frustrated lattices by using the complete graphs as basic building blocks. We show that the construction rule for the invariant energy levels leads to the existence condition of localized states, if the model is defined on the triangle-based line graphs, such as a kagome lattice. We further propose a procedure of engineering isolated flat bands with broken time-reversal symmetry, which are separated from other dispersive bands with finite energy gaps.

preprint2013arXiv

Linear Spin Wave Analysis for General Magnetic Orders in the Kondo Lattice Model

We extend the formulation of the spin wave theory for the Kondo lattice model, which was mainly used for the ferromagnetic metallic state, to general magnetic orders including complex noncollinear and noncoplanar orders. The 1/S expansion is reformulated in the matrix form depending on the size of the magnetic unit cell. The noncollinearity and noncoplanarity of the localized moments are properly taken into account by the matrix elements of the para-unitary matrix used in the diagonalization of the Bogoliubov-de Gennes type Hamiltonian for magnons. We apply the formulation within the linear spin wave approximation to a typical noncollinear case, the $120^{\circ}$ N{é}el order on a triangular lattice at half filling. We calculate the magnon excitation spectrum and the quantum correction to the magnitude of ordered moments as functions of the strength of the Hund's-rule coupling, $J_{\rm H}/t$. We find that the magnon excitation shows softening at $J_{\rm H}/t \simeq 2.9$, which indicates that the $120^{\circ}$ order is destabilized for smaller $J_{\rm H}/t$. On the other hand, we show that the $120^{\circ}$ order is stable in the entire range of $J_{\rm H}/t \gtrsim 2.9$, and, in the limit of $J_{\rm H}/t \to \infty$, the form of the spin wave spectrum approaches that for the antiferromagnetic Heisenberg model, while the bandwidth is proportional to $t^2/J_{\rm H}$. The reduction of the ordered moment is smaller than that for the spin-only model, except in the vicinity of the softening.

preprint2013arXiv

Loop Liquid in an Ising-Spin Kondo Lattice Model on a Kagome Lattice

Phase diagram of an Ising-spin Kondo lattice model on a kagome lattice is investigated by a Monte Carlo simulation. We find that the system exhibits a peculiar ferrimagnetic state at a finite temperature, in which each triangle is in a two-up one-down spin configuration but the spin correlation does not develop any superstructure. We call this state the loop liquid, as it is characterized by the emergent degree of freedom, self-avoiding up-spin loops. We elucidate that the system shows phase transitions from the loop liquid to ferrimagnetically ordered states and a crossover to a partially ferromagnetic state by changing the electron density and temperature. These can be viewed as crystallization and cohesion of the loops, respectively. We demonstrate that the loop formation is observed in the optical conductivity as a characteristic resonant peak.

preprint2013arXiv

Monte Carlo Study of an Effective Ising Model for the Spin-ice type Kondo Lattice Model

An effective Ising model for the spin-ice type Kondo lattice model is investigated by the classical Monte Carlo simulation. We clarify the magnetic phase diagram with four phases: ice-ferro, ice-(0,0,2π), 32-sublattice, and all-in/all-out ordered states. The result well reproduces the phase diagram of the Kondo lattice model studied previously [H. Ishizuka et al.: J. Phys. Soc. Jpn. 81 (2012) 113706], which suggests that the RKKY interactions up to third neighbors are sufficient to describe the magnetic properties of the itinerant electron model. We discuss the peculiar nature of phase transitions: the suppression of the critical temperatures down to zero temperature between two ice phases and the presence of the tricritical points on the phase boundary of the 32-sublattice ordered state.

preprint2013arXiv

On the Stability of Quantum Hall Kagome-ice Insulator

The stability of the kagome-ice insulator against the three dimensionality and the canting angle of the spin anisotropy axes is studied numerically using the exact diagonalization method. The kagome-ice insulator, which was recently found at 2/3 electron filling in the double-exchange model on a kagome lattice, is a Chern insulator showing a quantized Hall response despite the absence of magnetic long-range order. We find that this peculiar state remains stable when the lattice is extended to three dimensions in an anisotropic pyrochlore lattice in the weak interlayer coupling region. On the other hand, we show that the kagome-ice insulator is also stable against the change of the canting angle. Furthermore, we find another kagome-ice insulator at 1/3 filling in a different range of the canting angle. Our results indicate that the quantum anomalous Hall response is widely observed in the variants of the kagome-ice double-exchange systems.

preprint2013arXiv

Polynomial Expansion Monte Carlo Study of Frustrated Itinerant Electron Systems: Application to a Spin-ice type Kondo Lattice Model on a Pyrochlore Lattice

We present the benchmark of the polynomial expansion Monte Carlo method to a Kondo lattice model with classical localized spins on a geometrically frustrated lattice. The method enables to reduce the calculation amount by using the Chebyshev polynomial expansion of the density of states compared to a conventional Monte Carlo technique based on the exact diagonalization of the fermion Hamiltonian matrix. Further reduction is brought by a real-space truncation of the vector-matrix operations. We apply the method to the model with spin-ice type Ising spins on a three-dimensional pyrochlore lattice, and carefully examine the convergence in terms of the order of polynomials and the truncation distance. We find that, in a wide range of electron density at a relatively weak Kondo coupling compared to the noninteracting bandwidth, the results by the polynomial expansion method show good convergence to those by the conventional method within reasonable numbers of polynomials. This enables us to study the systems up to 4x8^3 = 2048 sites, while the previous study by the conventional method was limited to 4x4^3 = 256 sites. On the other hand, the real-space truncation is not helpful in reducing the calculation amount for the system sizes that we reached, as the sufficient convergence is obtained when most of the sites are involved within the truncation distance. The necessary truncation distance, however, appears not to show significant system size dependence, suggesting that the truncation method becomes efficient for larger system sizes.

preprint2013arXiv

Spin-orbital frustration in pyrochlores A2Mo2O7

Electronic and magnetic properties of molybdenum pyrochlores A2Mo2O7 are studied by the fully relativistic density-functional theory plus on-site repulsion (U) method, with focusing on the spin-glass insulating material Y2Mo2O7. We find that the system exhibits peculiar competition in energy between different magnetic states in the large-U insulating region. The magnetic competition cannot be explained by the conventional picture based on the geometrical frustration of isotropic Heisenberg antiferromagnetic exchange interactions. Through an analysis by using a generalized spin model, we find that the effective spin interactions are distinct from the simple Heisenberg form and strongly anisotropic in spin space. We also reveal that they give rise to keen competition between antiferromagnetic and ferromagnetic states. The complex form of the magnetic interactions indicates a crucial role of the orbital degree of freedom. Analyzing a three-orbital Hubbard model, we clarify that the magnetic competition is tightly connected with orbital frustration in the 4d2 electronic configuration through the spin-orbital interplay. The results challenge the conventional picture of the spin-glass behavior that attributes the origin to the geometrical frustration of purely antiferromagnetic exchange interactions.

preprint2013arXiv

Spontaneous Spatial Inversion Symmetry Breaking and Spin Hall Effect in a Spin-ice Double-exchange Model on a Pyrochlore Lattice

A formation of tetrahedral spin clusters is discovered by Monte Carlo simulation for a spin-ice type double-exchange model on a pyrochlore lattice. The spin-cluster phase is magnetically disordered, but breaks spatial inversion symmetry spontaneously by developing noncoplanar four-spin molecules periodically on the pyrochlore lattice. We find that the system exhibits a nonzero spin Hall conductivity in the spin-cluster phase. The result suggests that an intersite-multipole order induces the unconventional spin Hall state without the spin-orbit interaction.

preprint2013arXiv

Thermally-induced Phases in an Ising Kondo Lattice Model on a Triangular Lattice: Partial Disorder and Kosterlitz-Thouless State

Magnetic and electronic properties of a Kondo lattice model with Ising localized spins are studied on an isotropic triangular lattice. By using Monte Carlo simulation, we present that the model shows a rich phase diagram with four dominant states: two-sublattice stripe, three-sublattice ferrimganetic, partially disordered, and Kosterlitz-Thouless like quasi-long-range ordered states. Among them, the partially disordered state and Kosterlitz-Thouless like state are intermediate phases induced by thermal fluctuations in the phase competing regime; they are present only at finite temperatures and eventually taken over by another phases as the temperature is further lowered. Although the Kosterlitz-Thouless like state was found also in triangular Ising antiferromagnets with further-neighbor interactions, the partially disordered state has not been reported in the localized spin only models in two dimensions. Interestingly, the partially disordered phase is also peculiar in the charge degree of freedom of itinerant electrons; it is insulating and accompanied by charge disproportionation. From a combined analysis of a mean-field calculation of the band structure and Monte Carlo simulation, we conclude that the partial disorder in the present model is stabilized by the Slater mechanism.

preprint2012arXiv

Dirac half-metal in a triangular ferrimagnet

An idea is proposed for realizing a fully spin-polarized Dirac semimetal in frustrated itinerant magnets. We show that itinerant electrons on a triangular lattice exhibit the Dirac cone dispersion with half-metallic behavior in the presence of a three-sublattice ferrimagnetic order. The Dirac nodes have the same structure as those of graphene. By variational calculation and Monte Carlo simulation, we demonstrate that the ferrimagnetic order with the Dirac node spontaneously emerges in a simple Kondo lattice model with Ising anisotropy. The realization will be beneficial for spintronics as a candidate for spin-current generator.

preprint2012arXiv

Ground-state phase diagram of the Kondo lattice model on triangular-to-kagome lattices

We investigate the ground-state phase diagram of the Kondo lattice model with classical localized spins on triangular-to-kagome lattices by using a variational calculation. We identify the parameter regions where a four-sublattice noncoplanar order is stable with a finite spin scalar chirality while changing the lattice structure from triangular to kagome continuously. Although the noncoplanar spin states appear in a wide range of parameters, the spin configurations on the kagome network become coplanar as approaching the kagome lattice; eventually, the scalar chirality vanishes for the kagome lattice model.

preprint2012arXiv

Hidden Multiple-spin Interactions as an Origin of Spin Scalar Chiral Order in Frustrated Kondo Lattice Models

We reveal the significance of kinetic-driven multiple-spin interactions hidden in geometrically-frustrated Kondo lattice models. Carefully examining the perturbation in terms of the spin-charge coupling up to the fourth order, we find that a positive biquadratic interaction is critically enhanced and plays a crucial role on stabilizing a spin scalar chiral order near 1/4 filling in a triangular lattice case. This is a generalized Kohn anomaly, appearing only when the second-order perturbation is inefficient because of the degeneracy under frustration. The mechanism is potentially common to frustrated spin-charge coupled systems, leading to emergence of unusual magnetic orders.

preprint2012arXiv

Magnetic Order and Charge Disproportionation in a Spin-ice type Kondo Lattice Model: Large Scale Monte Carlo Study

Phase diagram of a spin-ice type Kondo lattice model, potentially relevant to metallic pyrochlore oxides, is obtained by the Monte Carlo simulation implementing the polynomial expansion technique up to the system size with 2048 sites. We identified a new 32-sublattice magnetic phase with concomitant charge disproportionation, along with other phases such as two-in two-out and all-in/all-out orders. The spin and charge pattern can be switched by external magnetic field to a different one accompanied by a half magnetization plateau.

preprint2012arXiv

Non-Kondo mechanism for resistivity minimum in spin ice conduction systems

We present a mechanism of resistivity minimum in conduction electron systems coupled with localized moments, which is distinguished from the Kondo effect. Instead of the spin-flip process in the Kondo effect, electrons are elastically scattered by local spin correlations which evolve in a particular way under geometrical frustration as decreasing temperature. This is demonstrated by the cellular dynamical mean-field theory for a spin-ice type Kondo lattice model on a pyrochlore lattice. Peculiar temperature dependences of the resistivity, specific heat, and magnetic susceptibility in the non-Kondo mechanism are compared with the experimental data in metallic Ir pyrochlore oxides.

preprint2012arXiv

Partial Disorder and Metal-Insulator Transition in the Periodic Anderson Model on a Triangular Lattice

Ground state of the periodic Anderson model on a triangular lattice is systematically investigated by the mean-field approximation. We found that the model exhibits two different types of partially disordered states: one is at half filling and the other is at other commensurate fillings. In the latter case, the kinetic energy is lowered by forming an extensive network involving both magnetic and nonmagnetic sites, in sharp contrast to the former case in which the nonmagnetic sites are rather isolated. This spatially extended nature of nonmagnetic sites yields a metallic partially-disordered state by hole doping. We discuss the mechanism of the metal-insulator transition by the change of electronic structure.

preprint2012arXiv

Partial disorder in an Ising-spin Kondo lattice model on a triangular lattice

Phase diagram of an Ising-spin Kondo lattice model on a triangular lattice near 1/3-filling is investigated by Monte Carlo simulation. We identify a partially disordered phase with coexistence of magnetic order and paramagnetic moments, which was unstable in two-dimensional Ising models with localized spins only. The partial disorder emerges in the competing regime between a twosublattice stripe phase and three-sublattice ferrimagnetic phase, at finite temperatures above an electronic phase separation. The peculiar magnetic structure accompanies a charge order and develops a gap in the electronic structure. The results manifest a crucial role of the nonperturbative interplay between spin and charge degrees of freedom in stabilizing the partial disorder.

preprint2012arXiv

Quantum Anomalous Hall Effect in Kagome Ice

An anomalous Hall insulator without magnetic long-range ordering is theoretically reported in the absence of the relativistic spin-orbit coupling. It is realized in itinerant electrons coupled with the Ising spins on a <111> kagome plane of pyrochlore spin ice in applied magnetic field. We find that the kagome-ice type local spin correlation in the magnetization plateau state opens a charge gap without magnetic ordering, which results in quantization of the Hall conductivity. By Monte Carlo simulation, we identify the anomalous Hall insulating region in the magnetic phase diagram, in addition to another anomalous Hall insulator in a fully-saturated state.

preprint2012arXiv

Thermally-induced magnetic phases in an Ising spin Kondo lattice model on a kagome lattice at 1/3-filling

Numerical investigation on the thermodynamic properties of an Ising spin Kondo lattice model on a kagome lattice is reported. By using Monte Carlo simulation, we investigated the magnetic phases at 1/3-filling. We identified two successive transitions from high-temperature paramagnetic state to a Kosterlitz-Thouless-like phase in an intermediate temperature range and to a partially disordered phase at a lower temperature. The partially disordered state is characterized by coexistence of antiferromagnetic hexagons and paramagnetic sites with period $\sqrt3 \times \sqrt3$. We compare the results with those for the triangular lattice case.

preprint2011arXiv

Application of polynomial-expansion Monte Carlo method to a spin-ice Kondo lattice model

We present the results of Monte Carlo simulation for a Kondo lattice model in which itinerant electrons interact with Ising spins with spin-ice type easy-axis anisotropy on a pyrochlore lattice. We demonstrate the efficiency of the truncated polynomial expansion algorithm, which enables a large scale simulation, in comparison with a conventional algorithm using the exact diagonalization. Computing the sublattice magnetization, we show the convergence of the data with increasing the number of polynomials and truncation distance.

preprint2011arXiv

Carrier doping to a partially disordered state in the periodic Anderson model on a triangular lattice

We investigate the effect of hole and electron doping to half-filling in the periodic Anderson model on a triangular lattice by the Hartree-Fock approximation at zero temperature. At half-filling, the system exhibits a partially disordered insulating state, in which a collinear antiferromagnetic order on an unfrustrated honeycomb subnetwork coexists with nonmagnetic state at the remaining sites. We find that the carrier doping destabilizes the partially disordered state, resulting in a phase separation to a doped metallic state with different magnetic order. The partially disordered state is restricted to the half-filled insulating case, while its metallic counterpart is obtained as a metastable state in a narrow electron doped region.

preprint2011arXiv

Cluster-size dependence in cellular dynamical mean-field theory

We examine the cluster-size dependence of the cellular dynamical mean-field theory (CDMFT) applied to the two-dimensional Hubbard model. Employing the continuous-time quantum Monte Carlo method as the solver for the effective cluster model, we obtain CDMFT solutions for 4-, 8-, 12-, and 16-site clusters at a low temperature. Comparing various periodization schemes, which are used to construct the infinite-lattice quantities from the cluster results, we find that the cumulant periodization yields the fastest convergence for the hole-doped Mott insulator where the most severe size dependence is expected. We also find that the convergence is much faster around (0,0) and (pi/2,pi/2) than around (pi,0) and (pi,pi). The cumulant-periodized self-energy seems to be close to its thermodynamic limit already for a 16-site cluster in the range of parameters studied. The 4-site results remarkably agree well with the 16-site results, indicating that the previous studies based on the 4-site cluster capture the essence of the physics of doped Mott insulators.

preprint2011arXiv

Critical property of spin-glass transition in a bond-disordered classical antiferromagnetic Heisenberg model with a biquadratic interaction

Motivated by puzzling spin-glass behaviors observed in many pyrochlore-based magnets, effects of magnetoelastic coupling to local lattice distortions were recently studied by the authors for a bond-disordered antiferromagnet on a pyrochlore lattice [Phys. Rev. Lett. 107, 047204 (2011)]. Here, we extend the analyses with focusing on the critical property of the spin-glass transition which occurs concomitantly with a nematic transition. Finite-size scaling analyses are performed up to a larger system size with 8192 spins to estimate the transition temperature and critical exponents. The exponents are compared with those in the absence of the magnetoelastic coupling and with those for the canonical spin-glass systems. We also discuss the temperature dependence of the specific heat in comparison with that in canonical spin-glass systems as well as an experimental result.

preprint2011arXiv

Loop algorithm for classical antiferromagnetic Heisenberg models with biquadratic interactions

Monte Carlo simulation using the standard single-spin flip algorithm often fails to sample over the entire configuration space at low temperatures for frustrated spin systems. A typical example is a class of spin-ice type Ising models. In this case, the difficulty can be avoided by introducing a global-flip algorithm, the loop algorithm. Similar difficulty is encountered in O(3) Heisenberg models in the presence of biquadratic interaction. The loop algorithm, however, is not straightforwardly applied to this case, since the system does not have a priori spin-anisotropy axis for constructing the loops. We propose an extension of the loop algorithm to the bilinear-biquadratic models. The efficiency is tested for three different ways to flip spins on a loop in Monte Carlo simulation. We show that the most efficient method depends on the strength of the biquadratic interaction.

preprint2011arXiv

Mean-field Study of Charge, Spin, and Orbital Orderings in Triangular-lattice Compounds ANiO2 (A=Na, Li, Ag)

We present our theoretical results on the ground states in layered triangular-lattice compounds ANiO2 (A=Na, Li, Ag). To describe the interplay between charge, spin, orbital, and lattice degrees of freedom in these materials, we study a doubly-degenerate Hubbard model with electron-phonon couplings by the Hartree-Fock approximation combined with the adiabatic approximation. In a weakly-correlated region, we find a metallic state accompanied by \sqroot3x\sqroot3 charge ordering. On the other hand, we obtain an insulating phase with spin-ferro and orbital-ferro ordering in a wide range from intermediate to strong correlation. These phases share many characteristics with the low-temperature states of AgNiO2 and NaNiO2, respectively. The charge-ordered metallic phase is stabilized by a compromise between Coulomb repulsions and effective attractive interactions originating from the breathing-type electronphonon coupling as well as the Hund's-rule coupling. The spin-orbital-ordered insulating phase is stabilized by the cooperative effect of electron correlations and the Jahn-Teller coupling, while the Hund's-rule coupling also plays a role in the competition with other orbital-ordered phases. The results suggest a unified way of understanding a variety of low-temperature phases in ANiO2. We also discuss a keen competition among different spin-orbital-ordered phases in relation to a puzzling behavior observed in LiNiO2.

preprint2011arXiv

Metal-insulator transition caused by the coupling to localized charge-frustrated systems under ice-rule local constraint

We report the results of our theoretical and numerical study on electronic and transport properties of fermion systems with charge frustration. We consider an extended Falicov-Kimball model in which itinerant spinless fermions interact repulsively by U with localized particles whose distribution satisfies a local constraint under geometrical frustration, the so-called ice rule. We numerically calculate the density of states, optical conductivity, and inverse participation ratio for the models on the pyrochlore, checkerboard, and kagome lattices, and discuss the nature of metal-insulator transitions at commensurate fillings. As a result, we show that the ice-rule local constraint leads to several universal features in the electronic structure; a charge gap opens at a considerably small U compared to the bandwidth, and the energy spectrum approaches a characteristic form in the large U limit, that is, the noninteracting tight-binding form in one dimension or the $δ$-functional peak. In the large U region, the itinerant fermions are confined in the macroscopically-degenerate ice-rule configurations, which consist of a bunch of one-dimensional loops: We call this insulating state the charge ice. On the other hand, transport properties are much affected by the geometry and dimensionality of lattices; e.g., the pyrochlore lattice model exhibits a transition from a metallic to the charge-ice insulating state by increasing U, while the checkerboard lattice model appears to show Anderson localization before opening a gap. Meanwhile, in the kagome lattice case, we do not obtain clear evidence of Anderson localization. Our results elucidate the universality and diversity of phase transitions to the charge-ice insulator in fully frustrated lattices.

preprint2011arXiv

Noncoplanar spin canting in lightly-doped ferromagnetic Kondo lattice model on a triangular lattice

Effect of the coupling to mobile carriers on the 120$^\circ$ antiferromagnetic state is investigated in a ferromagnetic Kondo lattice model on a frustrated triangular lattice. Using a variational calculation for various spin orderings up to a four-site unit cell, we identify the ground-state phase diagram with focusing on the lightly-doped region. We find that an electron doping from the band bottom immediately destabilizes a 120$^\circ$ coplanar antiferromagnetic order and induces a noncoplanar three-sublattice ordering accompanied by an intervening phase separation. This noncoplanar phase has an umbrella-type spin configuration with a net magnetic moment and a finite spin scalar chirality. This spin-canting state emerges in competition between the antiferromagnetic superexchange interaction and the ferromagnetic double-exchange interaction under geometrical frustration. In contrast, a hole doping from the band top retains the 120$^\circ$-ordered state up to a finite doping concentration and does not lead to a noncolpanar ordering.

preprint2011arXiv

Novel Magnetic Phases Revealed by Ultra-High Magnetic Field in the Frustrated Magnet ZnCr2O4

The Faraday rotation technique is used to map out the finite-temperature phase diagram of the prototypical frustrated magnet ZnCr2O4, in magnetic fields of up to 190 T generated by the single-turn coil method. We find evidence for a number of magnetic phase transitions, which are well-described by the theory based on spin-lattice coupling. In addition to the 1/2 plateau and a 3:1 canted phase, a 2:1:1 canted phase is found for the first time in chromium spinel oxides, which has been predicted by a theory of Penc et al. to realize in a small spin-lattice coupling limit. Both the new 2:1:1 and the 3:1 phase are regarded as the supersolid phases according to a magnetic analogy of Matsuda and Tsuneto, and Liu and Fisher.

preprint2011arXiv

Orbital degeneracy and Mott transition in Mo pyrochlore oxides

We present our theoretical results on an effective two-band double-exchange model on a pyrochlore lattice for understanding intricate phase competition in Mo pyrochlore oxides. The model includes the twofold degeneracy of $e_g'$ orbitals under trigonal field splitting, the interorbital Coulomb repulsion, the Hund's-rule coupling between itinerant $e_g'$ electrons and localized $a_{1g}$ spins, and the superexchange antiferromagnetic interaction between the $a_{1g}$ spins. By Monte Carlo simulation with treating the Coulomb repulsion at an unrestricted-type mean-field level, we obtain the low-temperature phase diagram as functions of the Coulomb repulsion and the superexchange interaction. The results include four dominant phases with characteristic spin and orbital orders and the metal-insulator transitions among them. The insulating region is characterized by a `ferro'-type orbital ordering of the $e_g'$ orbitals along the local $<111>$ axis, irrespective of the spin ordering.

preprint2011arXiv

Partial Disorder in the Periodic Anderson Model on a Triangular Lattice

We report our theoretical results on the emergence of a partially-disordered state at zero temperature and its detailed nature in the periodic Anderson model on a triangular lattice at half filling. The partially-disordered state is characterized by coexistence of a collinear antiferromagnetic order on an unfrustrated honeycomb subnetwork and nonmagnetic state at the remaining sites. This state appears with opening a charge gap between a noncollinear antiferromagnetic metal and Kondo insulator while changing the hybridization and Coulomb repulsion. We also find a characteristic crossover in the low-energy excitation spectrum as a result of coexistence of magnetic order and nonmagnetic sites. The result demonstrates that the partially-disordered state is observed distinctly even in the absence of spin anisotropy, in marked contrast to the partial Kondo screening state found in the previous study for the Kondo lattice model.

preprint2011arXiv

Quantum Monte Carlo study of the formation of molecular polarizations and the antiferroelectric ordering in squaric acid crystals

Effects of geometrical frustration and quantum fluctuation are theoretically investigated for the proton ordering in a quasi-two-dimensional hydrogen-bonded system, squaric acid crystal. We elucidate the phase diagram for an effective model, the transverse-field Ising model on a frustrated checkerboard lattice, by using quantum Monte Carlo simulation. A crossover to liquidlike paraelectric state with well-developed molecular polarizations is identified, distinguishably from long-range ordering. Emergence of long-range order from the liquidlike state exhibits peculiar aspects originating from the lifting of quasi-macroscopic degeneracy, such as colossal enhancement of the transition temperature and a vanishingly small anomaly in the specific heat.

preprint2011arXiv

Spin-glass transition in bond-disordered Heisenberg antiferromagnets coupled with local lattice distortions on a pyrochlore lattice

Motivated by puzzling characteristics of spin-glass transitions widely observed in pyrochlore-based frustrated materials, we investigate effects of coupling to local lattice distortions in a bond-disordered antiferromagnet on the pyrochlore lattice by extensive Monte Carlo simulations. We show that the spin-glass transition temperature $\TSG$ is largely enhanced by the spin-lattice coupling, and furthermore, becomes almost independent of $Δ$ in a wide range of the disorder strength $Δ$. The critical property of the spin glass transition is indistinguishable from that of the canonical Heisenberg spin glass in the entire range of $Δ$. These peculiar behaviors are ascribed to a modification of the degenerate manifold from continuous to semidiscrete one by the spin-lattice coupling.

preprint2011arXiv

Theory of pseudogap and superconductivity in doped Mott insulators

Underdoped Mott insulators provide us with a challenge of many-body physics. Recent renewed understanding is discussed in terms of the evolution of pole and zero structure of the single-particle Green's function. Pseudogap as well as Fermi arc/pocket structure in the underdoped cuprates is well reproduced from the recent cluster extension of the dynamical mean-field theory. Emergent coexisting zeros and poles set the underdoped Mott insulator apart from the Fermi liquid, separated by topological transitions. Cofermion proposed as a generalization of exciton in the slave-boson framework accounts for the origin of the zero surface formation. The cofermion-quasiparticle hybridization gap offers a natural understanding of the pseudogap and various unusual Mottness. Furthermore the cofermion offers a novel pairing mechanism, where the cofermion has two roles: It reinforces the Cooper pair as a pair partner of the quasiparticle and acts as a glue as well. It provides a strong insight for solving the puzzle found in the dichotomy of the gap structure.

preprint2011arXiv

Trimer Formation and Metal-Insulator Transition in Orbital Degenerate Systems on a Triangular Lattice

As a prototypical self-organization in the system with orbital degeneracy, we theoretically investigate trimer formation on a triangular lattice, as observed in LiVO2. From the analysis of an effective spin-orbital coupled model in the strong correlation limit, we show that the previously-proposed orbital-ordered trimer state is not the lowest-energy state for a finite Hund's-rule coupling. Instead, exploring the ground state in a wide range of parameters for a multiorbital Hubbard model, we find an instability toward a different orbital-ordered trimer state in the intermediately correlated regime in the presence of trigonal crystal field. The trimer phase appears in the competing region among a paramagnetic metal, band insulator, and Mott insulator. The underlying mechanism is nesting instability of the Fermi surface by a synergetic effect of Coulomb interactions and trigonal-field splitting. The results are compared with experiments in triangularlattice compounds, LiVX2 (X=O, S, Se) and NaVO2.

preprint2011arXiv

Variational Monte Carlo Study of the Kondo Necklace Model with Geometrical Frustration

We investigate the ground state of the Kondo necklace model on geometrically-frustrated lattices by the variational Monte Carlo simulation. To explore the possibility of a partially-ordered phase, we employ an extension of the Yosida-type wave function as a variational state, which can describe a coexistence of spin-singlet formation due to the Kondo coupling and magnetic ordering by the Ruderman-Kittel-Kasuya-Yosida interaction. We show the benchmark of the numerical simulation to demonstrate the high precision brought by the optimization of a large number of variational parameters. We discuss the ground-state phase diagram for the model on the kagome lattice in comparison with that for the triangular-lattice case.

preprint2010arXiv

Chirality-driven mass enhancement in the kagome Hubbard model

We investigate the quasiparticle-mass enhancement in the Hubbard model on the frustrated kagome lattice by using a cluster extension of the dynamical mean-field theory. By analyzing the cluster density matrix, we find a hierarchy of energy scale among charge, spin, and chirality degrees of freedom. Large amount of entropy associated with the chirality is released at a much lower temperature than other energy scales for spin and charge fluctuations, leading to a sharp peak in the specific heat and the single-particle spectrum. The results manifest a generic mechanism of mass enhancement driven by an emergent composite degree of freedom under geometrical frustration.

preprint2010arXiv

Chirality-spin separation in the Hubbard model on the kagome lattice

Effect of geometrical frustration in strongly-correlated metallic region is studied for the Hubbard model on the kagome lattice at half filling by a cluster extension of the dynamical mean-field theory combined with a continuous-time auxiliary-field quantum Monte Carlo method. We find that the electron correlation enhances the spin chirality in both vector and scalar channels. The chirality grows as decreasing temperature and exhibits a peak at a low temperature, indicating a new energy scale under strong correlation. The peak temperature is considerably lower than that for the local spin moment, namely, the characteristic temperatures for the chirality and the local moment are well separated. This is a signature of separation between spin and chiral degrees of freedom in the correlated metallic regime under geometrical frustration.

preprint2010arXiv

Cluster dynamical mean-field study of the Hubbard model on a 3D frustrated hyperkagome lattice

We study the Hubbard model on a geometrically-frustrated hyperkagome lattice by a cluster extension of the dynamical mean field theory. We calculate the temperature ($T$) dependences of the specific heat ($C$) and the spin-lattice relaxation time ($T_1$) in correlated metallic region. $C/T$ shows a peak at $T=T_{p1}$ and rapidly decreases as $T->0$. On the other hand, $1/T_1T$ has a peak at a higher temperature $T_{p2}$ than $T_{p1}$, and largely decreases below $T_{p2}$, followed by the Korringa law $1/T_1 propto T$ as $T->0$. Both peak temperatures are suppressed and the peaks become sharper as electron correlation is increased. These behaviors originate from strong renormalization of the energy scales in the peculiar electronic structure in this frustrated system; a pseudo-gap like feature, the van-Hove singularity, and the flat band. The results are discussed in comparison with the experimental data in the hyperkagome material, Na$_4$Ir$_3$O$_8$.

preprint2010arXiv

Doped high-Tc cuprate superconductors elucidated in the light of zeros and poles of electronic Green's function

We study electronic structure of hole- and electron-doped Mott insulators in the two-dimensional Hubbard model to reach a unified picture for the normal state of cuprate high-Tc superconductors. By using a cluster extension of the dynamical mean-field theory, we demonstrate that structure of coexisting zeros and poles of the single-particle Green's function holds the key to understand Mott physics in the underdoped region. We show evidence for the emergence of non-Fermi-liquid phase caused by the topological quantum phase transition of Fermi surface by analyzing low-energy charge dynamics. The spectra calculated in a wide range of energy and momentum reproduce various anomalous properties observed in experiments for the high-Tc cuprates. Our results reveal that the pseudogap in hole-doped cuprates has a d-wave-like structure only below the Fermi level, while it retains non-d-wave structure with a fully opened gap above the Fermi energy even in the nodal direction due to a zero surface extending over the entire Brillouin zone. In addition to the non-d-wave pseudogap, the present comprehensive identifications of the spectral asymmetry as to the Fermi energy, the Fermi arc, and the back-bending behavior of the dispersion, waterfall, and low-energy kink, in agreement with the experimental anomalies of the cuprates, do not support that these originate from (the precursors of) symmetry breakings such as the preformed pairing and the d-density wave fluctuations, but support that they are direct consequences of the proximity to the Mott insulator. Several possible experiments are further proposed to prove or disprove our zero mechanism.

preprint2010arXiv

Electronic phase separation in the pyrochlore double-exchange model

Electronic phase separation and related inhomogeneity is ubiquitously seen in strongly-correlated systems. A typical example is found between ferromagnetic metal and antiferromagnetic insulator in CMR manganese oxides. Here we demonstrate that the geometrical frustration brings distinctive aspects into the phase separation phenomena. From Monte Carlo simulation and a simple energy comparison for the pyrochlore double-exchange model, we show that such phase separation takes place between ferromagnetic and paramagnetic metals. We discuss the relevance of our results to a spin-glassy metallic phase found in Mo pyrochlore oxides under external pressure.

preprint2010arXiv

Ferromagnetic transition in the double-exchange model on the pyrochlore lattice

The double-exchange model, which has been extensively studied in the context of colossal magneto-resistance in perovskite manganese oxides, is known to exhibit a ferromagnetic metallic state at low temperatures because of the interplay between localized moments and itinerant electrons through the Hund's-rule coupling. Here we investigate numerically the ferromagnetic transition in the double-exchange model defined on the frustrated pyrochlore lattice as a simple model for ferromagnetic pyrochlore oxides. We demonstrate that the finite-size corrections are largely reduced by implementing averages over the twisted boundary conditions in the Monte Carlo simulation, which enables to estimate the ferromagnetic transition temperature in relatively small size clusters. The estimate is compared with that for the non-frustrated cubic lattice system.

preprint2010arXiv

Loop algorithm for classical Heisenberg models with spin-ice type degeneracy

In many frustrated Ising models, a single-spin flip dynamics is frozen out at low temperatures compared to the dominant interaction energy scale because of the discrete "multiple valley" structure of degenerate ground-state manifold. This makes it difficult to study low-temperature physics of these frustrated systems by using Monte Carlo simulation with the standard single-spin flip algorithm. A typical example is the so-called spin ice model, frustrated ferromagnets on the pyrochlore lattice. The difficulty can be avoided by a global-flip algorithm, the loop algorithm, that enables to sample over the entire discrete manifold and to investigate low-temperature properties. We extend the loop algorithm to Heisenberg spin systems with strong easy-axis anisotropy in which the ground-state manifold is continuous but still retains the spin-ice type degeneracy. We examine different ways of loop flips and compare their efficiency. The extended loop algorithm is applied to the following two models, a Heisenberg antiferromagnet with easy-axis anisotropy along the z axis, and a Heisenberg spin ice model with the local <111> easy-axis anisotropy. For both models, we demonstrate high efficiency of our loop algorithm by revealing the low-temperature properties which were hard to access by the standard single-spin flip algorithm. For the former model, we examine the possibility of order-from-disorder and critically check its absence. For the latter model, we elucidate a gas-liquid-solid transition, namely, crossover or phase transition among paramagnet, spin-ice liquid, and ferromagnetically-ordered ice-rule state.

preprint2010arXiv

Nematic, vector-multipole, and plateau-liquid states in the classical O(3) pyrochlore antiferromagnet with biquadratic interactions in applied magnetic field

The classical bilinear-biquadratic nearest-neighbor Heisenberg antiferromagnet on the pyrochlore lattice does not exhibit conventional Neel-type magnetic order at any temperature or magnetic field. Instead spin correlations decay algebraically over length scales r ~ \sqrt{T}, behavior characteristic of a Coulomb phase arising from a strong local constraint. Despite this, its thermodynamic properties remain largely unchanged if Neel order is restored by the addition of a degeneracy-lifting perturbation, e.g., further neighbor interactions. Here we show how these apparent contradictions can be resolved by a proper understanding of way in which long-range Neel order emerges out of well-formed local correlations, and identify nematic and vector-multipole orders hidden in the different Coulomb phases of the model. So far as experiment is concerned, our results suggest that where long range interactions are unimportant, the magnetic properties of Cr spinels which exhibit half-magnetization plateaux may be largely independent of the type of magnetic order present.

preprint2010arXiv

Nonequilibrium relaxation study of the anisotropic antiferromagnetic Heisenberg model on the triangular lattice

Effect of exchange anisotropy on the relaxation time of spin and vector chirality is studied for the antiferromagnetic classical Heisenberg model on the triangular lattice by using the nonequilibrium relaxation Monte Carlo method. We identify the Berezinskii-Kosterlitz-Thouless (BKT) transition and the chiral transition in a wide range of the anisotropy, even for very small anisotropy of $\sim$ 0.01%. As the anisotropy decreases, both the critical temperatures steeply decrease, while the BKT critical region becomes divergently wide. We elucidate a sharp "V shape" of the phase diagram around the isotropic Heisenberg point, which suggests that the isotropic case is exceptionally singular and the associated Z vortex transition will be isolated from the BKT and chiral transitions. We discuss the relevance of our results to peculiar behavior of the spin relaxation time observed experimentally in triangular antiferromagnets.

preprint2010arXiv

Novel spin-liquid states in the frustrated Heisenberg antiferromagnet on the honeycomb lattice

Recent experiment on a honeycomb-lattice Heisenberg antiferromagnet (AF) Bi$_3$Mn$_4$O$_{12}$(NO$_3$) revealed a novel spin-liquid-like behavior down to low temperature, which was ascribed to the frustration effect due to the competition between the AF nearest- and next-nearest-neighbor interactions $J_1$ and $J_2$. Motivated by the experiment, we study the ordering of the $J_1$ -$J_2$ frustrated classical Heisenberg AF on a honeycomb lattice both by a low-temperature expansion and a Monte Carlo simulation. The model has been known to possess a massive degeneracy of the ground state, which, however, might be lifted due to thermal fluctuations leading to a unique ordered state, the effect known as 'order-by-disorder'. We find that the model exhibits an intriguing ordering behavior, particularly near the AF phase boundary. The energy scale of the order-by-disorder is suppressed there down to extremely low temperatures, giving rise to exotic spin-liquid states like a "ring-liquid" or a "pancake-liquid" state accompanied by the characteristic spin structure factor and the field-induced antiferromagnetism. We argue that the recent experimental data are explicable if the system is in such exotic spin-liquid states.

preprint2010arXiv

Partial Kondo screening in frustrated Kondo lattice systems

We investigate the effect of geometrical frustration on the competition between the Kondo coupling and the Ruderman-Kittel-Kasuya-Yosida interaction in Kondo lattice systems. By variational Monte Carlo simulations, we reveal an emergent quantum phase with partial ordering in which the frustration is relieved by forming a magnetic order on a sublattice and leaving the rest in the Kondo screening with spin-singlet formation. The role of quantum fluctuations and spin-charge interplay is elucidated.

preprint2010arXiv

Phase competition in the double-exchange model on the frustrated pyrochlore lattice

Competition between the ferromagnetic double-exchange interaction and the super-exchange antiferromagnetic interaction is theoretically studied in the presence of geometrical frustration. As increasing the super-exchange interaction, the ferromagnetic metal becomes unstable, and is taken over by a cooperative paramagnetic metal, in sharp contrast with a discontinuous transition to the antiferromagnetic insulator in the absence of frustration. In the critical region, the system exhibits a peculiar temperature-independent behavior with highly incoherent transport, suggesting a large residual entropy at low temperatures. We discuss the relevance of the results to the pressure-induced behaviors in Mo pyrochlore oxides [S. Iguchi et al., Phys. Rev. Lett. 102, 136407 (2009)].

preprint2010arXiv

Quantum melting of charge ice and non-Fermi-liquid behavior: An exact solution for the extended Falicov-Kimball model in the ice-rule limit

An exact solution is obtained for a model of itinerant electrons coupled to ice-rule variables on the tetrahedron Husimi cactus, an analogue of the Bethe lattice of corner-sharing tetrahedra. It reveals a quantum critical point with the emergence of non-Fermi-liquid behavior in melting of the "charge ice" insulator. The electronic structure is compared with the numerical results for the pyrochlore-lattice model to elucidate the physics of electron systems interacting with the tetrahedron ice rule.

preprint2010arXiv

Spin chirality ordering and anomalous Hall effect in the ferromagnetic Kondo lattice model on the triangular lattice

We investigate the ground-state phase diagram for the ferromagnetic Kondo lattice model on a triangular lattice by a variational calculation for various spin orderings up to a four-site unit cell. We find that a noncoplanar four-sublattice ordering with a finite scalar spin chirality emerges at and around 1/4 filling, in addition to the 3/4-filled case, which was predicted to be induced by the perfect nesting of the Fermi surface [I. Martin and C. D. Batista: Phys. Rev. Lett. {\bf 101} (2008) 156402]. The 1/4-filling phase is stable in a wider range of parameters than the 3/4-filling one, and includes a large region of gapped insulating state characterized by a Chern number. We also compute the Hall conductivity in the chiral-ordered phases.

preprint2006arXiv

Spin Singlet State in Heptamers Emerging in Spinel Oxide AlV$_2$O$_4$

We present our theoretical results on the electronic state in vanadium spinel oxide AlV$_2$O$_4$. The material is a mixed-valent system with the average valence V$^{2.5+}$, and V cations constitute a frustrated pyrochlore structure. It shows a structural transition at $\sim 700$ K, leading to the formation of seven V-sites clusters -- heptamers. We study the electronic state of the heptamer by explicitly taking account of orbital degree of freedom as well as electron correlations. We show that the ground state of the heptamer for realistic parameters becomes spin-singlet because of strong $σ$-type bonding states of $t_{2g}$ orbitals. The temperature dependence of the magnetic susceptibility in experiments is naturally understood by this singlet formation in heptamers. Our results indicate that in AlV$_2$O$_4$ orbital physics is relevant to stabilize a cluster-type singlet state instead of a previously proposed charge-ordered state with valence skipping.

preprint2005arXiv

Colossal magnetoresistance and quenched disorder in manganese oxides

We give an overview on several recent topics of colossal magnetoresistive manganites in both experiments and theories, focusing on the effect of quenched disorder. The disorder is intrinsically involved since the compounds are solid solutions, and its importance has been pointed out in several experiments of transport and magnetic properties. Recent progress in the experimental control of the strength of disorder is also reviewed. Theoretically, the effect of the disorder has been explored within the framework of the double-exchange mechanism. Several efforts to understand the phase diagram and the electronic properties are reviewed. We also briefly discuss a recent topic on the effect of disorder on competing phases and the origin of colossal magnetoresistance.

preprint2004arXiv

Phase transitions in spin-orbital coupled model for pyroxene titanium oxides

We study the competing phases and the phase transition phenomena in an effective spin-orbital coupled model derived for pyroxene titanium oxides ATiSi2O6 (A=Na, Li). Using the mean-field-type analysis and the numerical quantum transfer matrix method, we show that the model exhibits two different ordered states, the spin-dimer and orbital-ferro state and the spin-ferro and orbital-antiferro state. The transition between two phases is driven by the relative strength of the Hund's-rule coupling to the onsite Coulomb repulsion and/or by the external magnetic field. The ground-state phase diagram is determined. There is a keen competition between orbital and spin degrees of freedom in the multicritical regime, which causes large fluctuations and significantly affects finite-temperature properties in the paramagnetic phase.

preprint2004arXiv

Spin Frustration and Orbital Order in Vanadium Spinels

We present the results of our theoretical study on the effects of geometrical frustration and the interplay between spin and orbital degrees of freedom in vanadium spinel oxides $A$V$_2$O$_4$ ($A$ = Zn, Mg or Cd). Introducing an effective spin-orbital-lattice coupled model in the strong correlation limit and performing Monte Carlo simulation for the model, we propose a reduced spin Hamiltonian in the orbital ordered phase to capture the stabilization mechanism of the antiferromagnetic order. Orbital order drastically reduces spin frustration by introducing spatial anisotropy in the spin exchange interactions, and the reduced spin model can be regarded as weakly-coupled one-dimensional antiferromagnetic chains. The critical exponent estimated by finite-size scaling analysis shows that the magnetic transition belongs to the three-dimensional Heisenberg universality class. Frustration remaining in the mean-field level is reduced by thermal fluctuations to stabilize a collinear ordering.

preprint2000arXiv

A Monte Carlo Method for Fermion Systems Coupled with Classical Degrees of Freedom

A new Monte Carlo method is proposed for fermion systems interacting with classical degrees of freedom. To obtain a weight for each Monte Carlo sample with a fixed configuration of classical variables, the moment expansion of the density of states by Chebyshev polynomials is applied instead of the direct diagonalization of the fermion Hamiltonian. This reduces a cpu time to scale as $O(N_{\rm dim}^{2} \log N_{\rm dim})$ compared to $O(N_{\rm dim}^{3})$ for the diagonalization in the conventional technique; $N_{\rm dim}$ is the dimension of the Hamiltonian. Another advantage of this method is that parallel computation with high efficiency is possible. These significantly save total cpu times of Monte Carlo calculations because the calculation of a Monte Carlo weight is the bottleneck part. The method is applied to the double-exchange model as an example. The benchmark results show that it is possible to make a systematic investigation using a system-size scaling even in three dimensions within a realistic cpu timescale.

Yukitoshi Motome

What is connected

Connect this record

See the researcher in context

Building this map preview

111 published item(s)

Asymmetric modulation of Majorana excitation spectra and nonreciprocal thermal transport in the Kitaev spin liquid under a staggered magnetic field

Feasibility of Kitaev quantum spin liquids in ultracold polar molecules

Magnetic field-temperature phase diagrams for multiple-$Q$ magnetic orderings: Exact steepest descent approach to long-range interacting spin systems

Magnetic Hedgehog Lattice in a Centrosymmetric Cubic Metal

Phase degree of freedom and topology in multiple-$Q$ spin textures

Spin dynamics in the square-lattice cupola system Ba(TiO)Cu$_4$(PO$_4$)$_4$

Electrical conductivity in helical and conical magnetic states

Noncoplanar multiple-$Q$ spin textures by itinerant frustration: Effects of single-ion anisotropy and bond-dependent anisotropy

Phase shift, ellipticity, angle, and topological number in skyrmion lattices

Anomalous Hall effect in $κ$-type organic antiferromagnets

Antiferromagnet-semiconductor van der Waals heterostructures: interlayer interplay of exciton with magnetic ordering

Magnetic hedgehog lattices in noncentrosymmetric metals

Majorana-Magnon Crossover by a Magnetic Field in the Kitaev Model: Continuous-Time Quantum Monte Carlo Study

Optical Hall response in spin-orbit coupled metals: Comparative study of magnetic cluster monopole, quadrupole, and toroidal orders

Perovskite as a spin current generator

Phase transitions between helices, vortices, and hedgehogs driven by spatial anisotropy in chiral magnets

Square skyrmion crystal in centrosymmetric itinerant magnets

The range of non-Kitaev terms and fractional particles in RuCl$_3$

Thermodynamic and Transport Properties in Disordered Kitaev Models

Thermodynamic classification of three-dimensional Kitaev spin liquids

Tracing Monopoles and Anti-monopoles in a Magnetic Hedgehog Lattice

Antiferromagnetic Kitaev Interactions in Polar Spin-Orbit Mott Insulators

Chiral spin liquids with crystalline $\mathbb{Z}_2$ gauge order in a three-dimensional Kitaev model

Asymmetric Magnon Excitation by Spontaneous Toroidal Ordering

Emergent spin-valley-orbital physics by spontaneous parity breaking

Engineering chiral density waves and topological band structures by multiple-$Q$ superpositions of collinear up-up-down-down orders

Fractional Spin Fluctuation as a Precursor of Quantum Spin Liquids: Majorana Dynamical Mean-Field Study for the Kitaev Model

Topological Insulators from Electronic Superstructures

A strong coupling theory for electron-mediated interactions in double-exchange models

Analysis of Charge-spin-orbital Fluctuations by Ab Initio Calculation and Random Phase Approximation: Application to Non-coplanar Antiferromagnet Cd2Os2O7

Charge-spin-orbital Fluctuations in Mixed Valence Spinels: Comparative Study of AlV2O4 and LiV2O4

Magnetic three states of matter: A quantum Monte Carlo study of spin liquids

Phase Diagram of the Kitaev-type Model on a Decorated Honeycomb Lattice in the Isolated Dimer Limit

Spin correlation and Majorana spectrum in chiral spin liquids in a decorated-honeycomb Kitaev model

Spontaneous formation of kagome network and Dirac half-semimetal on a triangular lattice

Spontaneous multipole ordering by local parity mixing

Theory of Valence Transition in BiNiO$_3$

Thermal fractionalization of quantum spins in a Kitaev model: $T$-linear specific heat and coherent transport of Majorana fermions

Thermodynamics of Chiral Spin Liquids with Abelian and Non-Abelian Anyons

Vortex Crystals with Chiral Stripes in Itinerant Magnets

3D Dirac Electrons on a Cubic Lattice with Noncoplanar Multiple-Q Order

Effect of magnetoelastic coupling on spin-glass behavior in Heisenberg pyrochlore antiferromagnets with bond disorder

Entanglement Spectrum in Cluster Dynamical Mean-Field Theory

Low-energy Majorana states in spin-liquid transitions in a three-dimensional Kitaev model

Magnetic Field Effect in One-Dimensional Charge Ordering Systems

Multiple-Q Instability by (d-2)-dimensional Connections of Fermi Surfaces

Quantum spin Hall effect in a two-orbital model on a honeycomb lattice

Reconstruction of chiral edge states in a magnetic Chern insulator

Spin-orbit coupling in octamers in a spinel sulfide CuIr$_2$S$_4$: Competition between spin-singlet and quadrupolar states, and its relevance to remnant paramagnetism

Spontaneous Parity Breaking in Spin-Orbital Coupled Systems

Surface and interface effects on a magnetic Chern insulator

Topological semimetal-to-insulator phase transition between noncollinear and noncoplanar multiple-Q states on a square-to-triangular lattice

Toroidal order in a partially disordered state on a layered triangular lattice: implication to UNi4B

Toroidal order in metals without local inversion symmetry

Vaporization of Kitaev spin liquids

Charge Order in a Two-Dimensional Kondo Lattice Model

Charge Order with a Noncoplanar Triple-Q Magnetic Order on a Cubic Lattice

Effect of Quantum Spin Fluctuation on Scalar Chiral Ordering in the Kondo Lattice Model on a Triangular Lattice

Invariant Energy Levels and Flat Band Engineering in a Kondo Lattice Model on Geometrically Frustrated Lattices

Linear Spin Wave Analysis for General Magnetic Orders in the Kondo Lattice Model

Loop Liquid in an Ising-Spin Kondo Lattice Model on a Kagome Lattice

Monte Carlo Study of an Effective Ising Model for the Spin-ice type Kondo Lattice Model

On the Stability of Quantum Hall Kagome-ice Insulator

Polynomial Expansion Monte Carlo Study of Frustrated Itinerant Electron Systems: Application to a Spin-ice type Kondo Lattice Model on a Pyrochlore Lattice

Spin-orbital frustration in pyrochlores A2Mo2O7

Spontaneous Spatial Inversion Symmetry Breaking and Spin Hall Effect in a Spin-ice Double-exchange Model on a Pyrochlore Lattice

Thermally-induced Phases in an Ising Kondo Lattice Model on a Triangular Lattice: Partial Disorder and Kosterlitz-Thouless State

Dirac half-metal in a triangular ferrimagnet

Ground-state phase diagram of the Kondo lattice model on triangular-to-kagome lattices

Hidden Multiple-spin Interactions as an Origin of Spin Scalar Chiral Order in Frustrated Kondo Lattice Models

Magnetic Order and Charge Disproportionation in a Spin-ice type Kondo Lattice Model: Large Scale Monte Carlo Study

Non-Kondo mechanism for resistivity minimum in spin ice conduction systems

Partial Disorder and Metal-Insulator Transition in the Periodic Anderson Model on a Triangular Lattice

Partial disorder in an Ising-spin Kondo lattice model on a triangular lattice