Solution of the heavily stacking faulted crystal structure of the honeycomb iridate H3LiIr2O6

A powder sample of pure H3LiIr2O6 was synthesized from α-Li2IrO3 powder by a soft chemical replacement of Li+ with H+. The crystal structure of H3LiIr2O6 consists of sheets of edge sharing LiO6- and IrO6-octahedra forming a honeycomb network with layers stacked in a monoclinic distorted HCrO2 type pattern. Heavy stacking faulting of the sheets is indicated by anisotropic peak broadening in the X-ray powder diffraction (XRPD) pattern. The ideal, faultless crystal structure was obtained by a Rietveld refinement of the laboratory XRPD pattern while using the LiIr2O63--layers of α-Li2IrO3 as a starting model. The low radial distances of the PDF function, derived from synchrotron XRPD data, as constraints to stabilize the structural refinement. DIFFaX-simulations, structural considerations, high radial distances of the PDF function and a Rietveld compatible global optimization of a supercell were employed to derive a suitable faulting model and to refine the microstructure using the experimental data. We assumed that the overall stacking pattern of the layers in the structure of H3LiIr2O6 is governed by interlayer O-HO contacts. From the constitution of the layers, different stacking patterns with similar amounts of strong O-HO contacts are considered. Random transitions among these stacking patterns can occur as faults in the crystal structure of H3LiIr2O6, which quantitatively describe the observed XRPD.


Introduction
H 3 LiIr 2 O 6 is a layered double hydroxide oxide salt, hence the formula can be expressed alternatively as LiIr 2 O 3 (OH) 3 .The compound was recently revealed to be the pH sensing phase in an IrO x -based solid state electrode. 13][4] IrO x -based solid state electrodes can be produced by the carbonate melt oxidation process [5][6][7] in which iridium wires are oxidized in an excess of molten lithium carbonate.This process leads to the formation of α-Li 2 IrO 3 . 1,8In a subsequent step the excess of lithium carbonate is dissolved by acid treatment and α-Li 2 IrO 3 is transformed into H 3 LiIr 2 O 6 by cation exchange (eqn (1)).
Profound crystallographic knowledge of the α-Li 2 IrO 3 -H 3 LiIr 2 O 6 system is essential for understanding and controlling the process presented in (eqn (1)) and in consequence also for the production of stable and reliable pH sensing IrO xbased solid state electrodes.
The crystal structure of α-Li 2 IrO 3 was determined from X-ray powder diffraction (XRPD) data. 9It consists of layers of edge sharing IrO 6 -and LiO 6 -octahedra and lithium ions that are situated in the interlayer spacing (Fig. 1, left).Accordingly the formula can be alternatively written as Li I 3 Li II Ir II 2 O 6 with atoms situated in the interlayer space indicated be "I" and atoms within the layers indicated by "II".O'Malley et al. 9 observed the occurrence of stacking faults in the crystal structure of α-Li 2 IrO 3 , which is a common phenomenon in layered alkali-and earth alkali metal iridates. 10,11Hence the slight occupational disorder between Li + and Ir 4+ in the layers, which was introduced to the structural model during the refinement, could be an artefact from daubed reflection intensities due to diffraction line broadening by the appearance of stacking faults. 10,11The overall structural motif of α-Li 2 IrO 3 is closely related, but not identical to LiCoO 2 . 12Due to the honeycomb ordering of the 2 Ir-sites and the Li-site within the layer of Li 2 IrO 3 and the associated monoclinic distortion of the unit cell of α-Li 2 IrO 3 the anion sublattice exhibits an ABCA′B′C′A″B″ C″…, stacking (Fig. 1, left), with anion positions indicated by capital Latin letters, that slight differs from the ABC anion stacking pattern in the structure of LiCoO 2 .
By acid treatment only the lithium ions situated in-between the layers are exchanged by protons which lead to the formation of H 3 LiIr 2 O 6 (Fig. 1, right).The occurrence of stacking faults in α-Li 2 IrO 3 is crucial for the production of H 3 LiIr 2 O 6 , as the degree of cation exchange is the higher, the higher the degree of faulting of the precursor material. 1O'Malley et al. 1 expected H 3 LiIr 2 O 6 to crystallize in a HCrO 2 like structure with an AABBCC stacking pattern of the anion sublattice, where oxygen sites of neighbouring layers directly oppose each other, which results in strong H-bonds.Because of the pronounced occurrence of stacking faults they only determined a rough estimate of the crystal structure.Neither the exact type nor the degree of stacking faulting in the H 3 LiIr 2 O 6 structure, which has an impact of the pH sensing properties, 1 has been understood, yet.
H 3 LiIr 2 O 6 attracts considerable attention not only as PH sensor materials but as a promising candidate for a topological quantum spin liquid described by Kitaev model. 13In particular iridates with a honeycomb like structural motif, like α-Li 2 IrO 3 9 and β-Li 2 IrO 3 14 are promising candidates for the materialization of the Kitaev model.In a honeycomb lattice strong spin frustration effects stabilize a quantum spin liquid state, in which novel excitations, like Majorana fermions and fluxes are apparent. 15α-Na 2 IrO 3 16 and α-Li 2 IrO 3 [17][18][19][20][21] were first expected to be a materialization of Kitaev model but a magnetically ordered state was found to be the ground state of the two honeycomb iridium oxides.H 3 LiIr 2 O 6 , a modified honeycomb iridate, was visited as an alternative candidate and was discovered not to show any magnetic ordering down to 1 K. 22 This is indicative for the realization of a spin liquid as the ground state.To clarify the reason why quantum spin liquid can be stabilized in H 3 LiIr 2 O 6 , a detailed structural understanding, including the microstructural effects like stacking faults, is necessary.
In the current study, we re-evaluate the structure model of H 3 LiIr 2 O 6 given by O'Malley et al. 1 and after we show that this model is not suitable, we describe the redetermination of the ideal crystal structure of H 3 LiIr 2 O 6 by using XRPD data and PDF-analysis.As the XRPD data indicate a strongly stacking faulted crystal lattice, possible structural defects are derived from the crystal structure and investigated by systematic DIFFaX simulations.A Rietveld compatible approach was used to determine the real structure, i.e. the degree of faulting by using a supercell approach.The obtained superstructure is confirmed by PDF-analysis.The approach used in this study should pave a way to determine the crystal structure and the degree of stacking faulting of the layered materials of interest.

Experimental section
Synthesis of H 3 LiIr 2 O 6 H 3 LiIr 2 O 6 was synthesized by using an α-Li 2 IrO 3 precursor.The precursor was obtained by a solid state reaction of Li 2 CO 3 (Chempur, 99.995+%) and IrO 2 (Tanaka Kikinzoku Kogyo, 86 wt% Ir) in which a 5% excess of Li 2 CO 3 was added to the stoichiometric mixture in order to compensate for volatilization of Li at high temperatures.Tempering of the mixture was conducted at 1000 °C for 24 h.During that period the powder was ground twice to provide a sufficient mixing.Afterwards the mixture was cooled down to room temperature in the oven.The obtained powder was confirmed to be a single phase of α-Li 2 IrO 3 by laboratory XRPD.For cation exchange the α-Li 2 IrO 3 powder was into in a Teflon-lined steel autoclave and 4 mol L −1 H 2 SO 4 aqueous solution was added.The mixture was heated in the sealed vessel at 120 °C for 96 hours.Afterwards the product was washed with distilled water and dried at room temperature.The obtained powder is insensitive against air.

Solid phase characterization
Laboratory X-ray powder diffraction (XRPD).The XRPDpattern of the solid phase was collected at room temperature on a laboratory powder diffractometer in Debye-Scherrer geometry (Stadi P-Diffractometer (Stoe), Ag-K α1 radiation from primary Ge(111)-Johann-type monochromator, Mythen 1 K detector (Dectris)).The sample was sealed in a 0.5 mm diameter borosilicate glass capillary (Hilgenberg glass no.14), which was spun during the measurement.
Synchrotron X-ray powder diffraction.Synchrotron X-ray total scattering pattern of the sample was recorded in Debye-Scherrer mode at the National Synchrotron Light Source (NSLS), Brookhaven National Laboratory, at a wavelength of λ = 0.1847 Å (67.13 keV) on beamline X28B1 equipped with a Si(311) sagittal focusing double Laue crystal monochromator.Diffracted X-rays were detected using a PerkinElmer detector (2048 × 2048 pixels and 200 × 200 μm 2 pixel size) mounted orthogonal to the beam path with a sample-to-detector distance of 205.80 mm.The measured samples were sealed in kapton tubes of 1.0 mm diameter.The collected Debye-Scherrer rings were subsequently azimuthally integrated with the program FIT2D 23 to one-dimensional powder diffraction patterns in Q [nm −1 ] and 2θ [°] versus intensity.Parameters for intensity integration were determined from a LaB 6 reference sample.Using the software TOPAS, 24 Rietveld refinements 25 were carried out, applying the fundamental parameter method 26,27 for peak modelling and Chebyshev polynomials for modelling the background.
Pair distribution function (PDF) analyses.PDFGETX2 28 was used to correct and normalize the diffraction data and the Fourier transform in order to obtain the experimental PDF, G(r).The Q range used in the Fourier transform was limited to 30 Å −1 .Nickel was previously measured as a standard material to determine the Q-damp and Q-broad parameters which are the parameters that correct the PDF envelope function for instrumental resolution effects. 29,30For structural modelling TOPAS 6.0 31 was used.Atomic positions of the compounds under investigation were taken from the corresponding structural models presented here and from the assumed crystal structure of H 3 LiIr 2 O 6 of O'Malley et al. 1 However, certain parameters were adjusted to improve the agreement.These were scale factor, Q-damp, Q-broad, the two trigonal lattice parameters and an overall isotropic displacement parameter.

Results and discussion
Crystal structure solution Evaluation of the HCrO 2 -like structure model.The measured XRPD pattern of the synthesized H 3 LiIr 2 O 6 appeared to be very similar to the reported pattern of O'Malley et al. 1 Accordingly, in a first attempt the trigonal HCrO 2 like crystal structure model of H 3 LiIr 2 O 6 which was created by them was used as a starting model for a Rietveld refinement 25 of the measured diffraction pattern using the program TOPAS 6.0. 31This procedure, however, led to an insufficient fit as several sharp peaks couldn't be described by the crystal structure model (Fig. 2, green circles).The misfit of the reflection intensities is most likely caused by the insufficient modelling of peak broadening, which is caused by the pronounced occurrence of stacking faults.This indicates that the R3 ˉm1 space group symmetry, as well as the unit cell of the predicted crystal structure is not suitable.
In addition, in the trigonal HCrO 2 like structure model the cation sublattice is completely occupationally disordered between lithium and iridium and all oxygen sites are situated on an identical plane.If this model was appropriate, the cation exchange by acid treatment would have caused a vast change in the constitution of the layers.The solid-state NMR data, however, point to an ordered cation sublattice (ESI, Fig. S1 †).
Derivation of the idealized, faultless crystal structure of H 3 LiIr 2 O 6 .As H 3 LiIr 2 O 6 was produced by cation exchange via acid treatment of α-Li 2 IrO 3 (eqn (1)), it was assumed that the crystal structure of the resulting solid phase shows some similarities to the one of the starting material, especially in terms of the constitution of the layers.By comparing a simulated XRPD-pattern of α-Li 2 IrO 3 with the measured pattern H 3 LiIr 2 O 6 , the latter seems to exhibit fewer diffraction lines (Fig. 3).A close inspection of the line profile, however, exhibits characteristic triangular, Warren-type 32   starting point for indexing the XRPD pattern of H 3 LiIr 2 O 6 .
Accordingly the intense basal reflection at ≈7°2θ was indexed as 001.The upshift of this reflection from α-Li 2 IrO 3 to H 3 LiIr 2 O 6 (Fig. 3) indicates that the interlayer distance ( product, c•sin(β − 0.5•π), of the lattice parameters c and β) of the latter phase is significantly smaller than in the starting material.In contrast some reflections like 130, 1 ˉ31, 131 and 3 ˉ3 ˉ1 show a downshift, which points to an elongation of the aand/or b-axis after the cation exchange.The 020 reflection seems to be missing in the diffraction pattern of H 3 LiIr 2 O 6 .A close inspection of the 001 reflection (Fig. 4), however, reveals an unusual peak shape.This reflection exhibits small shoulders both on the high and low angle side.The shoulder on the high angle side can be explained by an overlap with the Warren type peak, which, however, doesn't explain the shoulder on the low angle side.This shoulder also can't be explained by the instrumental profile, as the 130 and 1 ˉ31 reflections do not exhibit such a peak shape (Fig. 4).Accordingly the unusual profile of the 001 reflection must be caused by an overlap with another peak that is situated at a slightly lower diffraction angle.This was confirmed by a single line fit using two peaks.The additional reflection was indexed as 020 and its position pointed to a length of the b-axis of ≈9.24 Å.The information derived from the XRPD pattern was used for global optimization of the lattice parameters of H 3 LiIr 2 O 6 by a series of Pawley fits 34 employing the LP-search routine implemented into TOPAS. 35he resulting global minimum was used as a starting model for the determination of precise lattice parameters by using a LeBail 36 fit and applying the fundamental parameter approach of TOPAS. 27or the determination of the layer constitution of H 3 LiIr 2 O 6 the atomic positions of the LiIr 2 O 6 3− -layers of α-Li 2 IrO 3 were used as starting values.The inappropriate 10,11 occupational disorder between Li and Ir on the metal site (2) in the structural model of O'Malley et al. 9 was removed, as this is con-sidered as an artefact of the refinement, caused by the occurrence of stacking faults.In order to avoid a correlation between the diffractional effects caused by the stacking faults and the refinement of the layer constitution, regions in the XRPD pattern that are affected by peak broadening (7.18-11.88°2θ,13.20-13.76°2θ,14.30-15.00°2θ,15.46-16.00°2θ, 17.24-20.25°2θ,21.44-21.72°2θand 22.76-24.00°2θ)were excluded.Due to the exclusion of various regions of the XRPD pattern, an unconstrained refinement of the atomic coordinates and the lattice parameters was not possible, as the result strongly depends on the sequence in which the parameters were released to refinement.In addition, some resulting Ir-O and Li-O distances were either unreasonable long (>2.20 Å) or short (<1.80 Å).In order to introduce reasonable con-and restraints, the low distance region (<2.5 Å) of the measured PDF-curve was used, which refers to intralayer cation-oxygen distances (Fig. 5).The measured PDF-curve exhibits some noise and artefacts appearing as sine-type modulation in the region of low distances.Despite this noise a well-defined peak is clearly visible at 2.01 Å (Fig. 5, red font colour).This peak is mainly attributed to Ir-O pairs as more iridium than lithium is present in the solid and as iridium is a much stronger scatterer than lithium.According to the peak width and shape (small shoulder at higher distances) the distribution of the metal-oxygen distances is not necessarily unimodal and those distances can range from ≈1.90 Å to ≈2.10 Å.In consequence artificial penalty function depending on the metal-oxygen distances with minimums at 2.0 Å both for the Ir-O and Li-O distances were included into the Rietveld refinements.This led to robust and reproducible results.During the refinement all atomic and lattice parameters were released iteratively and could be refined without using further con-or restraints.The crystallographic data and the refined atomic coordinates of H 3 LiIr 2 O 6 are given in Tables 1 and  2, the graphical result of the refinement and the agreement factors are presented in Fig. 6.
The obtained structural model of the LiIr 2 O 6 3− -layers of H 3 LiIr 2 O 6 can be evaluated as a reasonable approximation (see

Dalton Transactions Paper
This journal is © The Royal Society of Chemistry 2017 A comparison of the crystal structure of the starting material and the product after the cation exchange by acid treatment is presented in Table 1 and Fig. 7.For this comparison it was assumed that α-Li 2 IrO 3 does not exhibit any occupational disorder between Li and Ir in the LiIr 2 O 6 3− -layers.The layers of H 3 LiIr 2 O 6 consist of edge linked, distorted LiO 6 -and IrO 6 -octahedra that form sheets within the ab-plane of the unit cell (Fig. 7a, on the right).Within these layers each LiO 6/3 3− -octahedron (Fig. 7, green) is surrounded by six IrO 6/3 -octahedra (yellow) and each oxygen atom is coordinated by two iridium ions and one lithium ion.Hence, the cation exchange did not produce significant changes in the cation sublattice (Fig. 7a).
Due to the exchange of Li + by less space filling H + , the interlayer distance was shortened from 4.82 Å in α-Li 2 IrO 3 to 4.54 Å in H 3 LiIr 2 O 6 (Fig. 7c) and the arrangement of the layers, i.e. the stacking order, was altered.The LiIr 2 O 6 3− -layers in α-Li 2 IrO 3 are stacked in a monoclinic distorted CdCl 2 (C19type) like fashion with an ABCA′B′C′A″B″C″…, stacking pattern of the anions, whereas the layers in H 3 LiIr 2 O 6 are stacked in a monoclinic distorted HCrO 2 (3R-type) like way with an ABBCCA′A′B′B′C′C′… stacking pattern of the anions (Fig. 7b  and d).In addition, slight changes in the anion sublattice are obvious.By the partial protonation of oxide ions during the cation exchange, hydroxide ions are formed.As the stacking    order of H 3 LiIr 2 O indicated that anions of neighboring layers directly oppose each other, attractive O⋯H-O interactions, most likely hydrogen bonds, between anion sites of neighboring layers can be assumed.In consequence, the curly character of the sheets, i.e. the distance between the planes in which the two oxygen sites are situated, is much less pronounced in H 3 LiIr 2 O 6 than in α-Li 2 IrO 3 (Fig. 7c and d).
A comparison of the distorted LiO 6/3 3− -and IrO 6/3 -octahedra in α-Li 2 IrO 3 and H 3 LiIr 2 O 6 is given in Fig. 8 and the metal-oxygen distances are listed in Table 3.In general the coordination polyhedra are only little affected by the cation exchange, only some bond distances in the LiO 6 -octahedra in H 3 LiIr 2 O 6 are slightly elongated compared to α-Li 2 IrO 3 , whereas the IrO 6 -octahedra exhibits almost no change.All metal-oxygen distances are in reasonable range.The distances between the cations (Table 3, bottom) are slightly elongated.The obtained values (3.06 Å-3.10 Å), however, are in very good agreement with the PDF-data (Fig. 5).As no penalty function was applied on the distances between the cations, these results indicate that a suitable model of the layer constitution was derived by Rietveld refinement.
The cation sublattice of H 3 LiIr 2 O 6 reveals a pronounced pseudo symmetry.A pseudo trigonal lattice as described by O'Malley et al. 1 can be found with a ≈3.1 Å (ESI, Fig. S2, † grey lines), which symmetry is broken by the ordering between Li and Ir.In addition the Ir 6 O 18 12− honeycombs (Fig. 7b, blue octahedra) exhibit a very pronounced pseudo hexagonal symmetry (ESI, Fig. S2, † green lines) that is only broken by the stacking order of the layers (see next sections).Derivation of the stacking pattern and of potential structural defects.Recently the real structures of related stacking faulted iridates with honeycomb like lattices were successfully described by super positioning of different stacking types. 10,11ased on these approaches possible stacking orders of the layers are derived from structural considerations.
The arrangement of the layers, i.e. the stacking order, in H 3 LiIr 2 O 6 can be described by a stacking vector S1 that is identical with the c-axis of the monoclinic unit cell.The stacking pattern of the layers results almost in a CrOOH like stacking, (AγB)(BαC)(CβA)(AγB) i (BαC) i (CβA) i (AγB) i+1 (BαC) i+1 …, with anion positions indicated by capital Latin letters, anion positions indicated by small Greek letters and layers indicated by parenthesis (Fig. 9a).The stacking order exhibits a monoclinic distortion, i.e. layer i + 3 is not identical with layer i, which is indicated by superscripted I and I + 1. Anions of neighboring layers are situated almost in direct opposition to each other According to the limits of the XRPD-method no attempt was made to determine and refine the hydrogen positions of the crystal structure.Due to the low interlayer spacing and the short O-O distances (Fig. 9a 3.    10b).Hence hydrogen bonds between these sites can only be formed if the bond between oxygen and hydrogen in a hydroxide ion is canted.This S2-stacking pattern with 4 strong hydrogen bonds and 2 bonds with potentially canted hydroxide ions should be energetically a little less favored with respect to the S1-stacking pattern with 6 strong and hydrogen bonds (Fig. 10a).Nevertheless, within the S1-stacking pattern of LiIr 2 O 3 (OH) 3 layers stacked in an alternative, S2-like way, can appear as stacking faults.
Alternative stacking vectors can be derived directly from the constitution of the layers.Due to the sharp 00l reflection in the diffraction pattern of H 3 LiIr 2 O 6 (Fig. 12, blue line) the interlayer spacing is not affected by stacking faults, i.e. the z-component of each stacking vector must be 1.0.The stacking vector S1 (Fig. 11a, magenta) directs an oxygen site, e.g.O2, in direct opposition to an identical site (red ball) of the preceding layer.Alternative stacking vectors bring different oxygen atoms in the surrounding of this O2-site in direct opposition to the O2 site of the preceding layer.Thus the x-and y-components of alternative stacking vectors can be derived from a projection of an oxygen layer onto the ab-plane (Fig. 11b).Therefore, the position of one O2 site is used as the origin for fractional xand y-coordinates using a-and b-lattice parameters.There are four O2 sites in direct surrounding of an O2 site, each of these positions can be reached by a shift, denoted as S1-1, of x = ±1/2 and y = ±1/2.Due to layer symmetry, with symmetry centers at (0,0); (±1/4, ±1/4) and (±1/2, ±1/2), each of these shifts is symmetry equivalent to S1 and therefore transitions among S1 and S1-1 stacking won't produce any additional diffraction effects.Six O1 sites are located in direct surrounding of one O2 site, as well.The shift from O2 to any of these sites produces a S2 like stacking described above.Accordingly, six alternative stacking vectors (S2-1 to S2-6) can be derived from the relative positions of O1 sites surrounding a central O2 site (Table 4).
In a hexagonal lattice all stacking vectors would be equal.The symmetry is broken, as one stacking vector, S1, seems to be preferred.

Systematic DIFFaX-studies
The DIFFaX-routine 37 was used to simulated the diffraction patterns of stacking faulted LiIr 2 O 3 (OH) 3 , containing transitions among S1 and S2-i-stacking.The simulations were carried out in recursive mode, i.e. an average diffraction pattern of all possible stacking sequences, weighted by transition probabilities is calculated in each simulation.The number of layers and the   layer extension in a-and b-direction was to infinite.A Pseudo-Voigt function, derived from a LeBail-fit 36 of the measured diffraction pattern without using crystalline size parameters, was used to describe the peak profile.
Transformation of the unit cell and the atomic coordinates.As only unit cells with α, β = 90°and γ = variable can be used in the DIFFaX-routine, a transformation of the monoclinic unit cell of H 3 LiIr 2 O 6 into an orthorhombic one with α, β, γ = 90°, was necessary.Therefore, a new, β-dependent, c′-axis had to be defined (Table 5).In addition all atomic coordinates (x, y, z), as well as the x-, y-and z-components of the stacking vectors (Table 4) were adopted to the transformed c′ and β′ lattice parameters using (eqn (1)).Due to the translational symmetry, the simulations can be limited to 4 stacking vectors: S1, S2-1, S2-2 and S2-3.A 4 × 4 matrix of transition probabilities, P ij , (Table 6) was used to describe different degrees of faulting in the crystal structure of H 3 LiIr 2 O 6 .
Diffraction effects of stacking faults within the structure of LiIr 2 O 3 (OH) 3 .As a starting point the diffraction pattern of faultless H 3 LiIr 2 O 6 using pure S1-stacking with P 1j = 1.0 was simulated (Fig. 12 The observed peak broadening and triangular peak shape in combination with structural considerations (Fig. 11) were used to introduce constraints for the transitions probabilities, P ij , to simplify the 4 × 4 transition probability matrix for systematic DIFFaX-simulations.As each alternative stacking vector, S2-i leads to the same number of strong O-H-O contacts, there shouldn't be any preference for any of these vectors after a shift from S1 stacking.Accordingly, the probability for a shift from S1-type stacking to S2-type stacking must be equal for all S2-i stacking vectors.In consequence the probability of a fault in an S1 stacking pattern can be described with the parameter x (Table 7).As each line of the transition probability matrix has to sum up to 1.0, each transition probability from S1 to S2-i can be expressed as x/3.With respect to the constitution of the layers, a shift from S1-stacking to S2-i does not implement any preference for a continuation of the S2-i stacking pattern, i.e. stacking faults do not have any range.This is indicated by the absence of additional, sharp reflections in the measured diffraction pattern of H 3 LiIr 2 O 6 .Hence each line of the transition probability matrix     The results of the systematic DIFFaX-study are summarized in Fig. 13.As predicted, increasing faulting in the S1-stacking of the crystal structure of LiIr 2 O 3 (OH) 3 leads to broadening of all non 00l, 13l, 1 ˉ3l and 3 ˉ3 ˉl reflections.The higher the degree of faulting, x, the greater the similarity between simulated and measured diffraction patterns.In addition, triangular peaks shapes evolve by broadening of the 110, 1 ˉ11, 021, 2 ˉ4 ˉ1 and 240 reflections.At x = 0.4 the greatest similarity between simulated (Fig. 13, dark magenta line) and measured (blue line) XRPD pattern is reached.A further increase of the degree of faulting, which is identical with a complete random stacking using S1 and S2-i vectors, will lead to disappearance of the remnants of the 1 ˉ11, 021, 041 and 240 reflections.As only a finite number of crystals, having a limited crystalline size and therefore a limited number of layer-to-layer-transitions, contributed to the measured diffraction pattern and as recursive DIFFaX-simulations led to diffraction patterns, produced by an infinite number of crystals, there can't be a complete match between simulated and measured patterns.In addition in the DIFFaXsimulation only idealized stacking vectors were used, that were created on the assumption, that the O-H bond is orientated perpendicular to the layers, which is not necessarily the case, taking the asymmetric coordination spheres of the oxygen sites (2x Ir, 1x Li) into account.Nevertheless, a good agreement was reached, which confirms that a realistic microstructural model has been developed.The degree of faulting in the S1-stacking pattern, x = 0.4, that created the best agreement between the simulated and measured pattern, however, is most likely slightly different in the real microstructure of the sample.Hence a Rietveld compatible approach was used, that allowed the refinement of the stacking vectors.

Rietveld compatible approach
In order to confirm the microstructural model derived from structural considerations and DIFFaX-simulations the Rietveld compatible approach [38][39][40] was used for global optimization of a limited supercell representing the microstructure of the sample.In this approach each layer is treated as a rigid body.
During global optimization the layers are allowed to move unconstrained and unrestrained in their ab-plane.A 12 c supercell containing 12 layers was subjected to global optimization using the simulated annealing approach 41 implemented into TOPAS.The resulting stacking vectors were calculated from each layer position in relation to the preceding layer with the first layer of the unit preceded by the last one.
Although the movement of the layer was not restrained each stacking vector within the supercell refined very close (≈±0.03x,±0.02y) to an idealized pendant derived from structural considerations (Table 8).This means that the free unconstrained movement of the layers led to a stacking Table 8 Comparison of stacking vectors obtained by global optimization of a 12 c supercell of H 3 LiIr 2 O 6 using uncon-and unrestrained movement of 12 rigid body like layers in xand y-direction with their idealized pendants derived from structural considerations pattern in which all oxygens of neighbouring layers are almost on opposition to each other.S1-like stacking is dominant and small S1-stacked sections are present, which breaks the hexagonal symmetry, whereas transitions to and among S2-like stacking occur unsystematically.The obtained supercell provides good agreement between calculated and measured diffraction pattern with good agreement factors (Fig. 14).

PDF-analysis
In order to support the microstructure model of H 3 LiIr 2 O 6 obtained by a Rietveld compatible approach using low energy (22.16 keV ≡ 0.5594 Å) laboratory X-ray diffraction data PDFanalyses were performed using high energy (67.13 keV ≡ 0.1847 Å) synchrotron XRPD data collected at the X28 beamline of NSLS-II (Brookhaven Nat.Lab.).PDF-curves were calculated (Fig. 15 1 and 2) and from the 12 c (12 layer) supercell structure (Fig. 15c) by a Rietveld compatible global optimization using rigid body like layers (Table 8, Fig. 14) and compared to the measured curve (Fig. 15, blue circles).Atomic and lattice parameters were kept fixed for all calculations.
PDF-analysis clearly demonstrates that the HCrO 2 like structure model of O'Malley et al. 1 is not suitable as the match between measured and calculated curve is poor both in low and high distance region (Fig. 15a).The exact position of the first peak at ≈2.0 Å, referring mainly to Ir-O pairs cannot be described by this structure model, which means the constitution of the layers differs from the expected HCrO 2 like structure model.In addition the misfit at higher distances indicates a different stacking pattern, as well.In contrast for the monoclinic structure model with pure S1-stacking there is a good match between calculated and measured curve at low distances (r < 7.0 Å) indicating that the constitution of the layers was refined properly (Fig. 15b).At high distances referring to interlayer pairs, however, there are serious deviations of the calculated PDF-curve from the measured one (green ellipses), which can be attributed to the occurrence of stacking faults. 40The usage of the structural data of the faulted 12 c supercell leads to a substantial improvement of the fit especially in the region of high interatomic distances (r > 7.0 Å).The match in this region is almost perfect.Thus the microstructure model which was obtained by structural considerations (Fig. 11, Table 4), DIFFaX-simulations (Table 6, Fig. 13) and Rietveld compatible global optimization using rigid body like layers (Table 8, Fig. 14) describes essential features of disorder in the crystal structure of H 3 LiIr 2 O 6 .8, Fig. 14).

Fig. 1
Fig. 1 Packing diagram in the crystal structure of α-Li 2 IrO 3 (left) that can be transformed into LiIr 2 O 3 (OH) 3 ( packing diagram, right) by acid treatment.The stacking order of the anion sublattice is indicated by capital Latin letters.† Electronic supplementary information (ESI) available.See DOI: 10.1039/ c7dt02978k

Fig. 2
Fig. 2 Excerpt of the final Rietveld refinement of diffraction pattern of H 3 LiIr 2 O 6 using the predicted crystal structure by O'Malley et al. 1 as starting model.

Fig. 3
Fig. 3 Comparison of the measured XRPD-pattern of H 3 LiIr 2 O 6 including selected reflection indices assuming C2/m space group symmetry, reflections that are vastly broadened due to faulting are indicated by grey font color and a simulated diffraction pattern of well-ordered Li 2 IrO 3 , for the simulation the instrumental function of the diffractometer used for the XRPD measurement, was applied.

Fig. 4
Fig. 4 Excerpts of the measured XRPD pattern of H 3 LiIr 2 O 6 including selected reflections indices.

Fig. 5
Fig. 5 Low r-region of the measured PDF-curve (blue circles) of synthesized H 3 LiIr 2 O 6 , contributions of Ir-O, Ir-Ir and Ir-Li pairs are highlighted.

Fig. 6
Fig. 6 Scattered X-ray intensities of H 3 LiIr 2 O 6 , at ambient conditions as a function of diffraction angle 2θ.The observed pattern (circles) measured in Debye-Scherrer geometry, the best Rietveld fit profiles (line) using the LiIr 2 O 6 3− -layers of α-Li 2 IrO 3 as starting model and the difference curve between the observed and the calculated profiles (below) are shown.The high angle part starting at 31.0°2θ is enlarged for clarity.

Fig. 11
Fig. 11 Derivation of potential stacking vectors in the microstructure of H 3 LiIr 2 O 6 by considering all potential O⋯H⋯O contacts of adjacent layers (a), projection of the oxygen sites at the bottom site of subsequent layer onto the ab-plane (b).

Table 4
Overview on potential stacking vectors in the microstructure of H 3 LiIr 2 O 6 derived by considering all potential O⋯H⋯O contacts of adjacent layers , black line) und compared to the measured pattern (blue line) of the sample.In the measured diffraction pattern only sharp 00l, 13l, 1 ˉ3l and 3 ˉ3 ˉl reflections are apparent (Fig. 12, blue indices).All other reflections are broadened in such a way that they almost completely merge with the background.Broadening of the 110, 1 ˉ11, 021, 2 ˉ4 ˉ1 and 240 reflections leads to a characteristic triangular peak shape 32 that indicates heavy faulting by random transitions among several stacking vectors.In addition, no further, sharp reflections are present in the measured diffraction pattern, hence the existence of domains of homogenous, S2-i stacked sections in the microstructure of H 3 LiIr 2 O 6 can be excluded.There are some tiny remnants of the 1 ˉ11, 021, 041, 022 and 240 reflections in the measured diffraction pattern, which indicate that despite heavy faulting S1-stacking is still present in the microstructure and that small S1-stacked sections are apparent.

Table 5
Transformation of the unit cell of H 3 LiIr 2 O 6 into a DIFFaXcompatible setting c′ = c•cos(β − 0.
, red line) from the HCrO 2 like crystal structure model of O'Malley et al. 1 (a), from a faultless monoclinic structure model with pure S1-stacking (b) obtained by a constraint Rietveld refinement (Fig. 6, Tables

Fig. 14
Fig. 14 Scattered X-ray intensities of H 3 LiIr 2 O 6 , at ambient conditions as a function of diffraction angle 2θ.The observed pattern (circles) measured in Debye-Scherrer geometry, the best Rietveld fit profiles (line) using a 12 c supercell and the difference curve between the observed and the calculated profiles (below) are shown.The high angle part starting at 21.0°2θ is enlarged for clarity.
Open Access Article.Published on 25 September 2017.Downloaded on 12/23/2018 11:04:21 AM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.below).According to the limits of the XRPD method no attempt was made to determine and refine the hydrogen positions.The stacking pattern of the LiIr 2 O 6 3− -layers of H 3 LiIr 2 O 6 in the structural model is determined by the lattice parameters and the space group symmetry.According to the pronounced occurrence of stacking faults, the obtained structural model cannot be considered as the crystal structure of H 3 LiIr 2 O 6 , but it will serve as a starting point for a detailed investigation of the microstructure of H 3 LiIr 2 O 6 .Crystal structure description of H 3 LiIr 2 O 6 in relation to α-Li 2 IrO 3 Dalton Trans., 2017, 46, 15216-15227 | 15219

Table 2
Atomic coordinates of H 3 LiIr 2 O 6 at room temperature

Table 1
Space group symmetry and lattice parameters of α-Li 2 IrO 3 1 and H 3 LiIr 2 O 6 (obtained by the LP-search routine)

Table 3
Comparison of the metal-oxygen and metal-metal distances in the crystal structures of α-Li 2 IrO 3 9 and H 3 LiIr 2 O 6 Open Access Article.Published on 25 September 2017.Downloaded on 12/23/2018 11:04:21 AM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.

Table 7
Constrained transition probability matrix, used for a systematic DIFFaX-study, with increasing probability, x, of a fault within S1-stacking Open Access Article.Published on 25 September 2017.Downloaded on 12/23/2018 11:04:21 AM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.canbe expressed by the same set of transitions probabilities (Table In consequence a systematic DIFFaX study could be carried out by varying only one parameter, x, describing the probability of faulting in the S1-stacking pattern of H 3 LiIr 2 O 6 .