Robust singlet dimers with fragile ordering in two-dimensional honeycomb lattice of Li2RuO3.

When an electronic system has strong correlations and a large spin-orbit interaction, it often exhibits a plethora of mutually competing quantum phases. How a particular quantum ground state is selected out of several possibilities is a very interesting question. However, equally fascinating is how such a quantum entangled state breaks up due to perturbation. This important question has relevance in very diverse fields of science from strongly correlated electron physics to quantum information. Here we report that a quantum entangled dimerized state or valence bond crystal (VBC) phase of Li2RuO3 shows nontrivial doping dependence as we perturb the Ru honeycomb lattice by replacing Ru with Li. Through extensive experimental studies, we demonstrate that the VBC phase melts into a valence bond liquid phase of the RVB (resonating valence bond) type. This system offers an interesting playground where one can test and refine our current understanding of the quantum competing phases in a single compound.

Systems with strongly correlated electrons usually harbour rich magnetic properties1. Most often these are different types of long-range magnetic ordering. However, other options are also possible. One of them, widely discussed and still being extensively studied, is the formation of different types of spin liquid states, which is generally expected to be realised in frustrated systems2. Yet another possibility is that singlet bonds are formed in the system, leading eventually to the emergence of valence bond crystals or valence bond solids. These are, for example, Peierls and spin-Peierls states in one-dimensional systems, but such states can also exist in higher dimensions.

The exact conditions in which such a VBC can be formed are not well known, although some general trends have been noted. One likely possibility is their formation in low-dimensional systems3. Orbital degrees of freedom may also help to stabilize such states45 via a particularly intriguing mechanism of orbital-selective dimerization6. Although interesting in its own right, the details of how the VBC are formed and stabilized remain largely unexplored. Equally fascinating questions are how the crossover between the usual magnetic ordering and the VBC occurs, and what the possible metastable states are close to such a crossover. There are in principle two possibilities: either there exist a quantum phase transition between these states, or the transition between them could be of first order. In the latter case, one could expect a possible coexistence of both states, which can imply, in particular, the existence of local dimers close to such transitions.

Another equally interesting and yet poorly understood question is how the VBC with a charge gap responds to low levels of doping. We are not aware of any studies of such a crossover or of the detailed doping effects in a real material under controlled conditions. We now identify one system which seems to be ideal for this kind of study: Li2RuO3, a layered material with a honeycomb lattice.

Systems with a honeycomb lattice have recently attracted special attention. The most celebrated case is of course graphene, with its Dirac points in the electronic spectrum. Correlated systems with honeycomb lattices are also generating considerable interest. The best known examples are Li2IrO3 and Na2IrO3, for which the applicability of a Kitaev-Heisenberg model was proposed78, with the eventual formation of a special spin-liquid state. Unfortunately, however, different magnetic ordered states were found in real materials: a zigzag-type ordering in Na2IrO39 and an incommensurate magnetic ordering in γ-Li2IrO310. We note that other types of magnetic ordering were also found in related compounds with similar structures: Li2MnO31112 and Li2RhO313. However, at least in one system of this class, Li2RuO3, the situation is drastically different. It was shown by Miura et al.1415 that in polycrystalline Li2RuO3 there is a transition of the singlet Ru-Ru dimer formation at Tc ~ 540 K, below which these dimers order in a herringbone fashion. Thus this material could be a classic example of the VBC state. As argued in ref. 16, orbital ordering seems to play a crucial role in the formation of this VBC state.

The situation in Li2RuO3 is far from trivial and much richer than originally thought. On the one hand, the singlet dimers with short Ru-Ru bonds seem to be very stable. As demonstrated by recent total scattering experiments and PDF (Pair Distribution Function) analysis17, the Ru-Ru dimers survive up to temperatures much higher than its Tc. On the other hand, a recent study of Li2RuO3 single crystals reported quite different behaviour. According to this latter study18, depending on the exact preparation conditions some single crystals show only a much weaker transition at about the same temperatures as reported in ref. 14, whereas other crystals show no transition at all, apart from a weak magnetic ordering at ~ 5 K18. Exact reasons for the sample dependent behaviour are not known at the moment, but one can guess that it may be due to small deviations in stoichiometry, which is more difficult to control in single crystals grown using the flux method, than in powders.

To have a better understanding of the situation, one has to answer the questions of how the VBC ground state evolves upon doping or introducing disorder in a real 2D material and how the magnetic ground state eventually emerges. This is a nontrivial problem and warrants careful study. In this work, we carried out a detailed study of this system with controlled changes of disorder and employing a comprehensive set of experimental techniques: resistivity, magnetization, specific heat, and both high-resolution elastic and inelastic neutron scattering (see Methods). On the basis of these studies, we conclude that even small levels of disorder, which may be expected to have a negligible effect on the ground state, are found to drastically modify the long-range ordered spin dimer state and thereby influence many of the physical properties of this system. Of particular interest is that all this happens despite the fact that the singlet dimers seem to survive locally even at higher doping. Thus the resulting state may be visualized as some kind of a dimer liquid, but with certain dimers broken. This, in particular, produces unusual spatial modulations in the hopping integrals leading to a variable-range hopping (VRH)-like conduction at low temperatures. Another nontrivial effect observed experimentally is the appearance of low-energy excitations giving rise to linear specific heat at low temperatures: , with a large coefficient γ. The broken dimers also seem to lead to enhanced magnetic signals, observed in the uniform susceptibility measured by bulk magnetization and in the dynamic susceptibility obtained from inelastic neutron scattering experiments. Our ab-initio GGA (generalized gradient approximation) calculations support this picture and, in particular, show the formation of a magnetic cloud close to impurities (in our case extra Li on the Ru sites).

At high temperatures Li2RuO3 with Ru4+ forms in the C2/m space group, one of several monoclinic phases common to this type of transition metal oxides, with the Ru occupying the symmetric honeycomb lattice with Li at the centre of Ru hexagons (Fig. 1a,b). These Ru honeycomb layers are separated by another layer of Li atoms (Fig. 1b). Upon cooling below 550 K, Li2RuO3 undergoes a structural transition into another monoclinic phase of P21/m symmetry by losing the face centring of the monoclinic plane. As shown in Fig. 1c,d, this transition is accompanied by a strong off-centring of Li atom breaking the symmetric honeycomb networks, and by splitting of three otherwise equal Ru−Ru nearest neighbour distances on the honeycomb lattice into one short and one intermediate and one long bond per unit cell (Supplementary Fig. 1 and Table 1). This off-centring and lowering of the symmetry was originally interpreted as being due to Ru dimerization caused by orbital ordering shown in Fig. 1e16. Also notable is the change of the unit cell parameters, and in particular the shrinkage of the lattice in the direction of the short Ru-Ru dimers, thus making the a (b) lattice constant shorter (longer), respectively (Supplementary Fig. 1 and Table 1).

Figure 1

Schematic diagram and in-plane view of the structure.

(a) Top and (b) side view of the crystal structure of P21/m space group. The Ru honeycomb network changes from (c) a structure at 550 K with the C2/m space group to (d) a structure with Ru dimer formation with the P2/m space group with the Ru orbital wave function as shown in (e).

All the physical properties of our samples prepared as described in the Methods section exhibit drastic changes at the structural transition as shown in Fig. 2. For example, the resistivity shows a marked increase, reminiscent of a metal-insulator transition seen in other Ru oxides19. At the same time, the magnetic susceptibility displays a sudden drop at almost the same temperature as reported previously1415. What is intriguing is that the low temperature resistivity of all our samples does not follow the usual Arrhenius law although the high temperature phase seems to be more consistent with the activation behaviour with a charge gap energy of around 200 meV (see Fig. 2a). The low temperature behaviour is clearly closer to variable-range hopping behaviour. Another interesting feature is that the magnetic susceptibility shows persistent paramagnetic signals of 2.5×10−4 emu/mole over a very wide temperature range of almost 500 K below the transition, implying a nonzero density of states in the spin susceptibility (see Fig. 2b). This then appears to be at variance with the picture of a complete dimerization as suggested earlier1415.

Figure 2

Bulk properties of our samples.

(a) Resistivity and (b) magnetization, and (c) the low-temperature specific heat for seven Li2RuO3 samples with different doping ratio (x) values, together with the data for Li3RuO4 and Li2TiO3. The dashed lines in (a) represent the fitting results using the activation formula. The insets in (a,b) show the resistivity and susceptibility data versus temperature for three representative samples, respectively.

The Ru honeycomb lattice can be perturbed by introducing more Li atoms into the Ru honeycomb lattice: for example, one can replace Ru by Li and vice versa Li2+Ru1–O3. An extreme example of such disruption of the Ru honeycomb lattice is Li3RuO4, which has one-dimensional zig-zag chains of Ru with Ru5+ valence, and is known to order antiferromagnetically at 40 K20. In order to investigate disorder effects on the metal-insulator transition (MIT) of Li2RuO3, we have made careful experiments, controlling the levels of disorder of the Ru honeycomb lattice by replacing Ru with Li ions and vice versa. In order to determine the exact amount of mixing (x) we carried out high-resolution neutron/X-ray diffraction experiments (Supplementary Table 4) and compared it with our best sample (LRO1 with a nominal value of x ≅ 0) to find that there is a monotonous variation in the unit cell volume with the amount of disorder for the x values of our interest. This picture of disorder at the Ru honeycomb lattice by Li is confirmed by the observation that as we increase the doping ratio (x) between Ru and Li in the Ru honeycomb lattice with increasing doping, the lattice parameter a increases while the lattice parameter b decreases (Supplementary Table 3).

As shown in the resistivity and magnetization data in Fig. 2a,b, the transition is progressively suppressed with increasing x. Several things are noteworthy here. First, the transition temperature does not change much while the transition itself becomes significantly broadened with doping. Then, the variable-range hopping behaviour seen below the transition temperature of Li2RuO3: with the exponent α ~ 1/3 and the charge excitation gap of Δ = 30 K for x = 0, transforms into a more insulating state with an activation energy increasing with x (Fig. 2a). For comparison, there is a clear activation behavior with Δ = 320 K for Li3RuO4 (Fig. 2a). We note in passing that the magnetoresistance of Li2RuO3 measured up to 14 tesla is always positive. Simultaneously, the drop in the susceptibility becomes less evident with increasing disorder, before recovering behaviour more reminiscent of Curie-Weiss localized moments for the Li3RuO4 sample, which shows the antiferromagnetic transition at 40 K as reported previously. We comment that according to our high-resolution X-ray diffraction data there is a sign of a miscibility gap for samples with higher disorder than reported here.

As regards the disorder effects on the MIT, it is worth mentioning that all our samples, including a sample (x = 0.13) labelled DTA that was prepared under slightly different conditions (see Methods), show the metal-insulator transition with more or less the same activation behaviour above the transition while there are clear variations in their low temperature behaviour (Fig. 2a). This sample dependence of the low temperature resistivity indicates that the low temperature phase is rather sensitive to low levels of disorder, while the transition itself appears to be more robust to modest doping. This experimental observation is rather remarkable given the high transition temperature.

What is more striking is the low temperature behaviour of the heat capacity shown in Fig. 2d. For our most stoichiometric sample (LRO1 with x = 0) with the highest transition temperature, the low temperature heat capacity includes a very small but finite linear contribution, with a γ value of 0.87 mJ/mole-K2. With increasing the mixing ratio (x), the γ value increases monotonically and reaches a large value of 40 mJ/mole-K2 for x = 0.22. For comparison, we measured Li2TiO3 prepared under similar conditions and found that it has the γ value of 0.09 mJ/mole-K2. Therefore, this measured γ value of 40 mJ/mole-K2 for x = 0.22 suggests that the low temperature phase should be a highly correlated phase. To put this value into perspective, let us assume that this γ value directly scales with the degree of disorder (which is confirmed by our data, Fig. 3) and that for x = 0.22 it comes from ~20% of the total volume of the sample. To stress that this γ value is indeed very large, we can consider a hypothetical case of 100% doping: of course, not realizable in practice. For such a hypothetical 100%-doped case, the rescaled γ value would be 200 mJ/mole-K2, which is on par with the typical values for heavy fermion systems21.

Figure 3

Doping dependence of key experimental parameters.

(a) the MIT transition temperature, (b) the charge gap estimated from the resistivity data above the transition, (c) the linear temperature dependence to the specific heat and (d) the paramagnetic contribution of the susceptibility at 300 K.

In Fig. 3, we summarize the main experimental findings (transition temperature, the gap value estimated from the high temperature resistivity data, the γ value, and the paramagnetic susceptibility) as a function of the unit cell volume: the unit cell volume is found to increase almost linearly with the doping ratio (x), thus these plots can be taken to show the dependence of these parameters on x. This summary reinforces our view that the transition temperature is only slightly suppressed with doping while the γ value and the room-temperature susceptibility increase noticeably.

The detailed origin of the most striking observation: linear specific heat, deserves special discussion. As noted above, the specific heat at low temperatures in modestly doped samples (up to x ~ 0.2) shows a linear temperature dependence: C = γT, with a rather large value of γ. This kind of linear temperature dependence in specific heat is usually ascribed to the presence of the Fermi-surface, with low-energy (electronic) excitations22 with a finite density of states: γ ~ N(EF), where EF is the Fermi energy. The observed linear specific heat data at low temperatures could be due to the presence of electronic excitations, e.g. the formation of inhomogeneous (phase-separated) state with metallic droplets immersed in the insulating matrix. We cannot exclude this possibility entirely at the moment, although it seems rather unlikely considering all the experimental observations, in particular the huge value of γ anticipated for such an imaginary metallic state occupying the whole sample as discussed above.

However, this is not the only possibility. Some other excitations with this property (a finite density of states at zero energy) could also give rise to such linear contributions to specific heat. This is indeed, for example, the case in some disordered systems23. In this sense, we can offer an alternative explanation, which is to attribute this linear specific heat to some other excitations in the magnetic subsystem. The ground state of Li2RuO3 consists of the ordered arrangement of singlet dimers, i.e. it is a valence bond crystal (VBC). The total scattering experiments with the PDF analysis17 demonstrate that at T > Tc, i.e. above the structural transition, dimers still exist locally, forming something like a (classical) dimer liquid. One can assume that similar state can also be generated at low temperatures by certain local disorder, i.e. by extra Li replacing some of the Ru ions in the honeycomb layers. In such a state, there should exist real magnetic (singlet-triplet) excitations, contributing to the magnetic susceptibility, but these may have a rather large energy gap (the singlet binding energy). But it is plausible that random and dynamic distribution of the dimers also allow for singlet excitations. It is known, for example, that in some frustrated magnetic systems, e.g. in Kagome magnets, there exist a lot of low-energy singlet excitations, which are accumulated at zero energy with increasing system size24. We can assume that similar excitations may also be possible in our disordered Li2RuO3 samples with suppressed long-range dimer ordering, but with dimers surviving locally, and that such excitations can make a significant contribution to the linear specific heat (see also the discussion below, after presentation of neutron scattering data in Fig. 4). But, the details of this picture are still unclear, and we cannot exclude that at least part of this specific heat comes from real electronic excitations.

Figure 4

Spin dynamics measured by inelastic neutron scattering of two Li2RuO3 samples.

(a) DTA (x = 0.13) & (b) LRO2 (x = 0.07) samples, (c) the difference (DTA-LRO2) plot and (d) the momentum average scattering response as a function of energy for both samples and (e) the difference plot of the momentum average data. We fitted the difference data in (e) using two Lorentzian functions (dash-dot line) with the sum of the two given in the solid line.

Further insight into the nature of the disorder-induced state can be obtained from measuring the excitations directly by using techniques such as inelastic neutron scattering on Li2RuO3 samples with different amount of disorder. For that purpose, we have measured the spin dynamics of two samples with different Li contents, i.e. a different amount of disorder: one is a sample with less disorder (LRO2 with x ~ 0.07), and the second is a slightly more disordered sample with x ~ 0.13 (DTA) (see Methods). As shown in Fig. 4a, the inelastic neutron scattering data of the DTA sample measured at 5 K exhibits strong scattering over the energy range from 2 to 6 meV. On the other hand, this scattering is strongly suppressed in the data taken on the LRO2 sample (x ~ 0.07) with less disorder as shown in Fig. 4b. That is, the low energy excitations observed in our inelastic neutron scattering experiments are clearly induced by a small amount of disorder, i.e. Ru on the Li site. This conclusion is further supported by the difference taken between the two data sets (Fig. 4c), which is obtained by directly subtracting the LRO2 data from the DTA data.

We note that our subsequent measurement on another sample (LRO5 with x ~ 0.16) reproduces exactly the same behaviour of strong magnetic scattering, as shown in the supporting information (Fig. SI5): the LRO5 sample was prepared under a more control protocol as described in the Methods section. Integrating the inelastic neutron scattering data of the DTA sample over the first Brilliouin zone suggests that approximately 0.5 μB/f.u. of magnetic moments are involved in the low energy excitations. It ought to be noted that as shown in Fig. SI6, these magnetic excitations are significantly weakened with increasing temperature, although they are still visible even in the data taken at room temperature.

Of further interest is that the low-energy excitations demonstrate a clear momentum modulation, which can be explained by the nearest neighbour correlation as shown in the supporting information (Fig. SI7). These correlations are probably connected with dimer correlations as seen in the total scattering measurements17. The other point worth mentioning is that the uniform susceptibility calculated from the inelastic neutron scattering data is in good agreement with the bulk data as shown in the supporting information (Fig. SI8). This latter observation reinforces our view that the unusually enhanced low-temperature susceptibility and the γ value are intrinsic and arise from the low-energy excitations measured by our inelastic neutron scattering experiments. A further confirmation can be found in the so-called Wilson plot as shown in Fig. SI9, the Sommerfeld-Wilson ratio is less than one for our Li2RuO3 materials, implying that strong correlations are present in our samples.

The picture emerging out of the magnetic behaviour of Li2RuO3 with disorder resembles that of the lightly doped canonical spin-Peierls system CuGeO3. Doping of CuGeO3 by nonmagnetic Zn breaks some of singlet pairs, and the resulting unpaired spins polarize the remaining singlet dimers, eventually leading to an inhomogeneous magnetic order2526. We can expect similar behaviour to occur in Li2RuO3, with the key difference being that because of a very large binding energy of singlet dimers (in our case Tc ~ 500–600 K, instead of Tc ~ 14 K in CuGeO3), the extension of the magnetic cloud around impurities would be much smaller in Li2RuO3.

To check this hypothesis we performed ab-initio GGA calculations and simulated the 4% Li/Ru interchange by constructing an appropriate supercell (see Methods). Such an interchange results in two types of defects. The first kind of defects are hexagons with a Ru instead of a Li atom at the centre, and the second kind of defects are Ru dimers broken by the substitution of Li atoms. We found that this breaking up of the Ru dimers leads to significant changes in the magnetic properties of the system. In an ionic model, unpaired magnetic Ru4+ ions would have S = 1 and hence a magnetic moment of 2 μB, while our GGA calculations show the magnetic moment of ~1.2 μB on this unpaired Ru atom. In addition, there are two other Ru ions next to the Li ion breaking the dimer, which themselves are magnetized by this defect with the induced magnetic moments of ~0.7 μB. We note that the total change of the magnetic moment with the Li/Ru interchange found in our calculations is ~0.25 μB/f.u., consistent with our experimental findings. The difference in the spin densities for pure and 4% doped Li2RuO3, shown in Fig. SI10, clearly demonstrates a formation of the magnetic cloud in a vicinity of defects. Thus the obtained theoretical results support the picture of short-scale magnetic cloud close to the impurity, as described above. A passing comment, we came across recently a report on the disorder effect in Li2RuO327.

Summary

To summarize, our detailed and extensive studies on the disorder effects of Li2RuO3 paint two apparently contradicting, yet rather revealing pictures of the spin dimerization. The Ru−Ru singlet dimers that form a long-range ordered state (valence bond crystal) for pure Li2RuO3 below Tc ~ 540 K, are quite robust. And yet, at the same time, such valence bond crystal itself appears to be very fragile and can be suppressed by rather low levels of disorder. Thus the state of Li2RuO3 with a finite level of disorder seems to correspond not to a valence bond crystal, but rather to a valence bond liquid (or valence bond glass), similar in spirit to the resonating valence bond state proposed originally for frustrated magnets28. Our work shows that Li2RuO3 is a very convenient and interesting material, providing good playground, on which one can test and improve our understanding of the unique transition from a valence bond crystal state to a valence bond liquid state, a quite nontrivial quantum state of matter, and eventually to a magnetically ordered state. Thus, it provides a rare window of opportunity to study the question of the destruction of quantum-entangled states in a real material.

Methods

Sample preparation & physical properties

In order to make comprehensive studies of doping experiments on the dimerized ground state of Li2RuO3, we prepared 7 samples with different levels of disorder using a solid state reaction method. All our samples were made with the starting materials better than 99.99% purity by mixing stoichiometric amounts of Ru or RuO2 and Li2CO3, preheated at 600 °C in air for overnight. The powder was then pressed into 10 mm diameter pellets and fired at 900 °C for 15 hr, before being further sintered in air at 1000 °C from 48 to 200 hr with intermediate grinding. The exact final sintering condition of our samples, all labelled with the prefix LRO, is summarized in Table SI2. In addition, one sample labelled DTA was also synthesized at ISIS, UK, by mixing commercially available powders of RuO2 and Li2CO3 following previously described procedures14. The starting raw materials were sintered in air at 1000 °C for 24 h in an alumina crucible, after which the product was reground and pelletized and heated at 900 °C for 48 h.

As well as noting the nominal starting composition of the six (LRO) samples, the exact disorder or doping (x) in these samples was determined from an analysis of our high-resolution neutron diffraction data (Supplementary Fig. 4 and Table 4). The values of x obtained suggest that at these low levels of doping the unit cell volume is approximately linear in x as shown in Fig. 3. The sample purity was monitored by collecting powder x-ray diffraction patterns (using a Miniflex II, Rigaku) with Cu Κα radiation. The local structure of our samples was examined by using the XANES (X-ray absorption near edge structure) technique at the 3C1 beam line of the Pohang Accelerator Laboratory and further high-resolution x-ray diffraction experiments were conducted at 8C2 beam line of the Pohang Accelerator Laboratory and using a high-resolution X-ray diffractometer (or Bruker XRD D8 Discover diffractometer) (Supplementary Figs 2 and 3). We also carried out microprobe chemical analysis by using ICP (inductively coupled plasma) & EPMA (electron probe micro-analyzer) techniques, and confirmed the chemical variations as discussed in the text.

We carried out low and high temperature resistivity measurements using a home-made setup covering the temperature range from 3 to 650 K. We also measured magnetization and heat capacity measurements using commercial set-ups: MPMS-5XL (Quantum Design, USA), PPMS9 (Quantum Design, USA), and VSM (Lakeshore Instrument). Magnetoresistance measurements were made up to 14 tesla using a home-made set-up at Tohoku Univ.

Neutron scattering

In order to investigate the structural transition in detail, we have also undertaken high-resolution neutron diffraction studies on two samples (LRO2 with x ~ 0.07 & DTA with x ~ 0.13) with slightly different levels of disorder by employing two neutron powder diffraction beam lines (HRPD and GEM) at the ISIS, UK, from 300 to 590 K. We also measured the spin dynamics of three samples (LRO2, DTA & LRO5 with x ~ 0.16) at the MERLIN and MARI time-of-flight inelastic neutron scattering beam lines of ISIS, from 5 to 580 K.

GGA calculations

The ab-initio calculations were carried out by pseudopotential method in the Quantum Espresso code29. We used the Generalized Gradient Approximation (GGA) with exchange-correlation potential as proposed by Perdew, Burke, and Ernzerhof 30. The cut-off energies for the wave functions and charge density were chosen to be 40 and 180 Ry, respectively. The crystal structure was taken from ref. 14 for T = 300 K. In order to simulate the Li/Ru interchange we constructed supercell, consisting of 144 atoms in the same hexagonal plane.

Additional Information

How to cite this article: Park, J. et al. Robust singlet dimers with fragile ordering in two-dimensional honeycomb lattice of Li2RuO3. Sci. Rep. 6, 25238; doi: 10.1038/srep25238 (2016).