Next Article in Journal
Dynamical Thermalization of Interacting Fermionic Atoms in a Sinai Oscillator Trap
Previous Article in Journal
Instabilities of the Vortex Lattice and the Peak Effect in Single Crystal YBa2Cu4O8
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Thermodynamic, Dynamic, and Transport Properties of Quantum Spin Liquid in Herbertsmithite from an Experimental and Theoretical Point of View

by
Vasily R. Shaginyan
1,2,*,
Alfred Z. Msezane
2,
Miron Ya. Amusia
3,4,
John W. Clark
5,6,
George S. Japaridze
2,
Vladimir A. Stephanovich
7 and
Yulya S. Leevik
8
1
Petersburg Nuclear Physics Institute of NRC “Kurchatov Institute”, 188300 Gatchina, Russia
2
Clark Atlanta University, Atlanta, GA 30314, USA
3
Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel
4
Ioffe Physical Technical Institute, RAS, 194021 St. Petersburg, Russia
5
McDonnell Center for the Space Sciences & Department of Physics, Washington University, St. Louis, MO 63130, USA
6
Centro de Investigação em Matemática e Aplicações, University of Madeira, 9020-105 Funchal, Madeira, Portugal
7
Institute of Physics, Opole University, Oleska 48, 45-052 Opole, Poland
8
Department of Finance, National Research University Higher School of Economics, 194100 St. Petersburg, Russia
*
Author to whom correspondence should be addressed.
Condens. Matter 2019, 4(3), 75; https://doi.org/10.3390/condmat4030075
Submission received: 15 June 2019 / Revised: 31 July 2019 / Accepted: 31 July 2019 / Published: 7 August 2019

Abstract

:
In our review, we focus on the quantum spin liquid (QSL), defining the thermodynamic, transport, and relaxation properties of geometrically frustrated magnet (insulators) represented by herbertsmithite ZnCu 3 ( OH ) 6 Cl 2 . The review mostly deals with an historical perspective of our theoretical contributions on this subject, based on the theory of fermion condensation closely related to the emergence (due to geometrical frustration) of dispersionless parts in the fermionic quasiparticle spectrum, so-called flat bands. QSL is a quantum state of matter having neither magnetic order nor gapped excitations even at zero temperature. QSL along with heavy fermion metals can form a new state of matter induced by the topological fermion condensation quantum phase transition. The observation of QSL in actual materials such as herbertsmithite is of fundamental significance both theoretically and technologically, as it could open a path to the creation of topologically protected states for quantum information processing and quantum computation. It is therefore of great importance to establish the presence of a gapless QSL state in one of the most prospective materials, herbertsmithite. In this respect, the interpretation of current theoretical and experimental studies of herbertsmithite are controversial in their implications. Based on published experimental data augmented by our theoretical analysis, we present evidence for the the existence of a QSL in the geometrically frustrated insulator herbertsmithite ZnCu 3 ( OH ) 6 Cl 2 , providing a strategy for unambiguous identification of such a state in other materials. To clarify the nature of QSL in herbertsmithite, we recommend measurements of heat transport, low-energy inelastic neutron scattering, and optical conductivity σ ¯ in ZnCu 3 ( OH ) 6 Cl 2 crystals subject to an external magnetic field at low temperatures. Our analysis of the behavior of σ ¯ in herbertsmithite justifies this set of measurements, which can provide a conclusive experimental demonstration of the nature of its spinon-composed quantum spin liquid. Theoretical study of the optical conductivity of herbertsmithite allows us to expose the physical mechanisms responsible for its temperature and magnetic field dependence. We also suggest that artificially or spontaneously introducing inhomogeneity at nanoscale into ZnCu 3 ( OH ) 6 Cl 2 can both stabilize its QSL and simplify its chemical preparation, and can provide for tests that elucidate the role of impurities. We make predictions of the results of specified measurements related to the dynamical, thermodynamic, and transport properties in the case of a gapless QSL.

1. Introduction

The frustrated magnet (insulator) herbertsmithite ZnCu 3 ( OH ) 6 Cl 2 is one of the best candidates for identification as a material that hosts a quantum spin liquid (QSL), thereby determining the nature of its thermodynamic, relaxation, and transport properties. The insulating nature of ZnCu 3 ( OH ) 6 Cl 2 has been established: There is a 3.3 eV charge gap (see e.g., [1,2]). In our review, we focus on a spin gap of QSL, and for brevity will call it “spin gap”. At low temperatures T, a QSL may have or may not have a gap in its excitation spectrum of spinons, which are fermion quasiparticles of zero charge that occupy the corresponding Fermi sphere with Fermi momentum p F . We note that, in contrast to metals, spinons cannot support charge current, but these, for instance, can carry heat, as electrons of metals do. The influence of a gap on the properties is huge, for at low T the properties of herbertsmithite are similar to those of common insulators, while the properties resemble those of metals, provided that the gap is absent. Thus, it is a challenge to experimentally establish the presence of the gap and its value, for current theoretical and experimental studies of herbertsmithite are controversial and cannot give a definite answer to the challenge. In the present review, based primarily on our own theoretical work, we offer a number of experiments that can discern unambiguously the presence of a gapless QSL in ZnCu 3 ( OH ) 6 Cl 2 and other geometrically frustrated magnets. Our aim here is to let the readers to make their own qualified judgments about the presented theoretical approach and its ability to deliver an adequate description of many seemingly different experiments in a uniform manner.
In a magnet with geometric frustration, where the spins cannot be ordered even at temperatures close to absolute zero, they stay in a liquid quantum spin state. Herbertsmithite is an antiferromagnet with a kagome lattice of spins S = 1 / 2 . Recent experimental studies have shown its unusual properties [3,4,5,6,7,8].
The balance of electrostatic forces for the Cu 2 + ions in the kagome structure is such that they occupy distorted octahedral sites. The magnetic planes formed by the Cu 2 + S = 1 / 2 ions are interspersed with nonmagnetic Zn 2 + layers. In samples, Cu 2 + defects occupy the nonmagnetic Zn 2 + sites between the kagome layers with x 15 % probability, thus introducing randomness and inhomogeneity into the lattice [9]. As we shall see, the starring role of impurities in the formation of QSL is not clearly understood and is now under thorough investigation [10,11,12,13,14,15]. However, we suggest that the influence of impurities induced by the homogeneity on the properties of ZnCu 3 ( OH ) 6 Cl 2 can be tested by varying x. We note that the impurities randomly located at the nanoscale level can support QSL as it is observed in measurements on Zn x Cu 4 x OH 6 Cl 2 with 1 x 0.8 , while at x > 0.8 , glassy dynamics emerge [5]. Moreover, the impurities randomly located at the nanoscale level can stabilize QSL as it is observed in measurements on the verdazyl-based complex Zn ( hfac ) 2 ( A x B 1 x ) [15] and make QSL stable in highly magnetic fields, as in the case of Mg -doped SrCu 2 ( BO 3 ) 2 [14].
The experiments done on ZnCu 3 ( OH ) 6 Cl 2 have not found any traces of magnetic order in it. Nor they have found the spin freezing down to temperatures of around 50 mK. In these respects, herbertsmithite is the best candidate among quantum magnets to contain QSL described above [3,4,5,6,7,8]. These assessments are supported by model calculations indicating that an antiferromagnet on a kagome lattice has a gapless spin liquid ground state [16,17,18,19,20,21,22,23]. At the same time, recent suggestions [9,24,25] that there can exist a small spin-gap in the kagome layers may stand in conflict to this emerging picture (see also refs. [26,27,28] for a recent review). The obtained results are a combination of experimental and theoretical ones. The latter have been obtained in the framework of the model, which takes the presence of Cu impurities into account. The experimental data has been obtained by inelastic neutron scattering on ZnCu 3 ( OH ) 6 Cl 2 crystals. It is assumed that the influence of the Cu impurity ensemble on the observed properties of herbertsmithite may be disentangled from that of the kagome lattice geometry [9,24,25]. It is further assumed that the impurity ensemble may be represented as a result of the dilution of some prototypic cubic lattice. The model supposes that the spin gap occurs in magnetic fields up to 9 T. We note that measurements of the local spin susceptibility of ZnCu 3 ( OH ) 6 Cl 2 show that the kagome spins exhibit a spin gap [24], while the measurements of the bulk spin susceptibility evidence the absence of the gap [3]. We shall see below that it is not possible to separate contributions coming from the local and bulk susceptibility, for both the kagome and the impurities spins form an integral system.
At the same time, without a magnetic field, the spin susceptibility χ shows behavior in many respects similar to the Curie law. The latter behavior demonstrates that the copper spin ensemble plays the role of weakly interacting impurities [9,24,25]. The same behavior was recently reported in a new kagome quantum spin liquid candidate, Cu 3 Zn ( OH ) 6 FBr [29]. As a result, we observe a challenging contradiction between two sets of experimental data when some of them state the absent of a gap while the others present evidence in favor of a gap. Subsequently, we shall demonstrate that the model based on the Cu spin ensemble is rather synthetic. In other words, this model does not have the possibility of distinguishing the Cu ensemble and kagome lattice contributions. This is because the impurities, being embedded in the kagome host lattice, form a single integral entity at nanoscale. The above model may contradict the accumulated knowledge about the physical properties of ZnCu 3 ( OH ) 6 Cl 2 . This knowledge is the result of vast theoretical and experimental efforts related to static and dynamic properties of herbertsmithite [3,4,5,6,7,8,16,17,18,19,20,21,22,23]. Thus, when analyzing the physics of quantum spin liquid in the herbertsmithite, it is of crucial importance to verify the existence of a spin gap by experimental and theoretical studies, for the gap strongly influences all its thermodynamic, transport, and relaxation properties. To analyze QSL behavior theoretically, we employ the strongly correlated quantum spin liquid (SCQSL) [16,17,18,20,23] model. A simple kagome lattice may host a dispersionless topologically protected branch of the quasiparticle spectrum with zero excitation energy, known as a flat band [16,18,23,30,31,32]. In that case, the topological fermion condensation quantum phase transition (FCQPT) can be considered as a quantum critical point (QCP) of the ZnCu 3 ( OH ) 6 Cl 2 spinon-composed QSL. Spinons have zero charge, occupying the Fermi sphere up to the Fermi momentum p F [16,17,32]. Taking into account that we are dealing with the real chemical compound ZnCu 3 ( OH ) 6 Cl 2 rather than with an ideal kagome lattice, we have to bear in mind that the actual magnetic interaction in the compound can shift QSL away from FCQPT, before or beyond QCP. Thus, the actual critical point location has to be established by experimental data analysis. The real part of the optical conductivity σ ¯ , measured at low frequency on geometrically frustrated magnetic insulators, can yield important experimental clues to the nature of spinon-based QSL [33,34], especially at a wide range of temperatures T and magnetic fields B. Note that the consistent interpretation of the above data is a difficult theoretical problem [33].
The main aim of the present review is to expose QSL as a new state of matter, also formed by heavy fermion (HF) metals, and to attract attention to experimental studies of ZnCu 3 ( OH ) 6 Cl 2 that have the potential to reveal both the underlying physics of QSL and the presence or absence of a gap in spinon excitations that form the thermodynamic, transport, and relaxation properties. To clarify the nature of QSL in herbertsmithite, we recommend measurements of heat transport, low-energy inelastic neutron scattering, and optical conductivity σ ¯ in ZnCu 3 ( OH ) 6 Cl 2 single crystals subject to an external magnetic field at low temperatures. We show that QSL is situated near FCQPT, which stems from the dependence of σ ¯ on the external magnetic field. We argue that the general physical picture inherent in QSL can be “purified” of the microscopic contributions coming both from phonons and the impurities polluting each specific sample being studied. We suggest that the influence of impurities on the properties of ZnCu 3 ( OH ) 6 Cl 2 , induced by the inhomogeneity, can be tested by varying x.
In the following sections on thermodynamic properties, relaxation, and optical conductivity, we shall demonstrate that the existence of a gap within the impurity model of herbertsmithite may contradict recent experimental data collected on this material; instead, the impurities and kagome planes can form a genuine SCQSL. Accordingly, the observed gap [24] related to the kagome planes may not be a real one as it is not a physical mechanism for the observed thermodynamic, relaxation, and conductivity properties of ZnCu 3 ( OH ) 6 Cl 2 .

2. Frustrated Insulator Herbertsmithite

If the quasiparticles forming QSL are approximately dispersionless spinons (uncharged fermions), this state of matter becomes an SCQSL, [17,18,20,35,36]. If an insulating chemical compound has a 2D lattice, an SCQSL can emerge as a result of the lattice geometrical frustration. Most frequently, such frustration results in a dispersionless topologically protected spectral branch with zero excitation energy, called a flat band. This scenario realizes in herbertsmithite with its kagome planes [18,30,32,37,38,39], shown in Figure 1a. In this case, SCQSL has its quantum critical point at the FCQPT. As we have mentioned above, the SCQSL consists of chargeless S = 1 / 2 spinons having the effective mass M * . Their momenta, as usual for fermions, reside in a Fermi sphere with momentum p F . Note that one frustrated valence spin is taken away to populate the approximately flat spinon band, as is displayed in Figure 1b. Such a behavior strongly resembles that of HF metals whose valence electrons form conduction flat bands, while the charges e cannot form any band because of the large charge gap of herbertsmithite or another geometrically frustrated insulating magnet. This result is in agreement with the experimental data [40]. Therefore, the above insulating compounds have physical properties resembling those of HF metals. There is, however, a notable exception. Namely, while HF metals are good conductors, the typical insulator does not sustain the electric current [18,19,20]. Thus, both frustrated insulators and HF metals, having universal behavior, form a new state of matter [35].

3. Thermodynamic Properties

If QSL-forming quasiparticles are approximately dispersionless spinons (uncharged fermions), this state of matter becomes an SCQSL [17,18,20,35,36]. As we have mentioned above, the geometrical frustration in herbertsmithite promotes the SCQSL formation [18,19,20], while the presence of randomly scattered impurities can facilitate the frustration (see e.g., [15]). Note that in real chemical compounds with many lattice imperfections, the actual SCQSL emergence point is shifted from the theoretically predicted FCQPT. This means that the mutual location of the SCQSL and FCQPT can only be extracted from the experimental data.
Famous Fermi liquid theory, proposed by Landau (the so-called Landau Fermi liquid (LFL) theory) [41], has de facto been the universal tool to describe the itinerant fermionic systems. It maps the ensemble of strongly interacting electrons in a solid to the effective quasiparticle gas with not so strong an interaction. In such an approach, the excitations are represented in terms of the above quasiparticles so that the low-temperature properties of the system under consideration depend on the latter. The quasiparticles possess certain effective mass M * , which is a parameter of the theory [41,42,43], being approximately independent of external stimuli including temperature, pressure, or an electromagnetic field. The LFL theory cannot, however, explain the experimental results related to strong temperature and/or magnetic field dependence of the effective mass M * , as observed in strongly correlated Fermi systems [16,44]. At the same time, deviations from LFL behavior are observed in the vicinity of an FCQPT [16,17,44]. The above peculiarities are usually addressed as non-Fermi-liquid (NFL) properties. They are due to large effective mass, which actually becomes infinite at the FCQPT point (see [16,17] for details).
Let us describe the physical mechanism yielding the temperature and field dependences of the effective mass of the Landau quasiparticle M * ( B , T ) . Again, the key point is that upon approach to the FCQPT from the LFL regime, the effective mass becomes strongly external-stimuli-dependent. The aforementioned stimuli are the temperature, the magnetic field, and the external pressure P, to name a few [16,17]. This is indeed a consequence of an additional instability channel of a normal Fermi liquid (to those of Pomeranchuk, [42]). Note that the new channel is activated when the effective mass approaches infinity.
To avoid unphysical situations related, for instance, to the negativity of the effective mass, the system alters the topology of its Fermi surface [32,38,39] so that the effective mass starts to depend on the above external stimuli [16,17]. For our studies of the effective mass properties in SCQSL, we adopt the so-called homogeneous Fermi liquid model. Such a description avoids consideration of nonuniversal (and highly irrelevant in our analysis) features like the exact structure of a specific sample [16,17]. In such a model, the LFL equation for M * ( B , T ) reads [16,17,41]
1 M * ( B , T ) = 1 M + σ 1 p F · p p F 3 F ( p F , p ) × n σ 1 ( p , T , B ) p d p ( 2 π ) 3 .
In this expression, M is the free electron mass, F ( p F , p ) is the Landau interaction function, depending on p (momentum), p F (constant Fermi momentum), and n σ ( p , T , B ) , which is the quasiparticle distribution function for the spin projection σ . The quasiparticle interaction F ( p , p 1 ) , assumed here to be spin-independent, is phenomenological. Without loss of generality, here we assume that F ( p , p 1 ) is independent of temperature so that n σ ( p , T ) has the Fermi–Dirac form
n σ ( p , T ) = 1 + exp ( ε σ ( p , T ) μ σ ) T 1 ,
where ε σ ( p , T ) is a single-particle energy spectrum. In addition, μ is a chemical potential, which is spin-dependent via Zeeman splitting μ σ = μ ± μ B B , where μ B is the Bohr magneton. As usual, the spectrum ε σ ( p , T ) is obtained variationally from the system energy E [ n σ ( p , T ) ] ,
ε σ ( p , T ) = δ E [ n ( p ) ] δ n σ .
In describing herbertsmithite as a strongly correlated quantum spin liquid, the choice of the function F ( p , p 1 ) is dictated by its possession of FCQPT [16]. Thus, the sole role of the Landau interaction is to drive the system to the FCQPT point, at which the topology of the Fermi surface is altered in such a way that the effective mass ceases to be a constant parameter obtaining the aforementioned dependence on external stimuli [16,17]. Performing the variation (3), we arrive at the expression for a single-particle energy spectrum
ε σ ( p , T ) p = p M F ( p , p 1 ) p n σ ( p 1 , T ) d 3 p 1 ( 2 π ) 3 .
Equations (2) and (4) constitute the closed set to find ε σ ( p , T ) and n σ ( p , T ) . In this case, the effective mass enters through the expression p F / M * = ε ( p ) / p | p = p F . At the FCQPT point, the analytical solution of Equation (1) is possible [16,17]. Namely, in contrast to the LFL approach with M * being a constant parameter, here at zero magnetic field, M * becomes temperature dependent. The latter feature comprises the strong deviation from the LFL picture, determining the NFL regime
M * ( T ) a T T 2 / 3 .
At elevated temperatures, the system undergoes a transition to the LFL region of the phase diagram, and being subjected to the magnetic field, exhibits the behavior
M * ( B ) a B B 2 / 3
of the effective mass.
The introduction of “internal” (or natural) scales greatly simplifies the analysis of the problem under consideration. We first observe that near the FCQPT, the effective mass M * ( B , T ) (i.e., the solution to Equation (1)) has a maximum M M * at a temperature T M B [16,17]. This means that it is convenient to measure the effective mass and temperature in the units M M * and T M , respectively. Thus, we arrive at normalized effective mass M N * = M * / M M * and temperature T N = T / T M . Near FCQPT, the dependence M N * ( T N ) can be rendered as a universal interpolating function [16,17]. This function describes the transition from LFL to NFL states, given by Equations (5) and (6), and presents the universal scaling behavior of M N * [16,17]
M N * ( y ) c 0 1 + c 1 y 2 1 + c 2 y 8 / 3 .
Here, y = T N = T / T M and c 0 = ( 1 + c 2 ) / ( 1 + c 1 ) , where c 1 and c 2 are free parameters. The magnetic field B enters Equation (1) only in the combination μ B B / T , making T M μ B B . It follows from Equation (7) that
T M a 1 μ B B ,
where a 1 is a dimensionless quantity. In this case, the variable y becomes y = T / T M T / μ B B . The Equation (8) permits us to assert Equation (7), giving the scaling properties of temperature and magnetic field dependence of the effective mass. This implies, in turn, that at different fields B, the curves M N * form a single one as a function of the normalized variable y = T / T M . As T and B enter symmetrically in Equation (7), it also manifests the scaling behavior of M N * ( B , T ) as a function of T / B :
T N = T T M = T a 1 μ B B T B B T .
The schematic phase diagram is portrayed in Figure 2. We assume for simplicity that at T = 0 and B = 0 , the system is near FCQPT. The external magnetic field and temperature are indeed the parameters controlling the system position on the phase diagram relatively to the FCQPT point. The same parameters are responsible for the transitions between the NFL and LFL regions of the phase diagram (see the arrows in Figure 2). The horizontal arrow corresponds to fixed temperatures so that the system motion along this arrow from the NFL to LFL region corresponds to a magnetic field increase. On the contrary, the vertical arrow fixes the magnetic field so that the motion along it signifies the elevated temperatures. The shaded area reports the region where the NFL state transits to a weakly polarized LFL one. The temperature T * ( B ) of the transition is defined by the expression
T * ( B ) T M ( B ) ,
which directly follows from Equation (8). The line Equation (10) is indeed the function T * μ B B , so that the width W ( B ) of the NFL region is proportional to the temperature. It can be shown similarly that the transition region width T W ( B ) is also T [16,17]. We note here that in essence, the transition region represents the crossover between LFL and NFL phases. In our case, the NFL phase is formed by quasiparticles that occupy the so-called fermion condensate (FC) state, in analogy to the Bose condensate for particles obeying Bose–Einstein statistics. (See [16,17] for a comprehensive explanation.) In a “pure” FC state, all fermions (quasiparticles) having momenta in a finite interval embracing the Fermi surface have energies pinned to the chemical potential. In reality, this state cannot be reached because of the Nernst theorem [45], and the NFL features arise from “traces” of the FC state manifested at finite temperatures. Moreover, at low but finite temperatures, the magnetic field acts to suppress NFL behavior (i.e., the “FC traces”) and, upon growing sufficiently strong, restores the LFL phase. On the other hand, thermal fluctuations destroy LFL behavior and generate NFL features related to the FC state.
To examine the impurity model and possible gap in the spinon single particle spectrum, we consider the measured properties of herbertsmithite magnetic susceptibility χ . It is seen from Figure 3 that the magnetic susceptibility diverges as χ ( T ) T 2 / 3 at B 1 T (full line). For weakly interacting impurities, it has been suggested that at low temperatures, the dependence χ ( T ) can be approximated by a Curie–Weiss law [9,24,25], i.e., χ CW ( T ) 1 / ( T + θ ) , where θ is a Curie temperature, which turns out to be very small. At the same time, with respect to firmly established behavior χ ( T ) T 2 / 3 at low B, the above Curie–Weiss approximation is in discord with both theory [17,18,20] and experiment [3]. Moreover, as seen in Figure 4, the normalized spin susceptibility χ N behaves like the normalized specific heat C e l / T extracted from measurements on YbRh 2 Si 2 in high magnetic fields [46] and displayed in Figure 5a. Note that C e l / T displayed in Figure 5a coincides approximately with C m a g ( B , T ) / T C e l ( B , T ) / T shown in Figure 6, pointing to similarity of the electron and spinon Fermi spheres. Figure 5b reports the normalized specific heat ( C e l / T ) N extracted from the data [46,47]. It is seen that at low magnetic fields, the electronic system is not polarized as it is at high ones (black line) [17,48]. We note that the same behavior is exhibited by both χ N shown in Figure 6 and ( C m a g / T ) N shown in Figure 7. These facts support our conclusion that the SCQSL of ZnCu 3 ( OH ) 6 Cl 2 behaves like both the HF electron liquid of YbRh 2 Si 2 and the insulator Cu ( C 4 H 4 N 2 ) ( NO 3 ) 2 [49,50].
In fact, at T = 0 , the FC state is represented by the superconducting state with the superconducting gap Δ = 0 , while the superconducting order parameter κ = n ( p ) ( 1 n ( p ) ) is finite in the region ( p i p f ) [16,31,51], for in the region n ( p ) < 1 .
The above observations point to a possible absence of the spin gap in ZnCu 3 ( OH ) 6 Cl 2 , for it is impossible to definitely separate the contributions of the impurities and of kagome planes. The impurity model implies that the intrinsic spin susceptibility of the kagome plane is decomposed as χ kag ( T ) = χ ( T ) χ CW ( T ) , leading to χ kag ( T 0 ) 0 and the result about the existence of a putative gap [25]. This shows that there is a problem with the impurity model as it cannot describe the firmly established behavior χ ( T ) T 2 / 3 [3]. Thus, to explain the observed behavior of χ , one should view the impurity ensemble embedded in the kagome host crystal lattice as an integral system [16,17,18,19,20,21,23,49,50,52] that acts coherently at nanoscale to produce additional frustration of the kagome planes, so as to make QSL robust at lower temperatures. Based on this analysis, we predict that QSL of quantum magnets can be stabilized by introducing a random distribution of impurities. We also note that artificially or spontaneously introducing impurities into ZnCu 3 ( OH ) 6 Cl 2 can both stabilize its QSL and simplify its chemical preparation.
It is instructive to plot the functions T M ( B ) and χ M ( B ) M M * ( B ) of ZnCu 3 ( OH ) 6 Cl 2 . These functions are displayed in Figure 8 with the data extracted from experimental facts [4]. It is seen from Figure 8a that the behavior of T M ( B ) B is in accordance with Equation (6). It is also seen from Figure 8a that ZnCu 3 ( OH ) 6 Cl 2 is not located exactly at the FCQPT, for T M ( B = 0 ) 0.4 K and the system demonstrate LFL behavior at T 0.4 K [4]. Similarly, it is seen from Figure 8b that χ m a x ( B ) M M * ( B ) , in accordance with Equations (6) and (7). Thus, the effect of the impurities on QSL can be analyzed basing on Figure 8: the LFL temperature T L F L grows if QSL is shifted by impurities from the FCQPT, while the gap is absent. It also follows from Figure 8 that ZnCu 3 ( OH ) 6 Cl 2 exhibits behavior typical of HF metals under the application of a magnetic field [16,17]. As a result, we conclude that impurities form the single integral entity with a kagome lattice at nanoscale.
The impurity model has been utilized by the authors of reference [9] to derive an intrinsic scattering measure S kag ( ω ) = S tot ( ω ) a S imp ( ω ) . In the latter expression, S tot ( ω ) is a total (i.e., host lattice plus impurities) scattering rate, while S imp ( ω ) is the impurity contribution with fitting parameter a. On finding that S kag ( ω ) goes to zero as ω < 0.7 meV (see Figure 4b of reference [9]), they assert the existence of a gap. However, as we have demonstrated above, such a procedure gives just a spurious gap. Indeed, the latter conclusion relies completely on the assumption of a weak inter-impurity interaction, which is unjustified empirically since the subtraction hypothesis is negated by the experimental behavior summarized in Figure 3.
We now examine the impurity model in further detail and put its conclusion about the existence of a spin gap under scrutiny. It is seen from Figure 3 that normal low-temperature Fermi liquid properties of the magnetic susceptibility χ are confirmed experimentally at least for B 3 T. In such strong magnetic fields and low temperatures, the impurity spins are aligned completely along the field direction. They do not have degrees of freedom, permitting them to demonstrate Curie–Weiss behavior. Thus, assuming that impurity spins are fully aligned along the magnetic field direction and hence do not contribute to χ , one remains with χ kag ( T ) = χ ( T ) . Similar properties of the heat capacity follow from Figure 6. In the above magnetic fields, C mag / T also demonstrates ordinary Fermi liquid behavior. This shows that at low temperatures and B 3 T, the impurities show minute contributions to χ and C mag / T . It becomes obvious that the main contribution in such a case comes from the host kagome lattice [9,24,25]. Moreover, the impurity model implies that both χ ( T ) and C mag ( T ) / T approach zero at T 0 and B 3 T. It is clear from Figure 3, Figure 4, Figure 6, and Figure 7 that this is not true. Namely, both χ and C mag / T do not go to zero at T 0 up to B 14 T. Moreover, the scaling properties of C mag / T from Figure 7 and the behavior of both T M ( B ) and χ m a x ( B ) displayed in Figure 8 confirm clearly the lack of a gap. The results of C mag measurements [4,5,6,7,8] are the same for powder and bulk samples (see Figure 6).
The aforementioned experimental results are consistent with the hypothesis that the vast majority of physical properties of herbertsmithite are due to stable SCQSL. First, there is no substantial gap in the spinon excitation spectrum. Moreover, the application of high magnetic fields 18 T does not trigger such a gap. This implies that the impurity model is inadequate. These findings are in conformity with recent measurements stating that the low-temperature plateau in the local susceptibility identifies the spin-liquid ground state as a gapless one [53], while recent theoretical analysis confirms the absence of a gap [22].

4. Dynamic, Relaxation, and Heat Transport Properties

Neutron scattering measurements are one more confirmation of our hypothesis’ validity. Recently, the dynamic spin susceptibility χ ( q , ω , T ) = χ ( q , ω , T ) + i χ ( q , ω , T ) has been measured by neutrons scattering as a function of q (momentum), ω (frequency), and T. At low temperatures, the results are consistent with the idea that these quasiparticles are spinons, which form an approximately flat band [54].
The imaginary part χ ( T , ω 1 ) satisfies the equation [17,19]
T 2 / 3 χ ( T , ω 1 ) a 1 ω 1 1 + a 2 ω 1 2 ,
where a 1 and a 2 are constants and ω 1 = ω / ( T ) 2 / 3 . Equation (11) demonstrates that T 2 / 3 χ ( T , ω 1 ) has a maximum ( T 2 / 3 χ ( T , ω 1 ) ) max at some ω max . Equation (11) describes the scaling properties of χ T 0.66 established in measurements in reference [3]. As in Equation (7), we introduce the dimensionless function ( T 2 / 3 χ ) N = T 2 / 3 χ / ( T 2 / 3 χ ) max and the (dimensionless) variable E N = ω 1 / ω max . In this case, Equation (11) is modified to read
( T 2 / 3 χ ) N b 1 E N 1 + b 2 E N 2 ,
where b 1 and b 2 are fitting parameters. They should be chosen from the condition that the right-hand side of Equation (12) should reach the maximum at E N = 1 . This means that the expression ( T 2 / 3 χ ) N = T 2 / 3 χ / ( T 2 / 3 χ ) max exhibits scaling as a function of E N [17,19]. In other words, we find that B 2 / 3 χ ( ω ) (cf. Equation (6)) exhibits the scaling behavior with E N = ω 1 / ω max :
( B 2 / 3 χ ) N d 1 E N 1 + d 2 E N 2 .
Similarly, d 1 and d 2 are fitting parameters chosen from the condition that the function ( B 2 / 3 χ ) N = 1 at E N = 1 . In the FCQPT point, the discussed scaling is valid almost to T = 0 .
Figure 9 reports the function ( T 2 / 3 χ ) N extracted from neutron-scattering measurements on the HF metal Ce 0.925 La 0.075 Ru 2 Si 2 [55]. The corresponding data collected on two systems, ZnCu 3 ( OH ) 6 Cl 2 [3] and ( D 3 O ) Fe 3 ( SO 4 ) 2 ( OD ) 6 [56], are reported in Figure 10 and Figure 11. The figures show pretty good agreement between theoretical [19] (solid curves) and experimental results on all considered chemical compounds over almost three decades in the scaled variable E N . Hence ( T 2 / 3 χ ) N also has scaling properties. This shows that the spinons in both ZnCu 3 ( OH ) 6 Cl 2 and ( D 3 O ) Fe 3 ( SO 4 ) 2 ( OD ) 6 demonstrate the same itinerant behavior as the conduction electrons in the heavy fermion compound Ce 0.925 La 0.075 Ru 2 Si 2 . The detection of the above itinerant behavior is very important as it convincingly shows the presence of a gapless SCQSL in herbertsmithite [17,19,20].
Under the assumption that a fermion condensate (FC) is present in the above HF metals, the imaginary part χ ( T , ω ) of the susceptibility reads [16,17]
T χ ( T , ω ) a 5 E 1 + a 6 E 2 ,
where E = ω / k B T , while a 5 and a 6 are adjustable parameters. Expression (14) shows that T χ ( T , ω ) is a function of a single variable E = ω / k B T . In this case, the expressions (12) and (14) define two types of scaling in χ ( ω , T ) . The dynamic susceptibility ( T χ ) taken from inelastic neutron scattering measurements on the HF metal YbRh 2 Si 2 [57] is reported in Figure 12. The scaling in ( T χ ) is clearly seen in both this function and variable E. This confirms both the feasibility of Equation (14) and the similarity of the HF metals and frustrated magnets properties. The scaled data acquired in measurements on such structurally different strongly correlated systems as ZnCu 3 ( OH ) 6 Cl 2 , Ce 0.925 La 0.075 Ru 2 Si 2 , ( D 3 O ) Fe 3 ( SO 4 ) 2 ( OD ) 6 , and YbRh 2 Si 2 merge in a single curve over almost three decades in the scaled variables, thus confirming that these strongly correlated Fermi systems exhibit universal scaling behavior and are symptomatic of the existence of a new state of matter [17,58]. This observation strengthens the credibility of the FC approach as it can reliably explain the experimental data concerned and has demonstrable predictive power [17,35,58].
Indeed, it is apparent from Figure 9, Figure 10, Figure 11 and Figure 12 that the calculations within the FC approach are in conformity with th experimental results, providing strong evidence that SCQSL is the underlying mechanism defining the properties of ZnCu 3 ( OH ) 6 Cl 2 and ( D 3 O ) Fe 3 ( SO 4 ) 2 ( OD ) 6 . We conclude that the spin gap apprehension is a rather artificial construction, which may contradict the acquired experimental knowledge about the ZnCu 3 ( OH ) 6 Cl 2 properties. The consistency of this description based on a FCQPT as the driving mechanism for the properties of herbertsmithite at low temperature can also be taken as strong evidence against the spin gap existence in this system. Note that the presence or absence of a spin gap is not that important as it does not define the physical properties of ZnCu 3 ( OH ) 6 Cl 2 . In particular, neither χ nor C mag / T vanishes in the high magnetic fields, eliminating the contribution coming from the impurities (see Figure 3, Figure 4, Figure 6, and Figure 7). Moreover, Figure 10 shows that such a gap does not contribute to the imaginary part of the magnetic susceptibility, which does not vanish at the lowest temperatures. Rather, all above physical properties are due to the SCQSL. This fact can be tested by magnetic field heat transport measurements similar to organic insulators EtMe 3 Sb [ Pd ( dmit ) 2 ] 2 and κ ( BEDT TTF ) 2 Cu 2 ( CN ) 3 [21,59,60]. Heat transport measurements are notably important as they probe the SCQSL excitations in ZnCu 3 ( OH ) 6 Cl 2 . They can also detect the itinerant spinons, which are the primary reason for heat transport. It is obvious that the above heat transport cannot do without phonons. On the other hand, the phonon contribution is barely affected by the magnetic field. To summarize, we anticipate that magnetic field dependence of heat transport measurements may be an important step forward to identify the SCQSL nature in ZnCu 3 ( OH ) 6 Cl 2 [17,20,21].
Note that if we set the electronic charge to zero, the SCQSL in herbertsmithite becomes identically similar to the itinerant electrons ensemble in HF metals. In this case, the SCQSL thermal resistivity w reads [17,20,21]
w w 0 = W r T 2 ρ ρ 0 ( M * ) 2 T 2 ,
where W r T 2 comes from spinon–spinon scattering, which is similar to the contribution A T 2 from electron–electron scattering to charge transport. Here, M * is the effective mass and ρ is the longitudinal magnetoresistivity (LMR), while w 0 and ρ 0 denote the residual thermal and electric resistivity, respectively.
Finally, we consider the magnetic field influence on the spin-lattice relaxation rate 1 / ( T 1 T ) . Figure 13A shows the normalized quantity 1 / ( T 1 T ) N as a function of the magnetic field. It shows that with a magnetic field increase, 1 / ( T 1 T ) decays progressively. Contrary to the previous cases, where the normalization has been fulfilled in the point of maximum, here we normalize our curves in the inflection point, B = B inf , marked by the arrows in the main panel and inset. The same procedure has been performed with the normalized magnetoresistance shown in the inset to Figure 13. The relation 1 / ( T 1 T ) N ( M * ) 2 shows that our system located near its QCP would exhibit similar behavior of 1 / ( T 1 T ) N [16,17,20,21]. Significantly, Figure 13A shows that the normalized spin-lattice relaxation rate of herbertsmithite [61] and the HF metal YbCu 5 x Au x [62] do exhibit the same behavior. As seen from Figure 13A for B B inf (or B N 1 ), the quantity 1 / ( T 1 T ) N is almost magnetic-field-independent. At the same time, at elevated magnetic fields it decays [16,17,20,21] according to
W r 1 / ( T 1 T ) N = ρ N ( M * ) 2 B 4 / 3 .
Thus, we hypothesize that magnetic field B yields a crossover between NFL and LFL regimes. It also reduces substantially the relaxation rate as well as the thermal resistivity, similar to the case of the normalized LMR of YbRh 2 Si 2 , (see Figure 13B). Experimental data for the LMR is taken from [63]. Also of high priority are measurements of low-energy inelastic neutron scattering on ZnCu 3 ( OH ) 6 Cl 2 crystalline samples under the application of a magnetic field, driving the system into the LFL sector of the phase diagram (see Figure 2). The latter measurements permit direct observation of a possible gap, since in this case the impurity contribution is insignificant.
Recent neutron-scattering measurements on YbMgGaO 4 reveal broad spin excitations covering a wide region of the Brillouin zone, thereby favoring the existence of a spinon Fermi surface [64,65]. At the same time, the measurements of heat transport do not show any substantial magnetic excitation contributions to thermal conductivity, thus raising doubt as to the presence of a QSL [66]. We speculate that the observed behavior can be attributed to the emergence of a Mott insulator (see e.g., [67]), placing the system beyond the FCQPT point in the phase diagram. At the same time, in the case of herbertsmithite, the transport and thermodynamic properties suggest that the gapless state represented by the SCQSL is situated before the FCQFT. We note that recent theoretical studies of possible gaps in the ground state of herbertsmithite lead to quite different conclusions about the presence of a gap [18,22,68,69,70]. Recent experimental studies also indicate that the spin-liquid ground state in kagome lattice is gapless [53], and that two distinct types of defects in herbertsmithite are found [71]. This observation makes the impurity model vulnerable. Therefore, to probe QSL in herbertsmithite reliably, it is of crucial importance to carry out the experimental studies suggested above.

5. Charge Transport and Optical Conductivity

Our next step is the analysis of the herbertsmithite optical conductivity σ ¯ at low frequencies. To avoid the contribution of phonon absorption to the conductivity, we consider low temperatures T and frequencies ω [33]. Under such assumptions, we can neglect the lattice symmetry (kagome vs. triangular) as the wavelength is much larger than the typical crystal size. In atomic units = c = 1 , the Hamiltonian of a particle having momentum p = i , spin s , and charge e is given by
H ^ = 1 2 m p e A 2 + e ϕ μ B s s · B ,
where A and ϕ are the vector and scalar potentials, respectively. The nature of the vector potential A = 0 implies that p and A operators commute. Since for spinons, e = 0 , only the last term on the right-hand side of Equation (17) contributes to σ ¯ .
Equations (5) and (11) show that at small ω , the imaginary part of the spin susceptibility reads [19]
χ ( ω ) ω ( M * ) 2 .
Observing that the energy transfer ε B [36,43] between magnetic field B ( ω ) and our system is due to the last term in Equation (17), we obtain
ε B = 2 π ω μ B s s · B χ ( ω ) 2 π ω 2 μ B 3 s ( M * ) 2 s · B .
The same quantity, ε E , coming from the electric field E ( ω ) , reads
ε E = E 2 ( ω ) σ ¯ ( ω ) .
Comparison of Equations (19) and (20) yields
σ ¯ ( ω ) ω χ ( ω ) ω 2 ( M * ) 2 .
It follows from Equation (21) that σ ¯ ( ω ) ω 2 , and that behavior is consistent with experimental facts obtained in measurements on ZnCu 3 ( OH ) 6 Cl 2 and EtMe 3 Sb [ Pd ( dmit ) 2 ] 2 , representing the best candidates for identification as materials that host QSL [2,33,72]. We note that χ (18) at low ω coincides Equation (11) and provides a good description of the experimental data (see Figure 9 and Figure 10). Equations (5) and (21) inform us that at elevated temperatures, σ ¯ ( T ) is a decreasing function, which conforms with the experimental data [33]. It also follows from Equations (6) and (21) that σ ¯ ( ω , B ) is a decreasing function of B. This observation seems to be contradictory as no systematic B dependence is observed experimentally [33].
To elucidate the σ ¯ ( B ) behavior, we note that corresponding experiments have been performed at T = 3 K and B 7 T [33]. This means that the system is still in the transition region of the phase diagram and does not yet have the normal Fermi liquid properties, where the effective mass M * is given by Equation (6). In this case, the effective mass properties is due to Equation (5), rather than Equation (6). This means that the dependence σ ¯ ( B ) cannot be observed. Accordingly, we think that substantial σ ¯ ( B ) dependence can be observed at B 7 T in the case of T 1 K. In that case, as one sees from Figure 6, the effective mass M * C mag / T is a diminishing function of the applied magnetic field. Thus, we predict that σ ¯ ( B ) diminishes with growing magnetic fields, as follows from Equations (6) and (21). Since the contribution coming from phonons does not depend on the magnetic field, we propose that measurements of the variation δ σ , i.e., δ σ = σ ¯ ( ω , B ) σ ¯ ( ω , B = 0 ) , can reveal both the physics of SCQSL and the ground state of ZnCu 3 ( OH ) 6 Cl 2 , as well as the ground state of other materials hosting a QSL. The above experiments on measurements of the heat transport and optical conductivity can be carried out on samples with varying x. Such experiments yield information on the influence of impurities on the value of the gap. We predict that at moderate x 20 %, the SCQSL remains robust, for both inhomogeneity and randomness facilitate frustration.

6. Summary

The central message of this review is that to achieve a satisfactory understanding of the quantum spin-liquid physics of herbertsmithite at the nanoscale, it is essential to perform the targeted measurements on ZnCu 3 ( OH ) 6 Cl 2 that we have discussed. These focus on heat transport, low-energy inelastic neutron scattering, and optical conductivity σ ¯ , in the presence of magnetic fields at low temperatures. Moreover, we have suggested that the increasing x, i.e., the percentage of Zn sites that are occupied by Cu , can facilitate the frustration of the lattice and thereby act to stabilize the quantum spin liquid state. This conjecture can be tested in experiments on samples of herbertsmithite with different x values under the application of a magnetic field. To be specific, the results of the measurements might yield an unambiguous answer to the question of whether the gap in spinon excitations really exists, and if so, how it depends on impurities. In addition, such experiments will help to separate the universal effects from those coming both from phonons and from the impurities that pollute each specific sample of the material, and to represent frustrated magnets with SCQSL as the new state of matter that exhibits the properties of HF metals, with one exception: these are insulators and do not support the charge current.
Other results reported in our review paper are based on the description of low-frequency optical conductivity experimental data on herbertsmithite. For that, we implied that ZnCu 3 ( OH ) 6 Cl 2 is a system of strongly correlated fermions with properties defined primarily by a quantum spin liquid of chargeless spinons. We have also predicted the optical conductivity dependence on the magnetic field and pointed out the explicit conditions under which such a dependence can be experimentally observed. While making a step towards the confirmation of the existence of the SCQSL state in ZnCu 3 ( OH ) 6 Cl 2 , our considerations may also pave the way to studies of this state in other magnetic frustrated insulators at the nanoscale.

Author Contributions

V.R.S., V.A.S., A.Z.M., and M.Y.A. designed the project and directed it with the help of J.W.C., G.S.J. and Y.S.L. Y.S.L. performed the calculations. V.R.S., V.A.S. and J.W.C. wrote the manuscript and all authors commented on it. The manuscript reflects the contributions of all authors.

Funding

This research received no external funding.

Acknowledgments

This work was partly supported by U.S. DOE, Division of Chemical Sciences, Office of Basic Energy Sciences, Office of Energy Research. JWC acknowledges support from the McDonnell Center for the Space Sciences and expresses gratitude to the University of Madeira and its branch of Centro de Investigação em Matemática e Aplicações (CIMA) for their gracious hospitality during periods of extended residence.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pustogow, A.; Li, Y.; Voloshenko, I.; Puphal, P.; Krellner, C.; Mazin, I.I.; Dressel, M.; Valentí, R. Nature of optical excitations in the frustrated kagome compound herbertsmithite. Phys. Rev. 2017, 96, 241114(R). [Google Scholar] [CrossRef]
  2. Sushkov, A.B.; Jenkins, G.S.; Han, T.-H.; Lee, Y.S.; Drew, H.D. Infrared phonons as a probe of spin-liquid states in herbertsmithite ZnCu3(OH)6Cl2. Phys. Condens. Matter 2017, 29, 095802. [Google Scholar] [CrossRef] [PubMed]
  3. Helton, J.S.; Matan, K.; Shores, M.P.; Nytko, E.A.; Bartlett, B.M.; Qiu, Y.; Nocera, D.G.; Lee, Y.S. Dynamic Scaling in the Susceptibility of the Spin-1/2 Kagome Lattice Antiferromagnet Herbertsmithite. Phys. Rev. Lett. 2010, 104, 147201. [Google Scholar] [CrossRef] [PubMed]
  4. Helton, J.S.; Matan, K.; Shores, M.P.; Nytko, E.A.; Bartlett, B.M.; Yoshida, Y.; Takano, Y.; Suslov, A.; Qiu, Y.; Chung, J.-H.; et al. Spin Dynamics of the Spin-1/2 Kagome Lattice Antiferromagnet ZnCu3(OH)6Cl2. Phys. Rev. Lett. 2007, 98, 107204. [Google Scholar] [CrossRef] [PubMed]
  5. De Vries, M.A.; Kamenev, K.V.; Kockelmann, W.A.; Sanchez-Benitez, J.; Harrison, A. Magnetic Ground State of an Experimental S = 1/2 Kagome Antiferromagnet. Phys. Rev. Lett. 2008, 100, 157205. [Google Scholar] [CrossRef] [PubMed]
  6. Han, T.H.; Helton, J.S.; Chu, S.; Prodi, A.; Singh, D.K.; Mazzoli, C.; Müller, P.; Nocera, D.G.; Lee, Y.S. Synthesis and characterization of single crystals of the spin- 1/2 kagome-lattice antiferromagnets ZnxCu4−x(OH)6Cl2. Phys. Rev. B 2011, 83, 100402(R). [Google Scholar] [CrossRef]
  7. Han, T.H.; Chu, S.; Lee, Y.S. Refining the Spin Hamiltonian in the Spin-1/2 Kagome Lattice Antiferromagnet ZnCu3(OH)6Cl2 Using Single Crystals. Phys. Rev. Lett. 2012, 108, 157202. [Google Scholar] [CrossRef] [PubMed]
  8. Han, T.H.; Chisnell, R.; Bonnoit, C.J.; Freedman, D.E.; Zapf, V.S.; Harrison, N.; Nocera, D.G.; Takano, Y.; Lee, Y.S. Thermodynamic Properties of the Quantum Spin Liquid Candidate ZnCu3(OH)6Cl2 in High Magnetic Fields. arXiv 2014, arXiv:1402.2693. [Google Scholar]
  9. Han, T.H.; Norman, M.R.; Wen, J.-J.; Rodriguez-Rivera, J.A.; Helton, J.S.; Broholm, C.; Lee, Y.S. Correlated impurities and intrinsic spin-liquid physics in the kagome material herbertsmithite. Phys. Rev. B 2016, 94, 060409(R). [Google Scholar] [CrossRef]
  10. Paddison, J.A.M.; Daum, M.; Dun, Z.; Ehlers, G.; Liu, Y.; Stone, M.B.; Zhou, H.; Mourigal, M. Continuous excitations of the triangular-lattice quantum spin liquid YbMgGaO4. Nat. Phys. 2017, 13, 117. [Google Scholar] [CrossRef]
  11. Saito, Y.; Minamidate, T.; Kawamoto, A.; Matsunaga, N.; Nomura, K. Site-specific 13C NMR study on the locally distorted triangular lattice of the organic conductor κ − (BEDTTTF)2Cu2(CN)3. Phys. Rev. B 2018, 98, 205141. [Google Scholar] [CrossRef]
  12. Lazic, P.; Pinteric, M.; Góngora, D.R.; Pustogow, A.; Treptow, K.; Ivek, T.; Milat, O.; Gumhalter, B.; Clić, N.; Dressel, M.; et al. Importance of van der Waals interactions and cation-anion coupling in an organic quantum spin liquid. Phys. Rev. B 2018, 97, 245134. [Google Scholar] [CrossRef] [Green Version]
  13. Pinterić, M.; Gongora, D.R.; Rapljenović, Z.; Ivek, T.; Culo, M.; Korin-Hamzić, B.; Milat, O.; Gumhalter, B.; Lazić, P.; Alonso, M.S.; et al. Electrodynamics in Organic Dimer Insulators Close to Mott Critical Point. Crystals 2018, 8, 190. [Google Scholar] [CrossRef]
  14. Shi, Z.; Steinhardt, W.; Graf, D.; Corboz, P.; Weickert, F.; Harrison, N.; Jaime, M.; Marjerrison, C.; Dabkowska, H.A.; Mila, F.; et al. Emergent bound states and impurity pairs in chemically doped Shastry-Sutherland system. Nat. Commun. 2019, 10, 2439. [Google Scholar] [CrossRef]
  15. Yamaguchi, H.; Okada, M.; Kono, Y.; Kittaka, S.; Sakakibara, T.; Okabe, T.; Iwasaki, Y.; Hosokoshi, Y. Randomness-induced quantum spin liquid on honeycomb lattice. Sci. Rep. 2017, 7, 16144. [Google Scholar] [CrossRef]
  16. Shaginyan, V.R.; Amusia, M.Y.; Msezane, A.Z.; Popov, K.G. Scaling behavior of heavy fermion metals. Phys. Rep. 2010, 492, 31–109. [Google Scholar] [CrossRef] [Green Version]
  17. Amusia, M.Y.; Popov, K.G.; Shaginyan, V.R.; Stephanovich, V.A. Theory of Heavy-Fermion Compounds; Springer Series in Solid-State Sciences; Springer: Berlin/Heidelberg, Germany, 2014; Volume 182. [Google Scholar]
  18. Shaginyan, V.R.; Msezane, A.Z.; Popov, K.G. Thermodynamic properties of the kagome lattice in herbertsmithite. Phys. Rev. B 2011, 84, 060401(R). [Google Scholar] [CrossRef]
  19. Shaginyan, V.R.; Msezane, A.Z.; Popov, K.G.; Khodel, V.A. Scaling in dynamic susceptibility of herbertsmithite and heavy-fermion metals. Phys. Lett. A 2012, 376, 2622–2626. [Google Scholar] [CrossRef] [Green Version]
  20. Shaginyan, V.R.; Msezane, A.Z.; Popov, K.G.; Japaridze, G.S.; Stephanovich, V.A. Identification of strongly correlated spin liquid in herbertsmithite. Europhys. Lett. 2012, 97, 56001. [Google Scholar] [CrossRef]
  21. Shaginyan, V.R.; Msezane, A.Z.; Popov, K.G.; Japaridze, G.S.; Khodel, V.A. Heat transport in magnetic fields by quantum spin liquid in the organic insulators EtMe3Sb[Pd(dmit)2]2 and κ − (BEDTTTF)2Cu2(CN)3. Europhys. Lett. 2013, 103, 67006. [Google Scholar] [CrossRef]
  22. Liao, H.J.; Xie, Z.Y.; Chen, J.; Liu, Z.Y.; Xie, H.D.; Huang, R.Z.; Normand, B.; Xiang, T. Gapless Spin-Liquid Ground State in the S = 1/2 Kagome Antiferromagnet. Phys. Rev. Lett. 2017, 118, 137202. [Google Scholar] [CrossRef]
  23. Shaginyan, V.R.; Amusia, M.Y.; Msezane, A.Z.; Popov, K.G.; Stephanovich, V.A. Heavy fermion spin liquid in herbertsmithite. Phys. Lett. A 2015, 379, 2092–2096. [Google Scholar] [CrossRef] [Green Version]
  24. Fu, M.; Imai, T.; Han, T.H.; Lee, Y.S. Evidence for a gapped spin-liquid ground state in a kagome Heisenberg antiferromagnet. Science 2015, 350, 655–658. [Google Scholar] [CrossRef] [Green Version]
  25. Imai, T.; Fu, M.; Han, T.H.; Lee, Y.S. Local spin susceptibility of the S = 1/2 kagome lattice in ZnCu3(OD)6Cl2. Phys. Rev. B 2011, 84, 020411(R). [Google Scholar] [CrossRef]
  26. Norman, M.R. Herbertsmithite and the search for the quantum spin liquid. Rev. Mod. Phys. 2016, 88, 041002. [Google Scholar] [CrossRef]
  27. Zhou, Y.; Kanoda, K.; Ng, T.-K. Quantum spin liquid states. Rev. Mod. Phys. 2017, 89, 025003. [Google Scholar] [CrossRef]
  28. Savary, L.; Balents, L. Quantum spin liquids: A review. Rep. Prog. Phys. 2017, 80, 016502. [Google Scholar] [CrossRef]
  29. Feng, Z.; Li, Z.; Yi, X.M.W.; Wi, Y.; Zhang, J.; Wang, Y.-C.; Jiang, W.; Liu, Z.; Li, S.; Liu, F.; et al. Gapped Spin-1/2 Spinon Excitations in a New Kagome Quantum Spin Liquid Compound Cu3Zn(OH)6FBr. Chin. Phys. Lett. 2017, 34, 077502. [Google Scholar] [CrossRef]
  30. Khodel, V.A.; Shaginyan, V.R. Superfluidity in system with fermion condensate. JETP Lett. 1990, 51, 553–555. [Google Scholar]
  31. Khodel, V.A.; Shaginyan, V.R.; Khodel, V.V. New approach in the microscopic Fermi system theory. Phys. Rep. 1994, 249, 1–134. [Google Scholar] [CrossRef]
  32. Volovik, G.E. Quantum Phase Transitions from Topology in Momentum Space. In Quantum Analogues: From Phase Transitions to Black Holes and Cosmology; Springer: Berlin/Heidelberg, Germany, 2007; Volume 718, pp. 31–73. [Google Scholar]
  33. Pilon, D.V.; Lui, C.H.; Han, T.-H.; Shrekenhamer, D.; Frenzel, A.J.; Padilla, W.J.; Lee, Y.S.; Gedik, N. Spin-Induced Optical Conductivity in the Spin-Liquid Candidate Herbertsmithite. Phys. Rev. Lett. 2013, 111, 127401. [Google Scholar] [CrossRef] [Green Version]
  34. Potter, A.C.; Senthil, T.; Lee, P.A. Mechanisms for sub-gap optical conductivity in Herbertsmithite. Phys. Rev. B 2013, 87, 245106. [Google Scholar] [CrossRef] [Green Version]
  35. Shaginyan, V.R.; Stephanovich, V.A.; Msezane, A.Z.; Schuck, P.; Clark, J.W.; Amusia, M.Y.; Japaridze, G.S.; Popov, K.G.; Kirichenko, E.V. New State of Matter: Heavy Fermion Systems, Quantum Spin Liquids, Quasicrystals, Cold Gases, and High-Temperature Superconductors. J. Low Temp. Phys. 2017, 189, 410–450. [Google Scholar] [CrossRef] [Green Version]
  36. Shaginyan, V.R.; Msezane, A.Z.; Stephanovich, V.A.; Popov, K.G.; Japaridze, G.S. Universal Behavior of Quantum Spin Liquid and Optical Conductivity in the Insulator Herbertsmithite. J. Low Temp. Phys. 2018, 191, 4–13. [Google Scholar] [CrossRef]
  37. Green, D.; Santos, L.; Chamon, C. Isolated flat bands and spin-1 conical bands in two-dimensional lattices. Phys. Rev. B 2010, 82, 075104. [Google Scholar] [CrossRef]
  38. Heikkilä, T.T.; Kopnin, N.B.; Volovik, G.E. Flat bands in topological media. JETP Lett. 2011, 94, 233. [Google Scholar] [CrossRef]
  39. Volovik, G.E. A new class of normal Fermi liquids. JETP Lett. 1991, 53, 222. [Google Scholar]
  40. Kelly, Z.A.; Gallagher, M.J.; McQueen, T.M. Electron Doping a Kagome Spin Liquid. Phys. Rev. X 2016, 6, 041007. [Google Scholar] [CrossRef] [Green Version]
  41. Landau, L.D. The theory of a Fermi liquid. Sov. Phys. JETP 1956, 3, 920. [Google Scholar]
  42. Lifshitz, E.M.; Pitaevskii, L.P. Statisticheskaya Fizika (Statistical Physics); Nauka: Moscow, Russia, 1978; Pt. 2; Translated into English, Pergamon Press: Oxford, UK, 1980. [Google Scholar]
  43. Pines, D.; Noziéres, P. Theory of Quantum Liquids; Benjamin: New York, NY, USA, 1966. [Google Scholar]
  44. Coleman, P. Lectures on the Physics of Highly Correlated Electron Systems VI; Mancini, F., Ed.; American Institute of Physics: New York, NY, USA, 2002. [Google Scholar]
  45. Clark, J.W.; Zverev, M.V.; Khodel, V.A. Entropy excess in strongly correlated Fermi systems near a quantum critical point. Ann. Phys. 2012, 327, 3063–3083. [Google Scholar] [CrossRef] [Green Version]
  46. Gegenwart, P.; Tokiwa, Y.; Westerkamp, T.; Weickert, F.; Custers, J.; Ferstl, J.; Krellner, C.; Geibel, C.; Kerschl, P.; Muller, K.-H. High-field phase diagram of the heavy-fermion metal YbRh2Si2. New J. Phys. 2006, 8, 171. [Google Scholar] [CrossRef]
  47. Oeschler, N.; Hartmann, S.; Pikul, A.; Krellner, C.; Geibel, C.; Steglich, F. Low-temperature specific heat YbRh2Si2. Phys. B 2008, 403, 1254. [Google Scholar] [CrossRef]
  48. Shaginyan, V.R.; Popov, K.G.; Stephanovich, V.A.; Fomichev, V.I.; Kirichenko, E.V. High-magnetic-fields thermodynamics of the heavy-fermion metal YbRh2Si2. Europhys. Lett. 2011, 93, 17008. [Google Scholar] [CrossRef]
  49. Shaginyan, V.R.; Stephanovich, V.A.; Popov, K.G.; Kirichenko, E.V. Quasi-one-dimensional quantum spin liquid in the Cu(C4H4N2)(NO3)2 insulator. JETP Lett. 2016, 103, 32. [Google Scholar] [CrossRef]
  50. Shaginyan, V.R.; Stephanovich, V.A.; Popov, K.G.; Kirichenko, E.V.; Artamonov, S.A. Magnetic quantum criticality in quasi-one-dimensional Heisenberg antiferromagnet Cu(C4H4N2)(NO3)2. Ann. Phys. 2016, 528, 483–492. [Google Scholar] [CrossRef]
  51. Shaginyan, V.R.; Stephanovich, V.A.; Msezane, A.Z.; Japaridze, G.S.; Popov, K.G. The influence of topological phase transition on the superfluid density of overdoped copper oxides. Phys. Chem. Chem. Phys. 2017, 19, 21964. [Google Scholar] [CrossRef]
  52. Shaginyan, V.R.; Amusia, M.Y.; Clark, J.W.; Japaridze, G.S.; Msezane, A.Z.; Popov, K.G.; Stephanovich, V.A.; Zverev, M.V.; Khodel, V.A. Comment on “Topological excitations and the dynamic structure factor of spin liquids on the kagome lattice”, Nature Physics 2014, 10, 289–293. arXiv 2014, arXiv:1409.4034. [Google Scholar]
  53. Gomilsek, M.; Klanjsek, M.; Pregelj, M.; Luetkens, H.; Li, Y.; Zhang, Q.M.; Zorko, A. μSR insight into the impurity problem in quantum kagome antiferromagnets. Phys. Rev. B 2016, 94, 024438. [Google Scholar] [CrossRef]
  54. Han, T.H.; Helton, J.S.; Chu, S.; Nocera, D.G.; Rodriguez-Rivera, J.A.; Broholm, C.; Lee, Y.S. Fractionalized excitations in the spin-liquid state of a kagome-lattice antiferromagnet. Nature 2012, 492, 406–410. [Google Scholar] [CrossRef] [Green Version]
  55. Knafo, W.; Raymond, S.; Flouquet, J.; Fåk, B.; Adams, M.A.; Haen, P.; Lapierre, F.; Yates, S.; Lejay, P. Anomalous scaling behavior of the dynamical spin susceptibility of Ce0.925La0.075Ru2Si2. Phys. Rev. B 2004, 70, 174401. [Google Scholar] [CrossRef]
  56. Fåk, B.; Coomer, F.C.; Harrison, A.; Visser, D.; Zhitomirsky, M.E. Spin-liquid behavior in a kagome antiferromagnet: Deuteronium jarosite. Europhys. Lett. 2008, 81, 17006. [Google Scholar] [CrossRef]
  57. Stock, C.; Broholm, C.; Demmel, F.; Duijn, J.V.; Taylor, J.W.; Kang, H.J.; Hu, R.; Petrovic, C. From Incommensurate Correlations to Mesoscopic Spin Resonance in YbRh2Si2. Phys. Rev. Lett. 2012, 109, 127201. [Google Scholar] [CrossRef] [PubMed]
  58. Shaginyan, V.R.; Msezane, A.Z.; Japaridze, G.S.; Popov, K.G.; Khodel, V.A. Strongly correlated Fermi systems as a new state of matter. Front. Phys. 2016, 11, 117103. [Google Scholar] [CrossRef] [Green Version]
  59. Yamashita, M.; Nakata, N.; Senshu, Y.; Nagata, M.; Yamamoto, H.M.; Kato, R.; Shibauchi, T.; Matsuda, Y. Highly Mobile Gapless Excitations in a Two-Dimensional Candidate Quantum Spin Liquid. Science 2010, 328, 1246–1248. [Google Scholar] [CrossRef] [PubMed]
  60. Yamashita, M.; Shibauchi, T.; Matsuda, Y. Probing non-equilibrium vibrational relaxation pathways in Qunatum Spin Liquids. Chem. Phys. 2012, 13, 74. [Google Scholar]
  61. Imai, T.; Nytko, E.A.; Bartlett, B.M.; Shores, M.P.; Nocera, D.G. 63Cu, 35Cl, and 1H NMR in the S = 1/2 Kagome Lattice ZnCu3(OH)6Cl2. Phys. Rev. Lett. 2008, 100, 077203. [Google Scholar] [CrossRef] [PubMed]
  62. Carretta, P.; Pasero, R.; Giovannini, M.; Baines, C. Magnetic-field-induced crossover from non-Fermi to Fermi liquid at the quantum critical point of YbCu5−xAux. Phys. Rev. B 2009, 79, 020401(R). [Google Scholar] [CrossRef]
  63. Gegenwart, P.; Westerkamp, T.; Krellner, C.; Tokiwa, Y.; Paschen, S.; Geibel, C.; Steglich, F.; Abrahams, E.; Si, Q. Multiple energy scales at a quantum critical point. Science 2007, 315, 969–971. [Google Scholar] [CrossRef] [PubMed]
  64. Shen, Y.; Li, Y.-D.; Wo, H.; Li, Y.; Shen, S.; Pan, B.; Wang, Q.; Walker, H.C.; Steffens, P.; Boehm, M.; et al. Evidence for a spinon Fermi surface in a triangular-lattice quantum-spin-liquid candidate. Nature 2016, 540, 559–562. [Google Scholar] [CrossRef] [Green Version]
  65. Luo, Z.-X.; Lake, E.; Mei, J.-W.; Starykh, O.A. Spinon Magnetic Resonance of Quantum Spin Liquids. Phys. Rev. Lett. 2018, 120, 037204. [Google Scholar] [CrossRef] [Green Version]
  66. Xu, Y.; Zhang, J.; Li, Y.S.; Yu, Y.J.; Hong, X.C.; Zhang, Q.M.; Li, S.Y. Absence of Magnetic Thermal Conductivity in the Quantum Spin-Liquid Candidate YbMgGaO4. Phys. Rev. Lett. 2016, 117, 267202. [Google Scholar] [CrossRef] [Green Version]
  67. Nakayama, M.; Kondo, T.; Tian, Z.; Ishikawa, J.J.; Halim, M.; Bareille, C.; Malaeb, W.; Kuroda, K.; Tomita, T.; Ideta, S.; et al. Slater to Mott Crossover in the Metal to Insulator Transition of Nd2Ir2O7. Phys. Rev. Lett. 2016, 117, 056403. [Google Scholar] [CrossRef]
  68. Yan, S.; Huse, D.A.; White, S.R. Spin-Liquid Ground State of the S = 1/2 Kagome Heisenberg Antiferromagnet. Science 2011, 332, 1173–1176. [Google Scholar] [CrossRef]
  69. Depenbrock, S.; McCulloch, I.P.; Schollwock, U. Nature of the Spin-Liquid Ground State of the S = 1/2 Heisenberg Model on the Kagome Lattice. Phys. Rev. Lett. 2012, 109, 067201. [Google Scholar] [CrossRef]
  70. Jiang, H.-C.; Wang, Z.; Balents, L. Identifying topological order by entanglement entropy. Nat. Phys. 2012, 8, 902. [Google Scholar] [CrossRef]
  71. Zorko, A.; Herak, M.; Gomilŝek, M.; van Tol, J.; Velàzquez, M.; Khuntia, P.; Bert, F.; Mendels, P. Symmetry Reduction in the Quantum Kagome Antiferromagnet Herbertsmithite. Phys. Rev. Lett. 2017, 118, 017202. [Google Scholar] [CrossRef] [PubMed]
  72. Pustogow, A.; Saito, Y.; Zhukova, E.; Gorshunov, B.; Kato, R.; Lee, T.-H.; Dobrosavljević, V.; Dressel, M. Low-Energy Excitations in Quantum Spin Liquids Identified by Optical Spectroscopy. Phys. Rev. Lett. 2018, 121, 056402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Kagome lattice and flat band. Panel (a) Crystal structure of ZnCu 3 ( OH ) 6 Cl 2 along the hexagonal c axis. The blue spheres are Cu sites. Arrows show Cu spins in the frustrated configuration. The kagome planes with frustrated spins, occupying a flat band (panel (b)), are shown by two triangles and highlighted by the dashed lines, displaying the hexagon. Groups of OH , Zn , and Cl atoms are shown schematically. Zn atoms may be situated below or above Cu planes. Panel (b) Calculated single particle spectrum ε ( p / p F , T / E F ) as a function of dimensionless momentum k = p / p F and temperature T / E F , with E F being the Fermi level [16,17]. The cases of almost zero ( T = 10 4 E F ; black curve) and finite ( T = 10 2 E F ; red curve) temperatures are reported. It is seen that the dispersionless part of the spectrum vanishes at finite temperatures, making the situation almost indistinguishable from normal Fermi liquid, while at T 0 the band becomes almost flat at k i < k < k f , where k i = p i / p F and k f = p f / p F .
Figure 1. Kagome lattice and flat band. Panel (a) Crystal structure of ZnCu 3 ( OH ) 6 Cl 2 along the hexagonal c axis. The blue spheres are Cu sites. Arrows show Cu spins in the frustrated configuration. The kagome planes with frustrated spins, occupying a flat band (panel (b)), are shown by two triangles and highlighted by the dashed lines, displaying the hexagon. Groups of OH , Zn , and Cl atoms are shown schematically. Zn atoms may be situated below or above Cu planes. Panel (b) Calculated single particle spectrum ε ( p / p F , T / E F ) as a function of dimensionless momentum k = p / p F and temperature T / E F , with E F being the Fermi level [16,17]. The cases of almost zero ( T = 10 4 E F ; black curve) and finite ( T = 10 2 E F ; red curve) temperatures are reported. It is seen that the dispersionless part of the spectrum vanishes at finite temperatures, making the situation almost indistinguishable from normal Fermi liquid, while at T 0 the band becomes almost flat at k i < k < k f , where k i = p i / p F and k f = p f / p F .
Condensedmatter 04 00075 g001
Figure 2. (Color online). Schematic strongly correlated quantum spin liquid (SCQSL) phase diagram in the “temperature–magnetic field” representation. The magnetic field B is the independent variable (control parameter). Vertical and horizontal arrows show Landau Fermi liquid (LFL)–non-Fermi-liquid (NFL) and NFL-LFL transitions at fixed B and T, respectively. The dependences M * ( T ) and M * ( B ) are given by Equations (5) and (6), respectively. The hatched area represents the transition region at T = T * ( B ) (see Equation (10)). The transitions occur according to the directions of the blue arrows. The solid line in the hatched area represents the function T * ( B ) T M ( B ) given by Equation (8). The functions W ( B ) T T * and T W ( B ) T T * shown by two-headed arrows define the width of the NFL state and the transition area, respectively. At the fermion condensation quantum phase transition (FCQPT) indicated by the arrow, the effective mass M * diverges, and both W ( B ) and T W ( B ) tend to zero.
Figure 2. (Color online). Schematic strongly correlated quantum spin liquid (SCQSL) phase diagram in the “temperature–magnetic field” representation. The magnetic field B is the independent variable (control parameter). Vertical and horizontal arrows show Landau Fermi liquid (LFL)–non-Fermi-liquid (NFL) and NFL-LFL transitions at fixed B and T, respectively. The dependences M * ( T ) and M * ( B ) are given by Equations (5) and (6), respectively. The hatched area represents the transition region at T = T * ( B ) (see Equation (10)). The transitions occur according to the directions of the blue arrows. The solid line in the hatched area represents the function T * ( B ) T M ( B ) given by Equation (8). The functions W ( B ) T T * and T W ( B ) T T * shown by two-headed arrows define the width of the NFL state and the transition area, respectively. At the fermion condensation quantum phase transition (FCQPT) indicated by the arrow, the effective mass M * diverges, and both W ( B ) and T W ( B ) tend to zero.
Condensedmatter 04 00075 g002
Figure 3. (Color online) The dependence χ ( T ) , measured experimentally on ZnCu 3 ( OH ) 6 Cl 2 , taken from [3]. Magnetic fields are coded by symbols and reported in the legend. The example of χ max and T max at B = 3 T is shown. The result of calculation at B = 0 (full line) shows the dependence χ ( T ) T 2 / 3 [17,18,21].
Figure 3. (Color online) The dependence χ ( T ) , measured experimentally on ZnCu 3 ( OH ) 6 Cl 2 , taken from [3]. Magnetic fields are coded by symbols and reported in the legend. The example of χ max and T max at B = 3 T is shown. The result of calculation at B = 0 (full line) shows the dependence χ ( T ) T 2 / 3 [17,18,21].
Condensedmatter 04 00075 g003
Figure 4. Normalized susceptibility χ N = χ / χ max = M N * versus normalized temperature T N (see Equation (7)) of ZnCu 3 ( OH ) 6 Cl 2 [3] as collected in Figure 3. We normalize specific heat ( C e l / T ) N = M N * taken from the experimental results on C e l / T in YbRh 2 Si 2 in magnetic fields B shown in Figure 5 (b) [46]. The corresponding values of B (legends) are coded by symbols. The full curve marks our theoretical calculations at B B * when the quasiparticle band is fully polarized. [17,20].
Figure 4. Normalized susceptibility χ N = χ / χ max = M N * versus normalized temperature T N (see Equation (7)) of ZnCu 3 ( OH ) 6 Cl 2 [3] as collected in Figure 3. We normalize specific heat ( C e l / T ) N = M N * taken from the experimental results on C e l / T in YbRh 2 Si 2 in magnetic fields B shown in Figure 5 (b) [46]. The corresponding values of B (legends) are coded by symbols. The full curve marks our theoretical calculations at B B * when the quasiparticle band is fully polarized. [17,20].
Condensedmatter 04 00075 g004
Figure 5. (Color online) Electronic specific heat C e l / T of the archetypical heavy fermion (HF) metal YbRh 2 Si 2 . Panel (a) reports the temperature dependence of the electronic specific heat C e l / T of YbRh 2 Si 2 at high magnetic fields [46] shown in the legend. The electronic specific heat C e l / T of YbRh 2 Si 2 strongly resembles C mag / T of ZnCu 3 ( OH ) 6 Cl 2 , shown in Figure 6. Panel (b): The normalized specific heat C e l / T at high [46] and low magnetic fields [47] extracted from the specific heat ( C / T ) measurements on the YbRh 2 Si 2 . The low-field calculations are depicted by the short-dash curve tracing the scaling behavior of M N * (see Equation (7)). Our high-field calculations (solid line) are taken at high magnetic field B, at which the quasiparticle band becomes fully polarized [48].
Figure 5. (Color online) Electronic specific heat C e l / T of the archetypical heavy fermion (HF) metal YbRh 2 Si 2 . Panel (a) reports the temperature dependence of the electronic specific heat C e l / T of YbRh 2 Si 2 at high magnetic fields [46] shown in the legend. The electronic specific heat C e l / T of YbRh 2 Si 2 strongly resembles C mag / T of ZnCu 3 ( OH ) 6 Cl 2 , shown in Figure 6. Panel (b): The normalized specific heat C e l / T at high [46] and low magnetic fields [47] extracted from the specific heat ( C / T ) measurements on the YbRh 2 Si 2 . The low-field calculations are depicted by the short-dash curve tracing the scaling behavior of M N * (see Equation (7)). Our high-field calculations (solid line) are taken at high magnetic field B, at which the quasiparticle band becomes fully polarized [48].
Condensedmatter 04 00075 g005
Figure 6. (Color online) Specific heat C mag / T measured on powder [4,5] and single-crystal [6,7,8] herbertsmithite samples as a function of temperature. Magnetic fields (legend) are coded by symbols. It is clearly seen that powder and single-crystal samples have approximately the same C mag / T . The example of C mag / T max and T max at B 3 T is shown.
Figure 6. (Color online) Specific heat C mag / T measured on powder [4,5] and single-crystal [6,7,8] herbertsmithite samples as a function of temperature. Magnetic fields (legend) are coded by symbols. It is clearly seen that powder and single-crystal samples have approximately the same C mag / T . The example of C mag / T max and T max at B 3 T is shown.
Condensedmatter 04 00075 g006
Figure 7. (Color online) Normalized specific heat ( C mag / T ) N versus normalized temperature T N at different magnetic fields (legend) [17,23] extracted from specific heat C mag / T measured on powder [4,5] and single-crystal [6,7,8] herbertsmithite samples (see Figure 6). The theoretical result from [17,18], represented by the solid curve, traces the scaling behavior of the normalized specific heat.
Figure 7. (Color online) Normalized specific heat ( C mag / T ) N versus normalized temperature T N at different magnetic fields (legend) [17,23] extracted from specific heat C mag / T measured on powder [4,5] and single-crystal [6,7,8] herbertsmithite samples (see Figure 6). The theoretical result from [17,18], represented by the solid curve, traces the scaling behavior of the normalized specific heat.
Condensedmatter 04 00075 g007
Figure 8. (Color online) Panel (a): The temperatures T M ( B ) at which the maxima of χ m a x ( B ) M M * ( B ) (see Figure 3) occur. The solid line represents the function T M B (see Equation (8)). Panel (b): the maxima χ m a x ( B ) of χ ( T ) versus B (see Figure 3). The solid curve is given by Equation (6).
Figure 8. (Color online) Panel (a): The temperatures T M ( B ) at which the maxima of χ m a x ( B ) M M * ( B ) (see Figure 3) occur. The solid line represents the function T M B (see Equation (8)). Panel (b): the maxima χ m a x ( B ) of χ ( T ) versus B (see Figure 3). The solid curve is given by Equation (6).
Condensedmatter 04 00075 g008
Figure 9. (Color online) Normalized dynamic spin susceptibility ( T 2 / 3 χ ) N , showing scaling properties. The data are taken from those on the HF metal Ce 0.925 La 0.075 Ru 2 Si 2 [55]. The full line is defined by Equation (12) [19].
Figure 9. (Color online) Normalized dynamic spin susceptibility ( T 2 / 3 χ ) N , showing scaling properties. The data are taken from those on the HF metal Ce 0.925 La 0.075 Ru 2 Si 2 [55]. The full line is defined by Equation (12) [19].
Condensedmatter 04 00075 g009
Figure 10. (Color online) Same as in Figure 9 but for herbertsmithite ZnCu 3 ( OH ) 6 Cl 2 [3]. The full line is once more defined by Equation (12) [17,19].
Figure 10. (Color online) Same as in Figure 9 but for herbertsmithite ZnCu 3 ( OH ) 6 Cl 2 [3]. The full line is once more defined by Equation (12) [17,19].
Condensedmatter 04 00075 g010
Figure 11. (Color online) Same as in Figure 9 but for deuteronium jarosite ( D 3 O ) Fe 3 ( SO 4 ) 2 ( OD ) 6 [56]. The full line is once more defined by Equation (12) [17,19].
Figure 11. (Color online) Same as in Figure 9 but for deuteronium jarosite ( D 3 O ) Fe 3 ( SO 4 ) 2 ( OD ) 6 [56]. The full line is once more defined by Equation (12) [17,19].
Condensedmatter 04 00075 g011
Figure 12. (Color online) Normalized dynamic spin susceptibility T χ against E = ω / k B T . Data are taken from those on YbRh 2 Si 2 [57]. The full line is defined by Equation (14) [17,19].
Figure 12. (Color online) Normalized dynamic spin susceptibility T χ against E = ω / k B T . Data are taken from those on YbRh 2 Si 2 [57]. The full line is defined by Equation (14) [17,19].
Condensedmatter 04 00075 g012
Figure 13. (Color online) Main panel (A). Magnetic field dependence of normalized spin-lattice relaxation rate ( 1 / T 1 T ) N . Data for different substances are shown by different markers: solid squares correspond to ZnCu 3 ( OH ) 6 Cl 2 [61] and solid triangles to YbCu 5 x Au x at x = 0.4 [62]. The normalization is done in the inflection point, shown by the arrow. Inset (B). Same as in main panel (A), but for the normalized magnetoresistance ρ N , for YbRh 2 Si 2 at different temperatures [63] (legend) coded by symbols. The results of calculations are shown in both panels by the full line, defining the scaling behavior of W r ( M * ) 2 (see Equations (15) and (16)).
Figure 13. (Color online) Main panel (A). Magnetic field dependence of normalized spin-lattice relaxation rate ( 1 / T 1 T ) N . Data for different substances are shown by different markers: solid squares correspond to ZnCu 3 ( OH ) 6 Cl 2 [61] and solid triangles to YbCu 5 x Au x at x = 0.4 [62]. The normalization is done in the inflection point, shown by the arrow. Inset (B). Same as in main panel (A), but for the normalized magnetoresistance ρ N , for YbRh 2 Si 2 at different temperatures [63] (legend) coded by symbols. The results of calculations are shown in both panels by the full line, defining the scaling behavior of W r ( M * ) 2 (see Equations (15) and (16)).
Condensedmatter 04 00075 g013

Share and Cite

MDPI and ACS Style

Shaginyan, V.R.; Msezane, A.Z.; Amusia, M.Y.; Clark, J.W.; Japaridze, G.S.; Stephanovich, V.A.; Leevik, Y.S. Thermodynamic, Dynamic, and Transport Properties of Quantum Spin Liquid in Herbertsmithite from an Experimental and Theoretical Point of View. Condens. Matter 2019, 4, 75. https://doi.org/10.3390/condmat4030075

AMA Style

Shaginyan VR, Msezane AZ, Amusia MY, Clark JW, Japaridze GS, Stephanovich VA, Leevik YS. Thermodynamic, Dynamic, and Transport Properties of Quantum Spin Liquid in Herbertsmithite from an Experimental and Theoretical Point of View. Condensed Matter. 2019; 4(3):75. https://doi.org/10.3390/condmat4030075

Chicago/Turabian Style

Shaginyan, Vasily R., Alfred Z. Msezane, Miron Ya. Amusia, John W. Clark, George S. Japaridze, Vladimir A. Stephanovich, and Yulya S. Leevik. 2019. "Thermodynamic, Dynamic, and Transport Properties of Quantum Spin Liquid in Herbertsmithite from an Experimental and Theoretical Point of View" Condensed Matter 4, no. 3: 75. https://doi.org/10.3390/condmat4030075

APA Style

Shaginyan, V. R., Msezane, A. Z., Amusia, M. Y., Clark, J. W., Japaridze, G. S., Stephanovich, V. A., & Leevik, Y. S. (2019). Thermodynamic, Dynamic, and Transport Properties of Quantum Spin Liquid in Herbertsmithite from an Experimental and Theoretical Point of View. Condensed Matter, 4(3), 75. https://doi.org/10.3390/condmat4030075

Article Metrics

Back to TopTop