In 5$d^2$ Mott insulators with strong spin-orbit coupling, the lowest pseudospin states form a non-Kramers doublet, which carries quadrupolar and octupolar moments. A family of double-perovskites where magnetic ions form a face-centered cubic (FCC) l
attice, was suggested to unveil an octupolar order offering a rare example in d-orbital systems. The proposed order requires a ferromagnetic octupolar interaction, since the antiferromagnetic (AFM) Ising model is highly frustrated on the FCC lattice. A microscopic model was recently derived for various lattices: for an edge sharing octahedra geometry, AFM Ising octupolar and bond-dependent quadrupolar interactions were found when only dominant inter- and intra-orbital hopping integrals are taken into account. Here we investigate all possible intra- and inter-orbital exchange processes and report that interference of two intra-orbital exchanges generates a ferromagnetic octupolar interaction. Applying the strong-coupling expansion results together with tight binding parameters obtained by density functional theory, we estimate the exchange interactions for the Osmium double-perovskites, Ba$_2$BOsO$_6$ (B = Mg, Ca). Using classical Monte-Carlo simulations, we show these double-perovskites exhibit type-I AFM quadrupolar order followed by an intriguing partial quadrupole order above the transition temperature. Implications of our theory and a way to generate the octupolar order are discussed.
One-dimensional gapped phases that avoid any symmetry breaking have drawn enduring attention. In this paper, we study such phases in a bond-alternating spin-1 $K$-$Gamma$ chain built of a Kitaev ($K$) interaction and an off-diagonal $Gamma$ term. In
the case of isotropic bond strength, a Haldane phase, which resembles the ground state of a spin-$1$ Heisenberg chain, is identified in a wide region. A gapped Kitaev phase situated at dominant ferromagnetic and antiferromagnetic Kitaev limits is also found. The Kitaev phase has extremely short-range spin correlations and is characterized by finite $mathbb{Z}_2$-valued quantities on bonds. Its lowest entanglement spectrum is unique, in contrast to the Haldane phase whose entanglement spectrum is doubly degenerate. In addition, the Kitaev phase shows a double-peak structure in the specific heat at two different temperatures. In the pure Kitaev limit, the two peaks are representative of the development of short-range spin correlation at $T_h simeq 0.5680$ and the freezing of $mathbb{Z}_2$ quantities at $T_l simeq 0.0562$, respectively. By considering bond anisotropy, regions of Haldane phase and Kitaev phase are enlarged, accompanied by the emergence of dimerized phases and three distinct magnetically ordered states.
We have used resonant inelastic x-ray scattering to reveal optical magnons in a honeycomb lattice iridate $alpha$-Li$_{2}$IrO$_{3}$. The spectrum in the energy region 20-25 meV exhibits momentum dependence, of which energy is highest at the location
of the magnetic Bragg peak, ($textit{h}, textit{k}$) = ($pm$0.32, 0), and lowered toward (0, 0) and ($pm$1, 0). We compare our data with a linear spin-wave theory based on a generic nearest-neighbor spin model. We find that a dominant bond-directional Kitaev interaction of order 20 meV is required to explain the energy scale observed in our study. The observed excitations are understood as stemming from optical magnon modes whose intensity is modulated by a structure factor, resulting in the apparent momentum dependence. We also observed diffuse magnetic scattering arising from the short-range magnetic correlation well above $textit{T}_{N}$. In contrast to Na$_{2}$IrO$_{3}$, this diffuse scattering lacks the $C_3$ rotational symmetry of the honeycomb lattice, suggesting that the bond anisotropy is far from negligible in $alpha$-Li$_{2}$IrO$_{3}$.
The appearance of nontrivial phases in Kitaev materials exposed to an external magnetic field has recently been a subject of intensive studies. Here, we elucidate the relation between the field-induced ground states of the classical and quantum spin
models proposed for such materials, by using the infinite density matrix renormalization group (iDMRG) and the linear spin wave theory (LSWT). We consider the $K Gamma Gamma$ model, where $Gamma$ and $Gamma$ are off-diagonal spin exchanges on top of the dominant Kitaev interaction $K$. Focusing on the magnetic field along the $[111]$ direction, we explain the origin of the nematic paramagnet, which breaks the lattice-rotational symmetry and exists in an extended window of magnetic field, in the quantum model. This phenomenon can be understood as the effect of quantum order-by-disorder in the frustrated ferromagnet with a continuous manifold of degenerate ground states discovered in the corresponding classical model. We compute the dynamical spin structure factors using a matrix operator based time evolution and compare them with the predictions from LSWT. We, thus, provide predictions for future inelastic neutron scattering experiments on Kitaev materials in an external magnetic field along the $[111]$ direction. In particular, the nematic paramagnet exhibits a characteristic pseudo-Goldstone mode which results from the lifting of a continuous degeneracy via quantum fluctuations.
A family of spin-orbit coupled honeycomb Mott insulators offers a playground to search for quantum spin liquids (QSLs) via bond-dependent interactions. In candidate materials, a symmetric off-diagonal $Gamma$ term, close cousin of Kitaev interaction,
has emerged as another source of frustration that is essential for complete understanding of these systems. However, the ground state of honeycomb $Gamma$ model remains elusive, with a suggested zigzag magnetic order. Here we attempt to resolve the puzzle by perturbing the $Gamma$ region with a staggered Heisenberg interaction which favours the zigzag ordering. Despite such favour, we find a wide disordered region inclusive of the $Gamma$ limit in the phase diagram. Further, this phase exhibits a vanishing energy gap, a collapse of excitation spectrum, and a logarithmic entanglement entropy scaling on long cylinders, indicating a gapless QSL. Other quantities such as plaquette-plaquette correlation are also discussed.
The Coulombic quantum spin liquid in quantum spin ice is an exotic quantum phase of matter that emerges on the pyrochlore lattice and is currently actively searched for. Motivated by recent experiments on the Yb-based breathing pyrochlore material Ba
$_3$Yb$_2$Zn$_5$O$_{11}$, we theoretically study the phase diagram and magnetic properties of the relevant spin model. The latter takes the form of a quantum spin ice Hamiltonian on a breathing pyrochlore lattice, and we analyze the stability of the quantum spin liquid phase in the absence of the inversion symmetry which the lattice breaks explicitly at lattice sites. Using a gauge mean-field approach, we show that the quantum spin liquid occupies a finite region in parameter space. Moreover, there exists a direct quantum phase transition between the quantum spin liquid phase and featureless paramagnets, even though none of theses phases break any symmetry. At nonzero temperature, we show that breathing pyrochlores provide a much broader finite temperature spin liquid regime than their regular counterparts. We discuss the implications of the results for current experiments and make predictions for future experiments on breathing pyrochlores.
We theoretically study three-dimensional topological semimetals (TSMs) with nodal lines protected by crystalline symmetries. Compared with TSMs with point nodes, e.g., Weyl semimetals and Dirac semimetals, where the conduction and the valence bands t
ouch at discrete points, in these new TSMs the two bands cross at closed lines in the Brillouin zone. We propose two new classes of symmetry protected nodal lines in the absence and in the presence of spin-orbital coupling (SOC), respectively. In the former, we discuss nodal lines that are protected by the combination of inversion symmetry and time-reversal symmetry; yet unlike any previously studied nodal lines in the same symmetry class, each nodal line has a $Z_2$ monopole charge and can only be created (annihilated) in pairs. In the second class, with SOC, we show that a nonsymmorphic symmetry (screw axis) protects a four-band crossing nodal line in systems having both inversion and time-reversal symmetries.
We study a classical ferromagnetic Heisenberg model in the presence of Dzyaloshinskii-Moriya interactions on the corner-shared triangle lattice formed by the Mn sites of MnSi. We show that a sizable spin helicity can be obtained only when the microsc
opic Moriya vectors lie parallel to the Mn-Mn bonds. Further, such vectors are shown to produce an unpinned helical order characterized by a particular ordering wavevector magnitude but unpinned direction, dubbed partial order, at physically realizable temperatures. A consequence of such an unpinned helical ordering is that the neutron scattering intensity is sharply peaked at this wavevector magnitude. The surface formed by connecting these wavevectors is a sphere, around which the neutron scattering weight is spread. We further show that the observed neutron scattering intensity can be anisotropic along this surface and that this anisotropy is dependent on the experimentalists choice of lattice Bragg peak through a geometric factor. A neutron scattering measurement near the Bragg point ($frac{2pi}{a}$,$frac{2pi}{a}$,0) naturally leads to a highest intensity along the (1,1,0) direction consistent with the observed anisotropy in MnSi [Pfleiderer {it{et al.}} Nature {bf{427}} 227 (2004)]. A possible mechanism for pinning the helical order, and a way to distinguish an ordered and a partially ordered state in the context of neutron scattering are discussed.
It was suggested that the two consecutive metamagnetic transitions and the large residual resistivity discovered in Sr$_3$Ru$_2$O$_7$ can be understood via the nematic order and its domains in a single layer system. However, a recently reported aniso
tropy between two longitudinal resistivities induced by tilting the magnetic field away from the c-axis cannot be explained within the single layer nematic picture. To fill the gap in our understanding within the nematic order scenario, we investigate the effects of bilayer coupling and in-plane magnetic field on the electronic nematic phases in a bilayer system. We propose that the in-plane magnetic field in the bilayer system modifies the energetics of the domain formation, since it breaks the degeneracy of two different nematic orientations. Thus the system reveals a pure nematic phase with a resistivity anisotropy in the presence of an in-plane magnetic field. In addition to the nematic phase, the bilayer coupling opens a novel route to a hidden nematic phase that preserves the x-y symmetry of the Fermi surfaces.
This paper consists of two important theoretical observations on the interplay between l = 2 condensates; d-density wave (ddw), electronic nematic and d-wave superconducting states. (1) There is SO(4) invariance at a transition between the nematic an
d d-wave superconducting states. The nematic and d-wave pairing operators can be rotated into each other by pseudospin SU(2) generators, which are s-wave pairing and electron density operators. The difference between the current work and the previous O(4) symmetry at a transition between the ddw and d-wave superconducting states (Nayak 2000 Phys. Rev. B 62 R6135) is presented. (2) The nematic and ddw operators transform into each other under a unitary transformation. Thus, when a Hamiltonian is invariant under such a transformation, the two states are exactly degenerate. The competition between the nematic and ddw states in the presence of a degeneracy breaking term is discussed.
