Next Article in Journal
Approach to the Fatigue and Cellular Behavior of Superficially Modified Porous Titanium Dental Implants
Previous Article in Journal
Selective Laser Melting of CuSn10: Simulation of Mechanical Properties, Microstructure, and Residual Stresses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Strongly Correlated Quantum Spin Liquids versus Heavy Fermion Metals: A Review

by
Vasily R. Shaginyan
1,2,*,
Alfred Z. Msezane
2,
George S. Japaridze
2,
Stanislav A. Artamonov
1 and
Yulya S. Leevik
3
1
Petersburg Nuclear Physics Institute, NRC Kurchatov Institute, 188300 Gatchina, Russia
2
Clark Atlanta University, Atlanta, GA 30314, USA
3
National Research University Higher School of Economics, 194100 St. Petersburg, Russia
*
Author to whom correspondence should be addressed.
Materials 2022, 15(11), 3901; https://doi.org/10.3390/ma15113901
Submission received: 3 April 2022 / Revised: 15 May 2022 / Accepted: 25 May 2022 / Published: 30 May 2022
(This article belongs to the Special Issue Exploration of Novel Quantum Spin Liquid Materials)

Abstract

:
This review considers the topological fermion condensation quantum phase transition (FCQPT) that explains the complex behavior of strongly correlated Fermi systems, such as frustrated insulators with quantum spin liquid and heavy fermion metals. The review contrasts theoretical consideration with recent experimental data collected on both heavy fermion metals (HF) and frustrated insulators. Such a method allows to understand experimental data. We also consider experimental data collected on quantum spin liquid in Lu 3 Cu 2 Sb 3 O 14 and quasi-one dimensional (1D) quantum spin liquid in both YbAlO 3 and Cu ( C 4 H 4 N 2 ) ( NO 3 ) 2 with the aim to establish a sound theoretical explanation for the observed scaling laws, Landau Fermi liquid (LFL) and non-Fermi-liquid (NFL) behavior exhibited by these frustrated insulators. The recent experimental data on the heavy-fermion metal α YbAl 1 x Fe x B 4 , with x = 0.014 , and on its sister compounds β YbAlB 4 and YbCo 2 Ge 4 , carried out under the application of magnetic field as a control parameter are analyzed. We show that the thermodynamic and transport properties as well as the empirical scaling laws follow from the fermion condensation theory. We explain how both the similarity and the difference in the thermodynamic and transport properties of α YbAl 1 x Fe x B 4 and in its sister compounds β YbAlB 4 and YbCo 2 Ge 4 emerge, as well as establish connection of these (HF) metals with insulators Lu 3 Cu 2 Sb 3 O 14 , Cu ( C 4 H 4 N 2 ) ( NO 3 ) 2 and YbAlO 3 . We demonstrate that the universal LFL and NFL behavior emerge because the HF compounds and the frustrated insulators are located near the topological FCQPT or are driven by the application of magnetic fields.

1. Introduction

Strongly correlated Fermi systems are represented by numerous heavy fermion (HF) compounds characterized by their diverse microscopic properties. To study these strongly correlated Fermi systems, we employ a topological symmetry representing a powerful method for gaining knowledge about physical systems spanning from solids to galaxies and their clusters in the Universe [1,2]. Understanding of such symmetry and conditions for its violation allows one to obtain a general information about physical systems. The low-temperature universal properties of strongly correlated systems, including HF metals and frustrated insulators, can be unveiled within the fermion condensation (FC) theory [1,2,3,4,5,6,7]. This universality suggests that strongly correlated systems, or HF compounds, represent a new state of matter. It means that this new state is independent of the atomic composition of HF compounds, exhibits universal properties, and is defined by the formation of flat or approximately flat bands [1,2,3,4,5,6,7]. These bands, predicted many years ago [3,5,6] and discovered recently in graphene, see, e.g., [8,9,10], originate from a specific quantum phase transition known as the topological fermion-condensation quantum phase transition (FCQPT) that rearranges the Fermi surface into the Fermi volume, generating a flat band. Thus, for very different substances and under very different external conditions the universal topological FCQPT occurs at microscopic level, determining the macroscopic properties and universal behavior of HF compounds. These compounds proliferated and they include HF metals, quantum spin liquids, quasicrystals and two dimensional systems like 3 He [1]. These HF compounds represent the new state of matter, since their behavior near the topological FCQPT acquires important similarities that make them universal [1,6,7,11,12,13,14].
In our brief review we consider recent experimental data collected on the frustrated insulator like Lu 3 Cu 2 Sb 3 O 14 and quasi-1D quantum spin liquid (1DQSL) in both YbAlO 3 and Cu ( C 4 H 4 N 2 ) ( NO 3 ) 2 and family of HF metals α YbAl 1 x Fe x B 4 , β YbAlB 4 and YbCo 2 Ge 4 . We show that these unexplained experimental data can be explained within the framework of fermion condensation (FC) theory based on the topological FCQPT. As a result, we demonstrate that these HF compounds belong to the new state of matter. Thus, the recent experimental data support our conclusion expressed in recent reviews [13,14] that HF compounds form the new state of matter.
Currently, numerous quantum spin liquids (QSLs) with various types of ground states are proposed [14,15,16,17,18,19,20,21,22,23,24,25]. These QSLs define the thermodynamic, transport and relaxation properties of frustrated insulators and represent the new state of matter formed by HF compounds [1,13,14]. QSLs are formed with fermionic quasiparticles with the effective mass M * which are called spinons. Spinons carry spin σ = 1 / 2 and no charge. At temperature T = 0 the Fermi sphere is shaped from spinons with the Fermi momentum p F . Thus, frustrated insulators can be viewed as a spinon metal, which differs from HF metals in that it cannot support the electric current.
The Fermi sphere of spinons can be located near the topological Fermion condensation phase transition (FCQPT) that forms the FC state and the corresponding flat band [1,3,4,7,11,26]. In the FC state, at T = 0 the corresponding flat band is given by the equation
ε ( p , T = 0 ) = μ , p i p F p f ; 0 n ( p ) 1 .
where p i and p f stand for initial and final momenta, where the flat band reside. At T > 0 the quasiparticle occupation numbers n ( p ) is given by the Fermi–Dirac distribution function which is represented in the form [5,27]
ε ( p , T ) μ ( T ) = T ln 1 n ( p , T ) n ( p , T ) .
Taking into account that T 0 , the distribution function satisfies the inequality 0 < n ( p ) < 1 for p i p F p f , we see that the logarithm is finite; the right hand side of Equation (2) vanishes, and lead to Equation (1). Near FCQPT flat band takes place and the notion of the strongly correlated quantum spin liquid (SCQSL) emerges that allows one to describe numerous data related to the thermodynamic, relaxation and transport properties of frustrated magnetic insulators [1,7,11,14,22,28,29,30].
We consider recent measurements in magnetic field B of the quantum spin liquid [23] that forms the thermodynamic properties of Lu 3 Cu 2 Sb 3 O 14 and the thermodynamic of 1DQSL that defines behavior of YbAlO 3 [24] and Cu ( C 4 H 4 N 2 ) ( NO 3 ) 2 [22,25], see Section 2. We analyze recently obtained measurements on the heavy-fermion (HF) metals β YbAlB 4 [31,32], YbCo 2 Ge 4 [33] and α YbAl 1 x Fe x B 4 [34], that have been performed under the application of magnetic field B, as well as on β YbAlB 4 under the application of hydrostatic pressure P [31,32]. These have received substantial theoretical analysis [31,32,35,36,37,38,39,40,41,42,43,44]. We explain these results within the framework of the FC theory, and show that the mentioned above HF metals exhibit the same behavior as that of SCQSL, forming properties of Lu 3 Cu 2 Sb 3 O 14 , YbAlO 3 and Cu ( C 4 H 4 N 2 ) ( NO 3 ) 2 , see Section 6. Our results are summarized in Section 8.

2. Universal Scaling Behavior of Quantum Spin Liquid

The ground state energy of QSL depends weakly on the spins configuration, since the spinons of the triangular lattice compounds form symmetric positions. Therefore, the triangular lattice is near to a topologically protected flat band of the spectrum with zero excitation energy [6,7,28,45,46,47]. As a result, the topological FCQPT can be considered as a quantum critical point (QCP) of the Lu 3 Cu 2 Sb 3 O 14 quantum spin liquid. In that case the elementary magnetic quasiparticles, dubbed spinons, defining the relaxation, transport and thermodynamic properties, carry the effective mass M * , zero charge and spin σ = 1 / 2 . Spinons occupy the corresponding Fermi sphere with the Fermi momentum p F , representing HF quasiparticles of deconfined SCQSL. Spinon quasiparticles generate the excitation spectrum typical for HF compounds located near the topological FCQPT. The ground state energy E ( n ) is given by the Landau functional, depending on the spinon distribution function n σ ( p ) , where p is the momentum. Near the FCQPT point, the spinon effective mass M * is defined by the Landau Equation [6,48]
1 M * ( T , B ) = 1 M * ( T = 0 , B = 0 ) + 1 p F 2 σ 1 p F p 1 p F F σ , σ 1 ( p F , p 1 ) δ n σ 1 ( p 1 ) p 1 d p 1 ( 2 π ) 3 .
In Equation (3) B is magnetic field and we rewrite the spinon distribution function as δ n σ ( p ) n σ ( p , T , B ) n σ ( p , T = 0 , B = 0 ) . Note that both functional E ( n ) and Equation (3) are exact [1,49]. This fact provides firm ground to construct the theory of HF compounds [1,6,7]. The geometric frustration of QSL in Lu 3 Cu 2 Sb 3 O 14 [23] is located near the topological FCQPT, therefore we employ the theory of HF compounds, that is the FC theory, to describe SCQSL of Lu 3 Cu 2 Sb 3 O 14 , see, e.g., [6,28]. This theory allows quantitative analysis of the thermodynamic, relaxation and transport properties of both HF compounds containing QSL and HF metals [1,6,7,11,14,28]. We will show that the thermodynamic properties of Lu 3 Cu 2 Sb 3 O 14 coincide with those of the frustrated magnet ZnCu 3 ( OH ) 6 Cl 2 and of HF metals including the archetypical HF metal YbRh 2 Si 2 [45].
In Equation (3) the only role of the Landau interaction is to drive the system to the FCQPT point, where the Fermi surface changes its topology so that the effective mass acquires strong temperature and field dependences [1,6,50,51], as seen from the inset of Figure 1. Indeed, it is seen from the inset that the effective mass M * ( T , B ) C mag / T strongly depends on both T and B. At the topological FCQPT the term 1 / M * ( T = 0 , B = 0 ) vanishes and Equation (3) becomes homogeneous and therefore is solved analytically. We remark that at p = p F the single-particle spectrum ε ( p ) acquires an inflection point at which the corresponding flat band related to FC is formed [3,6]. In the simplest case the inflection point turns out to be ε ( p ) = ( p p F ) 3 [5,6]; the other types of the inflection points are considered in Section 6. At B = 0 , the effective mass depends on T exhibiting the non-Fermi liquid (NFL) behavior [6]
M * ( T ) a T T 2 / 3 .
At finite T magnetic field B drives the system to the Landau Fermi liquid (LFL) behavior with
M * ( B ) a B B 2 / 3 .
Here a T and a B are fitting parameters. Note that the exponent 2 / 3 corresponds to the above considered inflection point.
Figure 1. The normalized specific heat ( C mag / T ) N as a function of the normalized temperature T N measured under the application of magnetic field is shown in the legend. The normalized specific heat C mag / T is extracted from the measurement of the specific heat of Lu 3 Cu 2 Sb 3 O 14 [23] shown in the inset. The solid orange curve displaces the theoretical calculations based on Equation (3) [45]. The same curve is shown in Figure 2 and Figure 3, exhibiting the scaling of the thermodynamic properties of the quantum spin liquid of Lu 3 Cu 2 Sb 3 O 14 .
Figure 1. The normalized specific heat ( C mag / T ) N as a function of the normalized temperature T N measured under the application of magnetic field is shown in the legend. The normalized specific heat C mag / T is extracted from the measurement of the specific heat of Lu 3 Cu 2 Sb 3 O 14 [23] shown in the inset. The solid orange curve displaces the theoretical calculations based on Equation (3) [45]. The same curve is shown in Figure 2 and Figure 3, exhibiting the scaling of the thermodynamic properties of the quantum spin liquid of Lu 3 Cu 2 Sb 3 O 14 .
Materials 15 03901 g001
The universal scaling of the effective mass M * is shown in Figure 4. This behavior is given by Equation (3), provided that the system is located near FCQPT. At finite B and T, near the topological FCQPT, the solution of Equation (3) M * ( B , T ) links the LFL ( M * ( T ) c o n s t ) and NFL ( M * ( T ) T 2 / 3 ) regions [1,6,7]. As seen from Figure 4, the LFL behavior and the NFL one, given by Equations (4) and (5) are separated by the crossover region at which M * reaches its maximum value M M * at temperature T M . It is seen from Figure 4, representing the universal scaling behavior of the dimensionless normalized effective mass M N * = M * / M M as a function of the dimensionless normalized temperature T N = T / T M , that M N * exhibits the usual behavior of the experimental thermodynamic function like C / T or χ ( T ) , see, e.g., [6]. Indeed, the region T N 1 represents the crossover region between the LFL behavior with almost constant effective mass and the NFL behavior, exhibiting the M * T 2 / 3 dependence, see Equation (4), and T i n f separates the beginning point of the crossover region from the NFL region [6]. Note that both M M * and T M , occurring at M N * = T N = 1 , depend on the microscopic properties of the system in question [6], while their normalized values exhibit the universal scaling. Thus, the FC theory, based on the topological FCQPT and the corresponding flat bands, incorporates the inherent universal scaling behavior that is experimentally exhibited by numerous HF compounds [1,6,7].

2.1. Universal Behavior of L u 3 C u 2 S b 3 O 14

It is seen from the inset of Figure 1 that C mag / T reaches its maximum value ( C mag / T ) max ( B ) under the application of magnetic fields at some temperature T max ( B ) . To reveal the scaling, we introduce the dimensionless normalized specific heat ( C mag / T ) N as a function of the dimensionless normalized temperature T N = T / T max ( B ) [6]
( C mag / T ) N = C mag / T ( C mag / T ) max = M N * .
From Equation (6) and from both Figure 1 and Figure 4, we see that ( C mag / T ) N ( T N ) exhibits the universal, scaling as a function of only T N T / B , presented by a single curve.
To construct the interpolating equation revealing the universal scaling of the effective mass M * C mag / T , we employ both the dimensionless normalized effective mass M N * and the dimensionless normalized temperature T N , defined by dividing the effective mass M * ( T , B ) by its maximal values, M max * ( T , B ) , and temperature T by T max at which the maximum M N * occurs, T N = T / T max [6]. Magnetic field B emerges in Equation (3) as the combination μ B B / k B T . As a result, k B T max μ B B where k B is the Boltzmann constant and μ b is the Bohr magneton [6,50]. This observation allows us to conlcude have that [6,7,44]
T max B
and
T N T / B
Thus, we obtain that the normalized effective mass M N * = M * / M max * = ( C mag / T ) N is well defined by the interpolating function, approximating the solutions of Equation (3) [6]
M N * ( y ) c 0 1 + c 1 y 2 1 + c 2 y 8 / 3 .
Here c 0 = ( 1 + c 2 ) / ( 1 + c 1 ) , c 1 and c 2 are fitting parameters, and y = T / T max T / B . Clearly, from both Equations (5) and (9) that under the application of magnetic field M * becomes finite and at low temperatures the system exhibits the LFL behavior, C mag ( T ) / T M * ( T ) M * ( T = 0 ) + a 1 T 2 . As seen from the inset of Figure 1, at increasing temperatures M * C mag / T increases and enters the crossover region, reaching its maximum M max * ( C mag ( T ) / T ) max at T = T max , with subsequent diminishing given by Equations (4) and (9). Scaling behavior is manifested by Equation (9), exhibiting the superior quality of Equation (3) at the topological FCQPT: the function M * ( T , B ) of two variables transforms into the function M * of the single variable T N T / B . We employ Equation (9) to outline the universal scaling, verifying our calculations based on Equation (3).
The scaling of ( C mag / T ) N = M N , extracted from the experimental data C mag ( T , B ) / T [23], is reported in Figure 1. The data for a wide range of B vales up to 9 T merge well into a single curve. Figure 2 reports the normalized specific heat ( C el / T ) N = M N * of YbRh 2 Si 2 versus normalized temperature T N as a function of B. Indeed, at low T N 0.1 the normalized specific heat ( C el / T ) N 0.4 [1,29]. This value is determined by the polarization of the heavy electron band under the application of magnetic fields B > 4 T, and coincides with that of χ N = M N * obtained on ZnCu 3 ( OH ) 6 Cl 2 and shown in Figure 3. Results of our calculations are represented by the same solid curve in Figure 1, Figure 2 and Figure 3. We stress that at low normalized temperatures T N and at low magnetic fields B the polarization becomes small, making ( C el / T ) N 0.9 at the LFL region, as it is seen from Figure 2. Thus, the scaling behavior of C mag / T shown in Figure 1 is of universal character; indeed, ( C mag / T ) N = M N * of Lu 3 Cu 2 Sb 3 O 14 behaves like ( C el / T ) N = χ N = M N * shown in Figure 2 and Figure 3 The data shown are extracted from measurements on ZnCu 3 ( OH ) 6 Cl 2 and YbRh 2 Si 2 [16,52,53]. Thus, the quantum spin liquid of Lu 3 Cu 2 Sb 3 O 14 can be viewed as SCQSL that exhibits gapless behavior even in the presence of a strong magnetic field. From Figure the data at B = 0 and T N < 1 indicate that the quantum spin liquid exhibits the LFL behavior. Thus, we see that the quantum spin liquid in Lu 3 Cu 2 Sb 3 O 14 is located before the topological FCQPT. Otherwise, the spin liquid, being on the ordered side of the topological FCQPT, would have been consumed by phase transition, eliminating the corresponding finite value of the residual entropy S 0 at S ( T 0 ) S 0 [1,6]. In that case one can experimentally observe competition of the different phase transitions at T 0 that make the corresponding phase diagram very complicated, and the only reason of these complexity and competition of possible phase transitions is to vanish S 0 due to the Nernst law [5,6]. Therefore, we conclude that SCQSL without a gap should be close to the topological FCQPT, and is located on the disordered side of the topological FCQPT.
In Figure 5a, the solid squares denote the values of the maxima ( C mag / T ) max ( B ) versus magnetic field B, taken from the inset of Figure 1. Clearly the agreement between the theory (solid curve) and the experiment is good. At B = 0 the arrow shows the position of the maximum with the subtracted impurity Schottky contribution [23]. We believe that there is no reason to subtract the the impurity Schottky contribution, since it is not possible to differentiate the contribution from the impurities and from those coming from the pure crystal holding SCQSL [45], since both of them form the integral SCQSL [14]. The solid line in Figure 5b represents function T max ( B ) B , given by Equation (8). It is seen that the data are well approximated by the straight line. At B = 0 T max is finite, pointing to the fact that SCQSL exhibits the LFL behavior. Indeed, at B = 0 and T 0 the specific heat C mag / T demonstrates the LFL behavior, as evident from the inset of Figure 1. This behavior indicates that SCQSL of Lu 3 Cu 2 Sb 3 O 14 is placed before the topological FCQPT. Thus, at T 0 the system exhibits the LFL behavior and the absence of a gap or some phase transition, that has to eliminate the residual entropy S 0 . This conclusion is consistent with the general properties of the phase diagrams of HF metals and quantum insulators [14,54]. To clarify the above mentioned properties, one needs to carry out low temperature measurements of both the magnetic susceptibility χ and the thermal transport in presence of the magnetic field.
A few remarks related to the heat transport in frustrated insulators are in order here. Recent measurements of the low-temperature thermal conductivity κ have shown that the value of κ ( T 0 ) strongly depends on the disorder of quantum magnets (insulators) and at high disorder κ ( T 0 ) 0 , see, e.g., [55,56]. Measurements on the transition metal dichalcogenide 1 T TaS 2 demonstrates that it hosts QSL on the two-dimensional perfect triangular lattice [56]. Experiments show that the application of magnetic field B enhances κ / T and suppresses C m a g / T [56]. These observations agree with the predictions of the FC theory [14,44,57]. On the other hand, κ ( T 0 ) 0 could signal that QSL is not present, while the thermodynamic properties of quantum insulators with κ ( T 0 ) 0 demonstrate the typical behavior of HF metals. As a result, we propose that the thermodynamic properties of these insulators are defined by QSL [1,14,28]. Thus, we have to suggest that there are at least two types of QSL: one is represented by QSL with high resistance to the heat transport, that is κ ( T 0 ) 0 , and the other is characterized by κ ( T 0 ) being finite. In the latter case κ depends on magnetic field similarly to the magnetoresistance of HF metals [29,30]. We assume that in two dimensional systems, formed by the kagomè lattice, spinons form weakly bound states with impurities and that bound states strongly obstruct the heat transport, κ ( T ) ) 0 , but this obstacle does not influence the thermodynamic properties of SCQSL [45].
There are perspective insulators with QSL that are the Kitaev materials. They can be thought of as Mott insulators, exhibiting specific exchange interactions leading to unconventional forms of magnetism induced by QSLs [58]. There are a number of experimentally studied examples like Na 2 IrO 3 , α - Li 2 IrO 3 and α RuCl 3 where local magnetic moments are aligned in interacting hexagonal layers, see, e.g., [14,59,60,61]. Measurements of thermal conductivity κ ( B ) in magnetic fields B of the insulator α RuCl 3 have shown that κ ( B ) is finite at T 0 when the antiferromagnetic order is suppressed by magnetic field B = B c 7 T, while κ ( B ) is an increasing function at B > B c [60,61]. Such a behavior points to the fact that QSL of α RuCl 3 is located on the ordered side of the topological FCQPT, and the application of magnetic field shifts the system to the point of FCQPT, similarly to the case of YbRh 2 Si 2 [62]. Then, the elevated magnetic field B enhances κ , since κ ( B ) ( M * ( B ) ) 2 with M * ( B ) given by Equation (5) [57]. These observations are in good agreement with the behavior of κ ( B ) of SCQSL, allowing us to suggest that QSL of α RuCl 3 represents SCQSL, resembling the corresponding behavior of HF metals [1,14,30,57].
The universal scaling behavior exhibited by Lu 3 Cu 2 Sb 3 O 14 is shown on Figure 6 that displays T / B scaling of the HF metal CeCu 6 x Au x and SCQSL of the frustrated insulator herbertsmithite ZnCu 3 ( OH ) 6 Cl 2 [16,63]. This universal scaling demonstrated by the very different HF compounds, allows to conclude that HF compounds represent a new state of matter [1,7]. In contrast to ordinary quantum phase transition, the universal scaling induced by FCQPT occurs up to high temperatures T < T f , where T f 100 K, for both LFL and NFL behaviors are defined by quasiparticles (with M N * given by Equation (9)), rather than by some kind of fluctuations or Kondo lattice [1,6,7]. A few remarks are in order here. A HF compound can be placed before the topological FCQPT, exactly at FCQPT, and behind it on its ordered side, see Section 6, Figure 7. One may expect that the T / B scaling is defined by some phenomena that are not related to both the presence of FCQPT and the corresponding divergence of M * , see, e.g., [1,6,7]. On the other hand, if the HF compound in question is located before FCQPT, it exhibits the LFL behavior even without the application of magnetic field B at low T 0 . At elevated magnetic fields, as magnetic field becomes B B 0 , Equation (5) is valid and the scaling restores, see also Section 6, Equation (28). Thus, to observe both scaling and the divergency of the effective mass in measurements on HF compounds, measurements should be performed at sufficiently low temperatures and magnetic fields. For instance, the HF metal CeRu 2 Si 2 shows the NFL behavior down to lowest temperatures of 170 mK and very low magnetic fields ( B 0.02 mT) [64]. We note that interpretations of the measurements carried out in presence of magnetic field can lead to incorrect theoretical results that CeRu 2 Si 2 demonstrates the LFL behavior at low temperatures [64]. Thus, we have to conclude that a theory is an important tool that allows one to understand what is being measured. For example, a conclusion that the scaling behavior of the thermodynamic properties of a HF compounds without both QCP and the divergency of the effective mass could be caused by simple misinterpretation of the obtained experimental data, see Section 6.

2.2. Schematic Phase Diagram of L u 3 C u 2 S b 3 O 14

Now we are in a position to construct the schematic phase diagram of Lu 3 Cu 2 Sb 3 O 14 displayed in Figure 8. As seen from (9) and Figure 1, at T = 0 and B = 0 the system is located before the topological FCQPT, that is on its disordered side. Therefore, at T < T 0 the system exhibits the LFL behavior, ensuring the existence of SCQSL without gap, as discussed above. Both magnetic field B and temperature T play the role of the control parameters, shifting the system from its location at B = 0 and T = T 0 (close to the topological FCQPT) and driving it from the NFL to LFL region as shown by the vertical and horizontal arrows in Figure 8. At a fixed temperature increasing B drives the system from the NFL to the LFL region.
This behavior is seen from Figure 5b: T max increases with B increasing. At T < T max the system exhibits the LFL behavior [6]. On the contrary, at fixed B and increasing T, the system, following the vertical arrow direction, shifts from the LFL to NFL region. The inset to Figure 8 displays the behavior of the normalized effective mass M N * as a function of the normalized temperature T N T / B , as seen from Equation (9).The region T N 1 represents the crossover region between the LFL behavior with almost constant effective mass and the NFL behavior, exhibiting the T 2 / 3 dependence, see Equation (4) and inset in Figure 8. We note that in the framework of the Fermion condensation theory it is possible to explain the crossover from the NFL behavior to the LFL one in the presence of the small magnetic fields [1,6,7,14,43], while the applying the pressure does not alter the NFL behavior, provided that the system is located on the ordered side of the topological FCQPT, see, e.g., [32]. In that case the residual entropy S 0 is eliminated by some phase transition like the superconducting one that occurs in the heavy-fermion superconductor β YbAlB 4 , later representing a strange metal located away from a magnetic instability, is not accompanied by fluctuations [32,43]. Furthermore, one cannot invoke quantum fluctuations to explain the corresponding properties of the phase diagram Figrue Figure 8 and the dependencies displayed in Figure 5a,b [23]. Thus, the main features of the schematic phase diagram Figrue Figure 8 show that the thermodynamic properties of Lu 3 Cu 2 Sb 3 O 14 are similar to those of the HF metal YbRh 2 Si 2 and ZnCu 3 ( OH ) 6 Cl 2 [1,6,14]. As a result, the quantum spin liquid of Lu 3 Cu 2 Sb 3 O 14 is represented by SCQSL [45].

3. Quasi-One Dimensional Quantum Spin Liquids

The behavior of quasi-one dimensional quantum spin liquid (1DQSL) is the subject of ongoing intensive experimental research in condensed matter physics (see, e.g., [22,24,25] and references therein). Recently, searching for 1DQSL, the salient experiments were performed on the 1D Heisenberg antiferromagnet insulators Cu(C 4 H 4 N 2 )(NO 3 ) 2 (CuPzN) and YbAlO 3 in the presence of magnetic field and interpreted in terms of SCQSL and the Tomonaga-Luttinger liquid (TLL) [22,24,25]. The observed thermodynamic properties of both CuPzN and YbAlO 3 are atypical and it is expected that they do not belong to the class of HF compounds, including HF metals and quasicrystals insulators with quantum spin liquid and HF metals [6,7,12,22,28,65]. In this Section we show that, contrary to conventional wisdom, both CuPzN and YbAlO 3 can be considered as insulators belonging to HF compounds, while their thermodynamic properties are defined by weakly interacting 1DQSL formed by spinons, and are similar to those of the HF compounds.
One dimensional (1D) chain of half-odd-integer spins described by the Heisenberg model can be mapped on the fermionic system [66,67,68,69]. One of the hallmark features of geometrically frustrated insulators is the spin-charge separation; frustrated spin system is disconnected from the electron system characterized by the charge gap, and forms approximate flat band, as an electron system does in metals, e.g., formed by Na [1,14]. At T = 0 1DQSL survives up to the saturation field B s 2 J , with J being the exchange coupling constant, e.g., between Cu 2 + in the 1D chains [25]. At B = B s the QCP occurs, creating the gapped field-induced paramagnetic phase [25,67,70]. That is, at B B s both antiferromagnetic (AFM) sublattices align toward the field direction, and the magnetic field B s fully polarizes 1DQSL spins, forming flat bands [22]. Thus, in 1DQSL the topological FCQPT plays a role of QCP, at which the energy band for spinons becomes almost flat at B = B s due to the purely kinematic reasons, and the effective mass M * of spinons diverges [22]. Beyond FCQPT at B > B s 1DQSL is fully polarized, leading to the vanishing magnetic susceptibility χ ( T 0 ) . Therefore, CuPzN and YbAlO 3 can be considered as weakly interacting fermions with simplest possible spectrum ε = p 2 / ( 2 m 0 ) , where p is the momentum and m 0 is the bare mass (we use the atomic units = c = 1). In vicinity of FCQPT occurring at B = B s and T = 0 , the fermion spectrum becomes almost flat, since at FCQPT p F 0 . The spinon effective mass diverges, M * m 0 / p F , as it is seen from Figure 9 [22]. In case of weak repulsion between spinons the divergence is associated with the onset of a topological transition at finite value of p F signaling that M * ( T ) T 1 / 2 , see Section 6 and Refs. [12,22,71,72,73]. Following [74], we assume that the weakly interacting 1DQSL could be thought as QSL formed by fermionic spinons generating the Fermi sphere (line) with the Fermi finite momentum p F , and carrying spin 1/2 and no charge. This observation is supported by experimental facts collected on quasi-1D HF metal YbNi 4 P 2 [75]. These facts demonstrate that the spin-charge separation is not observed, while the thermodynamic properties of YbNi 4 P 2 resemble those of HF metals including the emergence of the LFL behavior under the application of magnetic field with the resistance ρ ( T ) T 2 [75]. Note that recently a new state of matter, quasi-Fermi liquid, has been introduced [76,77] in context of 1DQSL with the bare interaction of spinons being weak. In that case the original Tomonaga-Luttinger system can exactly be mapped on a system of free spinons, whose low-temperature behavior in magnetic fields exhibits the LFL one [74]. As a result, we shall see that the T B phase diagram of 1DQSL in both CuPzN and YbAlO 3 resembles that of HF compounds. Thus, CuPzN, YbAlO 3 and YbNi 4 P 2 represent a unique possibility to observe a new type of 1DQSL whose thermodynamic properties resemble that of HF compounds including HF metals.
As 3D AFM ordering is observed in CuPzN and YbAlO 3 [24,78] with very weak coupling J between Cu 2 + chains, the proper spin Hamiltonian is of the form
B = J i S i · S i + 1 + J < i j > S i · S j B i S i z ,
with J being the interchain exchange coupling constant and J < < J is the interchain coupling. The experimental values J 10.3 K [79] and J 0.046 K [78] allow to conclude that the criterion of smallness of the ratio is met, J / J 0.0045. The Holstein-Primakoff model of the bozonization of the spin Hamiltonian (10) up to first nonlinear terms shows that their contribution to magnetization turns to be about 5% of that stemming from noninteracting boson gas [66,68,80,81,82].
As a result, both CuPzN and YbAlO 3 are indeed represented by weakly interacting fermions. Therefore, magnetization in terms of fermion number per spin is given by
N / L = 0 D ( ε ) f ( ε μ ( B ) ) d ε .
Here L is the number of spins in 1D chain, D ( ε ) is the density of states, corresponding to free fermion spectrum ε = p 2 / ( 2 m 0 ) . The chemical potential μ ( B ) = B s B , and f ( x ) = ( e x + 1 ) 1 is the well-known Fermi distribution function [25,68,69,83,84,85]. The magnetization can be expressed as M = M s N ( M s is the saturation magnetization)
M ( B , T ) = M s 2 m 0 T π 0 d x e x 2 B s B T + 1 .
Equations (9) and (12) will be used below to calculate the differential magnetic susceptibility χ
χ ( T , B ) = M ( T , B ) B
Dimensionless normalized magnetic susceptibility χ N of 1DQSL versus dimensionless variable ( T / | B B s | ) N for magnetic field below and above the saturation field B s is displayed in Figure 10. In the fermion representation of the 1DQSL ground state energy E ( n ) , it can be viewed as the Landau functional depending on the spinon distribution function n σ ( p ) , where p is the momentum. Near the topological FCQPT point, the effective mass M * is governed by the Landau Equation (3). In that case the main role of the Landau interaction F ( p 1 , p 2 ) = δ 2 E / δ n ( p 1 ) δ n ( p 2 ) is to bring the system to the topological FCQPT point, where M * at T = 0 , and the Fermi surface alters its topology so that the effective mass acquires the temperature and the magnetic field dependences, while the proportionality of the specific heat C / T and the magnetic susceptibility χ to M * holds: C / T χ M * ( T , B ) [6,7,50,51]. This feature can be used to separate solutions of the Equation (3), corresponding to specific experimental situation, for details see Section 6. Namely, Namely, experiments on YbAlO 3 [24] and on Cu ( C 4 H 4 N 2 ) ( NO 3 ) 2 [22,25] show that near the topological FCQPT at B = B s the temperature T M at which the maximum value of χ occurs vanishes, T M 0 , see Figure 11a,b. In accordance with these observations, the magnetic susceptibility of CuPzN diverges as χ ( T ) T 1 / 2 at B = 13.55 T [25] meaning that the divergence of M * is responsible for the observed behavior, as seen from Figure 9. Again, we recognize QCP at B = B s as the topological FCQPT.
It has been shown that near the FCQPT effective mass can behave as M * ( T ) T 1 / 2 , see Section 6, while the application of B drives the system to the LFL region with M * ( B ) ( B s B ) 1 / 2 [7,12]. At finite B and T, near FCQPT, solutions of Equation (3) M * ( T , B ) can be well approximated by a simple universal interpolating function [6,7,12]. The interpolation occurs between the LFL ( M * a + b T 2 ) and NFL ( M * T 1 / 2 ) regimes and represents the universal scaling of M N * ( T N ) independent of the spatial dimension of the considered system
M N * = 1 + c 2 1 + c 1 1 + c 1 T N 2 1 + c 2 T N 5 / 2 ,
where c 1 and c 2 are fitting parameters, M N * = M / M M * and T N = T / T M are the normalized effective mass and temperature respectively. Here, [6,7,12]
M M * | B s B | 1 / 2 ,
T M | B s B | .
We remaind that M M * is the maximum value of the effective mass M * taking place at T M . It is seen from Equation (16) that the normalized temperature reads
T N T / | B s B |
As a result, we obtain that χ ( T , B ) / χ max ( T , B ) = M N * ( T N ) becomes a function of the single variable T N T / B as it is shown in Figure 4, see, e.g., [1,6,7]. Here χ max ( T , B ) is the maximum of χ ( T , B ) occurring at T M ( B ) . Below Equations (14) and (17) is used along with Equation (12) to describe the experimental facts collected on CuPzN and YbAlO 3 .

4. Experiment versus Theory

As seen from Figure 12a,b our calculations of χ N for CuPzN and YbAlO 3 represented by solid curves agree quite well with the experimental values from [24,25]. Magnetic susceptibility at B < B s exhibits the LFL behavior in magnetic fields at which T / ( B s B ) < 1 , and at rising temperatures, that is T / ( B s B ) > 1 , the NFL behavior occurs. In between of the LFL and NFL behavior, the crossover takes place with the maximum of χ N ( T / B ) .
Theory agrees similarly well with the experiment in case B > B s , as shown in Figure 13a,b. Now, since he system is fully polarized, χ N 0 at T / ( B s B ) 0 ; at temperature increasing polarization dissolves, leading to increasing χ N . Then the χ N reaches its maximum value and the growing is intercepted as soon as the NFL behavior sets in.
The comparison of the experimental results for the normalized magnetization ( M c / B ) N obtained on β YbAlB 4 and CuPzN show very good agreement between these very different compounds, as seen from Figure 14. This result is in a good agreement with the theoretical curve taken from [22,43], as seen from Figure 14 that reports the scaling behavior of the magnetization M c / B 0.5 = a + ( M M s ) / ( B s B ) 0.5 as a function of T / B = T / ( B s B ) , with a being a constant. Indeed, from Figure 14, the LFL behavior occurs at T B , the crossover at T B , and the NFL one at T B , as in the case of the HF compounds [6,7,43]. As a result, we conclude that HF metals and 1DQSL exhibit the behavior similar to that of β YbAlB 4 . In Section 6, we show that the HF metals β YbAlB 4 exhibits the same universal behavior that α YbAl 1 x Fe x B 4 and YbCo 2 Ge 4 do.

5. Phase Diagram of One Dimensional Quantum Spin Liquids

Thermodynamic properties reported in Figure 10, Figure 12 and Figure 13 allow us to construct the T B phase diagram of 1DQSL, shown in Figure 15. We see from Figure 10, that the peak temperature T M vanishes as B approaches B s . Furthermore, from Figure 9 the effective mass M * of spinons does diverge at B B s , since it takes place at FCQPT. Based on these observations we construct the T B phase diagram reported in Figure 15, demonstrating that the peak dependence T M takes place over the wide range of B, for T M ( B s B ) . Thus, we conclude that the curves T M ( B ) are straight lines, representing energy scales typical for HF metals located at their QCP [1,54,86]. Since FCQPT occurs at B = B s , the phase diagram is approximately symmetric with respect to the point B = B s , and consists of the LFL, gapped Fermi liquid, crossover and the NFL regions. Some asymmetry comes from the impediment that the LFL region may be occupied by some ordered phases marked by OP in Figure 15, as it happens for YbAlO 3 [24]. The crossover regions in Figure 15 are depicted by arrows, and are represented by the straight lines that represent the B-dependencies of temperatures of approximate LFL and NFL boundaries as well as by that of T M . It is seen that NFL state occurs at relatively high temperatures. At the same time LFL region are located at low temperatures, where the spinon effective mass M * is almost constant, characteristic to the LFL behavior. At B > B s the 1DQSL transforms into a gapped field-induced paramagnetic spin liquid, as shown in Figure 15. With temperature increasing and the fixed magnetic field B, 1DQSL transits through the crossover, and enters the NFL region. The crossover region shown by oliver circles becomes wider, as 1DQSL moves from the topological FCQPT displayed by the filled red circle. We conclude that 1DQSL exhibits the typical behavior of HF compounds [1,54] forming the corresponding T B phase diagram displayed in Figure 15.

6. Universal Scaling in Heavy Fermion Metals

To address the universal scaling within the context of the topological FCQPT, we begin with examination of the scaling of the thermodynamic functions of α YbAl 1 x Fe x B 4 . The Landau functional E ( n ) representing the ground-state energy depends on the quasiparticle momentum distribution n σ ( p ) . Near the topological FCQPT, the effective mass m * is governed by the Landau equation Equation (3). In this section, in order reserve the capital letter for magnetization M, we denote the effective mass by m * . Let us recall that the Landau functional E ( n ) and Equation (3) are exact expressions [6,49]. The Landau interaction F ( p 1 , p 2 ) = δ 2 E / δ n ( p 1 ) δ n ( p 2 ) brings the system to the FCQPT point when m * at T = 0 . At this point the topology of the Fermi surface is altered, and in contrast to the LFL theory, near this point the effective mass m * acquires strong temperature and field dependencies. However, the typical LFL theory relations
C / T χ m * ,
remain intact. Approaching the FCQPT, m * ( T = 0 , B = 0 ) and Equation (3) becomes homogeneous, i.e., m * ( T = 0 , B ) B z and m * ( T , B = 0 ) T z , with z depending on the analytical properties of F [5,6,7,12]. On the ordered side of FCQPT at T = 0 , the single-particle spectrum ε ( p ) becomes flat in some interval p i < p F < p f surrounding the Fermi surface at p F . Thus, under the influence of the topological FCQPT, the two dimensional (2D) Fermi surface transforms into 3D Fermi volume:
ε ( p ) = μ ,
where μ is the chemical potential. At FCQPT the flat interval shrinks, since p i p F p f , and ε ( p ) possesses an inflection point at p F , with ε ( p p F ) μ ( p p F ) 3 . This inflection point can also emerge in the case of a non-analytical Landau interaction F, with [43]
ε ( p ) μ ( p F p ) 2 , p < p F ε ( p ) μ ( p p F ) 2 , p > p F .
At the inflection point given by Equation (20) the effective mass diverges as m * ( T 0 ) T 1 / 2 [1,12,43]. These features of ε ( p ) can be used to specify the solutions of Equation (3) according to different experimental situations. In particular, the experimental results obtained for both β - YbAlB 4 and α YbAl 1 x Fe x B 4 show that near QCP at B 0 , the magnetization obeys M ( B ) B 1 / 2 [12,31,35,36,37,42]. This behavior corresponds to the spectrum ε ( p ) given by Equation (20) with ( p f p i ) / p F 1 . Near the FCQPT and at finite B and T, the solutions of Equation (3) determining the T and B dependencies of m * ( T , B ) can be well approximated by the universal interpolating function [1,6,7,12]. The interpolation used between the LFL regime ( m * a + b T 2 ) with
m * B 1 / 2 ,
and the NFL is given by Equation (4).
The regimes given by Equations (4) and (15) are separated by the crossover region at which m * reaches its maximum value m M * at temperature T M , as seen from Figure 4, representing the universal scaling of the dimensionless normalized effective mass m N * = m * / m M as a function of the dimensionless normalized temperature T N = T / T M . Evidently, the region T N 1 represents the crossover region between the LFL behavior with almost constant effective mass and the NFL behavior, exhibiting the m * T 1 / 2 dependence, see Equation (4). As we shall see below, the inflection point T i n f can be used to reveal the universal scaling behavior. In Figure 4, T i n f shows approximately the beginning point of the crossover region [6]. Note that both m M * and T M depend on the microscopic properties of the system in question [6], while the normalized values exhibit universal scaling given by the equation:
m N * ( T N ) = m * ( T , B ) m M * = 1 + c 2 1 + c 1 1 + c 1 T N 2 1 + c 2 T N 5 / 2 .
Here c 1 and c 2 are fitting parameters. Clearly, from Equations (4) and (15), see, e.g., [1,44],
T M T i n f B ; T N = T M / T T / T i n f T / B .
From Equations (14) and (16) m * is seen to exhibit the universal scaling [1,6,44]
m * ( T , B ) = c 3 1 B m N * ( T / B ) ,
where c 3 is a constant [6,7,12]. We describe the experimental observations on α YbAl 1 x Fe x B 4 using Equations (14) and (24). Note that the scaling occurs at temperatures T T f , with T f being the temperature at which the influence of the FCQPT becomes negligible [6,7]. Based on Equation (24), we conclude that magnetization M as described within the theory of fermion condensation does exhibit the empirical scaling behavior, given by
d M ( T , B ) d T = d χ ( y ) d T d B 1 1 B d m N * ( y ) d y d y y ,
where y = T / B . Indeed, as seen from Equation (25), B d M / d T is a function of the only variable y.
To confirm the validity of Equation (25) and to demonstrate the universal scaling, we show in Figure 16 our calculated dimensionless normalized magnetization ( B 1 / 2 d M ( T , B ) / d T ) N versus the dimensionless normalized ratio ( T / B ) N . The normalization is obtained by dividing B 1 / 2 d M ( T , B ) / d T and T / B by their maximum values ( B 1 / 2 d M ( T , B ) / d T ) M and ( T / B ) M respectively. We recall that it is the normalization that reveals the universal behavior of HF compounds, since it allows one to get rid of microscopic properties of HF compounds, thereby elucidating their universal properties [1,6].
Evidently, as it is seen from Figure 16, the calculated scaling function of the ratio ( T / B ) N , taken from [43], tracks the data of the normalized quantity ( B 1 / 2 d M ( T , B ) / d T ) N well. It also follows from Equation (25) that the calculated function ( B 1 / 2 d M ( T , B ) / d T ) N exhibits the scaling as a function of ( B / T ) N . The NFL, crossover and LFL behavior are indicated in Figure 16a–c by arrows. The theory is represented by the solid curve [43], describing very well the scaling of ( B 1 / 2 d M ( T , B ) / d T ) N for the HF metals β YbAlB 4 , α YbAl 1 x Fe x B 4 and YbCo 2 Ge 4 . It is evident from Figure 16 that our calculations, not involving fitting parameters and ad hoc functions, are in good agreement with the experimental data [32,33,34,37,38].

7. Schematic Temperature—Doping Phase Diagram

To construct the schematic phase diagram of α YbAl 1 x Fe x B 4 , we define the location of α YbAl 1 x Fe x B 4 with respect to the topological FCQPT. The locations of the HF metals β YbAlB 4 and YbCo 2 Ge 4 are defined: β YbAlB 4 is located beyond FCQPT, and YbCo 2 Ge 4 before FCQPT [43,44]. When applying magnetic field B and at sufficiently low temperatures, α YbAl 1 x Fe x B 4 is driven to the LFL state having with the resistivity [34]
ρ ( T ) = ρ 0 + A ( B ) T 2 .
Measurements of the coefficient A ( B ) provides information on the location of the corresponding HF metal with respect to the topological FCQPT. Being proportional to the quasiparticle-quasiparticle scattering cross section, A ( B ) obeys the relation A ( m * ( B ) ) 2 [6,62,87], provided that system is located at the point of FCQPT. According to Equation (15) with B s = 0 , this implies that
A ( B ) 1 B .
If the system is located before FCQPT, then at low temperatures and as B 0 , the coefficient A ( B ) acquires the LFL behavior and can be approximated by the interpolating function, see, e.g., [6,44]
A ( B ) = a 1 B 2 + a 2 ,
where a 1 and a 2 are fitting parameters. From Equation (28), as B 0 the coefficient A c o n s t , similarly to the case of LFL behavior; at elevated magnetic field B one observes the behavior given by Equation (27). Figure 17 presents the fit of A ( B ) to the data extracted from the experimental data [34].
It is seen from Figure 17, that the theoretical dependence (28) agrees very well with the experimental data leading to a conclusion that the physics underlying the field-induced re-entrance into the LFL behavior under the application of magnetic field B is the same as for the HF metals and is defined by Equations (22) and (27). It is important to note here that deviations of the theoretical curve at low values of B from that given by Equation (27) are due to the fact that α - YbAl 1 x Fe x B 4 is located before the topological FCQPT, and since now the system exhibits the LFL behavior, m * does not diverge. To confirm this conclusion, we analyze the behavior of the resistivity ρ a ( T ) at low temperatures, as shown in Figure 18a,b.
From Figure 18a it is clearly seen that for 1 < T < 5 K the resistivity demonstrates the typical NFL behavior characterized by the linear dependence ρ a ( T ) T . At low temperatures T 0 this behavior is violated. Figure 18b shows that as T 0 the violation is defined by the LFL behavior, since the T-dependence of the resistivity is given by Equation (26). This observed agreement is consistent with that derived above from Figure 17.
We now turn to Figure 19 that displays the magnetization M ( T ) as a function of T for different values of the magnetic field B. In Figure 5, the inset displays the inflection point T i n f versus magnetic field B, and shows that at B = 0 the inflection point has a finite value T i n f 3 mK.
From Equation (18) it follows that at the crossover region, starting at T i n f , magnetization behaves as M / B χ m * ( T / B ) , as seen from Figure 4. Since from Equation (23), T i n f B , the finite value of T i n f as B 0 signals that the system is located before the topological FCQPT, exhibiting the LFL behavior at B = 0 and T 0 . This observation is consistent with the LFL behavior of both A ( B 0 ) shown in Figure 17 and the resistivity shown in Figure 18b.
Based on the results of the above analysis, we can construct the schematic T x / x c phase diagram of α YbAl 1 x Fe x B 4 with x c = 0.014 depicted in Figure 7, with the doping x / x c as a control parameter. At T 0 , B = 0 and x / x c > 1 the system exhibits the LFL behavior, and, therefore, is located before the topological FCQPT, that is on its disordered side, as it is shown by the blue arrow. Thus, at T 0 the system exhibits the LFL behavior, and we expect that the scaling of α YbAl 1 x Fe x B 4 (with x c = 0.042 and x c = 0.014 ) to be violated at low temperatures. At increasing temperatures T and fixed magnetic field B the system enters the crossover region and continues into the NFL one, displaying the restoration of the scaling behavior. We suggest that a fine tuning of x can place α YbAl 1 x Fe x B 4 at FCQPT with the T / B scaling down to lowest temperatures, while the doping x = 0.042 drives the system from FCQPT, and the system still exhibits the LFL behavior at relatively high temperatures [34]. Now we employ the phase diagram Figure 8 to demonstrate the similarity between frustrated insulators and HF metals. As seen from Figure 8, both the magnetic field B and the temperature T play the role of the control parameters, shifting the system from its location close to the topological FCQPT and driving it from the NFL to LFL regions as shown by the vertical and horizontal arrows. At fixed temperatures T > T 0 increasing B drives the system from the NFL to the LFL region. On the contrary, at fixed B and growing temperatures T, the system goes along the vertical arrow from the LFL to the NFL region. The same behavior is exhibited by the HF metal YbCo 2 Ge 4 that does not demonstrate scaling down to lowest temperatures, since at x / x c it is located before the quantum critical point associated with the topological FCQPT, see Figure 7. For the same reason the effective mass does not diverge at the lowest temperatures [44]. It is worth noting that within the framework of the theory of fermion condensation it is possible to explain the crossover from the NFL behavior to LFL when the small magnetic field is applied [1,6,43]. However, application of pressure does not change the NFL behavior, see, e.g., [32]. Such a behavior observed in the heavy-fermion superconductor β YbAlB 4 , a strange metal located away from a magnetic instability, is not accompanied by fluctuations [32]. Therefore, at B = 0 , β YbAlB 4 acquires a flat band, implying the presence of a fermion condensate in a strongly degenerate state of matter that becomes susceptible to transition into a superconducting state [43]. Accordingly, the HF metal β YbAlB 4 is located behind FCQPT. Thus, the general features of the schematic phase diagram Figure 8 demonstrate that at elevated T / B the thermodynamic properties of α YbAl 1 x Fe x B 4 become close to those of the HF metals YbCo 2 Ge 4 [44] and β YbAlB 4 [43], as it is seen from Figure 16. Note that at the doping x = 0.014 and relatively high temperatures [34] and elevated magnetic field B the system enters the LFL region from the transition region, and the behavior ρ ( T ) T 2 emerges, see, e.g., [6,7]. To carefully locate the position of the HF metal α YbAl 1 x Fe x B 4 with respect to the topological FCQPT one needs to carry out measurements of the resistivity at low temperatures at which the system is placed at the LFL region, as shown by the blue arrow in Figure 7. We suggest that this procedure can allow to tune the HF metal α YbAl 1 x Fe x B 4 to the topological FCQPT point by doping x.

8. Summary

We have demonstrated that the quantum spin liquid of Lu 3 Cu 2 Sb 3 O 14 and quasi-one dimensional quantum spin liquid of both YbAlO 3 and Cu ( C 4 H 4 N 2 ) ( NO 3 ) 2 can be considered as strongly correlated Fermi systems whose thermodynamic properties are defined by SCQSL located near FCQPT. Our calculations of thermodynamic properties and the constructed phase diagrams are in a good agreement with the experimental data. Thus, the quantum spin liquid and 1DQSL are well represented by SCQSL and well described within the theory of fermion condensation [1,6,7,14,45]. We remark that the observed universal scaling can hardly be explained within theories based on some kinds of fluctuations.
We have also analyzed the thermodynamic and transport properties of the heavy-fermion metal α YbAl 1 x Fe x B 4 and explained its enigmatic scaling behavior within a topological description based on the topological FCQPT. The similarity between the three different HF metals α YbAl 1 x Fe x B 4 , YbCo 2 Ge 4 and β YbAlB 4 has been explained, and have shown that their T / B scaling can be described by the same universal function. We predict that a fine tuning of x around x c = 0.014 can place α YbAl 1 x Fe x B 4 at FCQPT with the T / B scaling down to lowest temperatures, namely T 0 0 , with the effective mass exhibiting divergent behavior m * ( T ) T 1 / 2 down to lowest temperatures. We have also demonstrated that the fermion condensation theory provides good description of the scaling for various HF compounds. Our results are in a good agreement with experimental observations and allow us to conclude that both the HF metals α YbAl 1 x Fe x B 4 , β YbAlB 4 and YbCo 2 Ge 4 and the frustrated insulators Lu 3 Cu 2 Sb 3 O 14 , Cu ( C 4 H 4 N 2 ) ( NO 3 ) 2 and YbAlO 3 with SCQSL form the new state of matter [1].

Author Contributions

V.R.S., A.Z.M. and G.S.J. designed the project and directed it with the help of S.A.A. and Y.S.L. Y.S.L. and S.A.A. performed the calculations. V.R.S., A.Z.M. and G.S.J. wrote the manuscript and all authors commented on it. The manuscript reflects the contributions of all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was partly supported by U.S. DOE, Division of Chemical Sciences, Office of Basic Energy Sciences, Office of Energy Research. We thank Podlesnyak for valuable comments and generosity for providing us with some experimental results on YbAlO 3 .

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this review paper.

References

  1. Amusya, M.Y.; Shaginyan, V.R. Strongly Correlated Fermi Systems. A New State of Matter; Springer Tracts in Modern Physics; Springer Nature: Cham, Switzerland, 2020; Volume 283. [Google Scholar]
  2. Shaginyan, V.R.; Msezane, A.Z.; Japaridze, G.S.; Stephanovich, V.A. Violation of the Time-Reversal and Particle-Hole Symmetries in Strongly Correlated Fermi Systems. Symmetry 2020, 12, 1596. [Google Scholar] [CrossRef]
  3. Khodel, V.A.; Shaginyan, V.R. Superfluidity in systems with fermion condensate. JETP Lett. 1990, 51, 553. [Google Scholar]
  4. Volovik, G.E. A new class of normal Fermi liquids. JETP Lett. 1991, 53, 222. [Google Scholar]
  5. Khodel, V.A.; Shaginyan, V.R.; Khodel, V.V. New approach in the microscopic Fermi system theory. Phys. Rep. 1994, 249, 1. [Google Scholar] [CrossRef]
  6. Shaginyan, V.R.; Amusia, M.Y.; Msezane, A.Z.; Popov, K.G. Scaling behavior of heavy fermion metals. Phys. Rep. 2010, 492, 31. [Google Scholar] [CrossRef] [Green Version]
  7. Amusia, M.Y.; Popov, K.G.; Shaginyan, V.R.; Stephanovich, V.A. Theory of Heavy–Fermion Compounds; Series in Solid-State Sciences; Springer: Cham, Switzerland; Heidelberg, Germany; New York, NY, USA; Dordrecht, The Netherlands; London, UK, 2014; Volume 182. [Google Scholar]
  8. Cao, Y.; Fatemi, V.; Fang, S.; Watanabe, K.; Taniguchi, T.; Kaxiras, E.; Jarillo-Herrero, P. Unconventional superconductivity in magic-angle graphene superlattices. Nature 2018, 556, 43. [Google Scholar] [CrossRef]
  9. Regnault, N.; Xu, Y.; Li, M.-R.; Ma, D.-S.; Jovanovic, M.; Yazdani, A.; Parkin, S.S.P.; Felser, C.; Schoop, L.M.; Ong, N.P.; et al. Catalogue of flat-band stoichiometric materials. Nature 2022, 603, 824. [Google Scholar] [CrossRef]
  10. Lu, X.; Stepanov, P.; Yang, W.; Xie, M.; Aamir, M.A.; Das, I.; Urgell, C.; Watanabe, K.; Taniguchi, T.; Zhang, G.; et al. Superconductors, orbital magnets and correlated states in magic-angle bilayer graphene. Nature 2019, 574, 653. [Google Scholar] [CrossRef]
  11. Shaginyan, V.R.; Msezane, A.Z.; Popov, K.G.; Stephanovich, V.A. Universal Behavior of Two-Dimensional 3He at Low Temperatures. Phys. Rev. Lett. 2008, 100, 096406. [Google Scholar] [CrossRef] [Green Version]
  12. Shaginyan, V.R.; Msezane, A.Z.; Popov, K.G.; Japaridze, G.S.; Khodel, V.A. Common quantum phase transition in quasicrystals and heavy-fermion metals. Phys. Rev. B 2013, 87, 245122. [Google Scholar] [CrossRef] [Green Version]
  13. 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. [Google Scholar] [CrossRef] [Green Version]
  14. Shaginyan, V.R.; Stephanovich, V.A.; Msezane, A.Z.; Japaridze, G.S.; Clark, J.W.; Amusia, M.Y.; Kirichenko, E.V. Theoretical and experimental developments in quantum spin liquid in geometrically frustrated magnets: A review. J. Mater. Sci. 2020, 55, 2257. [Google Scholar] [CrossRef]
  15. Anderson, P.W. Resonating valence bonds: A new kind of insulator? Mater. Res. Bull. 1973, 8, 153. [Google Scholar] [CrossRef]
  16. 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 kagomè lattice antiferromagnet herbertsmithite. Phys. Rev. Lett. 2010, 104, 147201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. 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 kagomè-lattice antiferromagnets ZnxCu4-x(OH)6Cl2. Phys. Rev. B 2011, 83, 100402. [Google Scholar] [CrossRef] [Green Version]
  18. Balents, L. Spin liquids in frustrated magnets. Nature 2010, 464, 199. [Google Scholar] [CrossRef]
  19. 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] [Green Version]
  20. Mendels, P.; Bert, F. Quantum kagomè antiferromagnet: ZnCu3(OH)6Cl2. J. Phys. Conf. Ser. 2011, 320, 012004. [Google Scholar] [CrossRef]
  21. 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 kagomè-lattice antiferromagnet. Nature 2012, 492, 406. [Google Scholar] [CrossRef] [Green Version]
  22. 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)(NO 3)2. Ann. Phys. 2016, 528, 483. [Google Scholar] [CrossRef]
  23. Yang, Y.; Li, X.; Tan, C.; Zhu, Z.H.; Zhang, J.; Ding, Z.F.; Wu, Q.; Chen, C.S.; Shiroka, T.; MacLaughlin, D.E.; et al. Discovery of an ultra-quantum spin liquid. arXiv 2021, arXiv:2102.09271. [Google Scholar]
  24. Wu, L.S.; Nikitin, S.E.; Wang, Z.; Zhu, W.; Batista, C.D.; Tsvelik, A.M.; Samarakoon, A.M.; Tennant, D.A.; Brando, M.; Vasylechkoll, L.; et al. Tomonaga-Luttinger liquid behavior and spinon confinement in YbAlO3. Nat. Commun. 2019, 10, 698. [Google Scholar] [CrossRef] [PubMed]
  25. Kono, Y.; Sakakibara, T.; Aoyama, C.P.; Hotta, C.; Turnbull, M.M.; Landee, C.P.; Takano, Y. Field-induced quantum criticality and universal temperature dependence of the magnetization of a spin-1/2 Heisenberg Chain. Phys. Rev. Lett. 2015, 114, 037202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Khodel, V.A.; Clark, J.W.; Zverev, M.V. Topological disorder triggered by interaction-induced flattening of electron spectra in solids. Phys. Rev. B 2020, 102, 201108. [Google Scholar] [CrossRef]
  27. Nozières, P. Properties of Fermi liquids with a finite range interaction. J. Phys. 1992, 2, 443. [Google Scholar] [CrossRef]
  28. Shaginyan, V.R.; Msezane, A.Z.; Popov, K.G. Thermodynamic properties of the kagomè lattice in herbertsmithite. Phys. Rev. B 2011, 84, 060401. [Google Scholar] [CrossRef] [Green Version]
  29. Shaginyan, V.R.; Popov, K.G.; Stephanovich, V.A.; Fomichev, V.I.; Ki-richenko, E.V. High-magnetic-fields thermodynamics of the heavy-fermion metal YbRh2Si2. Europhys. Lett. 2011, 93, 17008. [Google Scholar] [CrossRef] [Green Version]
  30. 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] [Green Version]
  31. Matsumoto, Y.; Nakatsuji, S.; Kuga, K.; Karaki, Y.; Horie, N.; Shimura, Y.; Sakakibara, T.; Nevidomskyy, A.H.; Coleman, P. Quantum criticality without tuning in the mixed valence compound β-YbAlB4. Science 2011, 331, 316. [Google Scholar] [CrossRef] [Green Version]
  32. Tomita, T.; Kuga, K.; Uwatoko, Y.; Coleman, P.; Nakatsuji, S. Strange metal without magnetic criticality. Science 2015, 349, 506. [Google Scholar] [CrossRef] [Green Version]
  33. Sakai, A.; Kitagawa, K.; Matsubayashi, K.; Iwatani, M.; Gegenwart, P. T/B scaling without quasiparticle mass divergence. Phys. Rev. B 2016, 94, 041106. [Google Scholar] [CrossRef] [Green Version]
  34. Kuga, K.; Matsumoto, Y.; Okawa, M.; Suzuki, S.; Tomita, T.; Sone, K.; Shimura, Y.; Sakakibara, T.; Nishio-Hamane, D.; Karaki, Y.; et al. Quantum valence criticality in a correlated metal. Sci. Adv. 2018, 4, eaao3547. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Nakatsuji, S.; Kuga, K.; Machida, Y.; Tayama, T.; Sakakibara, T.; Karaki, Y.; Ishimoto, H.; Yonezawa, S.; Maeno, Y.; Pearson, E.; et al. Superconductivity and quantum criticality in the heavy-fermion system β-YbAlB4. Nat. Phys. 2008, 4, 603. [Google Scholar] [CrossRef]
  36. Matsumoto, Y.; Kuga, K.; Tomita, T.; Karaki, Y.; Nakatsuji, S. Anisotropic heavy-Fermi-liquid formation in valence-fluctuating α-YbAlB4. Phys. Rev. B 2011, 84, 125126. [Google Scholar] [CrossRef] [Green Version]
  37. Matsumoto, Y.; Nakatsuji, S.; Kuga, K.; Karaki, Y.; Shimura, Y.; Sakakibara, T.; Nevidomskyy, A.H.; Coleman, P. T/B scaling of magnetization in the mixed valent compound β-YbAlB4. J. Phys. Conf. Ser. 2012, 391, 012041. [Google Scholar] [CrossRef] [Green Version]
  38. Matsumoto, Y.; Kuga, K.; Karaki, Y.; Shimura, Y.; Sakakibara, T.; Tokunaga, M.; Kindo, K.; Nakatsuji, S. Field Evolution of Quantum Critical and Heavy Fermi-Liquid Components in the Magnetization of the Mixed Valence Compound β-YbAlB4. J. Phys. Soc. Jpn. 2015, 84, 024710. [Google Scholar] [CrossRef] [Green Version]
  39. Nevidomskyy, A.H.; Coleman, P. Layered Kondo Lattice Model for Quantum Critical β-YbAlB4. Phys. Rev. Lett. 2009, 102, 077202. [Google Scholar] [CrossRef] [Green Version]
  40. Ramires, A.; Coleman, P.; Nevidomskyy, A.H.; Tsvelik, A.M. Evolution of c-f Hybridization and Two-Component Hall Effect in β-YbAlB4. Phys. Rev. Lett. 2012, 109, 176404. [Google Scholar] [CrossRef] [Green Version]
  41. Watanabe, S.; Miyake, K. Robustness of Quantum Criticality of Valence Fluctuations. J. Phys. Soc. Jpn. 2013, 82, 083704. [Google Scholar] [CrossRef] [Green Version]
  42. Watanabe, S.; Miyake, K. T/B Scaling in β-YbAlB4. J. Phys. Soc. Jpn. 2014, 83, 103708. [Google Scholar] [CrossRef] [Green Version]
  43. Shaginyan, V.R.; Msezane, A.Z.; Popov, K.G.; Clark, J.W.; Khodel, V.A.; Zverev, M.V. Topological basis for understanding the behavior of the heavy-fermion metal β-YbAlB4 under application of magnetic field and pressure. Phys. Rev. B 2016, 93, 205126. [Google Scholar] [CrossRef] [Green Version]
  44. Shaginyan, V.R.; Msezane, A.Z.; Clark, J.W.; Japaridze, G.S.; Leevik, Y.S. Universal T/B scaling behavior of heavy fermion compounds. JETP Lett. 2020, 112, 657. [Google Scholar] [CrossRef]
  45. Shaginyan, V.R.; Msezane, A.Z.; Artamonov, S.A.; Japaridze, G.S.; Leevik, Y.S. Universal Ultra spin liquid in Lu3Cu2Sb3O14. Europhys. Lett. 2021, 136, 17003. [Google Scholar] [CrossRef]
  46. 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] [Green Version]
  47. Heikkilá, T.T.; Kopnin, N.B.; Volovik, G.E. Flat bands in topological media. JETP Lett. 2011, 94, 252. [Google Scholar] [CrossRef] [Green Version]
  48. Landau, L.D. The theory of a Fermi liquid. Sov. Phys. JETP 1956, 3, 920. [Google Scholar]
  49. Shaginyan, V.R. Density functional theory of fermion condensation. Phys. Lett. A 1998, 249, 237. [Google Scholar] [CrossRef]
  50. Clark, J.W.; Khodel, V.A.; Zverev, M.V. Anomalous low-temperature behavior of strongly correlated Fermi systems. Phys. Rev. B 2005, 71, 012401. [Google Scholar] [CrossRef] [Green Version]
  51. Khodel, V.A.; Clark, J.W.; Zverev, M.V. Topology of the Fermi surface beyond the quantum critical point. Phys. Rev. B 2008, 78, 075120. [Google Scholar] [CrossRef] [Green Version]
  52. Gegenwart, P.; Tokiwa, Y.; Westerkamp, T.; Weickert, F.; Custers, J.; Ferstl, J.; Krellner, C.; Geibel, C.; Kerschl, P.; Müller, K.H.; et al. High-field phase diagram of the heavy-fermion metal YbRh2Si2. New J. Phys. 2006, 8, 171. [Google Scholar] [CrossRef] [Green Version]
  53. Oeschler, N.; Hartmann, S.; Pikul, A.; Krellner, C.; Geibel, C.; Steglich, F. Low-temperature specific heat of YbRh2Si2. Phys. B 2008, 403, 1254. [Google Scholar] [CrossRef]
  54. Shaginyan, V.R.; Msezane, A.Z.; Popov, K.G.; Japaridze, G.S.; Khodel, V.A. General properties of phase diagrams of heavy-fermion metals. Europhys. Lett. 2014, 106, 37001. [Google Scholar] [CrossRef] [Green Version]
  55. Yamashita, M.; Sato, Y.; Tominaga, T.; Kasahara, Y.; Kasahara, S.; Cui, H.; Kato, R.; Shibauchi, T.; Matsuda, Y. Presence and absence of itinerant gapless excitations in the quantum spin liquid candidate EtMe3Sb[Pd(dmit)2]2. Phys. Rev. B 2020, 101, 140407. [Google Scholar] [CrossRef] [Green Version]
  56. Murayama, H.; Sato, Y.; Taniguchi, T.; Kurihara, R.; Xing, X.Z.; Huang, W.; Kasahara, S.; Kasahara, Y.; Kimchi, I.; Yoshida, M.; et al. Effect of quenched disorder on the quantum spin liquid state of the triangular-lattice antiferromagnet 1T-TaS2. Phys. Rev. Res. 2020, 2, 013099. [Google Scholar] [CrossRef] [Green Version]
  57. 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 κ-(BEDT-TTF)2Cu2(CN)3. Europhys. Lett. 2013, 103, 67006. [Google Scholar] [CrossRef] [Green Version]
  58. Kitaev, A. Anyons in an exactly solved model and beyond. Ann. Phys. 2006, 321, 2. [Google Scholar] [CrossRef] [Green Version]
  59. Savary, L.; Balents, L. Quantum spin liquids: A review. Rep. Prog. Phys. 2017, 80, 016502. [Google Scholar] [CrossRef]
  60. Kasahara, Y.; Ohnishi, T.; Mizukami, Y.; Tanaka, O.; Ma, S.; Sugii, K.; Kurita, N.; Tanaka, H.; Nasu, J.; Motome, Y.; et al. Majorana quantization and half-integer thermal quantum Hall effect in a Kitaev spin liquid. Nature 2018, 559, 227. [Google Scholar] [CrossRef]
  61. Czajka, P.; Gao, T.; Hirschberger, M.; Lampen-Kelley, P.; Banerjee, A.; Yan, J.; Mandrus, D.G.; Nagler, S.E.; Ong, N.P. Oscillations of the thermal conductivity in the spin-liquid state of α-RuCl3. Nat. Phys. 2021, 17, 915. [Google Scholar] [CrossRef]
  62. Gegenwart, P.; Custers, J.; Geibel, C.; Neumaier, K.; Tayama, T.; Tenya, K.; Trovarelli, O.; Steglich, F. Magnetic-Field Induced Quantum Critical Point in YbRh2Si2. Phys. Rev. Lett. 2002, 89, 056402. [Google Scholar] [CrossRef] [Green Version]
  63. Schroëder, A.; Aeppli, G.; Coldea, R.; Adams, M.; Stockert, O.; Loèhneysen, H.V.; Bucher, E.; Ramazashvili, R.; Coleman, P. Onset of antiferromagnetism in heavy-fermion metals. Nature 2000, 407, 351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Takahashi, D.; Abe, S.; Mizuno, H.; Tayurskii, D.; Matsumoto, K.; Suzuki, H.; Onuki, Y. ac susceptibility and static magnetization measurements of CeRu2Si2 at small magnetic fields and ultralow temperatures. Phys. Rev. B 2003, 67, 180407. [Google Scholar] [CrossRef] [Green Version]
  65. Shechtman, D.; Blech, I.; Gratias, D.; Cahn, J.W. Metallic phase with long-range orientational order and no translational symmetry. Phys. Rev. Lett. 1984, 53, 1951. [Google Scholar] [CrossRef] [Green Version]
  66. Mattis, D.C.; Lieb, E.H. Exact solution of a many-fermion system and its associated boson field. J. Math. Phys. 1965, 6, 304. [Google Scholar] [CrossRef] [Green Version]
  67. Haldane, F.D.M. General relation of correlation exponents and spectral properties of one-dimensional fermi systems: Application to the anisotropic S =1/2 Heisenberg chain. Phys. Rev. Lett. 1980, 45, 1358. [Google Scholar] [CrossRef]
  68. Haldane, F.D.M. Effective harmonic-fluid approach to low-energy properties of one-dimensional quantum fluids. Phys. Rev. Lett. 1981, 47, 1840. [Google Scholar] [CrossRef]
  69. Affleck, I. Bose condensation in quasi-one-dimensional antiferromagnets in strong fields. Phys. Rev. B 1991, 43, 3215. [Google Scholar] [CrossRef]
  70. Müller, G.; Thomas, H.; Beck, H.; Bonner, J.C. Quantum spin dynamics of the antiferromagnetic linear chain in zero and nonzero magnetic field. Phys. Rev. B 1981, 24, 1429. [Google Scholar] [CrossRef] [Green Version]
  71. Lifshitz, I.M. Anomalies of Electron Characteristics of a Metal in the High Pressure Region. Sov. Phys. JETP 1960, 11, 1130. [Google Scholar]
  72. Volovik, G.E. From Standard Model of particle physics to room-temperature superconductivity. Phys. Scr. T 2015, 164, 014014. [Google Scholar] [CrossRef] [Green Version]
  73. Khodel, V.A.; Clark, J.W.; Popov, K.G.; Shaginyan, V.R. Occurrence of Flat Bands in Strongly Correlated Fermi Systems and High-Tc Superconductivity of Electron-Doped Compounds. JETP Lett. 2015, 101, 413. [Google Scholar] [CrossRef] [Green Version]
  74. Rozhkov, A.V. Fermionic quasiparticle representation of Tomonaga-Luttinger Hamiltonian. Eur. Phys. J. B 2005, 47, 193. [Google Scholar] [CrossRef]
  75. Krellner, C.; Lausberg, S.; Steppke, A.; Brando, M.; Pedrero, L.; Pfau, H.; Tencé, S.; Rosner, H.; Steglich, F.; Geibel, C. Ferromagnetic quantum criticality in the quasi-one-dimensional heavy fermion metal YbNi4P2. New J. Phys. 2011, 13, 103014. [Google Scholar] [CrossRef]
  76. Rozhkov, A.V. One-dimensional fermions with neither Luttinger-liquid nor Fermi-liquid behavior. Phys. Rev. Lett. 2014, 112, 106403. [Google Scholar] [CrossRef] [Green Version]
  77. Lebed, A.G. Non-Fermi-liquid crossovers in a quasione- dimensional conductor in a tilted magnetic field. Phys. Rev. Lett. 2015, 115, 157001. [Google Scholar] [CrossRef] [Green Version]
  78. Lancaster, T.; Blundell, S.J.; Brooks, M.L.; Baker, P.J.; Pratt, F.L.; Manson, J.L.; Landee, C.P.; Baines, C. Magnetic order in the quasi-one-dimensional spin-1/2 molecular chain compound copper pyrazine dinitrate. Phys. Rev. B 2006, 73, 020410. [Google Scholar] [CrossRef] [Green Version]
  79. Hammar, P.R.; Stone, M.B.; Reich, D.H.; Broholm, C.; Gibson, P.J.; Turnbull, M.M.; Landee, C.P.; Oshikawa, M. Characterization of a quasi-one-dimensional spin-1/2 magnet which is gapless and paramagnetic for gμBHJ and kBT << J. Phys. Rev. B 1999, 59, 1008. [Google Scholar]
  80. Holstein, T.; Primakoff, H. Field Dependence of the Intrinsic Domain Magnetization of a Ferromagnet. Phys. Rev. 1940, 58, 1098. [Google Scholar] [CrossRef]
  81. Stephanovich, V.A.; Zhitomirsky, M.E. Role of dimensionality in spontaneous magnon decay: Easy-plane ferromagnet. Phys. Rev. B 2014, 89, 224415. [Google Scholar] [CrossRef] [Green Version]
  82. Gutfreund, H.; Schick, M. Momentum Distribution in the Tomonaga Model. Phys. Rev. 1968, 168, 418. [Google Scholar] [CrossRef]
  83. Maeda, Y.; Hotta, C.; Oshikawa, M. Universal temperature dependence of the magnetization of gapped spin chains. Phys. Rev. Lett. 2007, 99, 057205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Sachdev, S.; Senthil, T.; Shankar, R. Finite-temperature properties of quantum antiferromagnets in a uniform magnetic field in one and two dimensions. Phys. Rev. B 1994, 50, 258. [Google Scholar] [CrossRef] [Green Version]
  85. Nikuni, T.; Oshikawa, M.; Oosawa, A.; Tanaka, H. Bose-Einstein condensation of dilute magnons in TlCuCl3. Phys. Rev. Lett. 2000, 84, 5868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Klingner, C.; Krellner, C.; Brando, M.; Geibel, C.; Steglich, F. Magnetic behaviour of the intermetallic compound YbCo2Si2. New J. Phys. 2011, 13, 83024. [Google Scholar] [CrossRef]
  87. Khodel, V.A.; Schuck, P.Z. Universal behavior of the collision rate in strongly correlated Fermi systems. Z. Phys. B 1997, 104, 505. [Google Scholar] [CrossRef]
Figure 2. The normalized specific heat ( C el / T ) N = M N * of YbRh 2 Si 2 as a function of normalized temperature T N under the application of magnetic field B shown in the left hand legend (low B) and the right hand legend (high B) The experimental data is extracted from the measurement of C / T measurements on the archetypical HF metal YbRh 2 Si 2 [52,53]. The low-field calculations of M N * are depicted by the solid green curve. The solid blue curve representing high-field calculations ( B 18 ) is performed for the fully polarized quasiparticle band [29]. Starting from relatively high magnetic fields B 4 T, the specific heat demonstrates the same behavior as ( C mag / T ) = M N * of Lu 3 Cu 2 Sb 3 O 14 shown in Figure 1.
Figure 2. The normalized specific heat ( C el / T ) N = M N * of YbRh 2 Si 2 as a function of normalized temperature T N under the application of magnetic field B shown in the left hand legend (low B) and the right hand legend (high B) The experimental data is extracted from the measurement of C / T measurements on the archetypical HF metal YbRh 2 Si 2 [52,53]. The low-field calculations of M N * are depicted by the solid green curve. The solid blue curve representing high-field calculations ( B 18 ) is performed for the fully polarized quasiparticle band [29]. Starting from relatively high magnetic fields B 4 T, the specific heat demonstrates the same behavior as ( C mag / T ) = M N * of Lu 3 Cu 2 Sb 3 O 14 shown in Figure 1.
Materials 15 03901 g002
Figure 3. Normalized magnetic susceptibility χ N = χ / χ max = M N * as a function of the normalized temperature T N T / B . The values of magnetic field B are displayed in the legend. The data are extracted from the measurements of the magnetic susceptibility χ ( T , B ) on ZnCu 3 ( OH ) 6 Cl 2 [16]. The normalized data of ( C el / T ) N = M N * are obtained from the specific heat measurements C el / T of YbRh 2 Si 2 in the presence of magnetic field B (the legend) [52]. The solid curve represents theoretical calculations at B 18 T at the fully polarized quasiparticle band. It demonstrates universal scaling of M N * , and coincides with the universal scaling of the quantum spin liquid depicted in Figure 1. The crossover from the LFL behavior to the NFL one is displayed by the arrow.
Figure 3. Normalized magnetic susceptibility χ N = χ / χ max = M N * as a function of the normalized temperature T N T / B . The values of magnetic field B are displayed in the legend. The data are extracted from the measurements of the magnetic susceptibility χ ( T , B ) on ZnCu 3 ( OH ) 6 Cl 2 [16]. The normalized data of ( C el / T ) N = M N * are obtained from the specific heat measurements C el / T of YbRh 2 Si 2 in the presence of magnetic field B (the legend) [52]. The solid curve represents theoretical calculations at B 18 T at the fully polarized quasiparticle band. It demonstrates universal scaling of M N * , and coincides with the universal scaling of the quantum spin liquid depicted in Figure 1. The crossover from the LFL behavior to the NFL one is displayed by the arrow.
Materials 15 03901 g003
Figure 4. Schematic plot of the normalized effective mass M N * = M * / M M * versus the normalized temperature T N = T / T M . The crossover region, at which M N * reaches its maximum M N * = 1 at T N = 1 , is shown by the arrows. The inflection point T i n f at which the system enters the crossover region is displayed by the arrow. The LFL and NFL regions are labeled.
Figure 4. Schematic plot of the normalized effective mass M N * = M * / M M * versus the normalized temperature T N = T / T M . The crossover region, at which M N * reaches its maximum M N * = 1 at T N = 1 , is shown by the arrows. The inflection point T i n f at which the system enters the crossover region is displayed by the arrow. The LFL and NFL regions are labeled.
Materials 15 03901 g004
Figure 5. The properties of the specific heat. (a) The maximum values of ( C mag / T ) max of the specific heat C mag / T versus magnetic field B are shown by the solid squares, see the inset, Figure 1. The solid curve is approximated by M max * ( B ) B 2 / 3 in accordance with Equation (5). The arrow depicts the position of ( C mag / T ) max at B = 0 when the impurity Schottky contribution is subtracted [23]. (b) The temperature T max ( B ) , at which the maxima of ( C mag / T ) are located, see Figure 1. The solid straight line traces the function T max B , see Equation (7).
Figure 5. The properties of the specific heat. (a) The maximum values of ( C mag / T ) max of the specific heat C mag / T versus magnetic field B are shown by the solid squares, see the inset, Figure 1. The solid curve is approximated by M max * ( B ) B 2 / 3 in accordance with Equation (5). The arrow depicts the position of ( C mag / T ) max at B = 0 when the impurity Schottky contribution is subtracted [23]. (b) The temperature T max ( B ) , at which the maxima of ( C mag / T ) are located, see Figure 1. The solid straight line traces the function T max B , see Equation (7).
Materials 15 03901 g005
Figure 6. The universal scaling behavior of different strongly correlated Fermi systems versus B / T . The universal behavior of the HF metal CeCu 6 x Au x is extracted from data [63], and that of ZnCu 3 ( OH ) 6 Cl 2 is derived from data [16]. At B / T 1 the systems exhibit the NFL behavior, that is T 2 / 3 χ c o n s t . At B / T 1 the systems demonstrate the LFL behavior, with χ being a decreasing function of B / T .
Figure 6. The universal scaling behavior of different strongly correlated Fermi systems versus B / T . The universal behavior of the HF metal CeCu 6 x Au x is extracted from data [63], and that of ZnCu 3 ( OH ) 6 Cl 2 is derived from data [16]. At B / T 1 the systems exhibit the NFL behavior, that is T 2 / 3 χ c o n s t . At B / T 1 the systems demonstrate the LFL behavior, with χ being a decreasing function of B / T .
Materials 15 03901 g006
Figure 7. Schematic temperature T—doping x / x c phase diagrams of HF metals. The number density x is taken as the control parameter and depicted as x / x c . At x / x c < 1 the dashed arrow shows the the ordered phase of the topological FCQPT, when the system possesses flat bands. At any finite temperature T > 0 and at x / x c < 1 , the system exhibits the NFL. The shadowed area corresponds to the case x / x c > 1 and sufficiently low temperatures, where the system is in the LFL phase.
Figure 7. Schematic temperature T—doping x / x c phase diagrams of HF metals. The number density x is taken as the control parameter and depicted as x / x c . At x / x c < 1 the dashed arrow shows the the ordered phase of the topological FCQPT, when the system possesses flat bands. At any finite temperature T > 0 and at x / x c < 1 , the system exhibits the NFL. The shadowed area corresponds to the case x / x c > 1 and sufficiently low temperatures, where the system is in the LFL phase.
Materials 15 03901 g007
Figure 8. Schematic T B phase diagram of Lu 3 Cu 2 Sb 3 O 14 with magnetic field as the control parameter. The vertical and horizontal arrows depict LFL-NFL transitions at fixed B and T, respectively. The line separating LFL-NFL regions is shown by the arrow, representing the crossover region at T max ( B ) and displaying the function T max ( B ) . T 0 is the temperature at which the LFL behavior occurs. The inset demonstrates a chart of the normalized effective mass M N * as a function of T N T / B . Transition region is characterized by the maximum value M max * of M * at T N = T / T max = 1 .
Figure 8. Schematic T B phase diagram of Lu 3 Cu 2 Sb 3 O 14 with magnetic field as the control parameter. The vertical and horizontal arrows depict LFL-NFL transitions at fixed B and T, respectively. The line separating LFL-NFL regions is shown by the arrow, representing the crossover region at T max ( B ) and displaying the function T max ( B ) . T 0 is the temperature at which the LFL behavior occurs. The inset demonstrates a chart of the normalized effective mass M N * as a function of T N T / B . Transition region is characterized by the maximum value M max * of M * at T N = T / T max = 1 .
Materials 15 03901 g008
Figure 9. Temperature dependence of χ ( T ) for CuPzN. The experimental data are from [25]. The magnetic fields are shown in the legend. The three regions LFL, NFL and crossover are shown. The divergent behavior of χ ( T ) at B = B s is indicated by the black arrow.
Figure 9. Temperature dependence of χ ( T ) for CuPzN. The experimental data are from [25]. The magnetic fields are shown in the legend. The three regions LFL, NFL and crossover are shown. The divergent behavior of χ ( T ) at B = B s is indicated by the black arrow.
Materials 15 03901 g009
Figure 10. The normalized magnetic susceptibility χ N vs. the normalized ( T / | B B s | ) N , calculated from Equations (12) and (13) for B < B s (panel (a)) and B > B s (panel (b)). At B > B s the complete spin polarization occurs and χ vanish at ( T / B ) N 0 . Two different curves on panel (a) show an excellent scaling of χ N ( T / B ) N . As seen from panel (a), the LFL behavior holds for the weakly interacting quasi-1D quantum spin liquid [22].
Figure 10. The normalized magnetic susceptibility χ N vs. the normalized ( T / | B B s | ) N , calculated from Equations (12) and (13) for B < B s (panel (a)) and B > B s (panel (b)). At B > B s the complete spin polarization occurs and χ vanish at ( T / B ) N 0 . Two different curves on panel (a) show an excellent scaling of χ N ( T / B ) N . As seen from panel (a), the LFL behavior holds for the weakly interacting quasi-1D quantum spin liquid [22].
Materials 15 03901 g010
Figure 11. The magnetic field dependence of peak temperature T M ( χ ) . Panel (a) shows the B T M dependence of CuPzN, extracted from data [25]. Panel (b) demonstrate the B T M dependence for YbAlO 3 , extracted from experimental data [24]. The calculated straight lines T M = a | B s B | given by Equation (17). The excellent coincidence is seen, showing that Equation (17) holds for both CuPzN and YbAlO 3 .
Figure 11. The magnetic field dependence of peak temperature T M ( χ ) . Panel (a) shows the B T M dependence of CuPzN, extracted from data [25]. Panel (b) demonstrate the B T M dependence for YbAlO 3 , extracted from experimental data [24]. The calculated straight lines T M = a | B s B | given by Equation (17). The excellent coincidence is seen, showing that Equation (17) holds for both CuPzN and YbAlO 3 .
Materials 15 03901 g011
Figure 12. The normalized magnetic susceptibility χ N extracted from measurements at B < B s . Panel (a): χ N extracted from measurements on CuPzN [25]. Panel (b): χ N extracted from measurements on YbAlO 3 [24]. Our theoretical curve, plotted based on Equations (12) and (14), is represented by the solid line (shown in Figure 10a) tracing the scaling. It is seen that the dependence χ N for YbAlO 3 has three regions: LFL, crossover and NFL.
Figure 12. The normalized magnetic susceptibility χ N extracted from measurements at B < B s . Panel (a): χ N extracted from measurements on CuPzN [25]. Panel (b): χ N extracted from measurements on YbAlO 3 [24]. Our theoretical curve, plotted based on Equations (12) and (14), is represented by the solid line (shown in Figure 10a) tracing the scaling. It is seen that the dependence χ N for YbAlO 3 has three regions: LFL, crossover and NFL.
Materials 15 03901 g012
Figure 13. The normalized magnetic susceptibility χ N extracted from measurements in magnetic fields B > B s shown in the legend. Panel (a): χ N extracted from measurements on CuPzN [25]. Panel (b): χ N extracted from measurements on YbAlO 3 [24]. Our theoretical curves, plotted on the base of Equations (12) and (14), is taken from Figure 10b, and is reported by the solid line tracing the scaling behavior.
Figure 13. The normalized magnetic susceptibility χ N extracted from measurements in magnetic fields B > B s shown in the legend. Panel (a): χ N extracted from measurements on CuPzN [25]. Panel (b): χ N extracted from measurements on YbAlO 3 [24]. Our theoretical curves, plotted on the base of Equations (12) and (14), is taken from Figure 10b, and is reported by the solid line tracing the scaling behavior.
Materials 15 03901 g013
Figure 14. The scaling dependence of magnetization normalized to maximal values ( M c / B 1 / 2 ) N ( M c = M s M , B = B s B ) on ( T / B ) N for CuPzN and β YbAlB 4 . The experimental data are taken from [25,31]. The magnetic fields B (in T) are shown in the legends. The typical LFL and NFL behavior are shown by both the arrows and the straight lines, the crossover is displayed by the arrow. The theory is represented by the solid green curve [22,43].
Figure 14. The scaling dependence of magnetization normalized to maximal values ( M c / B 1 / 2 ) N ( M c = M s M , B = B s B ) on ( T / B ) N for CuPzN and β YbAlB 4 . The experimental data are taken from [25,31]. The magnetic fields B (in T) are shown in the legends. The typical LFL and NFL behavior are shown by both the arrows and the straight lines, the crossover is displayed by the arrow. The theory is represented by the solid green curve [22,43].
Materials 15 03901 g014
Figure 15. Schematic magnetic field—temperature phase diagram of 1DQSL. Straight lines on both sides of B s , which is a FCQPT point, indicate, respectively, the lines of LFL boundary (the lowest temperature), the temperatures of maxima (middle line, marked “ T M ” taken from Figure 11a) and the end of crossover region: The highest temperature at which the system enters the NFL regime. The left sector labeled as “LFL” and “OP” (ordered phase) displays the LFL behavior and possible ordered phase of spin liquid. The right sector labeled as “gapped FL” denotes the gapped field-induced ferromagnetic spin liquid.
Figure 15. Schematic magnetic field—temperature phase diagram of 1DQSL. Straight lines on both sides of B s , which is a FCQPT point, indicate, respectively, the lines of LFL boundary (the lowest temperature), the temperatures of maxima (middle line, marked “ T M ” taken from Figure 11a) and the end of crossover region: The highest temperature at which the system enters the NFL regime. The left sector labeled as “LFL” and “OP” (ordered phase) displays the LFL behavior and possible ordered phase of spin liquid. The right sector labeled as “gapped FL” denotes the gapped field-induced ferromagnetic spin liquid.
Materials 15 03901 g015
Figure 16. Scaling of dimensionless normalized magnetization ( B 1 / 2 d M ( T , B ) / d T ) N as a function of the dimensionless normalized ( T / B ) N at magnetic field values B given in the legends (a,b). Regions of LFL behavior, crossover, and NFL behavior are indicated by arrows. The theory is represented by the solid curve [14,43]. Data are extracted from [32,33,34]. Panel (a): Scaling of β YbAlB 4 [43]. Panel (b): Scaling of α YbAl 1 x Fe x B 4 . Panel (c): Scaling of YbCo 2 Ge 4 , measured at different field values B = 0.05 , 0.1 , 0.2 , 0.3 , 0.5 T [33].
Figure 16. Scaling of dimensionless normalized magnetization ( B 1 / 2 d M ( T , B ) / d T ) N as a function of the dimensionless normalized ( T / B ) N at magnetic field values B given in the legends (a,b). Regions of LFL behavior, crossover, and NFL behavior are indicated by arrows. The theory is represented by the solid curve [14,43]. Data are extracted from [32,33,34]. Panel (a): Scaling of β YbAlB 4 [43]. Panel (b): Scaling of α YbAl 1 x Fe x B 4 . Panel (c): Scaling of YbCo 2 Ge 4 , measured at different field values B = 0.05 , 0.1 , 0.2 , 0.3 , 0.5 T [33].
Materials 15 03901 g016
Figure 17. Experimental data for the coefficient A ( B ) , plotted as a function of magnetic field B (solid squares). Measured values of A ( B ) are extracted from the experimental data [34]. The solid curve is given by Equation (28).
Figure 17. Experimental data for the coefficient A ( B ) , plotted as a function of magnetic field B (solid squares). Measured values of A ( B ) are extracted from the experimental data [34]. The solid curve is given by Equation (28).
Materials 15 03901 g017
Figure 18. Experimental data for the resistivity ρ a ( T ) : in panel (a) resistivity is plotted as a function of T, solid circles and in panel (b) resistivity is plotted as a function of T 2 , solid squares. The data are taken from [34]. The solid lines are fits to the experimental data.
Figure 18. Experimental data for the resistivity ρ a ( T ) : in panel (a) resistivity is plotted as a function of T, solid circles and in panel (b) resistivity is plotted as a function of T 2 , solid squares. The data are taken from [34]. The solid lines are fits to the experimental data.
Materials 15 03901 g018
Figure 19. Magnetization M ( T ) versus T for B values shown in the legend. The data are taken from [34]. M ( T ) is displayed versus temperature on a logarithmic scale. The approximate location of the inflection point at temperature T i n f ( B ) (versus magnetic field B) is indicated by the arrow. At T = T i n f ( B ) the HF metal α YbAl 1 x Fe x B 4 enters the crossover region, separating the NFL behavior from the LFL one, see Figure 7. The inset displays T i n f versus magnetic field B.
Figure 19. Magnetization M ( T ) versus T for B values shown in the legend. The data are taken from [34]. M ( T ) is displayed versus temperature on a logarithmic scale. The approximate location of the inflection point at temperature T i n f ( B ) (versus magnetic field B) is indicated by the arrow. At T = T i n f ( B ) the HF metal α YbAl 1 x Fe x B 4 enters the crossover region, separating the NFL behavior from the LFL one, see Figure 7. The inset displays T i n f versus magnetic field B.
Materials 15 03901 g019
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shaginyan, V.R.; Msezane, A.Z.; Japaridze, G.S.; Artamonov, S.A.; Leevik, Y.S. Strongly Correlated Quantum Spin Liquids versus Heavy Fermion Metals: A Review. Materials 2022, 15, 3901. https://doi.org/10.3390/ma15113901

AMA Style

Shaginyan VR, Msezane AZ, Japaridze GS, Artamonov SA, Leevik YS. Strongly Correlated Quantum Spin Liquids versus Heavy Fermion Metals: A Review. Materials. 2022; 15(11):3901. https://doi.org/10.3390/ma15113901

Chicago/Turabian Style

Shaginyan, Vasily R., Alfred Z. Msezane, George S. Japaridze, Stanislav A. Artamonov, and Yulya S. Leevik. 2022. "Strongly Correlated Quantum Spin Liquids versus Heavy Fermion Metals: A Review" Materials 15, no. 11: 3901. https://doi.org/10.3390/ma15113901

APA Style

Shaginyan, V. R., Msezane, A. Z., Japaridze, G. S., Artamonov, S. A., & Leevik, Y. S. (2022). Strongly Correlated Quantum Spin Liquids versus Heavy Fermion Metals: A Review. Materials, 15(11), 3901. https://doi.org/10.3390/ma15113901

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop