Topological Phases emerging from Spin-Orbital Physics

We study the evolution of spin-orbital correlations in an inhomogeneous quantum system with an impurity replacing a doublon by a holon orbital degree of freedom. Spin-orbital entanglement is large when spin correlations are antiferromagnetic, while for a ferromagnetic host we obtain a pure orbital description. In this regime the orbital model can be mapped on spinless fermions and we uncover topological phases with zero energy modes at the edge or at the domain between magnetically inequivalent regions.


Introduction
Transition metal oxides are fascinating materials where several degrees of freedom (i.e., spin, orbital, charge, etc.) couple and, from a theoretical point of view, need to be treated on equal footing in order to provide reliable predictions. In undoped 3d Mott insulators large on-site Coulomb interactions localize electrons and the coupling between transition metal ions is controlled by a low-energy spin-orbital superexchange introduced first by Kugel and Khomskii [1]. As for spins, orbital degrees of freedom have a quantum character and can drive  [2] or lead to exotic novel types of magnetic order [3]. However, such cases are rare in e g systems and, typically, long range order in both spin and orbital sector develops in perovskite lattices [4], with the corresponding correlations following the Goodenough-Kanamori rules [5]. A well known example is the spin-orbital order in LaMnO 3 [6], with different energy scales for spin and orbital order [7]. However, there are numerous deviations from these rules caused either by superexchange on non-linear bonds [8], or by lattice frustration such as for instance in LiNiO 2 [9], or by spin-orbital entanglement [10], or, finally, by the presence of next nearest neighbor hopping [11]. In t 2g systems, orbital superexchange has leading contributions with SU(2) symmetry along a given cubic direction thus orbital fluctuations are much stronger [12] than in e g and a spin-orbital liquid emerging from intrinsic frustration is more likely to occur [13,14,15]. On the other hand, ordered states may be even stabilized by orbital fluctuations [16] as for instance in LaVO 3 [17] and Ca 2 RuO 4 [18] -in this latter case spin-orbit coupling also plays a role [19]. Quantum fluctuations and spin-orbital entanglement [20] are of great importance in this class of materials and may lead to novel phenomena as superconductivity in the pnictides driven by competing symmetries at orbital degeneracy [21], or spectacular topological structure of the excited states in the one-dimensional (1D) SU(2)⊗XY model [22], or, finally, dimerised phases [23].
Doping of Mott insulators adds another charge degree of freedom in spin-orbital systems and leads to several remarkable phenomena. Recently short range charge-density wave called stripe phase was reported in doped cuprates [24]. It has been suggested that the critical charge, orbital, and spin fluctuations near the quantum critical point provide the pairing interaction 2 Wojciech Brzezicki et al. [25]. As in doped cuprates, the holes doped in t 2g orbitals may be mobile due to three-site terms [26] or self-organization in stripe phases [27]. However, the formation of orbital molecules makes 1D insulating zigzag states kinetically more favorable than metallic stripes [28]. Insulating state is also found [29] when holes are confined near charge defects in Y 1−x Ca x VO 3 [30].
In contrast, neutral defects in spin-orbital systems lead to orbital dilution (with a local increase of spin to S = 3 2 ) [31] and the changes in spin-orbital order [33], or charge dilution [32] (with invariant spin S = 1 states). These phenomena are distinct from the orbital dilution in cuprates where holes remove simultaneously spin and orbital degree of freedom [34]. The t 2g systems with charge dilution are unexplored yet -they will likely play a major role in future functional materials and, possibly, in novel electronic devices. The purpose of the paper is to investigate the consequences of charge dilution in a t 2g system due to the substitution of a d 4 by a d 2 transition metal ion. Such type of doping allows to uniquely design a spin-orbital correlated environment with an orbital degree of freedom having an inequivalent charge character. Indeed, for d 2 and d 4 valence configurations, the empty orbital (i.e., holon) and the doubly occupied state (i.e., doublon) set the orbital degree of freedom, respectively. As an experimental motivation we mention, among the various emergent phenomena and the many possible hybrid oxides which could be designed, that: (i) dilute Cr doping for Ru reduces the temperature of the orthorhombic distortion, induces ferromagnetic (FM) order and anomalous negative thermal expansion in Ca 2 Ru 1−x Cr x O 4 (with 0 < x < 0.13) [35], and (ii) Mn-substituted single crystals of Sr 3 Ru 2−x Mn x O 7 rapidly drive an unusual metalinsulator transition and E-type antiferromagnetic (AF) order at low doping [36]. The theoretical search for the consequences of holon-doublon substitution is performed for a 1D ring and we analyze both spin and orbital correlations around the charge defect. We give reasons why the FM regime is well designed to search for topological aspects of the present model.

Spin-orbital physics and charge dilution
We consider a 1D ring made of d 4 transition metal ions in the insulating regime, with one d 2 charge defect, see Fig. 1. The physics of the undoped system is governed by a spin-orbital superexchange model which is equivalent, through an electron-hole transformation, to that introduced for vanadates [37]. It depends on two Kanamori's parameters: the intraorbital Coulomb element U 2 and Hund's exchange J 2 for t 2g electrons [4], responsible for the high spin states with spin S = 1 at the host d 4 ions. In the regime of strong electron interactions we obtain a spin-orbital model with spin S = 1 at every site and an orbital degree of freedom described by a pseudospin T = 1 2 . Since we work in one dimension and with t 2g orbitals, we select the cubic axis c with the active orbitals [16]: |a ≡ |yz and |b ≡ |xz .
However, the situation becomes less familiar when some of the d 4 ions are substituted by the d 2 ones. In the regime of low doping, all the bonds will be either between two d 4 ions, called host bonds, or between d 2 and d 4 ions around an impurity site -these we call hybrid bonds. The superexchange Hamiltonian for both kinds of bonds has a generic form (all the bonds are along the cubic axis c in the 1D chain) [32], where the label b = {0, h} stands for the type of bond and the operators K act in the orbital space at two sites, {i, i + 1}. These operators differ fundamentally for the hybrid and host bonds, i.e., for the host they take the U(1) symmetric form of, whereas for the hybrid bonds the symmetry is lowered, Here, τ α i are the Pauli operators describing doublon/holon fluctuating between |a and |b orbitals. They are defined by the Pauli matrices σ as where hardcore boson operators a † i and b † i create holon or doublon in the orbitals |a and |b , respectively.
Coefficient J b in the spin-orbital model (1) is a superexchange constant and is given by J 0 = 4t 2 0 /U 2 and J h = 2t 2 h /∆ (note that the excitations which provide the main contribution on the hybrid bond go in one direction only), where t b (b = 0, h) is a hopping amplitude along host-host or impurity-host bond, and U 2 is the Hubbard interaction for the host (d 4 ) atoms. ∆ is a typical excitation energy scale in the virtual process j (or charge transfer energy) given by, Here U 1 and J 1 are Hubbard and Hund's interactions at the d 2 impurity site. I e is the energy mismatch of the electronic levels at two ions; since ∆ must be positive and relatively large, this implies that I e > U i must be the largest energy scale in the system.
All these η's measure the relative strength of Hund's exchange with respect to typical excitation energyin case of host bonds it is U 2 whereas for hybrid bonds it is ∆ (7). The exact functional forms of these coefficients are complicated and will be reported elsewhere.
The properties of the host and hybrid bonds are the following: a single host's bond is always FM in spin and AF in orbital sector because of an orbital singlet which is formed on a bond [16]. This however is not stable when there are more than one bond -for a longer 1D system as for L = 8 chain considered here we find AF spin correlations for low η (8) turning FM in a high η limit. The case of a hybrid bond is much simpler: despite the complicated form of the Hamiltonian (1) it always gives AF spin correlations accompanied by FM τ z i τ z i+1 orbital correlations. Because of these intrinsic difference between host and hybrid bonds it is essential to check the ground state properties of a finite system with single impurity, see Fig. 1.
In Fig. 2 we show the ground state spin and orbital correlations obtained for a closed chain of L = 8 sites with a single d 2 impurity, see Fig. 1. The results are shown as functions of η for fixed values of J 0 , J h , η 1 , and η 2 which weakly influence the overall behavior. Due to translational invariance one finds four inequivalent bonds, see Fig. 1. There are two regimes: (i) AF with total spin S = 2 ( i S z i = 2) for η < 0.09, (ii) FM with S = 6 ( i S z i = 6) for η > 0.09 (however the hybrid impurity bonds are always AF). In the AF regime at η = 0 all the spin correlations are AF, but a level crossing occurs at η = 0.033 where the magnetic moment delocalizes from the impurity to its two neighbors, remaining nearly constant within these three sites.
Surprisingly, for increasing η < 0.09 the spin correlations between second and third neighbors of impurity become soon FM, due to spin-orbital entanglement, but the remaining spin correlations are AF. In the FM regime all the host bonds have almost saturated FM spin correlations, +1, while they tend to the classical value of −1 for increasing η on hybrid bonds, see Fig.  2(a). The orbital τ z i τ z i+1 correlations behave more regularly; they are AF for host bonds and FM for hybrid bonds in both regimes of η, see Fig. 2(b). For the offdiagonal orbital correlations we define the conventional are normalized to ±1). It turns out that τ + i τ − i+1 correlations are significant only for the host bonds and τ + i τ + i+1 only for hybrid bonds and they are always AF, see Fig. 2(c).
To investigate the spin-orbital entanglement we introduce covariances for the various correlators, (10) with σ = ±. In Fig. 3 we show the spin-orbital covariances in the AF and FM regime. One finds that both Ground state spin-orbital covariances C αβ i,i+1 on the bonds for a closed chain system of length L = 8 (Fig. 1), obtained via exact diagonalization as functions of η for orbital correlations: (a) τ z longitudinal (C zz i,i+1 ) and transverse (C +σ i,i+1 ) covariances are large in the AF regime. Moreover, as one could expect, they are much lower at higher η > 0.1 when the host spin correlations are FM, while they tend to zero as η increases, see Fig. 3. Interestingly, the transverse covariance for the hybrid bond is positive (C ++ i,i+1 > 0) in the entire regime of parameters which suggests that double orbital excitations are strong on hybrid bonds. Thus, we conclude that the factorization into spin and orbital operators is a good approximation only in the FM regime and for this case we set the spin-spin correlations as equal to ±1 for the host-impurity bonds.

Topological states in the orbital model
Factorization of spin and orbital degrees of freedom is allowed in the FM regime and leads to an effective orbital-only Hamiltonian, for the host and hybrid bonds, respectively. This purely orbital Hamiltonian can be mapped on spinless fermions by the Jordan-Wigner transformations. For symmetry reasons we find an exact relation A zz ≡ −A xx , and we also get that A xx and A yy almost compensate each other so their sum A xx + A yy = δ has a relatively small amplitude. It is however important to point out that δ = 0 because in the representation of (solid and dashed lines), with color convention as in Fig. 1; (b) magnified view of (a) for η 2 → 0; (c) orbital couplings Azz and sum Axx + Ayy (solid and dashed) for hybrid bonds, and (d) topologically relevant quantities, ImD π/2 and D 0 Dπ (solid and dashed). Parameters: Jordan-Wigner fermions δ is proportional to the pairing amplitude and its finite value can induce a topological non-trivial state. All the {A α,α } coefficients are functions of η 1 and η 2 , while we find that the dependence on η 1 is very weak. Thus, we fix η 1 = 4 (η is already fixed as η = 0.2) and we show the behavior of the A α,α coupling in Fig. 4(c). We note that η 2 = 0 is a high symmetry point where A α,α ≡ 0.5 for any η 1 .
Hence, by means of the Hartree-Fock decoupling we deal with fermion-interaction term τ z i τ z i+1 in a selfconsistent manner and we obtain the bond τ α i τ α i+1 orbital correlations for a periodic (and infinite) system with one d 2 impurity per every L = 8 sites. We find that in the present parameter regime the τ z i τ z i+1 vanish and one gets only the kinetic terms τ x i τ x i+1 and τ y i τ y i+1 . For the host bonds they are all AF while for the hybrid ones the xx correlations are AF and yy ones are FM, see Fig. 4(a). Interestingly, we obtain a discontinuous transition at η 2 = 0 + between anisotropic and isotropic phases -the difference between xx and yy correlations is triggered by any finite η 2 , see Fig. 4(b). Finally, at finite η 2 one always gets a regime with a non-trivial topological phase with respect to the Jordan-Wigner fermionic representation.
Indeed, the fermionic Hamiltonian in the momentum space is given by a matrix H k that belongs to the BDI Altland-Zirnbauer class [38]. Thus, it can have a non-trivial Z topological number. The topological invariant can be determined by looking at H k in the Topological Phases emerging from Spin-Orbital Physics eigen-basis of the chiral symmetry where it consists of two anti-diagonal blocks u k and u † k . The determinant of u k , D k ≡ det u k , is a complex number which yields a non-trivial topological number if it winds around the (0, 0) point in the complex plane as k changes from 0 to 2π. In present case this happens if: (i) the imaginary part of u π/2 is non-vanishing and (ii) the determinants D 0 and D π have opposite signs. In Fig. 4(d) we observe that indeed these conditions hold as long as η 2 > 0.
Recently, the topological phase diagram of a 1D tight-binding model of spinless electrons with an inhomogeneous distribution of pairing centers has been investigated [39]. The Hamiltonian includes inhomogeneities generated by diluted pairing centers with a given distribution profile in the unit cell of length L. For a periodic configuration with momentum k we get, with c L+1,k ≡ e ik c 1,k and {t p , ∆ p } being the nearest neighbor hopping and on-bond pairing amplitudes. There, we have found the topological invariant that can be generally expressed in terms of the physical parameters for any pairing center configuration [39].
Here, we emphasize the occurrence of edge states and present the spectra around zero energy for a closed and open system, see Figs. 5(a) and 5(b). We note that for an open system there are two zero-energy states appearing in the gap. These are Majorana end modes that arise as a consequence of the bulk-boundary correspondence in a topologically non-trivial configuration. In Figs. 5(c) and 5(d) the spatial occupation proba-bilities p i for the two zero energy states are explicitly shown in order to confirm their degree of localization on the right/left edges of the 1D chain. We also point out that the modification of the kinetic term with the inclusion of long-range hopping is expected to lead to multiple Majorana end modes both in spinless [40] and spinfull p-wave superconducting chains [41].

Discussion and Summary
In conclusion, we have studied a one-dimensional hybrid d 2 -d 4 system with a single d 2 impurity in a d 4 spin-orbital correlated host. Remarkably, the exact diagonalization analysis allows to single out regimes for which the orbitals and spins can be factorized if the host configuration is FM. By this decoupling one finds interacting orbital pseudospins exhibiting fully isotropic exchange for the host bonds and fully anisotropic for the hybrid ones. A Jordan-Wigner transformation and Hartree-Fock decoupling allow, then, to map the system on non-interacting fermions and to find topological non-trivial states. Unexpectedly, a topological nontrivial state occurs for any finite value of J 2 , i.e., the amplitude of the Hund's coupling at the host's ions. For a long chain we explicitly demonstrate that Majoranalike modes occur at the edge of the system. We argue that inhomogeneous topological patterns [42] can be achieved in the present spin-orbital scenario with Majorana modes occurring, for instance, at the boundary of the FM region if the impurities drive a magnetic configuration that has alternating FM with AF domains.