Arxiv:1308.3997V1 [Cond-Mat.Str-El] 19 Aug 2013 Ings Nw Stejh-Elr(T Feti Igemolecule
Total Page:16
File Type:pdf, Size:1020Kb
Vibronic Excitation Dynamics in Orbitally Degenerate Correlated Electron System Joji Nasu∗ and Sumio Ishihara Department of Physics, Tohoku University, Sendai 980-8578, Japan (Dated: August 12, 2018) Orbital-lattice coupled excitation dynamics in orbitally degenerate correlated systems are examined. We present a theoretical framework, where both local vibronic excitations and superexchange-type inter-site in- teractions are dealt with on an equal footing. We generalize the spin-wave approximation so as to take local vibronic states into account. Present method is valid from weak to strong Jahn-Teller coupling magnitudes. Two characteristic excitation modes coexist; a low-energy dispersive mode and high-energy multi-peak mode. These are identified as a collective vibronic mode, and Flanck-Condon excitations in a single Jahn-Teller cen- ter modified by the inter-site interactions, respectively. Present formalism covers vibronic dynamics in several orbital-lattice coupled systems. PACS numbers: 75.25.Dk, 75.30.Et,75.47.Lx I. INTRODUCTION multi-peak vibrational excitations with a broad envelop appear due to the Franck-Condon transitions. Orbital degree of freedom of an electron represents a di- Purpose in this paper is to present a theoretical framework rectional aspect of electronic wave function. It is widely rec- of vibronic excitations in orbitally degenerate correlated elec- ognized that the orbital degree of freedom influences signif- tron system; both the local vibronic excitations and the su- icantly magnetic, optical, and structural properties in corre- perexchange (SE)-type inter-site interaction between orbitals lated electron materials.1,2 A macroscopic symmetry breaking are taken into account on an equal footing. We set up a model of a degenerate orbital wave function, termed an orbital or- which consists of the SE interactions, the on-site JT coupling der, is often seen in several transition-metal compounds, rare- and the local lattice vibration. A generalized spin-wave ap- earth magnets, as well as molecular solids. A long-range or- proach where the local vibronic states are fully taken into ac- bital order is generally accompanied with a macroscopic lat- count is presented. Two characteristic excitation modes co- tice distortion which is compatible to a shape of the electronic exist; a low-energy dispersive vibronic mode interpreted as wave function. This is caused by an orbital-lattice interaction a renormalized “orbiton”, and high-energy multi-peaks origi- known as the Jahn-Teller (JT) effect in a single molecule. nating from the Franck-Condon excitation in a single JT cen- A collective orbital excitation in an orbital ordered state is ter modified by the SE interaction. The present formalism termed “orbital wave” and its quantized object is termed “or- does not only cover orbitally degenerate systems from weak biton”. This is an analogous to spin-wave excitation in a mag- to strong JT couplings, but also is applicable to several orbital- netically ordered state. Several experimental observations of lattice coupled models. orbiton by optical and Raman spectroscopies,3–7 x-ray scatter- In Sec. II, we introduce a model Hamiltonian for an orbital- ings8–12 and other experimental probes, as well as theoretical lattice coupled system. In Sec. III, a generalized spin-wave supports13–17 have been reported so far. Nonetheless, charac- approximation for vibronic excitations is presented. Before teristics of orbiton have not been revealed yet. This might be showing the detailed numerical calculations, results obtained attributed to a fact that coupling between orbiton and lattice by the present theory are compared with the results by the is not negligibly small, and an experimental assignment of or- exact diagonalization method in Sec. IV. The main part in biton is not so simple in comparison with that of magnon. this paper is Sec. V, where the detailed energy and momen- tum dependences of the vibronic excitations are presented. In The first theoretical examination of orbital wave was done by Cyrot and Lyon-Caen,18 and Komarov et al.,19 where Sec. VI, we focus on the low-energy excitation modes corre- sponding to the collective vibronic modes. In Sec. VII, results purely electronic orbital excitations as well as spin-orbital in the present E e system are compared with those in the coupled excitations were examined based on a correlated elec- ⊗ arXiv:1308.3997v1 [cond-mat.str-el] 19 Aug 2013 E b system. Section VIII is devoted to discussion and tron model. More realistic calculations of orbital wave were ⊗ 1 performed by one of the present authors and coworkers in summary. Refs. 20–22, where the lattice distortion is interpreted to be frozen. The adiabatic frozen-lattice treatment23 is justified in II. MODEL the limiting case where the orbital excitation energy is much higher than the phonon energy. A weak coupling approach for the JT effect was adopted in Ref. 24 where an anticrossing- In order to address an issue for the coupling between the type mixing between orbiton and phonon branches occurs. orbital excitation and the lattice dynamics, we introduce an Similar linear coupling between orbiton and phonon modes E e JT center at each lattice site and the SE-type interactions 25,26 ⊗ are examined in TmVO4 and related materials. On the between the nearest-neighbor (NN) eg orbitals. We adopt the other side, vibronic excitations in a single JT center have been following orbital-lattice coupled Hamiltonian, examined intensively.27 In particular, local vibronic excita- = + . (1) tions in an orbital ordered state were studied in Ref. 28, where H HJ HJT 2 The first term, J , represents the SE interactions and the sec- III. FORMULATION ond term, ,H is forthe local lattice vibrationand the JT cou- HJT pling. One of the prototypical SE interaction in an orbitally We present a formulation based on the generalized spin- degenerate magnet is the Kugel-Khomskii type spin-orbital wave (SW) approximation, where the local vibronic states 29,30 Hamiltonian which is derived from the multi-orbital Hub- are fully taken into account. We show later that this is valid bard model. Here, we focus on the orbital degree of freedom from weak to strong JT coupling regimes. A relation of the in the SE-type interaction, and consider the following spin- present formalism to the random-phase approximation (RPA) less orbital-only model, is given in Appendix A. Similar formalisms were presented in z z x x Refs. 32–36. J = JzTi Tj + JxTi Tj , (2) H − We assume that Jz Jx in J and the uniform orbital hiji X order for T z in the ground≥ state,H without loss of generality. where NN ij sites are represented by ij . The doubly- The z component of the PS operator is decomposed into the degenerate orbitals are described by the pseudo-spinh i (PS) op- z z z ordered moment and fluctuation as Ti = T + δTi where 1 † h i ′ ′ erator defined by Ti = 2 γγ′ diγ σγγ diγ , where diγ is an denotes the expectation value in the ground state. The annihilation operator for a spin-less fermion with orbital γ at Hamiltonianh· · · i in Eq. (1) is rewritten as site i, and σ are the Pauli matrices.P The eigen state of T z with z z x x MF the eigen value of +1/2 ( 1/2) corresponds to a state where = (JzδT δT + JxT T )+ . (6) − H − i j i j Hi the d 2 2 (d 2 2 ) orbital is occupied by an electron. The hiji i 3z −r x −y X X exchange constants, Jz and Jx, are set to be positive. Present formulation is able to be generalized easily to models where We define the on-site term z x other terms for the PS interactions, such as , and the MF z JT Ti Tj = h T + , (7) spin degree of freedom, are taken into account. This will be Hi − MF i Hi discussed in Sec. VIII. with The second term of the Hamiltonian is given by z hMF = zJz T , (8) = JT, (3) h i HJT Hi i where z is a coordination number. X z with The ordered moment T is determined by the following h i MF 2 2 2 way. The local Hamiltonian i is diagonalized numerically JT 1 ∂ ∂ Mω0 2 2 z H = + + (Q + Q ) under a given T . The eigen states Φn and the eigen Hi −2M ∂Q2 ∂Q2 2 iu iv h i {| i} iu iv energies En are obtained up to ( n) where the phonon- z x { } N ≥ +2A( Ti Qiu + Ti Qiv), (4) number is restricted to be less than Nph at each site. In the − present numerical calculations, we chose , which is where Q and Q represent the two vibrational modes at Nph = 80 iu iv enough to examine excitations of the present interest. The or- the i-th JT center with the E symmetry. The first two g dered moment is calculated by the ground-state wave-function terms describe the harmonic vibrations with frequency ω0 and z z Φ0 as T = Φ0 T Φ0 . This procedure is repeated un- a reduced mass M, and the third term describes the linear | i z h i h | | i MF til T converges. It is noted that i commutes with the JT coupling with a coupling constant A(> 0). For conve- h i v† v H z iπbi bi nience, we introduce the phonon coordinates for lattice vi- parity operator i = 2Ti e , and the eigen states are P brations at as follows. The creation and annihila- classified by the eigen values of i, i.e. i Φn = pn Φn . A = 0 P P | i | i tion operators for the mode γ(= u, v) phonons are defined by When pn = 1 (pn = 1), a parity of Φn is identified as γ† γ “even” (“odd”).