arxiv: v1 [cond-mat.str-el] 17 Apr 2014

Transcription

1 Evolution of Magnetism in the Na 3 δ (Na 1 x Mg x )Ir 2 O 6 Series of Honeycomb Iridates David C. Wallace a,b, Craig M. Brown c,d, Tyrel M. McQueen a,b arxiv: v1 [cond-mat.str-el] 17 Apr 2014 a Department of Chemistry, The Johns Hopkins University, Baltimore, MD 21218, USA b Institute for Quantum Matter, Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA c NIST Center For Neutron Research, Gaithersburg, MD, 20899, USA d Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, USA Abstract The structural and magnetic properties of a new series of iridium-based honeycomb lattices with the formula Na 3 δ (Na 1 x Mg x )Ir 2 O 6 (0 x 1) are reported. As x and δ are increased, the honeycomb lattice contracts and the strength of the antiferromagnetic interactions decreases systematically due to a reduction in Ir O Ir bond angles. Samples with imperfect stoichiometry exhibit disordered magnetic freezing at temperatures T f between 3.4 K and 5 K. This glassy magnetism likely arises due to the presence of non-magnetic Ir 3+, which are distributed randomly throughout the lattice, with a possible additional contribution from stacking faults. Together, these results demonstrate that chemical defects and non-stoichiometry have a significant effect on the magnetism of compounds in the A 2 IrO 3 materials family. Keywords: Honeycomb, Iridium, Magnetic Frustration, Disorder 1. Introduction 5d transition metal oxide materials have recently attracted significant interest due to the comparable energy scales of crystal field stabilization, electronic correlations (Hubbard U), and spin-orbit coupling [1]. In particular, iridium oxides are of great interest due to the ability of Ir to adopt a d 5 electronic configuration in its 4+ oxidation state, yielding a singly-occupied orbital with non-trivial spin texture in the strong spin-orbit coupling limit. This electronic configuration has been predicted to give rise to an unusual pattern of superexchange interactions, resulting in magnetic frustration and possible spin liquid behavior [2, 3, 4]. For this reason, the A 2 IrO 3 (A = Li, Na) family of layered honeycomb materials have been the subject of intense scrutiny [5, 6, 7, 8, 9]. The majority of experimental work on Na 2 IrO 3 and Li 2 IrO 3 has focused on singlecrystalline specimens grown via self-flux. Despite strong antiferromagnetic interactions in Na 2 IrO 3, evidenced by a reported large negative Weiss temperature θ W = -116 K, longrange antiferromagnetic order is only present below a Néel temperature of T N = 15 K, address: mcqueen@jhu.edu (Tyrel M. McQueen) Preprint submitted to Journal of Solid State Chemistry April 18, 2014

2 demonstrating that that Na 2 IrO 3 is in a magnetically frustrated regime [10, 11]. A study of the (Na 1 x Li x ) 2 IrO 3 solid solution revealed that θ W as well as T N vary significantly with x, and the magnetism reaches a peak in frustration at the intermediate value of x = 0.7 [12]. Isovalent substitution of Li for Na in the A 2 IrO 3 structure presumably has the effect of modulating the exchange interaction strengths between neighboring Ir sites, however this effect has not been systematically explored. Recently, a broad variety of new ternary Ir oxides with the formula Na 3 δ MIr 2 O 6 (M = Zn, Cu, Ni, Co, Fe and Mn) were reported [13], but no unifying picture of the magnetism demonstrated by this family can yet be drawn. Here, we report another new member of this chemical family, Na 3 δ (Na 1 x Mg x )Ir 2 O 6, in which Mg substitutes for Na in the honeycomb planes and the Na content between planes is variable(fig. 1(a)). Aliovalent substitution of Mg for Na in the honeycomb plane offers a convenient and controllable method by which to tune the lattice parameters of the unit cell and, in doing so, vary the exchange interactions between neighboring Ir sites. Additionally, we resolve the discrepancy in magnetic behavior between polycrystalline and single-crystalline samples of Na 2 IrO 3, where the former are spin glasses and the latter exhibit long range antiferromagnetic (AFM) order [10], by showing that Na 2 IrO 3 decomposes rapidly in air as previously reported, and that this decomposition destroys AFM order [14]. 2. Materials and Methods Polycrystalline Na 2 IrO 3 was prepared via a method similar to what was previously reported[10]: Na 2 CO 3 (NOAHTechnologies Corp., 99.9%) 1, MgO(NOAHTechnologiesCorp., 99.99%) and Ir black (J&J Materials Inc.) were intimately ground using an agate mortar and pestle and pelletized. The sample was heated in a covered alumina crucible to 750 C over a period of 4 hrs, held at that temperature for 30 hrs, then quenched in air. The pellet was then reground with an additional 5 mol % Na 2 CO 3, pelletized, and placed into a furnace preheated to 900 C where it dwelled for and additional 48 hrs. After the final heating, samples were removed from the furnace and allowed to cool in an argon filled glovebox to prevent reaction with laboratoryair. Polycrystalline samples of Na 3 δ (Na 1 x Mg x Ir 2 )O 6 were prepared via a similar method: Na 2 CO 3, MgO and Ir metal powder were intimately ground and pelletized. Samples were heated in covered alumina crucibles to 750 C over a period of four hrs, held at that temperature for 30 hrs, then quenched in air. Samples were then reground, pelletized, and placed into a furnace preheated to 900 C where they dwelled for between 48 and 96 hrs, followed by quenching in air. The 900 C heating was repeated with intermediate regrindings until the laboratory x-ray powder diffraction (XRPD) data showed no impurity phases and no change in lattice parameters or relative peak intensities between subsequent heatings. After the final heating, samples were processed and stored in an argon-filled glovebox to prevent reaction with air. Samples of suitable size for neutron powder diffraction (NPD) were prepared via a similar method, but with an additional two week heating at 900 C. 1 Certain commercial suppliers are identified in this paper to foster understanding. Such identification does not imply recommendation or endorsement by the National Institute of Standards and Technology, nor does it imply that the materials or equipment identified are necessarily the best available for the purpose. 2

3 Laboratory XRPD data were collected on a Bruker D8 Focus diffractometer using Cu-Kα radiation with λ 1 = Å and λ 2 = Å and equipped with a LynxEye CCD detector. Rietveld refinements to laboratory XRPD data were performed using TOPAS (Bruker AXS). NPD data were collected at the NIST Center for Neutron Research using the BT-1 high-resolution powder diffractometer with an incident wavelength of λ = Å using a Cu (311) monochromator and with 15 and 60 for primary and secondary collimation, respectively. Rietveld refinements to NPD data were performed in GSAS/EXPGUI [16, 17]. Samples were loaded into vanadium cells and sealed in an inert helium glovebox to prevent exposure to air. DC and AC magnetization measurements were carried out using a Quantum Design Physical Properties Measurement System. Elemental analyses were performed via inductively-coupled plasma optical emission spectroscopy (ICP-OES) by Evans Analytical Group. XRPD simulations were performed using DIFFaX [18]. 3. Results and Discussion 3.1. Structure The structure of the A 2 IrO 3 family of materials is derived from the layered triangular lattice compound α NaFeO 2 (Fig. 1(a)): when Na substitutes for one third of the Fe sites in NaFeO 2 in an ordered fashion, a layered honeycomb lattice is formed and the chemical formula becomes Na 2 FeO 3. For this reason, the structural formula of A 2 IrO 3 can be described more intuitively as A 3 A Ir 2 O 6, where electronically active A Ir 2 O 6 honeycomb layers, comprised of tilted, edge-sharing octahedra, are separated by electronically inert layers of A cations. The honeycomb order in these layers is energetically favorable due to the large differences in atomic radii and oxidation states of Ir in comparison to the A cation. Because electronically active honeycomb layers are separated by electronically inert triangular layers, the A 2 IrO 3 family are essentially layered two dimensional materials. This fact can give rise both to interesting physics, and to complex structural disorder. Fig. 2(a) shows simulated XRPD data that demonstrate how disordered stacking arrangements affect the observed XRPD pattern. When honeycomb lattices are stacked in a perfectly ordered manner, sharp reflections are observed in the XRPD data between scattering angle 2θ = 19 and 32 these are often referred to as supercell reflections, as the arise from the additional honeycomb ordering within the triangular layer. This structure (Fig. 1(a)) can be accounted for by a larger R3m cell where the a axis is expanded by 3 compared to the triangular R3m cell, or to lower symmetry variations such as the monoclinic C2/m or C2/c cells commonly observed (Fig. 1(b)) [11, 19]. If stacking faults are introduced, even with a low probability of occurrence, the supercell reflections become broadened and attenuated as shown in Fig. 2(a) with 5% stacking faults. Fig. 2(b) shows XRPD data collected on three different samples of Na 2 IrO 3, each with a different degree of structural disorder. Rietveld refinement to data collected on a sample with stacking faults can incorrectly account for the observed loss in supercell reflection intensity by introducing cation site mixing in the honeycomb plane. The energy cost associated with producing a stacking fault in Na 2 IrO 3 and structurally analogous systems is very small, and is thus significantly more likely to occur than antisite disorder [11, 15]. Similarly, a refinement can give the illusion of small particle size in order to account for the observed broadening of the supercell reflections. The same effect is observed in neutron diffraction, though the honeycomb order may be less apparent 3

4 if atoms in the honeycomb layer have similar neutron scattering cross sections. It is therefore quite difficult to make accurate structural characterizations of samples with significant stacking disorder using neutron or x-ray diffraction techniques. Fig. 3 shows NPD datasets collected on Na 3 Mg 0.5 Ir 2 O 6 at T = 5, 100, and 295 K, along with Rietveld refinements to the data. Small supercell reflections were observed that could not be indexed using the standard R3m cell for a triangular lattice. To account for this, a refinement was attempted using an R3m cell with a 3 expansion of the a-axis from the corresponding triangular cell, but this did not accurately index the observed supercell reflections. The reflections were properly indexed using the lower symmetry monoclinic space group C2/m (12). In order to determine cation site occupancies, it was assumed that no site mixing occurs between sites in the honeycomb plane (4g and 2a), and the 4g site was fixed at full Ir occupancy. The 2a site was assumed to have mixed Na/Mg occupancy. It was further assumed that Mg only substitutes for Na at the 2a site. The interplane Na sites (2d and 4h) were constrained to have the same occupancy, as were the O1 and O2 sites (8j and 4i). Occupancies were then allowed to refine freely within the constraints. The occupancies determined from refinement to the T = 5 K data were then fixed in refinements to data collected at higher temperatures. Due to the high correlation between isotropic displacement and occupancy, occupancies were determined with fixed U iso = 0.005, then U iso was allowed to refine freely. The results of these refinements are listed in Table 1. Fig. 4 shows NPD datasets collected on Na 2.4 MgIr 2 O 6 at T = 5 K and 100 K. Rietveld refinements to these datasets were carried out in the same way as for the previous sample, however, a small IrO 2 impurity was present in this sample, and was included as a separate phase in the refinements. The results of these refinements are listed in Table 2. Lattice parameters for Na 2.4 MgIr 2 O 6 at room temperature were determined from laboratory XRPD data collected using an air-free sample holder sealed under argon. To confirm the compositions of the samples used for NPD, elemental analyses were performed via ICP-OES. The results of these analyses are shown in Table 3. Initial attempts to dissolve the samples were unsuccessful, and the samples had to be recovered and re-dissolved using microwave digestion. While the Na/Mg ratios agree well with the compositions determined from NPD, the Ir content is anomalously low for both samples, which is likely due to the difficulty associated with fully dissolving Ir and its oxides. The structures of all other samples in the series were characterized via laboratory XRPD. Rietveld refinements were carried out using the reported structure for Na 2 IrO 3 (space group C2/m) [11]. However, because Ir is the dominant x-ray scatterer and stacking disorder is present, the a and b lattice parameters of the monoclinic cell were first obtained by performing a Le Bail fit to the data using space group R3m, where the a-axis is the nearest-neighbor distance in the honeycomb lattice. The resulting a lattice parameter was converted to the the a and b lattice parameters of the monoclinic cell, which were then fixed for the remainder of the refinement. In this way, site mixing between the 4g and 2a sites can be introduced in order to improve relative peak intensities without adversely affecting the lattice parameters. This is particularly important when comparing compounds within a series that may have varying degrees of stacking disorder and site mixing. For reporting purposes, we refer to the resulting site-mixing percentages as the percent of disorder, as they are likely due to stacking faults rather than site mixing. The lattice constants and percents of disorder determined for these compounds are listed in Table 4. 4

5 In order to distinguish between the effects of Mg substitution and Na deficiency, two series of compounds were prepared using the exact same heating schedule. The first series were prepared with the nominal formula Na 3 Na 1 x Mg x Ir 2 O 6, and the second series had the nominal formula Na 3 δ Na 0.75 Mg 0.25 Ir 2 O 6. Plots of lattice parameters versus nominal composition demonstrate the influence of both Mg substitution (Fig. 5(a)) and Na-deficiency (Fig. 5(b)) on the lattice parameters. Substitution of Mg for Na in the honeycomb plane causes the a- and b-axes to contract, and the c-axis to expand slightly. This result supports the assumption that Mg substitutes into the honeycomb plane, as substitution between the planes would cause the c-axis to contract due to electrostatics and the reduced size of Mg relative to Na. Similarly, Na-deficient samples show a contraction of the a- and b-axes and an expansion of the c-axis. Because both Na and Mg contents are variable, compounds in the Na 3 δ (Na 1 x Mg x )Ir 2 O 6 series have a broad range of possible unit cell parameters Magnetic Properties To investigate the influence of lattice geometry on the magnetic properties of the A 2 IrO 3 family, zero-field-cooled (ZFC) magnetization data were collected under an applied magnetic field of µ 0 H = 1 T on each of the compounds in the Na 3 δ (Na 1 x Mg x )Ir 2 O 6 series. Fig. 6 shows plots of inverse magnetic susceptibility (χ 1 H ) versus temperature for all compounds studied in this work. Datasets are staggered along the y-axis for visibility. Linear M fits to the data in the paramagnetic regime (T > 60 K) yield the Curie constant C, the Weiss temperature θ W, and the temperature-independent susceptibility χ 0 (Table 5). θ W, a measure of the average magnetic interaction strength, is negative for all members of the series, indicating that antiferromagnetic interactions dominate in these materials. Fig. 7 shows a plot of θ W versus the average nearest neighbor Ir Ir distance for all members of the Na 3 δ (Na 1 x Mg x )Ir 2 O 6 series. There is a clear trend in the magnitude of θ W as a function of Ir Ir distance, illustrated with a solid red line. Contraction of the honeycomb plane results in a decrease in the magnitude of θ W. The change in interaction strength is likely due to modulation of the Ir O Ir bond angles, which become narrower as the lattice contracts. Contraction of the lattice causes these angles to decrease, resulting in increased ferromagnetic exchange and a concomitant decrease in the magnitude of θ W. To illustrate this effect, the average Ir O Ir bond angles at T = 5 K are given on the plot for the two samples that were analyzed via NPD, along with the average of the two angles reported for Na 2 IrO 3 [11]. The Curie constant C is directly related to the number of Ir 4+ species in the sample, as Ir 3+, Na, and Mg are all diamagnetic. We therefore use both the nominal compositions and the calculated values for C to estimate the ratio of Ir 3+ to Ir 4+, using the compounds whose compositions are known (Na 2 IrO 3, Na 3 Mg 0.5 Ir 2 O 6 and Na 2.3 MgIr 2 O 6 ) as reference points (Table 5). It must be noted that in the strong spin-orbit coupling limit there is a substantial orbital contribution to the observed susceptibility of the Ir 4+ electrons. Nevertheless, these estimates show unambiguously that a significant quantity of low-spin Ir 3+ can exist as non-magnetic holes in A 2 IrO 3 honeycomb lattices. This is unsurprising from a chemical standpoint, as the Ir 3+ /Ir 4+ redox cycle is well documented in the literature, and has even been studied for use in water oxidation catalysis [20]. It is likely that the reduced Ir species are distributed in a disordered manner throughout the lattice, which would explain the origin of the spin glass-like behavior observed in many low quality Na 2 IrO 3 samples. 5

6 Samples in the Na 3 δ (Na 1 x Mg x )Ir 2 O 6 series exhibit spin glass-like behavior. Fig. 8 (a) shows AC magnetic susceptibility data collected on Na 3 Mg 0.5 Ir 2 O 6 and Na 2.3 MgIr 2 O 6. In both sets of data there is a broad peak in the real component of the AC susceptibilityχ at low temperatures, the position T f and height χ max of which vary as a function of the frequency ω of the AC field. Analysis of the data in the context of a canonical spin glass using the Vogel Fulcher law or standard theory for dynamical scaling yield unphysical parameters and therefore little insight [21]. The linear relationship between of ln(ω) and 1 T f indicates an Arrhenius type activation barrier, with activation energies E a between 100 K and 300 K and characteristic frequencies ω 0 on the scale of Hz to Hz [22] (Fig. 8 (b)). While these large values, located in Table 6, are clearly unphysical, they do imply that the frozen, disordered magnetic state observed in these and other iridates is closer to superparamagnetism than to spin glass. This implication is best understood by considering the two major sources of disorder in honeycomb iridates: stacking faults (Fig. 9(a)) and non-magnetic Ir 3+ holes (Fig. 9(b)). Stacking faults occur due to the negligible energy difference between all possible stacking arrangements, but may have little influence on magnetic order because exchange interactions between neighboring layers are small compared to the interactions within each two dimensional layer. Non-magnetic Ir 3+ holes are also likely to exist due to the ease with which Ir 4+ can be reduced. This type of defect, in contrast to stacking faults, can have a significant effect on magnetic order even at low levels, as it gives rise to the possibility of isolated antiferromagnetic domains [6]. Such a picture could explain the apparent proximity to a superparamagnetic state which in fact arises due to a distribution of local AF domain sizes. Fig. 10 (a) shows XRPD data collected over the 2θ range 10 to 50 in six minutes on a nearly pristine sample of Na 2 IrO 3, collected immediately upon removal from an argon-filled glovebox, where the sample was ground and prepared after heating at 900 C. XRPD data collected on powder from the same sample after eight hours of exposure to laboratory air show significant changes in the relative intensities of the Bragg peaks, as well as the development of broad new peaks, however there is only a minor change in the lattice parameters of the Na 2 IrO 3 unit cell over this time period. In contrast, the magnetic properties are significantly affectedbythereaction. Whenhandledproperly, polycrystalline Na 2 IrO 3 exhibits longrange AFM order, as evidenced by the local peak in the magnetic susceptibility data below T N = 15 K (Fig. 10(b)). After eight hours of exposure to air, there is no such local maximum. Curie-Weiss fitting to the linear region of the data above T = 60 K shows that θ W decreases substantially as a result of the reaction (Fig. 10(c)). The Curie constant also decreases significantly, suggesting that the average oxidation state of Ir is changing as a result of the reaction, resulting in fewer S = 1 2 Ir4+ species. The fact that the lattice parameters of the Na 2 IrO 3 cell do not change significantly as a result of exposure to air suggests that this is not simply a result of hydration, as one might expect from experience with structurally analogous systems [23]. The appearance of new peaks suggests that a chemical reaction is occurring between one or more components of laboratory air and Na 2 IrO 3, a possibility that was very recently explored in detail [14]. This result suggests the origin of the discrepancy between magnetic data collected on single crystalline and polycrystalline samples: the small surface area to volume ratio of single crystalline samples hinders the reaction, leaving much of the sample s bulk unreacted. Unfortunately, this implies that previous studies on the 6

7 magnetic properties and structure of single-crystalline Na 2 IrO 3 may have yielded inaccurate data due to sample inhomogeneity. 4. Conclusions AnewseriesofcompoundsbasedontheA 2 IrO 3 prototypewiththeformulana 3 δ (Na 1 x Mg x )Ir 2 O 6 (0 x 1) were synthesized and characterized using NPD, XRPD, ICP-OES, and magnetometry. Substitution of Mg for Na results in contraction of the a- and b-axes and expansion of the c-axis, suggesting that Mg substitution occurs primarily in the honeycomb plane. Similarly, Na deficiency between honeycomb planes contracts the a- and b-axes and expands the c-axis. All samples exhibit stacking disorder, which complicates structural characterization. Magnetic data collected on the compounds studied in this series suggest that there are two main variables that determine the magnetic properties of compounds in the A 2 IrO 3 family: (1) the strength of the magnetic superexchange interactions between neighboring S = 1 2 Ir4+ species is dependent on the angle of the Ir O Ir bonds, for which the in plane lattice parameters a and b are a good structural marker, and (2) the number of S = 1 2 Ir4+ sites is influenced by the chemical composition of the compound because the oxidation state of Ir is variable. This demonstrates that chemical defects can lead to significant numbers of non-magnetic holes, which are likely the root of the pseudo-superparamagnetic behavior observed in many polycrystalline samples. We further showed that Na 2 IrO 3 reacts quickly with laboratory air, producing significant changes in its magnetic behavior, and reported magnetic data on high-quality polycrystalline Na 2 IrO 3, which shows long-range AFM order. Together, our results demonstrate that defects and disorder have significant effects of the magnetism of the A 2 IrO 3 family, and that these materials can be chemically tuned to in order to explore experimentally what may prove to be a rich magnetic phase diagram. 5. Acknowledgements Acknowledgement is made to the donors of the American Chemical Society Petroleum Research Fund and the David and Lucille Packard Foundation for support of this research. References [1] G. Jackeli & G. Khaliullin, Mott Insulators in the Strong Spin-Orbit Coupling Limit: From Heisenberg to a Quantum Compass and Kitaev Models. Physical Review Letters 102, (2009) [2] A. Shitade, H. Katsura, J. Kunes, X.-L. Qi, S.-C. Zhang, N. Nagaosa, Quantum Spin Hall Effect in a Transition Metal Oxide Na 2 IrO 3. Physical Review Letters 102, (2009) [3] H. Gretarsson et al., Crystal-Field Splitting and Correlation Effect on the electronic structure of A 2 IrO 3. Physical Review Letters 110, (2013) [4] J. Chaloupka, G. Jackeli & G. Khaliullin, Kitaev-Heisenberg Model on a Honeycomb Lattice: Possible Exotic Phases in Iridium Oxides A 2 IrO 3. Physical Review Letters 105, (2010) 7

8 [5] I. Kimchi, J. J. Analytis, A. Vishwanath, Three dimensional quantum spin liquid in a hyperhoneycomb iridate model and phase diagram in an infinite D approximation. arxiv: (2013) [6] E. C. Andrade, M. Vojta, Magnetism in spin models for depleted honeycomb lattice iridates: Spin glass order towards percolation. arxiv: (2013) [7] J. G. Rau, E. K.-H. Lee, H.-Y. Kee, Generic spin model for the honeycomb iridates beyond the Kitaev limit. arxiv: (2013) [8] H. Lei, Z. Zhong, H. Hosono, Structural, Magnetic, and Electrical Properties of Li 2 Ir 1 x Ru x O 3. arxiv: (2013) [9] V. M. Katukuri, et al., Kitaev interactions between j = 1/2 moments in honeycomb Na 2 IrO 3 are large and ferromagnetic: Insights from ab initio quantum chemistry calculations. arxiv: (2013) [10] Y.Singh & P. Gegenwart, Antiferromagnetic Mott insulating state in single crystals of the honeycomb lattice material Na 2 IrO 3. Physical Review B 82, (2010) [11] S. K. Choi et al., Spin Waves and Revised Crystal Structure of Honeycomb Iridate Na 2 IrO 2. Physical Review Letters 108, (2012) [12] G. Cao et al., Evolution of Magnetism in Single-Crystal Honeycomb Iridates. arxiv: v1 (2013) [13] K. Baroudi et al., Structure and Properties of α-nafeo 2 -type Ternary Sodium Iridates. arxiv: (2013) [14] J. W. Krizan, J. H. Roudebush, G. M. Fox, R. J. Cava., The Chemical Instability of Na 2 IrO 3 in Air. arxiv: (2013) [15] J. Bréger et al., High-resolution X-ray diffraction, DIFFaX, NMR and first principles study ofdisorder inthe Li 2 MnO 3 Li[Ni 1/2 Mn 1/2 ]O 2 solid solution. Journal of Solid State Chemistry, 178, 2575 (2005) [16] A. C. Larson, R. B. Von Dreele, Los Alamos National Laboratory Report LAUR, (2000) [17] B. H. Toby, EXPGUI, a graphical user interface for GSAS, J. Appl. Crystallogr., 34, (2001) [18] M. M. J. Treacy, J. M. Newsam & M. W. Deem, A general recursion method for calculating diffracted intensities from crystals containing planar faults. Proceedings of the Royal Society, London A433, (1991) [19] E. Climente-Pascual, et. al, Spin 1 2 Delafossite Honeycomb Compound Cu 5SbO 6, Journal of Inorganic Chemistry, 51, 557 (2012) 8

9 [20] J. D. Bakemore, et al, Half-Sandwich Iridium Complexes for Homogeneous Water- Oxidation Catalysis. Journal of the American Chemical Society, 132, (2010) [21] J. A. Mydosh, Spin Glasses: An Experimental Introduction, Taylor and Francis Inc. (1993) ISBN [22] K. Binder, A. P. Young, Spin glasses: Experimental facts, theoretical concepts, and open questions. Rev. Mod. Phys., 58, 801 (1986) [23] J. W. Lynn et al., Structure and dynamics of superconducting Na x CoO 2 hydrate and its unhydrated analog., Physical Review B, 68, (2003) 6. Tables and Figures 9

10 Na 3 Mg 0.5 Ir 2 O 6 (C2/m (12)) 5 K 100 K 295 K a (Å) (3) (3) (5) b (Å) (5) (5) (9) c (Å) (3) (4) (4) β ( ) (7) (8) (1) Na1 (2d) 1 2,0,1 2 occ. 0.81(1) 0.81(1) 0.81(1) U iso 0.020(1) 0.021(1) 0.025(2) 1 Na2 (4h) 2,y,1 occ. 0.81(1) 0.81(1) 0.81(1) 2 y 0.328(2) 0.334(2) 0.326(3) U iso 0.020(1) 0.021(1) 0.025(2) Na3/Mg1 (2a) 0,0,0 occ. 0.54/ / /0.46 U iso (4) (1) (5) Ir1 (4g) 1 2,y,0 y (8) (6) (7) U iso (4) (1) (5) O1 (8j) x,y,z x 0.772(1) 0.754(2) 0.765(2) y (9) (7) (9) z (6) (8) 0.794(1) U iso (3) (6) (9) O2 (4i) x,0,z x 0.717(2) 0.728(2) 0.719(3) z 0.203(1) 0.204(2) 0.201(2) U iso (3) (6) (9) Ir O Ir 1( ) 99.1(4) 97.5(3) 100.7(6) Ir O Ir 2( ) 94.0(3) 95.0(3) 93.8(3) χ R wp R p Table 1: Results of Rietveld refinement to NPD data collected on Na 3 Mg 0.5 Ir 2 O 6 (sample 7). 10

11 Na 2.4 MgIr 2 O 6 (C2/m (12)) 5 K 100 K a (Å) (6) (5) b (Å) 9.197(1) (9) c (Å) (3) (3) β ( ) (7) (6) Na1 (2d) 1 2,0,1 2 occ. 0.77(1) 0.77(1) U iso 0.025(1) 0.041(9) 1 Na2 (4h) 2,y,1 occ. 0.77(1) 0.77(1) 2 y 0.331(3) 0.343(2) U iso 0.025(1) 0.041(9) Na3/Mg1 (2a) 0,0,0 occ. 0/1 0/1 U iso (5) 0.041(5) Ir1 (4g) 1 2,y,0 y (9) 0.171(1) U iso (5) 0.041(5) O1 (8j) x,y,z x 0.763(2) 0.766(2) y 0.162(1) 0.164(1) z 0.802(1) 0.798(1) U iso (4) 0.012(4) O2 (4i) x,0,z x 0.724(3) 0.726(2) z 0.198(2) 0.189(2) U iso (4) 0.012(4) Ir O Ir 1 ( ) 97.8(6) 99.5(5) Ir O Ir 2 ( ) 94.4(4) 93.9(3) χ R wp R p Table 2: Results of Rietveld refinement to NPD data collected on Na 2.4 MgIr 2 O 6 (sample 11). Sample No. NPD Formula Na : Mg : Ir molar ratio 7 Na 2.5 Na 0.5 Mg 0.5 Ir 2 O 6 ICP-OES 6.52 : 1.00 : 3.23 NPD 7.00 : 1.00 : Na 2.4 MgIr 2 O 6 ICP-OES 2.30 : 1.00 : 1.49 NPD 2.40 : 1.00 : 2.00 Table 3: results of ICP-OES elemental analysis performed on Na 3 Mg 0.5 Ir 2 O 6 and Na 2.4 MgIr 2 IrO 6 compared with compositions determined from Rietveld refinement to NPD data. Ir contents determined from ICP-OES are anomalously low due to complications which arose in the dissolution process. 11

12 Sample No. Target Stoichiometry a (Å) b (Å) c (Å) β ( ) % Disorder 1 Na 3 NaIr 2 O (1) (1) (2) (4) 30% 2 Na 3 NaIr 2 O 6 (8h in air) (1) (1) (2) (5) 20% 3 Na 3 Na 0.9 Mg 0.1 Ir 2 O (1) (1) (3) (5) 38% 4 Na 3 Na 0.8 Mg 0.2 Ir 2 O (3) (3) (3) (6) 50% 5 Na 3 Na 0.7 Mg 0.3 Ir 2 O (2) (2) (3) (6) 37% 6 Na 3 Na 0.5 Mg 0.5 Ir2O (1) (1) (3) (5) 32% 7 Na 2.5 Na 0.5 Mg 0.5 Ir 2 O 6 * (5) (9) (4) (1) 8 Na 3 Na 0.52 Mg 0.48 Ir 2 O (2) (2) (3) (5) 38% 9 Na 2.56 Na 0.5 Mg 0.5 Ir 2 O (3) (3) (1) (2) 30% 10 Na 2.46 Na 0.24 Mg 0.76 Ir 2 O (1) (1) (2) (4) 28% 11 Na 2.4 MgIr 2 O 6 * (1) (1) (3) (6) 32% Table 4: Target stoichiometries and unit cell parameters determined from XRPD data collected at room temperature of all compounds studied in this work. Compositions marked with an asterisk were determined by Rietveld refinement to NPD data collected on the final product and verified via ICP-OES. The 4g/2a site mixing percentages, referred to as % Disorder, are a rough measure of the amount of structural disorder present in the sample. emu Sample No. Target Stoichiometry χ 0 ( molir K molir Oe W (K) Ir 3+ :Ir 4+ (±0.1) 1 Na 3 NaIr 2 O 6 1.1(6) (3) (9) 0:1 2 Na 3 NaIr 2 O 6 (8h in air) -5(2) (8) -146(2) 0.3:0.7 3 Na 3 Na 0.9 Mg 0.1 Ir 2 O 6-9(1) (7) -87(1) 4 Na 3 Na 0.8 Mg 0.2 Ir 2 O 6-9(1) (4) -63(2) 5 Na 3.06 Na 0.7 Mg 0.3 Ir 2 O 6-6(1) (1) -48(3) 6 Na 3.26 Na 0.5 Mg 0.5 Ir2O 6 3(1) (5) -35(2) 0.4:0.6 7 Na 2.5 Na 0.5 Mg 0.5 Ir 2 O 6 * 3.7(8) (8) -30(1) 0.1:0.9 8 Na 3 Na 0.52 Mg 0.48 Ir 2 O 6-3.4(6) (6) -28(2) 0.3:0.7 9 Na 2.56 Na 0.5 Mg 0.5 Ir 2 O 6-9(1) (4) -16(2) 0.1: Na 2.46 Na 0.24 Mg 0.76 Ir 2 O 6-1.2(9) (9) -11(1) 11 Na 2.4 MgIr 2 O 6 * -1.1(7) (5) -1(1) 0.2:0.8 Table 5: Curie-Weiss parameters determined from linear least-squares fitting to inverse magnetic susceptibility data in the paramagnetic regime. Sample No. Target Stoichiometry T f (100 Hz) E a (K) ω 0 (Hz) 5 Na 3 Na 0.7 Mg 0.3 Ir 2 O (9) Na 3 Na 0.5 Mg 0.5 Ir2O (9) Na 2.5 Na 0.5 Mg 0.5 Ir 2 O 6 * (13) Na 3 Na 0.52 Mg 0.48 Ir 2 O (20) Na 2.56 Na 0.5 Mg 0.5 Ir 2 O (14) Na 2.4 MgIr 2 O 6 * (11) Table 6: Freezing temperatures T f, activation energies E a, and characteristic frequencies ω 0 determined from an Arrhenius analysis of AC magnetic susceptibility data collected on members of the Na 3 δ (Na 1 x Mg x )Ir 2 O 6 series. 12

13 Figure 1: (a) A simple triangular lattice can be described by an R-3m unit cell (green) with in-plane lattice parameter a. Formation of honeycomb order by substitution of one third of the atoms in the triangular lattice is accommodated by a 3 expansion of the triangular lattice constant yielding a larger R-3m unit cell (red). Different stacking arrangements(not shown) can require a lower symmetry unit cell, such as the C2/m cell shown (blue). (b) The structure of Na 3 δ (Na 1 x Mg x )Ir 2 O 6 is described by the C2/m cell shown in the ab-plane. Oxygen atoms above the honeycomb plane are bonded to Ir with thick black lines, while those below the plane are bonded to Ir with thin lines. IrO 6 octahedra share edges to form a honeycomb lattice, highlighted by a light blue hexagon. Two distinct Ir O Ir bond angles, labeled 1 and 2, are possible in this structure. The structure is also shown in the ac-plane to highlight the stacking arrangement corresponding to the C2/m cell. Figure 2: (a) XRPD patterns simulated using DIFFaX are shown for layered structures of the following types: triangular (green), ordered honeycomb (blue), and ordered honeycomb with 5% stacking faults (red). The influence of honeycomb ordering on the XRPD pattern is most apparent in the shaded region between 2θ = 18 and 22. (b) XRPD data collected on three different samples of Na 2 IrO 3 with varying degrees of structural disorder. The top sample has the fewest stacking faults, as illustrated by the sharp peaks in the shaded region between 2θ = 18 and

14 Figure 3: NPD data (black circles) collected on Na 3 Mg 0.5 Ir 2 O 6 at T = 5 K, 100 K, and 295 K are shown along with Rietveld refinements (red) and the difference between the data and the fit (green). Tick marks (gray) indicate the positions of expected Bragg reflections. Figure 4: NPD data (black circles) collected on Na 2.4 MgIr 2 IrO 6 at T = 5 K and 100 K are shown along with Rietveld refinements (red) and the difference between the data and the fit (green). A small ( 10%) IrO 2 impurity was present in the sample, and was included as a separate phase in the refinements. Tick marks (gray) indicate the positions of expected Bragg reflections. Tick marks corresponding to the IrO 2 impurity peak positions are located below the tick marks for the structure. 14

15 Figure 5: Plots of lattice parameters vs. nominal chemical composition for the cases of varied Mg substitution (Na 3 (Na 1 x Mg x )Ir 2 O 6 ) (a) and varied Na-content (Na 3 δ (Na 0.75 Mg 0.25 )Ir 2 O 6 ) (b). Substitution of Mg for Na in the honeycomb plane causes the a- and b-axes(black diamonds and red circles, respectively) to contract and the c-axis (blue squares) to expand. Decreasing the amount of sodium in the lattice has a similar effect. 15

16 Figure 6: Inverse magnetic susceptibility is plotted against temperature for each member of the series. Datasets are shifted along the y-axis for visibility. Linear fits to inverse magnetic susceptibility data over the range T = 60 K to 300 K yield the Curie constant C and the Weiss temperature θ W for each member of the series. Datasets are shifted along the y-axis for visibility. The Weiss temperatures θ W for each member of the series are indicated by red circles, corresponding to the x-intercept of the linear extrapolation from the Curie-Weiss regime. 16

17 Figure 7: Plot of the Weiss temperature θ W versus the average nearest neighbor Ir Ir distance for all compoundsstudiedinthis work. Negativevaluesofθ W indicatethat antiferromagneticinteractionsdominate the magnetic susceptibility of these layered honeycomb iridates, and the magnitude of θ W is a measure of the strength of these interactions. Expansion of the honeycomb plane increases the overall antiferromagnetic interaction strength as Ir O Ir bond angles deviate from 90. Average Ir O Ir bond angles, determined from analysis of NPD data, are shown for two points on the plot, as well as the reported value for Na 2 IrO 3, denoted by red boxes [11]. θ W is largest in Na 2 IrO 3 and is diminished by chemical substitution and sodium vacancies, as well as decomposition in air. The smooth curve through the data is a guide to the eye. 17

18 Figure 8: (a) AC susceptibility data collected on two representative samples from the Na 3 δ Na 1 x Mg x Ir 2 O 6 series under an applied field of H = 796 A/m (10 Oe) with amplitude of H = 398 A/m (5 Oe). Both samples show a peak in the real component χ of the AC susceptibility at the freezing temperature T f, and the magnitude of the peak decreases and shifts to higher temperatures as the frequency of the AC field is increased, consistent with a spin glass-like transition (inset). (b) A plot of ln(ω) versus 1 T f yields a roughly linear dependence, indicating proximity to a super paramagnetic regime, and the slope and intercept of the linear fit yield activation energy E a and the characteristic frequency ω 0 for spin reorientation. 18

19 Figure 9: Common types of disorder in A 2 IrO 3 honeycomb iridates. (a) In a fully ordered sample, honeycomb layers are stacked in a perfectly repeating pattern (ABCABCABC, for example). Stacking faults (ABCCACBCA, here) complicate structural characterization as they masquerade as (Na/Mg) mixing onto the Ir site. (b) Non-magnetic Ir 3+ also exist due to chemical defects. These holes perturb the magnetic order in A 2 IrO 3 materials, and can lead to the disordered freezing of spins commonly observed in polycrystalline samples. We speculate that this type of disorder leads to the formation of isolated islands of AFM order, which vary in size and interact weakly with one another, which is one possible explanation for the pseudo superparamagnetic behavior observed in AC magnetic susceptibility measurements. 19

20 Figure 10: (a) XRPD data collected over the 2θ range 10 to 50 in six minutes on Na 2 IrO 3 immediately after removal from an argon-filled glovebox (blue, bottom) and after eight hours of exposure to laboratory air (red, top). A reaction occurs between Na 2 IrO 3 and one or more components of the air that causes the relative intensities of the C2/m reflections to change, and new reflections emerge, highlighted with gray arrows. Small circles indicate reflections due to crystalline Silicon, which was used as an internal standard for the purpose of Rietveld refinement. (b) The magnetic susceptibilities of polycrystalline Na 2 IrO 3 (blue circles) and the same powder after eight hours of exposure to air at room temperature (red diamonds) are compared. When handled in air-free conditions, Na 2 IrO 3 exhibits long-range AFM order, as evidenced by the local maximum in the magnetic susceptibility at T N = 15 K. (c) A plot of inverse magnetic susceptibility versus temperature, along with linear fits to both datasets above T = 60 K yield the Weiss temperature θ W and Curie constant C for Na 2 IrO 3 before (blue circles) and after (red diamonds) exposure to air. 20