5d transition metal iridium (Ir)-based oxides, which are mostly insulators with exotic magnetic states.
Transition metal oxides, in particular, 3d or 4d perovskites, have provided diverse emergent physics that originates from the coupling of various degrees of freedom such as spin, lattice, charge, orbital, and also disorder. 5d perovskites form a distinct class because they have strong spin-orbit coupling that introduces to the system an additional energy scale that is comparable to bandwidth and Coulomb correlation. Consequent new physics includes novel Jeff = 1/2 Mott insulators, metal–insulator transitions, spin liquids, and topological insulators. After highlighting some of the phenomena appearing in the Ruddlesden–Popper iridate series Srn+1IrnO3n+1 (n = 1, 2, and ∞), we focus on the transport properties of perovskite SrIrO3. Using epitaxial thin films on various substrates, we demonstrate that metal–insulator transitions can be induced in perovskite SrIrO3 by reducing its thickness or by imposing compressive strain. The metal–insulator transition driven by thickness reduction is due to disorder, but the metal–insulator transition driven by compressive strain is accompanied by peculiar non-Fermi liquid behaviors, possibly due to the delicate interplay between correlation, disorder, and spin-orbit coupling. We examine various theoretical frameworks to understand the non-Fermi liquid physics and metal–insulator transition that occurs in SrIrO3 and offer the Mott–Anderson–Griffiths scenario as a possible solution.
- Transition metal oxides
- metal–insulator transitions
- non-Fermi liquid behaviors
Transition metal oxides (TMOs), ranging from simple binary oxides to more complex ternary or quaternary compounds, have been a subject of intense activities in condensed matter physics and materials science [1–5]. Of various transition metal oxides, perovskites and their variations are a particularly important class because they display a rich spectrum of various competing phases and physical properties. Structurally, in perovskites with the chemical formula ABO3, large rare-earth or alkali-metal cations occupy the A-sites (tetrahedral oxygen interstices) and transition metal cations occupy the B-sites (octahedral oxygen interstices). Each transition metal ion is thus surrounded by six oxygen anions to form a BO6 octahedron, and the octahedra are connected to each other three-dimensionally. This metal–oxygen network transition mostly determines the electronic properties, while the close-packed A cations and oxygen anions provide structural stability. One characteristic feature of the perovskite family is that numerous transition metal ions of different size and valence are allowed in the same structure. Along with this elemental diversity, mutually interacting quantum degrees of freedom such as lattice, spin, charge, and orbital lead to numerous emergent physical properties in perovskites. For example, the conducting properties of perovskites range from insulating to semiconducting to metallic and finally to superconducting. The electrical properties of perovskites are equally diverse; they also display various magnetic and multifunctional properties [6–14]. These unparalleled varieties in the electronic properties of perovskites hold tremendous application potential and thus provide a continuous impetus for research in the field.
The network of BO6 octahedra in transition metal perovskites strongly influences their electronic properties because d orbitals form bands with characteristic effective masses and different strengths of on-site Coulomb interactions. Indeed, diverse phenomena have been treated in terms of various relevant energy scales such as bandwidth (W) and Coulomb repulsion (U), and the relevant efforts can be termed d orbital physics [15,16]. Disorder (D), which is unavoidable in real solids also, has an important influence on the physics of transition metal oxides.
In materials science, numerous studies of 3d- and 4d-based transition metal perovskites have been conducted to find ways to exploit their functionalities that depend sensitively on structural distortion and crystal chemistry. However, perovskite studies have been limited mainly to those that contain 3d or 4d transition metal elements. In this regard, 5d perovskites form a special class because the strong spin-orbit coupling (SOC, Λ) introduces to the system an additional energy scale that is comparable to bandwidth and Coulomb correlation . For example, metal oxide compounds containing 5d iridium (e.g., Sr2IrO4, Sr3Ir2O7, SrIrO3, BaIrO3, Na2IrO3, Eu2Ir2O7, Na4Ir3O8, and Sr2GdIrO6) display a variety of emerging properties (Table 1). Of the various exotic properties of 5d perovskites, the transport properties and, in particular, the metal–insulator transitions (MIT) in iridate perovskites are the subject of the present chapter.
Simple extrapolation of the characteristics of 3d or 4d systems fails to predict the characteristics of 5d systems. Moving down the periodic table from 3d 4d 5d, the orbitals in the solids that contain the corresponding d orbitals become increasingly extended and so does the bandwidth (W3d < W4d < W5d). Along with an increase in the bandwidth, the associated on-site Coulomb repulsion would decrease correspondingly (U3d > U4d > U5d). From these considerations, one may expect higher metallicity and less magnetic instability for the materials with more extended 5d orbitals, compared with the systems containing 3d or 4d elements, because the Stoner criterion in these compounds favors paramagnetic metallic states. Surprisingly, however, the properties of the 5d-based transition metal oxides exhibit extremely rich behaviors (Table 1). The Ruddlesden–Popper (RP) series of iridium-based Srn+1IrnO3n+1 (n = 1, 2, and ∞) as well as other structural families such as pyrochlores, kagome-type lattices, honeycomb-type structures, and double perovskites are mostly insulating (few are metallic) but may have exotic physical properties, including novel Mott insulator, lattice-driven magnetoresistance, giant magnetocaloric effect, quantum criticality, charge density wave, geometrically frustrated magnetism, possible topological insulator, and Weyl semimetal [17–19]. These diverse behaviors in iridates require that a new area of physics must be considered; this is SOC (Λ3d < Λ4d < Λ5d). In iridates, the energy scales of W, U, and Λ are all comparable and thus compete with each other. In addition, D is also important, and the emergent properties are the results of strong interplay among Coulomb interaction, SOC, and disorder.
|J eff = 1/2 Mott insulator,|
|J eff = 1/2 Mott insulator,|
non-Fermi liquid, magnetism,
strong magnetic instability,
zig-zag magnetic order
non-Fermi liquid, semi-metal
spin wave excitations
Electrical conductivity is one of the properties that used to characterize or to classify solids, so a transition from a metal to insulator or vice versa as a function of a control parameter such as temperature, pressure, or magnetic field has been a topic of interest to the condensed matter physics community for several decades [20–24]. However, despite enormous efforts, understanding of the MIT at the microscopic level is still under debate. In transition metal oxides and strongly correlated systems, the MIT is often accompanied by a change in structural or magnetic symmetry. In contrast, in the presence of a sufficient amount of disorder, MIT is often not associated with any uniform ordering or change in symmetry. Perhaps the most important mechanisms that underlay MITs in correlated systems are interaction-driven or correlation-driven Mott localization, magnetic-order-driven Slater insulator, and disorder-driven Anderson localization. In perovskites that involve heavy 5d transition metal elements (e.g., Re, Os, Ir), SOC comes as another parameter and becomes comparable in strength to other relevant energy scales. All of these different parameters in the system compete with each other so their interplay stabilizes new exotic ground states. In this chapter, our foremost aim is to provide a brief description of current research on MITs that occur in Ir-based perovskites. Our emphasis will be on a model system (SrIrO3) to provide experimental evidence that the interplay among correlation, SOC, and disorder is important in achieving different MITs. We also wish to present our understanding about the observed non-Fermi liquid physics and MITs to provide insight and motivation for further activities in this rapidly developing, yet poorly understood, field of strong spin-orbit-coupled 5d-based oxide physics.
2. Metal–insulator transitions and representative types of insulators
2.1. Metal–insulator transitions
Electrical resistance R (or resistivity ρ if intrinsic quantity is used) is a key property that characterizes materials, which are typically classified as metals or insulators. A very good metal (e.g., Cu) can have an electrical resistivity as small as ρ ∼10-10 Ω cm, whereas a good electrical insulator (e.g., quartz) has a resistivity as high as ρ ∼1020 Ω cm. The causes of this huge difference in resistivity are now well understood in modern solid-state physics. The difference between the resistivities of metals and those of insulators implies that MITs may be accompanied by a large change in ρ (up to several orders of magnitude). However, many MITs of current interest are often accompanied not by a large change in magnitude but by a qualitative change in conducting behaviors.
Metals can be distinguished from insulators by their different responses of R to temperature T. Metals are defined as materials in which R decreases as T decreases (dR/dT > 0), whereas insulators are materials in which R increases as T decreases (dR/dT < 0) (Figure 1). Strictly, the ultimate distinction between the metal and the insulator can be made at absolute zero; a metal would continue to be conductive, whereas an insulator would lose its conductivity. MITs of fundamental interest change a system from a phase with dR/dT < 0 to another with dR/dT > 0, or vice versa.
The conventional theory of electrical transport was first formulated by Drude , immediately after the discovery of the electron. The semiclassical Drude theory of electronic conductivity was built on the idea of the kinetic theory of gases, considering a metal as a gas of electrons. The key concepts in that description are the mean free path ℓ (i.e., the average length that an electron travels between successive collisions) and relaxation time τ (i.e., the average time between successive collisions). According to the Drude theory, the electronic conductivity () is directly proportional to τ as , where n is the density of electrons, e is the electron charge, and m is electron mass. For good metals, ℓ ∼100 nm and τ ~10-14 s. This semiclassical Drude model is still used today as a quick way to estimate a material’s property.
Quantum mechanics has been used to clarify the transport properties of solids, and now metals, semiconductors, and insulators are classified according to the band theory of solids (Figure 2). In metals, Fermi energy EF is within the conduction band, whereas in insulators, EF separates the electronic band into an empty upper and filled lower band. For example, quartz is an insulator with a band gap Eg ∼8.9 eV, whereas Cu is a metal. Si is also an insulator but has a small Eg ∼1.1 eV and thus is called a semiconductor. Although band theory successfully describes the conducting nature of many materials, it fails to account for the behaviors of many transition metal oxides that have relatively narrow separation between conduction and valence band. Many materials with partially filled d orbitals and an odd number of electrons per lattice site, which should be metals according to band theory, are actually insulators (e.g., NiO, V2O3, and Fe3O4). The band theory considers single electron and thus may not adequately consider many-body effects such as electron–electron interaction that occur in these complex materials. In addition, D, which is unavoidable in real solids, must also be considered. Many theoretical explanations have been proposed to account for the metallic and insulating properties of so-called correlated materials; we briefly summarize these ideas in this section.
2.2. Mott insulator
In 1937, de Boer and Verwey found that many metal–oxides (e.g., NiO, MnO, and FeO) show insulating features despite the fact that these oxides have partially filled d-bands . Soon after their discovery, Peierls suggested that strong Coulomb repulsion between electrons could be the origin of the insulating behaviors . He remarked that ‘‘it is quite possible that the electrostatic interaction between the electrons prevents them from moving at all. At low temperatures the majority of the electrons are in their proper places in the ions. The minority which have happened to cross the potential barrier find therefore all the other atoms occupied, and in order to get through the lattice have to spend a long time in ions already occupied by other electrons. This needs a considerable addition of energy and so is extremely improbable at low temperatures’’ . Peierl’s speculation aroused interest in so-called strongly correlated systems. In 1949, Mott offered a theoretical explanation of how electron–electron correlation could yield an insulating state, now called a Mott Insulator .
Consider a lattice model with a single electron orbital at each site. If the electron–electron interactions are not considered, a single band would be formed from the overlap of the atomic orbitals: when two electrons, one with spin up () and the other with spin down (), occupy each site, the band becomes full. However, when two electrons occupy the same site, they would feel a large Coulomb repulsion, which can be explained using the Hubbard Hamiltonian (Eq. (1)):
where are the creation and the annihilation operators, respectively, for electrons on the site i and j with spin σ. Electron hopping between nearest neighbor sites occurs with hopping constant −t. U is the amount of energy for each pair of electrons that occupy the same lattice site and represents on-site Coulomb correlation. is the number operator and μ is the chemical potential.
The ratio between t and U determines whether the system is a metal or an insulator. When electron–electron correlation is negligible compared to hopping, t/U >> 1, the electrons tunnel between the sites without hindrance and the system is metallic. When t/U << 1, electron–electron correlation is strong, the electrons are localized due to Coulomb interaction, and the system becomes an insulator. In the Hubbard model, at half filling , the Mott–Hubbard insulating phase with one electron per site would appear.
In transition metal oxides, often oxygen p-bands remain unchanged when interaction (U) occurs (Figure 3(a)). In contrast, the d-band splits into two subbands (upper Hubbard band and lower Hubbard band), and the half-filled state becomes an insulator with the opening of a charge gap (Figure 3(b)). The Hubbard model allows two types of MIT: the filling-control metal–insulator transition (FC-MIT), which originates from variation of electron concentration or chemical potential (μ/U); and a bandwidth-control metal–insulator transition (BC-MIT), which originates from variation of hopping energy or bandwidth (t/U).
2.3. Slater insulator
Mott insulators in the previous section are electron systems with an odd number of electrons per unit cell; despite the fact that although band theory suggests that these systems would be metallic, they are insulators. Mott insulators belong to the U >> t regime in the Hubbard model. However, another aspect must be considered to explain the insulating property even in the Hubbard model. In fact, many insulating systems, especially transition metal oxides (e.g., Na2IrO3 ), often have antiferromagnetic ground states. Slater focused on this point. He proposed that formation of spin density waves (long-range magnetic order) due to electron–electron interaction may be an origin of the insulating phase itself . The difference between a Mott insulator and a Slater insulator is that for a Mott insulator, the system remains insulating even above Néel temperature TN, whereas a Slater insulator the system should be metallic above TN.
The principle behind a magnetic-order-driven Slater insulator may be easily understood from the band picture. Slater insulators belong to the t >> U regime, which is opposite to the case for Mott insulators. Suppose for simplicity that the lattice is half filled and thus on average each site holds one electron. A low-energy state for a half-filled system ensues if a periodicity doubling of the lattice occurs (i.e., the Brillouin zone is halved). At the new Brillouin zone boundary, an energy gap for charge excitations occurs and the system becomes insulating. According to Slater, such an insulating behavior is closely connected with the appearance of magnetic order at TN, for example, electrons in a bipartite lattice , i.e., one that consists of two interpenetrate sublattices A and B in such a way that the nearest neighbors of any site are members of the opposite sublattice, for example, a rock salt arrangement in a simple cubic lattice. Thus, the nearest neighbors of an electron from A are those from B, and vice versa. Because and electrons mutually repel due to Coulomb interaction, they become preferentially arranged on alternating A and B sublattice sites. Hence, spins form a spin-density wave (SDW) whose wave vector is commensurate with the lattice. Because the electrons tend to avoid each other, the potential energy increases, and this gain is balanced by a loss in kinetic energy because of localization of electrons. The spin-density wave provides a necessary periodicity-doubling potential. With an increase in temperature, thermal fluctuations affect the ordering and narrow the energy gap. Eventually, at TN, which is typically ~102 K, the ordering is destroyed and the insulating property disappears.
2.4. Anderson insulator
Defects and impurities are unavoidable in real materials, and disorder can never be neglected in reality. Thus, disorder-induced MIT is a subject of continuing interest in condensed matter physics. In 1958, Anderson initiated the field of so-called localization by arguing that sufficiently strong randomness will localize all the electronic states within a given band and that diffusion may be completely suppressed, thereby leading to an MIT at T = 0 K . In the absence of disorder, even a small amount of hopping can delocalize the electrons. However, in the presence of sufficient disorder, the hopping process only spreads an initially localized state over a finite distance; defined as the localization length ξ0. As a result, the states at the band tail become localized and sometimes even a whole band becomes localized. Mott argued that there must exist a critical energy EC called mobility edge, which separates the localized states from the extended ones (Figure 4(a)) . If EF < EC, the system is an insulator; otherwise, it is a metal.
The basic idea behind the localization phenomenon is rather straightforward. Suppose an electron propagates in a disordered medium from point A to point B. To obtain the probability that an electron reaches point B, the probabilities of all possible paths from A to B must be summed. In a disordered medium, the phases of the interference terms may be so random that on average they mutually cancel; the resultant state can be explained by the diffusion model of conductivity, but this simple argument may not be always applicable. Imagine a wave that travels from point A along a random path to point B and then goes back to A. Two possible paths are illustrated (Figure 4(b)): a randomly chosen path and the same path traversed in time reversed sense. The two paths interfere constructively and should be treated coherently as long as time reversal symmetry is preserved. Then the probability of the electron’s return to A is twice as large as it would have been if probabilities were added (first squared and then added). The enhanced back-scattering, known as weak localization, reduces the conductance between A and B by increasing the electron’s likelihood of returning to its starting point, and eventually leads to localization . At high temperatures, coherence is lost due to thermal vibrations and back-scattering effects diminish. However, at low temperatures, thermal vibrations and inelastic scattering cease and other channels through which electrons can exchange energy become interrupted, so quantum interference immobilizes electrons to induce localization transitions. In 1979, Abrahams et al. developed the phenomenological scaling theory of localization [33, 34]. In metallic systems in two dimensions, resistivity often increases at low temperatures when T is decreased. This variation is due to weak localization and was first experimentally demonstrated by Dolan and Osheroff .
In 1960, Ioffe and Regel realized that the electron mean free path ℓ can never be shorter than the lattice spacing a, because at ℓ ≤ a, coherent quasi-particle motion would vanish and the system would be an insulator . In 1972, Mott proposed a similar argument and formed a criterion ℓmin ≈ a for minimum metallic conductivity . This criterion was named the Mott–Ioffe–Regel (MIR) limit for resistivity saturation [38, 39]. Unfortunately, the MIR limit is not universally observed and disagrees with the now-widely accepted scaling theory of localization. Nevertheless, the MIR limit is still quoted frequently and is now generalized to more complex media; the criteria ranging from kFℓmin ≈ 1 through ℓmin ≈ a to kFℓmin ≈ 2π, where the wave vector kF = 2π/λ. In some materials MIT occurs at the MIR limit. For two-dimensional materials, kFℓ is related to the ratio of quantum resistance h/e2 to sheet resistance Rsheet as .
3. Recent results on 5d perovskite iridates
5d perovskites, in particular iridates, show many new interesting phenomena; examples include novel insulating states, exotic magnetism, spin-liquid behaviors, Weyl semimetals, and topological insulators. In this regard, RP series Srn+1IrnO3n+1 have been widely investigated and are probably the most widely studied iridate system. In this section, we review some of the recent results on 5d perovskite iridates.
The term perovskite was named for the discovery of CaTiO3 in honor of Count Lev Alexander Von Perovski, a Russian mineralogist. An ideal perovskite structure has an ABX3 stoichiometry. Although most of the perovskite are oxides (X = Oxygen), other forms like fluorides, halides, and sulfides also exist in literature. Perovskite ABO3 and, more generally, perovskite families including its variations accept a large number of transition metal elements of various size and valence into the B-sites (e.g., , so the variety of transition metal oxides is unlimited: examples include cuprates, manganites, ruthenates, nickelates, titanates, and recently iridates [41, 42]. Along with this elemental diversity, mutually interacting quantum degrees of freedom such as lattice, spin, charge, and orbital leads to numerous emergent physical properties in perovskite families [1–14]. The perovskite families of current interest, iridates, can be classified as follows:
General formula of ABO3 where A or B are metal cations and O as anions. B cations are surrounded by six oxygen anions forming BO6 octahedra (e.g., SrIrO3) (Figure 5).
General formula of A2BB’O6 (ABO3+AB’O3), where the BO6 and B′O6 octahedra are alternatively arranged in two sublattices. For double perovskites, alkaline or rare earth ions occupy the A-sites, while the B-sites represent transition metal ions (e.g., Sr2RuIrO6 constitutes alternative unit cells of SrRuO3 and SrIrO3).
Layered perovskites consist of stacked two-dimensional slabs of the ABO3 layer. Three subcategories have been recognized.
The general formula of Ruddlesden–Popper (RP) series is An+1BnO3n+1, where n represents the number of octahedral layers in the repeating unit and can be visualized as repeated stacking of AO(ABO3)n [43, 44]. n = 1 corresponds to one BO6 octahedron (Sr2IrO4), and n = 2 to two BO6 octahedra (Sr3Ir2O7) (Figure 5). The differentiating characteristics for the layered perovskites are the motif (Sr2+) that separates the layers, and the offsetting of the layers from each other. As n increases from n = 1 to n = ∞, the dimensionality of the compounds changes from two to three.
The general formula of the Aurivillius phase is (Bi2O2) (A n+1BnO3n+1) . This structure consists of layers of Bi2O2 separated by n layers of perovskites.
Ideal perovskite ABO3 has symmetric cubic structure with a lattice constant of ∼4 Å. However, most perovskites, deviate from the cubic structure if the Goldschmidt tolerance factor , given by (, , and represents the ionic radii of ions A, B, and O, respectively), deviates from 1 . When , (i.e., radius of cation A is small), the O anions move toward the A cation, so BO6 octahedra tilt to shrink the available volume for A cations that is empty. Because B-O-B bonds are highly flexible and BO6 octahedra are flexible in shape and size depending upon the cationic size, valence, and position, the overall deformation reduces the cubic symmetry, and result in structural transitions to orthorhombic, tetragonal, or hexagonal states that have lower symmetry than the cubic state.
3.2. d orbitals
The critical factor in the orbital physics involving d orbitals is their anisotropic charge distributions that arise from the wave functions, which take different shapes depending on energy when electrons are bound to atomic nuclei by Coulomb force. Electronic properties of perovskites would be severely affected by the chemistry of the transition metals at the center of BO6 octahedra that have corner-sharing oxygen anions [15, 16]. Consider a transition metal element which is surrounded by six O2- in the BO6 octahedron. This configuration gives rise to a crystal field which hinders the free motion of d-electrons; consequently, orbital angular momentum is usually quenched, and the d orbitals due to the crystal field energy ∆oct split in energy into (1) dx2-y2 and dz2 states, which form twofold degenerate eg orbitals with higher energy, and (2) threefold degenerate t2g orbitals dxy, dyz, and dzx states at lower energy (Figure 6). As the degeneracy of spherical symmetry of an isolated atom is removed, the d orbitals begin to fill, starting from a low-energy state and continuing to higher energy states. The actual filling arrangement depends on the competition between the crystal field and on-site exchange interaction described by Hund’s rule . For example, for 5d perovskite SrIrO3, Ir4+ has five electrons in the d orbitals (5d5); the electrons distributions from basic viewpoint are illustrated (Figure 7). For transport properties, the band structure would depend sensitively on an overlap between the d orbital of the B-site transition metal element and the p orbitals of the surrounding oxygen’s.
Generally, in transition metal oxides, the electronic properties are further complicated by the interaction of various degrees of freedom surrounding the d-electrons, i.e., charge, orbital, spin, and lattice. These degrees of freedom would produce relevant energy scales of similar magnitude such as bandwidth, Coulomb repulsion, and SOC. Because of the large number of transition metals, the number of possible metal-based perovskite transition is enormous. These transition metals include those of which the order orbital is 3d (e.g., Fe, Co, Ni), 4d (e.g., Mo, Ru, Rh), or 5d (e.g., W, Re, Ir). 3d orbitals are well localized and thus form a narrow band (W) with a large on-site Coulomb interaction (U). 4d orbitals are spatially more extended than their 3d counterparts. As 5d orbitals are more spatially extended than 3d or 4d ones, as a result, nearest-neighbor orbitals overlap significantly, and therefore W is wider in 5d orbitals than in 3d or 4d cases, i.e., W3d < W4d < W5d.
Because 5d orbitals are extended, on-site Coulomb repulsion or correlation U is weaker for 5d orbitals than form 3d or 4d ones, i.e., U5d < U4d < U3d. Thus, naively one would expect 5d systems to be more metallic and less magnetic than those based on 3d and 4d transition metal oxides. Indeed, 5d perovskite SrIrO3 is a correlated paramagnetic metal according to the expectation, but this is an exception. Surprisingly, many other 5d perovskites such as Sr2IrO4, and Sr3Ir2O7 are insulators. This unexpected fact can be explained by the high SOC in 5d systems. In general, 5d systems have larger SOC than do 3d or 4d systems, i.e., Λ3d < Λ4d < Λ5d because SOC is proportional to the fourth power of the atomic number Z (i.e., ) . In iridate compounds, Z = 77 for Ir, so the SOC strength is very high, even comparable to on-site Coulomb repulsion (∼0.5 eV). This high SOC strength leads to modification of the electronic structure and gives rise to novel emergent phenomena.
3.3. Ruddlesden–Popper series Srn+1IrnO3n+1
3.3.1. Spin-orbit coupling and band structure evolution
In this section, we present some unusual phenomena that are observed in RP series Srn+1IrnO3n+1 (n = 1, 2, and ∞). Because the SOC strength is so much larger in Ir than in a typical 3d system, SOC contributes to lifting the fivefold degeneracy of the atomic d-levels. In crystals, the crystal field and SOC act together to split the t2g levels into Jeff = 3/2 and Jeff = 1/2 levels (Figure 8). Two equivalent views of the splitting of 5d orbitals due to ∆oct and Λsoc are illustrated .
Now, considering electron hopping in solids, the Jeff levels would become bands. Ir4+ has five 5d electrons; four of them fill the lower Jeff = 3/2 band and one partially fills the Jeff = 1/2 band with the Fermi level in the Jeff = 1/2 band. The band structure evolves across the RP series Srn+1IrnO3n+1 (Figure 9). In the n = 1 case, Sr2IrO4, SOC splits the t2g band into two Jeff bands, so the bandwidth of the conduction band (Jeff = 1/2) is effectively reduced and correlation U can split the conduction band again into upper Hubbard bands (UHB) and lower Hubbard bands (LHB). Thus, Sr2IrO4 becomes a Mott insulator (Figure 9(a)). As n increases in the RP series, the bandwidth of the Jeff = 1/2 band also increases, possibly due to increase in the coordination number. Sr3Ir2O7, a barely insulator, is the intermediate case (Figure 9(b)). In SrIrO3 with n = ∞, U cannot split the conduction band with relatively wide bandwidth, so this compound is a correlated metal (Figure 9(c)).
The first material in the RP series of Srn+1IrnO3n+1 is a layered compound Sr2IrO4 (n = 1). The crystal structure of Sr2IrO4 consists of alternating layered stacking of SrO-IrO2-SrO perovskite units with ideal K2NiF4-type tetragonal cell (Figure 5). The lattice constants are ∼5.494 Å in the ab-plane and ∼25.796 Å along the c-axis with space group I41/acd . Most importantly, it has a highly insulating nature (Figure 10(a))  in contrast to a naive expectation that the extended nature of 5d orbitals would lead to a significant overlap of the nearest-neighbor orbitals and thus a broad electronic bandwidth. This insulating state arises from strong SOC and subsequent Coulomb repulsion (Figure 9(a)) and is termed a “Mott Jeff = 1/2 insulating state” of the 5d electron system [53, 54]. Magnetically, Sr2IrO4 exhibits weak ferromagnetism (canted antiferromagnetism) with TN ~ 240 K (Figure 10(b)) with very small ferromagnetic moment (0.023 μB/Ir) . The Mott insulating state due to strong SOC was first confirmed by Kim et al. by x-ray absorption (Figure 10(c)), based on the selection rules associated with the 2p to 5d transitions. Kim et al. also used theoretical calculations to identify the unusual nature of the Ir4+ state in accord with . Angle Resolved Photoemission Spectroscopy (ARPES) (Figure 10(d)) also confirmed the Mott Jeff = 1/2 state .
Considering the origin of the Mott Jeff = 1/2 insulating state, it should be very sensitive to external perturbations; if so, the relevant energy parameters can be tuned to modify the electronic states. In Sr2IrO4 thin films under tensile (compressive) strain, the correlation energy is affected by in-plane lattice strain with increase (decrease) in bandwidth . Sr2IrO4 remains an insulator even under pressure up to 55 GPa ; this stability illustrates the robustness of this insulating state. Also, a pressure-induced, fully reversible, giant piezoresistance was detected at room temperature . The electronic band gap could be tuned electrically, and this characteristic demonstrates potential application in next-generation electronic devices . Alkali-metal doping induces a close resemblance of the electronic state to that of high-temperature cuprates and therefore represents a step toward high-temperature superconductivity . Many other exciting properties have been observed in this compound, including magnetic structural change, spin-orbit tuned MITs, lattice-driven magnetoresistivity, electron-doped tuned electronic structure, anisotropic magnetoresistance, and excitonic quasi-particle [61–66]. The interplay of crystal field splitting, SOC, and correlation effects in layered Sr2IrO4 determines the 5d electronic structure and leads to realization of a completely new class of materials with a novel quantum state.
The crystal structure of Sr3Ir2O7 (n = 2) in this RP series of Srn+1IrnO3n+1, is of tetragonal cell with a = 3.896 Å and c = 20.879 Å and space group I4/mmm . It consists of strongly coupled bilayers of Ir-O octahedra which are separated by Sr-O interlayer (Figure 5). With the increase of number of octahedral layers (n), the electronic bands progressively broaden, and in particular, the bandwidth of the Jeff = 1/2 band increases from 0.48 eV for n = 1 to 0.56 eV for n = 2 . Still the transport measurement shows a well-defined, barely insulating Jeff = 1/2 states (Figure 11(a)) . Theoretical calculations based on LDA + U + SOC also provide evidence for the existence of the barely insulating band structure in which Fermi energy is between the Jeff = 1/2 bands (Figure 9(b)) . The onset of weak ferromagnetism occurs at TC ~ 285 K and is closely associated with the rotation of IrO6 octahedra about the c-axis. Indeed, the temperature dependence M(T) of magnetization closely tracks the rotation of the octahedra, as characterized by Ir-O-Ir bond angle. Sr3Ir2O7 also exhibits an intriguing M reversal for in-plane magnetization below 20 K with the onset of a rapid reduction at TD ~ 50 K (Figure 11(b)). The barely insulating nature related to splitting of the Jeff = 1/2 band was again observed with x-ray scattering and absorption (Figure 11(c))  and supported by ARPES measurements (Figure 11(d)) .
The ground state of bilayered Sr3Ir2O7 is highly sensitive to small external perturbations such as chemical doping, high pressures, and magnetic field. By replacing Sr2+ with La3+, electrons can be doped into bulk samples and can lead to an insulator-to-metal transition . Also, the application of a high hydrostatic pressure leads to a drastic reduction in the electrical resistivity; this observation suggests that the system is near an MIT . Although this system has not been fully explored yet, some studies such as resonant inelastic x-ray scattering, scanning tunneling spectroscopy, and optical conductivity have been performed [72–74]. These observations indicated that Sr3Ir2O7 is a good model system to explore the mechanism for novel Mott states near an MIT boundary at which competitive interplay between SOC and Coulomb interactions persists.
The last compound in the RP series (n = ∞) is SrIrO3; it is a correlated metal. SrIrO3 is a three-dimensional system and has the largest coordination number in the RP series. An increase in the coordination number would lead to the increase of the bandwidth, and thus correlation-driven band splitting would not occur (Figure 9(c)). In fact, in SrIrO3, the bandwidth of the Jeff = 1/2 band reaches 1.01 eV, so no gap appears within the Jeff = 1/2 band or between the Jeff = 1/2 band and Jeff = 3/2 band .
In 1971, Longo et al. first synthesized polycrystalline SrIrO3 . The stable ambient structure of SrIrO3 is a monoclinic distortion of hexagonal BaTiO3 structure [75, 76]. However, the structure transforms to an orthorhombic perovskite (space group Pnma) at 40 kbar and T = 1000 °C (Table 2). If a perovskite sample obtained at high temperature and high pressure is quenched, it remains a perovskite at room temperature. Because this book is about perovskites, our main focus remains on the properties of perovskite SrIrO3. Zhao et al.  and Blanchard et al.  again synthesized orthorhombic perovskite samples under high pressure and performed electric and magnetic measurements which showed that perovskite SrIrO3 is truly a paramagnetic metal (Figure 12).
|Synthesis conditions||Structure||Lattice parameters|
T = 900 °C
|a = 5.604 Å|
b = 9.618 Å
c = 14.17 Å
β = 93.26°
|PO2 = 40 kbar,|
T = 1000 °C
|Orthorhombic||a = 5.60 Å|
b = 5.58 Å
c = 7.89 Å
In situ ARPES showed that perovskite SrIrO3 is an exotic narrow-band semimetal . The bandwidth was surprisingly narrower than other two-dimensional RP phases, and the semimetallic nature is caused by the unusual coexistence of heavy hole-like and light electron-like bands contrary to the coordination number argument given previously. The observed unusual property may originate from the interplay of strong SOC, dimensionality, and both in- and out-of-plane IrO6 octahedral rotations. Recent theoretical calculations also suggest an extremely interesting possibility that the interplay of the lattice structure and large SOC produces Dirac nodes in the Jeff = 1/2 band, and engineering topological phases at interfaces and in superlattices would alter the system to be close to a topological crystalline metal [80–82]. It is obvious that SrIrO3 is an intriguing system in its own right and further studies are warranted.
4. Non-Fermi liquid physics and metal–insulator transitions in SrIrO3 films
Among RP series compounds Srn+1IrnO3n+1, the end member SrIrO3 is of particular interest because it is a correlated metal that exhibits unusual electronic transport properties that deviate from the normal Fermi liquid behaviors. The MITs and associated non-Fermi liquid physics that occur in SrIrO3 with strong SOC would provide an opportunity to extend the limit of our current knowledge of the physics of MITs as presented in the previous section. Thus, the present section constitutes the core of this chapter, and we first summarize salient features of the transport properties of SrIrO3, particularly in thin film form. Then we proceed to present the relevant theoretical frameworks to understand the non-Fermi liquid physics that underlay the MITs in the system.
4.1. Salient features of the metal–insulator transitions in SrIrO3
SrIrO3 is believed to be close to an MIT as evidenced by the evolution of the RP series Srn+1IrnO3n+1 from being an insulator to a correlated metal with increasing n = 1 2 ∞. The transport properties of correlated SrIrO3 can be anticipated to be susceptible to external perturbations and MITs, and that the associated unusual properties could be induced in SrIrO3 if, for example, an external stress is applied to the system. As introduced in the earlier section, the two most important mechanisms for MITs in correlated transition metal oxides are correlation-driven Mott localization and disorder-driven Anderson localization. When the system is under variable external stress, one may be able to tune W, U, and D to some extent and expose the interesting physics that is controlled by the parameters known as effective correlation (U/W) and effective disorder (D/W). For these reasons, we attempted to synthesize perovskite SrIrO3. Perovskite SrIrO3, however, is metastable at room temperature and is obtainable only by applying an elevated pressure (~40 kbar) at high temperature (~1,000 °C) and subsequent quenching. While it is not easy to obtain single crystals of perovskite SrIrO3 due to technical difficulties dealing with high pressures, the crystals can be stabilized by using thin film synthesis to produce them. In this case, the underlying substrates provide compressive strain, which replaces pressure, and epitaxial perovskite SrIrO3 thin films are easily obtained. More importantly, compressive strains can be imposed on SrIrO3 films by choosing substrates with appropriate lattice parameters. By depositing SrIrO3 films on GdScO3, DyScO3, SrTiO3, and NdGaO3, one can impose progressively larger compressive strain in the films (Figure 13). Biswas et al.  provides further details on thin-film synthesis.
The various scenarios on MITs suggest that an MIT may be obtained by varying the thickness of SrIrO3 films. Indeed in films deposited on GdScO3, which has a lattice constant well matched with that of SrIrO3 (Figure 13), an MIT occurs as the thickness is reduced from 4 nm to 3 nm [83, 84]; a 4-nm film is metallic (with a resistivity upturn at low temperatures), but a 3-nm film is insulating. The resistivity upturn at low temperatures in the 4-nm film is well described by the weak localization in two dimensions, and the insulating behavior of the 3-nm film can be explained by the variable-range hopping. Thus, the thickness-driven MIT for SrIrO3 (Figure 14(a)) is due to disorder and falls in the class of Anderson localization . Reducing the thickness of the film seems to increase the effectiveness of given disorder and of the grain size effect to scatter the charge carriers. The transition from being metallic with a low temperature upturn in resistivity for small disorder to fully insulating over the whole temperature range with the increase of disorder is a realization of the disorder-driven MIT. However, the temperature variation of the resistance of all the metallic films (thickness 4 nm, 10 nm, and 35 nm) follows ρ ∝ T4/5. This nontrivial exponent indicates that although the MIT itself is driven by disorder, the underlying transport mechanism not so simple. Also, in some materials (e.g., SrVO3), thickness-dependent MIT is caused by a decrease of the coordination number or by an increase in the correlation effect .
MIT also occurs in perovskite SrIrO3 thin films when compressive strain is imposed on the films when the thickness is kept constant (Figure 14(b)) . The imposed strain changes the Ir-O bond length d and the Ir-O-Ir bond angle θ; these changes affect the bandwidth as , where ψ = (π – θ)/2 is the buckling deviation of the Ir-O-Ir bond angle θ from π . d is relatively difficult to change, but θ can be readily affected by lattice strain. When θ is decreased by compressive strain, the electronic hopping integral between Ir 5d orbitals is reduced, so the buckling drives the system toward the insulating state. In addition, compression would also slightly increase correlation (U). Thus, the overall change in effective correlation (U/W) would seem to induce the MIT; however, as compressive strain increases, the effective disorder (D/W) in the system also increase. Disorder was already shown to play a decisive role in the thickness-driven MIT and has a strong influence in the strain-driven MIT, as exposed by optical conductivity measurements . Optical absorption spectroscopy provides insight into electronic band structure and free-carrier dynamics; optical absorption spectroscopy measurements of compressively strained SrIrO3 films showed Drude-like, metallic responses without an optical gap opening. This result indicates that localization has a measurable effect on strain-induced MIT in perovskite SrIrO3 thin films. In fact, extensive transport measurements in the compressively strained SrIrO3 revealed is that the MIT is not simply due to either disorder or correlation [83, 87]. Thus, SrIrO3 seems to provide a rare example in which the interplay of correlation and disorder in the presence of SOC causes the MIT.
The most remarkable feature of the electrical transport of perovskite SrIrO3 under compression is that the temperature variation of resistivity deviates from the Fermi liquid behavior in a peculiar way. For SrIrO3 films of 35-nm thickness on GdScO3, DyScO3, SrTiO3, and NdGaO3 substrates, the electrical resistivity not only shows non-Fermi liquid behaviors (ρ ∝ Tε with ε ≠2), but ε evolves from 4/5 to 1 to 3/2 as the compressive strain is increased, specifically, ρ ∝ T4/5 for films on GdScO3, ρ ∝ T for films on DyScO3, and ρ ∝ T3/2 for films on SrTiO3 (Figure 15) . Films on NdGaO3 are subjected to the largest strain and become insulating (Figure 14(b)). The present strain-driven MIT is clearly contrasted to the thickness-driven MIT, during which ε remains constant at 4/5 as the thickness is reduced and the system approaches the MIT .
The cause of these peculiar non-Fermi liquid behaviors and consequent MIT in the strain-driven case is not clear. Arguably, because SrIrO3 is paramagnetic without long-range magnetic ordering, localized states might induce the formation of local magnetic moments without collective magnetic fluctuations because they occur near the MIT. Such localized moments or small magnetic clusters can influence the electronic transport significantly in the presence of disorder. Indeed, the presence of disorder in the strained films is indicated by the increases in resistivity at low temperatures (Figure 14). Also, disorder and the evolution of non-Fermi liquid physics are possibly inter-connected. In a correlated metal, when disorder is sufficiently high, the system can enter a so-called Griffiths phase, which consists of a mixture of islands of Fermi liquid and Mott insulating regions, and that has non-Fermi liquid behaviors. In the following sections, we will try to propose a model that incorporates all these factors and that may explore a new paradigm for non-Fermi liquid physics.
4.2. UV physics: Emergence of localized magnetic moments and dynamical mean-field theory
Essential experimental features for SrIrO3 films of thickness 35-nm grown on GdScO3, DyScO3, SrTiO3, and NdGaO3 are non-Fermi liquid transport phenomena near MIT, where electrical resistivity shows anomalous temperature dependences with ε = 4/5, 1, 3/2, and –1/4 respectively. (The resistivity for the insulating phase can be described with a negative exponent.) The continuous change of the temperature exponent implies that a particular type of interplay between correlations of electrons and disorders is expected to have an important influence on physics near the MIT and raises fundamental questions about the nature of the non-Fermi liquid. The first question is whether this non-Fermi liquid physics is involved with either UV (ultraviolet) or IR (infrared) physics. Here, UV physics means that localized magnetic moments appear near an MIT and are regarded to be the source of strong inelastic scattering events due to their extensive entropy and to be responsible for non-Fermi liquid transport phenomena. The theoretical framework of dynamical mean-field theory was designed to simulate this local-moment physics quite well .
Dynamical mean-field theory describes Mott quantum criticality successfully and suggests that it appears at high temperatures, where not only MITs but also the so-called bad metal physics have been revealed [89, 90]. A noticeable point is that this local-moment UV physics seems to be universal, regardless of IR physics in which such local moments are expected to disappear by forming either singlets or magnetic orders and thereby reducing the entropy dramatically at low temperatures. Because incompletely screened local moments have important effects in non-Fermi liquid physics, the appearance of negative magnetoresistance (MR) (MR ) can be expected. However, the experiment on 35-nm SrIrO3 films on various substrates confirms positive MR, which is less than 1% up to magnetic field strength of 9 Tesla (Figure 16) . Although this result does not necessarily mean that the local-moment physics may not be important in the MIT of SrIrO3 thin films, the positive MR leads us to focus on IR physics.
Another important phenomenon is that spectral weight transfer in optical conductivity from the Drude part to a mid-infrared region as temperature increases . This change is a characteristic feature of bad metals and occurs at finite temperatures in the metallic side near an MIT. Dynamical mean-field theory reveals that local-moment UV physics describes this bad-metal Mott physics nicely. If local-moment physics are not considered, the spectral weight transfer in the optical conductivity cannot be easily reproduced without symmetry breaking. Unfortunately, no experiment has been performed yet in determining the importance of UV local-moment physics on non-Fermi liquid states of SrIrO3 thin films. With the contradictory picture of UV physics and observed positive MR, in the next section, we focus on IR physics, which assumed to be responsible for the observed non-Fermi liquid physics near the MIT of SrIrO3 thin films.
4.3. IR physics: Slater quantum criticality vs. Mott–Anderson–Griffiths scenario
Here, IR physics means that long-wave length and low-energy fluctuations determine the non-Fermi liquid physics near the MIT of SrIrO3 thin films. Then, quantum criticality would be the first choice, where quantum critical fluctuations involved with symmetry breaking, which cause to strong inelastic scattering between low-energy electrons. This scattering is responsible for non-Fermi liquid transport phenomena . Unfortunately, we failed to figure out quantum criticality involved with any kinds of orders, in particular, those associated with magnetism, although we do not exclude more-delicate symmetry breaking near the MIT of SrIrO3 thin films. Furthermore, quantum criticality itself cannot readily explain the continuous change of ε in the relationship of electrical resistivity to temperature on the metallic side of the MIT. One way to understand this continuous change is to consider the effect of disorder on quantum criticality. The Harris criterion is on the stability of quantum criticality against weak randomness, in the sense of average, involved with a space dimension and a critical exponent of correlation length . When the Harris criterion is violated, the clean quantum critical point becomes destabilized. As a result, a novel disordered quantum critical point can emerge, respecting the Harris criterion. However, the continuous evolution of the non-Fermi liquid physics is difficult to understand even in this situation. However, the strength of randomness can grow indefinitely, resulting in the so-called infinite randomness fixed point; when this happens, samples become extremely inhomogeneous due to effectively enhanced disorders, with ordered regions coexisting with disordered islands. The statistical distribution of the ordered islands shows a power-law tail, which implies that rare events that correspond to the power-law tail of the distribution function have important influences in non-Fermi liquid physics, particularly thermodynamics . For example, an antiferromagnetic Heisenberg model with random exchange coupling lies at an infinite randomness fixed point, at which the antiferromagnetic quantum critical point disappears and is replaced by a random singlet state . An attractive feature of the infinite randomness fixed point is that non-Fermi liquid physics still survive near it, i.e., away from the quantum critical point, where dynamics of rare regions associated with the long tail part of the distribution function governs singular behaviors of the extremely inhomogeneous state. These phenomena are referred to as quantum Griffith effects . The continuous change of non-Fermi liquid transport exponents reminds us of the quantum Griffith phase.
We suspect that the MIT in 35-nm SrIrO3 films on GdScO3, DyScO3, SrTiO3, and NdGaO3 can be achieved by electron correlations, in which lattice mismatches between substrates and SrIrO3 thin films are expected to control the ratio between interactions and hopping integrals. An essential question is on the role of disorder in this Mott transition, combined with electron correlations. Considering that the films show positive MR at magnetic field strength up to 9 Tesla, we suggest that magnetic correlations may not have important influence in this MIT. In fact, the positive MR suggests that the resulting insulating phase may be paramagnetic.
Furthermore, the residual resistivity is close to the Mott–Ioffe–Regel (MIR) limit; this similarity implies that the concentration of disorders is not low. Therefore, we conjecture that the interplay between electron correlations and disorders can give rise to a Griffith-type phase between Landau’s Fermi liquid state and the Mott–Anderson insulating phase; the Griffith-type phase allows the continuous change of transport exponents. We call this physics the Mott–Anderson–Griffith scenario. Until now, the Griffith scenario has been realized near the infinite randomness fixed point , at which extreme inhomogeneity and associated rare events are responsible for non-Fermi liquid physics that have varying critical exponents. Although the mechanism by which such an infinite randomness fixed point can appear in the Mott–Anderson transition has not been identified, fluctuations between metallic and insulating islands as rare events are expected to allow development of the Mott–Anderson–Griffith phase.
4.4. Model Hamiltonian and tentative Global phase diagram
Recalling the interplay between the spin-orbit coupling and the Hubbard interaction, we start from an Anderson–Hubbard model with one band (Eq. (2)):
where represents an electron field at site , and expresses a Kramers doublet state given by total angular momentum, is a hopping integral, is the strength of on-site Coulomb interaction, and is a random potential introduced by disorder. When electron correlation becomes negligible, the Hamiltonian reduces to the Anderson model, showing a continuous phase transition from a diffusive Fermi liquid state to an Anderson insulating phase in three dimensions. Local density of states (LDOS) occurs (Figure 17), in which localized eigenstates are given by discrete energy levels. This condition occurs below the mobility edge in the diffusive Fermi liquid state, and over the whole range of energy in the insulating phase. When contributions of random impurities can be neglected, the Hamiltonian becomes the Hubbard model, which shows a Mott transition from Landau’s Fermi liquid state to a paramagnetic Mott insulating phase. Regarding this Mott transition based on the UV local-moment physics yields first-order MIT . In contrast, using IR physics to explain the Mott transition yields a continuous MIT, which will be discussed below. In the Landau’s Fermi liquid state, a coherent peak occurs at the Fermi energy in the electron spectral function besides incoherent humps, and it disappears to transfer into an incoherent background, resulting in upper and lower Hubbard bands in the paramagnetic Mott insulating phase. A further question is whether the Griffith-type phase appears near the Mott–Anderson transition in the middle region of the phase diagram (Figure 17).
As discussed, we try to describe IR physics to explain the Mott–Anderson MIT in SrIrO3 thin films. Incompletely screened local moments can be expected to be screened completely at low temperatures. One mechanism for this screening is the Kondo effect, in which localized magnetic moments form singlets with itinerant electrons, well described by dynamical mean-field theory. The other mechanism entails localized magnetic moments form singlets with themselves as an RKKY-type (Ruderman, Kittel, Kasuya, Yosida) interactions [97–99]. The resulting paramagnetic Mott insulating phase is called a spin-liquid state, in which charge-neutral spinons may emerge as low-energy elementary excitations and may interact among themselves through abundant singlet fluctuations, referred to as gauge fields . Recalling the positive MR, we discuss the Mott–Anderson transition based on the spin-liquid physics.
4.5. Generalization of Finkelstein’s nonlinear sigma model approach near the Mott–Anderson transition
To describe the Mott transition involved with the spin-liquid physics, we take the U(1) slave-rotor representation  as , in which an electron field is expressed as a composite operator of charge and spin degrees of freedom. accounts for the dynamics of collective charge fluctuations (sound modes), and expresses a charge-neutral spinon field for collective dynamics of spin degrees of freedom. Then, an effective mean-field theory for such variables of and can be easily formulated from the Hubbard model without randomness; in this model, the Hubbard Hamiltonian is decomposed into two sectors describing the dynamics of spinons and zero-sound modes, respectively, given by (Eq. (3) and Eq. (4)):
Here, the conventional saddle-point approximation has been performed for a spin-liquid-type Mott insulating phase. describes band renormalization for electrons, and approximately expresses the width of incoherent bands. is a Lagrange multiplier field to control the spin-liquid to Fermi liquid phase transition and is regarded to be the chemical potential of bosons. This results from a nonlinear model description, in which the rotor variable is replaced with and the unimodular constraint is considered. is the nearest coordinate number of our lattice, and is the size of the system.
The Mott transition in this spin-liquid approach is realized by the condensation transition of , by which zero-sound modes are gapless in the Fermi liquid state, i.e., but become gapped in the spin-liquid phase, i.e., . As a result, the height of the coherent peak in the electron spectral function, proportional to the condensation amplitude , decreases gradually to disappear toward the spin-liquid Mott insulating phase; during this process, the spectral weight is transferred to the incoherent background of upper and lower Hubbard bands (Figure 18). Our problem is to introduce a random potential into this spin-liquid Mott transition. As long as the strength of the random potential remains smaller than the spinon bandwidth, we are allowed to deal with the role of the random potential perturbatively. This renormalization group analysis has been performed to reveal that the clean spin-liquid Mott critical point becomes unstable as soon as the random potential is turned on. As a result, a disorder critical point appears to be identified with a transition from diffusive spin-liquid glass insulator to diffusive Fermi liquid metal, in which the diffusive spin-liquid glass state consists of a diffusive spin-liquid phase of spinons and a charge glass phase of sound modes (Figure 19) .
In this renormalization group analysis, the spinon conductivity has been used as an input parameter, rather than being self-consistently determined, and the diffusive dynamics of spinons has been assumed. As a result, the renormalization group study reached the conclusion that the weak-disorder quantum critical point depends on the residual spinon conductivity, which means that the universality may not appear around this disorder quantum critical point. Suppose that SrIrO3 thin films on the substrate of GdScO3 are near this disordered MIT. Since SrIrO3 thin films on DyScO3, SrTiO3, and NdGaO3 substrates are also near this quantum critical point and their spinon conductivities differ from each other, we may observe continuous change of the temperature exponent of the electrical resistivity. Unfortunately, however, this previous study does not evaluate the spinon conductivity self-consistently, so the discussion cannot be beyond our speculation. Therefore, a theoretical framework must be developed to determine both transport coefficients of spinons and sound modes self-consistently.
One problem is that the previous renormalization group analysis does not consider effective interactions between diffusions and Cooperons, i.e., weak-localization corrections . The replica nonlinear model approach serves a natural theoretical framework to introduce such quantum corrections [104, 105]. By introducing the replica trick into the effective action for spinons and sound modes, and taking the diffusive spin-liquid fixed point as a saddle point for the spinon dynamics, one can derive an effective field theory as follows (Eq. (5)) :
Dynamics of diffusions and Cooperons of the spinon sector are described by the replica nonlinear model of the first line. is a matrix field, located on the Grassmannian manifold of , which spans retarded and advanced, Kramers doublet, and N replica spaces. is a constant diagonal matrix, given by for the retarded part and for the advanced sector. is the spinon density of states at the spinon Fermi-energy and is the diffusion coefficient. a gauge field to describe singlet excitations, more precisely, spin-chirality fluctuations that are neglected in the previous mean-field analysis but introduced here beyond the saddle-point approximation , in which the dynamics of gauge fluctuations is described by the last term, which is given by the polarization function of the diffusive dynamics of spinons. When the effects of gauge fluctuations in the dynamics of diffusions and Cooperons can be neglected, the replica nonlinear model reduces to that for the Anderson localization [104, 105]. The boson sector of sound modes is essentially the same as before, but gauge fluctuations are introduced beyond the mean-field analysis and the time sector is modified by the renormalization of the diffusive dynamics of spinons. Nonlocal derivative terms for both space and time should be defined in the momentum and frequency space, where the above expression is just for our formal writing. As discussed before, the boson sector describes the spin-liquid Mott transition, whereas the nonlinear model part describes Anderson localization.
We claim that this replica nonlinear model approach [107–110] generalizes the Finkelstein’s nonlinear model approach into the region of Mott transitions. Finkelstein’s approach does not incorporate Mott physics associated with strong correlations of electrons. Instead, this approach considers the influences of effective interactions of both singlet and triplet channels in the dynamics of diffusions and Cooperons, considers instabilities of interacting diffusive Fermi liquids, and reveals the nature of the Anderson MIT. In contrast, the existing nonlinear model approach considers instabilities of interacting diffusive spin-liquids and the nature of the Mott–Anderson MIT. In diffusive Fermi liquids, effective interactions of the triplet channel have an important influence on the Anderson MIT of two spatial dimensions [107–110]. In our problem, we do not consider triplet-channel interactions because the positive magnetoelectrical resistivity suggests that the magnetic correlations may not be important. However, we speculate that gauge fluctuations, regarded to be singlet-channel interactions, affect the dynamics of Cooperons seriously, thereby suppressing weak-localization corrections and stabilizing the metallic state of spinons. This would be a novel mechanism to produce metallicity near the Mott–Anderson transition. Renormalization group analysis has not been performed yet for this effective field theory.
Transition metal perovskites have been the platform for numerous emergent physics that originate from the coupling of the fundamental degrees of freedom such as spin, lattice, charge, and orbital as well as of disorder, which is unavoidable in solids. The physics of 5d perovskite iridates, in particular, has attracted considerable attention after the discovery of the novel Jeff = 1/2 Mott insulating state and the evolution of dimensionality controlled MIT in the RP series Srn+1IrnO3n+1, originating due to strong SOC of 5d element which is comparable to its Coulomb correlation or bandwidth. SrIrO3, the end member of the RP series, shows many unusual phenomena, including different kinds of MITs. We fabricated epitaxial thin films of various thickness on various substrates to induce MIT in perovskite SrIrO3 by thickness reduction or imposed compressed strain. The MIT driven by thickness reduction occurs due to disorder, but the MIT driven by compressive strain in the films on different substrates is accompanied by peculiar non-Fermi liquid behaviors with an evolving temperature exponent in the electrical resistivity relationship. The latter MIT and associated non-Fermi liquid behaviors are probably due to the delicate interplay between correlation, SOC, and disorder, and thus pose a theoretical challenge to our understanding of non-Fermi liquid physics and MIT.
To reveal the nature of this non-Fermi liquid physics near the MIT of SrIrO3 thin films on various substrates, we first discussed the influence of emergent localized magnetic moments, referred to as UV physics. The observed positive MR in the whole temperature range led us to pursue another direction, referred to as IR physics, in which long-wave length and low-energy fluctuations would be important. Because quantum criticality itself cannot explain the continuous evolution of the non-Fermi liquid physics, we speculate that the interplay between strong correlation of electrons and not-so-weak disorder may contribute to the continuously varying non-Fermi liquid physics. One possible scenario within IR physics is to consider quantum Griffiths effects with extreme inhomogeneity, in which local fluctuations between metallic and insulating phases, referred to as rare events, and may dominate the non-Fermi liquid physics. To realize this so-called Mott–Anderson–Griffiths scenario, we try to combine spin-liquid Mott physics with Anderson localization that describes a Mott transition from Landau’s Fermi liquid state to a spin-liquid Mott insulating phase. Theoretical understanding of this special class of materials is in its early stages, and many new emergent phenomena in iridates are yet to be explored theoretically or experimentally.
YHJ acknowledges the support by NRF (grant no. 2011-0009231). KSK was supported by NRF (grant no. 2012R1A1B3000550) and also by TJ Park Science Fellowship of the POSCO TJ Park Foundation. YHJ and KSK were jointly supported by CTM at POSTECH (grant nos. 2011-0030785 and 2011-0030786).