Signature of a randomness-driven spin-liquid state in a frustrated magnet

Collective behaviour of electrons, frustration induced quantum fluctuations and entanglement in quantum materials underlie some of the emergent quantum phenomena with exotic quasi-particle excitations that are highly relevant for technological applications. Herein, we present our thermodynamic and muon spin relaxation measurements, complemented by ab initio density functional theory and exact diagonalization results, on the recently synthesized frustrated antiferromagnet Li4CuTeO6, in which Cu2+ ions (S = 1/2) constitute disordered spin chains and ladders along the crystallographic [101] direction with weak random inter-chain couplings. Our thermodynamic experiments detect neither long-range magnetic ordering nor spin freezing down to 45 mK despite the presence of strong antiferromagnetic interaction between Cu2+ moments leading to a large effective Curie-Weiss temperature of − 154 K. Muon spin relaxation results are consistent with thermodynamic results. The temperature and magnetic field scaling of magnetization and specific heat reveal a data collapse pointing towards the presence of random-singlets within a disorder-driven correlated and dynamic ground-state in this frustrated antiferromagnet. Frustrated magnetic systems are characterised by spin-spin interactions, which are mediated by strong quantum fluctuations, and can potentially lead to magnetically disordered ground states such as a quantum spin liquid. Here, the authors experimentally investigate the frustrated magnet, Li4CuTeO6 and present evidence to suggest that this material may exhibit a random-singlet ground state.

G eometrically frustrated magnets with incompatible magnetic interactions between spins are characterized by strong quantum fluctuations and the absence of trivial order parameters related to broken symmetry. They are, therefore, extremely sensitive to perturbations such as quenched disorder, nextnearest-neighbor interactions, and magnetic anisotropy 1,2 . Strong quantum fluctuations preclude symmetry breaking phase transitions and lead to novel magnetically disordered ground states such as quantum spin liquids (QSLs) 1,3-6 . QSLs are characterized by quantum fluctuations, absence of long-range order even at T → 0, and exotic fractionalized excitations in contrast to spin-wave excitations found in conventional magnets 7,8 . These properties of QSLs are fingerprints of strong quantum entanglement. Triangular-lattice compound κ-(ET) 2 Cu 2 (CN) 3 9 and 1-TaS 2 10 , kagome-lattice compound ZnCu 3 (OH) 6 Cl 2 11 , and hyperkagome-lattice compound PbCuTe 2 O 6 12,13 are some of the most investigated candidate spinliquid materials with unusual ground-state properties. Recently, enormous efforts have been devoted towards experimental realizations of QSLs on the honeycomb lattice ever since Alexei Kitaev proposed an exactly solvable highly anisotropic Kitaev model for S = 1/2 spins on this two-dimensional spin lattice 14 . In sharp contrast to an isotropic nearest-neighbor exchange model on the honeycomb lattice, which exhibits long-range order 15 , the ground state of the Kitaev model is a QSL with fractionalized Majorana fermions and gauge flux excitations [16][17][18][19] . The realization of this exotic state of matter is interesting not only from a fundamental physics point of view but also holds immense promise for potential applications in robust quantum computing technology 20 .
Despite tremendous experimental efforts, the observed signatures of QSL are far away from those of theoretically proposed QSLs states for a clean periodic system 2,14,21 . However, it has been suggested that some highly frustrated magnets show quantum disordered ground state due to randomness that might play an important role in destabilizing conventional Néel order at low T. There are several origins of randomness in magnetic insulators and the interplay between quantum fluctuation and randomness may lead to novel quantum ground states. The triangular-lattice YbMgGaO 4 with antisite mixing between nonmagnetic Mg 2+ and Ga 3+ ions 22,23 and square lattice Sr 2 CuTe 1−x W x O 6 with a few percent of W on the Te site are some of examples of frustrated magnets with a randomnessinduced quantum disordered ground state [24][25][26] . Interestingly, while most of the reported Kitaev magnets including 5d iridates A 2 IrO 3 (where A = Na, Li) and ruthenate α-RuCl 3 with 4d electrons show long-range order at sufficiently low temperature, the iridate H 3 LiIr 2 O 6 and ruthenate α-Ru 1−x Ir x Cl 3 , which are derived from a parent ordered magnet by substitution of nonmagnetic ions, do not exhibit long-range magnetic order [27][28][29][30] . Furthermore, some honeycomb-lattice magnets, like A 3 LiIr 2 O 6 (A = H, D, Ag etc.), show scaling behavior in the field dependence of the specific heat 27,31,32 , which is ascribed to the existence of disorder, making these quantum materials interesting for studies of the effect of disorder on a QSL state 32 . A similar scenario is observed in threedimensional spin lattice where quantum fluctuation is less pronounced. For instance, the oxynitride Lu 2 Mo 2 O 5 N 2 with S = 1/ 2 spins on a pyrochlore lattice, which is derived from the spin-glass material Lu 2 Mo 2 O 7 via random substitution of O 2− and N 3− anions, shows a dynamic ground state 33,34 .
Disorder in quantum magnets is usually unavoidable, yet very often it can act as a new prism in revealing many interesting quantum phenomena of the host material [35][36][37][38] . Frustrated magnets with the quenched disorder in the form of material defects or a broad distribution of exchange interaction strengths can exhibit a randomness-induced spin-liquid state. Understanding the impact of quenched disorder on the quantum fluctuations that drive quantum disordered states in antiferromagnets is an interesting area of study 22,[39][40][41] . In the so-called random-singlet phase, singlets form due to a random distribution of antiferromagnetic exchange interaction between spin-1/2 defect sites, leading to a power-law distribution of exchange energies and density of states, which lead to unconventional scaling behavior in magnetic susceptibility and specific heat at low temperatures 22,32,42 . The signature of such novel random singlet state was initially observed in one-dimensional spinchain organic compounds due to the presence of random antiferromagnetic Heisenberg exchange interaction 43,44 . Furthermore, there is also experimental evidence of a random singlet state in inorganic spin-chain compound BaCu 2 (Si 1−x Ge x ) 2 O 7 45 wherein the random substitution of Ge in place of Si introduce the bond disorder in antiferromagnetic long-range ordered magnet BaCu 2 Si 2 O 7 46 . Theoretically, it has been proposed that such exchange randomness in the spin-chain compound is one of the key ingredients for quantum information processors which deals with properties of quantum entanglement in quantum materials [47][48][49][50] . In recent years, higher-dimensional quantum materials with intrinsic disorder either from spin vacancies, anti-site disorder or nonmagnetic impurities and dislocations or grain boundaries have drawn enormous attention to realize novel low-temperature phases including random singlet that is expected to be unique in disordered quantum many-body systems 25,33,34 . In frustrated magnets with quenched disorder, the majority of intrinsic S = 1/2 spins may constitute a dynamic liquid-like state while the remaining S = 1/2 moments act as defects and form a random network with exotic low-lying excitations that destabilize a glassy state and instead lead to a dynamic low-energy ground state 32 . Theoretically, it has been proposed that randomness in the frustrated lattice can induce spinliquid-like states in a wide region of parameter space 51 . Experimentally, this phenomenon is much less explored due to a lack of suitable frustrated magnets, which leads one to look for candidate materials wherein the degree of randomness can be tuned in a controlled manner. In this context, 3d transition-metal frustrated S = 1/2 spin lattices are promising potential hosts of a randomnessdriven quantum spin-liquid ground state.
Herein, we report crystal structure, neutron diffraction, magnetic properties, specific heat, and muon spectroscopy, as well as complementary density functional theory (DFT) and exact diagonalization (ED) results of the antiferromagnet Li 4 CuTeO 6 (henceforth LCTO), where Cu 2+ (S = 1/2) ions are arranged on a frustrated spin lattice with unavoidable anti-site disorder between Li + and Cu 2+ . In LCTO, the majority of Cu 2+ ions (84%) at the 2d crystallographic site form random-length strongly antiferromagnetic spin chains while a minority of Cu 2+ ions (~7%) at defect 4g sites strongly antiferromagnetically couple to these spin chains, leading to significant frustration. The large and negative value of Curie-Weiss temperature reflects the presence of strong antiferromagnetic interaction between Cu 2+ spins. Furthermore, we find that weaker interchain interactions effectively couple the chains in a random, disordered network at low T (low temperature). As a consequence, this compound neither undergoes a phase transition to long-range magnetic order nor spin-freezing down to at least 45 mK. Furthermore, specific heat and magnetization results reveal data collapse behavior, which suggests the presence of a random-singlet state 32 . Muon spin relaxation measurements corroborate a dynamic ground state in this frustrated magnet. This is attributed to the presence of subdominant interchain interactions that couple the chains in a disordered network. Our results thus establish that LCTO hosts a randomness-induced spin-liquid-like state in a frustrated magnet.

Results
Rietveld refinement and crystal structure. The roomtemperature powder X-ray diffraction data were taken to confirm the phase purity of the polycrystalline samples of LCTO. The Rietveld refinement of the X-ray diffraction data was performed using GSAS software 52 , which confirms that LCTO crystallizes in a monoclicnic structure with the space group C2/m (No. 12) 53 . To deepen the analysis of the crystalline structure and in particular to determine the occupation of the mixed sites by the light elements a neutron diffractogram was measured at 300 K at an incident wavelength 1.28 Å. Figure 1a depicts the Rietveld refinement of neutron-diffraction data performed using fullprof which yields lattice parameters a = 5.2752(5) Å, b = 8.8163(8) Å, c = 5.2457(5) Å and α = γ = 90°, β = 113.172(8)°, in agreement with powder X-ray diffraction data (see Supplementary Table I) and previously reported values 53 . A unit cell of LCTO is shown in Fig. 1b. We found that among three crystallographic sites of Li, the sites 2d and 4g are partially shared with Cu 2 and Cu 3 sites, respectively, with Cu 2 site occupancy p 2 = 84% and Cu 3 site occupancy p 3 = 7%. This is likely due to the similar ionic radii of Li + and Cu 2+ ions. Cu 2+ ions form distorted CuO 6 octahedra with nearest-neighbor oxygen ions with Cu 2 and Cu 3 sites connecting to each other directly by corner-and edge-sharing, and indirectly via corner-and edge-sharing TeO 6 octahedra in an a priori 3D lattice. However, a critical look at the crystal structure reveals only two plausibly strong exchange pathways: (i) a nearly-linear Cu-O-Cu superexchange (SE) bridge between Cu 2 and Cu 3 sites with a bond angle 173.3°and bond distance 4.12 Å, which is expected to yield a strong antiferromagnetic exchange J 0 by Goodenough-Kanamori rules 54,55 , and (ii) a double Cu-O ⋯ O-Cu super-superexchange (SSE) bridge between two Cu 2 sites going around a Te 6+ ion with a large symmetric Cu-O ⋯ O bond angle of 138.6°and a small O ⋯ O distance 2.87 Å (well below the van der Waals radii sum~3.04 Å), which is expected to yield a strong antiferromagnetic exchange J 56-60 (Fig. 1b). The resulting 1D spin model is shown in Fig. 1c and consists of random-length linear spin chains of Cu 2 (occupation probability p 2 ) with exchange J running along the [101] direction, with randomly occupied Cu 3 defect sites (occupation probability p 3 ) displaced along ± b from each main-chain bond J introducing strong frustration by coupling to both neighboring Cu 2 main-chain sites via a J 0 exchange. All other SE and SSE pathways, which, if present, introduce 3D interchain couplings in general, are expected to have much smaller strength as they either proceed via very nearly 90°C u-O-Cu and Cu-O ⋯ O angles, or involve O ⋯ O distances well above the O 2− van der Waals radii sum. The powder ESR spectrum consists of two components, both featuring uniaxial symmetry, in agreement with the crystal symmetry. A fit of the spectrum with a powder-averaged Lorentzian shape (see Fig. 1d) yields principal g factors g 1 x;y = 2.06 and g 1 z = 2.31 for the component with larger intensity and g 2 x;y = 2.06 and g 2 z = 2.20 for the component with smaller intensity. The ratio of the intensities is 85:15, which is close to the ratio p 2 :2p 3 of the occupied Cu 2 and Cu 3 sites as confirmed by the structural refinement.
Ab initio and exact diagonalization calculations. The expected interaction between Cu 2+ moments are fully born out by ab initio DFT exchange-coupling calculations on a 2 × 1 × 2 supercell of LCTO using the local (spin) density approximation (LSDA) in an LSDA+U scheme with an effective Hubbard U eff = 9 meV (see Supplementary Note 6 for details). By first assuming full occupancy of Cu 2 sites and zero occupancy of Cu 3 sites the mainchain Cu 2 -Cu 2 antiferromagnetic exchange strength J DFT = 177 K was extracted, and then by introducing a single occupied Cu 3 site the frustrating Cu 2 -Cu 3 antiferromagnetic exchange strength J 0 DFT = 824 K was obtained. As expected, all other exchange interactions, which introduce 3D couplings between the chains, were found to have strengths below~6 K, making them difficult to resolve by DFT. For T ≫ 6 K the frustrated 1D randomoccupancy model of Fig. 1c is thus expected to hold, while at low T further frustrating interchain exchanges might be expected. Temperature-dependent thermodynamic quantities of the high-T random model obtained from DFT were calculated using ED by probability weighing ED results for different random-length chains of Cu 2 spins with a different distribution of occupied Cu 3 sites around them (see Supplementary Note 7 for details). This way, ED predictions for magnetic susceptibility χ ED , magneticspecific heat C mag,ED and magnetic entropy S mag,ED as a function of T, applied magnetic field μ 0 H (μ 0 is the permeability of free space), g factor, and J and J 0 exchange strengths at any given occupation probabilities p 2 and p 3 could be calculated to high accuracy.
Magnetic susceptibility. Figure 2a presents the experimental temperature dependence of the magnetic susceptibility (χ(T)) in an applied magnetic field μ 0 H = 1 T in the temperature range 1.9 K ≤ T ≤ 300 K. The magnetic susceptibility data do not indicate the presence of a phase transition down to 1.9 K. In order to estimate dominant magnetic interactions, the high-temperature 1/χ(T) data were fitted with a Curie-Weiss law (right y axis; Fig. 2a) In Eq. (1), θ CW and C are the Curie-Weiss (CW) temperature and the Curie constant, respectively. The CW fitting in the temperature range 180 K ≤ T ≤ 300 K yields θ CW = −154 K and C = 0.39 cm 3 Kmol −1 . The estimated effective magnetic moment μ eff = ffiffiffiffiffi ffi 8C p = 1.76 μ B is close to the expected value for a Cu 2+ ion with S = 1/2 spin. The large negative θ CW confirms the presence of dominant antiferromagnetic interactions between Cu 2+ spins. To obtain a more quantitative estimate of the J and J 0 exchange strengths, the χ(T) data were next fitted with the results of ED calculations on the random spin model of Fig. 1c. Fixing the ratio of experimental J=J 0 to be the same as that found in DFT, J DFT =J 0 DFT , which is expected to be a more robust DFT prediction than raw exchange strengths, and fixing the main-chain Cu 2 site occupancy to the experimental value p 2 = 84%, a good fit at T > 20 K is obtained (Fig. 2b) with only 3 free, sensible parameters: g = 2.28, J=J DFT ¼ J 0 =J 0 DFT ¼ 1:38, and p 3 = 20%. We note that while fit parameters g and J=J DFT ¼ J 0 =J 0 DFT are quite robust to perturbations like magnetic anisotropies or unwanted extrinsic contributions due to any impurities, the fitted value of p 3 is substantially less so since the experimental data is limited to T ( J 0 . Thus, the experimental value of Cu 3 occupancy from neutron diffraction should be taken as definitive. The final estimates for the dominant exchanges are thus J = 244 K and J 0 = 1140 K, with the discrepancy from DFT likely arising from the choice of Hubbard U eff , which acts to rescale all exchanges. Importantly, for T < 20 K the high-T model of Fig. 1c starts to systematically deviate from the experiment, indicating the presence of further, most likely interchain, exchanges J inter ≲ 10 K in strength acting between the random spin-chain fragments of the high-T model and dominating the physical response of LCTO at low T < 20 K ( J; J 0 . In this regime most spin degrees of freedom of spin-chain fragments are expected to freeze out, leaving only their ground-state degrees of freedom to couple via J inter exchanges in an emergent ≥ 2 dimensional (likely 3D) random The temperature dependence of magnetic susceptibility at 1 T with a high-T (high-temperature) fit using exact diagonalization (ED) results, and a low T fit using a power-law. The inset shows the zero-field-cooled (ZFC) and filed cooled (FC) susceptibility in a magnetic field of 100 Oe. c The temperature dependence of the real part of ac susceptibility in different frequencies down to 45 mK. d Data collapse for MT −0.15 , where M is the magnetization and T is the temperature, as a function of μ 0 H/T. lattice. Remarkably, though LCTO exhibits such quenched disorder, ultimately arising from partial Cu 2+ site occupancy, that could lead to a spin-glass state at low T, the absence of glass-like freezing behavior is witnessed by the absence of any separation between zero-field-cooled and field-cooled susceptibility measurements in 100 Oe (inset of Fig. 2b). Furthermore, to confirm the absence of spin-freezing, we have performed ac susceptibility measurements down to 45 mK at different frequencies. It is observed that there is no peak or frequency dependence of the real part of the ac susceptibility ( Fig. 2c and see Supplementary  Fig. 3) in the temperature range 45 mK ≤ T ≤ 50 K, which strongly rules out the presence of spin-glass transition in Li 4 CuTeO 6 61 . Quenched Li/Cu anti-site disorder in the host lattice could, alternatively, also lead to unconventional scaling behavior of thermodynamic observables. Such a scenario has been proposed in several frustrated quantum materials 22,25,32,62 .
As shown in Fig. 2d, a data collapse of MT −0.15 , where M is the magnetization and T is temperature, as a function of μ 0 H/T is found in LCTO in a broad temperature range between 2 and 50 K (see Supplementary Note 3), which could be a sign of the presence of random singlets 62 . Furthermore, theoretically, it has been suggested that in the absence of the spin-orbit interaction, low-temperature magnetic susceptibility data of random-singlet states show a power-law χ ∝ T −γ behavior with γ < 1 24 . The susceptibility data were fitted with a power-law χ ∝ T −γ in the low-temperature region, which yielded γ = 0.82 ± 0.01 (Fig. 2b). A similar power-law fit was also observed in the quasi-twodimensional magnet Sr 2 CuTe x W 1−x O 6 , which is a promising candidate for realizing a random-singlet state on a square lattice 25 . A power-law fit of low-temperature magnetic susceptibility thus suggests the realization of a random-singlet state attributed to the presence of isolated magnetic moments in the host lattice 25,32 . The data collapse is a remarkable feature found also in other QSL candidates with a quenched disorder such as the honeycomb lattice H 3 LiIr 2 O 6 27 and the triangular-lattice Y 2 CuTiO 6 62 . The data collapse behavior in the low-temperature region suggests the presence of random singlets in LCTO induced by quenched disorder as predicted theoretically 32 .
Specific heat. Figure 3a depicts the temperature dependence of specific heat in zero magnetic field. The absence of any anomaly in the entire temperature range of investigation rules out any phase transition down to 52 mK, which is much below both the CW temperature and the dominant J and J 0 exchange interactions and thus indicates the presence of strong spin frustration that promotes a dynamic ground state in LCTO. The specific heat results in several magnetic fields upto 7 T are shown in Fig. 3b. It is observed that upon increasing the magnetic field the broad maximum shifts toward higher temperatures. Though similar to a Schottky contribution, which could arise from quasi-free spins, it, in fact, describes intrinsic spin physics due to both spins that appear isolated in the high-T spin model (Fig. 1c), and due to finite-length spin-chain fragments with an odd total number of spins, as the total spin of those cannot be lower than 1/2 by the rules for quantum addition of spin. In both cases, an applied field Zeeman splits the spin-1/2 ground state introducing a local-fielddependent energy gap that results in a specific heat peak whose position is proportional to the applied field. However, all spins, both those in finite-length spin-chain fragments, as well as nominally isolated spins, in fact couple together to form a frustrated 3D random spin lattice via weaker J inter interchain exchanges, together shaping the intriguing, liquid-like low T behavior of LCTO at T ≲ J inter .
In order to gain further insights into the magnetic properties relevant to this frustrated magnet and understand the ground-state degeneracy, it is important to extract the magnetic-specific heat. Therefore, we have considered the zero-field specific heat data (C p ) as a sum of intrinsic magnetic-specific heat (C mag ) due to exchange-coupled Cu 2+ ions, and lattice specific heat (C lattice ) due to phonons. After subtracting the lattice contribution (see Supplementary Note 5), the obtained magnetic-specific heat exhibits a clear broad maximum (see Supplementary Fig. 6) around 0.4 K releasing~13% of entropy (inset in Fig. 3c), which is a typical signature of short-range spin correlations in spin liquid candidates 27,62 . Comparing this with ED results on the high-T model with no interchain coupling, which shows no such peak in zero field, and instead exhibits residual entropy at absolute zero due to unavoidable Kramers ground-state degeneracy of individual spin-chain fragments (see Supplementary Note 7), this clearly shows the dominant impact of frustrated interchain couplings J inter at low T. Furthermore, it is observed that the temperature dependence of C mag /T shows a T −0.5 power-law behavior (Fig. 3c) at low temperature, which suggest the presence of a random singlet state with strong low-lying excitations 27,63 . Next, the entropy was obtained by integrating C mag /T in the temperature range 0.052 K ≤ T ≤ 2 K and is shown in the inset of Fig. 3c, with the contribution to entropy from specific heat below at T below the lowest measured T = 52 mK, estimated as S(52 mK) = C mag (52 mK)/(0.5). The saturation entropy amounts to just 13% of expected total entropy for a spin-1/2 system, which could suggest the presence of a highly degenerate ground state, but also naturally arises just from the fact that many spin degrees of freedom inside individual spin-chain fragments already freeze out at these low T ( J; J 0 , reducing the total number of available spin degrees of freedom and making individual spin-chain fragments behave as single effective spins in a frustrated random 3D spin lattice formed by weaker interchain interactions J inter . The presence of these frustrated interactions is also seen from the shift of specific heat to lower T than predicted by theory if they were not present (see Supplementary Note 7 for a detailed comparison with ED calculations). The power-law behavior and the absence of any long-range ordering down to 52 mK in zero field data suggest that the antiferromagnet Li 4 CuTeO 6 is a quantum spin liquid candidate.
Recently, Kimchi et al. 32 suggested that quenched disorder in spin-liquid candidates can form a network of local moments in a random-singlet phase 39,64 and the low-temperature-specific heat data recorded in various magnetic fields should collapse to a universal curve characterized by a parameter γ, where γ is defined by the relation C p =T / 1 ðμ 0 HÞ γ Fð T μ 0 H Þ and F is a scaling function 32 . We scaled the magnetic-specific heat data of LCTO accordingly (see Fig. 3d) and we indeed find data collapse for γ = 0.85, which could be attributed to the presence of random singlets in the spin-liquid ground state. A similar data collapse is reported in a triangularlattice Y 2 CuTiO 6 with γ = 0.70, where 50% of Cu 2+ -sites are diluted by Ti 4+ ions that yield quenched disorder 62 .
Muon spin relaxation. In order to gain microscopic insights into the intrinsic susceptibility and the ground-state spin dynamics, we performed zero-field (ZF) and transverse-field (TF) muon spin-relaxation (μSR) measurements on polycrystalline samples. μSR is a unique local-probe technique with a broad time window to track electron spin fluctuations. Therefore, ZF-μSR is one of the best tools to gain comprehensive microscopic information on magnetic ordering and spin dynamics of frustrated quantum materials. In Fig. 4a, we present ZF-μSR spectra at a few representative temperatures showing that muon asymmetry remains more or less unchanged for all temperatures, which confirms the absence of a phase transition down to 1.55 K 65 . The time dependence of the muon asymmetry provides information concerning internal field distributions at the muon stopping site. In multidomain or polycrystalline samples, static local field of nuclear origin B loc is randomly oriented with respect to the initial muon spin direction. The corresponding local-field distribution can be approximated by a Gaussian shape. In this case, the muon relaxation function P z (t) (also known as the Gaussian Kubo-Toyabe function (KT) is of the form where σ gives the local-field distribution width. In LCTO, the presence of competing magnetic fields arising from electron spins provides dynamic muon spin relaxation in addition to relaxation due to quasi-static nuclear fields. Namely, a damping factor exp(− λ(T)t) characteristic of dynamical fields has to be added to the KT function to account for the experimental results. Therefore, the zero-field asymmetry spectra were fitted by the damped Gaussian Kubo-Toyabe relaxation function where λ(T) is the dynamical muon spin-relaxation rate due to electronic spins. The fitting of the observed ZF asymmetry data (Fig. 4a) by Eq. (3) yields the nuclear field distribution width of approximately 3.66 G and the dynamical relaxation rate λ = 0.023 μs −1 . In addition, the strong change of asymmetry spectra between zero field and a small longitudinal field of 50 G as shown in Fig. 4b also supports this model, as the applied field decouples the relaxation due to small static nuclear fields while it does not affect the relaxation due to larger but dynamical fields of electronic origin. The obtained temperature dependence of the muon spin-relaxation rate due to dynamical fields in ZF remains constant down to 1.55 K (Fig. 4c) without any signs of magnetic ordering, which suggests a dynamic ground state. We note that the presence of a disordered static local field at the muon site generally leads to the "1/3" tail in zero-field muon decay asymmetry of powder samples that is a hallmark of static magnetism. However, zero-field muon asymmetries do not show the "1/3" tail, which confirms the absence of spin-freezing 66 . In LCTO, the absence of any oscillation in ZF spectra (see Fig. 4a) also corroborates the absence of long-range magnetic ordering down to 1.55 K. Next, we carried out muon rotation measurements in the presence of a large transverse field (TF) B T = 4900 G, which is much greater than the local fields. The observed experimental TF spectra were fitted using an oscillatory Gaussian decaying function Here ϕ~π/3 is the initial phase of the oscillation, ν(T) is the frequency of muon precession originating from the local magnetic field and external transverse magnetic field, whereas λ T (T) is the transverse muon spin-relaxation rate. In the presence of an external transverse field (B T ), the total magnetic field experienced by muon can be defined as the sum of B T and B loc . So, the mean precession frequency can be defined as ν(T) = γ μ (B T + B loc (T))/2π, where only the local magnetic field B loc (T) changes with temperature (γ μ = 2π × 135.53 MHzT −1 is the muon gyromagnetic ratio). The muon Knight shift (K) that is directly proportional to magnetic susceptibility (χ) via K = (ν(T) − ν 0 )/ν 0 = Aχ, where ν 0 = (γ μ /2π) B T is the reference frequency and A is the coupling constant between the muon magnetic moment and the electron magnetic moments, thus captures intrinsic susceptibility relevant to the spin lattice. Inset of Fig. 4d depicts the muon Knight shift as a function of bulk  Fig. 3 Specific heat measurements reveal a randomness-induced spin-liquid state in Li 4 CuTeO 6 . a The temperature dependence of the specific heat (C p ) in zero magnetic field in a broad temperature range 0.052 K ≤ T ≤ 250 K. b The temperature-dependent specific heat in different magnetic fields upto 7 T. c The temperature dependence of magnetic-specific heat divided by temperature in zero field in a log-log scale where the red line is the power-law fit. Inset shows the temperature dependence of associated entropy change. d The scaling behavior of the magnetic-specific heat with C mag ðμ 0 HÞ 0:85 =T as a function of the scaled temperature T/μ 0 H for several fields.
magnetic susceptibility (χ) with temperature as an implicit parameter. It is seen that K as a function of χ exhibits two linear regions that are consistent with the observed crossover region in the temperature dependence of λ T . Below T < 75 K, the linear fit (orange line) yields a coupling constant A = −114 Oeμ À1 B whereas we find A = −211 Oeμ À1 B from the linear fit (olive line) above 75 K. Furthermore, it is observed that the Knight-shift data follows K~T −0.75 power-law behavior at low T (see Fig. 4d), which is consistent with bulk magnetic susceptibility and confirms that such susceptibility is intrinsic to the system. This behavior further corroborates the presence of "random-singlets" in LCTO. μSR results do not indicate any static local fields of electronic origin down to 1.55 K. The measured muon asymmetry decay in zero field is, in fact, dominated by static nuclear fields, while relaxation due to dynamical fields of electronic origin remains hindered behind the nuclear relaxation due to exchange narrowing down to 1.55 K. This reveals spin fluctuations are fast even at T/θ CW ≈ 0.01, which again suggests a dynamic ground state in LCTO. The temperature dependence of λ T (T) is shown in Fig. 4e. Upon decreasing the temperature, the relaxation rate increases and a rather strong change of λ T is observed in the temperature range 1.55 K < T < 10 K, which suggests the presence of a crossover region 12 . Above T > 20 K, λ T varies very slowly with temperature as would be expected for a paramagnetic region. We note that the zero-field muon spin-relaxation rate due to electronic fields (Fig. 4c) is much lower than the transverse muon spin-relaxation rate (Fig. 4e) and its temperature dependence is different. As the former relaxation is only due to field fluctuations, while the latter is affected also by static localfield distributions induced by the applied transverse field, it is the induced local fields that are responsible for the increase of the transverse relaxation rate with decreasing temperature. This increase is characterized by a power-law T −0.2 at the lowest temperatures (Fig. 4e).

Discussion
Our investigation reveals the presence of a strong antiferromagnetic exchange interactions between spins of Cu 2+ ions constituting a random frustrated spin-chain model with weaker interchain exchanges couplings forming an effective frustrated random 3D lattice of spin-chain fragments at low T. Remarkably, the material does not show any signature of a phase transition or spin-glass freezing down to at least 45 mK, which suggests a dynamical liquid-like ground state in this antiferromagnet. The scaling behavior of magnetization and specific heat data at low T reflects the presence of "random-singlets" in a quantum disordered ground state of LCTO. Magnetic materials with a dilute random network of S = 1/2 sites can be considered as a combination of two sub-systems where the majority of constituent spins develop spin-liquid or valence bond crystal, and the minority forms random singlets, and leaves out a random network of S = 1/2 sites. For such a system, in the ground state, spins do not undergo any long-range antiferromagnetic order. The spins that are strongly coupled rather constitute singlet pairs and the remaining spins are almost free due to the development of effective antiferromagnetic exchange interaction (J) with a powerlaw probability distribution, [P](J)~J −γ , between them in the low-T renormalization group flow of the random-singlet phase, where γ is a positive quantity related to the spin correlation length 22 . The small sub-system of spin-1/2 moments in the lowenergy limit leads to unconventional scaling behavior in specific heat C[T]~T 1−γ and spin susceptibility χ[T]~T −γ32 . In the presence of an applied magnetic field, specific heat captures the modified distribution of energy through a field-dependent coefficient, i.e., C $ T=ðμ 0 HÞ γ at low temperatures. It is worth to mention that the temperature and field dependence of magneticspecific heat displays T/μ 0 H data collapse for γ = 0.85 in LCTO, which is a bit higher than the observed γ = 0.5 for the honeycomb lattice H 3 LiIr 2 O 6 and triangular-lattice LiZn 2 Mo 3 O 8 but it is close to the observed γ = 0.7 for triangular-lattice Y 2 CuTiO 6 . The absence of strong spin-orbit interaction in Y 2 CuTiO 6 and LCTO possibly plays a vital role resulting in data collapse with higher value of γ compared to H 3 LiIr 2 O 6 and Li 2 Zn 2 Mo 3 O 8 . In addition, we also found a divergence in the low-temperature susceptibility of LCTO, which scales as T −γ with γ = 0.82 for bulk measurements and γ = 0.75 for local-probe muon Knight-shift measurements. This value of γ is comparable with the reported γ = 0.68 and 0.76 for Y 2 CuTiO 6 and Sr 2 CuTe 1−x W x O 6 , respectively 24,25,62 . Further theoretical and experimental studies are desired to shed additional insights into the magnetism and spin dynamics of this promising spin-liquid candidate and to clarify the degree to which the random-singlet picture might be modified due to the presence of an underlying spin-chain fragment structure with random connectivity. The chain fragments could conceivably contribute further ingredients to the low-temperature physics of LCTO, leading to novel phenomena outside the paradigm of the usually considered random-singlet states.
Conclusion. We have investigated the structure, thermodynamic and local magnetic properties of Li 4 CuTeO 6 by employing X-ray, neutron diffraction, magnetization, specific heat, ESR, and muon spin-relaxation measurements. These were supported by state-ofthe-art ab initio DFT and ED calculations to elucidate the high-T spin model and its thermodynamics. Li 4 CuTeO 6 crystallizes in the monoclinic space group C2/m, wherein the Cu 2+ site constitute spin chains with exchange J along the [101] direction of random length due to partial occupancy of the Cu 2 site, with random additional Cu 2+ ions on partially occupied Cu 3 sites surrounding the main chains at positions displaced along ± b from the main-chain bond centers and coupling to two Cu 2 sites on the main chain through an exchange J 0 . The large and negative CW temperature, θ CW = −154 K, suggests the strong antiferromagnetic nature of the J and J 0 exchange interactions between Cu 2+ moments, which is confirmed by DFT and fits of ED calculations to experimental data for T > 20 K. These gradual deviations from these fits at low T reveals further frustrating 3D interactions between the random spin-chain fragments, which dominate the low T response. Specific heat experiments reveal the absence of long-range magnetic ordering down to 52 mK despite the presence of a large CW temperature, implying a highly frustrated spin lattice. The absence of spin-freezing despite the Li/Cu anti-site disorder reflects that Cu 2+ moments remain fluctuating down to 45 mK. The zero-field and transverse-field μSR measurements support a dynamic ground state down to 1.55 K. Moreover, the existence of T/μ 0 H scaling as well as magnetization and specific heat data collapse indicate the presence of a spin-1/2 network of random singlets due to unavoidable anti-site disorder in the host lattice and frustrated weak 3D couplings. Our results demonstrate a randomness-induced spin-liquid state in a frustrated magnet. These results offer an exciting ground to explore novel quantum states with unconventional low-energy excitations in novel frustrated quantum materials with quenched disorder.

Methods
Polycrystalline samples of Li 4 CuTeO 6 were synthesized by the conventional solidstate reaction route as outlined in ref. 53 . To obtain a single phase of LCTO, stoichiometric amounts of Li 2 CO 3 (Alfa Aesar, 99.0%), CuO (Alfa Aesar, 99.995%), and TeO 2 (Alfa Aesar, 99.9995%) were mixed and the mixture was finely ground and pressed into a pellet, and after several intermediate steps, finally sintered at 850°C for 30 h. The X-ray diffraction measurements were performed at 300 K by employing a Rigaku smart LAB X-ray diffractometer with Cu Kα radiation (λ = 1.54 Å). Due to the light mass of Li and O atoms, the analysis of X-ray powder diffraction data is not sufficient to determine the structural disorder and the mixed site occupancy. Therefore, neutron-diffraction measurements were performed at room temperature for the incident wavelength of 1.28 Å using crg-D1B Two-Axis Powder Diffractometer at the Institut of Laue-Langevin, France. Our neutrondiffraction studies detects neither spurious phases nor non-stoichiometry of oxygen. The neutron-diffraction results at 1.28 Å are shown in Fig. 1a, while the neutron-diffraction pattern, XRD pattern, and resulting refinements parameters are presented in Supplementary Notes 1 and 2. This suggests that the polycrystalline samples used in this study are of very high quality. The refined occupancy of all the atomic sites are within the experimental uncertainty of about 1%, thus ruling out the presence of non-stoichiometric oxygen.
Magnetization measurements were carried out using a Quantum Design, SQUID (MPMS) in the temperature range 1.9 K ≤ T ≤ 300 K under magnetic fields upto 5 T. The electron spin resonance (ESR) spectrum was measured at 250 K at the irradiation frequency of 256.3 GHz on a custom-built spectrometer with homodyne detection at the NHMFL, Tallahassee, USA. Thermoremanent measurements (see Supplementary Note 4) were performed using a Quantum Design, SQUID (MPMS) and the sample was cooled down from 100 to 2 K in either 1000 or 8 Oe, and turned down at 2 K. After 1 h, the signal was measured again (at t = 0). These results rule out the presence of a spin-glass transition and are shown in Supplementary Fig. 4.
Specific heat measurements were performed in a Quantum Design, Physical Properties Measurement System (PPMS) by the thermal relaxation method, in the temperature range 1.9 K ≤ T ≤ 240 K and in magnetic fields upto 7 T. Furthermore, specific heat measurements were carried out separately in the temperature range 0.052 K ≤ T ≤ 4 K in zero field and 3 T using a dilution refrigerator which was also used to measure ac susceptibility in the temperature range 0.045 K ≤ T ≤ 4 K at four different frequencies using a Dynacool PPMS instrument from Quantum Design, San Diego, USA.
μSR measurements were performed using the GPS spectrometer at the Paul Scherrer Institute, Villigen, Switzerland, on a polycrystalline sample with mass of 1 g down to 1.55 K in zero field and in a transverse field of 4900 G. DFT+U calculations on a 2 × 1 × 2 supercell of LCTO were carried out using the plane-wave DFT code 67 using a local (spin) density approximation (LSDA) exchange-correlation functional in a LSDA+U scheme 68 with an effective Hubbard U eff = 9 eV, a 2000 plane-wave energy cutoff, and a 2 × 3 × 2 Monkhorst-Pack grid reciprocal-space sampling 69 . To extract exchanges total energies of ≥120 random collinear spin configurations were calculated and fitted to a symmetry-adapted spin model 70 . See Supplementary Note 6 for further details. ED calculations with upto N = 18 spins were carried out on > 13,000 distinct configurations of occupied Cu 2 and Cu 3 sites around finite-length spin-chain fragments in the high-T random 1D model of Fig. 1c, achieving a total probability between 60.2% and more than 98.8%, depending on parameters, that a randomly-chosen Cu 2+ spin in LCTO is part of one of the ED calculated spin-chain fragments. Details of these ED calculations are presented in Supplementary Note 7.

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request. The neutron-diffraction data can be found at https://doi. org/10.5291/ILL-DATA.EASY-1043.