Chemical and hydrostatic-pressure effects on the Kitaev honeycomb material Na2IrO3 1, 2,† 3 2 4 1 G. Simutis, ∗ N. Barbero, K. Rolfs, P. Leroy-Calatayud, K. Mehlawat, R. Khasanov, H. Luetkens,1 E. Pomjakushina,5 Y. Singh,4 H.-R. Ott,2, 6 J. Mesot,2, 6 A. Amato,1 and T. Shiroka2, 6 1Laboratory for Muon Spin Spectroscopy, Paul Scherrer Institut, Villigen PSI, Switzerland 2Laboratorium für Festkörperphysik, ETH Zürich, CH-8093 Zürich, Switzerland 3Laboratory for Developments and Methods, Paul Scherrer Institut, CH-5232 Villigen, Switzerland 4Indian Institute of Science Education and Research, Sector 81, Mohali 140306, India 5Laboratory for Multiscale Materials Experiments, Paul Scherrer Institut, CH-5232 Villigen, Switzerland 6Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland (Dated: October 11, 2018) The low-temperature magnetic properties of polycrystalline Na2IrO3, a candidate material for the real- ization of a quantum spin-liquid state, were investigated by means of muon-spin relaxation and nuclear magnetic resonance methods under chemical and hydrostatic pressure. The Li-for-Na chemical substitution promotes an inhomogeneous magnetic order, whereas hydrostatic pressure (up to 3.9 GPa) results in an enhancement of the ordering temperature TN. In the first case, the inhomogeneous magnetic order suggests 4 either short- or long-range correlations of broadly distributed j = ½ Ir + magnetic moments, reflecting local disorder. The increase of TN under applied pressure points at an increased strength of three dimensional interactions arising from interlayer compression. I. INTRODUCTION O α β Magnetic frustration resulting from bond-dependent ex- change interactions is a possible route to the formation of Ir a quantum spin-liquid (QSL) state,1 in which quantum ef- γ fects prevent a long-range ordering of magnetic moments even at zero temperature. The Kitaev-Heisenberg model on a single-layer honeycomb lattice with bond-dependent b Ising interactions (see Fig.1), which can be solved exactly a Na in certain cases, is known to host QSL as one of the possible states. The latter has been predicted to exhibit a range of un- conventional features, such as emergent Majorana fermions FIG. 1. Schematic sketch of the honeycomb layer in Na2IrO3. The and gauge fluxes as effective excitations.2 effective spin-½ 5d5 Ir4+ ions (at the center of oxygen octahedra) interact via three distinct NN (nearest neighbor) exchange cou- Real systems featuring the interactions assumed in the plings, here denoted by α, β, and γ. The J and K parameters in Kitaev model may be found in spin-orbit coupled materials Eq. (1) represent the Heisenberg- and Kitaev exchange integrals, 3 with edge-shared octahedra. However, experiments have respectively, taken over all the NN links in the honeycomb. shown that only a few compounds exhibit the relevant Ki- taev interactions required to realize this model. The most prominent examples are the layered honeycomb-lattice iri- 4,5 6,7 dates Na2IrO3 and α-Li2IrO3, as well as the recently 8 identified H3LiIr2O6. In addition, a few three-dimensional level substitution, well below the phase-separation threshold, systems that realize bond-frustrated lattices have been dis- may transform the original zigzag magnetic order into a spi- 9,10 covered, including β and γ polytypes of Li2IrO3. While ral one.13 Subsequently, optical-spectroscopy measurements all the above materials have bond-dependent anisotropic showed that Li-substitution reduces the metal-metal hop- (i.e. Kitaev-type) interactions, also non-negligible Heisen- ping integral t compared to the Coulomb repulsion term U, berg and symmetric off-diagonal exchange interactions, re- thus enhancing the magnitude of the Mott insulating gap.14 sulting in competing ground states, have to be taken into Despite these interesting results, to date a microscopic in- account. In particular, all the known candidate materials, vestigation of the magnetic properties of Na2IrO3 under except H3LiIr2O6, adopt long-range ordered magnetic struc- chemical or hydrostatic pressure is still missing. tures at low temperatures. Since the character of these states depends on the exact details of the underlying Hamiltonian, In this report, we use muon-spin relaxation (µSR) and nu- arXiv:1806.10766v4 [cond-mat.str-el] 10 Oct 2018 different materials exhibit distinct magnetic structures. clear magnetic resonance (NMR) techniques as local probes Considering that subtle modifications may lead to rather and examine the evolution of magnetism in Na2IrO3 upon different ground states, it is intuitively tempting to deliber- chemical doping (Li for Na substitution) and under applied ately induce variations of states or phases, e.g., by changing hydrostatic pressure. In the former case, we find that at the chemical composition or by varying external parameters low substitution levels the magnetic order is rather robust, such as pressure. From this perspective, Na2IrO3 represents whereas above a 5%-substitution threshold, an inhomo- a particularly interesting case. Initial measurements have geneous static magnetic∼ order sets in. In case of applied shown that substituting Na by Li leads to a substantial de- pressure, up to 3.9 GPa, we find a linear increase of the or- 11 crease of the magnetic ordering temperature TN. Later dering temperature with pressure, but no qualitative changes studies pointed out that, at high levels of Li substitution, in the magnetic ground state. The increased TN indicates an phase separation unfortunately prevents a wider range of enhancement of effective exchange interactions, most likely, tunability.12 Nevertheless, it was suggested that even low- arising from a reduction of the interlayer distances.14 2 II. THE KITAEV-HEISENBERG MODEL In addition to µSR, we also performed complementary high-pressure NMR measurements. By employing a hybrid piston-clamped zirconia-anvil cell,22 smaller samples could To date, Na2IrO3 is the only system to have been di- rectly shown to host Kitaev-like interactions.15 To fully de- be probed via NMR at higher pressures (up to 3.9 GPa). The scribe the interactions between the effective spin one-half NMR investigations under pressure included line-shape and Ir4+ magnetic moments (see Fig.1), a so-called extended spin-lattice T1 relaxation-time measurements in an applied Kitaev-Heisenberg Hamiltonian on a honeycomb lattice is magnetic field of 7.057 T. The inner part of the pressure required:16 cell was filled with Daphne oil 7575, acting as a pressure- transmitting medium. The applied hydrostatic pressure was monitored via the pressure-dependent NQR (nuclear X X γ γ J p Si Sj 2K p Si Sj 63 23 = [ ( ) + ( ) + quadrupole resonance) signal of Cu in Cu2O. The most H <i j> α,β,γ · (1) 23 suitable nucleus for our study was Na, an I = 3=2 nucleus α β β α Γ (p)(Si Sj + Si Sj )]. with a 100% abundance and 79.47 MHz Larmor frequency in the chosen applied magnetic field. The Heisenberg exchange coupling J and the Ising-like Ki- taev coupling K define the plain Kitaev-Heisenberg model Hamiltonian. Here the double summation runs over the IV. EXPERIMENTAL RESULTS bond- (α, β, and γ) and NN- (i, j) indexes, respectively. Depending on the ratio between the J and K interaction pa- A. Weak transverse-field µSR measurements rameters, the zero-field ambient-pressure solution predicts 17 six possible ground states. These include two QSL states, 24 a ferromagnetic phase, and three antiferromagnetic phases, Weak transverse-field (wTF) µSR can be used to deter- the latter exhibiting three possible spin arrangements (Néel, mine the magnetic ordering temperature and the magnetic stripy, and zigzag).17 The extended model, represented in volume fraction. In the paramagnetic phase muon spins pre- Eq. (1), includes also a Γ parameter, which captures the cess coherently around the externally applied field. Upon symmetric off-diagonal exchanges and requires a numerical cooling, the electronic magnetic moments order and the solution. Finally, since all the coupling parameters depend resulting magnetic-field distribution at the muon stopping on orbital hybridization, in turn depending on structural sites dephases the muon spins. The ensuing damped os- details, a possible pressure dependence of the coupling pa- cillation can be fitted with a simple harmonic model and rameters is included to adapt the Hamiltonian to our case. the resulting metadata serve to determine the paramagnetic volume fraction Fpara. The latter is related to the parameters of a magnetic-phase transition via:25 III. EXPERIMENTAL DETAILS TN T Fpara(T) = Fnm + (1 Fnm)= exp − + 1 . (2) − ∆ Na2 x Lix IrO3 powder samples with x = 0, 0.05, 0.1, and 0.15 were− synthesized following the procedure described Here Fnm is a residual fraction, corresponding to that part in Ref. 13 and used for the µSR measurements. The solid- of the sample which is not magnetically ordered even at the lowest temperatures, T is the inflection point of each state reaction synthesis used IrO2, Na2CO3, and Li2CO3 as N starting materials, mixed in stoichiometric ratios and heated Fpara(T) curve, interpreted as the transition temperature, 13 26 up to 1000◦C. The original characterization showed a lat- and ∆ is the broadening parameter. The evolution of the tice contraction upon Li doping, with Li uniformly replacing paramagnetic volume-fraction with temperature is shown in the Na ions. Bulk magnetic measurements indicated that Fig.2, both as a function of Li content (a) and as a function both the Curie-Weiss and the transition temperature stays of hydrostatic pressure (b). approximately the same, independent of Li concentration, The best-fit parameters using Eq. (2) are summarized in 11,12 in contrast with reports on single crystals. For the NMR Fig.3. While at low Li-substitution levels ( x = 0.05 and measurements, polycrystalline Na2IrO3 samples were pro- 0.10) the transition temperature and the magnetically-frozen duced using a similar solid-state synthesis, as reported in sample-fraction remain virtually unchanged, for x = 0.15 Ref.