Article

# Probing the Magnetic Ground State of Ba₂YIrO₆: Impact of Nonmagnetic Dopants and Spin-Orbit Coupling

Shuvajit Halder <fnref n="1" />, Md Salman Khan <fnref n="1" />, Fabrice Bert <fnref n="2" />, Payel Aich <fnref n="3" />,<fnref n="4" />, Carlo Meneghini <fnref n="3" />,* and Sugata Ray <fnref n="1" />,<fnref n="5" />

**Abstract:** Strong spin-orbit coupling (SOC) in iridates has long been predicted to lead to exotic electronic and magnetic ground states. Ba₂YIrO₆ (BYIO) has attracted particular attention due to the expectation of a *J*<sub>eff</sub> = 0 state for Ir<sup>5+</sup> ions under the *jj*-coupling scheme. However, controversies surround the actual realization of this state, as finite magnetic moments are consistently observed experimentally. We present a multi-physics study of this system by progressively introducing nonmagnetic Sb<sup>5+</sup> ions in place of Ir<sup>5+</sup> (Ba₂YIr<sub>1-y</sub>Sb<sub>y</sub>O<sub>6</sub>, BYISO). Despite similar charge and ionic radii, Sb<sup>5+</sup> doping appears highly inhomogeneous, coexisting with a fraction of nearly pure BYIO regions, as confirmed by X-ray diffraction (XRD). This aligns with observations in related compounds. While inhomogeneity creates uncertainty, the doped majority phases offer valuable insights. It is relevant that the inclusion of even small amounts of Sb<sup>5+</sup> (10–20%) leads to a rise in magnetization. This strengthens our previous suggestion that magnetic Ir ions form dynamic singlets in BYIO, resulting in a near-nonmagnetic background. The observed moment enhancement with nonmagnetic doping supports the breakdown of these singlets. Furthermore, the magnetization steadily increases with an increasing Sb<sup>5+</sup> content, contradicting the anticipated approach towards the *J*<sub>eff</sub> = 0 state with increased SOC due to reduced hopping between Ir<sup>5+</sup> ions. This reinforces the presence of individual Ir<sup>5+</sup> moments. Overall, our findings suggest that Ba₂YIrO₆ might not possess sufficiently strong SOC to be solely described within the *jj*-coupling picture, paving the way for further investigation.

**Keywords:** iridates; spin-orbit coupling; local-structure; *jj*-coupling; doping; magnetism

<figure><img src="image_1.png" /></figure>

check for
updates

**Citation:** Halder, S.; Khan, M.S.; Bert, F.; Aich, P.; Meneghini, C.; Ray, S. Probing the Magnetic Ground State of Ba₂YIrO₆: Impact of Nonmagnetic Dopants and Spin-Orbit Coupling. *Materials* **2024**, *17*, 1766. https://doi.org/10.3390/ma17081766

Academic Editor: George Kioseoglou

Received: 3 March 2024

Revised: 29 March 2024

Accepted: 6 April 2024

Published: 11 April 2024

<figure><img src="image_2.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Exotic electronic and magnetic ground states due to the influence of strong spin-orbit coupling have been the subject of intensive study for more than a decade now. As the spin-orbit coupling is proportional to the lower power of the atomic number (proportional to Z<fnref n="2" />-Z<fnref n="4" />), the 4d and 5d transition metals have become the central focus of study, especially the 5d iridium-based compounds [1–5]. The interplay between the strong spin-orbit coupling (SOC) and comparable onsite Coulomb interaction (U), along with the crystal field effect (Δ<sub>CFE</sub>), inter-site hopping (*t*<sub>ij</sub>), Hund's coupling (*J*<sub>H</sub>), and superexchange interaction energy 4t<fnref n="2" />/U, drives the system into many rich quantum mechanical states [3,6], such as Mott insulators [7], Weyl semimetals [8], quantum spin liquids [9,10], topological insulators [3,8], etc. In the strong spin-orbit coupling regime, *m*<sub>j</sub> becomes the only good quantum number instead of *m*<sub>1</sub> (orbital) and *m*<sub>s</sub> (spin), where the total angular momentum *J* determines the multiplates and degeneracy in the system.

<footnote n="1">¹ School of Materials Sciences, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, India; mssr@iacs.res.in</footnote>

<footnote n="2">² Laboratoire de Physique des Solides, UMR CNRS 8502, Université Paris-Sud, 91405 Orsay, France</footnote>

<footnote n="3">³ Dipartimento di Scienze, University Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy; payelaich2104@gmail.com</footnote>

<footnote n="4">⁴ LASR3 Surface Analysis Laboratory, University Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy</footnote>

<footnote n="5">⁵ Technical Research Center, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, India</footnote>

<footnote>* Correspondence: carlo.meneghini@uniroma3.it</footnote>
The study of iridates increased manifold after the reports of an SOC-driven Mott insulating ground state in layered Sr₂IrO₄ [6] and Na₂IrO₃ [11], which contradicts the uncorrelated band metallicity in iridates.

In the single-particle picture, large SOC splits the six-fold degenerate $t_{2g}$ orbitals of tetravalent iridium (Ir⁴⁺, 5d⁵ with one hole) into four-fold degenerate fully filled $J_{eff} = 3/2$ and doubly degenerate half-filled $J_{eff} = 1/2$ states. The $J_{eff} = 1/2$ band undergoes a Mott transition due to a relatively small (compared to the bandwidth) on-site Coulomb repulsion U [2,3,6,12,13]. This contrasts with pentavalent iridates (Ir⁵⁺, 5d⁴ with two $t_{2g}$ holes), where SOC leads to various (15) states, with a $J_{eff} = 0$ nonmagnetic state as the ground state. Surprisingly, to date, such an unusual nonmagnetic state with two unpaired electrons in the $t_{2g}$ band has never been realized in any kind of Ir⁵⁺ compounds [10,14–19]. Such deviations are often explained using different existing solid-state effects, such as the non-cubic crystal field ($\Delta_{CEF}^{NC}$), which modifies the effective SOC [11,14,20,21], the ligand–metal charge transfer, inter-site hopping, or ionic disorder [22,23], which can modify the SOC description at the atomic level and introduce small magnetic moments [7,10,14,17,19,24–26]. Particularly, the proposal of condensation of Van Vleck excitons, when the SOC strength and the superexchange energy scale become comparable, has been widely accepted as a feasible mechanism for moment development [24].

In this context, the perfectly cubic perovskite Ba₂YIrO₆ (BYIO) containing Ir⁵⁺ generated quite a bit of curiosity. Surprisingly, despite being free of any non-cubic crystal field ($\Delta_{CEF}^{NC} = 0$) effect, various studies have consistently reported finite magnetic moments for BYIO ranging from 0.16 µB/Ir to 1.44 µB/Ir with no long-range magnetic order observed down to the lowest measured temperature of 60 mK [17,19,26–28]. However, the origin of these observed magnetic moments is a subject of intense debate. A significant body of research, encompassing both experimental and theoretical investigations, has claimed that BYIO is truly a nonmagnetic $J_{eff} = 0$ system. These researchers attribute the observed magnetization solely to sample-related issues [1, 22, 23, 27, 29–32]. Contrarily, others argue that BYIO is intrinsically magnetic like other iridates, with each Ir⁵⁺ ion retaining a finite moment. [17,26,28,33–35]. Our group previously concluded the existence of hopping-induced finite intrinsic Ir moments in BYIO, forming fluctuating nonmagnetic Ir-Ir singlets with no long-range order (Figure 1a), resulting in a net nonmagnetic state between 60 mK and 10 K, as confirmed by µ-SR measurements [19]. In order to validate this hypothesis, earlier we attempted to weaken inter-site hopping and potentially drive the system towards the predicted atomic SOC-driven nonmagnetic $J_{eff} = 0$ ground state, by diluting Ir⁵⁺ ions in BaYSbO₆ (Ba₂YIr₁⁻ₓSbₓO₆ with y = 0.7–0.9), which has the same cubic structure as BYIO, but instead of Ir⁵⁺ there is Sb⁵⁺, which is not magnetic [16]. In this work we denote the content of Sb in BYIO by y, being complementary to the previous study [16], which focused on the dilution of Ir ions into a nonmagnetic structure, in which we denoted the content of Ir by x = 1 – y). However, this previous study [16] revealed an unexpected trend: the magnetic moments, instead of decreasing, exhibited a systematic increase with an increasing Sb⁵⁺ concentration (y). This unanticipated result raises concerns about the true strength of Ir⁵⁺ SOC in BYIO and necessitates a deeper understanding.

In the present study, we probe the low Sb doping regime, i.e., y = 0.1, 0.2, and 0.5, with an aim to progressively disrupt the nonmagnetic singlet pairs [19] in Ba₂YIrO₆ (BYIO) (Figure 1b), in order to confirm their existence. Clearly, the inclusion of nonmagnetic Sb⁵⁺ in BYIO is expected to break proportionate numbers of Ir⁵⁺-Ir⁵⁺ antiferromagnetically coupled singlet pairs, leaving an increasing number of unpaired Ir⁵⁺ ions around the dopant Sb⁵⁺ and consequently increasing the net magnetization by the free Ir⁵⁺ moments, if at all. For the sake of completeness, we have also included a few experimental results of our earlier reported y = 0.8 compound (x = 1 – y = 0.2 in ref. [16]) in the present manuscript. Consistent with our previous observation [16], the partial immiscibility problem of Sb⁵⁺ within the BYIO matrix has been observed here too, especially when the doping is taken above 10%. However, even then we actually obtain a major Ba₂YIr₁⁻ₓSbₓO₆ phase with y < y', which can still be utilized for the stated purpose. Consequently, we have carried out
detailed structural, electronic, and magnetic studies on the series of compounds, and our
detailed experiments reveal that the effective paramagnetic moments continuously increase
with the addition of nonmagnetic Sb<sup>5+</sup> ions in the system (from ~0.4 µ<sub>B</sub>/Ir to ~0.7 µ<sub>B</sub>/Ir),
with Θ<sub>CW</sub> being always negative [16], and, here again, similar to the parent Ba<sub>2</sub>YIrO<sub>6</sub> and
high Sb-doped end systems [16], no long-range magnetic order develops. Valence band
photoemission spectroscopic data established that the bandwidth of hybridized Ir 5d-O 2p
decreases with increased Sb doping, as expected, and the consequent increase in magnetic
moments only points towards the applicability of the localized moment LS coupling model
instead of the proposed *jj*-coupling interaction.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Crystal structures of (a) BYIO and (b) Sb-doped BYISO, showcasing the disruption caused by Sb doping (colour code: Green ball: Barium, Deep Green: Yttrium, Brown: Iridium, Deep Pink: Antimony, and Maroon: Oxygen). Schematic diagrams of (c) resonating valence bond (RVB) singlets (Cyan and Dark Cyan colour gradient) in BYIO and (d) antiferromagnetic Ir<sup>5+</sup> (Cyan ball) singlets disrupted by nonmagnetic Sb dopants (Red ball).</figcaption></figure>

**2. Experimental Details**

Poly crystalline samples of the Ba₂YIr₁₋₃Sb{Y}O₆ with doping y = 0.1, 0.2, 0.5, and 0.8 (abbreviated as BYISO-10, BYISO-20, BYISO-50, and BYISO-80, respectively) were synthesized by the conventional solid-state reaction technique. Stoichiometric amounts of high purity (>99%, Sigma-Aldrich, St. Louis, MO, USA) BaCO₃, Y₂O₃, IrO₂, and Sb₂O₅ powders were thoroughly grounded in an agate mortar. Initially, the mixture was calcined at 1173 K for 12 h to decompose carbonates. The mixture was then pressed into pellets and annealed in air at 1623 K for 72 h with a few intermittent grindings. The structural characterization and phase purity of all the compounds were checked in a Rigaku SmartLab X-ray Diffractometer (Tokyo, Japan) with a Cu Kα (λ = 1.5406 Å) X-ray source at room temperature (300 K). The powder X-ray diffraction (PXRD) data were analyzed to extract the structural information through Rietveld refinement using the Fullprof software [36,37]. The X-ray photoemission spectroscopy (XPS) experiments were carried out in an OMICRON electron spectrometer (Taunusstein, Germany), equipped with a SCIENTA OMICRON SPHERA analyzer and an Al Kα monochromatic X-ray source with an energy resolution of 0.5 eV. The in situ argon sputtering was used to clean the surface of the pellets. The Ir L₃ edge (11,215 eV) X-ray absorption fine structure (XAFS) experiments were performed at the XAFS beamline of an Elettra synchrotron radiation facility in Italy [38]. The Si(111) double
crystal was used to scan the X-ray beam energy across the Ir L3 edge. The absorption spectra were measured at room temperature in transmission geometry, using two gas-filled ionization chambers to measure incident and transmitted X-ray fluxes. The BYISO samples were grounded, mixed with boron nitride (BN) matrix in an approximately 1/10 weight ratio, pressed in thin solid pellets, and mounted on the beamline measurement chamber. The absorption edge discontinuity was approximately 0.5 for all the samples. The absorption spectra from a pure Ir foil placed after the second ionization chamber were measured at the same time and used to precisely monitor the X-ray beam energy calibration. The extended X-ray absorption fine structure (EXAFS) data were analyzed quantitatively using the open-source DEMETER (Athena and Artemis) [39,40] and ESTRA-FitEXA [41] software packages. The magnetization measurements in the temperature range of 2–300 K and ±5 T were performed in a superconducting quantum interference device (SQUID) magnetometer (Quantum Design). The muon spin resonance (μ-SR) experiments were performed using the muon spectrometer at the ISIS Neutron and Muon source facility in the United Kingdom.

### 3. Results and Discussions

#### 3.1. *Crystal Structure from X-ray Diffraction*

Room temperature X-ray diffraction patterns and the best-fitted Rietveld analysis curves for all the polycrystalline Ba₂YIr₁⁻ₓSbᵧO₆ (y = 0.1, 0.2, 0.5, and 0.8) compounds are shown in Figure 2a-d along with the high 2θ angle peak (inset to Figure 2a-d), which signify the influence of Sb on crystallographic phases as a function of doping concentrations, and the refined crystallographic lattice parameters and the phase percentages are listed in Table 1. The structural analysis reveals that the BYISO-10 has a single cubic phase, the space group is *Fm$\bar{3}$m, and as the doping percentage of Sb increases, the signature of two phases becomes more prominent (see Figure 2). Above 10% Sb doping, a dominating phase Ba₂YIr₁⁻ₓSbᵧO₆ (y < y', phase I) and a minor phase Ba₂YIrO₆ (phase II), which both have the cubic space group *Fm$\bar{3}$m [42], have been found to coexist in all the compounds. As the slightly large cation Sb⁵⁺ (< r >Sb⁵⁺ = 0.60 Å) is doped in place of Ir⁵⁺ (< r >Ir⁵⁺ = 0.57 Å), the cubic lattice constants are increased accordingly [43] (Table 1). The reason behind this inhomogeneous replacement of Ir⁵⁺ by Sb⁵⁺ is not clear. Given the fact that Y³⁺ is substantially larger (< r >Y³⁺ = 0.90 Å) than Ir⁵⁺/Sb⁵⁺, along with charge differences (see Table 2), B/B' anti-site disorder is found to be only marginal in both the phases. However, the nature of Sb⁵⁺-O-Y³⁺-O-Sb⁵⁺ bonding and Ir⁵⁺-O-Y³⁺-O-Ir⁵⁺ bonding may differ because of differences in *p*-block Sb⁵⁺ and Ir⁵⁺ *d*-block ions and could create certain local preferences giving rise to this partial immiscibility.

<table><caption>**Table 1.** Room temperature crystallographic information of Ba₂YIr₁⁻ₓSbᵧO₆ (y = 0.1, 0.2, 0.5, and 0.8): two cubic phases of space group *Fm$\bar{3}$m (space group no. 225, *a* = *b* = *c*, *α* = *β* = *γ* = 90°) are used to refine the powder XRD patterns. BYISO-10: Rₚ = 10.4, R<sub>wp</sub> = 7.8, R<sub>exp</sub> = 4.76, χ² = 2.68; BYISO-20: Rₚ = 11.1, R<sub>wp</sub> = 8.08, R<sub>exp</sub> = 5.9, χ² = 1.87. BYISO-50: Rₚ = 13.1, R<sub>wp</sub> = 8.33, R<sub>exp</sub> = 6.81, χ² = 1.49. BYISO-80: Rₚ = 19.6, R<sub>wp</sub> = 16.3, R<sub>exp</sub> = 7.19, χ² = 5.14 [16]. Standard uncertainty on the last digit of refined parameters is reported in parentheses. Values of the fixed or constrained parameters have no uncertainty.</caption><thead><tr><th>Sample</th><th>Phase (%)</th><th>a (Å)</th><th>Atom</th><th>Occupancy</th><th>x</th><th>y</th><th>z</th><th>B (Å)<sup>2</sup></th></tr></thead><tbody><tr><td rowspan="5">BYISO-10</td><td rowspan="5">Ba₂YIr₀.₉Sb₀.₁O₆ (100%)</td><td rowspan="5">8.356 (1)</td><td>Ba</td><td>1</td><td>0.25</td><td>0.25</td><td>0.25</td><td>0.357 (4)</td></tr><tr><td>Y</td><td>1</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.276 (1)</td></tr><tr><td>Ir</td><td>0.90</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.225 (3)</td></tr><tr><td>Sb</td><td>0.10 (1)</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.225</td></tr><tr><td>O</td><td>1</td><td>0.263 (5)</td><td>0.00</td><td>0.00</td><td>0.205 (2)</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th>Sample</th><th>Phase (%)</th><th>a (Å)</th><th>Atom</th><th>Occupancy</th><th>x</th><th>y</th><th>z</th><th>B (Å)<sup>2</sup></th></tr></thead><tbody><tr><td rowspan="10">BYISO-20</td><td rowspan="5">Ba<sub>2</sub>YIr<sub>0.78</sub>Sb<sub>0.22</sub>O<sub>6</sub> (91%)</td><td rowspan="5">8.364 (1)</td><td>Ba</td><td>1</td><td>0.25</td><td>0.25</td><td>0.25</td><td>0.388 (2)</td></tr><tr><td>Y</td><td>1</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.306 (2)</td></tr><tr><td>Ir</td><td>0.78</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.295 (1)</td></tr><tr><td>Sb</td><td>0.22 (2)</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.295</td></tr><tr><td>O</td><td>1</td><td>0.262 (9)</td><td>0.00</td><td>0.00</td><td>0.191 (4)</td></tr><tr><td rowspan="5">Ba<sub>2</sub>YIrO<sub>6</sub> (9(1)%)</td><td rowspan="5">8.357 (1)</td><td>Ba</td><td>1</td><td>0.25</td><td>0.25</td><td>0.25</td><td>0.188</td></tr><tr><td>Y</td><td>1</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.316</td></tr><tr><td>Ir</td><td>1</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.195</td></tr><tr><td>O</td><td>1</td><td>0.262 (1)</td><td>0.00</td><td>0.00</td><td>0.191</td></tr><tr><td>Ba<sub>2</sub>YIr<sub>0.41</sub>Sb<sub>0.59</sub>O<sub>6</sub> (85%)</td><td>8.384 (1)</td><td>Ba</td><td>1</td><td>0.25</td><td>0.25</td><td>0.25</td><td>0.363 (1)</td></tr><tr><td rowspan="10">BYISO-50</td><td rowspan="5">Ba<sub>2</sub>YIr<sub>0.41</sub>Sb<sub>0.59</sub>O<sub>6</sub> (85%)</td><td rowspan="5">8.384 (1)</td><td>Y</td><td>1</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.288 (3)</td></tr><tr><td>Ir</td><td>0.41</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.257 (3)</td></tr><tr><td>Sb</td><td>0.59 (3)</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.257</td></tr><tr><td>O</td><td>1</td><td>0.262 (2)</td><td>0.00</td><td>0.00</td><td>0.211 (5)</td></tr><tr><td>Ba</td><td>1</td><td>0.25</td><td>0.25</td><td>0.25</td><td>0.363</td></tr><tr><td rowspan="5">Ba<sub>2</sub>YIrO<sub>6</sub> (15(1)%)</td><td rowspan="5">8.359 (1)</td><td>Y</td><td>1</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.288</td></tr><tr><td>Ir</td><td>1</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.257</td></tr><tr><td>O</td><td>1</td><td>0.262 (3)</td><td>0.00</td><td>0.00</td><td>0.211</td></tr><tr><td>Ba<sub>2</sub>YIr<sub>0.12</sub>Sb<sub>0.88</sub>O<sub>6</sub> (91%)</td><td>8.408 (1)</td><td>Ba</td><td>1</td><td>0.25</td><td>0.25</td><td>0.25</td><td>0.284 (9)</td></tr><tr><td rowspan="10">BYISO-80</td><td rowspan="5">Ba<sub>2</sub>YIr<sub>0.12</sub>Sb<sub>0.88</sub>O<sub>6</sub> (91%)</td><td rowspan="5">8.408 (1)</td><td>Y</td><td>1</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.369 (7)</td></tr><tr><td>Ir</td><td>0.12</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.378 (8)</td></tr><tr><td>Sb</td><td>0.88 (1)</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.378</td></tr><tr><td>O</td><td>1</td><td>0.263 (3)</td><td>0.00</td><td>0.00</td><td>0.378 (8)</td></tr><tr><td>Ba</td><td>1</td><td>0.25</td><td>0.25</td><td>0.25</td><td>0.284</td></tr><tr><td rowspan="5">Ba<sub>2</sub>YIrO<sub>6</sub> (9(1)%)</td><td rowspan="5">8.356 (1)</td><td>Y</td><td>1</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.369</td></tr><tr><td>Ir</td><td>1</td><td>0.50</td><td>0.50</td><td>0.50</td><td>0.378</td></tr><tr><td>O</td><td>1</td><td>0.267 (1)</td><td>0.00</td><td>0.00</td><td>0.378</td></tr><tr><td>Ba<sub>2</sub>YIr<sub>0.12</sub>Sb<sub>0.88</sub>O<sub>6</sub> (91%)</td><td>8.408 (1)</td><td>Ba</td><td>1</td><td>0.25</td><td>0.25</td><td>0.25</td><td>0.284 (9)</td></tr><tr><td>Ba<sub>2</sub>YIr<sub>0.12</sub>Sb<sub>0.88</sub>O<sub>6</sub> (91%)</td><td>8.408 (1)</td><td>Y</td><td>1</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.369 (7)</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Room temperature X-ray diffraction (black open circle) along with Rietveld refined patterns (red solid line) of Ba<sub>2</sub>YIr<sub>1-y</sub>Sb<sub>y</sub>O<sub>6</sub> are shown in (a) BYISO-10, (b) BYISO-20, (c) BYISO-50, and (d) BYISO-80 panels, respectively. The expanded view of the higher angle contribution of *hkl* (620) of two phases (phase I: green line; phase II: pink line) is shown in the inset of corresponding panels.</figcaption></figure>
<table><caption>**Table 2.** Observed distances for Ir-O, Sb-O, and Y-O from XRD analysis of phase I in BYIO, BYISO-10, BYISO-20, BYISO-50, BYISO-80, and BYSO [16,19,44].</caption><thead><tr><th rowspan="2">Compounds</th><th colspan="2">Bond Length in Å</th></tr><tr><th>Ir(Sb)-O</th><th>Y-O</th></tr></thead><tbody><tr><td>BYIO</td><td>1.98</td><td>2.19</td></tr><tr><td>BYISO-10</td><td>1.97</td><td>2.20</td></tr><tr><td>BYISO-20</td><td>1.98</td><td>2.20</td></tr><tr><td>BYISO-50</td><td>1.99</td><td>2.19</td></tr><tr><td>BYISO-80</td><td>1.99</td><td>2.21</td></tr><tr><td>BYSO</td><td>1.96</td><td>2.24</td></tr></tbody></table>

## 3.2. Local Structure from Extended X-ray Absorption Fine Structure (EXAFS)

The Ir L$_{3}$ edge EXAFS data analysis has been carried out with the aim of revealing the details of the local coordination chemistry around the average Ir absorber and, in particular, the mid-range chemical order details [45], which is mandatory to shed light on the Ir-X (X = Ir, Y, Sb) correlations. To this aim, the quantitative analysis has been carried out by applying a multi-shell refinement procedure [46] to reproduce the main structural features around the average Ir absorber in the samples till around 4.5 Å. The EXAFS data analysis necessitates the careful identification of statistically relevant and physically meaningful structural signals. This task involves a trial-and-error process to establish meaningful constraints on fitting parameters, thereby minimizing correlations and improving the result reliability. Following the approach outlined in Ref. [16], we utilized local atomic clusters around Ir ions (derived from XRD data) to identify key single and multiple scattering contributions. These were then employed to calculate the photoelectron amplitude and scattering functions needed for theoretical EXAFS simulations [41].

The *k*-weighted EXAFS spectra for all the investigated samples are presented in Figure 3a, along with the best fit curves for the sake of comparison. The moduli of the Fourier transforms (|FTs|) of data and best fits are presented in Figure 3b, providing a more intuitive description of the average local atomic structure around Ir, which is a pseudoradial distribution function in which peaks represent the average interatomic distances (or photoelectron half path length for MS terms), which are roughly 0.5 Å squeezed by the phase shift effect. The most intense peak around 2 Å is the signature of the first Ir coordination shell, consisting of the six surrounding oxygen atoms (IrO$_{6}$ octahedra) with a coordination number (N$_{Ir-O}$) of 6. The very equal peak intensity and shape across all four doped samples suggest a highly similar Ir-O environment with minimal distortions in the IrO$_{6}$ octahedra, consistent with the XRD data analysis. The next evident peak in the |FT| at around 4 Å originates from the Ir next-nearest neighbors located along the perovskite cube edges. Notably, the multiple scattering (MS) contributions to this peak are significantly enhanced by the aligned Ir-O-X configurations. We accounted for Y/Ir chemical disorder, including the Ir and Y contributions whose multiplicities were by *x* and 1 − *x*, respectively, and *x* was refined. It should be noted that for all samples the best fit has *x* = 0, as expected for the ideal double perovskite, establishing the high degree of chemical order. Only the analysis of BYISO-80 reveals a fraction of antisite defects *x* = 0.09 ± 0.02, corresponding to averagely 5.5 Ir-Y and 0.5 Ir-Y next neighbors. This signifies that a marginal amount of antisite disorder between Ir and Y occurs only for the higher percentage of Sb doping (Table 3). We have also checked the possibility of Sb/Y antisite defects, but any attempt to include Ir-O-Sb contributions degrades the best fit quality, also in the the highest doped sample BYISO-80. The accuracy of EXAFS data analysis diminishes for shells further away either because structural disorder and the finite mean free path of the photoelectron attenuate the XAFS structural signals or because of the increasing number of interfering contributions, which lead to a complex overlap of signals from various atomic arrangements. Therefore, it is hard to comment specifically on eventual Sb substitution at the Ir position, which is around 6 Å (i.e., along the perovskite cube diagonal).
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Stacked plot of Ir L$_{3}$ edge EXAFS experimental and analysis data of BYISO-10, BYISO-20, BYISO-50, and BYISO-80; each graph is labeled accordingly: (a) *k*-weighted experimental data and the corresponding fits in the *k* range 3–12 Å$^{-1}$; (b) the Fourier transforms (moduli) of *k*-weighted EXAFS experimental data, magnitudes (FT), imaginary parts (Imm), and fitted curves.</figcaption></figure>

<table><caption>**Table 3.** Local structure parameters as obtained from the EXAFS analysis of the Ir L$_{3}$ edge for the four samples. The absolute mismatches between the experimental data and the best fit are R$^{2}$ = 0.038, 0.042, 0.047, and 0.05 for BYISO-10, BYISO-20, BYISO-50, and BYISO-80, respectively. Uncertainty on the last digit of the refined parameters is reported in parentheses. The values of fixed or constrained parameters have no uncertainties.</caption><thead><tr><th>Sample</th><th>Shell</th><th>N</th><th>σ<sup>2</sup> × 10<sup>-3</sup> (Å<sup>2</sup>)</th><th>R (Å)</th></tr></thead><tbody><tr><td rowspan="6">BYISO-10</td><td>Ir-O</td><td>6</td><td>0.18 (2)</td><td>1.961 (6)</td></tr><tr><td>Ir-O-O</td><td>24</td><td>0.9 (1)</td><td>3.31 (3)</td></tr><tr><td>Ir-Ba</td><td>8</td><td>13 (2)</td><td>3.62 (3)</td></tr><tr><td>Ir-Y (SS)</td><td>6</td><td>1.1 (1)</td><td>4.17 (3)</td></tr><tr><td>Ir-O-Y (MS-3 legs)</td><td>12</td><td>6.6</td><td>4.17</td></tr><tr><td>Ir-O-Y-O (MS-4 legs)</td><td>6</td><td>12.1 (1)</td><td>4.17</td></tr><tr><td rowspan="6">BYISO-20</td><td>Ir-O</td><td>6</td><td>0.27 (3)</td><td>1.968 (6)</td></tr><tr><td>Ir-O-O</td><td>24</td><td>1.1 (2)</td><td>3.36 (3)</td></tr><tr><td>Ir-Ba</td><td>8</td><td>14 (2)</td><td>3.64 (2)</td></tr><tr><td>Ir-Y (SS)</td><td>6</td><td>2.7 (2)</td><td>4.20 (2)</td></tr><tr><td>Ir-O-Y (MS-3 legs)</td><td>12</td><td>5.8</td><td>4.20</td></tr><tr><td>Ir-O-Y-O (MS-4 legs)</td><td>6</td><td>8.8 (4)</td><td>4.20</td></tr><tr><td rowspan="6">BYISO-50</td><td>Ir-O</td><td>6</td><td>0.15 (2)</td><td>1.962 (6)</td></tr><tr><td>Ir-O-O</td><td>24</td><td>1.2 (2)</td><td>3.31 (3)</td></tr><tr><td>Ir-Ba</td><td>8</td><td>12 (1)</td><td>3.63 (2)</td></tr><tr><td>Ir-Y (SS)</td><td>6</td><td>2.7 (2)</td><td>4.21 (2)</td></tr><tr><td>Ir-O-Y (MS-3 legs)</td><td>12</td><td>5.0</td><td>4.21</td></tr><tr><td>Ir-O-Y-O (MS-4 legs)</td><td>6</td><td>7.4 (3)</td><td>4.21</td></tr><tr><td rowspan="8">BYISO-80</td><td>Ir-O</td><td>6</td><td>0.13 (1)</td><td>1.958 (6)</td></tr><tr><td>Ir-O-O</td><td>24</td><td>2.0 (2)</td><td>3.25 (3)</td></tr><tr><td>Ir-Ba</td><td>8</td><td>10 (1)</td><td>3.53 (2)</td></tr><tr><td>Ir-Y (SS)</td><td>5.5 (1)</td><td>3.8 (2)</td><td>4.24 (2)</td></tr><tr><td>Ir-O-Y (MS-3 legs)</td><td>11</td><td>4.5</td><td>4.24</td></tr><tr><td>Ir-O-Y-O (MS-4 legs)</td><td>5.5</td><td>5.2 (2)</td><td>4.24</td></tr><tr><td>Ir-Ir (SS) (ASD)</td><td>0.5</td><td>9.5 (2)</td><td>4.24</td></tr><tr><td>Ir-O-Ir (MS-3 legs)</td><td>1.0</td><td>7.8</td><td>4.24</td></tr><tr><td></td><td>Ir-O-Ir-O (MS-4 legs)</td><td>1.0</td><td>5.3 (2)</td><td>4.24</td></tr></tbody></table>
### 3.3. Oxidation State and Valence Band Spectra from X-ray Photoelectron Spectroscopy

In order to discuss the predicted nonmagnetic ($J_{eff} = 0$) state in Ir$^{5+}$ ions under the influence of SOC (jj-coupling), it is important to confirm the oxidation state first [6,47,48]. X-ray photoelectron spectra of the Ir 4f core level from all three compounds could be fitted using a single spin-orbit doublet, as shown in Figure 4. The energy positions of 4f$_{7/2}$ (63.18 eV, 63.21 eV, 63.15 eV, and 63.13 eV) and 4f$_{5/2}$ (66.23 eV, 66.25 eV, 66.19 eV, and 66.17 eV) and their spin-orbit separations of around 3.05 eV, 3.04 eV, 3.04 eV, and 3.04 eV (listed in Table 4) for BYISO-10, BYISO-20, BYISO-50, and BYISO-80, respectively, confirm the presence of pure Ir$^{5+}$ only in all the compounds. The valence band (VB) spectra of these four compounds (shown in Figure 5) shows the absence of the density of states (DOS) at the Fermi level, indicating the insulating behavior in all of them that is exactly like the end members, Ba$_2$YIrO$_6$ and Ba$_2$YSbO$_6$. It is clearly seen from the VB spectra that as the Sb doping increases there is a gradual decrease in the DOS near to the Fermi level, indicating that this DOS is contributed mostly by the Ir 5d-band, and such a depletion makes the system more insulating, similar to Ba$_2$YSbO$_6$, which is a highly insulating dielectric compound [19,44,49].

<figure><img src="image_6.png" /><figcaption>**Figure 4.** X-ray photoelectron spectroscopic spectra of 4f core level (green circle) with corresponding fitting (pink solid line): (a) BYISO-10, (b) BYISO-20, and (c) BYISO-50, respectively.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Valence band spectra of BYISO-10 (violet circle), BYISO-20 (arctic circle), BYISO-50 (pink circle), and BYISO-80 (green circle), respectively.</figcaption></figure>
<table><caption>**Table 4.** X-ray photoelectron spectroscopy (XPS) data for the binding energy (eV) of *f*<sub>7/2</sub> orbital and spin-orbit separation of Ir<sup>5+</sup> element.</caption><thead><tr><td>Sample</td><td>Energy (eV)</td><td>SO Splitting (eV)</td></tr></thead><tbody><tr><td>BYISO-10</td><td>63.18</td><td>3.05</td></tr><tr><td>BYISO-20</td><td>63.21</td><td>3.04</td></tr><tr><td>BYISO-50</td><td>63.15</td><td>3.04</td></tr><tr><td>BYISO-80</td><td>63.13</td><td>3.04</td></tr></tbody></table>

### 3.4. Magnetic Susceptibility

The temperature-dependent (2 K to 300 K) dc magnetic susceptibility of BYISO-10, BYISO-20, BYISO-50, and BYISO-80 has been measured in zero field-cooled (ZFC) and field-cooled (FC) protocols with a 2000 Oe magnetic field, as shown in Figure 6a-c, respectively. The $\chi$ vs. T curves were fitted at the higher-temperature region using Curie-Weiss (CW) law: $\chi = \chi_0 + \frac{C}{T - \Theta_{CW}}$, where $\chi_0$, C, and $\Theta_{CW}$ represent the temperature-independent paramagnetic susceptibility, Curie constant, and Curie-Weiss temperature, respectively [50]. The linear fittings $1/(\chi - \chi_0)$ are satisfactory down to 150 K for BYISO-10, BYISO-20, and BYISO-50, below which they deviate from the paramagnetic behavior (inset Figure 6a-c). The dc magnetic susceptibility curves for all the compounds show no long-range ordering down to 2 K, and the $\Theta_{CW}$ values are largely negative, hence the frustration parameter $f = |\Theta_{CW}|/T_N$ values are high. The negative $\Theta_{CW}$ values also indicate antiferromagnetic interaction between the Ir<sup>5+</sup> ions in all the samples. However, contrary to the expectations, the net effective magnetic moments extracted from the CW fitting (listed in Table 5) are found to increase gradually with an increasing nonmagnetic Sb content, and the most isolated Ir<sup>5+</sup>-ion (in BYISO-80) possesses the highest moment (Figure 7a) [16]. This observation clearly indicates the following: (a) each Ir<sup>5+</sup> ion possesses a finite moment, and, as a result of breaking the singlets with Sb doping, the Ir free spins and net magnetization increase; (b) increasing the Ir dilution expectedly narrows down the Ir bandwidth, but instead of moving towards the atomic $J_{eff} = 0$ limit, it takes the system more and more away from it, further confirming the fact that the system is still better described by a moderate SOC picture.

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Magnetic susceptibility $\chi$ vs. T in zero field-cooled (black open circle) and field-cooled (red solid circle) contexts along with corresponding $1/(\chi - \chi_0)$ vs. T in the inset are plotted with Curie-Weiss fitting (red solid line): (a) BYISO-10, (b) BYISO-20, and (c) BYISO-50.</figcaption></figure>
<table><caption>**Table 5.** Effective paramagnetic moments from CW fitting and phase I moments calculated from XRD phase percentage (considering phase II moment ~0.44 µB/Ir) of BYIO [17], BYISO-10, BYISO-20, BYISO-50, and BYISO-80 [16].</caption><thead><tr><th rowspan="2">Sample</th><th colspan="2">Overall Compound</th><th colspan="2">Phase I</th></tr><tr><th>Θ<sub>CW</sub> (K)</th><th>μ<sub>eff</sub> (μ<sub>B</sub>/Ir)</th><th>Sb (%)</th><th>μ<sub>Phase I</sub> (μ<sub>B</sub>/Ir)</th></tr></thead><tbody><tr><td>BYISO-10</td><td>-96</td><td>0.45</td><td>10</td><td>0.45</td></tr><tr><td>BYISO-20</td><td>-76</td><td>0.47</td><td>22</td><td>0.47</td></tr><tr><td>BYISO-50</td><td>-44</td><td>0.57</td><td>59</td><td>0.59</td></tr><tr><td>BYISO-80</td><td>-31</td><td>0.67</td><td>88</td><td>0.69</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** (a) Net magnetic moment of Ba₂YIr₁⁻ₓSbᵧO₆ with respect to the stoichiometric Sb doping percentage; (b) magnetic moment of phase I as a function of the effective Sb doping percentage in phase I.</figcaption></figure>

Effective paramagnetic moments of phase I (listed in Table 5) are calculated using the phase percentage observed from the XRD analysis (Table 1) and considering the effective paramagnetic moment of BYIO, i.e., phase II as 0.44 µB/Ir [17].

### 3.5. Muon Spin Relaxation—μSR

From the macroscopic magnetic measurements, such as dc magnetic susceptibility measurement, it is hard to understand the complex magnetic ground state of these cubic double perovskite systems. Hence, to understand the magnetic nature of the BYISO system it becomes important to perform a μSR experiment, which is highly sensitive to the tiny internal magnetic field due to the large gyromagnetic ratio of the muon ($γμ = 851.615$ MHz/T). Polarization $P(t)$ variations with time at different temperatures under a zero applied field from BYISO-20 are shown in Figure 8a. It is observed that till the lowest measured temperature of 270 mK, the system is not magnetically ordered, almost similar to Ba₂YIrO₆ [19]. Under the application of zero external fields, the polarization at different temperatures follows the stretched exponential function, $P(t) = e^{-(λ't)β}$. Fitted values of the relaxation rate ($λ'$) and stretched exponent ($β$) as a function of temperature are plotted in Figure 8b and are almost similar to the parent compound Ba₂YIrO₆ [19]. As the temperature decreases, the $λ'$ increases gradually, and below 1 K it becomes constant till 270 mK. The characteristic spin fluctuation frequency $ν$, defined as $ν ∝ 1/λ'$, becomes nearly constant below 1 K till 270 mK. The spin dynamics signifies no magnetic freezing until the base temperature (270 mK). To explore the origin of the relaxation at low temperatures, the longitudinal field ($B_{LF}$)-dependent polarization has been measured at 270 mK. A significant variation in relaxation is observed with changes in the external longitudinal field $B_{LF}$ (shown in Figure 8c). To estimate the local field ($B_μ$) that the muons experience and the fluctuation frequency ($ν$), the Redfield formula (Equation (1)) is used to fit the relaxation ($λ'$) as a function of the external applied longitudinal field $B_{LF}$ (see Figure 8d, plotted $1/λ'$ as a function of $B_{LF}^2$).
$$ \lambda' = \lambda'_{0} + \frac{\nu\gamma_{\mu}^{2}B_{\mu}^{2}}{(\nu^{2} + \gamma_{\mu}^{2}B_{LF}^{2})} \quad (1) $$

The outcome from the Redfield formula yields $B_{\mu} \sim 2.2$ mT and $\nu \sim 0.1$ MHz, which strongly support the magnetization measurement data that 20% Sb doping in BYIO increases the internal magnetic field almost 10 times higher compared to the parent compound BYIO with less fluctuation [19]. Hence, it is clear that some resonating valence bond singlets are broken due to doping, giving rise to free Ir$^{5+}$ spins, which contribute to this enhanced local magnetic moment [51].

<figure><img src="image_10.png" /><figcaption>Figure 8. (a) Time evolution of the muon polarization of Ba$_{2}$YIr$_{0.8}$Sb$_{0.2}$O$_{6}$ in zero fields with fits to a stretched exponential function (continuous lines) with temperature variation; (b) fitted relaxation parameter $\lambda'$ as a function of temperature (T); (c) time evolution of the muon polarization with different applied longitudinal fields with fits to a stretched exponential function (continuous lines) at T = 270 mK; (d) $1/\lambda'$ vs. $B_{LF}^2$ for BYIO [19] (cyan color) and BYISO-20 (magenta color) data extracted using Redfield formula (Equation (1)).</figcaption></figure>

## 4. Conclusions

We conducted a comprehensive multi-physical characterization of Sb-doped Ba$_{2}$YIr$_{1-y}$Sb$_{y}$O$_{6}$ (BYISO) samples to investigate the nonmagnetic Ir$^{5+}$-Ir$^{5+}$ singlets and shed light on the magnetic state of isolated Ir$^{5+}$. Our structural investigations revealed the presence of defects, which are challenging to eliminate entirely in any real samples. Notably, XRD analysis showed a weak residual presence (around 10%) of a secondary BYIO phase (phase II). This results in a slight difference in the nominal composition of the Sb-doped phase (phase I), but such a difference does not exhibit a monotonic trend with the doping concentration that may motivate the evolution of magnetic properties (see below). The Ir-L$_{3}$ edge XAFS analysis did not reveal any defects in the Ir coordination (IrO$_{6}$ octahedra) or the presence of Ir-Sb correlations from Sb/Y substitutional defects, even at the highest Sb concentrations. A weak fraction (<10%) of Ir/Y antisite defects providing some Ir-O-Ir units was observed only at the highest Sb doping (BYISO-80).

The magnetic properties exhibit a clear progressive increase in the total sample magnetization, as well as the magnetization attributed solely to phase I (Sb-doped). Notably,
the magnetization is highest in the sample with the highest Sb concentration (BYISO-80), where the Ir ions should be completely isolated. This finding clearly indicates that each Ir⁵⁺ ion possesses an intrinsic magnetic moment. Breaking the singlets with Sb doping leads to an increase in free Ir spins and so in the net magnetization. As expected, increasing the Ir dilution narrows the Ir bandwidth. However, contrary to the expectation of approaching the atomic *J*<sub>eff</sub> = 0 limit, this effect pushes the system further away from it, reinforcing the idea that the system is better described by a moderate SOC picture. The μSR measurements further support this hypothesis. They show an order of magnitude increase in the internal magnetization field for the BYISO-20 sample compared to the pure BYIO sample. This effect evidently arises from the free Ir⁵⁺ spins originating from the Sb broken singlets, which enhance the local magnetic moments, against the *J*<sub>eff</sub> = 0 ground state model,

Overall our results clearly demonstrate the inadequacy of a *J*<sub>eff</sub> = 0 state for Ir<sup>5+</sup> and suggest that the SOC strength in BYIO is better described within the framework of LS coupling.

**Author Contributions:** S.H. synthesized and characterized the samples and largely carried out all the spectroscopic measurements. M.S.K. contributed in discussion and shared a part of the results related to Ref. [16]. F.B. helped to perform muon spectroscopy measurement and analysis. P.A. contributed to EXAFS analysis and manuscript editing. C.M. performed the XAS measurement and analysis and supervised the project. S.R. conceptualized the problem, supervised the overall work, and finalized the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** The authors acknowledge the financial support for XAFS experiments at the ELETTRA synchrotron radiation facility (exp. 20200283, 20205397). C.M. and P.A. acknowledge the financial support from the Grant of Excellence Departments 2023–2027, MIUR (ARTICOLO 1, COMMI 314–337 LEGGE 232/2016), and ECS 0000024 Rome Technopole (CUP B83C22002820006, PNRR Missione 4 Componente 2 Investimento 1.5, financed by NextGenerationEU). S.R. acknowledges SERB, DST, for financial support (Project No. CRG/2019/003522).

**Data Availability Statement:** XAFS data will be available at the ELETTRA-XAFS beamline after an embargo (see https://www.elettra.eu/userarea/scientific-data-policy.html (accessed on 1 January 2024) for the ELETTRA synchrotron data policy).

**Acknowledgments:** C.M. gratefully acknowledges the XAFS beamline staff at the ELETTRA Synchrotron Radiation Facility for their valuable technical support. S.H. acknowledges CSIR, India, for fellowship. The authors also thank the Department of Central Scientific Services (CSS) and Technical Research Center (TRC) of IACS for providing experimental facilities. The authors also thank P. Biswas and A. D. Hiller in ISIS, RAL, UK, for help in µSR experiments.

**Conflicts of Interest:** The authors declare no conflicts of interest.

**References**

1. Kusch, M.; Katukuri, V.M.; Bogdanov, N.A.; Büchner, B.; Dey, T.; Efremov, D.V.; Hamann-Borrero, J.E.; Kim, B.H.; Krisch, M.; Maljuk, A.; et al. Observation of heavy spin-orbit excitons propagating in a nonmagnetic background: The case of (Ba,Sr)₂YIrO₆. *Phys. Rev. B* **2018**, *97*, 064421. [CrossRef]

2. Rau, J.G.; Lee, E.K.H.; Kee, H.Y. Spin-Orbit Physics Giving Rise to Novel Phases in Correlated Systems: Iridates and Related Materials. *Annu. Rev. Condens. Matter Phys.* **2016**, *7*, 195–221. [CrossRef]

3. Witczak-Krempa, W.; Chen, G.; Kim, Y.B.; Balents, L. Correlated Quantum Phenomena in the Strong Spin-Orbit Regime. *Annu. Rev. Condens. Matter Phys.* **2014**, *5*, 57–82. [CrossRef]

4. Khomskii, D. *Transition Metal Compounds*; Cambridge University Press: Cambridge, UK, 2014.

5. Shanavas, K.; Popović, Z.S.; Satpathy, S. Theoretical model for Rashba spin-orbit interaction in d electrons. *Phys. Rev. B* **2014**, *90*, 165108. [CrossRef]

6. Kim, B.J.; Jin, H.; Moon, S.J.; Kim, J.Y.; Park, B.G.; Leem, C.S.; Yu, J.; Noh, T.W.; Kim, C.; Oh, S.J.; et al. Novel $J_{eff} = 1/2$ Mott State Induced by Relativistic Spin-Orbit Coupling in $Sr_2IrO_4$. *Phys. Rev. Lett.* **2008**, *101*, 076402. [CrossRef]

7. Svoboda, C.; Randeria, M.; Trivedi, N. Effective magnetic interactions in spin-orbit coupled $d^4$ Mott insulators. *Phys. Rev. B* **2017**, *95*, 014409. [CrossRef]

8. Wan, X.; Turner, A.M.; Vishwanath, A.; Savrasov, S.Y. Topological semimetal and Fermi-arc surface states in the electronic structure of pyrochlore iridates. *Phys. Rev. B* **2011**, *83*, 205101. [CrossRef]
9. Corredor, L.T.; Aslan-Cansever, G.; Sturza, M.; Manna, K.; Maljuk, A.; Gass, S.; Dey, T.; Wolter, A.U.B.; Kataeva, O.; Zimmermann, A.; et al. Iridium double perovskite Sr<sub>2</sub>YIrO<sub>6</sub>: A combined structural and specific heat study. *Phys. Rev. B* **2017**, *95*, 064418. [CrossRef]

10. Nag, A.; Middey, S.; Bhowal, S.; Panda, S.K.; Mathieu, R.; Orain, J.C.; Bert, F.; Mendels, P.; Freeman, P.G.; Mansson, M.; et al. Origin of the Spin-Orbital Liquid State in a Nearly *J* = 0 Iridate Ba<sub>3</sub>ZnIr<sub>2</sub>O<sub>9</sub>. *Phys. Rev. Lett.* **2016**, *116*, 097205. [CrossRef] [PubMed]

11. Kim, C.H.; Kim, H.S.; Jeong, H.; Jin, H.; Yu, J. Topological Quantum Phase Transition in 5d Transition Metal Oxide Na<sub>2</sub>IrO<sub>3</sub>. *Phys. Rev. Lett.* **2012**, *108*, 106401. [CrossRef] [PubMed]

12. Moon, S.J.; Jin, H.; Kim, K.W.; Choi, W.S.; Lee, Y.S.; Yu, J.; Cao, G.; Sumi, A.; Funakubo, H.; Bernhard, C.; et al. Dimensionality-Controlled Insulator-Metal Transition and Correlated Metallic State. *Phys. Rev. Lett.* **2008**, *101*, 226402. [CrossRef] [PubMed]

13. Kim, B.; Ohsumi, H.; Komesu, T.; Sakai, S.; Morita, T.; Takagi, H.; Arima, T.H. Phase-sensitive observation of a spin-orbital Mott state in Sr<sub>2</sub>IrO<sub>4</sub>. *Science* **2009**, *323*, 1329–1332. [CrossRef] [PubMed]

14. Cao, G.; Qi, T.F.; Li, L.; Terzic, J.; Yuan, S.J.; DeLong, L.E.; Murthy, G.; Kaul, R.K. Novel Magnetism of Ir<sup>5+</sup>(5d<sup>4</sup>) Ions in the Double Perovskite Sr<sub>2</sub>YIrO<sub>6</sub>. *Phys. Rev. Lett.* **2014**, *112*, 056402. [CrossRef] [PubMed]

15. Bandyopadhyay, A.; Chakraborty, A.; Bhowal, S.; Kumar, V.; Sala, M.M.; Efimenko, A.; Bert, F.; Biswas, P.K.; Meneghini, C.; Büttgen, N.; et al. Breakdown of atomic spin-orbit coupling picture in an apparently isolated pseudo-one-dimensional iridate: Sr<sub>3</sub>NaIrO<sub>6</sub>. *Phys. Rev. B* **2022**, *105*, 104431. [CrossRef]

16. Khan, M.S.; Meneghini, C.; Bert, F.; Sala, M.M.; Ray, S. Phase inhomogeneity and inapplicability of *J*<sub>eff</sub> = 0 description in dilute Ir<sup>5+</sup> doped Ba<sub>2</sub>YSb<sub>1-x</sub>Ir<sub>x</sub>O<sub>6</sub> (*x* = 0.1, 0.2, 0.3). *Phys. Rev. B* **2021**, *104*, 214414. [CrossRef]

17. Dey, T.; Maljuk, A.; Efremov, D.V.; Kataeva, O.; Gass, S.; Blum, C.G.F.; Steckel, F.; Gruner, D.; Ritschel, T.; Wolter, A.U.B.; et al. Ba<sub>2</sub>YIrO<sub>6</sub>: A cubic double perovskite material with Ir<sup>5+</sup> ions. *Phys. Rev. B* **2016**, *93*, 014434. [CrossRef]

18. Bremholm, M.; Dutton, S.; Stephens, P.; Cava, R. NaIrO<sub>3</sub>—A pentavalent post-perovskite. *J. Solid State Chem.* **2011**, *184*, 601–607. [CrossRef]

19. Nag, A.; Bhowal, S.; Chakraborty, A.; Sala, M.M.; Efimenko, A.; Bert, F.; Biswas, P.K.; Hillier, A.D.; Itoh, M.; Kaushik, S.D.; et al. Origin of magnetic moments and presence of spin-orbit singlets in Ba<sub>2</sub>YIrO<sub>6</sub>. *Phys. Rev. B* **2018**, *98*, 014431. [CrossRef]

20. Khomskii, D.; Kugel, K.; Sboychakov, A.; Streltsov, S. Role of local geometry in the spin and orbital structure of transition metal compounds. *J. Exp. Theor. Phys.* **2016**, *122*, 484–498. [CrossRef]

21. Dodds, T.; Choy, T.P.; Kim, Y.B. Interplay between lattice distortion and spin-orbit coupling in double perovskites. *Phys. Rev. B* **2011**, *84*, 104439. [CrossRef]

22. Hammerath, F.; Sarkar, R.; Kamusella, S.; Baines, C.; Klauss, H.H.; Dey, T.; Maljuk, A.; Gaß, S.; Wolter, A.U.B.; Grafe, H.J.; et al. Diluted paramagnetic impurities in nonmagnetic Ba<sub>2</sub>YIrO<sub>6</sub>. *Phys. Rev. B* **2017**, *96*, 165108. [CrossRef]

23. Chen, Q.; Svoboda, C.; Zheng, Q.; Sales, B.C.; Mandrus, D.G.; Zhou, H.D.; Zhou, J.S.; McComb, D.; Randeria, M.; Trivedi, N.; et al. Magnetism out of antisite disorder in the *J* = 0 compound Ba<sub>2</sub>YIrO<sub>6</sub>. *Phys. Rev. B* **2017**, *96*, 144423. [CrossRef]

24. Khaliullin, G. Excitonic Magnetism in Van Vleck-type *d*<sup>4</sup> Mott Insulators. *Phys. Rev. Lett.* **2013**, *111*, 197201. [CrossRef] [PubMed]

25. Laguna-Marco, M.A.; Kayser, P.; Alonso, J.A.; Martínez-Lope, M.J.; van Veenendaal, M.; Choi, Y.; Haskel, D. Electronic structure, local magnetism, and spin-orbit effects of Ir(IV)-, Ir(V)-, and Ir(VI)-based compounds. *Phys. Rev. B* **2015**, *91*, 214433. [CrossRef]

26. Bhowal, S.; Baidya, S.; Dasgupta, I.; Saha-Dasgupta, T. Breakdown of *J* = 0 nonmagnetic state in *d*<sup>4</sup> iridate double perovskites: A first-principles study. *Phys. Rev. B* **2015**, *92*, 121113. [CrossRef]

27. Pajskr, K.; Novák, P.; Pokorný, V.; Kolorenč, J.; Arita, R.; Kuneš, J. On the possibility of excitonic magnetism in Ir double perovskites. *Phys. Rev. B* **2016**, *93*, 035129. [CrossRef]

28. Terzic, J.; Zheng, H.; Ye, F.; Zhao, H.D.; Schlottmann, P.; De Long, L.E.; Yuan, S.J.; Cao, G. Evidence for a low-temperature magnetic ground state in double-perovskite iridates with Ir<sup>5+</sup>(5d<sup>4</sup>) ions. *Phys. Rev. B* **2017**, *96*, 064436. [CrossRef]

29. Gong, H.; Kim, K.; Kim, B.H.; Kim, B.; Kim, J.; Min, B. Is the ground state of 5d<sup>4</sup> double-perovskite Iridate Ba<sub>2</sub>YIrO<sub>6</sub> magnetic or nonmagnetic? *J. Magn. Magn. Mater.* **2018**, *454*, 66–70. [CrossRef]

30. Fuchs, S.; Dey, T.; Aslan-Cansever, G.; Maljuk, A.; Wurmehl, S.; Büchner, B.; Kataev, V. Unraveling the Nature of Magnetism of the 5 d 4 Double Perovskite Ba 2 YIrO 6. *Phys. Rev. Lett.* **2018**, *120*, 237204. [CrossRef] [PubMed]

31. Kim, B.H.; Efremov, D.V.; van den Brink, J. Spin-orbital excitons and their potential condensation in pentavalent iridates. *Phys. Rev. Mater.* **2019**, *3*, 014414. [CrossRef]

32. Aczel, A.A.; Chen, Q.; Clancy, J.; dela Cruz, C.; Reig-i Plessis, D.; MacDougall, G.; Pollock, C.; Upton, M.; Williams, T.; LaManna, N.; et al. Spin-orbit coupling controlled ground states in the double perovskite iridates A<sub>2</sub>BIrO<sub>6</sub> (A = Ba, Sr; B = Lu, Sc). *Phys. Rev. Mater.* **2022**, *6*, 094409. [CrossRef]

33. Paramekanti, A.; Singh, D.J.; Yuan, B.; Casa, D.; Said, A.; Kim, Y.J.; Christianson, A.D. Spin-orbit coupled systems in the atomic limit: Rhenates, osmates, iridates. *Phys. Rev. B* **2018**, *97*, 235119. [CrossRef]

34. Singh, B.; Cansever, G.; Dey, T.; Maljuk, A.; Wurmehl, S.; Büchner, B.; Kumar, P. Orbiton-phonon coupling in Ir<sup>5+</sup> (5d<sup>4</sup>) double perovskite Ba<sub>2</sub>YIrO<sub>6</sub>. *J. Phys. Condens. Matter* **2018**, *31*, 065603. [CrossRef] [PubMed]

35. Schnait, H.; Bauernfeind, D.; Saha-Dasgupta, T.; Aichhorn, M. Small moments without long-range magnetic ordering in the zero-temperature ground state of the double perovskite iridate Ba<sub>2</sub>YIrO<sub>6</sub>. *Phys. Rev. B* **2022**, *106*, 035132. [CrossRef]

36. Rodríguez-Carvajal, J. Recent advances in magnetic structure determination by neutron powder diffraction. *Phys. B Condens. Matter* **1993**, *192*, 55–69. [CrossRef]
37. Rodríguez-Carvajal, J. Recent developments of the program FULLPROF. *Comm. Powder Diffr. (IUCr)* **Newsl. 2001**, *26*, 12–19.

38. Cicco, A.D.; Aquilanti, G.; Minicucci, M.; Principi, E.; Novello, N.; Cognigni, A.; Olivi, L. Novel XAFS capabilities at ELETTRA synchrotron light source. *J. Phys. Conf. Ser.* **2009**, *190*, 012043. [CrossRef]

39. Ravel, B.; Newville, M. *ATHENA, ARTEMIS, HEPHAESTUS*: Data analysis for X-ray absorption spectroscopy using *IFEFFIT*. *J. Synchrotron Radiat.* **2005**, *12*, 537–541. [CrossRef] [PubMed]

40. Newville, M. *IFEFFIT*: Interactive XAFS analysis and *FEFF* fitting. *J. Synchrotron Radiat.* **2001**, *8*, 322–324. [CrossRef]

41. Meneghini, C.; Bardelli, F.; Mobilio, S. ESTRA-FitEXA: A software package for EXAFS data analysis. *Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At.* **2012**, *285*, 153–157. [CrossRef]

42. Fu, W.; IJdo, D. Re-examination of the structure of Ba₂MIrO₆ (M=La, Y): Space group revised. *J. Alloys Compd.* **2005**, *394*, L5–L8. [CrossRef]

43. Shannon, R.D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. *Acta Crystallogr. Sect. A* **1976**, *32*, 751–767. [CrossRef]

44. Wu, L.; Mei, X.; Zheng, W. Hydrothermal synthesis and characterization of double perovskite Ba₂YSbO₆. *Mater. Lett.* **2006**, *60*, 2326–2330. [CrossRef]

45. Bandyopadhyay, A.; Neogi, S.K.; Atanu, P.; Meneghini, C.; Daasgupta, I.; Bandyopadhyay, S.; Ray, S.R. Atomic-scale chemical fluctuation in LaSrVMoO₆, a proposed half-metallic antiferromagnet. *Phys. Rev. B* **2017**, *95*, 024432. [CrossRef]

46. Jana, S.; Meneghini, C.; Sanyal, P.; Sarkar, S.; Saha-Dasgupta, T.; Karis, O.; Ray, S. Signature of an antiferromagnetic metallic ground state in heavily electron-doped Sr₂FeMoO₆. *Phys. Rev. B* **2012**, *86*, 054433. [CrossRef]

47. Cao, G.; Durairaj, V.; Chikara, S.; DeLong, L.E.; Parkin, S.; Schlottmann, P. Non-Fermi-liquid behavior in nearly ferromagnetic SrIrO₃ single crystals. *Phys. Rev. B* **2007**, *76*, 100402. [CrossRef]

48. Cao, G.; Subedi, A.; Calder, S.; Yan, J.Q.; Yi, J.; Gai, Z.; Poudel, L.; Singh, D.J.; Lumsden, M.D.; Christianson, A.D.; et al. Magnetism and electronic structure of La₂ZnIrO₆ and La₂MgIrO₆: Candidate *J*<sub>eff</sub> = 1/2 Mott insulators. *Phys. Rev. B* **2013**, *87*, 155136. [CrossRef]

49. Konopka, J.; Jose, R.; Wołcyrz, M. Structural and dielectric properties of Ba₂YbTaO₆, Ba₂YSbO₆ and Ba₂EuZrO₅.₅. *Phys. C Supercond. Appl.* **2006**, *435*, 53–58. [CrossRef]

50. Nag, A.; Ray, S. Misjudging frustrations in spin liquids from oversimplified use of Curie-Weiss law. *J. Magn. Magn. Mater.* **2017**, *424*, 93–98. [CrossRef]

51. Uemura, Y.J.; Yamazaki, T.; Harshman, D.R.; Senba, M.; Ansaldo, E.J. Muon-spin relaxation in AuFe and CuMn spin glasses. *Phys. Rev. B* **1985**, *31*, 546–563. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.