## Abstract

We report a comprehensive muon spectroscopy study of the Zn-barlowite series of \(S=\frac{1}{2}\) kagomé antiferromagnets, Zn_{x}Cu_{4−x}(OH)_{6}FBr, for *x* = 0.00 to 0.99(1). By combining muon spin relaxation and rotation measurements with state-of-the-art density-functional theory muon-site calculations, we observe the formation of both *μ*–F and *μ*–OH complexes in Zn-barlowite. From these stopping sites, implanted muon spins reveal the suppression of long-range magnetic order into a possible quantum spin liquid state upon the increasing concentration of Zn-substitution. In the parent compound (*x* = 0), static long-range magnetic order below *T*_{N} = 15 K manifests itself in the form of spontaneous oscillations in the time-dependent muon asymmetry signal consistent with the dipolar fields expected from the calculated muon stopping sites and the previously determined magnetic structure of barlowite. Meanwhile, in the *x* = 1.0 end-member of the series—in which antiferromagnetic kagomé layers of Cu^{2+}\(S=\frac{1}{2}\) moments are decoupled by diamagnetic Zn^{2+} ions—we observe that dynamic magnetic moment fluctuations persist down to at least 50 mK, indicative of a quantum disordered ground state. We demonstrate that this crossover from a static to dynamic magnetic ground state occurs for compositions of Zn-barlowite with *x* > 0.5, which bears resemblance to the dynamical behaviour of the widely studied Zn-paratacamite series that contains the quantum spin liquid candidate herbertsmithite.

## Introduction

Quantum effects play a significant role in the low-temperature physics of magnetic systems in which antiferromagnetic \(S=\frac{1}{2}\) moments decorate a two-dimensional kagomé array of corner-sharing triangles^{1}. The competing interactions resulting from the geometric frustration in such systems combined with quantum fluctuations may give rise to elusive states of matter, such as quantum spin liquids (QSLs); a possibility that continues to capture the interest of the quantum materials research community^{2}. Indeed, QSLs represent a unique state of matter that evades long-range magnetic order despite the sometimes considerable exchange interactions between its magnetic moments, which would typically drive a symmetry-breaking magnetic phase transition in a conventional magnet following Landau theory^{3}. Instead, QSLs are quantum superpositional states formed from entangled pairs of magnetic moments, or spins, that can, in certain cases, be envisioned within a resonating valence bond model^{4}. QSLs are of fundamental importance as their highly entangled nature should give rise to exotic physical phenomena, such as topological phases, fractionalised excitations and emergent gauge fields^{5}. However, the unambiguous realisation of a QSL remains an outstanding challenge.

In particular, while from a theoretical perspective it is generally accepted that the ground state of the \(S=\frac{1}{2}\) kagomé antiferromagnet is a QSL^{6}, the true nature of this putative QSL state is still under debate, and new material realisations of it are needed to help resolve the ongoing disconnect between theory and experiment^{7,8,9}. In this regard, synthetic analogues of Cu(II) hydroxyl halide minerals have been central to the experimental exploration into the physics of the \(S=\frac{1}{2}\) kagomé antiferromagnet in recent years^{10,11,12}. To date, the most widely studied example of this class of materials is the Zn-paratacamite series, Zn_{x}Cu_{4−x}(OH)_{6}Cl_{2}, of which the *x* = 1 end-member is the QSL candidate herbertsmithite, ZnCu_{3}(OH)_{6}Cl_{2}. The crystal structure of herbertsmithite contains kagomé layers of Cu^{2+}\(S=\frac{1}{2}\) moments separated by diamagnetic Zn^{2+} ions^{13,14,15}, which give rise to a dynamic, disordered magnetic ground state with magnetic moment fluctuations persisting to the lowest measurable temperatures despite strong nearest-neighbour antiferromagnetic exchange (*J* = 200 K)^{14,16}. Precisely characterising the instrinsic physics of the \(S=\frac{1}{2}\) kagomé layers of herbertsmithite, however, is extremely challenging due to the presence of Cu^{2+}/Zn^{2+} site occupancy disorder within its crystal structure^{15,17,18}.

A closely related Cu(II) hydroxyl halide that has garnered considerable attention in the recent literature is barlowite, Cu_{4}(OH)_{6}FBr^{19}, a naturally occurring mineral which has been synthetically produced by several groups^{20,21,22,23,24,25}. A recent powder neutron diffraction (PND) study contributed to the ongoing debate around the nuclear and magnetic structures of this material^{23}: like herbertsmithite, barlowite is formed from Cu^{2+}\(S=\frac{1}{2}\) kagomé layers but with a second interlayer Cu^{2+} site, which is disordered at room-temperature in the widely reported *P*6_{3}/*m**m**c* model^{20,21,24,25}. Below ~250 K, one reported class of sample undergoes a subtle structural distortion to an orthorhombic *P**n**m**a* phase^{23} (see Fig. 1a), which appears to relieve magnetic frustration and allows for the onset of magnetic order at *T*_{N} = 15 K^{26}. Intriguingly, computational studies have suggested that barlowite may be tuned towards a QSL phase by substituting the interlayer Cu^{2+} ions with diamagnetic Zn^{2+} or Mg^{2+}^{27,28} thus reducing the magnetic coupling between the frustrated \(S=\frac{1}{2}\) kagomé layers. Furthermore, the calculated energy of formation for cation site occupancy defects is an order of magnitude greater in substituted barlowite than for herbertsmithite^{27}, suggesting Mg-barlowite or Zn-barlowite may be highly promising candidates in which to explore the intrinsic properties of a Cu^{2+}-based \(S=\frac{1}{2}\) kagomé antiferromagnet. Indeed, the end-member of the Zn-barlowite series, nominally ZnCu_{3}(OH)_{6}FBr, has been studied using inelastic neutron scattering and ^{19}F NMR and seems to display evidence for a gapped Z_{2} QSL^{29,30}, a theoretically predicted ground state for the \(S=\frac{1}{2}\) Heisenberg model on an antiferromagnetic kagomé net^{31}.

Despite this promise, the evolution of the magnetic ground state and characteristic magnetic moment correlation dynamics across the Zn-barlowite series are yet to be explored in any great detail using local probe techniques, such as muon spin relaxation and rotation (*μ*SR). In these techniques, spin-polarised, positively charged muons are implanted into a sample of interest, in which they thermalise in regions of high electron density. From there, the spin polarisation of implanted muons will evolve over time, depending on the presence of local magnetic fields at the muon stopping site, which may originate from both nuclear and electronic moments within the sample. Following the time-dependence of the muon spin polarisation—via the asymmetry of the emitted muon-decay positrons—thus provides a unique insight into the local magnetic properties of condensed matter systems^{32}. As the muon spin is sensitive to very small magnetic fields (~10^{−5} T) and fluctuating field dynamics on a microsecond timescale, *μ*SR has proved to be invaluable in revealing dynamical magnetic correlations in herbertsmithite and many other frustrated quantum magnets^{10,16,33,34,35,36,37}. However, limitations in the interpretation of *μ*SR data often stem from the lack of knowledge of the possible muon stopping sites within complex crystal structures and the extent to which the implanted muons distort the surrounding lattice^{38}. Crucially, the development of density-functional theory (DFT) to calculate muon stopping sites in crystalline solids is increasingly enabling the determination of these factors to advance the quantitative analysis of experimental *μ*SR data^{39}. Here, we present such a study—combining comprehensive *μ*SR measurements on the Zn-barlowite series with supporting DFT muon-site calculations—to reveal the onset of the possible QSL phase within this family of quantum materials. In this case, we take advantage of the formation of both *μ* − F and *μ* − OH complexes upon muon implantation in Zn-barlowite which—in combination with our DFT calculations—gives a strong constraint for the accurate determination of muon stopping sites and the local magnetic properties of this series.

## Results and discussion

### Combining zero-field *μ*SR measurement and DFT to identify muon stopping sites in barlowite

We first consider the parent compound barlowite, Cu_{4}(OH)_{6}FBr, which—in accordance with previous studies^{20,21,22,24}—orders at *T*_{N} = 15 K as indicated by the sharp upturn in the magnetic susceptibility shown in Fig. 2a. Curie-Weiss fitting the high-temperature inverse susceptibility (shown in Supplementary Fig. 1) yields an antiferromagnetic Weiss temperature, *θ* = −111.4(3) K, which is manifested in the antiferromagnetic alignment of Cu^{2+} moments within the kagomé planes of barlowite in its magnetic structure previously determined by neutron diffraction^{23}. This magnetic structure includes a ferromagnetic component due to the canting of the magnetic moments at the interlayer Cu^{2+} site. Figure 3a shows the muon decay asymmetry measured in zero-field (ZF) on the MuSR spectrometer for such a sample. Above *T*_{N}, the long-time oscillations observed in the data are characteristic of the local nuclear magnetic fields at the muon stopping sites within the sample. In particular, such oscillations are indicative of dipolar coupling between the muon and nuclear spins that leads to the formation of entangled state complexes. DFT muon-site calculations reveal two distinct classes of muon stopping site complexes in barlowite, which we show are consistent with the experimental ZF asymmetry signal shown in Fig. 3a. The first class of muon stopping site localises ~1.0 Å away from the oxygen atoms in the hydroxide groups that connect the Cu^{2+} ions within the kagomé layers of barlowite. This forms a triangular *μ*–OH complex, with the muon-proton distance found to be 1.54 Å in the lowest-energy sites, as shown in Fig. 1c. In the second class of muon stopping site, muons are found to localise near the fluoride anions in between the kagomé layers of barlowite, with a *μ*–F separation of 1.1 Å, as depicted in Fig. 1d. Such *μ*–F sites lie at substantially higher energies above the lowest energy *μ*–OH sites in our calculations (≈1 eV). Although this suggests—on purely energetic grounds—that the formation of *μ*–F complexes in barlowite is unlikely compared to *μ*–OH, it is possible that during the stopping process, muons can be captured in these local potential minima.

Based on the results of these DFT calculations, we have devised a model that describes the ZF asymmetry of barlowite in the paramagnetic regime, incorporating the two distinct classes of muon stopping sites,

*A*_{bg} accounts for the background contribution of muons that stop outside the sample, which is weakly relaxing, and the expressions *P*_{F}(*t*) and *P*_{OH}(*t*) describe the dipolar interactions between the nuclear magnetic moments of fluoride and hydroxide groups and the muon spin at each class of stopping site within the sample, given by^{40},

The oscillation frequencies *ω*_{F} and *ω*_{OH} are related to the respective *μ*–F and *μ*–H distances, *d*_{F(H)}, through,

where *μ*_{0} is the permeability of free space and *γ*_{μ} = 2*π* × 135.5 MHz T^{−1}, *γ*_{F} = 2*π* × 40.1 MHz T^{−1} and *γ*_{H} = 2*π* × 42.6 MHz T^{−1} are the gyromagnetic ratios of the muon spin and the fluorine and hydrogen nuclear spins, respectively. A fit of this model to ZF data collected for barlowite is shown in Fig. 3a with the fitting parameters given in Table 1. The experimentally determined *μ*–F and *μ*–H distances are in good agreement with those calculated from DFT. The Gaussian damping term in Equation 1 phenomenologically describes a possible distribution of nuclear fields surrounding the muon spin, with Δ_{F} and Δ_{OH} describing the field distributions at each class of stopping site. Indeed, DFT calculations also suggest that there are many possible stopping sites within the *μ*–OH class, separated in total energy by <0.1 eV. This reflects the possibility that muons localise near the oxygen atoms in a number of similar configurations. There is some variation in the *μ*–H distances calculated for these similar stopping sites, with 1.63 Å found in a few of the higher-energy sites in this class. In all calculated cases of *μ*–OH complex, however, the muon causes a distortion to the local crystal structure in the vicinity of its stopping site. In the lowest-energy sites of this class, the muon pushes the closest Cu^{2+} ion away from the muon site by ~0.45 Å, as demonstrated in Fig. 1c and d. In contrast, little distortion of the local environment occurs close to the *μ*–F sites.

### Resolving magnetic order in barlowite through high-resolution *μ*SR

Below *T*_{N}, the internal electronic magnetic fields arising from the ordered Cu^{2+} moments in barlowite rapidly depolarise the implanted muon spins. As a consequence, much of the initial asymmetry in our MuSR datasets is lost, as can be seen in Fig. 3a, owing to the limited time-resolution of the pulsed muon beam at the ISIS Neutron and Muon Source. We do, however, observe a ‘one-third tail’ in the baseline asymmetry, a feature attributed to the presence of static magnetic moments in a polycrystalline sample, where on average, one-third of the implanted muon spins will be initially parallel to the local magnetic field, and whose polarisation will thus remain unaffected in the absence of dynamics^{41}. To overcome this limitation, we turn to data collected for barlowite on the GPS instrument at PSI, shown in Fig. 3b and Supplementary Fig. 4, for which the continuous nature of the muon source provides sufficient time-resolution to capture the spontaneous oscillations apparent below 0.25 μs in the ZF asymmetry signal. These oscillations are indicative of the muon spin precession around the static internal magnetic fields of electronic origin at a frequency given by \({\nu }_{i}=({\gamma }_{\mu }/2\pi ){B}_{i}^{{\mathsf{loc}}}\), where \({B}_{i}^{{\mathsf{loc}}}\) is the local magnetic field at each muon stopping site, *i*, arising from the neighbouring ordered Cu^{2+} moments. Consistent with the DFT calculations and the analysis of the ZF asymmetry data collected above *T*_{N}, the spontaneous oscillations in the ZF asymmetry data collected below *T*_{N} in Fig. 3b can also be modelled in terms of two distinct classes of muon stopping site with the following expression,

The two dominant oscillating frequencies, *ν*_{1} = 6.3(2) MHz and *ν*_{2} = 13.7(2) MHz, extracted from this fit can also be clearly resolved in the Fourier transform of the GPS data collected at 4.2 K, as shown in Fig. 3c. Comparing the fitted amplitudes of each component, *a*_{1} = 64(2)% and *a*_{2} = 34(2)%, with the stopping site fractions, *f*_{F} and *f*_{OH}, obtained from fitting the ZF data above *T*_{N} in Fig. 3a (see Table 1), we assign the former to the *μ*–F and the latter to the *μ*–OH classes of muon stopping site. From this assignment, the data imply that the local magnetic fields at each class of muon stopping site arising from neighbouring ordered Cu^{2+} moments are thus \({B}_{{\mathsf{F}}}^{{\mathsf{loc}}}\,\approx \,46\) mT and \({B}_{{\mathsf{OH}}}^{{\mathsf{loc}}}\,\approx \,101\) mT, respectively. It should be noted that the weak relaxation of the one-third tail in the asymmetry data should not be mistaken for residual spin dynamics in barlowite below *T*_{N}. In the case of Fig. 3b, it is due to the small fraction of muons that stop in the aluminium foil used to contain the sample, whose spins thus experience the nuclear fields of ^{27}Al.

As a final consistency check for our combination of DFT methods and muon experiment, we estimate the local electronic magnetic moment dipolar fields at each of the calculated muon stopping site classes based on the known magnetic structure for barlowite below *T*_{N}^{23}. For the lowest-energy *μ*–OH stopping sites, and assuming an undistorted crystal structure, these calculations suggest typical local fields around 120 mT, along with a few around 75 mT. If the muon-induced distortion is taken into account we obtain still smaller fields on the order of 50 mT, owing to the movement of the closest Cu^{2+} moment. The \({B}_{{\mathsf{OH}}}^{{\mathsf{loc}}}\) fields obtained are, therefore, highly sensitive to this distortion, but it is probable that this effect is over estimated in our computations owing to the relatively small size of the supercell used. The average \({B}_{{\mathsf{F}}}^{{\mathsf{loc}}}\) for the *μ*–F class of stopping sites is smaller, at typically 55 mT, and these values are insensitive to the inclusion of the distortions. Overall, the agreement of these calculations with the fields determined from our PSI data (summarised in Table 1) reinforces the likely formation of *μ*–F stopping sites in barlowite despite their overall higher energy compared with *μ*–OH.

### Following the evolution towards quantum disorder in Zn-barlowite through transverse-field (TF) *μ*SR

We next consider the effect of substituting Cu\({}^{2+}\,S=\frac{1}{2}\) ions in barlowite for diamagnetic Zn^{2+} in the series, Zn_{x}Cu_{4−x}(OH)_{6}FBr. Our analysis of PND data collected for Zn-substituted samples, shown in Supplementary Fig. 2 and Supplementary Table 3, yields two important observations. First, for highly substituted samples (*x* > 0.5) the crystal structure retains hexagonal *P*6_{3}/*m**m**c* symmetry down to 1.5 K—the base temperature of our diffraction experiments—in which the kagomé layers remain undistorted. Second, refinement of the cation site occupancies within this hexagonal structure indicate that the substituted Zn^{2+} ions predominately occupy the trigonal prismatic interplane site, with only 0.6(2) % of the intralayer sites occupied by Zn^{2+}, an order of magnitude less than reported for herberstmithite^{17}. This is in keeping with theory^{28} and suggests that Zn-substitution effectively reduces the magnetic coupling between the kagomé layers of Cu^{2+} moments in barlowite. Correspondingly, we observe a progressive suppression of the magnetic ordering transition in the temperature-dependent magnetic susceptibility for Zn_{x}Cu_{4−x}(OH)_{6}FBr with increasing *x*, as shown in Fig. 2a. Figure 2b shows the fitting parameters obtained from the Curie-Weiss analysis of these data, revealing a reduction in the effective magnetic moment, *μ*_{eff}, per formula unit across the series as Cu^{2+} is replaced with Zn^{2+}, and an increase in the absolute value of the Weiss constant, *θ*, since the net magnetic exchange becomes increasingly negative as the positive ferromagnetic contribution from the interlayer coupling is reduced across the series.

In order to follow the evolution of the magnetic ground state in Zn-barlowite, we performed a set of TF *μ*SR measurements across the series, in which a magnetic field of 2 mT is applied perpendicular to the initial muon spin polarisation. In this TF geometry for a sample in the paramagnetic regime, implanted muon spins precess around the applied field with a frequency, *ν*_{TF}. Below *T*_{N}, however, the local internal magnetic fields within the magnetically ordered state can dominate the muon spin relaxation, resulting in a loss of asymmetry and a dephasing of the oscillating signal measured in the TF experiment. Both of these effects can be seen in the TF data collected for an *x* = 0.00 sample above and below *T*_{N} shown in Fig. 4a, and can be modelled by,

where *A*_{p} provides a measure of the paramagnetic (dynamic) volume fraction of the sample, *A*_{bg} gives the background contribution of muons that stop in the silver sample holder and *B* accounts for the fraction of muons with their spin polarisation aligned with the local magnetic fields in the ordered state of our polycrystalline samples. Therefore, in the *x* = 0.00 sample below *T*_{N}, *B* tends to one-third of the asymmetry we observe in the paramagnetic state, in a similar manner to the one-third tail observed in our ZF data discussed above. By fitting this model to TF data collected over a range of temperatures, we may plot the temperature dependence of *A*_{p} or, in this case, the normalised inverse of this value (i.e. 1 − *A*_{p}(*T*)/*A*_{p}(20 K)), which represents the evolution of the frozen volume fraction of static magnetic moments within a sample. For example, for the *x* = 0.00 sample, a sudden drop in *A*_{p} at *T*_{N} = 15 K corresponds to a dramatic increase in the static volume fraction, shown in Fig. 4b, a clear signal of the magnetic ordering transition, in agreement with our magnetic susceptibility (Fig. 2a) and PND measurements^{23}.

The same analysis can be applied to samples in the Zn-barlowite series upon increasing Zn^{2+} content, and Fig. 4b shows the temperature dependence of the frozen volume fraction for Zn_{x}Cu_{4−x}(OH)_{6}FBr with *x* = 0.159(3) and 0.32(1). For *x* = 0.159(3) we observe an extended magnetic ordering regime, while for the *x* = 0.32(1) sample, the static moment fraction at base temperature is vastly reduced compared with the parent *x* = 0.00 compound. Moreover, there is also evidence to support the coexistence of static and dynamic magnetic correlations in the *x* = 0.159(3) and 0.32(1) members of the Zn-barlowite series, as in both ZF and TF measurements of these samples we do not recover the one-third asymmetry expected for a fully ordered magnetic ground state in a polycrystalline sample, as shown in Supplementary Fig. 3. As the Zn^{2+} content of Zn_{x}Cu_{4−x}(OH)_{6}FBr is increased further to *x* = 0.52(1), dynamic magnetic correlations dominate the ZF muon asymmetry at low-temperature—likely indicating a static to dynamic crossover region in the magnetic phase diagram—and by *x* = 0.66(1), the magnetic moments of the Cu^{2+}\(S=\frac{1}{2}\) ions remain dynamically fluctuating down to the lowest temperatures of our *μ*SR experiment. This is most clearly evidenced by the low-temperature ZF data collected for the *x* = 0.66(1) sample shown in Fig. 5a, which can be modelled using the same expression developed for the *x* = 0.00 compound in the paramagnetic regime (Eq. 1 and Supplementary Table 4). An important result from our DFT calculations in this regard is that the inclusion of Zn^{2+} ions and the hexagonal space group symmetry of the Zn-substituted samples does not yield any notably different classes of low-energy muon stopping sites compared with the parent orthorhombic structure. Together, our observations from ZF and TF *μ*SR measurements from across the Zn-barlowite series allow us to map out the low-temperature magnetic phase diagram shown in Fig. 4c, which reveals a critical reduction in the static volume fraction from below *x* = 0.5 as the system enters a possible QSL state. This is strikingly similar to the magnetic phase diagram known for the Zn-paratacamite series^{16}, in which the static volume fraction drops to zero near *x* = 0.6, and is reproduced well by a recent study of the Zn-barlowite phase diagram by Wei et al.^{42}.

### Probing spin fluctuations in Zn-barlowite through longitudinal-field *μ*SR

Finally, we turn to explore the dynamical correlations within the quantum disordered ground state of Zn_{x}Cu_{4−x}(OH)_{6}FBr for *x* ≥ 0.66. In order to isolate and probe the dynamics of the electronic magnetic moment fluctuations in these highly substituted samples, we applied a longitudinal field (LF) of 10 mT in our *μ*SR experiments. In doing so, the implanted muon spins are decoupled from the nuclear magnetic moments that dominate the ZF signal, meaning any remaining depolarisation of the muon spin arises from its interaction with the local electronic magnetic fields (see Supplementary Fig. 5). Figure 5a shows representative LF data for the *x* = 0.66(1) sample, which can be modelled at all temperatures using a simple exponential relaxation function,

where *λ* is the muon spin relaxation rate, *A*_{r} is the relaxing amplitude and *A*_{bg} is the background contribution from muons stopping outside the sample. Such an exponential relaxation of the LF signal of a sample is characteristic of electronic magnetic moment fluctuations within a motional narrowing regime, with the corresponding relaxation rate described by the Redfield equation^{43}:

where *γ*_{μ} is the muon spin gyromagnetic ratio, *B*_{LF} is the applied longitudinal-field and Δ and *τ* are the internal electronic magnetic field distribution and its fluctuation time—characteristic to a particular sample—respectively. A common feature that emerges from the analysis of LF *μ*SR data of many QSL candidates of diverse chemical and structural varieties is the presence of a plateau in the muon spin relaxation rate, *λ*, at low-temperatures^{16,37,44,45,46,47}.

Qualitatively, this can be taken to indicate the presence of a QSL ground state characterised by dynamic electronic magnetic moment fluctuations with a temperature-independent fluctuation time, which is unlike any system with a static magnetic ground state, for example, a spin glass, which is associated with a maximum in *λ* near its spin freezing temperature^{48}. Quantitatively, however, it is challenging to meaningfully compare muon spin relaxation rates observed in the plateau states across different classes of QSL candidate materials because, as Eq. 7 highlights, there are several key parameters that will determine the overall magnitude of *λ* for any given system, as well as experimental factors, such as the applied LF and signal background. Even across a single family of materials—here Zn_{x}Cu_{4−x}(OH)_{6}FBr—substituting Cu^{2+} for diamagnetic Zn^{2+} will affect the size and distribution of the internal electronic magnetic field, which can have pronounced consequences for the magnitude of *λ*. For instance, in this particular case, one might expect that as *x* increases across the Zn_{x}Cu_{4−x}(OH)_{6}FBr series at a given temperature, Δ will decrease as magnetic Cu^{2+} is replaced by diamagnetic Zn^{2+}, but that *τ* will increase as the system approaches a dynamic spin liquid state. Together, this implies that the magnitude of *λ* will decrease with increasing concentration of Zn-substitution, which has been observed in the Zn-paratacamite family of herbertsmithite^{16}, but does not appear to be the case for samples of Zn-barlowite with *x* > 0.5 (see Supplementary Fig. 6). However, another compelling viewpoint to consider in this regard is that at low-temperatures, muon spins implanted within highly substituted samples of Zn-barlowite may be predominantly influenced by their dipolar interaction with defect spins of Cu^{2+} moments remaining at the interlayer sites^{49}. If this is the case, then one would expect a decrease in the dominant dipolar fields coupling to the muon spin as the concentration of defect Cu^{2+} moments decreases with increasing Zn-subsitution. Such an argument could thus account for the observed variation in the magnitude of *λ* in Zn-barlowite for *x* > 0.5, and further supports the notion that cation site occupancy defects are reduced in substituted barlowite compared with herbertsmithite owing to the differences in their crystal strucures^{27}.

However, for members of the Zn-barlowite series with *x* > 0.5—for which we have demonstrated that magnetic order is suppressed—one would not expect the internal field distribution, Δ, within a particular sample to vary significantly with temperature. Therefore, by normalising the muon spin relaxation rates obtained for highly substituted members of the Zn-barlowite series, shown in Supplementary Fig. 6, we can directly compare their dynamical magnetic moment correlations across the series, which will be primarily governed by their temperature-dependent field fluctuation time, *τ*. Figure 5b shows the comparison of the normalised relaxation rate, *λ*_{norm}, for Zn_{x}Cu_{4−x}(OH)_{6}FBr with *x* = 0.66(1) and 0.99(1). Both samples have a prominent plateau in the muon spin relaxation that persists to the lowest measured temperatures—in the case of the *x* = 0.99(1) sample, to at least 50 mK—indicating an absence of any ordering or freezing of Cu^{2+} magnetic moments despite the strong exchange interactions between them (∣*θ*∣ ≈ 190 K for *x* = 0.99(1), (see Fig. 2b). Our observations are thus in support of the recent proposals of a QSL ground state for fully subsituted ZnCu_{3}(OH)_{6}FBr from both experiment and theory^{30,50}. Another important result from our present study, however, is the striking similarity of the temperature dependence of *λ*_{norm} for *x* = 0.66(1) and 0.99(1) samples, as seen in Fig. 5b. This suggests that the nature of the dynamical field fluctuations in highly substituted Zn-barlowite—and its potential QSL state—is remarkably robust to the presence of interlayer exchange between the kagomé layers of its hexagonal crystal structure due to partial substitution at the interlayer sites or that the combination of chemically randomised magnetic interactions and geometric frustration results in a theoretically proposed disorder-induced QSL^{51}. In either case, we conclude that *x* ≈ 0.5 marks a critical threshold in the Zn-barlowite series, above which dynamic magnetic moment correlations characterise the low-temperature magnetic behaviour. We note that this is in agreement with a recent study which proposes that magnetic moments in single crystals of Zn-barlowite with *x* = 0.56 remain dynamic to low-temperatures via bulk magnetometry measurements^{50}.

In summary, we have demonstrated that the Zn-barlowite series offers a promising route to a new QSL phase in an \(S=\frac{1}{2}\) kagomé antiferromagnet. By combining *μ*SR experiment with DFT, we find that the electronic magnetic moments in the parent material of this series, barlowite, Cu_{4}(OH)_{6}FBr, are static below *T*_{N} = 15 K and that *μ*SR spectra can be quantitatively interpreted by the formation of both *μ*–F and *μ*–OH complexes. We note that the application of DFT methods may be especially important in the *μ*SR study of QSL candidates, as they allow us to assess not only the location of implanted muons but also the degree of distortion that they cause, thus providing confidence that the muon spin probes the intrinsic properties of the quantum disordered ground state^{52}. In this regard, barlowite offers a valuable benchmark for DFT muon-site calculation and its application to other materials in which there is neither long-range magnetic order nor muon complex formation to give any indication as to where muons might implant. We find that incorporating as little as *x* = 0.16 of Zn^{2+} into the crystal structure of barlowite leads to a suppression of its ordered, frozen fraction, while for *x* > 0.66 the electronic magnetic moments remain dynamically fluctuating at all measurable temperatures. Whether the magnetic ground states of such partially substituted compositions of Zn-barlowite correspond to a QSL state—as proposed for the *x* = 1.0 end-member^{21,29,30,50}—awaits further investigation, for instance, by inelastic neutron scattering.

Also of note is that we find no evidence in our *μ*SR measurements for the coexistence of long-range magnetic order and dynamic magnetic moment fluctuations in barlowite, which has been previously observed by other groups via NMR measurements^{53} and in the structurally related clinoatacamite via *μ*SR measurements^{33}. Indeed, the reduced ordered moments obtained from the magnetic structure refinement of barlowite also indicate the presence of persistent dynamics in its ordered state^{23}. However, an additional consideration here is the potential sample dependence on the magnetic properties of barlowite. Whilst the subtle structural differences depending on synthetic route are unlikely to have any significant effect on the muon stopping sites, some studies report several magnetic transitions in their samples synthesised by alternative routes to us which may eventually be insightful to explore using *μ*SR^{20,22,24,50}.

To conclude this present study, we finish by noting the striking similarity of the temperature dependence of the normalised muon spin relaxation rates, *λ*_{norm}, not only between members of the Zn-barlowite series, but also with members of the related Zn-paracamite series, and in particular, the QSL candidate herbertsmithite^{16}, as shown in Fig. 5b. The apparent equivalence of the dynamics of the magnetic moment fluctuations that drive both of these related families of \(S=\frac{1}{2}\) kagomé antiferromagnets into the plateau state is perhaps indicative of a universality to the way in which such hydroxyl halides evolve towards a QSL regime upon chemical substitution. On the other hand, it could rather indicate the role of the interlayer Cu^{2+} moments present within the crystal structures of both systems, which has recently been proposed as an in situ probe of their QSL ground states via a Kondo-like interaction between the magnetic moments at the interlayer sites and the spinon excitations inherent to the possible QSL phases formed by the moments within their kagomé layers^{54}. Further direct comparison of the evolution of magnetic moment correlations in the Zn_{x}Cu_{4−x}(OH)_{6}FBr and Zn_{x}Cu_{4−x}(OH)_{6}Cl_{2} series across complementary time and length scales by, for example, muon Knight shift, NMR spectroscopy or neutron spin-echo techniques, could, therefore, provide much-needed further experimental insight into the roles of interlayer coupling, cation occupancy disorder and local structure in the ground state selection of the \(S=\frac{1}{2}\) kagomé antiferromagnet.

## Methods

### Hydrothermal synthesis

Polycrystalline samples of barlowite, Cu_{4}(OH)_{6}FBr, were synthesised via a hydrothermal reaction as previously reported^{23}. Zn-substituted samples, Zn_{x}Cu_{4−x}(OH)_{6}FBr, were achieved through the addition of excess ZnBr_{2} whilst varying the amount of CuBr_{2} accordingly as detailed in Supplementary Table 1. For our fully Zn-substituted sample, *x* = 0.99(1), the autoclave was heated at 1 K/min to 483 K, held for 24 h and cooled at a rate of 0.1 K/min to room temperature. The resulting products, coloured turquoise to pale blue depending on the level of Zn-substitution, were filtered and washed several times with distilled water and dried in air.

### XRD, inductively coupled plasma optical emission spectrometry (ICP-OES) and superconducting quantum interference device (SQUID) magnetometry

Sample purity was confirmed via powder X-ray diffraction (PXRD) using a Bruker d8 diffractometer in Debye-Scherrer geometry and utilising Cu K_{α1} (*λ* = 0.154056 Å) radiation. The Cu and Zn content was determined using ICP-OES on an Agilent 5110 SVDV ICP-OES. Samples were digested in dilute HNO_{3} and raw data are detailed in Supplementary Table 2. Time-of-flight PND measurements were performed on a sample containing *x* = 0.60(1) on the General Materials (GEM) diffractometer at the ISIS Facility of the Rutherford Appleton Laboratory^{55}. The 1.4 g sample was loaded into a 6 mm vanadium can and measured at 2 and 200 K. Magnetic susceptibility data were measured on a Quantum Design magnetic properties measurement system (MPMS) with a SQUID magnetometer. Samples were measured in VSM mode in an applied field of 1 T between 2 and 300 K in 2 K intervals.

### Muon spectroscopy measurements

*μ*SR data were collected on the MuSR spectrometer at the ISIS Neutron and Muon Source^{56,57}. Samples ranging from 60 to 500 mg and *x* = 0.00 to 0.99(1) were packed into silver foil sachets with sizes ranging from 2.25 to 4 cm^{2}, attached to a silver backing plate with vacuum grease and loaded into a ^{4}He cryostat. For *x* = 0.99(1) two additional silver sheets were placed in front of the sachet due to the small sample size of 60 mg. For *x* = 0.66(1) and 0.99(1), samples were additionally measured in a ^{3}He/^{4}He dilution fridge. Data were collected in zero-field (ZF), longitudinal-field (LF) and transverse-field (TF) geometries and analysed using the MANTID software^{58}. A separate sample of Cu_{4}(OH)_{6}FBr was measured on the General Purpose Surface-Muon (GPS) instrument of the Swiss Muon Source at Paul Scherrer Institut. The 90 mg sample was contained in an aluminium foil packet, suspended in the muon beam using a silver fork-type holder and measured in VETO mode to reduce background signal.

### DFT calculations

DFT muon-site calculations were carried out using the MuFinder software^{59} and the plane-wave-based code CASTEP^{60} using the local density approximation. A supercell consisting of 1 × 1 × 2 unit cells of the *P**n**m**a* structural model of barlowite was used in order to minimise the effects of muon self-interaction resulting from the periodic boundary conditions. Muons, modelled by an ultrasoft hydrogen pseudopotential, were initialised in low-symmetry positions and the structure was allowed to relax (keeping the unit cell fixed) until the change in energy per ion was <1 × 10^{−5} eV. A cutoff energy of 544 eV and a 1 × 2 × 1 Monkhorst-Pack grid^{61} were selected for *k*-point sampling.

## Data availability

Data supporting the findings within this study are available from the corresponding author upon request. Raw data sets from ISIS experiments can be accessed via links provided in refs. ^{55,56,57}.

## Code availability

The MuFinder program used to classify and analyse muon stopping sites can be downloaded, along with installation and usage instructions, via the link provided in ref. ^{59}.

## References

- 1.
Greedan, J. E. Geometrically frustrated magnetic materials.

*J. Mater. Chem.***11**, 37–53 (2001). - 2.
Broholm, C. et al. Quantum spin liquids.

*Science***367**, eaay0668 (2020). - 3.
Knolle, J. & Moessner, R. A field guide to spin liquids.

*Annu. Rev. Condens. Matter Phys.***10**, 451–472 (2019). - 4.
Anderson, P. W. Resonating valence bonds: a new kind of insulator?

*Mat. Res. Bull.***8**, 153–160 (1973). - 5.
Savary, L. & Balents, L. Quantum spin liquids.

*Rep. Prog. Phys.***80**, 016502 (2017). - 6.
Yan, S., Huse, D. A. & White, S. R. Spin-liquid ground state of the \(S=\frac{1}{2}\) kagome Heisenberg antiferromagnet.

*Science***332**, 1173–1176 (2011). - 7.
Han, T.-H. et al. Fractionalized excitations in the spin-liquid state of a kagome-lattice antiferromagnet.

*Nature***492**, 406–410 (2012). - 8.
Fu, M., Imai, T., Han, T.-H. & Lee, Y. S. Evidence for a gapped spin-liquid ground state in a kagome Heisenberg antiferromagnet.

*Science***350**, 655–658 (2015). - 9.
Khuntia, P. et al. Gapless ground state in the archetypal quantum kagome antiferromagnet ZnCu

_{3}(OH)_{6}Cl_{2}.*Nat. Phys.***16**, 469–474 (2020). - 10.
Fåk, B. et al. Kapellasite: a kagome quantum spin liquid with competing interactions.

*Phys. Rev. Lett.***109**, 037208 (2012). - 11.
Boldrin, D., Knight, K. & Wills, A. S. Orbital frustration in the \(S=\frac{1}{2}\) kagome magnet vesignieite kagome magnet vesignieite, BaCu

_{3}V_{2}O_{8}(OH)_{2}.*J. Mater. Chem.***4**, 10315–10322 (2016). - 12.
Hiroi, Z., Ishikawa, H., Yoshida, H., Yamaura, J.-i. & Okamoto, Y. Orbital transitions and frustrated magnetism in the kagome-type copper mineral volborthite.

*Inorg. Chem.***58**, 11949–11960 (2019). - 13.
Shores, M. P., Nytko, E. A., Bartlett, B. M. & Nocera, D. G. A structurally perfect \(S=\frac{1}{2}\) kagomé antiferromagnet.

*J. Am. Chem. Soc.***127**, 13462–13463 (2005). - 14.
Helton, J. S. et al. Spin dynamics of the spin-1/2 kagome lattice antiferromagnet ZnCu

_{3}(OH)_{6}Cl_{2}.*Phys. Rev. Lett.***98**, 107204 (2007). - 15.
de Vries, M. A., Kamenev, K. V., Kockelmann, W. A., Sanchez-Benitez, J. & Harrison, A. Magnetic ground state of an experimental \(S=\frac{1}{2}\) kagome antiferromagnet.

*Phys. Rev. Lett.***100**, 157205 (2008). - 16.
Mendels, P. et al. Quantum magnetism in the paratacamite family: towards an ideal kagomé lattice.

*Phys. Rev. Lett.***98**, 077204 (2007). - 17.
Freedman, D. E. et al. Site specific X-ray anomalous dispersion of the geometrically frustrated kagomé magnet, herbertsmithite, ZnCu

_{3}(OH)_{6}Cl_{2}.*J. Am. Chem. Soc.***132**, 16185–16190 (2010). - 18.
Nilsen, G. J., De Vries, M. A., Stewart, J. R., Harrison, A. & Rønnow, H. M. Low-energy spin dynamics of the \(S=\frac{1}{2}\) kagome system herbertsmithite.

*J. Phys.:Condens. Matter***25**, 106001 (2013). - 19.
Elliott, P., Cooper, M. A. & Pring, A. Barlowite, Cu

_{4}FBr(OH)_{6}, a new mineral isotructural with claringbullite: description and crystal structure.*Mineral. Mag.***78**, 1755–1762 (2014). - 20.
Han, T.-H., Singleton, J. & Schlueter, J. A. Barlowite: a spin-1/2 antiferromagnet with a geometrically perfect kagome motif.

*Phys. Rev. Lett.***113**, 227203 (2014). - 21.
Feng, Z. et al. Effect of Zn doping on the antiferromagnetism in kagome Cu

_{4−x}Zn_{x}(OH)_{6}FBr.*Phys. Rev. B***98**, 155127 (2018). - 22.
Pasco, C. M. et al. Single-crystal growth of Cu

_{4}(OH)_{6}BrF and universal behavior in quantum spin liquid candidates synthetic barlowite and herbertsmithite.*Phys. Rev. Mater.***2**, 044406 (2018). - 23.
Tustain, K., Nilsen, G. J., Ritter, C., da Silva, I. & Clark, L. Nuclear and magnetic structures of the frustrated quantum antiferromagnet barlowite, Cu

_{4}(OH)_{6}FBr.*Phys. Rev. Mater.***2**, 111405 (2018). - 24.
Smaha, R. W., He, W., Sheckelton, J. P., Wen, J. & Lee, Y. S. Synthesis-dependent properties of barlowite and Zn-substituted barlowite.

*J. Solid State Chem.***268**, 123–129 (2018). - 25.
Henderson, A. et al. Order-disorder transition in the \(S=\frac{1}{2}\) antiferromagnets barlowite and claringbullite.

*Chem. Commun.***55**, 11587–11590 (2019). - 26.
Han, T.-H., Isaacs, E. D., Schlueter, J. A. & Singleton, J. Anisotropy: spin order and magnetization of single-crystalline Cu

_{4}(OH)_{6}FBr barlowite.*Phys. Rev. B***93**, 214416 (2016). - 27.
Liu, Z., Zou, X., Mei, J. W. & Liu, F. Selectively doping barlowite for quantum spin liquid: a first-principles study.

*Phys. Rev. B***92**, 220102(R) (2015). - 28.
Guterding, D., Valentí, R. & Jeschke, H. O. Reduction of magnetic interlayer coupling in barlowite through isoelectronic substitution.

*Phys. Rev. B***94**, 125136 (2016). - 29.
Feng, Z. et al. Gapped spin-1/2 spinon excitations in a new kagome quantum spin liquid compound Cu

_{3}Zn(OH)_{6}FBr.*Chin. Phys. Lett.***34**, 077502 (2017). - 30.
Wei, Y. et al. Evidence for a Z

_{2}topological ordered quantum spin liquid in a kagome-lattice antiferromagnet. Preprint at http://arxiv.org/abs/1710.02991 (2020). - 31.
Depenbrock, S., McCulloch, I. P. & Schollwöck, U. Nature of the spin-liquid ground state of the \(S=\frac{1}{2}\) Heisenberg model on the kagome lattice.

*Phys. Rev. Lett.***109**, 067201 (2012). - 32.
Blundell, S. J. Spin-polarized muons in condensed matter physics.

*Contemp. Phys.***40**, 175–192 (1999). - 33.
Zheng, X. G. et al. Coexistence of long-range order and spin fluctuation in geometrically frustrated clinoatacamite Cu

_{2}ClOH_{3}.*Phys. Rev. Lett.***95**, 057201 (2005). - 34.
Kermarrec, E. et al. Spin-liquid ground state in the frustrated kagome antiferromagnet MgCu

_{3}(OH)_{6}Cl_{2}.*Phys. Rev. B***84**, 100401 (2011). - 35.
Barthélemy, Q. et al. Local study of the insulating quantum kagome antiferromagnets YCu

_{3}(OH)_{6}O_{x}Cl_{3−x}(*x*= 0, 1/3).*Phys. Rev. Mater.***3**, 074401 (2019). - 36.
Zorko, A. et al. Coexistence of magnetic order and persistent spin dynamics in a quantum kagome antiferromagnet with no intersite mixing.

*Phys. Rev. B***99**, 214441 (2019). - 37.
Clark, L. et al. Gapless spin liquid ground state in the \(S=\frac{1}{2}\) vanadium oxyfluoride kagome antiferromagnet [NH

_{4}]_{2}[C_{7}H_{1}4N][V_{7}O_{6}F_{1}8].*Phys. Rev. Lett.***110**, 207208 (2013). - 38.
Foronda, F. R. et al. Anisotropic local modification of crystal field levels in Pr-based pyrochlores: a muon-induced effect modeled using density functional theory.

*Phys. Rev. Lett.***114**, 017602 (2015). - 39.
Möller, J. S. Playing quantum hide-and-seek with the muon: localizing muon stopping sites.

*Phys. Scr.***88**, 068510 (2013). - 40.
Lord, J. S., Cottrell, S. P. & Williams, W. G. Muon spin relaxation in strongly coupled systems.

*Phys. B Condens. Matter***289**, 495–498 (2000). - 41.
Yaouanc, A. & Dalmas de Réotier, P.

*Muon Spin Rotation, Relaxation, and Resonance*(Oxford University Press, 2010). - 42.
Wei, Y. et al. Magnetic phase diagram of Cu

_{4−x}Zn_{x}(OH)_{6}FBr studied by neutron-diffraction and*μ*SR techniques. Preprint at http://arxiv.org/abs/2007.11816 (2020). - 43.
Slichter, C. P.

*Principles of Magnetic Resonance*(Springer Verlag, 1978). - 44.
Gardner, J. S. et al. Cooperative paramagnetism in the geometrically frustrated pyrochlore antiferromagnet Tb

_{2}Ti_{2}O_{7}.*Phys. Rev. Lett.***82**, 1012–1015 (1999). - 45.
Li, Y. et al. Muon spin relaxation evidence for the U(1) quantum spin-liquid ground state in the triangular antiferromagnet YbMgGaO

_{4}.*Phys. Rev. Lett.***117**, 097201 (2016). - 46.
Mustonen, O. et al. Spin-liquid-like state in a spin-1/2 square-lattice antiferromagnet perovskite induced by

*d*^{10}-*d*^{0}cation mixing.*Nat. Commun.***9**, 1085 (2018). - 47.
Kenney, E. M. et al. Coexistence of static and dynamic magnetism in the Kitaev spin liquid material Cu

_{2}IrO_{3}.*Phys. Rev. B***100**, 094418 (2019). - 48.
Dunsiger, S. R. et al. Muon spin relaxation investigation of the spin dynamics of geometrically frustrated antiferromagnets Y

_{2}Mo_{2}O_{7}and Tb_{2}Mo_{2}O_{7}.*Phys. Rev. B***54**, 9019–9022 (1996). - 49.
Gomilšek, M. et al.

*μ*SR insight into the impurity problem in quantum kagome antiferromagnets.*Phys. Rev. B***94**, 024438 (2016). - 50.
Smaha, R. W. et al. Materializing rival ground states in the barlowite family of kagome magnets: quantum spin liquid, spin ordered, and valence bond crystal states.

*npj Quant. Mater.***5**, 23 (2020). - 51.
Kawamura, H., Watanabe, K. & Shimokawa, T. Quantum spin-liquid behavior in the spin-\(\frac{1}{2}\) random-bond Heisenberg antiferromagnet on the kagome lattice.

*J. Phys. Soc. Jpn.***83**, 103704 (2014). - 52.
Lancaster, T. et al. Quantum magnetism in molecular spin ladders probed with muon-spin spectroscopy.

*N. J. Phys.***20**, 103002 (2018). - 53.
Ranjith, K. M. et al. Magnetic resonance as a local probe for kagome magnetism in barlowite Cu

_{4}(OH)_{6}FBr.*Sci. Rep.***8**, 10851 (2018). - 54.
Gomilšek, M. et al. Kondo screening in a charge-insulating spinon metal.

*Nat. Phys.***15**, 754–758 (2019). - 55.
Tustain, K. et al. STFC ISIS Facility, https://doi.org/10.5286/ISIS.E.95671117 (2018).

- 56.
Tustain, K. et al. STFC ISIS Facility, https://doi.org/10.5286/ISIS.E.RB1720027 (2018).

- 57.
Tustain, K. et al. STFC ISIS Facility, https://doi.org/10.5286/ISIS.E.RB1910387 (2019).

- 58.
Arnold, O. et al. Mantid-Data analysis and visualization package for neutron scattering and

*μ*SR experiments.*Nucl. Instrum. Methods Phys. Res.***764**, 156–166 (2014). - 59.
Huddart, B. M. https://gitlab.com/BenHuddart/mufinder/ (2020).

- 60.
Clark, S. J. et al. First principles methods using CASTEP.

*Z. Kristallogr.***220**, 567–570 (2005). - 61.
Monkhorst, H. J. & Pack, J. D. Special points for Brillouin-zone integrations.

*Phys. Rev. B***13**, 5188 (1976).

## Acknowledgements

L.C. and K.T. acknowledge the University of Liverpool for financial support and a studentship to K.T. B.M.H. thanks the STFC for the provision of a studentship. T.-H.H. acknowledges the support of the Grainger Fellowship provided by the University of Chicago. Work at the ISIS Neutron and Muon Source was supported by the STFC. Part of this work is supported by the EPSRC grant no EP/N024028/1. F.B. acknowledges support from the French Agence Nationale de la Recherche under Grant No. ANR-18-CE30-0022 LINK and from the Université Paris-Sud (MRM MAGMAG grant). The authors are grateful to S. Moss and the University of Liverpool Analytical Services for performing ICP-OES measurements, H. Niu and B. Le Pennec for assistance with susceptibility measurements, as well as I. da Silva for assistance on GEM. We acknowledge computing resources provided by Durham Hamilton HPC.

## Author information

### Affiliations

### Contributions

K.T., B.W.-O’ B. and T.-H.H. prepared the samples. K.T., F.B., P.J.B. and L.C. performed the MuSR experiments and K.T. analysed the data. F.B. and H.L. performed the GPS experiments and F.B. analysed the data. T.L. and B.M.H. performed the DFT calculations. K.T. performed supplementary experiments (magnetic susceptibility and PND) and analysed the data. L.C. conceived of and supervised the study. L.C. and K.T. wrote the manuscript with input from all of the authors.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

**Publisher’s note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Supplementary information

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Tustain, K., Ward-O’Brien, B., Bert, F. *et al.* From magnetic order to quantum disorder in the Zn-barlowite series of *S* = 1/2 kagomé antiferromagnets.
*npj Quantum Mater.* **5, **74 (2020). https://doi.org/10.1038/s41535-020-00276-4

Received:

Accepted:

Published: