ORIGINAL RESEARCH article
Sec. Condensed Matter Physics
Volume 10 - 2022 | https://doi.org/10.3389/fphy.2022.843627
Critical Scaling of Entropy and Thermal Drude Weight in Anisotropic Heisenberg Antiferromagnets: A Thermodynamic Quest for Quantum Criticality
- Department of Physics, ChongQing Three Gorges University, Wanzhou, China
Up to now, probing the quantum phase transition (QPT) and quantum critical (QC) phenomena at finite temperatures in one-dimensional (1D) spin systems still lacks an in-depth understanding. Herein, we study the QPT and thermodynamics of 1D spin-1/2 anisotropic Heisenberg antiferromagnetic chains by Green’s function theory. The quantum phase diagram is renormalized by the anisotropy (∆), which manifests a quantum critical point (QCP) hc = 1 + ∆ signaling the transition from gapless Tomonaga–Luttinger liquid (TLL) to gapped ferromagnetic (FM) state, demonstrated by the magnetic entropy and thermal Drude weight. At low temperatures, it is shown that two crossover temperatures fan out a QC regime and capture the QCP via the linear extrapolation to zero temperature. In addition, around QCP, the QC scaling is performed by analyzing the entropy and thermal Drude weight to extract the critical exponents (α, δ, and β) that fulfill the Essamm–Fisher scaling law, which provides a novel thermodynamic means to detect QPT for experiment. Furthermore, scaling hypothesis equations with two rescaled manners are proposed to testify the scaling analysis, for which all the data points fall on a universal curve or two independent branches for the plot against rescaled field or temperature, implying the self-consistency and reliability of the obtained critical exponents.
A quantum phase transition (QPT) takes place at zero temperature, as a result of the quantum fluctuations arising from Heisenberg uncertainty relation . Usually, it manifests that the ground state of a quantum system changes upon tuning an external nonthermal parameter such as pressure or magnetic field to a critical value, which is marked by a quantum critical point (QCP) for a continuous transition. The QC fluctuations will result in an exotic behavior in stark contrast to the conventional gapped or gapless low-lying excitations of materials. Near the QCP, a QC regime emerges at finite temperatures in an extended parameter space attributed to the interplay between quantum and thermal fluctuations . This intriguing region is featured by the absence of energy scales other than temperature as well as the corresponding critical properties of quantum correlations or thermodynamic quantities, which culminate into scaling characteristic and universality [3–6]. In this regard, one-dimensional (1D) quantum spin systems offer field tunability for probing QPT and critical phenomena. A variety of magnetic compounds have been regarded as the 1D antiferromagnetic (AF) Heisenberg chain model [3–9], whose ground state resides in a Tomonaga–Luttinger liquid (TLL) characterized by spinon quasiparticles with gapless elementary excitations [10, 11]. However, the magnetic anisotropy (∆) would make a significant impact on its properties. For ∆ > 1, it turns into an Ising-type model, which is realized in experiment for CoNb2O6 and BaCo2V2O8 [12–14]. A transverse field induces an Ising QPT with gapless quantum criticality and self-duality of QCP . Differently, for ∆ < 1, it becomes an XY-type spin chain, the experimental realization of which is Cs2CoCl4 [16, 17], whose ground state lies in a TLL state. In a longitude magnetic field, it induces a QPT without self-duality of QCP. Herein, we focus on the AF Heisenberg chain for the anisotropy ∆ ≤ 1, whose ground state still resides in a TLL. Up to the critical field Bc = J (1+∆)/gμB, a QPT occurs from TLL to a spin polarized ferromagnetic (FM) state with an excitation gap opened up by the field .
As we know, the magnetic entropy shows anomaly close to the QCP, where two ground states compete with each other and it does not determine which state to be reside in . Besides, the thermal Drude weight Dth is a good indicator signaling gapped (Dth = 0) or gapless (Dth>0) low-lying excitations . In fact, due to that the absolute zero temperature cannot be attained experimentally, the field-induced quantum criticality has been intensively investigated by using the field dependence of magnetization measurement to determine the kink at finite temperature [3, 6, 9, 13, 16], which becomes rounded such that the QCP cannot be measured exactly. The nature of quantum fluctuations near the QCP remains enigmatic. Thus, it is urgent need to provide a new clue to capture QCP for diagnosing QPT at finite temperatures such that the QC scaling becomes rather important, which is one of the cornerstone concepts in modern physics and plays a key role in understanding quantum criticality, namely, the QPT. Near QCP, the physical quantities such as magnetization, magnetic susceptibility, and specific heat are featured by a set of critical exponents and scaling functions [3, 4]. Note that the critical behavior of entropy should be the same as the specific heat. However, whether the thermal Drude weight can be done critical scaling analysis instead of magnetization (ineffective in Ising model) or not? What critical scaling forms to feature quantum criticality? Although several thermodynamic quantities were carried out to do critical scaling analysis for detecting the QPT, it is still challenging to measure the thermal Drude weight as an effective detector of quantum criticality.
To fully assess the universality of quantum criticality, we will demonstrate the scaling behavior of entropy and thermal Drude weight divided by temperature to extract the critical exponents and capture the QCP. On the one hand, such scaling is a direct consequence of the scaled temperature to a certain universal power multiplied by the field that collapse onto a single curve for the plot against a scaling transformation for thermodynamic quantity. On the other hand, as the magnetic field is scaled to a certain universal power multiplied by temperature, the data points of scaling transformation for thermodynamic quantity will fall on two independent branches. In the forthcoming section, we present the model Hamiltonian and Green’s function theory. In Results and Discussion, the renormalized quantum phase diagram and phase crossover behavior are explored, and the field dependence of magnetization for different anisotropies has been tested, which are compared to the experimental observations; the QC scaling behavior is analyzed and discussed. Finally, a conclusion is drawn in Conclusion.
Model Hamiltonian and Method
where J > 0 denotes the AF coupling and ∆ represents the anisotropy. Hereafter, we define the reduced magnetic field
the Hamiltonian (1) becomes
where the operator
The method that we employ is the two-time Green’s function theory. The retarded Green’s function for JW fermions is defined as 
where the subscripts i and j label lattice sites. After the time Fourier transformation, the Green’s function is put into the equation of motion,
It is clearly shown that a rigorous calculation is not available as a result of the Ising interacting quartic terms. For the high-order Green’s function
For further Fourier transformation into k-space, the Green’s function can be expressed as
The integral of the wavevector k extends over the first Brillouin zone. Thus, the momentum space Green’s function
Then, the average magnetization M per site, specific heat, and thermal entropy are defined as
Hence, the isothermal magnetic entropy change is expressed as 
Results and Discussion
In what follows, kB = ħ = 1 and J = 1 is set as an energy unit; hereafter, all other parameters are scaled by it. At first sight, without loss of generality, the field dependence of magnetization at low temperatures is calculated for different anisotropies and done a comparison to the experimental results. For Δ = 0.25, the system is an XY-like antiferromagnet with gapless TLL. The magnetization increases gradually as the field ascends and reaches its saturation at Bc = J (1 + ∆)/gμB = 2.07T for J = 0.23 meV and g = 2.4, as shown in Figure 1A, which is in accordance with the experimental observation value Bc = 2.1T on compound Cs2CoCl4 . Nonetheless, it becomes an isotropic case with Δ = 1, for which a series of compounds have been identified as candidates in experiment [6–9]. For example, a copper-containing coordination polymer CuPzN (a good realization of a spin-1/2 AF chain) , whose magnetization is fitted by J = 10.6 K and g = 2.27 at low temperatures in Figure 1B, manifests a relatively sharp kink around Bc = 13.9T, close to the experimental, Bethe–Ansatz and QTM results. Beyond Bc, the magnetization reaches the saturation value 1.135 μB per Cu2+, in excellent agreement with 1.15 μB obtained from the experimental observation . As the temperature ascends, the sharp critical signatures broaden systematically due to the enhanced thermal fluctuations. As it is well known, the absolute zero temperature cannot be attained experimentally such that the critical field cannot be measured exactly. How to feature the universality of quantum criticality and phase crossover with gapped or gapless low-lying excitations? How to capture the exact critical field at finite temperature to diagnose the QPT? In addition to the abovementioned magnetization characterization, on the one hand, the magnetic entropy is a good quantity to characterize the quantum criticality by its maximum value at the lowest temperature because two quantum phases compete with each other at the QCP, where it is not determined which ground state to be resided in. Herein, the magnetic field is renormalized by the anisotropy λ = h/(1 + ∆). Figure 1C presents the ∆-λ phase diagram by the contour plot of entropy. It is clearly shown that the renormalized critical field λc = 1.0 associated with a sharp peak of entropy separates the TLL and FM phases for any anisotropy. At finite temperature, the different states are featured by the magnetic entropy change in Figure 1D, in which the inverse magnetocaloric effect (IMCE) (−∆S < 0) predominates in TLL at low temperature for λ < λc, while only conventional magnetocaloric effect (CMCE) (−∆S > 0) persists in FM state beyond λc. Furthermore, it has been pointed out that |ΔS| follows a power–law dependence of the field: |ΔS|∼λn with
demonstrated by the linear double logarithm with δ = 2.009 (see the inset in Figure 1F), which agrees well with the experimental value 1.98 obtained from the critical scaling of magnetization on CuPzN .
FIGURE 1. At low temperatures, the field dependence of magnetization for ∆ = 0.25 in (A) (J, g) = (0.23 meV, 2.4) and for ∆ = 1.0 in (B) (J, g) = (10.6 K, 2.27). (C) The renormalized quantum phase diagram by the contour plot of entropy at T = 0.005 with Tomonaga–Luttinger liquid (TLL) and spin polarized ferromagnetic (FM) phase being unveiled. For ∆ = 0.25, (D) the temperature dependence of magnetic entropy change under different fields. (E) The field dependence of magnetic entropy change; the inset is its log-log scale under different temperatures. (F) The thermal Drude weight landscape of renormalized quantum phase diagram; the inset is its scaling around the quantum critical point.
However, at finite temperatures, the QPT disappears, and it only shows a phase crossover behavior. Proceeding the same way as that observed in Sr3Ru2O7 , the field derivative of the entropy divided by temperature is employed to characterize the landscape of temperature–field phase diagram, as shown in Figure 2A. A peak-valley structure marks two crossover temperatures fanning out a QC regime, which constitute two branches of the phase crossover boundaries. In the QC regime, it is well known that the correlation length ξ diverges as
with α = 0.497 demonstrated by the slope of linear lnS/T–lnT curve in the inset, which turns out to be a T1/2 behavior of the entropy. In fact, one can obtain α = 1/2 from the scaling relation α = 2-(d + z)/z with dynamic exponent z = 2 and spatial dimension d = 1 [3, 4, 34].
FIGURE 2. Temperature–field phase diagram based on the contour plot of the field derivative of (A) entropy and (D) thermal Drude weight divided by temperature, fanning out a quantum critical (QC) regime. The temperature dependence of (B) entropy and (E) thermal Drude weight under different fields; the insets are the logarithm as a function of 1/T, revealing the spin excitation gap by the magnitude of slope. (C) S/T and (F) Dth/T as a function of temperature at the critical field; the insets are the corresponding double logarithm plot, unveiling the critical exponents.
For a comparison, the field derivative of thermal Drude weight divided by temperature is also employed to characterize the landscape of temperature–field phase diagram, as shown in Figure 2D, which shows the similar behavior as the entropy plotted in Figure 2A. A double-valley structure marks two crossover temperatures fanning out a QC regime, which also feature the phase crossover boundaries. Under different fields, the thermal Drude weight as a function of temperature behaves differently, as shown in Figure 2E. As T→0, one can find that the thermal Drude weight displays a T-linear relation for λ < λc, implying a gapless behavior, while it is exponentially activated:
Thus, one can find that α+β(1+δ) = 1.968 fulfills the Essam–Fisher relation 
In addition, we propose some scaling hypothesis equations to confirm the reliability and self-consistency of the obtained critical exponents. From the above analysis, it is clearly shown that the critical exponents δ and β describe the field and temperature dependence of thermal Drude weight
FIGURE 3. Scaling transformation of (A) entropy and (B) thermal Drude weight as a function of rescaled field at different temperatures, while they are plotted as a function of rescaled temperature around the critical field in (C) and (D), respectively. All the data points collapse onto a universal curve or two independent branches, indicating the self-consistency and reliability of the obtained critical exponents.
The other is the temperature rescaled, i.e., the temperature multiplied by a scaled magnetic field. Differently, in the asymptotic QC region, the thermal Drude weight equation is proposed as
in which f(x) is a regular function that behaves differently for λ < λc and λ > λc, respectively. Thus, as T→0, it allows an unambiguous determination of the critical exponent δ from the thermal Drude weight
In conclusion, the QPT and low-temperature properties of 1D spin-1/2 anisotropic Heisenberg AF chains are investigated by means of Green’s function theory. For different anisotropies, the field dependence of magnetization is calculated at low temperatures, which are in good agreement with the experimental results. We further renormalize the quantum phase diagram by the anisotropy that manifests a gapless TLL transition into gapped FM state at hc = 1 + ∆, which is demonstrated by the drops to zero from a finite value of thermal Drude weight. At low temperature, two crossover temperatures fan out a QC regime and capture the QCP from linear high-temperature extrapolation to zero temperature. The T-linear dependence of entropy and thermal Drude weight signals the gapless low-lying excitation in TLL, while it decays exponentially upon cooling down to zero temperature, suggesting the gapped behavior in FM state. At the QCP, it takes on a T1/2 or T3/2 behavior. Furthermore, we demonstrate the QC scaling via analyzing the entropy and thermal Drude weight around QCP to extract the critical exponents (α, δ, and β) that fulfill the Essamm–Fisher relation α+β(1+δ) = 2, which provides a novel thermodynamic means to detect QPT for experiment. Meanwhile, scaling hypothesis equations with two scaling transformation manners are proposed to check the scaling analysis. One is the plot against the rescaled magnetic field, for which all the data points collapse onto a universal single curve, whereas the other is the rescaled temperature such that the data points fall on two independent branches, indicating the self-consistency and reliability of the obtained critical exponents, which provide an explicit physical picture for understanding QC phenomena at finite temperatures.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
LD conceived the work and wrote the manuscript. YZ performed some theoretical calculation results. All the authors participated in the discussion and interpretation of the results.
This work was supported by the Natural Science Foundation of Chongqing (cstc2021jcyj-msxmX0878) and the Youth Project of Scientific and Technological Research Program of Chongqing Education Commission (KJQN202101210).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
3. Kono Y, Sakakibara T, Aoyama CP, Hotta C, Turnbull MM, Landee CP, et al. Field-induced Quantum Criticality and Universal Temperature Dependence of the Magnetization of a Spin-1/2 Heisenberg Chain. Phys Rev Lett (2015) 114:037202. doi:10.1103/physrevlett.114.037202
6. Breunig O, Garst M, Klümper A, Rohrkamp J, Turnbull MM, Lorenz T. Quantum Criticality in the Spin-1/2 Heisenberg Chain System Copper Pyrazine Dinitrate. Sci Adv (2017) 3:eaao3773. doi:10.1126/sciadv.aao3773
7. Wolf B, Tsui Y, Jaiswal-Nagar D, Tutsch U, Honecker A, Removic-Langer K, et al. Magnetocaloric Effect and Magnetic Cooling Near a Field-Induced Quantum-Critical point. Proc Natl Acad Sci (2011) 108:6862–6. doi:10.1073/pnas.1017047108
9. Chakraborty T, Singh H, Chaudhuri D, Jeevan HS, Gegenwart P, Mitra C. Investigation of Thermodynamic Properties of Cu(NH 3 ) 4 SO 4 ·H 2 O, a Heisenberg Spin Chain Compound. J Magnetism Magn Mater (2017) 439:101–6. doi:10.1016/j.jmmm.2017.05.020
12. Liang T, Koohpayeh SM, Krizan JW, McQueen TM, Cava RJ, Ong NP. Heat Capacity Peak at the Quantum Critical point of the Transverse Ising Magnet CoNb2O6. Nat Commun (2015) 6:7611. doi:10.1038/ncomms8611
13. Wang Z, Lorenz T, Gorbunov DI, Cong PT, Kohama Y, Niesen S, et al. Quantum Criticality of an Ising-like Spin- 1/2 Antiferromagnetic Chain in a Transverse Magnetic Field. Phys Rev Lett (2018) 120:207205. doi:10.1103/physrevlett.120.207205
14. Wang Z, Schmidt M, Loidl A, Wu JD, Zou HY, Yang W, et al. Quantum Critical Dynamics of a Heisenberg-Ising Chain in a Longitudinal Field: Many-Body Strings versus Fractional Excitations. Phys Rev Lett (2019) 123:067202. doi:10.1103/physrevlett.123.067202
16. Kenzelmann M, Coldea R, Tennant DA, Visser D, Hofmann M, Smeibidl P, et al. Order-to-disorder Transition in theXY-like Quantum magnetCs2CoCl4induced by Noncommuting Applied fields. Phys Rev B (2002) 65:144432. doi:10.1103/physrevb.65.144432
17. Breunig O, Garst M, Sela E, Buldmann B, Becker P, Bohatý L, et al. Spin-12XXZChain SystemCs2CoCl4in a Transverse Magnetic Field. Phys Rev Lett (2013) 111:187202. doi:10.1103/physrevlett.111.187202
18. Rost AW, Perry RS, Mercure J-F, Mackenzie AP, Grigera SA. Entropy Landscape of Phase Formation Associated with Quantum Criticality in Sr 3 Ru 2 O 7. Science (2009) 325:1360–3. doi:10.1126/science.1176627
24. Franco V, Blázquez JS, Ipus JJ, Law JY, Moreno-Ramírez LM, Conde A. Magnetocaloric Effect: From Materials Research to Refrigeration Devices. Prog Mater Sci (2018) 93:112–232. doi:10.1016/j.pmatsci.2017.10.005
29. Hess C, ElHaes H, Waske A, Büchner B, Sekar C, Krabbes G, et al. Linear Temperature Dependence of the Magnetic Heat Conductivity in CaCu2O3. Phys Rev Lett (2007) 98:027201. doi:10.1103/physrevlett.98.027201
30. Biswas A, Chandra S, Samanta T, Ghosh B, Datta S, Phan MH, et al. Universality in the Entropy Change for the Inverse Magnetocaloric Effect. Phys Rev B (2013) 87:134420. doi:10.1103/physrevb.87.134420
31. Biswas A, Chandra S, Samanta T, Phan MH, Das I, Srikanth H. The Universal Behavior of Inverse Magnetocaloric Effect in Antiferromagnetic Materials. J Appl Phys (2013) 113:17A902. doi:10.1063/1.4793768
32. Yang B, Chen Y-Y, Zheng Y-G, Sun H, Dai H-N, Guan X-W, et al. Quantum Criticality and the Tomonaga-Luttinger Liquid in One-Dimensional Bose Gases. Phys Rev Lett (2017) 119:165701. doi:10.1103/physrevlett.119.165701
33. Adelhardt P, Koziol JA, Schellenberger A, Schmidt KP. Quantum Criticality and Excitations of a Long-Range Anisotropic XY Chain in a Transverse Field. Phys Rev B (2020) 102:174424. doi:10.1103/physrevb.102.174424
Keywords: quantum criticality, one-dimensional spin system, thermodynamics, critical scaling, magnetic anisotropy
Citation: Ding LJ and Zhong Y (2022) Critical Scaling of Entropy and Thermal Drude Weight in Anisotropic Heisenberg Antiferromagnets: A Thermodynamic Quest for Quantum Criticality. Front. Phys. 10:843627. doi: 10.3389/fphy.2022.843627
Received: 26 December 2021; Accepted: 10 January 2022;
Published: 10 February 2022.
Edited by:Guoying Gao, Huazhong University of Science and Technology, China
Reviewed by:Jiyu Fan, Nanjing University of Aeronautics and Astronautics, China
Sun Zhao-Yu, Wuhan Polytechnic University, China
Copyright © 2022 Ding and Zhong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: L. J. Ding, email@example.com