Orbital reconstruction in nonpolar tetravalent transition-metal oxide layers.

A promising route to tailoring the electronic properties of quantum materials and devices rests on the idea of orbital engineering in multilayered oxide heterostructures. Here we show that the interplay of interlayer charge imbalance and ligand distortions provides a knob for tuning the sequence of electronic levels even in intrinsically stacked oxides. We resolve in this regard the d-level structure of layered Sr2IrO4 by electron spin resonance. While canonical ligand-field theory predicts g||-factors less than 2 for positive tetragonal distortions as present in Sr2IrO4, the experiment indicates g|| is greater than 2. This implies that the iridium d levels are inverted with respect to their normal ordering. State-of-the-art electronic-structure calculations confirm the level switching in Sr2IrO4, whereas we find them in Ba2IrO4 to be instead normally ordered. Given the nonpolar character of the metal-oxygen layers, our findings highlight the tetravalent transition-metal 214 oxides as ideal platforms to explore d-orbital reconstruction in the context of oxide electronics.

Their unique diversity of transport and magnetic properties endows transition-metal (TM) oxides with a long-term potential for applications in microelectronics and electrical engineering. Nowadays the search for new or superior properties goes beyond known bulk phases and includes oxide interfaces and stacked superlattices12. As compared with the bulk material, at interfaces the modification of the nearby surroundings can significantly affect the valence electronic structure, in particular, the occupation of the d-shell levels12345. This is often referred to as orbital reconstruction345 and brings to the fore the most basic aspect in electronic-structure theory: how energy levels in quantum matter are formed and populated.

A variety of intrinsically stacked crystalline oxides is presently known. The high-temperature cuprate superconductors6, for example, fall in this category but also iridates of the type A2IrO4 (A=Sr2+, Ba2+) that closely resemble undoped cuprates, both structurally and magnetically7891011. Sr2IrO4 has a rather simple crystalline structure displaying stacked, quasi two-dimensional (2D) IrO2 and double SrO layers. We shall demonstrate that in this system the occupation of the valence d electronic levels differs from what is expected in textbook ligand-field theory due to electrostatics that involves both types of metal-oxygen sheets. In particular, we show that, as compared with the isostructural cuprate La2CuO4, a different distribution of ionic charges between the TM-O2 and A-O layers modifies the sequence of energy levels within the t2g and eg manifolds and consequently very fundamental physical properties such as the magnetic g factors, which determine the relation between the magnetic moment and quantum number of a magnetic particle. Our findings are of direct relevance to the field of stacked oxide heterostructures and provide a guideline on how low-symmetry crystal fields at d-metal sites can be altered and potentially engineered through the appropriate design of successive ionic layers.

To show this we first use electron spin resonance (ESR) measurements to untangle the 5d-shell electronic structure of crystalline Sr2IrO4, in particular, the exact order of the Ir t2g levels. The single s=1/2 hole present in these t2g orbitals carries an angular moment leff=1 and is subject to a large spin-orbit coupling (SOC), which in first approximation results in an effective Ir4+ moment, or pseudospin, jeff=leff−s≈1/2 (refs 7, 12, 13). We compare the experimental properties of these pseudospins with the ones we have calculated by ab initio quantum chemistry methods. This combined approach, explored here on a strongly spin-orbit coupled material for the first time, provides direct access to the spatial anisotropies of the g factors and further to the detailed microscopic superexchange interactions. The ESR measurements and theory are found to agree on a quantitative level and moreover undoubtedly show that the d-level ordering in Sr2IrO4 is inverted with respect to the normal ordering in the sister iridate Ba2IrO4 or the isostructural 214 cuprate superconductors. The good agreement between the ESR data and the outcome of the computational methodology we describe and employ here establishes the latter as a reliable tool for the investigation of nontrivial electronic structures and magnetic couplings.

Results

Pseudospins and effective Hamiltonian

Mott-Hubbard physics in d-metal compounds has been traditionally associated with first-series (3d) TM oxides. However, recently, one more ingredient entered the TM-oxide ‘Mottness' paradigm—large SOC's in 5d systems. SOC in 5d and to some extent 4d anisotropic oxides modifies the very nature of the correlation hole of an electron, by admixing the different t2g components1213, changes the conditions for localization7, the criteria of Mottness and further gives rise to new types of magnetic ground states and excitations911. While various measurements indicate that indeed spin-orbit-coupled jeff≈1/2 states form in A2IrO4 (refs 7, 9, 14), it has been also pointed out that off-diagonal SOC's may mix into the ground state (GS) wavefunction substantial amounts of character15161718. Such many-body interactions were shown to produce remarkable effects in X-ray absorption and X-ray magnetic circular dichroism (XMCD): the branching ratio between the L3 and L2 Ir 2p absorption edges reaches values as large as 4, nearly 50% higher than the 2.75 value for a ‘pure' jeff=1/2 system19. In addition, low-symmetry noncubic fields produce sizeable splittings of the 5d t2g levels, in some cases close to or even larger than ∼1/2 eV (refs 20, 21, 22), and therefore admix the jeff=1/2 and components1315. The structure of the spin-orbit GS depends on both the strength and sign of these splittings. Interestingly, the best fits of the X-ray absorption and XMCD data are achieved in Sr2IrO4 with a negative t2g tetragonal splitting19, although the oxygen octahedra in this material display a distinct positive tetragonal distortion—the IrO6 octahedra are substantially elongated23 (a negative tetragonal splitting should occur when the IrO6 octahedra are compressed1324, see Fig. 1). This is already a first indication of the level inversion that our ESR measurements and quantum chemistry calculations show to take place in Sr2IrO4.

Figure 1

TM t2g splittings for tetragonal distortions of the oxygen octahedron sans SOC.

(a) z axis compression of the octahedron corresponds to a tetragonal splitting δ<0, causes an orbital dublet to be lowest in energy and the g factors to be ordered as . (b) Elongation of the octahedron (δ>0) causes an orbital singlet to be lowest in energy and the g factors to be ordered as . Purple dashed lines indicate the conventional zero level used to define the sign of δ.

The interactions between a pair 〈ij〉 of nearest-neighbour (NN) 1/2 pseudospins in the presence of an external magnetic field h is given by the effective Hamiltonian

where , are pseudospin (jeff≈1/2) operators, J is the isotropic Heisenberg exchange, D=(0,0,D) defines the antisymmetric Dzyaloshinskii-Moriya (DM) coupling, is a symmetric traceless second-rank tensor describing the symmetric anisotropy and due to the staggered rotation of the IrO6 octahedra the tensor splits for each of the two sites into uniform and staggered components (see for example, refs 16, 25). This effective spin Hamiltonian is of direct relevance to the interpretation of the ESR data.

ESR measurements

For a single crystal of Sr2IrO4 we observe antiferromagnetic resonance (AFR) modes in the sub-THz frequency domain26 as displayed in Fig. 2. There are two modes if h||z: a gapless Goldstone mode v||1=0 and a gapped excitation

Figure 2

ESR data for Sr2IrO4.

(a) Out-of-plane magnetic-field geometries; the inset shows a representative AFR spectrum. (b) In-plane magnetic fields; the inset demonstrates the g-factor anisotropy as function of the tetragonal distortion parameter α (see text). Symbols denote experimental data points—solid lines are theoretical curves using equations (2, 3, 4) and the quantum chemically computed g factors g||=2.31, g⊥=1.76 (see Table 1); dashed lines are calculated assuming isotropic g factors.

where , and Γ couples the and components (along the c axis, perpendicular to the ab-plane IrO2 layers23) in the third term of equation (1) (see Methods for details). Experimental results are shown for v||2 in Fig. 2a. The data comprise at T≪TN=240 K a group of overlapping resonances (Fig. 2a, inset), possibly due to some distribution of internal fields in the sample. Though revealing some scatter, the AFR data follow approximately a parabolic dependence on h and, most importantly, they lie substantially above the curve corresponding to the free-electron Landé factor ge=2 (dashed line in Fig. 2a). The experimental dependence can be reasonably well modelled with g|| values of 2.3–2.45. The solid line in Fig. 2a is obtained by using g||=2.31, as derived from quantum chemistry calculations that will be discussed later on.

Sizeable deviations to values >2 of g|| is clear indication of the presence of low-symmetry, noncubic crystal fields. In the simplest approximation, that is, restricting ourselves to the Ir4+ manifold, the anisotropic g factors in axial noncubic environment can be expressed up to the sign as13 g||=gc=(2+2k)cos2α−2sin2α and where k is a covalency reduction factor, α=(1/2) arctan parameterizes the deviation from octahedral symmetry, λ is the SOC constant and δ the Ir t2g splitting. A plot for the dependence of the diagonal g factors on the distortion parameter α is shown in the inset to Fig. 2b. For simplicity, k=1 is for the moment assumed but smaller values of k do not bring qualitative changes. In cubic symmetry δ=0, αcub=35.26° and the g matrix is isotropic with . According to standard textbooks on ligand-field theory24, an elongation of the out-of-plane Ir-O bond induces a positive tetragonal splitting of the Ir t2g levels, with δ>0 and α>αcub, whereas a bond compression yields δ<0 and α<αcub (see Fig. 1). As in Sr2IrO4 the IrO6 octahedra are substantially elongated in the z direction23, the g factors are expected to correspond to the case of positive splitting α>αcub, see the area to the right of the crossing point shown in the inset to Fig. 2b. It follows from the plot that g||<2, which obviously contradicts our AFR data for h||z [Fig. 2a]. The value g||=2.31 used to draw the curve connecting the open circles in Fig. 2a in fact corresponds to α=32.05°<αcub (see the inset in Fig. 2b) and indicates that, despite the positive c-axis tetragonal distortion, a counterintuitive negative tetragonal splitting of the 5d t2g levels is present in Sr2IrO4.

It should be noted that there is no Goldstone mode for finite in-plane magnetic fields. The canting angle depends in this case both on the strength of the DM interaction and the applied field. The two modes are

and

where and m is the in-plane ferromagnetic component of the effective moments. To first order in magnetic field, m can be expressed as

which holds for weak fields where m remains small. The first term corresponds to the zero-field canting, arising from the DM interaction, while the second term shows how this canting evolves with increasing h⊥. Plots based on equations (3, 4) and the quantum chemically derived interaction parameters (see below) are displayed in Fig. 2b together with ESR data.

Quantum chemistry calculations of g factors

Results of ab initio quantum chemistry calculations for the g factors in Sr2IrO4 and in the structurally related material Ba2IrO4 are listed in Table 1. Our computational scheme follows the prescription of Bolvin27 and Vancoillie et al.28. It maps the matrix elements (MEs) of the ab initio Zeeman Hamiltonian onto the MEs of the effective pseudospin Hamiltonian , where μ, L and S are magnetic moment, angular-momentum and spin operators, respectively. The spin-orbit GS wave functions are computed either at the complete-active-space self-consistent-field (CASSCF) or multireference configuration-interaction (MRCI) level of theory29, as described in ref. 30 and using the MOLPRO quantum chemistry package31. All necessary angular-momentum MEs are calculated as well with MOLPRO (see Methods). In a first set of calculations, only the three t2g orbitals at a given Ir site and five electrons were considered in the active space. The self-consistent-field optimization was carried out for the corresponding state. We use here the more convenient notations associated to O symmetry, although the calculations were performed for the actual experimental geometry, with point-group symmetry lower than octahedral. Inclusion of SOC yields in this case a set of three Kramers doublets (KDs), see Table 1.

Table 1

g factors for Sr2IrO4 and Ba2IrO4.

States considered	CASSCF		MRCI
	g⊥	g||	g⊥	g||
Sr2IrO4 (δ=−155 meV):
(3 KDs)	1.67	2.25	1.60	2.35
(27 KDs)	1.81	2.27	1.76	2.31
Ba2IrO4 (δ=65 meV):
(3 KDs)	2.00	1.61	2.01	1.60
(27 KDs)	2.09	1.77	2.10	1.76

CASSCF=complete-active-space self-consistent-field; KDs=Kramers doublets; MRCI=multireference configuration-interaction. Results of many-body quantum chemistry calculations are shown. The left column displays the electron configurations entering the spin-orbit treatment. Only the high-spin sextet state is considered out of the manifold. Final results are indicated in bold face.

Subsequently we performed calculations with larger active spaces, including also the Ir eg orbitals. One plus four (2A2g, 2T1g, 2Eg and 2T2g) spin doublets, two spin quartets [ and ] and one spin sextet entered here the spin-orbit treatment. The orbitals were optimized for an average of all these terms.

The effect of enlarging the active space to include terms in the reference wavefunction is in the range of 10%, in line with earlier semi-empirical estimates for 4d5 and 5d5 systems15161718. Most importantly, the calculations yield a negative tetragonal splitting of the Ir t2g levels in Sr2IrO4, δ=−155 meV by MRCI, and positive t2g splitting in Ba2IrO4 (see Table 1 and Methods). Similar signs, negative in Sr2IrO4 and positive in Ba2IrO4, but much larger magnitudes (≈0.7 eV) are found for the computed Ir eg splittings (not shown in Table 1).

Taken together, the ESR and quantum chemistry results unequivocally point at an anomalous order of the split Ir 5d levels in Sr2IrO4, related to the important role of the extended crystalline surroundings in generating low-symmetry fields that compete with ‘local' distortions of the ligand cage. Similar effects were found by ab initio calculations on the 214 layered rhodate Sr2RhO4 (ref. 32) and the 227 pyrochlore iridates22. In contrast, in Ba2IrO4, the stretch of the apical Ir-O bonds is strong enough33 to overcome the longer-range electrostatics, turning the tetragonal t2g positive again, as discussed in more detail in the following. Consequently, the structure of the tensor in Ba2IrO4 is qualitatively different, with g⊥>2 and g||<2 (see Table 1), the ordering that one normally expects and encounters for elongated octahedra13.

Exchange couplings from quantum chemistry

To obtain ab initio quantum chemistry values for the inter-site effective magnetic couplings in Sr2IrO4 (see equation (1)), we carried out additional calculations on larger clusters that incorporate two 5d5 sites. The two-octahedra cluster has C2v symmetry (see Fig. 3), which implies a diagonal form for and in equation (1) (see Methods). By one-to-one correspondence between the MEs of the ab initio Hamiltonian

Figure 3

Planar IrO2 network in Sr2IrO4.

(a) Coordination of the Ir site. Dashed lines show the boundaries of the crystallographic unit cell within a given IrO2 layer. (b) The point-group symmetry of the [Ir2O11] block is C2v; associated symmetry elements are indicated in the figure.

and the MEs of the effective spin Hamiltonian (1) in the basis of the lowest four spin-orbit states defining the magnetic spectrum of two NN octahedra, we can derive in addition to the g factors the strengths of the Heisenberg and anisotropic intersite couplings. In equation (6), is the scalar-relativistic Born-Oppenheimer Hamiltonian, describes spin-orbit interactions30 and is the two-site Zeeman Hamiltonian.

Diagonalization of the spin Hamiltonian (1) provides the expected singlet and three (split) triplet components |t〉, |t〉 and . Owing to the DM interaction, and are admixtures of ‘pure' |0,0〉 and |1,0〉 spin functions. Our mapping procedure yields J≈48 meV, somewhat lower than J values of 55–60 meV derived from experiment934, and a ratio between the DM and Heisenberg couplings D/J=0.25, in agreement with estimates based on effective superexchange models123536 and large enough to explain the nearly rigid rotation of magnetic moments that is observed when the IrO6 octahedra revolve3738. Ab initio results for the NN anisotropic couplings , also relevant for a detailed understanding of the magnetic properties of Sr2IrO4, are shown as well in Table 2. In our convention the x axis is taken along the Ir-Ir link, that is, it coincides with the 〈110〉 crystallographic direction23, and z||c. We obtain Γ≈Γ, which then allows to recast the Heisenberg and symmetric anisotropic terms in (1) as , with Γ=Γ=−Γ/2. Equally interesting, for no rotation of the IrO6 octahedra and straight Ir-O-Ir bonds in Ba2IrO4, it is Γ and Γ which are approximately the same, providing a realization of the compass-Heisenberg model3639 since the DM coupling is by symmetry 0 in that case.

Table 2

Nearest-neighbour magnetic couplings in Sr2IrO4.

J	D	Γxx	Γyy	Γzz
47.8	±11.9	0.42	−0.84	0.42

Results of spin-orbit MRCI calculations on two-octahedra clusters are displayed (meV). Γ+Γ+Γ=0 since is traceless.

The two-site magnetic Hamiltonian (1) features in-plane symmetric-anisotropy couplings Γ and Γ, which were not considered in previous studies1226. In the presence of two-sublattice order, terms containing these couplings cancel each other in the mean-field energy but they are in general relevant for pseudospin fluctuations and excitations. Using spin-wave theory and effective parameters derived from the quantum chemistry calculations, we nicely reproduce the correct GS and character of the modes, as shown in Fig. 2. To reproduce the experimental zero-field gap, in particular, we used J, D and g-factor values as listed in Table 2 and a somewhat larger Γ parameter of 0.98 meV. To leading order, the dependence of and on h is linear, see equations (3, 4), and the slope is proportional to m. At low fields (≤1 T), m can be actually replaced with its field independent value26. Using the MRCI coupling constants, the first term in equation (5) then yields a moment m≈0.12μB, in good agreement with the experiment840.

Discussion

The exact d-level order is of fundamental importance in TM oxides, dictating for instance the symmetry of the quasiparticle states in photoemission103241 and the nature of the magnetic ordering4243. In 214 iridates specifically, it determines the various isotropic as well as anisotropic contributions to the magnetic exchange couplings12353639, the evolution of those magnetic interactions with strain44 and/or pressure19 and most likely the nature of the intriguing transition to a nonmagnetic phase in Sr2IrO4 under high pressure19. Having established that in Sr2IrO4 the d levels are inverted and that in the closely related Ba2IrO4 they are not raises the question what actually drives the inversion. To address this, we performed an additional set of calculations, in which we change the charges around the reference IrO6 octahedron. As a simple numerical experiment that preserves charge neutrality of the A2IrO4 system, we assigned the 4 NN iridium sites (in-plane, see Fig. 4) the charge QTM−2Δq and the 8 closest A-site cations (out of plane) the valence QA+Δq. In a fully ionic picture, QTM and QA are 4+and 2+, respectively. However, since in our calculations the NN TM and A sites are not modelled as just formal point charges (see Methods), the actual valence states depart from their formal values, with larger ‘deviations' for QTM. The way we introduce Δq in the computations is therefore by appropriately modifying the nuclear charge at the respective site. For variable Δq, this interpolates linearly between nearby surroundings corresponding to 5d 214 layered perovskites (with Δq=0 and TM4+, A2+ formal valence states) and their cuprate 214 equivalents (with Δq=1, TM2+/A3+ formal ionic charges and ‘normal' order of the TM t2g and eg levels45).

Figure 4

Effect of interlayer charge imbalance in A2IrO4 iridates.

(a) The nearby surroundings of TM sites in A2TMO4-layered perovskites. In test calculations one can assign the adjacent (in-plane) TM ions the formal charge QTM−2Δq, which is compensated by assigning the NN A sites the charge QA+Δq. (b) Tetragonal crystal-field energy splitting between t2g orbitals (δ) as a function of the charge redistribution Δq for Sr2IrO4 and Ba2IrO4.

As is illustrated in Fig. 4a, increasing Δq amounts to moving positive charge from the IrO2 plane to the adjacent A-O layers. The calculations show that upon moving charge in such a manner, the Ir t2g splitting δ increases, see Fig. 4b. In other words, this redistribution of charge counteracts the level inversion in Sr2IrO4 and further increases the already positive δ in Ba2IrO4. In Sr2IrO4 the cubic-like jeff=1/2 limit occurs for Δq=0.22. This effect can easily be understood: placing more positive charge out of the IrO2 plane stabilizes the out-of-plane t2g orbitals, corresponding to the (yz,zx) orbital doublet, and thus enhances δ. One can also do the opposite and drive Δq negative. In this case more positive charge piles up in the IrO2 plane, which one expects to lower the energy of the xy orbital singlet, thus enhancing the level inversion in Sr2IrO4. This is indeed what happens, see Fig. 4b. What is more, driving Δq negative even causes a level inversion in Ba2IrO4, when Δq−0.25. It is interesting to note that the slope of the δ versus Δq lines in Sr2IrO4 is much steeper than in Ba2IrO4, which is caused by the significantly smaller Ir–Ir distances in Sr2IrO4.

From these test calculations it is clear that low-symmetry crystal fields associated to neighbours beyond the first ligand coordination shell, in particular, the highly charged Ir4+ NNs, counteract the local tetragonal crystal field that is caused by the elongation of the IrO6 octahedra in both Sr2IrO4 and Ba2IrO4. In the case of Ba2IrO4 the local distortion is still strong enough to overcome these longer-range effects but in Sr2IrO4, with a slightly smaller tetragonal distortion, the longer-range electrostatics wins, causing the observed level inversion.

While the role of the high ionic charge of in-plane ions has been earlier invoked in the tetravalent Ru oxide compound Ca2RuO4 (ref. 43) and in mixed-valence manganites46, we here explicitly prove it by combined ESR measurements and many-body ab initio calculations on structurally and chemically simpler systems in which additional complications arising from further distortions4347 or the presence of multiple TM valence states46 are excluded. A reversed order of the Ir t2g levels in Sr2IrO4 has been also indirectly implied by fits of X-ray absorption19 and X-ray magnetic scattering3748 spectra. As a more direct and more sensitive experimental technique to such details of the valence electronic structure and with back up from truly ab initio many-body calculations, ESR now provides irrefutable evidence for such physics. The numerical ‘experiment' outlined in Fig. 4 further shows that at the heart of this effect is not the intersite exchange, as assumed in ref. 19, and not the t2g-eg orbital hybridization invoked in ref. 48, but basic interlayer electrostatics.

We have, in sum, provided an integrated picture on the d-level structure and magnetic anisotropies in Sr2IrO4, a prototype spin-orbit driven magnetic insulator. Both the single-site tensor and intersite effective exchange interactions are analysed in detail. To access the latter, we build on an earlier computational scheme for deriving intersite matrix elements in mixed-valence spin-orbit coupled systems49. While the ratio D/J of the antisymmetric Dzyaloshinskii-Moriya and isotropic Heisenberg couplings is remarkably large in Sr2IrO4 and concurs with an in-plane rotation pattern of the Ir magnetic moments that follows nearly rigidly the staggered rotation of the IrO6 octahedra3738, the most prominent symmetric anisotropic terms are according to the quantum chemistry data in-plane, perpendicular to the Ir-Ir links. The structure of the tensor, as measured by ESR and computed with first-principles electronic-structure methods, is such that and distinctly indicates a negative tetragonal-like splitting of the Ir t2g levels, in spite of sizable positive tetragonal distortions in Sr2IrO4. We further observe that a much stronger tetragonal distortion in Ba2IrO4 renders the tetragonal d-level splitting positive and g||Sr2IrO4. The d-level inversion in Sr2IrO4 and the surprisingly small splitting in Ba2IrO4 have to do with the way the positive ionic charge is distributed between adjacent Ir4+O2 and A2+O layers, having in contrast to the 214 cuprate superconductors, for example, more positive charge in the TM-O planes. This almost compensates the ‘local' tetragonal field arising from the z axis elongation of the IrO6 octahedra in Ba2IrO4 and overcompensates it in Sr2IrO4.

The subtle interplay between local distortions of the O ligand cage and additional uniaxial fields associated with the anisotropic extended surroundings opens new perspectives on strain44 and pressure19 experiments in square-lattice iridates, for example, in connection to the spin-flop transition earlier predicted in Sr2IrO4 (refs 12, 36). It also opens up the perspective of manipulating this way the d-level ordering in oxide heterostructures with highly charged, trivalent and tetravalent species. Compounds with tetravalent species within the TM-O2 layers, in particular, given the nonpolar character of the quasi 2D sheets, provide ideal playgrounds to explore the mechanism of d-level ordering pointed out here since that will not be hindered by ‘interface' charge redistribution and structural reconstruction occuring in polar heterostructures from polar discontinuities5051. A recent experimental realization of such mixed, tetravalent/divalent TM-oxide interfaces is for example the SrRuO3/NiO interface52, one system that requires in this respect closer theoretical examination.

Methods

Single-site magnetic properties

The g factors were obtained by computations on clusters which contain one central IrO6 octahedron, the four NN IrO6 octahedra and the nearby ten Sr/Ba ions. The solid-state surroundings were modelled as a large array of point charges fitted to reproduce the crystal Madelung field in the cluster region. To obtain a clear picture on crystal-field effects and spin-orbit interactions at the central Ir site, we cutoff the magnetic couplings with the adjacent Ir ions by replacing the tetravalent open-shell d5 NNs with tetravalent closed-shell Pt4+ species. This is a usual procedure in quantum chemistry studies on TM systems, see for example, refs 42, 45, 53, 54, 55, 56. We used energy-consistent relativistic pseudopotentials and valence basis sets of quadruple-zeta quality supplemented with f polarization functions for the central Ir ion57 and all-electron triple-zeta basis sets for the six adjacent ligands58. For the TM NNs, we applied energy-consistent relativistic pseudopotentials and triple-zeta basis functions57 along with minimal atomic-natural-orbital basis sets59 for the Os coordinating those TM sites but not shared with the central octahedron. The Sr and Ba species were modelled by divalent total-ion effective potentials supplemented with a single s function60. All O 2p and metal t2g electrons at the central octahedron were correlated in the MRCI calculations. The latter are performed with single and double substitutions with respect to the CASSCF reference (for technicalities, see refs 61, 62), which is referred to as MRCISD. To separate the metal 5d and O 2p valence orbitals into different groups, that is, central-octahedron and adjacent-octahedra orbitals, we used the Pipek-Mezey localization module63 available in MOLPRO. The computations with hypothetical (QTM−2Δq) and (QA+Δq) ionic charges at the TM and Sr/Ba sites next to the reference Ir ion were carried out as frozen-orbital multideterminant calculations (also referred to as CASCI) with three Ir t2g and five electrons in the active space and orbitals optimized for Δq=0.

The spin-orbit treatment was performed according to the procedure described in ref. 30. In a first step, the scalar relativistic Hamiltonian is used to calculate correlated wavefunctions for a finite number of low-lying states, either at the CASSCF level or at the MRCI level. In a second step, the spin-orbit part is added to the initial scalar relativistic Hamiltonian, matrix elements of the aforementioned states are evaluated for this extended Hamiltonian and the resulting matrix is finally diagonalized to yield SO wavefunctions.

The g factors were computed following the scheme proposed by Bolvin27 and Vancoillie et al.28 (for alternative formulations, see, for example, ref. 64). For the KD GS , the Abragam-Bleaney tensor1365, G=gg can be written in matrix form as

where

The MEs of are here provided by MOLPRO while those of are derived using the conventional expressions for the generalized Pauli matrices :

G is next diagonalized and the g factors are obtained as the possitive square roots of the three eigenvalues. The corresponding eigenvectors specify the rotation matrix to the main magnetic axes. In our case, the magnetic z axis is along the crystallographic c coordinate, while x and y are ‘degenerate' and can be any two perpendicular directions in the ab plane.

To cross-check the g-factor values computed with our subroutine, we further performed g-factor calculations using the module available within the ORCA quantum chemistry package66. We applied all-electron DKH (Douglas-Kroll-Hess) basis sets of triple-zeta quality for the TM ions67, triple-zeta basis functions for the ligands of the central octahedron58 and double-zeta basis functions for additional Os at the NN octahedra58. Dynamical correlation effects were accounted for by N-electron valence-state second-order perturbation theory (NEVPT2)6869. CASSCF and NEVPT2 results are listed in Table 3, for both Sr2IrO4 and Ba2IrO4. It is seen that the data in Tables 1 and 3 compare very well and indicate the same overall trends.

Table 3

Ir t 2g splittings and g factors for Sr2IrO4 and Ba2IrO4.

	δt2g (meV)	g⊥	g||
Sr2IrO4:
CASSCF	−127	1.66	2.23
NEVPT2	−199	1.55	2.41
Ba2IrO4:
CASSCF	30	1.93	1.74
NEVPT2	70	2.01	1.58

Results as obtained with the ORCA program66 are shown. Only the 2T2g() states were included in the CASSCF optimization and in the spin-orbit treatment.

Superexchange interactions in Sr2IrO4

NN magnetic coupling constants were obtained for Sr2IrO4 by calculations on an embedded cluster that includes two IrO6 octahedra as magnetically active units. As for the calculation of single-site magnetic properties, to accurately describe the charge distribution in the immediate neighbourhood, the adjacent IrO6 octahedra (six) and the closest Sr2+ ions (16) were also incorporated in the actual cluster. We used energy-consistent relativistic pseudopotentials along with quadruple-zeta basis sets for the valence shells of the two magnetically active Ir ions57, all-electron quintuple-zeta basis sets for the bridging ligand58 and triple-zeta basis functions for the other Os associated with the two reference octahedra58. We further employed polarization functions at the two central Ir sites and for the bridging anion, namely 2 Ir f and 4 O d functions5758. Additional ions defining the NN octahedra, the nearby Sr2+ species and the farther crystalline surroundings were modelled as in the single-site study, see above.

For two adjacent magnetic sites, the manifold entails nine singlet and nine triplet states. The CASSCF optimization was carried out for an average of these nine singlet and nine triplet eigenfunctions of the scalar relativistic Hamiltonian . In the subsequent MRCI treatment, only the Ir t2g and the O 2p electrons at the bridging ligand site were correlated. Results in good agreement with the experimental data were recently obtained with this computational approach for related 5d5 iridates2170.

Diagonalization of the Hamiltonian in the basis of the lowest nine singlet and nine triplet states provides a total of 36 spin-orbit-coupled eigenfunctions, namely, four , eight , eight and sixteen states. In the simplest picture, the lowest four roots imply either singlet or triplet coupling of the spin-orbit jeff=1/2 (or, more generally, pseudospin =1/2) on-site objects and are separated from higher-lying states by a gap of 0.5 eV, much larger than the strength of the intersite exchange. It is this set of lowest four spin-orbit MRCI roots that we map onto the eigenstates of the effective two-site (pseudo)spin Hamiltonian (1). The Zeeman interaction shows up on the quantum chemistry side as , where and are transformed to the spin-orbit-coupled basis using the spin-orbit wavefunctions as unitary transformation matrix. MEs of the ab initio model Hamiltonian are shown in Table 4. Diagonal components show the energies of the zero-field states, while the off-diagonal MEs describe the coupling to magnetic field.

Table 4

Matrix elements of the ab initio spin-orbit Hamiltonian.

		|tx〉	|ty〉
	0	0.2308iμBhy	−0.1768iμBhx	0
〈tx|	−0.2308iμBhy	48.3328	2.3083iμBhz	−1.6854iμBhy
〈ty|	0.1768iμBhx	−2.3083iμBhz	48.9626	1.6266iμBhx
	0	1.6854iμBhy	−1.6266iμBhx	49.0630

The latter is described by expression (6). Results of spin-orbit MRCI calculations are shown (meV). The two-site singlet and (split) triplet states are labelled and {|t〉, |t〉, }, respectively. Due to the antisymmetric exchange, and are admixtures of ‘pure' |0,0〉 and |1,0〉 spin functions.

For the experimentally determined crystal structure of Sr2IrO4 (ref. 23), the two-octahedra [Ir2O11] cluster displays C2 symmetry. Having the x axis along the 〈110〉 crystallographic direction23 and z||c, the effective anisotropic couplings read D=(0,0,D),

for Ir-Ir links along x and

for Ir-Ir links along y, with Γ+Γ+Γ=0. The uniform and staggered components of the tensor take for individual Ir sites the following form :

g=g|| while g and g are directly related to g⊥ but not restricted to be equal due to the lower symmetry of the two-octahedra cluster as compared with the IrO6 unit.

To solve now the actual problem, we need to transform the effective spin Hamiltonian (1) to the same form as the ab initio Hamiltonian shown in Table 4, that is, diagonal in zero magnetic field. The result of such a transformation is shown in Table 5. Direct correspondence between homologous MEs in the two arrays yields a set of eight independent equations that finally allow to derive hard values for all effective coupling constants that enter expression (1). While the results for the intersite exchange interactions in (1) and (10,11) are shown in Table 2, the -tensor data obtained from the two-octahedra calculations are g=1.64, g=1.70, g=2.31 and g=0.02. The way the additional parameters introduced in Table 5 are defined is explained in Table 6. The intermediate steps followed to arrive to the matrix form provided in Table 5 are outlined in Supplementary Tables I–III and the Supplementary Methods.

Table 5

Matrix elements of the effective spin Hamiltonian.

		|tx〉	|ty〉
	0			0
〈tx|			−igzzhz
〈ty|		igzzhz
	0

The explicit form of the latter is given by (1); for additional notations we use here, see Table 6.

Table 6

Notations used for anisotropic exchange coupling parameters in Sr2IrO4.

	J,Γxx=−Γyy−Γzz (present study)	Jx,Jy,Jz (ref. 36)	J,δJxy,δJz (x bond) (refs 123536)
Jx+y		(Jx+Jy)/2
Γx−y	(Γxx−Γyy)/2	(Jx−Jy)/2
Γz		Jz−Jx+y	δJz

The definitions on the left-hand side are applied in Table 5. Other conventions presently employed in the literature are also shown for comparison.

Spin-wave calculations

For understanding all details of the ESR spectrum, we carried out a spin-wave analysis using the Hamiltonian (1) and a site-factorized variational GS wave function

where for each sublattice index L∈{A,B},

For the 2D unit cell displayed in Fig. 3a, D2 symmetry is considered.

The spin components of the GS configuration depend on the α and variational parameters as

For magnetic fields parallel to the c axis, the GS energy only depends on the parameters α=α=α and . is the angle between neighbouring in-plane spins and α describes how much the spins are tilted away from the c axis. For α=π/2 the spins are lying within the ab plane, while α=0 corresponds to the fully polarized high-magnetic-field case. The GS energy

is minimized when

For zero field we find that α=0, the spins are confined to the ab plane (that is, to the IrO2 layer) and the angle is controlled by the strength of D. is not affected by fields along the c axis while α changes smoothly from π/2 to 0 with increasing the field strength.

The two magnons characteristic for spin-1/2 antiferromagnets are related to states orthogonal to and . The 4 × 4 Hamiltonian defining these magnons can be derived by using the well-known Holstein-Primakoff approach and diagonalized through a Bogoliubov transformation. For h||c, the two spin-wave modes are the gapless v||1=0 Goldstone mode, corresponding to U(1) symmetry breaking, and the gapped mode given by (2). When the magnetic field lies in the ab plane, the z component of Ir spins remains zero (with α=π/2) and the GS energy only depends on the angle :

For simplicity, we select x for the direction of the magnetic field. The result is, however, independent of how this choice is made as there is no anisotropy within the ab plane.

An infinitesimally small in-plane field fixes the direction of uniform magnetization as . At finite field both modes thus become gapped. To determine the uniform magnetization m one needs to minimize equation (21), which to leading order in magnetic field leads to equation (5). The first term of equation (5) corresponds to the zero-field moment, which arises due to the canting induced by the DM interaction. This ferromagnetic order-parameter is further enhanced in finite h⊥ field. As long as h⊥ is small equation (5) remains valid.

As discussed in the main text, in low fields, m can be approximated by its field independent value26. Using the quantum chemically derived interaction parameters (see Table 2), we then find m≈0.12μB, in good agreement with recent experiments840. Yet the zero-field gap comes out too large as compared to experiment. A good fit can nevertheless be reached with J, D and g values as obtained in the MRCI treatment and by increasing Γ from 0.42 to 0.98 meV.

Additional information

How to cite this article: Bogdanov, N. A. et al. Orbital reconstruction in nonpolar tetravalent transition-metal oxide layers. Nat. Commun. 6:7306 doi: 10.1038/ncomms8306 (2015).