Magnetic and electronic properties of La MO and 3 7 7 1 possible polaron formation in hole-doped La MO 0 3 7 2 (M=Ru and Os) n a J Bin Gao, Yakui Weng, Jun-Jie Zhang, Huimin Zhang, Yang 6 2 Zhang and Shuai Dong ] E-mail: [email protected] i c Department of Physics, Southeast University, Nanjing 211189,China s - l r t m . Abstract. Oxideswith4d/5dtransitionmetalionsarephysicallyinterestingfortheir t a particularcrystalline structures aswell as the spin-orbitcoupledelectronic structures. m Recentexperimentsrevealedaseriesof4d/5dtransitionmetaloxidesR3MO7 (R: rare - earth; M: 4d/5dtransitionmetal) with unique quasi-one-dimensionalM chains. Here d n first-principles calculations have been performed to study the electronic structures o of La3OsO7 and La3RuO7. Our study confirm both of them to be Mott insulating c antiferromagnets with identical magnetic order. The reduced magnetic moments, [ whicharemuchsmallerthantheexpectedvalueforidealhigh-spinstate(3t2g orbitals 1 occupied), are attributed to the strong p d hybridization with oxygen ions, instead v − 5 of the spin-orbit coupling. The Ca-doping to La3OsO7 and La3RuO7 can not only 5 modulate the nominal carrier density but also affect the orbital order as well as the 6 local distortions. The Coulombic attraction and particular orbital order would prefer 7 0 to form polarons, which might explain the puzzling insulating behavior of doped 5d . transition metal oxides. In addition, our calculation predict that the Ca-doping can 1 0 trigger ferromagnetism in La3RuO7 but not in La3OsO7. 7 1 : Keywords: 4d/5d transition metal oxides, polaron, antiferromagnetism v i X r 1. Introduction a Transition metal oxides have attracted enormous attentions for its plethoric members, divergent properties, novel physics, andgreatimpactsonpotential applicationsbasedon correlatedelectrons. Inpastdecades, theoverwhelming balanceofinterestsweredevoted to those compounds with 3d elements, which showed high-T superconductivity, colossal C magnetoresistivity, multiferroicity, andso on[1, 2]. However, the4dand5dcounterparts were much less concerned and only in very recent years a few of them, e.g. Sr IrO , have 2 4 beenfocusedon[3,4,5]. Inprinciple, for4d/5delectrons, theelectron-electronrepulsion, e.g. Hubbard U, is much weaker due to more extended wave functions, while the spin- orbit coupling (SOC) is much stronger due to the large atom number, comparing with the 3d electrons [6]. These characters may lead to non-conventional physics in 4d/5d Magnetic andelectronicpropertiesof La MO and possiblepolaronformationin hole-dopedLa MO (M=Ru 3 7 3 7 metal oxides, e.g. p-wave superconductors, spin-orbit Mott insulator, Kitaev magnets, topological materials, and possible high-T superconductors [5, 7, 8, 9, 10]. C Till now, the most studied 4d/5d metal oxides owns quasi-two-dimensional layer structures (e.g. Sr IrO and Na IrO ) or three-dimensional structures (e.g. SrIrO and 2 4 2 3 3 SrRuO ). Recently, those 4d/5d metal oxides with quasi-one-dimensional chains have 3 also been synthesized, which may lead to unique low-dimensional physics, e.g. charge density waves, spin-Peierls transitions, and novel magnetic excitations [11, 12, 13]. For example, recent experiments reported the basic physical properties of R MO , 3 7 which owns the weberite structure, as shown in Fig. 1(a) [14]. Since here the 4d/5d electrons are mostly confined in one-dimensional chains instead of two-dimensional plane or three-dimensional framework, their electronic and magnetic structures, may be markedly different from the higher-dimensional structural 4d/5d counterparts. Given the decreased electron correlations and increased SOC of the 4d/5d electrons, the physical behavior of these compounds may also show differences comparing with quasi- one-dimensional 3d metal oxides [15]. In fact, there is rare 3d metal oxide forming the weberite R MO structure. It is therefore physically interest to study these new 3 7 systems. Taking La OsO for example, recent experimental studies reported its structural, 3 7 transport, and magnetic properties, characterized by magnetic susceptibility, x-ray diffraction, as well as neutron diffraction [16]. The corner-shared OsO octahedra 6 form chains along the [001] direction of the orthorhombic framework. The nearest- neighbor distance of Os-Os is 3.81 ˚A within a chain, but 6.75 ˚A between chains. The ◦ Os-O-Os bond angle within a chain is about 153 , implying strong octahedra tilting, which is also widely observed in other oxides. Its ground state is an antiferromagnetic (AFM) insulator. The Ca-doped La OsO was also studied. Despite the change of 3 7 nominal carrier density, surprisingly, this hope-doped system remain an insulator (or a semiconductor), violating the rigid band scenario [16]. Similar robust insulating behavior was also found in some doped iridates [17, 18], which was expected to show superconductivity after doping [19, 10]. In this work, we have performed systematic first-principles calculations to understand the electronic and magnetic properties of La OsO as well as the 3 7 isostructural La RuO . The doping effect has also been studied, which may provide 3 7 a reasonable explanation to the insulating behavior, based on the polaron forming. To ourbestknowledge, therewereveryfewtheoreticalstudiesonthesetwomaterialsbefore. Only Khalifah et al. calculated several magnetic states of La RuO [20]. Even though, 3 7 their predicted ground state (see Fig. 1(b)) seems to be inaccurate, according to our results. 2. Model & methods All following calculations were performed using the Vienna ab initio Simulation Package (VASP) based on the generalized gradient approximation (GGA) [21, 22]. The new- Magnetic andelectronicpropertiesof La MO and possiblepolaronformationin hole-dopedLa MO (M=Ru 3 7 3 7 Figure 1. Structure and possible magnetism of La3MO7 (M=Os or Ru). (a) Schematic of the weberite R3MO7 structure. (b)-(e) The possible AFM orders considered in our calculations. A and B are the chain indexes of M chains. + and − denote the spin directions of M. developed PBEsol function is adopted [23], which can improve the accurate description of crystal structure comparing with the old-fashion PBE one. The plane-wave cutoff is 550 eV and Monkhorst-Pack k-points mesh centered at Γ points is adopted. Starting from the low-temperature experimental orthorhombic (No. 63 Cmcm) structures [16, 24], the lattice constants and inner atomic positions are fully optimized till the Hellman-Feynman forces are all less than 0.01 eV/˚A. The Hubbard repulsion U (= U J) is imposed on Ru’s 4d orbitals and Os’s 5d orbitals [25]. Various values eff − of U have been tested from 0 eV to 4 eV. It is found that U (Ru)=1 eV is the best eff eff choice to reproduce the experimental structure of La RuO , while for La OsO the bare 3 7 3 7 GGA without U is the best choice. Comparing with experimental values [16, 20], eff the deviation of calculated lattice constants are only < 0.8% for La OsO , and 0.5% for 3 7 La RuO , providing a goodstart point to study other physical properties. These choices 3 7 of U are quite reasonable considering the gradually deceasing Hubbard repulsion from eff 3d to 5d. Considering the fact of heavy atoms, the relativistic SOC is also taken into consideration, comparing with those calculations without SOC. 3. Results & discussion 3.1. Undoped La MO : magnetic orders and reduced moments 3 7 First, the magnetic ground state is checked by comparing several possible magnetic orders, including ferromagnetic (FM) state, and various AFM ones (AFM I-IV as shown in Fig. 1(b-e)). For AFM I, III, and IV states, the -up-down-up-down- ordering is Magnetic andelectronicpropertiesof La MO and possiblepolaronformationin hole-dopedLa MO (M=Ru 3 7 3 7 Table 1. The energy difference (in unite of meV) for a minimal unit cell (four f.u.’s), local magnetic moment per M within the default Wigner-Seitz sphere (in unit of µB), and band gap (in unit of eV) of La3MO7. Magnetism Energy Moment Gap La OsO FM 0 0.851,0.846 0.19 3 7 AFM I -689 1.670,-1.670 0.42 AFM II -228 1.435,-1.435 0 AFM III -709 1.652,-1.652 0.54 AFM IV -714 1.651,-1.651 0.53 La RuO FM 0 1.958,1.959 0.53 3 7 AFM I -23 1.902,-1.902 0.60 AFM II -22 1.948,-1.948 0.78 AFM III -38 1.902,-1.902 0.73 AFM IV -42 1.900,-1.899 0.70 adopted within each chain, but with different coupling between chains. Taking the FM state as the energy reference, the energies of all candidates are summarized in Table 1, which suggest the AFM IV to be the possible ground state for both M=Os and Ru. By mapping the system to a classical spin model, the exchange coefficients between neighbor (within each chain and between chains, as indicated in Fig. 1(a)) spins (normalized to S =1) can be extracted as: J =72.00 meV, J =8.34 meV, and J =2.96 1 2 3 | | meV for M=Os, J =3.45 meV, J =0.66 meV, and J =0.37 meV for M=Ru. Obviously, 1 2 3 the exchanges between Os chains are quite prominent even for the nearest neighbor chains (distance up to 6.75 ˚A), implying strongly coupled AFM chains, different from the one-dimensional intuition. And these exchanges are much stronger in La OsO than 3 7 the correspondences in La RuO . These characters of La OsO are benefited from the 3 7 3 7 more extended distribution of 5d orbitals. Experimentally, theAFMtransitiontemperaturesofLa OsO aremuchhigherthan 3 7 the corresponding ones of La RuO . For La OsO , the intrachain magnetic correlation 3 7 3 7 emerges near 100 K (mainly due to J ) and the fully three-dimensional AFM ordering 1 ∼ occurs at 45 K (also determined by J and J [20]. In contrast, the signal for magnetic 2 3 ordering in La RuO appears at 17 K with short-range characters [20]. Note Ref. [20] 3 7 ∼ once predicted the ground state of La RuO to be AFM I, which is ruled out according 3 7 to our calculation. More neutron experiments are needed to refine the subtle magnetic order of La RuO . 3 7 Second, the total density of states (DOS) and atomic-projected density of states (PDOS) of La OsO are displayed in Fig. 2(a). Clearly, the system is insulating with a 3 7 band gaps of 0.53 eV, even in the pure GGA calculation. Both the topmost valence ∼ band(s) and bottommost conducting band(s) of La OsO are from Os, in particular the 3 7 t orbitals. Since the 5d orbitals have a large SOC coefficient, we also calculate the 2g DOS and PDOS with SOC enabled, which are presented in Fig. 2(b) for comparison. Magnetic andelectronicpropertiesof La MO and possiblepolaronformationin hole-dopedLa MO (M=Ru 3 7 3 7 Total DOS La O Os (a) Os 30 5 0 0 -30 -5 90 SOC-enabled (b) Os (SOC-en abled) S O 5 45 D OS/P 0 Tota l DOS La O Ru (c) Ru 0 30 5 D 0 0 -30 -5 90 SOC-enabled (d) Ru (SOC-e nabled) 5 45 0 0 -1.0 -0.5 0.0 0.5 1.0 1.5 -4 -3 -2 -1 0 1 Energy (eV) Figure 2. Density ofstate (DOS) and projecteddensity ofstate (PDOS) ofLa3MO7 (M=OsorRu). (a-b)La3OsO7;(c-d)La3RuO7. (a)and(c)SOC-disabled;(b)and(d) SOC-enabled. Rightpanels: thecorrespondingnear-Fermi-levelPDOSofanindividual M ion. However, there is no qualitative difference between the SOC-enabled and SOC-disabled calculations. The quantitative differences include: 1) a shrunk band gap 0.37 ∼ eV (SOC-enabled); 2) a slightly reduced local magnetic moment from 1.661 µ /Os B (SOC-disabled) to 1.578 µ /Os (SOC-enabled). In particular, the magnitude of orbital B moment is only 0.087 µ . Noting this local moment is obtained by integrating the B ∼ wave function within the Wigner-Seitz radius of Os (0.58 ˚A) and thus not absolutely precise. Even though, the theoretical values are still quite close to the experimental one 1.71 µ /Os [16]. Such a local moment is significantly reduced from the high-spin B ∼ expectation (3 µ /Os) of three t electrons as in Os5+ here, but agrees with recent B 2g neutron diffraction results of Os5+ in several double perovskites [26, 27, 28]. According to the PDOS (insert of Fig. 2(a-b)), every Os seems to be in the high- spin state, i.e. only spin-up electrons within the Wigner-Seitz radius. Then how to understand the reduced local moment? Above SOC-enabled calculation has ruled out SOC as the main contribution, which can only slightly affect the value of moment. Instead, the real mechanism is the covalency between Os and O. As revealed in PDOS, there exists strong hybridization between Os’s 5d and O’s 2p orbitals around the Fermi energy level, owning to the spatial extended 5d orbitals. In fact, the previous neutron study also attributed the reduced moment to the hybridization between Os and O [16]. Furthermore, the same calculations have been done for La RuO and the 3 7 DOS/PDOS are shown in Fig. 2(c-d), which are qualitatively similar to La OsO . The 3 7 local magnetic moment of Ru5+ is 1.892 µ (SOC-enabled) or 1.900 µ (SOC-disabled), B B and such a negligible difference implies an weaker SOC effect comparing with La OsO . 3 7 Magnetic andelectronicpropertiesof La MO and possiblepolaronformationin hole-dopedLa MO (M=Ru 3 7 3 7 In particular, the magnitude of orbital moment is only 0.018 µ per Ru, even lower B ∼ than that of Os. The total moment is also lower than the ideal 3 µ but higher than the B moment of Os, which is also reasonable considering the more localized distribution of 4d orbitals than 5d. The reduced moment of Ru5+ is also due to the covalency between Ru and O, as indicated in Fig. 2(c). The calculated band gap of La RuO is 0.70 eV, which is higher than the 3 7 experimental value ( 0.28 eV) extracted from transposrt [20]. This inconsistent ∼ is probably due to their polycrystalline nature of samples and the presence of small amounts of the highly insulating La O , as admitted in Ref. [20]. More measurements, 2 3 especially the optical absorption spectrum, are needed to clarify the intrinsic band gap of La RuO . 3 7 The aforementioned weak SOC effects to magnetism and band structures in La RuO and La OsO seem to contradict with the intuitive expection of strong SOC 3 7 3 7 coefficients for 4d/5d electrons. This paradox can be understood as following. Since in La MO the low-lying t orbitals are half-filled (t3 ), the Hund coupling between t 3 7 2g 2g 2g electrons will prefer the high-spin state, in which the orbit moment is mostly quenched. Then the net effect of SOC is weak even if the SOC efficiency is large. Other 5d electronic systems with own more or less electrons than t3 , e.g. Sr IrO , can active the 2g 2 4 SOC effects. 3.2. Chemical doping and polaron forming Doping is a frequently used method to tune physical properties of materials. For Mott insulators, proper doping may result in superconductivity (e.g. for cuprates) or colossal magnetoresistivity (e.g. for manganites). One of the most anticipant doping effects on 5d metal oxides is the possible superconductivity, as predicted in Sr IrO [10, 19]. 2 4 However, till now, not only the superconductivity has not been found, but also there is an unsolved debate regarding the metallicity of doped Sr IrO . Some experiments 2 4 reported the metallic transport behavior upon tiny doping and observed Fermi arcs using angle-resolved photoelectron spectroscopy (ARPES) [29, 30, 31, 32], while some others reported robust insulating (or semiconducting) behavior even upon heavy doping by element substitution and field-effect gating [18, 33]. Similarly, for La OsO , the experiment found that the Ca-doping up to 6.67%could 3 7 reduce the resistivity but the system remained insulating [16]. Then it is interesting to investigate the doping effect. In our calculation, by using one Ca to replace one La in a unit cell, i.e. 8.33% doping, the crystal structure is re-relaxed with various magnetism. Then the ground state turns to be AFM III, a little different from the original AFM IV (see Table 2 for more details). Even though, the in-chain AFM order remains robust. In contrast, when the 8.33% Ca-doping is applied to La RuO , our calculations 3 7 predict that the ground state magnetism would probably transformed from AFM IV to FM, different from above Os-based counterpart. As summarized in Table 2, no matter the lowest energy FM state or the second lowest energy AFM II state, the in- Magnetic andelectronicpropertiesof La MO and possiblepolaronformationin hole-dopedLa MO (M=Ru 3 7 3 7 Table2. Theenergydifference(inunitofmeVfor4f.u.’s)andlocalmagneticmoments of M in unit of µB of doped La11/4Ca1/4MO7. The SOC is disabled in calculations except for those items with ”+SOC”. Magnetism Energy Moment M=Os FM 0 1.690,1.694,1.707,1.709 AFM I 79 1.559,1.622,-1.565,-1.219 AFM II -63 1.647,1.645,-1.647,-1.699 AFM III -315 1.393,-1.393,-1.528,1.511 AFM III (+SOC) 1.298,-1.298,-1.442,1.449 AFM IV -285 1.257,1.617,-1.620,-1.289 M=Ru FM 0 1.753,1.765,1.797,1.803 FM (+SOC) 1.749,1.753,1.791,1.792 AFM I 363 1.896,1.912,-1.887,1.265 AFM II 9 1.726,1.735,-1.789,-1.794 AFM III 40 -1.697,1.706,1.787,-1.749 AFM IV 42 1.568,1.891,-1.892,-1.531 chain FM order is unambiguous. This result is also reasonable considering the much weaker in-chain antiferromagnetism (i.e. J ) of La RuO . Thus the antiferromagnetism 1 3 7 of La RuO should be more fragile against chemical doping. Further experiments are 3 7 needed to verify our prediction. As shown in Fig. 3, the DOS’s of doped La OsO and La RuO own finite values 3 7 3 7 at Fermi levels, implying metallic behavior, which seems to be opposite to experimental observation of doped La OsO . However, a careful analysis finds that this finite DOS at 3 7 Fermi level should be due to a technical issue of calculation. The substitution of one La by one Ca will bring one hole into the system. However, the AFM state implies at least doubly degenerate bands (spin up and spin down). Thus one hole to doubly degenerate bands always leads to half-filling, as observed in our DOS. Thus, the finite DOS at Fermi level does not guarantee metallicity of Ca-doped La OsO , while the metallicity 3 7 of Ca-doped La RuO needs experimental verification. 3 7 The SOC-enabled calculations have also been performed for the doped La MO . 3 7 However, due to the partial hole concentration ( 1/4 per M), the SOC effect is not ∼ prominent. For example, the near-Fermi-level DOS (Fig. 3(g-h)) are similar to the corresponding non-SOC ones. The local moments for the ground states are also listed in Table 2, which are only slightly lower than the original one without SOC, especially for the Ru case. The PDOS’s of Ca-doped La OsO show that the Os ions can be classified into two 3 7 types: a) two near-Ca Os’s (one spin up and one spin down); b) other two Os’s. Their PDOS’s are slightly different, and the type-b Os’s are less affected by the Ca-doping, namely the doping effect has a tendency to be localized. In contrast, the PDOS’s of Ca-doped La RuO show that almost all four Ru ions are equally effected by the Ca- 3 7 Magnetic andelectronicpropertiesof La MO and possiblepolaronformationin hole-dopedLa MO (M=Ru 3 7 3 7 (g) (h) Figure 3. DOS/PDOS and corresponding band structures around the Fermi level of La11/4Ca1/4MO7 (M=OsorRu). (a-c)La11/4Ca1/4OsO7; (d-f)La11/4Ca1/4RuO7. In (e), spinup/downbandsaredistinguishedbyblue/redcolors. (g-h)The SOC-enabled DOS/PDOS which are only slightly changed compared with the non-SOC ones. doping. To clarify the effect of doping, the charge density distribution of hole in Ca-doped La OsO is visualized in Fig. 4(a). Here we only extract the wave function of the 3 7 above-Fermi-level partial bands, which can represent the hole (half spin-up hole plus half spin-down hole). Clearly, the orbitals of hole are the d type on Os site and the xy p type on O site (the chain direction is chosen as the z axis), implying an spatially x extended wave function. According to the Slater-Koster equation [34], the lying-down d orbitalhasavery weak hoppingamplitude alongthez-axis, if notideallyzero. Thus, xy the hole will be restricted near the Ca dopant by the Coulombic interaction, leading to the semiconducting behavior of Ca-doped La OsO . 3 7 Magnetic andelectronicpropertiesof La MO and possiblepolaronformationin hole-dopedLa MO (M=Ru 3 7 3 7 (a) (b) Os Ru Ca Ca Os Ru Q <0 Q >0 3 3 (c) Os Os Os Os Ca Os Os Os Os Figure 4. Plots of charge density distribution of hole in Ca-doped La3OsO7 (a) and La3RuO7 (b). Lower: schematic of corresponding Jahn-Teller Q3 mode. In La3OsO7, Q3 is negative, leading to a low-lying dxy orbitals for electron. In La3RuO7, Q3 is positive, leading to a higher-energy d orbitals for electron. However, the orbital xy shapes of hole after Ca-doping are very similar between these two systems, suggesting the Coulombic interaction between hole and dopant to be the common driven force. (c) The spatial hole distribution in Ca-doped La3OsO7 supercells doubled along the c axes. The corresponding doping concentration is 4.17%. Similar results are obtained for Ca-doped La3RuO7 supercells (not shown). According to the plenty experience of 3d electron systems, the lattice distortions, e.g. Jahn-Teller modes, will be always activated by partially occupied t orbitals (or 2g e orbitals) to split the energy degeneration between/among orbitals. By carefully g analyzing the bond lengths of oxygen octahedra, it is easy to verify the effect of hole modulated lattice distortions. First, the breathing mode Q can be defined as 1 (l +l +l )/√3 to characterize the size of oxygen octahedral cage, where l denote the x y z O-M-O bond length along a particular axis [35, 36]. After the doping, the changes of Q are 5.307 pm for the near-Ca Os’s and 3.636 pm for other two Os’s. These 1 − − shrunk octahedral cages are due to the Coulombic attraction between positive-charged hole on Os and negative charged oxygen ions. Second, the Jahn-Teller modes Q and Q 2 3 Magnetic andelectronicpropertiesof La MO and possiblepolaronformationin hole-dopedLa MO (M=Ru 3 7 3 7 can be defined as (l l )/√2 and ( l l +2l )/√6 respectively, which can split the x y x y z − − − degeneration among triplet t orbitals or between doublet e orbitals. ForLa OsO , the 2g g 3 7 original Jahn-Teller modes Q =0.100 pm and Q = 3.938 pm. This prominent Q mode 2 3 3 − prefers the d orbital for electrons. Therefore, the d hole after doping is not driven xy xy by this preseted lattice distortion, but can only be due to the Coulombic interaction from Ca2+ since the spatial distribution of d hole is more closer to dopant (see Fig. 4). xy Then this Coulombic-driven d hole will suppress the Q mode: Q = 2.414 pm for xy 3 3 − near-Ca Os’s and 1.007 pm for other two Os’s. Meanwhile, the Q mode are enhanced: 2 − Q =0.942 pm for near-Ca Os’s and 1.152 pm for the other two Os’s. 2 For doped La RuO , the changes of Q are 4.074 pm for the near-Ca Ru’s and 3 7 1 − 3.470 pm for the other two, similar to the case of doped La OsO . For original 3 7 − La RuO , the Jahn-Teller modes are: Q =0.110 pm, Q =5.351 pm. In contrast with 3 7 2 3 La OsO , this lattice distortion dislike the d orbital (for electron) in energy. Then the 3 7 xy Coulombic-driven d hole will further enhance this positive Q mode: Q =7.574 pm xy 3 3 for the near-Ca Ru’s and 7.913 pm for the other two Ru’s. Meanwhile, the Q mode 2 are enhanced as in La OsO : Q =1.227 pm for the near-Ca Ru’s and 1.129 pm for the 3 7 2 other two. The localization of hole can be further confirmed by calculations of supercells. As shown in Fig. 4(c-e), the hole occupancies on four near-Ca Os ions are much prominent than other Os ions which stay only one u.c. space from the dopant. Similar results exist for doped La RuO . This localized hole with distorted lattices form the polaron. 3 7 According previous studies, there are magnetic polarons in manganites, which are ferromagnetic clusters of several Mn sites embedded in AFM background [37]. This scenario is quite possible for La RuO considering the ferromagetically-aligned Ru ions 3 7 as revealed in Table 2. However, for La OsO , both the result of the minimal cell 3 7 (shown in Table 2) and the calculation of doubled-supercell (along the c-axis) dislike magnetic polaron, at least the small (up to three-site) magnetic polaron. It is also reasonable considering the differences among 3d/4d/5d electrons. The 3d electrons have strong Hubbard interaction which prefers localized magnetic moments. Thus the energy gain from forming magnetic polaron is large [37]. In contrast, the 5d electrons are more spatially extended andownweaker Hubbardinteraction, which aredisadvantage toform magnetic polaron. In fact, in Ref. [37], the second nearest-neighbor hopping (which equals to the spatially extended effect) can suppress the formation of magnetic polaron, which provides a hint to understand the difference between 3d polaron and 5d polaron. Of course, the current calculations can not fully exclude the possibility of magnetic polarons with larger sizes or higher dimensional, which are beyond our computational capability considering the fact that a minimal cell of La OsO already contains 44 ions. 3 7 In short, considering above results of hole restricted by Coulombic interaction and the lattice distortions followed the particular d -orbital hole, it is reasonable to argue xy that the carriers generate by Ca doping would be mostly localized near dopant to form polaron and contribute to the semiconducting behavior. This scenario may explain the puzzle why the expected metallicity is absent in some doped 5d-Mott-insulators.